Wastewater Based Epidemiology Enabled Surveillance of Antibiotic Resistance
===========================================================================

* M.V. Riquelme
* E. Garner
* S. Gupta
* J. Metch
* N. Zhu
* M.F. Blair
* G. Arango-Argoty
* A. Maile-Moskowitz
* A. Li
* C-F Flach
* D.S. Aga
* I. Nambi
* D.G.J. Larsson
* H. Bürgmann
* T. Zhang
* A. Pruden
* P.J. Vikesland

## ABSTRACT

Wastewater-based epidemiology (WBE) for disease monitoring is highly promising, but requires consistent methodologies that incorporate predetermined objectives, targets, and metrics. We demonstrate a comprehensive metagenomics-based approach for global surveillance of antibiotic resistance in sewage, enabling assessment of: 1) which antibiotic resistance genes (ARGs) are shared across regions/communities; 2) which ARGs are discriminatory; and 3) factors associated with overall trends including antibiotic concentrations in sewage. Across an internationally-sourced transect of sewage samples collected using a centralized, standardized protocol, ARG relative abundances (16S rRNA gene-normalized) were highest in Hong Kong and India and lowest in Sweden and Switzerland, reflecting national policy, measured antibiotic concentrations, and metal resistance genes. Asian versus European/US resistomes were distinct, with macrolide-lincosamide-streptogramin, phenicol, quinolone, and tetracycline versus multidrug resistance ARGs being discriminatory, respectively. Sales data were not predictive of antibiotics measured in sewage, emphasizing need for direct measurements. The WBE approach defined herein demonstrates multi-site comparability and sensitivity to local/regional factors.

Wastewater-based epidemiology (WBE) is a rapidly emerging framework for public health surveillance as it quickly and non-invasively provides anonymous population-scale information about human disease and anthropogenic chemical use.[1] Activity in this field has focused on xenobiotic and human biomarker measurements[2, 3] and the monitoring of infectious diseases such as polio[4] and typhoid.[5] Most recently, WBE has been proposed for surveillance of the spread of antibiotic resistance[6-8] as well as COVID-19.[9, 10] Within any given sewershed, the analysis of temporal or spatial changes in WBE targets (e.g., genes, microbes, chemicals) could provide an early warning of local and regional disease outbreaks,[11] while comparisons across sewersheds can enable insights into environmental and socioeconomic factors contributing to disease patterns.[12] To produce globally actionable information, however, data comparability across disparate systems is crucial and should be considered at the outset of commencing measurements. In this study, we applied a consistent sample collection and analysis protocol[13, 14] and demonstrated the capability of an array of data analysis approaches[15-17] for discriminating an international transect of sewage samples from Asia (Hong Kong, India, the Philippines), Europe (Sweden, Switzerland), and North America (US) with the aim of advancing a WBE framework for antibiotic resistance surveillance.

Sewage represents a composite of human-associated flora; including pathogens, antibiotic resistant bacteria (ARB), and antibiotic resistance genes (ARGs),[18-21] carried across the corresponding community and thus provides opportunities for environmental antibiotic resistance surveillance.[6, 22] The profiles of ARB and ARGs in sewage are expected to be dependent on current and historic antibiotic resistance management practices, including the types and quantities of antibiotics used and local attention to transmission control.[18, 23] A targeted and well-designed WBE program could help establish baseline levels of antibiotic resistance for a given community, inform effective antibiotic use policy and wastewater treatment practices, and serve as a point of comparison for assessing the impacts of such interventions.[12, 22]

Shotgun metagenomic sequencing provides a powerful means to characterize sewage microbiomes[24] and viromes,[25] producing nucleotide sequences that can be archived and compared to publicly-available databases to profile pathogens, pandemic viruses,[26, 27] and ARG composition (i.e., the “resistome”).[28-30] Recent studies have applied metagenomic sequencing to identify compelling sewage resistome trends. Hendriksen et al.[6] collected wastewater samples from 79 sites in 60 countries and found that ARG abundance correlated with socioeconomic, health, and environmental factors. Using the same dataset, Karkman et al. explored to what extent sewage metagenomes (alone or in combination with socioeconomic factors) predict local clinical resistance.[8] Pärnänen et al.[7] showed that sewage samples collected across a European North-South transect correlated with differences in antibiotic usage and average local temperatures. Going forward, precise surveillance objectives, targets, and metrics must be identified and validated if metagenomics enabled WBE is to be globally recommended and adopted.[22] In the context of antibiotic resistance, identification of informative features of resistomes and their relationship to other key sewage constituents, including concentrations of residual antibiotics, metals, metal resistance genes (MRGs), and mobile genetic elements (MGEs), is needed to gain insight to the factors driving resistance proliferation (**Figure 1**).

![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/06/03/2021.06.01.21258164/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2021/06/03/2021.06.01.21258164/F1)

Figure 1. 
General strategy for antibiotic resistance monitoring.

Here we assessed various dimensions of sewage metagenomes, antibiotic concentrations, and the associated metainformation, to identify discriminatory features and factors that might reflect the corresponding resistomes. Taxonomic composition and the diversity of the corresponding microbiomes were evaluated as key constraints on ARG composition.[31, 32] In addition, plasmids, transposons, and integrons were quantified as plausible carriers of resistance[28, 29, 33] (i.e., indicating the potential for associated ARGs to spread horizontally across species). Genes that could provide co-selection opportunities, such as MRGs, were also quantified, addressing concerns that even if antibiotic use is curbed, metals and other selective agents could still select for ARGs. Finally, “core resistome”[31, 33, 34] and “discriminatory resistome”[16] analyses were carried out to account for the ARG “background” of globally-distributed or naturally-prevalent ARGs (e.g., those found even in 30,000 year old Arctic soils[35] and permafost[36]), and to differentiate ARGs that result primarily from anthropogenic activity and are of clinical concern. Herein, we advance a WBE framework for sewage-based surveillance of antibiotic resistance, focusing monitoring towards targets relevant to local evolution and dissemination.

## RESULTS AND DISCUSSION

### Comparison of Resistomes Across a Global Transect of Sewage Samples

Fourteen influent sewage samples from twelve WWTPs located in six countries were systematically collected by our team using standardized protocols,[13, 14] profiled, and compared. To identify and quantify ARGs, we performed shotgun metagenomic sequencing (20-30 million reads per sample; **Supplementary Table S1**) and annotated reads using the CARD database.[37] We detected an average of 449 ARGs (Range: 309-489) at each site, with the Chao index estimating true ARG richness ranging from 501-683 ARGs/sample (Average: 577; **Supplementary Figure S1**). When considering 16S rRNA gene normalized relative abundances, a clear trend emerged, whereby Asian sewage samples contained higher relative abundances of ARGs than European/US samples (**Figure 2A, Supplementary Figures S2** and **S3, Supplementary Data 1**). Such a pattern is likely driven primarily by the relative inputs of ARGs and ARBs from human populations, although other factors could be at play, including differences in relative industrial/hospital inputs or differential selection pressures when sewage travels through the collection network.[38, 39]

![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/06/03/2021.06.01.21258164/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2021/06/03/2021.06.01.21258164/F2)

Figure 2. 
**A**: ARG distribution and relative abundance by corresponding antibiotic categories. ARGs were annotated via MetaStorm using CARD version 1.2.1. ARG categories were assigned in-house (**Supplementary Data 7**). Genes corresponding to two or more categories were labeled as “multidrug.” ARG abundances were normalized via MetaStorm to 16S rRNA gene abundances. **B:** Non-metric multidimensional scaling (NMDS) ordination of WWTPs according to ARG-based Bray-Curtis distance. Ellipses enclose sites of noted similarities, ranging from 0 to 100 (perfect similarity); thus, ellipses with similarity of 80, represent sites with highest similarity **C**: Relative abundance of top 20 bacterial phyla in WWTP influents. Genus-level annotations were done via MG-RAST using the Refseq database. **D:** NMDS ordination of WWTPs according to Euclidian distances between rlog-transformed reads to bacterial genera. Ellipses enclose sites of noted distances, where the shortest distance denotes the highest resemblance and longer distance denotes the least resemblance among the sites. Sample names refer to countries (IND: India, PHL: Philippines, USA: United States, CHE: Switzerland, HKG: Hong Kong, SWE: Sweden) along with the visit number (1,2) and WWTP plant number (P).

Hendriksen et al.[6] recently reported a global sewage survey and, although the analytical approach was distinct, the trend in ARG relative abundance was strikingly similar to the present study: higher in Africa (median: 2,100 fragments per kilo base per million fragments [FPKM]), Asia (1,200 FPKM), South America (1,900 FPKM), and the Middle East (1,100 FPKM), relative to Europe (750 FPKM), North America (900 FPKM), and Oceania (800 FPKM). Hendriksen et al.[6] applied FPKM normalization to address highly variable sequencing library depth (8-398 million reads per sample), while in the present study (20-30 million reads per sample) we normalized to 16S rRNA genes as a biomarker of total bacteria and correspondingly indicate the proportion of bacteria carrying ARGs.[40] Hendriksen et al. applied the ResFinder database[41] for ARG annotation, which resulted in detection of an apparently larger number of ARGs. This difference may result from the inclusion of multiple similar variants of the same ARG in ResFinder, but not in CARD, or from the inclusion of multiple ARGs that result from detection of a single ARG cluster. The differential number of ARGs included in different databases reflects the rapid evolution of the field of resistome bioinformatics[42] and the need to improve comparability across studies for sewage surveillance.[43]

