Estuarine salinity gradient governs sedimentary bacterial community but not antibiotic resistance gene profile

Antibiotic resistance gene (ARG) pollution in estuarine environment has drawn great attention, and it is not clear if the physical and chemical parameters such as salinity, total organic carbon, total nitrogen, total phosphorus and antibiotics affects the distribution of ARGs. Herein, we deciphered the ARG profiles and microbial community compositions in sediments from Jiulong River Estuary (JRE) and Min River Estuary (MRE) of China using high-throughput sequencing-based metagenomics analysis. Furthermore, we explored the influence of salinity on bacterial community and ARG profiles. The results showed that Proteobacteria, Chloroflexi, and Acidobacteria were the dominant phyla in these two estuaries. The abundance of ARGs ranged from 1.05 × 10−1 to 2.93 × 10−1 copy of ARG per copy of 16S rRNA gene in all the sediment samples and the profiles of ARGs presented similar patterns in two estuaries. Multidrug resistance genes were the dominant ARG types in both estuaries, with an overall abundance of 2.39 × 10−2–1.07 × 10−1 copy of ARG per copy of 16S rRNA gene, followed by genes conferring resistance to bacitracin and macrolide-lincosamide-streptogramin. Salinity was an important influencing factor on the bacterial community but not on the ARG profiles. Instead, stochastic processes exerted the main influence on the distribution of ARGs. The comparison of ARG profiles among estuary sediments, marine sediments, and samples from anthropogenic pollution environments revealed remarkable similarity of ARG profiles between samples from estuary sediments and those from municipal wastewater treatment plants. These results suggested that the complex emission of anthropogenic pollution could cause the stochastic ecological pattern of ARGs.