Supplementary data are available at Bioinformatics online.

Wastewater treatment plants (WWTPs) are implicated as hotspots for the dissemination of antibacterial resistance into the environment. However, the in situ processes governing removal, persistence, and evolution of resistance genes during wastewater treatment remain poorly understood. Here, we used quantitative metagenomic and metatranscriptomic approaches to achieve a broad-spectrum view of the flow and expression of genes related to antibacterial resistance to over 20 classes of antibiotics, 65 biocides, and 22 metals. All compartments of 12 WWTPs share persistent resistance genes with detectable transcriptional activities that were comparatively higher in the secondary effluent, where mobility genes also show higher relative abundance and expression ratios. The richness and abundance of resistance genes vary greatly across metagenomes from different treatment compartments, and their relative and absolute abundances correlate with bacterial community composition and biomass concentration. No strong drivers of resistome composition could be identified among the chemical stressors analyzed, although the sub-inhibitory concentration (hundreds of ng/L) of macrolide antibiotics in wastewater correlates with macrolide and vancomycin resistance genes. Contig-based analysis shows considerable colocalization between resistance and mobility genes and implies a history of substantial horizontal resistance transfer involving human bacterial pathogens. Based on these findings, we propose future inclusion of mobility incidence (M%) and host pathogenicity of antibiotic resistance genes in their quantitative health risk ranking models with an ultimate goal to assess the biological significance of wastewater resistomes with regard to disease control in humans or domestic livestock.

Antibiotic resistance genes (ARGs) are widespread among bacteria. However, not all ARGs pose serious threats to public health, highlighting the importance of identifying those that are high-risk. Here, we developed an ‘omics-based’ framework to evaluate ARG risk considering human-associated-enrichment, gene mobility, and host pathogenicity. Our framework classifies human-associated, mobile ARGs (3.6% of all ARGs) as the highest risk, which we further differentiate as ‘current threats’ (Rank I; 3%) - already present among pathogens - and ‘future threats’ (Rank II; 0.6%) - novel resistance emerging from non-pathogens. Our framework identified 73 ‘current threat’ ARG families. Of these, 35 were among the 37 high-risk ARGs proposed by the World Health Organization and other literature; the remaining 38 were significantly enriched in hospital plasmids. By evaluating all pathogen genomes released since framework construction, we confirmed that ARGs that recently transferred into pathogens were significantly enriched in Rank II (‘future threats’). Lastly, we applied the framework to gut microbiome genomes from fecal microbiota transplantation donors. We found that although ARGs were widespread (73% of genomes), only 8.9% of genomes contained high-risk ARGs. Our framework provides an easy-to-implement approach to identify current and future antimicrobial resistance threats, with potential clinical applications including reducing risk of microbiome-based interventions.

Antibiotic resistance has emerged as one of the most pressing, global threats to public health. In single-species experiments selection for antibiotic resistance occurs at very low antibiotic concentrations. However, it is unclear how far these findings can be extrapolated to natural environments, where species are embedded within complex communities. We competed isogenic strains of Escherichia coli, differing exclusively in a single chromosomal resistance determinant, in the presence and absence of a pig faecal microbial community across a gradient of antibiotic concentration for two relevant antibiotics: gentamicin and kanamycin. We show that the minimal selective concentration was increased by more than one order of magnitude for both antibiotics when embedded in the community. We identified two general mechanisms were responsible for the increase in minimal selective concentration: an increase in the cost of resistance and a protective effect of the community for the susceptible phenotype. These findings have implications for our understanding of the evolution and selection of antibiotic resistance, and can inform future risk assessment efforts on antibiotic concentrations.

Recent research has demonstrated that selection for antibiotic resistance occurs at very low antibiotic concentrations in single-species experiments, but the relevance of these findings when species are embedded in complex microbial communities is unclear. We show that the strength of selection for naturally occurring resistance alleles in a complex community remains constant from low subinhibitory to above clinically relevant concentrations. Selection increases with antibiotic concentration before reaching a plateau where selection remains constant over a 2-order-magnitude concentration range. This is likely to be due to cross protection of the susceptible bacteria in the community following rapid extracellular antibiotic degradation by the resistant population, shown experimentally through a combination of chemical quantification and bacterial growth experiments. Metagenome and 16S rRNA analyses of sewage-derived bacterial communities evolved under cefotaxime exposure show preferential enrichment for blaCTX-M genes over all other beta-lactamase genes, as well as positive selection and co-selection for antibiotic resistant, opportunistic pathogens. These findings have far-reaching implications for our understanding of the evolution of antibiotic resistance, by challenging the long-standing assumption that selection occurs in a dose-dependent manner.