Metagenomic analysis manifesting intrinsic relatedness between antibiotic resistance genes and sulfate- and iron-reducing microbes in sediment cores of the Pearl River Estuary.

Antibiotic resistance is an increasingly concerned hotspot of human health. Microbial determinants that may affect the profiles of antibiotic resistance genes (ARGs) in the environments need be explored. Here, sediment cores in the Pearl River Estuary (PRE) were analyzed using metagenomic approaches. ARGs were vertically stratified in the PRE sediment cores in terms of both diversity and absolute levels. Multidrug resistance genes could account for approximately 65.0% of the total ARGs, followed by sulfonamides (19.1%), aminoglycosides (5.9%), beta-lactams (4.5%), etc. ARGs related to aminoglycosides, lincosamides, macrolides, sulfonamides and tetracyclines were preferentially enriched in the surface layers of sediment cores. Sulfate-reducing microbes (SRMs) (e.g., Desulfocapsa and Desulfobulbus) and iron-reducing microbes (IRMs) (e.g., Pseudomonas and Sulfurospirillum) were consistently popular and dominant in the PRE sediment cores. The total levels of both SRMs and IRMs were significantly correlated with those of ARGs in the PRE sediment cores (p < 0.01). Network analysis showed that SRM and IRM genera (i.e., Pseudomonas, Shewanella, and Desulfovibrio) had the high co-occurrence with multiple ARG subtypes in the PRE sediment cores such as rsmA, mexK, and mexF. This study highlighted that anaerobic microbes could play significant roles in shaping vertical ARG distribution in the sediments of aquatic environments.