We visualized ARG compositional similarities via non-metric multidimensional scaling (NMDS) ordination derived from Bray-Curtis dissimilarity matrices (**Figure 2B**). Clustering within the NMDS plot was tightest for WWTP influents originating from the same country, with the exception of the two Philippines WWTPs. A second level of clustering was observed between the Asian and European/US samples. One exception was the HKG-P1 influent, which clustered separately from all other samples. This latter result is consistent with the saline nature of this sewage, where ocean water is used for toilet flushing. Analysis of similarity (ANOSIM) indicated that sewage samples strongly separated from one another when grouped by location (R=0.741, p=0.001). Clusters also strongly separated when grouped as Asia versus Europe/US (R=0.700; p=0.002). Three clusters, less strongly separated (R=0.586; p=0.001), were apparent when grouped by continent. Importantly, two sets of temporally independent influent samples collected from two Hong Kong WWTPs (HKG1-P1/HKG2-P1 vs. HKG1-P2/HKG2-P2) clustered separately, but reproducibly, thus indicating less variation with time relative to geographic location. Ju et al.[31] also noted a high level of repeatability when characterizing resistomes in sewage samples from different WWTPs in Switzerland. Other studies[6, 44] have similarly shown locational reproducibility, which is key for inferring that observed differences are primarily driven by geographical location rather than the particular time of sample collection.

A number of clinically relevant ARGs (**Supplementary Data 2)** were detected in the influent of all WWTPs. OXA-type carbapenemase ARGs were ubiquitous and were detected at relatively high abundance (0.03-0.16 copies/16S rRNA gene copies). A variety of other β-lactamase ARG types (GES, CARB, SHV, CTX-M, and TEM) were also broadly detected, although at lower abundances. The *qnr*S gene (quinolone resistance) was detected in all sewages (0.001-0.24 copies/16S rRNA gene copies) and tended towards higher abundance in Asia (0.085±0.099 copies / 16S rRNA gene copies) relative to Europe/US (0.012±0.011 copies/16S rRNA gene copies), but not significantly so (Wilcox, p=0.1649). Others have noted differences in *qnr* carriage across populations,[45, 46] but statistical power here may have been insufficient to detect a difference. Reads annotated as the mcr*-1* colistin ARG were found at all WWTPs, with the exception of one Hong Kong plant.

### Core and Discriminatory Aspects of Resistomes

