Next Article in Journal
Proposed Countermeasures against Woody Debris Damage Considering Runoff Characteristics
Next Article in Special Issue
Partial Nitrification Characteristics of an Immobilized Carrier in Municipal Wastewater under Low-Temperature Shock: The Role of the Nitrifying Bacterial Community Structure
Previous Article in Journal
Trend Analysis of Streamflows in Relation to Precipitation: A Case Study in Central Italy
Previous Article in Special Issue
Effect of Anammox Granular Sludge Type on the CANON Process with Immobilized Fillers Treating Domestic Wastewater
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Distribution of Antibiotic Resistance Genes and Their Association with Microbes in Wastewater Treatment Plants: A Metagenomics Analysis

1
Key Laboratory of Songliao Aquatic Environment, Ministry of Education, Jilin Jianzhu University, Changchun 130118, China
2
Tianjin Eco-Environmental Monitoring Center, Tianjin 300191, China
3
Yulin Municipal Ecology and Environment Bureau, Yulin 719000, China
*
Authors to whom correspondence should be addressed.
Water 2023, 15(8), 1587; https://doi.org/10.3390/w15081587
Submission received: 2 March 2023 / Revised: 13 April 2023 / Accepted: 15 April 2023 / Published: 19 April 2023
(This article belongs to the Special Issue Functional Microorganisms in Wastewater Treatment)

Abstract

:
The extensive use of antibiotics has resulted in the generation and accumulation of antibiotic resistance genes (ARGs) in the environment, and domestic wastewater is one of the main reservoirs of ARGs and resistant bacteria. In this study, Illumina high-throughput sequencing and network analysis were used to study the microbial community characteristics, ARGs’ occurrence status, and resistance mechanism in the influent and effluent of a domestic sewage treatment plant in Changchun, China. A total of 29 phyla were found in the influent and effluent, including 23 bacterial phyla, 3 archaeal phyla, and 3 eukaryotic phyla. In total, 112 ARG subtypes were detected in the samples, and the dominant ARG subtypes were Erm(35) and tet(W/N/W). In this study, ARGs related to tetracycline and macrolide accounted for a high proportion, and the resistance mechanisms of ARGs detected in the samples were mainly antibiotic inactivation and antibiotic efflux pumps. Co-occurrence maps of ARGs and microbes demonstrated by network analysis indicated that the resistance genes kdpE, GES-5, and tetX may easily bind to microbes, potentially making them antibiotic-resistant bacteria. Fifty-seven bacteria in the genera Cupriavidus, Escherichia, and Collinsell are potential hosts of multiple ARGs. The findings can increase our understanding of the distribution of ARGs and their association with microbes in wastewater treatment plants, and also provide a research foundation for controlling the diffusion of ARGs in the environment.

1. Introduction

Antibiotics are widely used to treat and control bacterial infections and have revolutionized medicine and saved countless lives. However, the abuse of antibiotics in agriculture, animal husbandry, and medical care has led to the presence of antibiotic residues in the environment, humans, animals, and plants [1]. Antibiotics released into the environment can cause microorganisms to develop antibiotic resistance genes (ARGs). ARGs may replicate and spread with microorganisms in the environment through horizontal gene transfer (HGT) and vertical gene transfer (VGT), which causes even more microorganisms to develop resistance. Clearly, this poses a serious threat to human health and ecological security, causing ARGs to be recognized as emerging environmental pollutants [1,2,3,4].
Environmental sample monitoring results have shown that antibiotics and ARGs are present in various environmental media, such as natural water [5], sediment, soil [6], and domestic sewage. Wastewater treatment plants (WWTPs) are the main places to collect and treat municipal domestic sewage. The biological treatment process can remove most pollutants, such as N, P, and organic pollutants, in wastewater [7]. However, the growth and reproduction of microorganisms during biological treatment provide convenient conditions for the spread of ARGs, resulting in the accumulation of a large number of antibiotics and ARGs in WWTPs [8,9]. Therefore, municipal wastewater treatment plants are important repositories of ARGs and sites for ARG propagation and diffusion [10,11]. The concentration, composition, and distribution characteristics of antibiotics and ARGs in WWTPs have become popular research topics. From an epidemiological perspective, the monitoring and analyzing of ARGs in WWTPs also are more beneficial to overcoming traditional surveillance limitations due to sample data only from hospitals [12]. Lee et al. [13] quantitatively detected the ARG concentration in the UV-disinfection effluent of two sewage treatment plants and found that up to 5.4 × 1016–4.2 × 1018 copies/L of ARGs were discharged into the environment every day. Rodríguez et al. [14] used the shotgun metagenome and culture-based approaches to determine bacterial communities, the resistome and mobilome in a WWTP in Colombia, the results showed that the genes encoding resistance to macrolide-lincosamide-streptogramin, β-lactams, and those conferring multidrug resistance were the most abundant in the samples. In addition, tetracycline, sulfonamide, and quinolone ARGs are also highly detected in the WWTPs [1,15]. The research results above show that there are abundant and diverse ARGs with high concentrations in sewage treatment plants. In general, ARGs decrease significantly after sewage treatment. However, residual ARGs in sewage may still cause ecological risks. A study on the characteristics and risk assessment of ARGs in WWTPs in Hong Kong found that the ARGs’ removal efficiency in WWTPs could reach a 2.34–2.43 log reduction rate. The investigated WWTPs had relatively high efficiency in eliminating ARGs compared with the global average. Despite this, ARG variants with high risk were still highly abundant in effluents, which deserves our attention [16]. Therefore, it is essential to study the distribution characteristics of ARGs during sewage treatment to analyze the horizontal transfer of ARGs and perform an accurate ecological risk assessment.
The species of ARGs have a strong co-occurrence relationship with the microbial community structure in the environment, in addition to extensive diffusion by microbial reproduction during pollutant removal. Zhao et al. [17] studied the correlation between bacteria and virus groups and ARGs subtypes, and found that in addition to bacteria, viruses, especially phages, also contributed to the spread of ARGs. They also found that rpoB, drfE, gyrA, and parC were significantly correlated with bacteria and phages. According to correlation analysis of microorganisms and ARGs in the activated sludge of sewage treatment plants, some environmental bacteria (such as Clostridium and Nitrosomonas) may be potential hosts for many ARGs [11]. Similarly, Aeromonadaceae and Enterobacteriaceae have also been found to have a significant association with some ARGs as blaOXA-58 and blaKPC genes [14]. However, there is less research on the correlation analysis of the ARGs’ subtype and microbes in terms of the resistance mechanism. Based on the resistance mechanism, the ARG subtype and microbial correlation can be analyzed to determine which resistance mechanism is common in each ARG subtype. This can be combined with the microbes to provide theoretical guidance for solving ARGs and antibiotic resistance bacteria (ARB) problems in different environments.
Previous studies have shown that PCR and metagenomics are effective tools for studying the abundance of ARGs [18]. However, PCR relies on available primers of known ARGs and is not suitable for the discovery of unknown ARGs in samples [19]. Metagenomics sequencing can overcome the shortcomings of PCR-based techniques [19,20] and provide comprehensive information about the diversity and abundance of ARGs. Although there have been some studies on ARGs in WWTPs using metagenomic sequencing, due to the differences in characteristics, composition, sources, and sewage treatment process in different regions, more extensive research is needed to obtain more comprehensive data. The wastewater plant in this study is in Changchun, China, which is located in the severe cold region of China with a permanent population of more than 9 million. Changchun is an important industrial base in northeast China, and there are many pharmaceutical enterprises. Therefore, it is necessary to explore the occurrence characteristics of ARGs in wastewater in this region. In this study, the influent and effluent water samples of a domestic sewage treatment plant in Changchun, China, were collected, and the microbial community composition and ARG characteristics of the water samples were studied by Illumina NovaSeq and macrogenetic sequencing. The distribution characteristics of microbial communities and ARGs in water samples were analyzed, and the correlation between ARG subtypes and microbial was analyzed based on the resistance mechanism. The results provided a basis for further research on the environmental behavior of ARGs in WWTPs.

