A metagenomic approach and network analysis was used to investigate the wide-spectrum profiles of antibiotic resistance genes (ARGs) and their co-occurrence patterns in 50 samples from 10 typical environments. In total, 260 ARG subtypes belonging to 18 ARG types were detected with an abundance range of 5.4 × 10⁻⁶–2.2 × 10⁻¹ copy of ARG per copy of 16S-rRNA gene. The trend of the total ARG abundances in environments matched well with the levels of anthropogenic impacts on these environments. From the less impacted environments to the seriously impacted environments, the total ARG abundances increased up to three orders of magnitude, that is, from 3.2 × 10⁻³ to 3.1 × 10⁰ copy of ARG per copy of 16S-rRNA gene. The abundant ARGs were associated with aminoglycoside, bacitracin, β-lactam, chloramphenicol, macrolide-lincosamide-streptogramin, quinolone, sulphonamide and tetracycline, in agreement with the antibiotics extensively used in human medicine or veterinary medicine/promoters. The widespread occurrences and abundance variation trend of vancomycin resistance genes in different environments might imply the spread of vancomycin resistance genes because of the selective pressure resulting from vancomycin use. The simultaneous enrichment of 12 ARG types in adult chicken faeces suggests the coselection of multiple ARGs in this production system. Non-metric multidimensional scaling analysis revealed that samples belonging to the same environment generally possessed similar ARG compositions. Based on the co-occurrence pattern revealed by network analysis, tetM and aminoglycoside resistance protein, the hubs of the ARG network, are proposed to be indicators to quantitatively estimate the abundance of 23 other co-occurring ARG subtypes by power functions.

Understanding environmental and biological influences on the dynamics of microbial communities has received great attention in microbial ecology. Here, utilizing large time-series 16S rRNA gene data, we show that in activated sludge of an environmentally important municipal wastewater treatment plant, 5-year temporal dynamics of bacterial community shows no significant seasonal succession, but is consistent with deterministic assemblage by taxonomic relatedness. Biological interactions are dominant drivers in determining the bacterial community assembly, whereas environmental conditions (mainly sludge retention time and inorganic nitrogen) partially explain phylogenetic and quantitative variances and indirectly influence bacterial assembly. We demonstrate a correlation-based statistical method to integrate bacterial association networks with their taxonomic affiliations to predict community-wide co-occurrence and co-exclusion patterns. The results show that although taxonomically closely related bacteria tend to positively co-occur (for example, out of a cooperative relationship), negative co-excluding correlations are deterministically observed between taxonomically less related species, probably implicating roles of competition in determining bacterial assembly. Overall, disclosures of the positive and negative species–species relations will improve our understanding of ecological niches occupied by unknown species and help to predict their biological functions in ecosystems.

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 co-localization 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.