We subjected the collected data to “core” and “discriminatory” resistome analysis to further delineate similarities and differences between sites. ARGs and resistance classes that were ubiquitous across all sites reflect the “core”, while those that enable locational differentiation reflect “discriminatory” resistomes. The core resistome of 216 ARGs (**Supplementary Data 3**) was comprised of ARGs conferring multidrug resistance and resistance to aminoglycosides, β-lactams, macrolide-lincosamide-streptogramins (MLS), polymyxins, tetracyclines, and trimethoprims. The 25 most abundant ARGs found across all locations included one polymyxin (*pmr*E), one trimethoprim (*dfr*E), two tetracycline (*tet*Q and *tet*C), seven MLS (*msr*E, *mph*D, *erm*F, *mel, mac*B, *mef*A, *erm*B), one sulfonamide (*sul*1), three multidrug (*CRP, msb*A, *ade*J), two quinolone (*qnrS2, qac*H), one β-lactam (*cfx*A6) and two aminoglycoside (*aph(3’’)-lb* and *aph(6)-ld*) ARGs (**Supplementary Table S2**). Notably, among the core ARGs for which information was available, many were “ancient” or “natural background” ARGs identified in prior studies of pristine or less human-impacted environments (e.g., isolated caves, glacial soil, permafrost); thus, these ARGs likely predate the antibiotic era **(Supplementary Data 3)**. However, other core ARGs are known to have become globally distributed more recently (e.g., *sul*1, *sul*2, *dfrA*3, *tet*(G), *aad*A8). Based on pairwise comparison, the number of ARGs shared between sites ranged from 296 to 401 (**Supplementary Figure S4** with 349 average shared annotations. The overall proportion of core ARGs averaged 69% (min=62%, IND1-P2; max=72%, HKG2-P1). Notably, the least number of reads in common were between one of the Swedish (SWE1-P1) versus one of the Hong Kong sewages (HKG2-P1), which aligns with the Asia versus Europe/US divide in sewage resistome similarity. Consistent with the repeatability of the sampling events, the greatest number of common ARGs was noted for two Hong Kong sewage samples collected six months apart.

To focus analysis on the “variable” resistome, we removed core ARGs from the profiles in **Figure 2A**. This accentuated the continental differences in the sewage resistomes (**Figure 3**). Specifically, the relative proportions of aminoglycoside, β-lactam, rifamycin, sulfonamide, trimethoprim, and quinolone ARGs tended to increase, while those of MLS, multidrug, polymyxin, and tetracycline ARGs tended to decrease. Interestingly, the proportion of core ARGs tended to increase from Asia to Europe/US, which hypothetically could reflect a longer history of antibiotic use and management in Europe/US.

![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/06/03/2021.06.01.21258164/F3.medium.gif)

[Figure 3.](http://medrxiv.org/content/early/2021/06/03/2021.06.01.21258164/F3)

Figure 3. 
Distribution and relative abundance of ARGs comprising the “variable” resistome, i.e., with “core” ARGs removed from the analysis, by corresponding antibiotic categories.

We next characterized the “discriminatory” resistomes (i.e., ARGs that are, in terms of presence or abundance, indicative of local population or region). We applied the ExtrARG machine learning-based algorithm[16] to identify ARGs that discriminate sewage samples based on *a priori* selected groups. Sewage samples were grouped at the continent level (i.e., Asia vs. Europe/US) and then the relative abundances of the 50 most discriminatory ARGs were visualized using a heatmap (**Figure 4**). Two aminoglycoside ARGs, *aad*A8 and *aad*A2, were prominent in Asia, while ant*(3”)-IIb, aac(3)-Ia* and *axy*X were more prominent in Europe/US. Various β-lactam ARGs (*bla*CphA5, *bla*OXA-363, *bla*OXA-309, *bla*OXA-371, *bla*FOX-8, *bla*FOX-1, *bla*CPS-1) were abundant in Europe/US, but were found at low concentrations in Asian sewage. Interestingly, *bla*CARB-2 and *bla*CARB-4 were not detected in European/US sewage, but were found in Asia. The MLS ARG *erm*T, phenicol ARGs *cml*A1 and *cml*A4, quinolone ARG *qnr*VC6, and the tetracycline ARG *tet*(L) were primarily detected in Asian sewage. The sulfonamide ARG, *sul*2, was 3× to 5× more abundant in Asian sewage, while a number of the multidrug ARGs (e.g., efflux pumps *ade*J, *mex*K, *abe*M), were 3× to 15× higher in sewage from Europe/US. Further, we also found that a number of ARGs identified as discriminatory were only discovered within the last decade (e.g., ant*(3”)-IIb, axy*X, *bla*CphA5, *bla*CPS-1, *bla*FOX-8, *qnr*VC6 and *qnr*VC6; **Supplementary Data 4**). We suggest that the location-specific distribution of discriminatory ARGs in sewage can be used to determine if their global distribution is widening. Such a result would suggest potential clinical concern.

![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/06/03/2021.06.01.21258164/F4.medium.gif)

[Figure 4.](http://medrxiv.org/content/early/2021/06/03/2021.06.01.21258164/F4)

Figure 4. 
Heatmap representing the relative abundance of discriminatory ARGs for influent samples grouped according to their sampling continent identified using the ExtrARG python package.

### Comparison of Microbiomes Across the Global Transect of Sewage Samples

Similar to the resistomes, NMDS analysis of the sewage microbiomes (i.e., bacterial taxonomic composition derived from metagenomic data) resulted in clustering of samples as a function of geographic location (**Figure 2C** and **2D; Supplementary Figure S5**). ANOSIM indicated that the sewage microbiomes were most strongly separated when grouped by country (R=0.701, p=0.001). When grouped by Asia versus Europe/US, the groups were better separated (R=0.421, p=0.001) than when grouped by continent (i.e., Asia vs. Europe vs. US; R=0.3; p=0.003). As was the case for the resistomes, Hong Kong WWTP1 (HKG-P1) clustered separately from the other sewage samples. Similar patterns in the clustering of the resistomes and microbiomes are expected, since phylogenetic constraints to ARG carriage are well known.[31]

In accordance with the ARG-based class separations for Asia versus Europe/US, LEfSe analysis was performed on the microbiome data (**Supplementary Figure S6**). Four and ten classifying genera were identified in European/US versus Asian sewages, respectively. European/US sewages were distinguished by genera belonging to *Betaproteobacteria*, while *Clostridia* or *Negativicute* classes of the *Firmicute* phylum were discriminatory in Asian sewage. Only two non-*Firmicute*s genera (*Dehalococcoides* and *Shigella*) were found to be discriminatory in Asian sewage. Typically, higher levels of *Proteobacteria* are thought to reflect environmental conditions within the sewer system, while *Firmicutes* reflect greater human fecal influence.[47, 48]

Potential parallels between the microbiome and resistome compositions were investigated. First, the Mantel test was applied to either full or continent-specific distance matrices representative of the microbiome (Euclidian distances between log-transformed reads to bacterial genera) and resistome (Bray-Curtis dissimilarity matrix). When comparing all the collection sites, a strong positive correlation was observed between the ARG and genus-based taxonomy distance matrices (r = 0.74, p < 0.001). When separated by continent, a stronger correlation was observed between the Asian distance matrices (r = 0.80, p < 0.001), than between the European/USA distance matrices (r = 0.48, p = 0.03). For comparative purposes, a Bray-Curtis dissimilarity matrix was generated based on MetaPhlAn[49] output of taxonomic assignments obtained from MetaStorm[15] for Mantel-based comparison to the ARG similarity matrix. Again, a strong correlation was observed, both overall (r = 0.76, p < 0.001) and when compared continentally within Asia (r = 0.80, p < 0.001) and Europe/US (r = 0.56, p < 0.01).

Common pathogenic hosts of the discriminatory ARGs identified in **Figure 4** were also examined (**Supplementary Table S3**). Based on the data reported in the CARD database,[37] eighty-two pathogenic bacterial species are known to host these ARGs. Interestingly, discriminatory ARGs most representative of Asian sewages tended to have a broader range of potential host pathogens, whereas discriminatory ARGs associated with European/US sewages tended to have a narrower range.

### Genes Providing Opportunities for Co-Selection and Horizontal Gene Transfer

Genes documented to sometimes be found on plasmids were investigated as indicators of horizontal gene transfer potential. The average relative abundance of such genes identified via alignment to the ACLAME database was 56.7 copies/16S rRNA gene copies, with an observed minimum of 33.3 copies/16S rRNA gene copies in the Philippines (PHL1-P1) and a maximum of 111.6 copies/16S rRNA gene copies in Hong Kong (HKG2-P1; **Supplementary Figure S7**). No clear country- or continent-specific patterns were observed in terms of abundance.

Potential associations between ARGs and MRGs known to sometimes occur on plasmids or other MGEs were investigated and visualized via network analysis of assembled scaffolds (**Figure 5**). While acknowledging that there is uncertainty in the accuracy of scaffolds assembled from metagenomic data,[50] we carried out the analysis for empirical comparison assuming a consistent error rate across a sample set generated from the same sequencing platform. A substantial portion of reads were assembled, averaging 34% and ranging from 24% (HKG1-P1) to 50% (PHL1-P2). For each sewage sample, the number of generated scaffolds ranged from 64,576 (USA1-P1) to 243,356 (CHE1-P2), with an average length of 719 bp. All of the major classes of ARGs were detected on the scaffolds, with the most frequently observed belonging to the multidrug, MLS, glycopeptide, and tetracycline categories (**Supplementary Figure S8**). The relative distribution of ARGs (based on antibiotic category) among the scaffolds was remarkably similar for all sewage samples. Across all locations, the most common ARGs with plasmid associations were *pmr*E (peptide ARG), *acr*B (multidrug efflux), *dfr*E (trimethoprim ARG), *Mux*B (multidrug efflux), and *ros*B (peptide ARG) with 301, 186, 133, 125, and 101 co-occurrences, respectively. *pmrE* and *dfr*E were present at high abundance at all of the sites, thus suggesting globalization of these genes. Other notable plasmid ARG co-occurrences include *sul*1 (sulfonamide ARG, 20 co-occurrences), *MCR-*3 (peptide, 32 co-occurrences), and the *bla*CARB, *bla*OXA, and *bla*TEM β-lactam ARGs (between 1-13 co-occurrences). *sul*1 is known to be highly associated with both class 1 integrons[51] as well as plasmids. Some of these co-occurrences correspond to resistance against carbapenems and colistin, which are critically-important antibiotics of last resort. This analysis illustrates that clinically-important classes of ARGs were readily found on assembled scaffolds that also contained genes commonly found on plasmids, an indicator that they are potentially in a mobile form in sewage.

![Figure 5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/06/03/2021.06.01.21258164/F5.medium.gif)

[Figure 5.](http://medrxiv.org/content/early/2021/06/03/2021.06.01.21258164/F5)

Figure 5. 
Co-occurrence of ARGs and plasmid gene markers on *de novo* assembled scaffolds generated by shotgun metagenomic sequencing reads pooled from all samples. This analysis highlights ARGs that are probable candidates for potential horizontal gene transfer (co-occurrences with plasmid associated genes) or co-selection (co-occurrences with ARGs of different classes). Proximity of nodes and width of lines indicate frequency of associations between genes. Node diameter is proportional to the number of co-occurrences for that gene. Co-occurences with fewer than 3 instances were excluded from the network analysis rendering. MLS = macrolide, lincosamide, and streptogramin resistance.

ARGs and MRGs are often subject to co-selection pressure when they are both present on a single genetic element, such as a plasmid, or cross-selection if the same gene is both an ARG and an MRG, as is commonly the case for multidrug efflux pumps capable of pumping both metals and antibiotics.[52] The ranking of total MRG relative abundances, identified via alignment to the BacMet database,[53] was strikingly similar to that of the total ARG relative abundances (Spearman r = 0.66, p = 0.013), except both Indian WWTPs and one Philippines plant ranked much lower than the other Asian WWTPs, while a Swiss WWTP had extremely low MRG relative abundance (**Figure 6**). This finding is suggestive that there are, or have been, common drivers, such as metals, antibiotics or other agents, either within human gastrointestinal tracts, in industrial inputs, or in sewers themselves, that exert (or have exerted) selection pressure for carriage of both ARGs and MRGs.

![Figure 6.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/06/03/2021.06.01.21258164/F6.medium.gif)

[Figure 6.](http://medrxiv.org/content/early/2021/06/03/2021.06.01.21258164/F6)

Figure 6. 
Distribution and relative abundance of metal resistance genes in WWTP influents. Metal resistance genes were annotated via MetaStorm using the BacMet database version 1.1. Meta resistance gene abundances were normalized to 16S rRNA gene abundances. Sample names refer to countries (IND: India, PHL: Philippines, USA: United States, CHE: Switzerland, HKG: Hong Kong, SWE: Sweden) along with the visit number and WWTP plant (P) number. Order is according to the ranked comparison of total ARG relative abundance shown in Figure 1.

Across all locations, multidrug ARGs were the most common type associated with MRGs (**Supplementary Figure S9**). The most common co-occurrences within the assembled scaffolds were noted between *acr*B and *Mex*B ARGs, co-located with copper/zinc MRGs (160 and 83 co-occurrences, respectively), and *mux*B ARG co-located with zinc MRGs (96 co-occurrences). These co-locations were consistently among the most abundant at all sample sites. A literature review was conducted to examine whether any of these presumed co-occurrences were actually the same gene, and thus indicative of cross-rather than co-selection. We identified nine genes that were consistently abundant across all the sample sites and confer resistance to both antibiotics and metals, namely, *mdt*A, *mdt*B, *mdt*C, *bae*S, *bae*R, *cme*B, *acr*D, *mex*I and *pmr*C. As anticipated, most of these genes are known to be associated with efflux pump systems (**Supplementary Data 5**).

### Sewage Antibiotic Concentrations

Antibiotic, pharmaceutical, and personal care product (PPCP) concentrations measured in these sewage samples were recently reported in detail by Singh et al.[14] and trends relevant to the present study are summarized in **Table 1** and **Supplementary Figure S10**. We note that chemical stability issues precluded quantification of the β-lactams and that, due to logistical constraints, the antibiotic measurements for the Indian sewage were made during a sampling trip that occurred one year later. We measured the highest antibiotic concentrations in Hong Kong sewage, where total concentrations of macrolides, quinolones, sulfonamides, tetracyclines, and trimethoprim reached levels >63,000 ng/L. The antibiotics with the highest concentrations in the sewage were ciprofloxacin (48,100 ng/L, HKG1-P1) and clarithromycin (39,551 ng/L, IND2-P2). In particular, Indian sewage contained very high levels of anhydro-erythromycin and norfloxacin. To contextualize these measured antibiotic concentrations, we first compared them to predicted no effect concentrations (PNECs; **Table 1**)[54] for resistance selection. All samples, except one from Sweden, were above the PNEC for ciprofloxacin, while norfloxacin was near or above the PNEC in both Swiss sewages and one Indian sewage sample. Clarithromycin was above the PNEC in one US sample. Two Hong Kong sewage samples had concentrations greater than the PNEC for the tetracyclines, while both the Hong Kong and Indian sewages contained the macrolides azithromycin and clarithromycin at levels exceeding the PNEC. It is noted that PNEC values provide a point of reference, but measured concentrations that either exceed or fall short of a PNEC do not necessarily delineate the presence or absence of selection pressure within sewage.[54] Resistance metrics quantified in sewage likely reflect both current and historic drivers, including past patterns of antibiotic use and selection pressures that may have occurred in the human gut, that collectively shape the resistome; however, antibiotic measurements can only reflect usage at the time of sampling.

View this table:
[Table 1:](http://medrxiv.org/content/early/2021/06/03/2021.06.01.21258164/T1)

Table 1: 
Proposed no effect concentrations (PNEC),52 representing the cutoff below which no selection of resistant bacteria is anticipated to occur and concentrations (ng/L) of antibiotics detected in the influent of each WWTP. Concentration data originally published in Singh et al. **13**

To further contextualize the measured antibiotic concentrations, we compared them to values predicted based upon antibiotic sales data. Prior studies have used sales data to explain global differences in resistomes in sewage[6] and other contexts.[55] Antibiotic consumption data between 2000-2015 for each site was obtained from ResistanceMap.[56] As shown in **Supplementary Figure S11**, total antibiotic consumption (defined daily doses (DDD) per 1000 population) was reported to be highest in the US, followed by Hong Kong, Sweden, and Switzerland.[55] Reported antibiotic use in India and the Philippines has historically been lower than that in Europe/US, but in recent years, reported use in India has increased and is estimated to be similar to that in Sweden and Switzerland as of the 2015. Reported antibiotic use in the Philippines lags relative to other countries and has remained fairly stable over the 2000-2015 period.

We calculated *expected* antibiotic concentrations within our sewage samples based upon sales data, data on pharmaceutical excretion rates, and published relationships between xenobiotic analyte concentrations and sewershed populations. As illustrated in **Supplementary Figure S10**, these estimated antibiotic concentrations overestimate the measured antibiotic concentrations by orders of magnitude. Such overestimates hold true regardless of the locality or the xenobiotic used. Similar overestimates were also found using estimated sewershed populations and sewage flow rates employed. These disparities are unsurprising given the many potential disconnects at the national, local, and regional level between antibiotic sales, consumption, and disposal patterns[57, 58] as well as the wide differences in biological and chemical decay rates exhibited by these antibiotics.[59] Considering that measured antibiotic data better reflected measured trends in the resistome, (i.e., Hong Kong and Indian sewage ranked highest in total ARG relative abundance; **Figure 2A**), we recommend that determination of actual antibiotic concentrations in sewage should take precedence over extrapolations from general sales or consumption reports. Because antibiotic use often follows seasonal patterns, ideally such measurements should take place regularly throughout the year. It is recognized, however, that some antibiotics (e.g., beta lactams) have short half-lives and/or recoveries and thus, while better, direct measurement of antibiotics in sewage is not a perfect solution.

We also explored if prior reports of links between antibiotic concentrations and resistant clinical isolates58,59 held true for sewage resistomes. Bonferroni-corrected Spearman rank order correlation analyses were conducted to identify potential relationships between concentrations of antibiotics or pharmaceuticals and PPCPs, individual ARGs, and ARG classes (**Supplementary Data 6**). At the individual ARG level, few correlations were observed when comparing against antibiotic compound concentrations: clarithromycin concentration showed a positive correlation with *tet*(L) (≥ 0.93, p<0.01) and a negative correlation with *mph*B (≥ -0.91, p<0.05); and norfloxacin showed a positive correlation with the *cat* gene (≥ 0.91, p<0.05). At the ARG class level, macrolide antibiotics positively correlated with quinolone ARGs (≥ 0.84, p<0.05). Of the four macrolide antibiotics tested, clarithromycin presented the strongest correlations with quinolone ARGs. A remarkably similar pattern of macrolide antibiotics and quinolone ARG correlation was previously observed in Swiss sewage.30 Overall, few correlations between ARGs and antibiotics were observed and those that were observed pertained to macrolides. Lack of simplistic correlations between antibiotics and ARGs is consistent with prior work[6, 19] and is to be expected considering that the subset of antibiotics analyzed do not capture the full range of antibiotics actually present, while certain dominant ARGs may not be subject to selection pressure and others may be negatively selected by certain antibiotics. Also, effective concentrations of antibiotics may be too low to exert measurable selection pressure in the sewage or selection pressure might have occurred historically in the human gut and thus there is a spatiotemporal disconnect between the two measures.

### Surveillance Effectively Reveals Key Sewage Resistome Trends

We observed that the relative abundances of total ARGs were elevated in Asian versus European/US sewages, consistent with general perceptions of rigor in antibiotic use policy across countries[19] and the generally higher antibiotic loadings in Asia. The lowest relative abundances of ARGs were found in the sewage from Sweden which has incorporated numerous proactive measures to avoid misuse and overuse of antibiotics.[60] Hong Kong, India, and the Philippines, on the other hand, are characterized by dense urban populations and either have problems with illegal antibiotic sales without a prescription[61] or do not require prescriptions for antibiotic use.[57]

Removal of core ARGs (i.e., those ubiquitously detected in urban sewage) from the resistome profiles accentuated the differences between Asia and Europe/US, while delineation of discriminatory ARGs identified those that uniquely circulate in each corresponding region. In general, a greater variety of ARGs across multiple resistance classes were identified as discriminatory in Asian sewage, while a subset of aminoglycoside and multi-drug ARGs were characteristic of European/US sewages. Such knowledge provides information about the ARGs within local regions in terms of their prevalence, novelty, or their potential carriage by pathogens. Surveillance is a key component of the global action plan to combat antimicrobial resistance.[62] Within that context, consistent methods and approaches are required to facilitate robust data comparability (**Figure 1**). The approach outlined herein provides such consistency, but it would costly to apply everywhere due to its reliance on a dedicated sampling and analysis team. The alternative approach,[6] relying upon centralized processing of large numbers of shipped samples, has the capacity to more broadly elucidate global trends, while more focused efforts, such as that delineated herein, provide the potential for greater granularity (e.g., comparing local sewersheds, sampling within a sewershed). Implementation of coordinated sewage resistome surveillance efforts will be key to improved understanding of baseline resistome characteristics.[42] Such information can be used as a reference point for future sampling to aid in early identification of previously unobserved ARGs that may pose clinical threats, including the emergence of new resistance types.[22]

### Implications for Wastewater-Based Epidemiology

WBE is rapidly gaining attention as a highly promising tool for infectious disease monitoring.[2, 6, 7, 12, 33, 63] Here we demonstrated the potential for metagenomic-based surveillance of ARGs and other key metainformation in sewage to provide a foundation for WBE of antibiotic resistance, and we have illustrated challenges that must be considered as the field develops. In particular, it is critical to define WBE surveillance objectives; including what will be monitored, how samples (and appropriate metainformation) will be collected and analyzed, and how the data will ultimately be interrogated prior to investing in extensive sampling campaigns. The hypothesis that the sewage resistome composition differs by geographical location could only be supported through the use of consistent experimental and data analysis protocols, which will be particularly vital in future efforts to further dissect the role of geographical sub-factors, such as socioeconomics, population density, antibiotic use, diet, sanitation status, and local climate/temperature.[19, 20] Comparison of results obtained using widely varying methods (e.g., nucleic acid extraction kits,[13] solid-phase extraction protocols[14]) or metagenomic platforms present inherent and substantial limitations to identifying key factors driving the spread of antibiotic resistance.

For WBE to be applied towards estimating the prevalence of infectious disease or antibiotic resistance for a given region, reasonably accurate estimates of the population contributing to a sewage sample are needed.[2, 64] We utilized xenobiotic analyte measurements to develop estimates of the population within a given sewershed. These measurements enabled estimation of the *de facto*[64] population (i.e., incorporating transient contributors to the sewershed) at the time of sampling. A required input to these calculations is the sewage flow rate. Such metainformation is not necessarily collected in all studies and the lack of such information for the sewage samples from the Philippines precluded estimation of antibiotic concentrations for that country. An additional consideration of such population estimates is the assumption that relationships developed within Australia relating PPCP concentrations to population are globally applicable. Given disparities in antibiotic and PPCP use across the world this latter assumption requires evaluation. Incorporation of antibiotic and PPCP measurements within the context of WBE requires *a priori* planning if such measurements are to be used to estimate study populations within the sewershed. Given these considerations, alternative approaches to estimate human contributions to sewage (e.g., quantification of pepper mottle mild virus[65] or crAssphage[21]) may be preferable.

The WBE approach defined herein demonstrates a means to identify geographic “hot-spots” where attention towards curbing the spread of antibiotic resistance may be particularly warranted. This is especially relevant as we emerge from the new era of pandemic awareness and concern, where multidrug resistant pathogens, ARGs, and infectious viruses have been observed to rapidly spread across international borders.[66] Estimates that several thousands of deaths per year globally can be attributed to antibiotic resistance amplifies the need for coordinated, global surveillance. Ultimately, widescale online implementation of metagenomic surveillance of ARGs and other disease determinants within sewage in the context of a WBE framework promises to yield new insights. Such insights will help illuminate the specific factors that drive the dissemination and spread not only of antibiotic resistance, but also other diseases of emerging concern.

## METHODS

### Sample locations

A total of 14 sewage samples were collected from 12 WWTPs located in India, Hong Kong, the Philippines, Sweden, Switzerland and the US, with 2 WWTPs per location and the Hong Kong WWTPs subjected to two separate sampling events. The WWTP capacities ranged from 2.6-66 MGD with a majority of the sewage consisting of domestic wastewater, with some industrial and hospital inputs (**Table 2**).

View this table:
[Table 2:](http://medrxiv.org/content/early/2021/06/03/2021.06.01.21258164/T2)

Table 2: 
Sewage sample characteristics

### Sample collection and processing

Sample collection and processing were standardized as described previously.[13, 14] In brief, influent grab samples were collected in sterile polypropylene containers and transported to a local laboratory on ice. Sewage biomass was concentrated in triplicate via vacuum filtration of equal influent volumes onto three separate 0.22 µm mixed cellulose ester membranes (MilliporeSigma Darmstadt, Germany), which were then preserved in 50% ethanol and stored at -20 °C. Ethanol-fixed membrane samples were shipped to Virginia Tech for DNA extraction and biomolecular analyses. Sewage samples (0.5 L) for chemical analysis 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 isotopically labeled surrogate standards (50 µL of 1000 µg/L surrogate mix solution) were added to each sample.

### DNA extraction, shotgun metagenomic sequencing and analyses

As described previously,[13] DNA extraction was carried out using the FastDNA Spin Kit for Soil (MP Biomedicals, Solon, OH). Composite samples were prepared by pooling equal masses of triplicate samples. TrueSeq libraries (Illumina, San Diego, CA) were prepared for shotgun metagenomic sequencing via Illumina HiSeq 2500 with 2 × 100 paired-end reads. Sequencing was performed at the Virginia Tech Biocomplexity Institute Genomic Sequencing Center (Blacksburg, VA). Library preparation and sequencing were duplicated independently for two samples to verify sequencing reproducibility. After quality trimming, sequencing depth ranged from 10 to 16 million paired-end reads per sample with an average of 13 million reads per sample (**Supplementary Table S1**).

### Taxonomic, ARG, plasmid, and metal resistance gene (MRG) analyses

Sequences were uploaded to the metagenomics RAST server (MG-RAST) and annotated against the RefSeq database to identify bacterial phyla and genera using default parameters.[67] Taxonomic beta diversity was examined via NMDS based ordination according to Euclidian distances between rlog-transformed reads to bacterial genera.[68] ARG annotation was carried out via MetaStorm[15] with read matching to CARD (protein homolog model) version 1.2.1.[69] The pipeline uses the DIAMOND BLASTX71 aligner with the representative hit approach (E-value<1e-10, identity>90%, and minimum length of 25 amino acids). The ARGs in CARD version 1.2.1 were assigned into resistance categories based on the updated ARGminer database.[17] Manual curation was subsequently performed to fill in missing categories or to update existing categories and to segregate the genes thought to confer resistance via point mutations. Gene names and the ARG category to which they were assigned are found in **Supplementary Data 7**. If an ARG conferred resistance to two or more antibiotics classes, it was categorized as “multidrug.” The core resistome is defined as the collection of ARGs that are ubiquitous and are an inherent part of the resistome. Core resistome analysis was carried out by subsetting the list of ARGs that were identified in all samples (**Supplementary Data 3**). Plasmid-associated sequences were annotated using the ACLAME database, version 0.4.[70] MRGs were annotated using the list of experimentally confirmed metal resistance genes from BacMet, antibacterial biocide and metal resistance genes database, version 1.1.[53] Gene abundances were normalized within MetaStorm to 16S rRNA gene abundances as previously described.[40] NMDS and analysis of similarities (ANOSIM) was conducted considering each individual ARG detected (i.e., not grouped by class) with square root transformation and Bray-Curtis dissimilarities. R-value cutoffs as defined by Clarke and Warwick[71] were used (R>0.75, well separated; 0.75>R>0.25, separated but overlapping; R<0.25, barely separated). Clinically-relevant ARGs were selected from a manually curated list prioritizing human health impacts of ARGs carried by infectious bacteria (**Supplementary Table S3**). Discriminatory ARGs were identified via the ExtrARG[16] machine learning based approach. ExtrARG is based on the extremely randomized trees classifier that uses a Bayesian strategy to optimize the parameters of the classifier and identify discriminatory ARGs based on the user-defined categorizing scheme or groups. LEfSe[72] enables biomarker discovery among class conditions within metagenomic samples. It first performs a non-parametric Kruskal-Wallis sum-rank test to detect features with differential abundance among the conditions of interest; it then applies an unpaired Wilcoxon rank-sum test to test for biological consistency. LEfSe uses linear discriminant analysis to estimate the effect size of each differentially abundant feature. Reads were assembled *de novo* in MetaStorm according to default parameters and scaffolds were mapped against the CARD and ACLAME databases. Network visualization was conducted using Gephi (version 0.8.2). Potential pathogenic ARG hosts were determined using the CARD prevalence tool, which examines the occurrence of each ARG of interest on NCBI archived pathogen chromosomes and assembled whole genome sequences, as well as plasmids.[37]

### Antibiotic analysis

Sample processing and analysis for antibiotics was conducted as previously described.[14] Solid phase extraction (SPE) was performed on 1 L filtered wastewater samples by conditioning Oasis HLB cartridges with acetonitrile and nanopure water before the water samples were loaded at a rate of 3-5 mL/min. Cartridges were dried under vacuum and then shipped to the University at Buffalo for elution and liquid chromatography with tandem mass spectrometry (LC-MS/MS) analysis. LC-MS/MS analysis was carried out using an Agilent 1200 LC system (Palo Alto, CA).

### Antibiotic loading

Anticipated antibiotic loads in sewage were calculated using the following expression: ![Formula][1]</img>  where C*i,sewage*(g/L) is the expected concentration of compound *i* in sewage, DDD (g/day) is the defined daily dose per person, P is the estimated sewershed population, F is the fraction of antibiotic excreted, and Q is the reported wastewater flow. Sewershed populations are notoriously challenging to obtain; however, they may be estimated based upon known excretion rates for widely used PPCPs such as carbamazepine, caffeine, and iopromide. For those systems where we detected one or more of these PPCPs, we estimated the sewershed population using: ![Formula][2]</img>  where C*PPCP,sewage* is the measured concentration of carbamazepine (CBZ), caffeine (CAF), or iopromide (IOP) and β*PPCP* is an empirical parameter relating sewage PPCP concentration to population. β*CBZ* = 8.6 (± 1.4 at 95% confidence level) × 10−8 kg/day-person; β*CAF* = 6.8 (± 3.2 at 95% confidence level) × 10−7 kg/day-person; β*IOP* = 1.1 (± 0.19 at 95% confidence level) × 10−5kg/day-person. To date, β*PPCP* values have only been reported for samples collected within Australia and thus this calculation inherently assumes it is universally applicable across all samples. Such an assumption requires future verification.

## Supporting information

Supplemental Figures and Tables [[supplements/258164_file02.pdf]](pending:yes)

Supplementary Data 1 [[supplements/258164_file03.xlsx]](pending:yes)

Supplementary Data 2 [[supplements/258164_file04.xlsx]](pending:yes)

Supplementary Data 3 [[supplements/258164_file05.xlsx]](pending:yes)

Supplementary Data 4 [[supplements/258164_file06.xlsx]](pending:yes)

Supplementary Data 5 [[supplements/258164_file07.xlsx]](pending:yes)

Supplementary Data 6 [[supplements/258164_file08.xlsx]](pending:yes)

Supplementary Data 7 [[supplements/258164_file09.xlsx]](pending:yes)

## Data Availability

The wastewater influent metagenomic datasets have been deposited in NCBI Short Read Archive (SRA) under bioproject PRJNA527877.

## Data availability

The wastewater influent metagenomic datasets have been deposited in NCBI Short Read Archive (SRA) under bioproject PRJNA527877.

## Supporting Information

**SI:** Supplementary Figures (S1-S11) and Tables (S1-S3).

**Supplementary Data 1**: 16S Normalized Gene Counts

**Supplementary Data 2**: Clinically Relevant ARGs

**Supplementary Data 3**: Common ARGs

**Supplementary Data 4**: Discriminatory ARGs

**Supplementary Data 5**: Cross Resistance Genes

**Supplementary Data 6**: Antibiotic correlations

**Supplementary Data 7**: Manually curated database of ARGs included in analyses.

## Footnotes

*   * apruden{at}vt.edu

*   Received June 1, 2021.
*   Revision received June 1, 2021.
*   Accepted June 3, 2021.


*   © 2021, Posted by Cold Spring Harbor Laboratory

The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission.

## References

1.  1.Sims N, Kasprzyk-Hordern B. Future perspectives of wastewater-based epidemiology: Monitoring infectious disease spread and resistance to the community level. Environ Int. 2020;139:105689. Epub 2020/04/14. doi: 10.1016/j.envint.2020.105689. PubMed PMID: 32283358; PubMed Central PMCID: PMCPMC7128895.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.envint.2020.105689&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32283358&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 

2.  2.Choi PM, Tscharke BJ, Donner E, O’Brien JW, Grant SC, Kaserzon SL, et al. Wastewater-based epidemiology biomarkers: Past, present and future. TrAC-Trends in Analytical Chemistry. 2018;105:453–69. doi: 10.1016/j.trac.2018.06.004. PubMed PMID: WOS:000440408800036.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.trac.2018.06.004&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=WOS:00044040&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 

3.  3.Zuccato E, Chiabrando C, Castiglioni S, Bagnati R, Fanelli R. Estimating community drug abuse by wastewater analysis. Environ Health Perspect. 2008;116(8):1027–32. Epub 2008/08/19. doi: 10.1289/ehp.11022. PubMed PMID: 18709161; PubMed Central PMCID: PMCPMC2516581.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1289/ehp.11022&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18709161&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000258270200023&link_type=ISI) 

4.  4.Hovi T, Shulman LM, van der Avoort H, Deshpande J, Roivainen M, Em DEG. Role of environmental poliovirus surveillance in global polio eradication and beyond. Epidemiol Infect. 2012;140(1):1–13. Epub 2011/08/19. doi: 10.1017/S095026881000316X. PubMed PMID: 21849095.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/S095026881000316X&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21849095&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000298547400001&link_type=ISI) 

5.  5.Sikorski MJ, Levine MM. Reviving the “Moore Swab”: a Classic Environmental Surveillance Tool Involving Filtration of Flowing Surface Water and Sewage Water To Recover Typhoidal Salmonella Bacteria. Applied and Environmental Microbiology. 2020;86(13):e00060–20.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1128/aem.00060-20&link_type=DOI) 

6.  6.Hendriksen RS, Munk P, Njage P, van Bunnik B, McNally L, Lukjancenko O, et al. Global monitoring of antimicrobial resistance based on metagenomics analyses of urban sewage. Nat Commun. 2019;10(1):1124. Epub 2019/03/10. doi: 10.1038/s41467-019-08853-3. PubMed PMID: 30850636; PubMed Central PMCID: PMCPMC6408512.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-019-08853-3&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30850636&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 

7.  7. Pärnänen Kmm, Narciso-da-Rocha C, Kneis D, Berendonk TU, Cacace D, Do TT, et al. Antibiotic resistance in European wastewater treatment plants mirrors the pattern of clinical antibiotic resistance prevalence. Science Advances. 2019;5(3):eaau9124. doi: 10.1126/sciadv.aau9124
    
    [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6MzoiUERGIjtzOjExOiJqb3VybmFsQ29kZSI7czo4OiJhZHZhbmNlcyI7czo1OiJyZXNpZCI7czoxMjoiNS8zL2VhYXU5MTI0IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDYvMDMvMjAyMS4wNi4wMS4yMTI1ODE2NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

8.  8.Karkman A, Berglund F, Flach CF, Kristiansson E, Larsson DGJ. Predicting clinical resistance prevalence using sewage metagenomic data. Communications Biology. 2020;3(1). doi: 10.1038/s42003-020-01439-6. PubMed PMID: WOS:000595962500002.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s42003-020-01439-6&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=WOS:00059596&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 

9.  9.Medema G, Been F, Heijnen L, Petterson S. Implementation of environmental surveillance for SARS-CoV-2 virus to support public health decisions: opportunities and challenges. Curr Opin Environ Sci Health. 2020. Epub 2020/10/08. doi: 10.1016/j.coesh.2020.09.006. PubMed PMID: 33024908; PubMed Central PMCID: PMCPMC7528975.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.coesh.2020.09.006&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33024908&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 

10. 10.Polo D, Quintela-Baluja M, Corbishley A, Jones DL, Singer AC, Graham DW, et al. Making waves: Wastewater-based epidemiology for COVID-19 - approaches and challenges for surveillance and prediction. Water Res. 2020;186:116404. Epub 2020/09/18. doi: 10.1016/j.watres.2020.116404. PubMed PMID: 32942178; PubMed Central PMCID: PMCPMC7480445.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.watres.2020.116404&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32942178&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 

11. 11.Daughton CG. Monitoring wastewater for assessing community health: Sewage Chemical-Information Mining (SCIM). Sci Total Environ. 2018;619–620:748-64. Epub 2017/11/22. doi: 10.1016/j.scitotenv.2017.11.102. PubMed PMID: 29161600; PubMed Central PMCID: PMCPMC6091531.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.scitotenv.2017.11.102&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29161600&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 

12. 12.Huijbers PMC, Flach CF, Larsson DGJ. A conceptual framework for the environmental surveillance of antibiotics and antibiotic resistance. Environ Int. 2019;130. doi: 10.1016/j.envint.2019.05.074. PubMed PMID: WOS:000477934800004.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.envint.2019.05.074&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=WOS:00047793&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 

13. 13.Li AD, Metch JW, Wang YL, Garner E, Zhang AN, Riquelme MV, et al. Effects of sample preservation and DNA extraction on enumeration of antibiotic resistance genes in wastewater. FEMS Microbiology Ecology. 2018;94(2). doi: 10.1093/femsec/fix189. PubMed PMID:WOS:000428122100016.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/femsec/fix189&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=WOS:00042812&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 

14. 14.Singh RR, Angeles LF, Butryn DM, Metch JW, Garner E, Vikesland PJ, et al. 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:361–9. doi: 10.1016/j.envint.2019.01.025. PubMed PMID: WOS:000457122700036.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.envint.2019.01.025&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=WOS:00045712&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 

15. 15.Arango-Argoty G, Singh G, Heath LS, Pruden A, Xiao WD, Zhang LQ. MetaStorm: A Public Resource for Customizable Metagenomics Annotation. Plos One. 2016;11(9). doi: 10.1371/journal.pone.0162442. PubMed PMID: WOS:000383706900048.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0162442&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=WOS:00038370&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 

16. 16.Gupta S, Arango-Argoty G, Zhang L, Pruden A, Vikesland P. Identification of discriminatory antibiotic resistance genes among environmental resistomes using extremely randomized tree algorithm. Microbiome. 2019;7(1):123. Epub 2019/08/31. doi: 10.1186/s40168-019-0735-1. PubMed PMID: 31466530; PubMed Central PMCID: PMCPMC6716844.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s40168-019-0735-1&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31466530&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 

17. 17.Arango-Argoty GA, Guron GKP, Garner E, Riquelme MV, Heath LS, Pruden A, et al. ARGminer: A web platform for crowdsourcing-based curation of antibiotic resistance genes. Bioinformatics. 2019;36(9):2966–73. doi: 10.1093/bioinformatics/btaa095.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btaa095&link_type=DOI) 

18. 18.Bürgmann H, Frigon D, Gaze WH, Manaia CM, Pruden A, Singer AC, et al. Water and sanitation: an essential battlefront in the war on antimicrobial resistance. FEMS Microbiology Ecology. 2018;94(9). doi: 10.1093/femsec/fiy101. PubMed PMID: WOS:000442373800002.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/femsec/fiy101&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=WOS:00044237&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 

19. 19.Collignon P, Beggs JJ, Walsh TR, Gandra S, Laxminarayan R. Anthropological and socioeconomic factors contributing to global antimicrobial resistance: a univariate and multivariable analysis. The Lancet Planetary Health. 2018;2(9):e398–e405. doi: 10.1016/s2542-5196(18)30186-4.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/s2542-5196(18)30186-4&link_type=DOI) 

20. 20.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. Epub 2019/03/09. doi: 10.1021/acs.accounts.8b00643. PubMed PMID: 30848890.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1021/acs.accounts.8b00643&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30848890&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 

21. 21.Karkman A, Parnanen K, Larsson DGJ. Fecal pollution can explain antibiotic resistance gene abundances in anthropogenically impacted environments. Nat Commun. 2019;10(1):80. Epub 2019/01/10. doi: 10.1038/s41467-018-07992-3. PubMed PMID: 30622259; PubMed Central PMCID: PMCPMC6325112.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-018-07992-3&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30622259&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 

22. 22.Aarestrup FM, Woolhouse M. Using sewage for surveillance of antimicrobial resistance. Science. 2020;367:630–2.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzNjcvNjQ3OC82MzAiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wNi8wMy8yMDIxLjA2LjAxLjIxMjU4MTY0LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

23. 23.Hutinel M, Huijbers PMC, Fick J, Ahren C, Larsson DGJ, Flach CF. Population-level surveillance of antibiotic resistance in Escherichia coli through sewage analysis. Euro Surveillance. 2019;24(37). doi: 10.2807/1560-7917.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2807/1560-7917&link_type=DOI) 

24. 24.Wu L, Ning D, Zhang B, Li Y, Zhang P, Shan X, et al. Global diversity and biogeography of bacterial communities in wastewater treatment plants. Nat Microbiol. 2019;4(7):1183–95. Epub 2019/05/16. doi: 10.1038/s41564-019-0426-5. PubMed PMID: 31086312.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41564-019-0426-5&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31086312&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 

25. 25.Bibby K, Crank K, Greaves J, Li X, Wu Z, Hamza IA, et al. Metagenomics and the development of viral water quality tools. npj Clean Water. 2019;2(1). doi: 10.1038/s41545-019-0032-3.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41545-019-0032-3&link_type=DOI) 

26. 26.Izquierdo-Lara R, Elsinga G, Heijnen L, Oude Munnink BB, Schapendonk CME, Nieuwenhuijse D, et al. Monitoring SARS-CoV-2 circulation and diversity through community wastewater sequencing. medRxiv. 2020. doi: 10.1101/2020.09.21.20198838.this.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1101/2020.09.21.20198838.this&link_type=DOI) 