2. Materials and Methods

2.1. Sampling and DNA Extraction

The sewage plant in this study is located in Changchun, China. It is a conventional third-level treatment sewage plant, with an improved A2/O + ozonation disinfection as the primary process. The designed daily treatment capacity was 100,000 m3/d, and the primary sewage source was urban domestic sewage. Wastewater samples were collected at the influent and effluent of the sewage plant. The water samples were filtered by a 0.45 μm filter membrane, which was placed in a sterile centrifuge tube for measurements.
Genomic DNA extraction and metagenomics sequencing were completed by Shanghai Personalbio Technology Co., Ltd. DNA samples were extracted using an OMEGA Soil DNA Kit (D5625-01). A Qubit™ 4 Fluorimeter (Qubit™ 4 Fluorimeter with Wi-Fi: Q33238; Qubit™ assay tubes: Q32856; Qubit™ 1X dsDNA HS Assay Kit: Q33231) was used to assay the DNA concentration, and mass was determined by 1% agarose gel electrophoresis. The level of DNA concentration was about 100 ng/μL.

2.2. Metagenomics Sequencing and Data Analysis

Illumina NovaSeq whole genome shotgun DNA library construction and sequencing were carried out as follows. After DNA fragmentation, inserted fragment libraries of approximately 400 base pairs were constructed, one library per sample. These libraries were double-ended sequenced (2 × 150 bp). Quality control was performed on the obtained original sequence data, and sequence lengths < 50 bp and sequences containing fuzzy bases were removed. Kranken2 [21] was used for species annotation to obtain the species abundance table of each sample, and then the microbial diversity was analyzed. The gene sequences were annotated in the comprehensive antibiotic resistance database (CARD) database to obtain the abundance of ARGs in each sample. The ARGs whose ratio of protein length predicted by the sample gene sequence to that of the CARD database was 100 were screened. TPM (transcripts per kilobase per million mapped reads) was used as the standard to measure the gene abundance of each sample (Equations (1) and (2)). To visualize the correlations between major bacteria, archaea, and viruses and ARGs, Pearson correlations between microbial genera and ARGs with the top 60 relative abundances were calculated to construct a correlation matrix:
T P M = r g × r l × 10 6 f l g × T
T = r g × r l f l g
where rg is the number of reads aligned to the gene sequence, rl is the length of the read, and flg is the length of the gene.

3. Results and Discussion

3.1. Diversity and Richness of Microbial Communities

