Biomarkers of antibiotic resistance genes during seasonal changes in wastewater treatment systems

To evaluate the seasonal distribution of antibiotic resistance genes (ARGs) and explore the reason for their patterns in different seasons and different systems, two wastewater treatment systems were selected and analyzed using high-throughput qPCR. Linear discriminant analysis (LDA) effect size (LEfSe) was used to discover the differential ARGs (biomarkers) and estimate the biomarkers' effect size. We found that the total absolute abundances of ARGs in inflows and excess sludge samples had no obvious seasonal fluctuations, while those in winter outflow samples decreased in comparison with the inflow samples. Eleven differentially abundant ARGs (biomarker genes, BmGs) (aadA5-02, aac-6-II, cmlA1-01, cmlA1-02, blaOXA10-02, aadA-02, tetX, aadA1, ereA, qacEΔ1-01, and blaTEM) in summer samples and 10 BmGs (tet-32, tetA-02, aacC2, vanC-03, aac-6-I1, tetE, ermB, mefA, tnpA - 07, and sul2) in winter samples were validated. According to 16S rRNA gene sequencing, the relative abundance of bacteria at the phylum level exhibited significant seasonal changes in outflow water (OW), and biomarker bacteria (BmB) were discovered at the family (or genus) level. Synechococcus and vadinCA02 are BmB in summer, and Trichococcus, Lactococcus, Pelosinus, Janthinobacterium, Nitrosomonadaceae and Sterolibacterium are BmB in winter. In addition, BmB have good correlations with BmGs in the same season, which indicates that bacterial community changes drive different distributions of ARGs during seasonal changes and that LEfSe is an acute and effective method for finding significantly different ARGs and bacteria between two or more classes. In conclusion, this study demonstrated the seasonal changes of BmGs and BmB at two wastewater treatment systems.