27. 27.Crits-Christoph A, Kantor RS, Olm MR, Whitney ON, Al-Shayeb B, Lou YC, et al. Genome Sequencing of Sewage Detects Regionally Prevalent SARS-CoV-2 Variants. mBio. 2021;12(1):e02703–20. doi: 10.1128/mBio.02703-20.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1128/mBio.02703-20&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33468686&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 

28. 28.Wang Z, Zhang XX, Huang K, Miao Y, Shi P, Liu B, et al. Metagenomic profiling of antibiotic resistance genes and mobile genetic elements in a tannery wastewater treatment plant. PLoS One. 2013;8(10):e76079. Epub 2013/10/08. doi: 10.1371/journal.pone.0076079. PubMed PMID: 24098424; PubMed Central PMCID: PMCPMC3787945.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0076079&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24098424&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 

29. 29.Guo JH, Li J, Chen H, Bond PL, Yuan ZG. Metagenomic analysis reveals wastewater treatment plants as hotspots of antibiotic resistance genes and mobile genetic elements. Water Research. 2017;123:468–78. doi: 10.1016/j.watres.2017.07.002. PubMed PMID: WOS:000410010500046.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.watres.2017.07.002&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=WOS:00041001&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 

30. 30.D’Costa VM, McGrann KM, Hughes DW, Wright GD. Sampling the Antibiotic Resistome. Science. 2006;311:374–7. doi: 10.1126/science.1120800.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzMTEvNTc1OS8zNzQiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wNi8wMy8yMDIxLjA2LjAxLjIxMjU4MTY0LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