After Illumina high-throughput sequencing was used to assemble and filter the data, the microbial community diversity of different samples was analyzed and compared. As shown in Figure 1, in both the influent and effluent samples, the sparse curves were close to the saturation trend, indicating that the sequencing depth in this study was sufficient for classification coverage (Figure 1A). At the species level, the alpha diversity of microbes in the influent and effluent was analyzed, including the ACE richness indices and Shannon diversity indices. The ACE showed that microbial species richness was lower in the effluent than influent (Figure 1B). Shannon analysis showed that the diversity of microbial species in the effluent was lower than that in the influent (Figure 1C).
To better understand the characteristics of microbial flora structures in the influent and effluent, Kranken2 was used to annotate the effective sequences of each sample from the phylum to species level (Table 1). Bacteria were the dominant microbes in both the influent and effluent, accounting for 99.12–99.61% in the influent and 97.45–99.36% in the effluent. Yasir et al. [22] studied the microbial community in a biological sewage treatment plant in Jeddah and also found that bacteria were the main microbial type in the influent and effluent (accounting for 98.9% of influent and 98.1% of effluent). Giwa et al. [23] analyzed the microbial community structure in the influent, activated, return sludge and effluent of the sewage plant, and found that bacteria were also dominant microbes in the influent. The analysis on the microbial community structure of the whole wastewater plant found that except bacteria, eukaryota had a high abundance, followed by archaea and viruses. In this study, except bacteria, archaea had a high abundance, followed by viruses and eukaryota, which was slightly different from the results of Giwa et al.
A total of 23 bacterial phyla, 3 archaeal phyla, 3 eukaryotic phyla, and 1 viral order were found in influent and effluent samples. At the species level, 2723 bacterial species were found in the influent samples, while only 971 bacterial species were found in the effluent samples. In the influent samples, 1 viral order, 4 viral families, 25 viral genera, and 15 viral species were identified. In the effluent samples, the number of taxonomic units for viruses was significantly reduced. The number of taxonomic units for archaea and eukaryota showed little difference between influent and effluent samples (Table 2).
Bacteria are important vectors of ARGs, and ARGs may spread, transfer, and mutate through the reproduction of bacterial hosts [24]. At phylum level, Proteobacteria were the most abundant bacteria in this study, followed by Bacteroidetes, Firmicutes, and Actinobacteria (Figure 2A). Proteobacteria accounted for 36.78% in the influent but increased to 52.85% in the effluent, in which it was dominant. Compared with the influent, the relative abundance of Firmicutes and Bacteroidetes decreased from 24.29% to 12.69% and from 34.34% to 17.24%, respectively. Proteobacteria, Bacteroidetes, and Actinomyces are considered to be common hosts of ARGs in the activated sludge of sewage treatment plants [25]. At the family level (Figure 2B), Bacteroidaceae (19.15%), Campylobacteraceae (14.91%), and Prevotellaceae (11.33%) were the dominant families in the influent, followed by Lachnospiraceae (9.04%) and Ruminococcaceae (7.23%). In the effluent samples, the dominant bacteria were Enterobacteriaceae (18.81%), Bacteroidaceae (11.27%), and Nitrospiraceae (9.61%). A previous study on Enterobacteriaceae in municipal sewage found that this type of bacteria could carry a variety of ARGs [26]. At the genus level, the top 6 genera were Bacteroides (accounting for 19.49% in the influent and 11.31% in the effluent), Prevotella (accounting for 11.23% in the influent and 3.53% in the effluent), Arcobacter (accounting for 11.22% in the influent and 3.99% in the effluent), Acinetobacter (accounting for 6.30% in the influent and 4.57% in the effluent), Aliiarcobacter (accounting for 3.63% in the influent and 1.34% in the effluent), and Aeromonas (accounting for 3.58% in the influent and 3.15% in the effluent) (Figure 2C).
Previous studies have shown that Bacteroides, Acinetobacter, Aeromonas, and Pseudomonas are all potential pathogens that may carry different types of multi-resistant ARGs [1,27]. This study showed that the influent of the WWTPs contained some potential host bacteria of ARGs. After treatment, the abundance of host bacteria such as Bacteroides in the sewage was reduced, thus reducing the risk of the spread of drug resistance. It is worth noting that ozonation was used in the wastewater plants sampled in this study, but in the effluent, some potential host bacteria of ARGs were still present. This may be due to the fact that ozonation disinfection often does not effectively inactivate bacteria in wastewater, and that cells in agglomerates, sludge flocs or particles of suspension may be shielded and protected against the action of a disinfectant [28]. Therefore, how to choose a more effective disinfection process may be the focus of future research on ARG removal. Previous studies have shown that incomplete mineralization by-products may be formed during ozonation of treated wastewater that is toxic to aquatic organisms. The toxicity class of treated wastewater may change from the completely non-toxic to very high hazard category [29]. Therefore, when selecting the disinfection process, in addition to considering the inactivation effect, the ecological safety of the disinfection method also needs our close attention.
The characteristics of archaeal communities in influent and effluent samples were also studied. The results showed that Euryarchaeota accounted for 99.6% and 98.73% of all Archaea in the influent and the effluent, respectively, and was the dominant archaeal phylum, followed by Thaumarchaeota and Crenarchaeota (Figure S1). Li et al. [30] found that the dominant archaea in the receiving water of a domestic sewage discharge were Euryarchaeota, Thaumarchaeota, and Crenarchaeota, successively, which was consistent with the results of this study. Euryarchaeota can appear in various ecological niches and contain a variety of methanogenic archaea [31] that play an indispensable role in anaerobic wastewater treatment [32]. Thaumarchaeota contain most of the ammoxidation archaea, which participate in ammoxidation. Euryarchaeota and Thaumarchaeota play an important role in the nitrogen cycle of sewage treatment systems [33]. The dominant archaea at the genus level were Methanothrix (accounting for 45.24% of the influent and 12.06% of the effluent), Methanosarcina (accounting for 26.44% of the influent and 79.87% of the effluent), and Methanobrevibacter (accounting for 12.87% of the influent and 2.41% of the effluent) (Figure 3). Compared with the influent, except for Methanosarcina, Methanothrix, Methanobrevibacter, and Methanomethylovorans were significantly lower in the effluent. A previous study showed that Methanothrix and Methanosarcina are dominant methanogens in the anaerobic digestion [34]. In this study, Methanothrix had higher abundance than Methanosarcina in the influent. However, in the effluent, this genus was completely outnumbered by the metabolically versatile methanogens Methanosarcinaceae. The change of this dominant genus may be related to the concentrations of volatile fatty acids (VAFs) [35]. Both methanogens require VFAs in the process of anaerobic digestion. Based on K/r hypothesis, there might be a competitive relationship between them. High concentrations of VFAs may be more conducive to reproducing Methanosarcina, while low concentrations of VFAs may enable Methanothrix to remain the dominant genus. Perhaps the change in these substrate concentrations caused the change in the dominant genus in the present study.
Next, the main types of viruses in sewage were studied. In the influent, at the viral family level, Myoviridae, Siphoviridae, and Podoviridae were dominant. At the genus level, Kp34virus (15.01%) and Kp15virus (10.14%) were dominant, followed by Kp36virus (7.67%), Pamx74virus (7.61%), and Hk578virus (6.85%), which belong to the order Caudovirales. Zhao et al. investigated the influent and effluent microbial communities in six decentralized facilities receiving livestock and household wastewater from rural areas of Changzhou City, Jiangsu Province, China. Consistent with this study, they found that the dominant viral structures in the influent and effluent were also Myoviridae, Siphoviridae, and Podoviridae [17]. Some studies have found that Caudovirales play a dominant role in virus-related sequences among sewage microbes [36], and Caudovirales are prone to horizontal gene transfer and may affect the ecological environment of sewage microorganisms through interactions with bacterial hosts [37]. The free ARGs in sewage may bind to Caudovirales through horizontal gene transfer, and viruses carrying AGRs may invade bacteria and eventually turn bacteria into ARBs.

