Abstract Background The dissemination of antibiotic resistance genes (ARGs) from anthropogenic activities into the environment poses an emerging public health threat. Water constitutes a major vehicle for transport of both biological material and chemical substances. The present study focused on putative antibiotic resistance and integrase genes present in the microbiome of agricultural, urban influenced and protected watersheds in southwestern British Columbia, Canada. A metagenomics approach and high throughput quantitative PCR (HT qPCR) were used to screen for elements of resistance including ARGs and integron-associated integrase genes ( intI ). Sequencing of bacterial genomic DNA was used to characterize the resistome of microbial communities present in watersheds over a one-year period. Results Data mining using CARD and Integrall databases enabled the identification of putative antibiotic resistance genes present in watershed samples. Antibiotic resistance genes presence in samples from various watershed locations was low relative to the microbial population (&lt;1 %). Analysis of the metagenomic sequences detected a total of 78 ARGs and intI1 across all watershed locations. The relative abundance and richness of antibiotic resistance genes was found to be highest in agriculture impacted watersheds compared to protected and urban watersheds. Gene copy numbers (GCNs) from a subset of 21 different elements of antibiotic resistance were further estimated using HT qPCR. Most GCNs of ARGs were found to be variable over time. A downstream transport pattern was observed in the impacted watersheds (urban and agricultural) during dry months. Urban and agriculture impacted sites had a higher GCNs of ARGs compared to protected sites. Similar to other reports, this study found a strong association between intI1 and ARGs (e.g., sul1 ), an association which may be used as a proxy for anthropogenic activities. Chemical analysis of water samples for three major groups of antibiotics was negative. However, the high richness and GCNs of ARGs in impacted sites suggest effects of effluents on microbial communities are occurring even at low concentrations of antimicrobials in the water column. Conclusion Antibiotic resistance and integrase genes in a year-long metagenomic study showed that ARGs were driven mainly by environmental factors from anthropogenized sites in agriculture and urban watersheds. Environmental factors accounted for almost 40% of the variability observed in watershed locations.

Antibiotic resistance genes in agriculture and urban influenced watersheds in southwestern British Columbia