31. 31.Ju F, Beck K, Yin X, Maccagnan A, McArdell CS, Singer HP, et al. Wastewater treatment plant resistomes are shaped by bacterial composition, genetic exchange, and upregulated expression in the effluent microbiomes. ISME Journal. 2019;13(2):346–60. Epub 2018/09/27. doi: 10.1038/s41396-018-0277-8. PubMed PMID: 30250051; PubMed Central PMCID: PMCPMC6331547.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41396-018-0277-8&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30250051&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 

32. 32.Narciso-da-Rocha C, Rocha J, Vaz-Moreira I, Lira F, Tamames J, Henriques I, et al. Bacterial lineages putatively associated with the dissemination of antibiotic resistance genes in a full-scale urban wastewater treatment plant. Environ Int. 2018;118:179–88. Epub 2018/06/09. doi: 10.1016/j.envint.2018.05.040. PubMed PMID: 29883764.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.envint.2018.05.040&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29883764&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 

33. 33.Su JQ, An XL, Li B, Chen QL, Gillings MR, Chen H, et al. Metagenomics of urban sewage identifies an extensively shared antibiotic resistome in China. Microbiome. 2017;5(1):84. Epub 2017/07/21. doi: 10.1186/s40168-017-0298-y. PubMed PMID: 28724443.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s40168-017-0298-y&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28724443&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 