3.2. Characteristics of ARGs

Based on the CARD database, this study identified 109 subtypes of ARGs in the influent samples and 97 subtypes of ARGs in the effluent samples after screening (the percentage length of reference sequence was 100%), with a total of 112 subtypes (Figure 4). The relative abundance of 26 ARG subtypes was >1%. The 26 ARG subtypes mainly included 7 macrolide resistance genes (Erm(35), mphE, mel, msrE, EreB, EreA2, and Mrx); 6 tetracycline resistance genes (tet(W/N/W), tet(39), tet(C), tet(40), tetX, and acrB); and 4 cephalosporin resistance genes (AAC(6′)-Ie-APH(2″)-Ia, OXA-333, acrB, and GES-5). The ARG abundance calculation based on TPM showed that the total abundance of ARGs in the influent was 305 TPM, while that in the effluent was only 15 TPM, indicating that the wastewater treatment plant (A2/O + ozonation) in this study had a good removal effect on ARGs. Hultman et al. [38] collected influent and effluent samples from two different municipal sewage treatment plants in Finland and found that the ARG content in the effluent decreased significantly after treatment. Pazda et al. [39] compared and analyzed the removal efficiency of quinolone, macrolide, sulfonamide resistance genes, and other ARGs by different sewage treatment processes and found that the number of ARG types in the treated effluent decreased significantly. Despite this significant decrease after treatment, Tang et al. [40] found that the emission of these low-dose “quasi-persistent” ARGs still posed ecological risks to the receiving water bodies.
In this study, we observed a high abundance of ARG subtypes, including Erm (35) (macrolide/lincosamide/streptogramin antibiotic), tet (W/N/W) (tetracycline antibiotic), mphE (macrolide antibiotic), and mel (macrolide/streptogramin antibiotic). As shown in Figure 5B, ARGs associated with tetracycline and macrolides accounted for a higher proportion in the sample. In the influent, the abundance of ARGs associated with tetracycline was 17.77% and the abundance of ARGs associated with macrolides was 12.66%. In the effluent, the abundance of two types of ARGs was 12.43% and 19.43%, respectively. This may be related to two factors: tetracycline and macrolide are still the most commonly used antibiotics in China [41]; these two antibiotics are difficult to metabolize in vivo, so most antibiotics entering biological bodies end up in wastewater as unchanged parent compounds [42]. As a result, the content of tetracycline and macrolides in urban domestic sewage is high, which has led to the abundance of related ARGs. Of further concern in this study is that macrolide antibiotic abundance related to ARGs was higher in the antibiotic effluent compared with antibiotic influent, suggesting that macrolide ARGs were more difficult to remove than other types.
Previous studies have shown that the main mechanisms of microbial resistance to antibiotics include antibiotic inactivation, antibiotic target protection, and antibiotic efflux pumps [1,39]. The 112 subtypes of ARGs found in this study included 7 drug resistance mechanisms. In the influent, the highest abundance of ARGs with antibiotic inactivation mechanism was 118 TPM, followed by target protection mechanism, efflux pump mechanism, target alteration mechanism, target replacement mechanism, ATA-AEF, and ATA-AEF-RPA. In the effluent, ARGs with an antibiotic inactivation mechanism also had the highest abundance (7 TPM) (Figure 5A). The results showed that ARGs present in the wastewater treatment plant mainly caused microbial resistance via an antibiotic inactivation mechanism. Microbes often develop resistance to aminoglycosides, β-lactams, and macrolides through an antibiotic inactivation mechanism [40]. Antibiotic target protection, another antibiotic resistance mechanism in this study, is often related to microbial resistance to erythromycin antibiotics [43]. Antibiotic efflux has long been considered to be the main mechanism responsible for the development of simultaneous resistances to multiple antibiotics (e.g., tetracyclines, macrolides, and fluoroquinolones) [44].

3.3. Co-Occurrence of ARGs and Microbes

To further investigate the relationships between microbes, ARGs, and resistance mechanisms using correlation analysis, we used Gephi to visualize the co-occurrence patterns between ARG subtypes and microbial communities. ARGs co-occurring with bacteria (top 80 genera in influent samples), archaea, and viruses (top 40 genera in influent samples) were analyzed based on network analysis by Pearson correlation analysis. The results in Figure 6A show that the network of bacteria and ARGs consisted of 100 nodes (62 bacterial genera and 38 ARGs subtypes) and 109 edges. A total of 57 out of 62 bacterial genera may carry more than 2 ARG subtypes. Most of the 38 ARG subtypes were associated with antibiotic inactivation and antibiotic efflux mechanisms. Among these microbes, 14 harbored kdpE (antibiotic efflux), 12 harbored GES-5 (antibiotic inactivation), and 11 harbored tetX (antibiotic efflux). Although the abundances of these three ARG subtypes were not high, many microbial species carried these three ARG subtypes, indicating that microorganisms may have been more likely to serve as hosts of these three ARG subtypes. The 62 bacterial genera were mainly classified into Proteobacteria and Firmicutes phyla, among which Anaerostipes, Fusicatenibacter, and Cupriavidus carried more diverse ARG isoforms than the other genera (5, 4, and 4 species, respectively). In Figure 6B, the network of archaea with viruses and ARGs contained 71 nodes (13 Archaeal genera, 22 viral genera, and 36 ARG subtypes) and 62 edges. Among the 36 ARG subtypes, the antibiotic inactivation and antibiotic efflux mechanisms were still dominant. Methanosarcina carried four species (msrE, tet(39), QnrS2, and cphA5), and three resistance mechanisms (antibiotic inactivation, antibiotic efflux, and antibiotic target protection) (Figure 6B). Analysis of the microbial community structure showed that the abundance of Methanosarcina was high in both influent and effluent, and its abundance in the effluent was higher than that in the influent, indicating that the removal effect of Methanosarcina by A2/O + ozonation process in this study was lower than in the case of other archaeal genera.
Therefore, Methanosarcina needs to be focused on as a potential ARG host. Bacteriophages, the most abundant organisms on Earth, play an important role in the composition and dynamics of microbial communities, biogeochemical cycles, and bacterial evolution in ecosystems around the world. Bacteriophages and their mediated transduction may play an important role in the enrichment and diffusion of ARGs. In this study, three kinds of phages (sp58virus, G7cvirus, and Np1virus) carried three kinds of ARGs. Among them, sp58virus may have carried GlpT, dfrF, and ANT(6)-Ia; G7cvirus may have carried VEB-9, OXA-1, and CMY-1; Np1virus may have carried tet(C), OXA-347, and dfrB4.

