Metagenomic profiles of the resistome in subtropical estuaries: Co-occurrence patterns, indicative genes, and driving factors

Estuaries are resistome hotspots owing to resistome accumulation and propagation at these locations from surrounding rivers, yet the large-scale biogeographic pattern of resistome, especially biocide and metal resistance genes (BMRGs) and its driving mechanisms in estuarine waters remains to be elucidated. Here, a metagenomics-based approach was firstly used to investigate resistome and mobilome profiles in waters from 30 subtropical estuaries, South China. The Pearl River estuaries had a higher diversity and abundance of antibiotic resistance genes (ARGs), BMRGs, and mobile genetic elements (MGEs) when compared with estuaries from east and west regions. Genes resistant to multiple antibiotics, metals, and biocides were the most abundant gene types in the resistome. The abundance of MGEs (e.g., intI1, IS91, and tnpA) was highly associated with the total abundance of resistance genes, suggesting their utility as potential indicators for quantitative estimations of the resistome contamination. Further, MGEs contributed more than bacterial communities in shaping the resistome in subtropical estuaries. Physicochemical factors (e.g., pH) regulated MGE composition and stochastic assembly, which mediated the co-selection of ARGs and BMRGs via horizontal gene transfer. Our findings have important implications and provide a reference on the management of ARGs and BMRGs in subtropical estuarine ecosystems.