34. 34.Munck C, Albertsen M, Telke A, Ellabaan M, Nielsen PH, Sommer MOA. Limited dissemination of the wastewater treatment plant core resistome. Nature Communications. 2015;6. doi: 10.1038/ncomms9452. PubMed PMID: WOS:000363143400002.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ncomms9452&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=WOS:00036314&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 

35. 35.Allen HK, Moe LA, Rodbumrer J, Gaarder A, Handelsman J. Functional metagenomics reveals diverse beta-lactamases in a remote Alaskan soil. Isme Journal. 2009;3(2):243-51. PubMed PMID: WOS:000263261800011.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=WOS:00026326&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 

36. 36.D’Costa VM, King CE, Kalan L, Morar M, Sung WW, Schwarz C, et al. Antibiotic resistance is ancient. Nature. 2011;477(7365):457–61. doi: 10.1038/nature10388. PubMed PMID: 21881561.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature10388&link_type=DOI) 
    
    [GeoRef](http://medrxiv.org/lookup/external-ref?access_num=2012003519&link_type=GEOREF) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21881561&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000295080500039&link_type=ISI) 

37. 37.Jia B, Raphenya AR, Alcock B, Waglechner N, Guo P, Tsang KK, et al. CARD 2017: expansion and model-centric curation of the comprehensive antibiotic resistance database. Nucleic Acids Res. 2017;45(D1):D566–D73. doi: 10.1093/nar/gkw1004. PubMed PMID: 27789705; PubMed Central PMCID: PMCPMC5210516.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkw1004&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27789705&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 

38. 38.Hu Y, Yang X, Qin J, Lu N, Cheng G, Wu N, et al. Metagenome-wide analysis of antibiotic resistance genes in a large cohort of human gut microbiota. Nature Communications. 2013;4:2151. Epub 2013/07/24. doi: 10.1038/ncomms3151. PubMed PMID: 23877117.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ncomms3151&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23877117&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 

39. 39.Pormohammad A, Nasiri MJ, Azimi T. Prevalence of antibiotic resistance in Escherichia coli strains simultaneously isolated from humans, animals, food, and the environment: a systematic review and meta-analysis. Infect Drug Resist. 2019;12:1181–97. Epub 2019/06/14. doi: 10.2147/IDR.S201324. PubMed PMID: 31190907; PubMed Central PMCID: PMCPMC6512575.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2147/IDR.S201324&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31190907&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 

40. 40.Li B, Yang Y, Ma L, Ju F, Guo F, Tiedje JM, et al. Metagenomic and network analysis reveal wide distribution and co-occurrence of environmental antibiotic resistance genes. Isme Journal. 2015;9(11):2490–502. doi: 10.1038/ismej.2015.59. PubMed PMID: WOS:000365091700015.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ismej.2015.59&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=WOS:00036509&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 

41. 41.Zankari E, Hasman H, Cosentino S, Vestergaard M, Rasmussen S, Lund O, et al. Identification of acquired antimicrobial resistance genes. Journal of Antimicrobial Chemotherapy. 2012;67(11):2640–4. doi: 10.1093/jac/dks261.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/jac/dks261&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22782487&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000309923800014&link_type=ISI) 