4. Conclusions

In this study, metagenomics sequencing and network analysis were used to explore the community structure characteristics of bacteria, archaea, and viruses in the influent and effluent of urban wastewater treatment plants and to determine correlations between microbial communities and ARGs. The results showed that bacteria were the dominant microorganisms in both the influent and effluent. Proteobacteria, Bacteroidetes, Firmicutes, and Actinobacteria were the most abundant bacterial phyla in samples. ARGs in WWTPS conferred microbial resistance mainly through an antibiotic inactivation mechanism. Among bacteria, Anaerostipes, Fusicatenibacter, and Cupriavidus carried more diverse ARG isoforms than other genera. Among archaea, Methanosarcina was the dominant genus in the effluent and possibly carried four ARGs, which requires special attention in future research. Among viruses, Siphoviridae and Myoviridae showed significant positive correlations with ARGs. Ozonation was used in the sewage plants sampled in this study, but there were still some potential host microbes of ARGs in the effluent, and the potential ecological risks caused by residual ARGs cannot be ignored. Therefore, how to choose a more effective disinfection process may be the focus of future research on ARGs’ removal. In the future, it is necessary to conduct long-term monitoring of ARGs in WWTPs, and to carry out tracking monitoring ARGs in receiving water to assess the risks caused by ARGs to human health at a long-term and large-scale level.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/w15081587/s1, Figure S1: The top 3 abundant archaeal phyla in influent and effluent, visualized by Origin 2022.

Author Contributions

Conceptualization, M.S. and C.L. (Chenyang Li); methodology, M.L., Y.H. and C.L. (Chen Lyu); validation, X.B., T.C. and H.C.; writing—original draft preparation, X.H.; writing—review and editing, J.L. and M.S.; supervision, J.L. and K.F.; funding acquisition, M.S. and C.L. (Chenyang Li). All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by Science and Technology Development Plane of Jilin Province (Grant no. 20210203177SF), Education Department of Jilin Province (Grant No. JJKH 20230350KJ), and the Bureau of Science and Technology of Yulin.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviation

AbbreviationFull Name
MLSmacrolide antibiotic; lincosamide antibiotic; streptogramin antibiotic
MSmacrolide antibiotic; streptogramin antibiotic
CPcephalosporin; penam
MCmonobactam; cephalosporin
FCGPTRPTfluoroquinolone antibiotic; cephalosporin; glycylcycline; penam; tetracycline antibiotic; rifamycin antibiotic; phenicol antibiotic; triclosan
SSsulfonamide antibiotic; sulfone antibiotic
GTglycylcycline; tetracycline antibiotic
CCPcarbapenem; cephalosporin; penam