42. 42.Hendriksen RS, Bortolaia V, Tate H, Tyson GH, Aarestrup FM, McDermott PF. Using Genomics to Track Global Antimicrobial Resistance. Frontiers in Public Health. 2019;7:242. Epub 2019/09/26. doi: 10.3389/fpubh.2019.00242. PubMed PMID: 31552211; PubMed Central PMCID: PMCPMC6737581.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fpubh.2019.00242&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31552211&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 

43. 43.Martinez JL, Coque TM, Baquero F. What is a resistance gene? Ranking risk in resistomes. Nature Reviews Microbiology. 2015;13(2):116–23. PubMed PMID: WOS:000348144800011.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nrmicro3399&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25534811&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 

44. 44.Joseph SM, Battaglia T, Maritz JM, Carlton JM, Blaser MJ. Longitudinal Comparison of Bacterial Diversity and Antibiotic Resistance Genes in New York City Sewage. Msystems. 2019;4(4). doi: 10.1128/mSystems.00327-19. PubMed PMID: WOS:000482629400020.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1128/mSystems.00327-19&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=WOS:00048262&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 

45. 45.Colomer-Lluch M, Calero-Cáceres W, Jebri S, Hmaied F, Muniesa M, Jofre J. Antibiotic resistance genes in bacterial and bacteriophage fractions of Tunisian and Spanish wastewaters as markers to compare the antibiotic resistance patterns in each population. Environ Int. 2014;73:167–75.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.envint.2014.07.003&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25127043&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 

46. 46.Rutgersson C, Fick J, Marathe N, Kristiansson E, Janzon A, Angelin M, et al. Fluoroquinolones and qnr genes in sediment, water, soil, and human fecal flora in an environment polluted by manufacturing discharges. Environ Sci Technol. 2014;48(14):7825–32. Epub 2014/07/06. doi: 10.1021/es501452a. PubMed PMID: 24988042.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1021/es501452a&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24988042&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 

47. 47.Bengtsson-Palme J, Hammaren R, Pal C, Ostman M, Bjorlenius B, Flach CF, et al. Elucidating selection processes for antibiotic resistance in sewage treatment plants using metagenomics. Science of the Total Environment. 2016;572:697–712. doi: 10.1016/j.scitotenv.2016.06.228. PubMed PMID: WOS:000387807200068.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=doi:10.1016/j.scitotenv.2016.06.228&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=WOS:00038780&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 

48. 48.McLellan SL, Huse SM, Mueller-Spitz SR, Andreishcheva EN, Sogin ML. Diversity and population structure of sewage-derived microorganisms in wastewater treatment plant influent. Environmental Microbiology. 2010;12(2):378–92. doi: 10.1111/j.1462-2920.2009.02075.x. PubMed PMID: WOS:000274234600008.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1462-2920.2009.02075.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19840106&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 

49. 49.Segata N, Waldron L, Ballarini A, Narasimhan V, Jousson O, Huttenhower C. Metagenomic microbial community profiling using unique clade-specific marker genes. Nature methods. 2012;9(8):811–4. doi: 10.1038/nmeth.2066.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nmeth.2066&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22688413&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000307015700022&link_type=ISI) 

50. 50.Brown CL, Keenum IM, Dai D, Zhang L, Vikesland PJ, Pruden A. Critical evaluation of short, long, and hybrid assembly for contextual analysis of antibiotic resistance genes in complex environmental metagenomes. Scientific reports. 2021;11(1):3753. doi: 10.1038/s41598-021-83081-8.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41598-021-83081-8&link_type=DOI) 

51. 51.Gillings MR, Gaze WH, Pruden A, Smalla K, Tiedje JM, Zhu YG. Using the class 1 integron-integrase gene as a proxy for anthropogenic pollution. ISME Journal. 2015;9(6):1269–79. doi: 10.1038/ismej.2014.226. PubMed PMID: 25500508; PubMed Central PMCID: PMCPMC4438328.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ismej.2014.226&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25500508&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 

52. 52.Baker-Austin C, Wright MS, Stepanauskas R, McArthur JV. Co-selection of antibiotic and metal resistance. Trends Microbiol. 2006;14(4):176–82. doi: 10.1016/j.tim.2006.02.006. PubMed PMID: 16537105.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.tim.2006.02.006&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16537105&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000237259100007&link_type=ISI) 

53. 53.Pal C, Bengtsson-Palme J, Rensing C, Kristiansson E, Larsson DGJ. BacMet: antibacterial biocide and metal resistance genes database. Nucleic Acids Research. 2014;42(D1):D737–D43. doi: 10.1093/nar/gkt1252. PubMed PMID: WOS:000331139800109.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkt1252&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24304895&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000331139800109&link_type=ISI) 

54. 54.Bengtsson-Palme J, Larsson DG. Concentrations of antibiotics predicted to select for resistant bacteria: Proposed limits for environmental regulation. Environmental International. 2016;86:140–9. doi: 10.1016/j.envint.2015.10.015. PubMed PMID: 26590482.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.envint.2015.10.015&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26590482&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 

55. 55.Klein EY, Van Boeckel TP, Martinez EM, Pant S, Gandra S, Levin SA, et al. Global increase and geographic convergence in antibiotic consumption between 2000 and 2015. Proc Natl Acad Sci U S A. 2018;115(15):E3463–E70. Epub 2018/03/28. doi: 10.1073/pnas.1717295115. PubMed PMID: 29581252; PubMed Central PMCID: PMCPMC5899442.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTE1LzE1L0UzNDYzIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDYvMDMvMjAyMS4wNi4wMS4yMTI1ODE2NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

56. 56.Joint Programming Initiative on Antimicrobial Resistance. Joint Programming Initiative on Antimicrobial Resistance Mapping Report. 2015.
    
    

57. 57.Morgan DJ, Okeke IN, Laxminarayan R, Perencevich EN, Weisenberg S. Non-prescription antimicrobial use worldwide: a systematic review. Lancet Infectious Diseases. 2011;11(9):692-701. PubMed PMID: WOS:000295239300024.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1473-3099(11)70054-8&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21659004&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000295239300024&link_type=ISI) 

58. 58.Auta A, Hadi MA, Oga E, Adewuyi EO, Abdu-Aguye SN, Adeloye D, et al. Global access to antibiotics without prescription in community pharmacies: A systematic review and meta-analysis. J Infect. 2019;78(1):8–18. Epub 2018/07/10. doi: 10.1016/j.jinf.2018.07.001. PubMed PMID: 29981773.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jinf.2018.07.001&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29981773&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 

59. 59.Patel M, Kumar R, Kishor K, Mlsna T, Pittman CU Jr..,, Mohan D. Pharmaceuticals of Emerging Concern in Aquatic Systems: Chemistry, Occurrence, Effects, and Removal Methods. Chemical reviews. 2019;119(6):3510–673. Epub 2019/03/05. doi: 10.1021/acs.chemrev.8b00299. PubMed PMID: 30830758.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1021/acs.chemrev.8b00299&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30830758&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 

60. 60.Molstad S, Lofmark S, Carlin K, Erntell M, Aspevall O, Blad L, et al. Lessons learnt during 20 years of the Swedish strategic programme against antibiotic resistance. Bull World Health Organ. 2017;95(11):764–73. Epub 2017/11/18. doi: 10.2471/BLT.16.184374. PubMed PMID: 29147057; PubMed Central PMCID: PMCPMC5677604.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2471/BLT.16.184374&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29147057&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 

61. 61.Fernandes N. Illegal sale of antibiotics still common in Hong Kong: Journalism and Media Studies Center - The University of Hong Kong; 2017 [cited 2019 October 7]. Available from: [https://jmsc.hku.hk/2017/03/illegal-sale-of-antibiotics-still-common-in-hong-kong/](https://jmsc.hku.hk/2017/03/illegal-sale-of-antibiotics-still-common-in-hong-kong/).
    
    

62. 62.World Health Organization. Global Antimicrobial Resistance and Use Surveillance System (GLASS) Report: Early Implementation 2020. 2020.Contract No.: 978-92-4-000558-7.
    
    

63. 63.Fahrenfeld N, Bisceglia KJ. Emerging investigators series: sewer surveillance for monitoring antibiotic use and prevalence of antibiotic resistance: urban sewer epidemiology. Environmental Science: Water Research & Technology. 2016;2(5):788–99. doi: 10.1039/c6ew00158k.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1039/c6ew00158k&link_type=DOI) 

64. 64.O’Brien JW, Thai PK, Eaglesham G, Ort C, Scheidegger A, Carter S, et al. A model to estimate the population contributing to the wastewater using samples collected on census day. Environ Sci Technol. 2014;48(1):517–25. Epub 2013/11/29. doi: 10.1021/es403251g. PubMed PMID: 24283359.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1021/es403251g&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24283359&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 

65. 65.Symonds EM, Nguyen KH, Harwood VJ, Breitbart M. Pepper mild mottle virus: A plant pathogen with a greater purpose in (waste)water treatment development and public health management. Water Res. 2018;144:1–12. Epub 2018/07/14. doi: 10.1016/j.watres.2018.06.066. PubMed PMID: 30005176; PubMed Central PMCID: PMCPMC6162155.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.watres.2018.06.066&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30005176&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 

66. 66.Berrazeg M, Diene SM, Medjahed L, Parola P, Drissi M, Raoult D, et al. New Delhi Metallo-beta-lactamase around the world: An eReview using Google Maps. Eurosurveillance. 2014;19(20):2-15. PubMed PMID: WOS:000336499500001.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=WOS:00033649&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 

67. 67.Meyer F, Paarmann D, D’Souza M, Olson R, Glass EM, Kubal M, et al. The metagenomics RAST server - a public resource for the automatic phylogenetic and functional analysis of metagenomes. BMC bioinformatics. 2008;9:386. doi: 10.1186/1471-2105-9-386. PubMed PMID: 18803844; PubMed Central PMCID: PMC2563014.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1471-2105-9-386&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18803844&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 

68. 68.Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biology. 2014;15(12):550. doi: 10.1186/s13059-014-0550-8.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s13059-014-0550-8&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25516281&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 

69. 69.McArthur AG, Waglechner N, Nizam F, Yan A, Azad MA, Baylay AJ, et al. The comprehensive antibiotic resistance database. Antimicrob Agents Chemother. 2013;57(7):3348–57. doi: 10.1128/AAC.00419-13. PubMed PMID: 23650175; PubMed Central PMCID: PMC3697360.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYWFjIjtzOjU6InJlc2lkIjtzOjk6IjU3LzcvMzM0OCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzA2LzAzLzIwMjEuMDYuMDEuMjEyNTgxNjQuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

70. 70.Leplae R, Lima-Mendez G, Toussaint A. ACLAME: A CLAssification of Mobile genetic Elements, update 2010. Nucleic Acids Research. 2010;38:D57–D61. doi: 10.1093/nar/gkp938. PubMed PMID: WOS:000276399100009.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkp938&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19933762&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000276399100009&link_type=ISI) 

71. 71.Clarke KR, Gorley RN, Somerfield PJ, Warwick RM. Change in Marine Communities: An Approach to Statistical Analysis and Interpretation. Plymouth: Primer-E; 2014.
    
    

72. 72.Segata N, Izard J, Waldron L, Gevers D, Miropolsky L, Garrett WS, et al. Metagenomic biomarker discovery and explanation. Genome Biology. 2011;12(6):R60. doi: 10.1186/gb-2011-12-6-r60.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/gb-2011-12-6-r60&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21702898&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F06%2F03%2F2021.06.01.21258164.atom)

 [1]: /embed/graphic-9.gif
 [2]: /embed/graphic-10.gif