References

  1. Manoharan, R.K.; Ishaque, F.; Ahn, Y.H. Fate of antibiotic resistant genes in wastewater environments and treatment strategies—A review. Chemosphere 2022, 298, 134671. [Google Scholar] [CrossRef]
  2. van der Donk, C.F.M.; van de Bovenkamp, J.H.B.; De Brauwer, E.I.G.B.; De Mol, P.; Feldhoff, K.-H.; Kalka-Moll, W.M.; Nys, S.; Thoelen, I.; Trienekens, T.A.M.; Stobberingh, E.E. Antimicrobial Resistance and Spread of Multi Drug Resistant Escherichia coli Isolates Collected from Nine Urology Services in the Euregion Meuse-Rhine. PLoS ONE 2012, 7, e47707. [Google Scholar] [CrossRef] [PubMed]
  3. Taylor, N.G.H.; Verner-Jeffreys, D.W.; Baker-Austin, C. Aquatic systems: Maintaining, mixing and mobilising antimicrobial resistance? Trends Ecol. Evol. 2011, 26, 278–284. [Google Scholar] [CrossRef]
  4. Aminov, R.I. The role of antibiotics and antibiotic resistance in nature. Environ. Microbiol. 2009, 11, 2970–2988. [Google Scholar] [CrossRef] [PubMed]
  5. Blasco, M.D.; Esteve, C.; Alcaide, E. Multiresistant waterborne pathogens isolated from water reservoirs and cooling systems. J. Appl. Microbiol. 2008, 105, 469–475. [Google Scholar] [CrossRef]
  6. Forsberg, K.J.; Reyes, A.; Wang, B.; Selleck, E.M.; Sommer, M.O.A.; Dantas, G. The Shared Antibiotic Resistome of Soil Bacteria and Human Pathogens. Science 2012, 337, 1107–1111. [Google Scholar] [CrossRef] [PubMed]
  7. Zhao, W.; Bi, X.; Peng, Y.; Bai, M. Research advances of the phosphorus-accumulating organisms of Candidatus Accumulibacter, Dechloromonas and Tetrasphaera: Metabolic mechanisms, applications and influencing factors. Chemosphere 2022, 307, 135675. [Google Scholar] [CrossRef]
  8. Łuczkiewicz, A.; Jankowska, K.; Fudala-Książek, S.; Olańczuk-Neyman, K. Antimicrobial resistance of fecal indicators in municipal wastewater treatment plant. Water Res. 2010, 44, 5089–5097. [Google Scholar] [CrossRef] [PubMed]
  9. Novo, A.; Manaia, C.M. Factors influencing antibiotic resistance burden in municipal wastewater treatment plants. Appl. Microbiol. Biot. 2010, 87, 1157–1166. [Google Scholar] [CrossRef]
  10. Ju, F.; Beck, K.; Yin, X.; Maccagnan, A.; McArdell, C.S.; Singer, H.P.; Johnson, D.R.; Zhang, T.; Bürgmann, H. Wastewater treatment plant resistomes are shaped by bacterial composition, genetic exchange, and upregulated expression in the effluent microbiomes. ISME J. 2019, 13, 346–360. [Google Scholar] [CrossRef]
  11. Guo, J.; Li, J.; Chen, H.; Bond, P.L.; Yuan, Z. Metagenomic analysis reveals wastewater treatment plants as hotspots of antibiotic resistance genes and mobile genetic elements. Water Res. 2017, 123, 468–478. [Google Scholar] [CrossRef] [PubMed]
  12. Sims, N.; Kasprzyk-Hordern, B. Future perspectives of wastewater-based epidemiology: Monitoring infectious disease spread and resistance to the community level. Environ. Int. 2020, 139, 105689. [Google Scholar] [CrossRef] [PubMed]
  13. Lee, J.; Jeon, J.H.; Shin, J.; Jang, H.M.; Kim, S.; Song, M.S.; Kim, Y.M. Quantitative and qualitative changes in antibiotic resistance genes after passing through treatment processes in municipal wastewater treatment plants. Sci. Total Environ. 2017, 605–606, 906–914. [Google Scholar] [CrossRef] [PubMed]
  14. Rodríguez, E.A.; Ramirez, D.; Bálcazar J., L.; Jiménez J., N. Metagenomic analysis of urban wastewater resistome and mobilome: A support for antimicrobial resistance surveillance in an endemic country. Environ. Pollut. 2021, 276, 116736. [Google Scholar] [CrossRef]
  15. Mao, D.; Yu, S.; Rysz, M.; Luo, Y.; Yang, F.; Li, F.; Hou, J.; Mu, Q.; Alvarez, P.J.J. Prevalence and proliferation of antibiotic resistance genes in two municipal wastewater treatment plants. Water Res. 2015, 85, 458–466. [Google Scholar] [CrossRef]
  16. Yin, X.; Yang, Y.; Deng, Y.; Huang, Y.; Li, L.; Chan, L.Y.L.; Zhang, T. An assessment of resistome and mobilome in wastewater treatment plants through temporal and spatial metagenomic analysis. Water Res. 2022, 209, 117885. [Google Scholar] [CrossRef] [PubMed]
  17. Zhao, J.; Li, B.; Lv, P.; Hou, J.; Qiu, Y.; Huang, X. Distribution of antibiotic resistance genes and their association with bacteria and viruses in decentralized sewage treatment facilities. Front. Environ. Sci. Eng. 2021, 16, 35. [Google Scholar] [CrossRef]
  18. Su, J.-Q.; Wei, B.; Ou-Yang, W.-Y.; Huang, F.-Y.; Zhao, Y.; Xu, H.-J.; Zhu, Y.-G. Antibiotic Resistome and Its Association with Bacterial Communities during Sewage Sludge Composting. Environ. Sci. Technol. 2015, 49, 7356–7363. [Google Scholar] [CrossRef]
  19. Schmieder, R.; Edwards, R. Insights into antibiotic resistance through metagenomic approaches. Future Microbiol. 2011, 7, 73–89. [Google Scholar] [CrossRef]
  20. Yang, Y.; Li, B.; Zou, S.; Fang, H.H.P.; Zhang, T. Fate of antibiotic resistance genes in sewage treatment plant revealed by metagenomic approach. Water Res. 2014, 62, 97–106. [Google Scholar] [CrossRef]
  21. Wood, D.E.; Lu, J.; Langmead, B. Improved metagenomic analysis with Kraken 2. Genome Biol. 2019, 20, 257. [Google Scholar] [CrossRef] [PubMed]
  22. Yasir, M. Analysis of Microbial Communities and Pathogen Detection in Domestic Sewage Using Metagenomic Sequencing. Diversity 2020, 13, 6. [Google Scholar] [CrossRef]
  23. Giwa, A.S.; Ali, N.; Athar, M.A.; Wang, K. Dissecting microbial community structure in sewage treatment plant for pathogens’ detection using metagenomic sequencing technology. Arch. Microbiol. 2020, 202, 825–833. [Google Scholar] [CrossRef]
  24. Sørensen, S.J.; Bailey, M.; Hansen, L.H.; Kroer, N.; Wuertz, S. Studying plasmid horizontal transfer in situ: A critical review. Nat. Rev. Microbiol. 2005, 3, 700–710. [Google Scholar] [CrossRef] [PubMed]
  25. Zhang, G.; Guan, Y.; Zhao, R.; Feng, J.; Huang, J.; Ma, L.; Li, B. Metagenomic and network analyses decipher profiles and co-occurrence patterns of antibiotic resistome and bacterial taxa in the reclaimed wastewater distribution system. J. Hazard. Mater. 2020, 400, 123170. [Google Scholar] [CrossRef]
  26. Ferreira da Silva, M.; Vaz-Moreira, I.; Gonzalez-Pajuelo, M.; Nunes, O.C.; Manaia, C.M. Antimicrobial resistance patterns in Enterobacteriaceae isolated from an urban wastewater treatment plant. FEMS Microbiol. Ecol. 2007, 60, 166–176. [Google Scholar] [CrossRef]
  27. Tong, J.; Tang, A.; Wang, H.; Liu, X.; Huang, Z.; Wang, Z.; Zhang, J.; Wei, Y.; Su, Y.; Zhang, Y. Microbial community evolution and fate of antibiotic resistance genes along six different full-scale municipal wastewater treatment processes. Bioresource Technol. 2019, 272, 489–500. [Google Scholar] [CrossRef]
  28. Affek, K.; Muszyński, A.; Załęska-Radziwiłł, M.; Doskocz, N.; Ziętkowska, A.; Widomski, M. Evaluation of ecotoxicity and inactivation of bacteria during ozonation of treated wastewater. Desalin. Water Treat. 2020, 192, 176–184. [Google Scholar] [CrossRef]
  29. Affek, K.; Muszyński, A.; Doskocz, N.; Załęska-Radziwiłł, M. Ecotoxicological effects of disinfection of treated wastewater. Desalin. Water Treat. 2021, 233, 190–198. [Google Scholar] [CrossRef]
  30. Li, D.; Jiang, X.; Wang, J.; Wang, K.; Zheng, B. Effect of Sewage and Industrial Effluents on Bacterial and Archaeal Communities of Creek Sediments in the Taihu Basin. Water 2017, 9, 373. [Google Scholar] [CrossRef]
  31. Xu, R.; Zhang, Y.; Xiong, W.; Sun, W.; Fan, Q.; Zhaohui, Y. Metagenomic approach reveals the fate of antibiotic resistance genes in a temperature-raising anaerobic digester treating municipal sewage sludge. J. Clean. Prod. 2020, 277, 123504. [Google Scholar] [CrossRef]
  32. Tabatabaei, M.; Rahim, R.A.; Abdullah, N.; Wright, A.-D.G.; Shirai, Y.; Sakai, K.; Sulaiman, A.; Hassan, M.A. Importance of the methanogenic archaea populations in anaerobic wastewater treatments. Process Biochem. 2010, 45, 1214–1225. [Google Scholar] [CrossRef]
  33. Robertson, C.E.; Harris, J.K.; Spear, J.R.; Pace, N.R. Phylogenetic diversity and ecology of environmental Archaea. Curr. Opin. Microbiol. 2005, 8, 638–642. [Google Scholar] [CrossRef] [PubMed]
  34. Wojcieszak, M.; Pyzik, A.; Poszytek, K.; Krawczyk, P.S.; Sobczak, A.; Lipinski, L.; Roubinek, O.; Palige, J.; Sklodowska, A.; Drewniak, L. Adaptation of methanogenic inocula to anaerobic digestion of maize silage. Front. Microbiol. 2017, 8, 1881. [Google Scholar] [CrossRef]
  35. Garlicka, A.; Umiejewska, K.; Nielsen, P.H.; Muszynski, A. Hydrodynamic disintegration of thickened excess sludge and maize silage to intensify methane production: Energy effect and impact on microbial communities. Bioresource Technol. 2023, 376, 128829. [Google Scholar] [CrossRef] [PubMed]
  36. Gulino, K.; Rahman, J.; Badri, M.; Morton, J.; Bonneau, R.; Ghedin, E. Initial Mapping of the New York City Wastewater Virome. mSystems 2020, 5, e00876-00819. [Google Scholar] [CrossRef] [PubMed]
  37. Ackermann, H.-W. Tailed Bacteriophages: The Order Caudovirales. Adv. Virus Res. 1998, 51, 135–201. [Google Scholar] [CrossRef] [PubMed]
  38. Hultman, J.; Tamminen, M.; Pärnänen, K.; Cairns, J.; Karkman, A.; Virta, M. Host range of antibiotic resistance genes in wastewater treatment plant influent and effluent. FEMS Microbiol. Ecol. 2018, 94, fiy038. [Google Scholar] [CrossRef]
  39. Pazda, M.; Kumirska, J.; Stepnowski, P.; Mulkiewicz, E. Antibiotic resistance genes identified in wastewater treatment plant systems – A review. Sci. Total Environ. 2019, 697, 134023. [Google Scholar] [CrossRef] [PubMed]
  40. Tang, J.; Bu, Y.; Zhang, X.-X.; Huang, K.; He, X.; Ye, L.; Shan, Z.; Ren, H. Metagenomic analysis of bacterial community composition and antibiotic resistance genes in a wastewater treatment plant and its receiving surface water. Ecotoxicol. Environ. Saf. 2016, 132, 260–269. [Google Scholar] [CrossRef] [PubMed]
  41. Versporten, A.; Zarb, P.; Caniaux, I.; Gros, M.-F.; Drapier, N.; Miller, M.; Jarlier, V.; Nathwani, D.; Goossens, H.; Koraqi, A.; et al. Antimicrobial consumption and resistance in adult hospital inpatients in 53 countries: Results of an internet-based global point prevalence survey. Lancet Global Health 2018, 6, e619–e629. [Google Scholar] [CrossRef] [PubMed]
  42. Nnadozie, C.F.; Kumari, S.; Bux, F. Status of pathogens, antibiotic resistance genes and antibiotic residues in wastewater treatment systems. Rev. Environ. Sci. Bio/Technol. 2017, 16, 491–515. [Google Scholar] [CrossRef]
  43. Kumar, N.; Radhakrishnan, A.; Wright, C.C.; Chou, T.-H.; Lei, H.-T.; Bolla, J.R.; Tringides, M.L.; Rajashankar, K.R.; Su, C.-C.; Purdy, G.E.; et al. Crystal structure of the transcriptional regulator Rv1219c of Mycobacterium tuberculosis. Protein Sci. 2014, 23, 423–432. [Google Scholar] [CrossRef] [PubMed]
  44. Lin, J.; Nishino, K.; Roberts, M.C.; Tolmasky, M.; Aminov, R.I.; Zhang, L. Mechanisms of antibiotic resistance. Front. Microbiol. 2015, 6, 34. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Analysis of metagenomics sequence diversity of samples, visualized by Origin 2022. (A) Sparsity curve analysis of influent and effluent samples; (B) ACE analysis of influent and effluent samples; and (C) Shannon analysis of influent and effluent samples.
Figure 1. Analysis of metagenomics sequence diversity of samples, visualized by Origin 2022. (A) Sparsity curve analysis of influent and effluent samples; (B) ACE analysis of influent and effluent samples; and (C) Shannon analysis of influent and effluent samples.
Water 15 01587 g001
Figure 2. Metagenomics analysis of bacterial communities in influent and effluent, visualized by Origin 2022. The top 10 abundant phyla (A); the top 15 abundant families (B); and the top 24 abundant genera (C). Relative abundance is expressed as a percentage of the total valid sequences for each sample for each bacterial taxonomic units.
Figure 2. Metagenomics analysis of bacterial communities in influent and effluent, visualized by Origin 2022. The top 10 abundant phyla (A); the top 15 abundant families (B); and the top 24 abundant genera (C). Relative abundance is expressed as a percentage of the total valid sequences for each sample for each bacterial taxonomic units.
Water 15 01587 g002
Figure 3. Metagenomics analysis of archaeal and viral communities at genus level, visualized by Origin 2022. The top 7 abundant archaeal genera (A) and the top 15 abundant viral genera (B). Relative abundance is expressed as a percentage of the total valid sequences for each sample for each archaeal/or viral taxonomic unit. Since viruses were not annotated at the genus level in the effluent, Figure 3B only shows the abundance of virus genera in the influent.
Figure 3. Metagenomics analysis of archaeal and viral communities at genus level, visualized by Origin 2022. The top 7 abundant archaeal genera (A) and the top 15 abundant viral genera (B). Relative abundance is expressed as a percentage of the total valid sequences for each sample for each archaeal/or viral taxonomic unit. Since viruses were not annotated at the genus level in the effluent, Figure 3B only shows the abundance of virus genera in the influent.
Water 15 01587 g003
Figure 4. The distribution structure of ARG subtypes (relative abundance > 1%) in the influent and effluent samples, visualized by Circos software. The names of influent and effluent samples and ARG subtypes are listed in the outside circle. The connecting lines within the circles relate the ARG subtypes to the sample, and the width of each line indicates the abundance of each ARG subtype in the influent and effluent samples. The color of each line was used to distinguish the association between the sample and the ARG subtypes.
Figure 4. The distribution structure of ARG subtypes (relative abundance > 1%) in the influent and effluent samples, visualized by Circos software. The names of influent and effluent samples and ARG subtypes are listed in the outside circle. The connecting lines within the circles relate the ARG subtypes to the sample, and the width of each line indicates the abundance of each ARG subtype in the influent and effluent samples. The color of each line was used to distinguish the association between the sample and the ARG subtypes.
Water 15 01587 g004
Figure 5. Abundance of antibiotic resistance genes, visualized by Origin 2022. (A) Abundance of ARGs classified by resistance mechanism in influent and effluent samples. ATA-AEF: antibiotic target alteration-antibiotic efflux; ATA-AEF-RPA: antibiotic target alteration-antibiotic efflux-reduced permeability to antibiotic; (B) Percent abundance (relative abundance > 1%) of influent and effluent samples classified by Drug Class in samples. MLS, MS, CP, MC, Fcgptrpt, SS, GT and CCP represent multi-resistance categories, as shown in Abbreviation.
Figure 5. Abundance of antibiotic resistance genes, visualized by Origin 2022. (A) Abundance of ARGs classified by resistance mechanism in influent and effluent samples. ATA-AEF: antibiotic target alteration-antibiotic efflux; ATA-AEF-RPA: antibiotic target alteration-antibiotic efflux-reduced permeability to antibiotic; (B) Percent abundance (relative abundance > 1%) of influent and effluent samples classified by Drug Class in samples. MLS, MS, CP, MC, Fcgptrpt, SS, GT and CCP represent multi-resistance categories, as shown in Abbreviation.
Water 15 01587 g005
Figure 6. ARGs and microbial co-occurrence network based on correlation analysis, visualized by Gephi 0.9. Nodes represent ARG subtypes or species, and lines represent p-values < 0.05 and a strong correlation. The size of a node is proportional to its number of connections. Nodes are colored by the type of resistance mechanism of ARG subtypes and by the type of microorganism. Bacterial and ARG co-occurrence network (A). (B) Archaeal co-occurrence network with viruses and ARGs.
Figure 6. ARGs and microbial co-occurrence network based on correlation analysis, visualized by Gephi 0.9. Nodes represent ARG subtypes or species, and lines represent p-values < 0.05 and a strong correlation. The size of a node is proportional to its number of connections. Nodes are colored by the type of resistance mechanism of ARG subtypes and by the type of microorganism. Bacterial and ARG co-occurrence network (A). (B) Archaeal co-occurrence network with viruses and ARGs.
Water 15 01587 g006
Table 1. Distribution of taxonomic units in influent and effluent samples a (%).
Table 1. Distribution of taxonomic units in influent and effluent samples a (%).
BacteriaArchaeaVirusEukaryotaUnassigned
InfluentPhylum99.12430.00867-0.00300.7860
Class99.34400.00867-0.00230.5670
Order99.45800.08670.01030.00270.4423
Family99.47730.08570.00200.00730.4277
Genus99.61930.08700.00230.00200.2893
Species99.59570.08370.00430.00200.3143
EffluentPhylum97.45400.7853-0.12971.6310
Class97.83330.7803-0.12931.2570
Order98.07370.78270.00630.13271.0047
Family99.13800.78370.00500.13200.9413
Genus99.36170.78770.00000.13430.7163
Species98.56070.77100.00030.18970.4783
Note(s): a: The abundance in percentage was based on the taxonomic units using the identified DNA gene tags.
Table 2. Number of taxonomic units in influent and effluent samples.
Table 2. Number of taxonomic units in influent and effluent samples.
BacteriaArchaeaVirusEukaryota
InfluentPhylum233-3
Class404-5
Order94617
Family1959411
Genus680142512
Species272326159
EffluentPhylum193-3
Class384-9
Order886111
Family1748314
Genus41012NA11
Species97117110
Note(s): NA: Not annotated.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Shen, M.; Hu, X.; Li, M.; Lyu, C.; Hu, Y.; Bu, X.; Chen, T.; Cai, H.; Li, C.; Liu, J.; et al. Distribution of Antibiotic Resistance Genes and Their Association with Microbes in Wastewater Treatment Plants: A Metagenomics Analysis. Water 2023, 15, 1587. https://doi.org/10.3390/w15081587

AMA Style

Shen M, Hu X, Li M, Lyu C, Hu Y, Bu X, Chen T, Cai H, Li C, Liu J, et al. Distribution of Antibiotic Resistance Genes and Their Association with Microbes in Wastewater Treatment Plants: A Metagenomics Analysis. Water. 2023; 15(8):1587. https://doi.org/10.3390/w15081587

Chicago/Turabian Style

Shen, Mengnan, Xiaowei Hu, Ming Li, Chen Lyu, Yi Hu, Xiaodan Bu, Tao Chen, Hang Cai, Chenyang Li, Jiahong Liu, and et al. 2023. "Distribution of Antibiotic Resistance Genes and Their Association with Microbes in Wastewater Treatment Plants: A Metagenomics Analysis" Water 15, no. 8: 1587. https://doi.org/10.3390/w15081587

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop