Next Article in Journal
Novel Criteria to Provide a Locality/Normality Degree in Molecules and Their Relevance in Physical Chemistry
Previous Article in Journal
Structural Unfolding of G-Quadruplexes: From Small Molecules to Antisense Strategies
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effect of Stepwise Exposure to High-Level Erythromycin on Anaerobic Digestion

1
School of Environmental and Material Engineering, Yantai University, Yantai 264005, China
2
State Key Laboratory of Pollution Control and Resource Reuse, School of the Environment, Nanjing University, Nanjing 210023, China
3
Key Laboratory of Coastal Environment and Resources of Zhejiang Province, School of Engineering, Westlake University, Hangzhou 310030, China
4
Department of Chemical and Biochemical Engineering, Technical University of Denmark, DK-2800 Kongens Lyngby, Denmark
*
Author to whom correspondence should be addressed.
Molecules 2024, 29(15), 3489; https://doi.org/10.3390/molecules29153489
Submission received: 8 June 2024 / Revised: 12 July 2024 / Accepted: 23 July 2024 / Published: 25 July 2024

Abstract

:
High-level erythromycin (ERY) fermentation wastewater will pose serious threats to lake environments. Anaerobic digestion (AD) has advantages in treating high-level antibiotic wastewater. However, the fate of antibiotic resistance genes (ARGs) and microbial communities in AD after stepwise exposure to high-level ERY remains unclear. In this study, an AD reactor was first exposed to 0, 5, 10, 50, 100 and 200 mg/L ERY and then re-exposed to 0, 50, 200 and 500 mg/L ERY to investigate the effect of ERY on AD. The results show that AD could adapt to the presence of high-level ERY (500 mg/L) and could maintain efficient CH4 production after domestication with low-level ERY (50 mg/L). The AD process could achieve higher removal of ERY (>94%), regardless of the initial ERY concentration. ErmB and mefA, conferring resistance through target alteration and efflux pumps, respectively, were dominant in the AD process. The first exposure to ERY stimulated an increase in the total ARG abundance, while the AD process seemed to discourage ARG maintenance following re-exposure to ERY. ERY inhibited the process of acetoclastic methanogenesis, but strengthened the process of hydrogenotrophic methanogenesis. This work provides useful information for treating high-level ERY fermentation wastewater by the AD process.

1. Introduction

Erythromycin (ERY), a broad-spectrum macrolide antibiotic, is produced by the fermentation of Streptomyces species [1]. Since ERY is frequently used in human medicine and veterinary medicine, its demand is high. For example, the usage of ERY in the UK is high, at 48 tonnes/year [2]. However, due to incomplete extraction during the ERY production process, its residue in fermentation wastewater is high, reaching the 100–500 mg/L level [3]. ERY has been detected globally in sewage (max. conc. of 1037 μg/kg dry wt.), surface water (max. conc. of 3847 ng/L) and sediment (max. conc. of 175.38 μg/kg dry wt.) [4]. The accumulation of ERY in the environment may pose a significant risk to the ecosystem. Previous studies found that ERY could cause toxic effects on algae, such as Microcystis flos-aquae [5] and Pseudokirchneriella subcapitata [6], as well as fish, such as Sparus aurata L. [7] and Oncorhynchus mykiss [8]. ERY is included on the US EPA’s Drinking Water Contaminant Candidate List 4 (CCL-4) [9] and the EU’s Water List [10]. Therefore, it is urgent to deal with high-level ERY fermentation wastewater.
Various treatment methods have been used to remove ERY from wastewater [11,12,13]. For example, Gholamiyan et al. used magnetic-activated carbon as an adsorbent for the removal of ERY and found that 95% of ERY could be adsorbed [11]. Mohammed et al. synthesized CuFe2O4@Chitosan for photocatalytic degradation of ERY and found a high ERY degradation efficiency (90%) [12]. Chu et al. investigated the removal of ERY by ionizing radiation and found that the ERY removal rate could reach 87% [13]. However, the above treatment methods are expensive. In addition, the high chemical oxygen demand (COD) and total suspended solids (TSSs) in ERY fermentation wastewater limited their application in ERY wastewater treatment. Anaerobic digestion (AD) is considered an economically feasible method to deal with recalcitrant pollutants, owing to its low cost and high energy efficiency [14,15]. Spielmeyer et al. [16] investigated the removal of tetracyclines and sulfonamides by anaerobic fermentation and found that their elimination rates ranged from 14% to 89%. It has been reported that the AD process could achieve 100% elimination of tetracycline [17].
Nowadays, there are, however, only limited studies investigating the treatment of high-level ERY fermentation wastewater by the AD process. Zeynep Cetecioglu et al. [18] found that methane generation was almost completely inhibited when the ERY dosage was above 500 mg/L. They also found that low-level ERY (1–100 mg/L) did not affect the removal of volatile fatty acids, and high-level ERY (250–1000 mg/L) impaired the utilization of propionate and reduced the removal of butyrate [19]. The above-mentioned research mainly focused on the transformation of organic matter. Actually, during the AD process, high-level ERY may hinder the growth of bacteria by binding to the 50S ribosomal subunits, further affecting the microbial community structure [20]. Moreover, the presence of high-level ERY may induce or affect the amplification of antibiotic resistance genes (ARGs) (four types, i.e., efflux genes, methylase genes, phosphorylase genes and esterase genes) [21]. Wang et al. reported that high-level ERY caused an increase in the relative abundance of Sedimentibacter and Methanosarcina and accelerated methylase gene (erm A/T) amplification [22]. The relative abundance of Syntrophomonas, Syntrophorhabdus and Methanosaeta increased under the pressure of high-level clarithromycin (an ERY derivative) [23]. Methylase genes (ermF and ermB) were enriched with the addition of high-level lincomycin (macrolide–lincosamide–streptogramin (MLS) antibiotic) [24]. However, the mentioned study only evaluated the influence of the first exposure to high-level ERY on the AD process. Up to now, the effect of stepwise exposure to high-level ERY on ARGs and microbial communities of the AD system has been little reported.
In order to investigate the influence of stepwise exposure to high-level ERY on anaerobic fermentation, an AD reactor was first exposed to 0, 5, 10, 50, 100 and 200 mg/L ERY, and then re-exposed to 0, 50, 200 and 500 mg/L ERY. The main objectives of this study were to (1) explore methane production during stepwise exposure to high-level ERY; (2) evaluate the removal efficiency of ERY in the AD system; (3) investigate the fate of key ARGs’ abundance and discuss the resistant mechanism; and (4) determine the diversity and structure of the microbial communities. Our study will provide a theoretical basis for the treatment of high-level ERY fermentation wastewater by the AD process.

2. Results and Discussion

2.1. Methane Generation Performance in the Existence of ERY

Methane production is usually used as an indicator of AD performance. Figure 1 displays the effects of ERY on methane production during AD. As shown in Figure 1a, in the presence of 5 and 10 mg/L ERY, CH4 production was higher than that of the glucose control group during the whole experiment, indicating that low levels of ERY would promote CH4 production. Low concentrations of ERY may have limited the antibacterial spectra of the bacteria community [25]. Mustapha et al. [26] found that macrolide antibiotics (ERY at 15 mg/L, clarithromycin at 15 mg/L and roxithromycin at 15 mg/L) promoted CH4 production. However, when the concentration of ERY increased further (50, 100 and 200 mg/L), decreases in CH4 production were observed during the experiment. Therefore, it could be speculated that high levels of ERY could inhibit the activity of CH4 production. In addition, the ultimate CH4 production remained at similar levels, indicating that the inhibition effect of ERY on AD progressively recovered with extended exposure. The recovery of CH4 production might be due to the domestication of certain bacteria or the removal of ERY.
According to Figure 1b, 50–50 and 200–200 showed no apparent inhibition effects on CH4 production. This might be because bacteria had already adapted to the same concentration of ERY during phase Ⅰ inoculation. In addition, the inhibition effect of 0–500 on CH4 production was similar to that of high-level ERY in phase Ⅰ. Interestingly, CH4 production remained relatively efficient in 50–500, while it was negatively affected by 200–500. A low concentration of the antibiotic is more likely to foster microbial resistance [27]. Fan et al. [28] studied the causal correlation of ERY with resistance proliferation. The result suggested that the presence of ERY at low concentrations could increase the establishment of antibiotic resistance. Thus, compared with 200, bacteria might find it easier to develop resistance under 50 in phase Ⅰ and could quickly adapt to 500 in phase II.

2.2. The Removal of ERY during AD

Figure 2 illustrates the variation in ERY concentration during AD. After the anaerobic reaction of phase Ⅰ, the removal efficiency of ERY in 50 and 200 reached 95.9% and 99.8%, respectively. The results imply that ERY could be effectively removed, even though ERY was difficult to biodegrade under anaerobic conditions. Feng et al. [29] examined the removal of ERY during anaerobic digestion of pig manure. In their study, the removal efficiency of ERY reached 99.9% at day 40, which indicated that this antibiotic could be effectively degraded in the AD system. In addition, the above results show that the increase in ERY concentration led to the improvement in removal efficiency. One possible explanation is that the antibiotic as a carbon source was more likely utilized by microorganisms at higher concentrations. Yang et al. [30] explored the effect of antibiotic concentration on the removal of antibiotics and found that the removal efficiency of antibiotics (sulfamethoxazole, chlorotetracycline and penicillin) improved with increasing antibiotic concentration in a sequencing batch reactor. This might be the reason that the inhibition effect of ERY gradually weakened and could recover when prolonging the exposure time.
After the anaerobic reaction of phase II, the removal efficiency of ERY was still high. The removal efficiency of ERY was 94.5% in 0–500, which was not consistent with the concept that the higher the antibiotic concentration, the better the removal efficiency. This was likely due to the microbial acclimation in phase II or the inhibition of microorganism activity under a higher concentration of ERY. In addition, the removal efficiency of ERY in 200–500 (96.4%) was higher than that in 50–500 (94.5%), but CH4 production in 200–500 was lower than that in 50–500. It could be speculated that the efficient CH4 production in 50–500 might not be caused by ERY removal, but by bacteria domestication under ERY pressure during phase I.

2.3. Variations in ARG Types and Relative Abundance

Ten ERY resistance genes, including two phosphorylase genes (mphA-01, mphB), two esterase genes (ereA, ereB), four methylase genes (ermA, ermB, ermC, ermX) and two efflux pump genes (mefA, mefE), were detected. As shown in Figure 3, ermB and mefA were abundant in 0. The former confers resistance to ERY by target alteration due to methylation, while the latter pumps ERY into the extracellular environment [31]. The results indicate that target alteration and efflux pumps were the dominant resistance mechanisms against ERY in the original sludge. The total ARG abundance was enhanced under ERY loading conditions in phase I, which might be owing to the increase in ERY selection pressure on the microbial communities. Meanwhile, the addition of ERY substantially increased the relative abundance of ermB and mefA, suggesting that ERY exposure improved the influence of methylation and efflux pumps. However, phosphorylase genes and esterase genes displayed insignificant changes, as well as the remaining efflux pump genes and methylase genes. The phosphorylase genes and esterase genes could inactivate ERY by encoding ERY phosphotransferase and esterase, respectively [19]. The results demonstrate that antibiotic inactivation played a minor role in ERY resistance.
It is clear that ermB and mefA were still dominant in phase II. The total ARG abundance in 0–0 was much higher than that in 0, which might be caused by the microbial acclimation. In addition, 0–500 presented the highest enrichment of total ARGs. Interestingly, the total ARG abundance in 50–50 and 200–200 in phase II decreased compared with those in 50 and 200 in phase I, respectively. This might be due to the fact that the first exposure of the antibiotic was more likely to stimulate an increase in the total ARG abundance. Moreover, in comparison with 200–200 and 200–500, the total ARG abundance was much lower in 50–50 and 50–500. The ARGs expressed in response to ERY could be adapted and reserved in 50. It was likely that these ARGs were less sensitive to ERY pressure in 50–500.

2.4. Variations in Microbial Community Diversity and Structure

2.4.1. Microbial Community Diversity

The richness and diversity of the microbial communities were estimated using the alpha diversity indices. As shown in Figure 4a, the Chao1, Observe, Shannon and Simpson indices of 50 were smaller than those of 0, suggesting that the microbial richness and diversity decreased. A possible explanation is the adaptation of certain bacteria and the inhibition of other bacteria under 50. Unlike 50, the microbial richness and diversity increased in 200, which could be explained by the removal of ERY. A previous study showed that bacteria diversity was positively related to the biotransformation rate of micro-pollutants [32]. In addition, the above four indices of 0–0 decreased compared with those of 0. The addition of ERY reduced the microbial richness and diversity in phase II, among which 0–500 presented the lowest levels, indicating that microorganism activity was inhibited under a higher concentration of ERY. On the other hand, microbial richness and diversity indices of 50–500 were higher than those of 50–50, 200–200 and 200–500, which demonstrated the adaptation of certain bacteria during phase I. The similarity in the microbial communities was revealed by PCoA (Figure 4b). PCoA1 and PCoA2 together explained 56.47% of the total variation. As we can see, 0–0 was deviating considerably from 0, which was caused by the microbial acclimation. The microbial communities of 0 were similar to those of 50, but different from those of 200. In addition, PCoA did not reveal any difference between 0–0 and 50–500, but revealed a difference between 0–0 and 0–500. These results were consistent with the changes in CH4 production.

2.4.2. Microbial Community Structure

The microbial community composition was analyzed at the phylum and genus levels (Figure 5). The phyla in all of the samples contained Bacteroidota (17.6–55.6%), Firmicutes (17.3–47.4%), Synergistota (7.9–28.2%), Caldatribacteriota (0.8–7.2%), Cloacimonadota, Halobacterota, Proteobacteria, Thermotogota and WS1. Bacteroidota could degrade glucose, glycerol and amino acids to produce volatile fatty acids in the anaerobic digestion process [33]. Moreover, Bacteroidota were correlated with the hydrolysis of macromolecular organic compounds [33]. Firmicutes played a major part in hydrolysis and acidogenesis [13]. Additionally, Firmicutes were involved in acetogenesis/syntrophic acetate oxidation in the AD system [13]. Synergistota contained a large number of organic matter-degrading bacteria, contributing to the production of methane [34]. Bacteroidota, Firmicutes and Synergistota were often found in anaerobic digesters owing to their ability to decompose a range of organic compounds [35]. The phylum Halobacterota was the dominant archaea species in all of the samples, which is known as a group of methanogens [36]. Most methanogens in the phylum Halobacterota were reported to use methanol, acetate and hydrogen as electron donors. As shown in Figure 5a, Bacteroidota increased in 50 but decreased in 200, which might be due to the fact that the hydrolytic process associated with Bacteroidota was promoted at low-level ERY but inhibited at high-level ERY. The addition of ERY and the microbial acclimation both affected the composition of phyla, which is in agreement with the alpha diversity results (Figure 4a).
At the genus level, Proteiniphilum and Acetomicrobium, typical hydrolytic/fermentative bacteria, were dominant in all of the samples. In phase I, the relative abundance of Lentimicrobium, WS1 and Methanosarcina displayed a decreasing trend with the increase in ERY concentration, which demonstrated the inhibition effects of ERY on these bacteria. Lentimicrobium is a strictly anaerobic bacterium that can degrade glucose to acetate, malate, propionate, formate and hydrogen [37]. Methanosarcina is one kind of acetoclastic methanogen and can metabolize acetate to produce CH4 [38]. Its decrease indicated that ERY inhibited the acetoclastic methanogenesis. The relative abundance of DTU014, Aminobacterium and Alkaliphilus maintained an increasing trend under ERY pressure, indicating that these bacteria could resist the toxicity of ERY. Aminobacterium, one kind of amino-acid-degrading bacterium, is closely related to hydrogenotrophic methanogens [39]. Alkaliphilus is capable of producing glycoside hydrolase under anaerobic fermentation [40], which may contribute to the removal of ERY. This is because ERY, composed of the hexose sugar cladinose, is a candidate target of glycosylase enzymes [41]. In addition, DTU014 and Alkaliphilus have shown potential as syntrophic acetate-oxidizing bacteria [42], which could couple to perform hydrogenotrophic methanogenesis. The increased relative abundance of the above genera might be the reason for the recovery of CH4 production under 200.
In phase II, Fermentimonas and Acholeplasma, typical hydrolytic/acetogenic bacteria, were enriched when glucose changed into avicel. The relative abundance of MgMjR-022, Ruminofilibacter, Syntrophomonas and Proteiniborus in 0–500 was higher than that in 0–0, which demonstrated that they could endure ERY biotoxicity. On the contrary, the relative abundance of Fermentimonas and Acholeplasma in 0–500 was lower than that in 0–0. Fermentimonas can use carbohydrates and complex proteinaceous to produce acetic acid, hydrogen and carbon dioxide [43]. Acholeplasma produces fatty acids that can be used for methanogenesis [44]. The reduced relative abundance of these genera suggests that hydrolysis and acidogenesis were affected by ERY, in turn leading to the inhibition of CH4 production. It should be noted that the relative abundance of MgMjR-022, Ruminofilibacter, Fermentimonas, Acholeplasma and Aminobacterium was still quite high in 50–50, 50–500 and 200–200. Hence, the high relative abundance of hydrolytic and acetogenic bacteria might contribute to the production of CH4. Furthermore, Methanoculleus showed the highest relative abundance in 50–500. Methanoculleus is responsible for hydrogenotrophic methanogenesis and is correlated with syntrophic acetate oxidation [45]. The result suggests that hydrogenotrophic methanogens were facilitated under 50–500.
Overall, the effect of stepwise exposure to high-level ERY on AD resulted in the enrichment of Aminobacterium, MgMjR-022 and Ruminofilibacter. Those enriched bacteria might be associated with the changes in ARGs to a certain extent. In addition, based on the above analysis, it was inferred that the hydrogenotrophic methanogenesis pathway might be less disturbed than the acetoclastic methanogenesis pathway under 50 in phase I. A previous study found that hydrogenotrophic methanogens were less sensitive to 50 mg/L of ciprofloxacin compared with acetoclastic methanogens [33]. Thus, hydrogenotrophic methanogens in 50 could quickly adapt to the high-level ERY in phase II and played an important role in CH4 production.

2.5. Assessment of the Influence of Microbial Community on ARGs

The microbial community is believed to influence the variation in ARGs [46]. Network analysis is an essential method to reveal the correlation between ARGs and the microbial community [47]. In this study, the relative abundance of mefA was high, but the correlation between mefA and the microbial community was not found through network analysis, which might be caused by the small sample size. However, it has been reported that Fermentimonas was a potential host of mefA [48]. Thus, the increase in the abundance of mefA in 0–0 might be related to the enrichment of Fermentimonas. In addition, the increased abundance of ermB could be explained by the proliferation of the potential hosts, such as Alkaliphilus and Ruminofilibacter, according to previous reports [49]. Moreover, Acholeplasma was identified as a co-host of ermB and mefA which is saprotrophic and pathogenic [49]. Alkaliphilus and Acholeplasma belong to Firmicutes, while Fermentimonas and Ruminofilibacter belong to Bacteroidota. Those two phyla were reported to be responsible for carrying and transferring ARGs in a variety of environments [50]. In this study, mefA and erm B, as well as most of the above-mentioned hosts, were not enriched in 50–500, suggesting that bacteria were better adapted to the presence of the high-level antibiotic in phase II after domestication with the low-level antibiotic in phase I.

3. Materials and Methods

3.1. Experiment Setup

The inoculum used in the experiment was obtained from an anaerobic reactor treating sewage sludge at 37 °C. After storage in the incubator for 1 month, its methane potential was removed. The total solid (TS) of the inoculum was 17.8 g/L; the volatile solid (VS) was 10.6 g/L; and the pH was 8.40. Erythromycin (C37H67NO13) was purchased from Aladdin Chemistry Co., Ltd. (Shanghai, China).
Batch tests were carried out in 118 mL serum bottles. Two phases were included in the experiment. The procedure of phase I was as follows: 10 mL inoculum was added into the serum bottle; glucose was used as a carbon source to achieve the required COD concentration (5000 mg/L); 20 mL ERY was injected into the bottle to obtain a range of 0, 5, 10, 50, 100 and 200 mg/L. A blank test was performed with only inoculum. Thereafter, those bottles were purged with nitrogen to exclude oxygen and then closed using butyl rubber stoppers and aluminum crimps. The bottles were named blank, 0, 5, 10, 50, 100 and 200. Then, the prepared bottles were placed in an incubator at 37 °C. The experiment was terminated when the gas production was insignificant. All of the experiments were carried out in triplicate. The digestate from phase I (0, 50 and 200) was used as inoculum in phase II. The procedure of phase II was as follows: Avicel was used as a carbon source to achieve the required COD concentration (5000 mg/L). Then, 10 mL products of 0, 50 and 200 were spiked with 20 mL ERY in new bottles to obtain a range of 0 and 500 mg/L (named 0–0 and 0–500), 50 and 500 mg/L (named 50–50 and 50–500) as well as 200 and 500 mg/L (named 200–200 and 200–500). The remaining steps were the same as in phase I. In addition, the concentration of ERY was determined at the end of the experiment. Cs and Ce represent the initial and final concentrations of ERY, respectively. As for ARGs and microbial analysis, the products of 0, 50, 200, 0–0, 0–500, 50–50, 50–500, 200–200 and 200–500 were taken from the reactor.

3.2. Microbial Analysis

Firstly, the sample was washed using phenol to improve DNA purification. Then, DNA was extracted with DNeasy PowerSoil® (QIAGEN GmbH, Hilden, Germany). DNA quality was evaluated with NanoDrop (Thermo Scientific, Waltham, MA, USA). After that, 16S rRNA genes were amplified using the primers 515F/806R. Sequencing was conducted using the Illumina high-throughput sequencing technology [51].

3.3. ARG Analysis

The abundance of the ERY resistance gene was examined through RT-qPCR using the StepOnePlus™ Real-Time PCR system (Wcgene Biotechnology, Shanghai, China). The target genes included ereA, ereB, mphA-01, mphB, ermA, ermB, ermC, ermX, mefA and mefE. The ERY resistance gene primers are shown in the Supplementary Data. The detailed composition of the PCR mixture and the PCR procedures were processed according to a previous study [52]. The efficiency of amplification should be within 1.8–2.2. The relative abundance of ARGs was normalized to 16S rRNA genes, which was measured according to a comparative CT method, 2−ΔCT, where ΔCT = (CTARG − CT16S rRNA gene).

3.4. Analytical Methods

TS and VS were analyzed according to [53]. The pH was measured by a pH meter (FiveEasy Plus Benchtop FP20, Mettler Toledo, Greifensee, Switzerland). The content of CH4 was quantified by using a gas chromatograph (Trace 1310 GC-TCD, Thermo Fisher, Lillerød, Denmark) equipped with a thermal conductivity detector (TCD) and Thermo (P/N 26004–6030, Thermo Scientific, Waltham, MA, USA) Column (30 m length, 0.320 mm inner diameter and film thickness 10 µm) with helium as the carrier gas. The concentration of ERY was determined using a high-performance liquid chromatography (Agilent 1290 Infinity, Santa Clara, CA, USA, HPLC) system coupled with a tandem mass spectrometer (Agilent 6470 series, USA, MS/MS). The system consisted of electrospray ionization (ESI) and a C18 column (2.1 × 50 mm, 1.7 μm). The positive ion mode (ESI+) was selected. ERY was detected by the prominent fragment ion at m/z 158. Mobile phase A was MeOH and mobile phase B was water plus 0.1% formic acid (at a flow rate of 0.2 mL/min for 8 min). The elution gradient mode was as follows: 45% A (0–2.5 min), 25% A (2.5–4 min), 10% A (4–7 min) and 25% A (7–8 min). Principal coordinate analysis (PCoA), a heatmap and a chord diagram were conducted using R software (4.2.0). OriginPro 2021 was used to analyze other figures.

4. Conclusions

This study investigated the influence of stepwise exposure to high-level ERY on anaerobic digestion, and the main conclusions are as follows: (1) After domestication with low-level ERY, AD could adapt to the presence of high-level ERY and could maintain efficient CH4 production. (2) The AD process could achieve higher removal of ERY. (3) The first exposure to ERY stimulated an increase in the total ARG abundance, while the AD process seemed to discourage ARG maintenance following re-exposure to ERY. (4) Under the influence of ERY, the hydrogenotrophic methanogenesis pathway replaced the acetoclastic methanogenesis pathway for CH4 production. There are still some limitations in this study: the effect of the transformation products of ERY on the fate of ARGs and microorganisms needs further research, and the application of AD for the treatment of actual antibiotic fermentation wastewater should be clarified.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/molecules29153489/s1. Table S1: qPCR primer sequences for antibiotic resistance genes analysis. References [54,55,56] are cited in the Supplementary Materials.

Author Contributions

Y.Z.: investigation, formal analysis, data curation, writing—original draft. C.L.: writing—review and editing. X.Z.: formal analysis, writing—review and editing. I.A.: supervision, writing—review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the funding of the China Scholarship Council (CSC), China, and the Natural Science Foundation of Shandong Province (ZR2022QE138).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data are contained within the article or the Supplementary Materials.

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

References

  1. Dinos, G.P. The macrolide antibiotic renaissance. Br. J. Pharmacol. 2017, 174, 2967–2983. [Google Scholar] [CrossRef] [PubMed]
  2. Ashton, D.; Hilton, M.; Thomas, K.V. Investigating the environmental transport of human pharmaceuticals to streams in the United Kingdom. Sci. Total Environ. 2004, 333, 167–184. [Google Scholar] [CrossRef] [PubMed]
  3. Aydin, S.; Ince, B.; Cetecioglu, Z.; Arikan, O.; Ozbayram, E.G.; Shahi, A.; Ince, O. Combined effect of erythromycin, tetracycline and sulfamethoxazole on performance of anaerobic sequencing batch reactors. Bioresour. Technol. 2015, 186, 207–214. [Google Scholar] [CrossRef] [PubMed]
  4. Schafhauser, B.H.; Kristofco, L.A.; de Oliveira, C.M.R.; Brooks, B.W. Global review and analysis of erythromycin in the environment: Occurrence, bioaccumulation and antibiotic resistance hazards. Environ. Pollut. 2018, 238, 440–451. [Google Scholar] [CrossRef] [PubMed]
  5. Wan, J.; Guo, P.; Peng, X.; Wen, K. Effect of erythromycin exposure on the growth, antioxidant system and photosynthesis of Microcystis flos-aquae. J. Hazard. Mater. 2015, 283, 778–786. [Google Scholar] [CrossRef] [PubMed]
  6. Gonzalez-Pleiter, M.; Gonzalo, S.; Rodea-Palomares, I.; Leganes, F.; Rosal, R.; Boltes, K.; Marco, E.; Fernandez-Pinas, F. Toxicity of five antibiotics and their mixtures towards photosynthetic aquatic organisms: Implications for environmental risk assessment. Water Res. 2013, 47, 2050–2064. [Google Scholar] [CrossRef] [PubMed]
  7. Rodrigues, S.; Antunes, S.C.; Correia, A.T.; Golovko, O.; Zlabek, V.; Nunes, B. Assessment of toxic effects of the antibiotic erythromycin on the marine fish gilthead seabream (Sparus aurata L.) by a multi-biomarker approach. Chemosphere 2019, 216, 234–247. [Google Scholar] [CrossRef] [PubMed]
  8. Rodrigues, S.; Antunes, S.; Correia, A.; Nunes, B. Acute and chronic effects of erythromycin exposure on oxidative stress and genotoxicity parameters of Oncorhynchus mykiss. Sci. Total Environ. 2016, 545, 591–600. [Google Scholar] [CrossRef] [PubMed]
  9. USEPA. Drinking Water Contaminant Candidate List 4; USEPA: Washington, DC, USA, 2016; Volume 2024.
  10. EU. Commission Implementing Decision (EU) 2018/840 of 5 June 2018 Establishing a Watch List of Substances for Union-Wide Monitoring in the Field of Water Policy Pursuant to Directive 2008/105/EC of the European Parliament and of the Council and Repealing Commission Implementing Decision (EU) 2015/495. 2018. Available online: https://eur-lex.europa.eu/legal-content/EN/TXT/PDF/?uri=CELEX:32018D0840 (accessed on 20 July 2024).
  11. Gholamiyan, S.; Hamzehloo, M.; Farrokhnia, A. RSM optimized adsorptive removal of erythromycin using magnetic activated carbon: Adsorption isotherm, kinetic modeling and thermodynamic studies. Sustain. Chem. Pharm. 2020, 17, 100309. [Google Scholar] [CrossRef]
  12. Mohammed, N.A.H.; Shamma, R.N.; Elagroudy, S.; Adewuyi, A. Copper ferrite immobilized on chitosan: A suitable photocatalyst for the removal of ciprofloxacin, ampicillin and erythromycin in aqueous solution. Catal. Commun. 2023, 182, 106745. [Google Scholar] [CrossRef]
  13. Chu, L.; Wang, J.; Chen, C.; He, S.; Wojnárovits, L.; Takács, E. Advanced treatment of antibiotic wastewater by ionizing radiation combined with peroxymonosulfate/H2O2 oxidation. J. Clean. Prod. 2021, 321, 128921. [Google Scholar] [CrossRef]
  14. Guo, Y.; Askari, N.; Smets, I.; Appels, L. A review on co-metabolic degradation of organic micropollutants during anaerobic digestion: Linkages between functional groups and digestion stages. Water Res. 2024, 256, 121598. [Google Scholar] [CrossRef] [PubMed]
  15. Wang, Y.; Du, B.; Wu, G. Powdered activated carbon facilitated degradation of complex organic compounds and tetracycline in stressed anaerobic digestion systems. Bioresour. Technol. 2024, 400, 130672. [Google Scholar] [CrossRef] [PubMed]
  16. Spielmeyer, A.; Breier, B.; Groißmeier, K.; Hamscher, G. Elimination patterns of worldwide used sulfonamides and tetracyclines during anaerobic fermentation. Bioresour. Technol. 2015, 193, 307–314. [Google Scholar] [CrossRef] [PubMed]
  17. Shi, J.C.; Liao, X.D.; Wu, Y.B.; Liang, J. Effect of antibiotics on methane arising from anaerobic digestion of pig manure. Anim. Feed Sci. Technol. 2011, 166, 457–463. [Google Scholar] [CrossRef]
  18. Cetecioglu, Z.; Ince, B.; Orhon, D.; Ince, O. Acute inhibitory impact of antimicrobials on acetoclastic methanogenic activity. Bioresour. Technol. 2012, 114, 109–116. [Google Scholar] [CrossRef] [PubMed]
  19. Cetecioglu, Z.; Ince, B.; Ince, O.; Orhon, D. Acute effect of erythromycin on metabolic transformations of volatile fatty acid mixture under anaerobic conditions. Chemosphere 2015, 124, 129–135. [Google Scholar] [CrossRef] [PubMed]
  20. Alighardashi, A.; Pandolfi, D.; Potier, O.; Pons, M.N. Acute sensitivity of activated sludge bacteria to erythromycin. J. Hazard. Mater. 2009, 172, 685–692. [Google Scholar] [CrossRef] [PubMed]
  21. Zeng, S.; Sun, J.; Chen, Z.; Xu, Q.; Wei, W.; Wang, D.; Ni, B.J. The impact and fate of clarithromycin in anaerobic digestion of waste activated sludge for biogas production. Environ. Res. 2021, 195, 110792. [Google Scholar] [CrossRef]
  22. Wang, M.; Ren, P.; Wang, Y.; Cai, C.; Liu, H.; Dai, X. Erythromycin stimulates rather than inhibits methane production in anaerobic digestion of antibiotic fermentation dregs. Sci. Total Environ. 2022, 807 Pt 3, 151007. [Google Scholar] [CrossRef] [PubMed]
  23. Huang, X.; Liu, X.; Chen, F.; Wang, Y.; Li, X.; Wang, D.; Tao, Z.; Xu, D.; Xue, W.; Geng, M.; et al. Clarithromycin affect methane production from anaerobic digestion of waste activated sludge. J. Clean. Prod. 2020, 255, 120321. [Google Scholar] [CrossRef]
  24. Zhu, W.; Bu, F.; Xu, J.; Wang, Y.; Xie, L. Influence of lincomycin on anaerobic digestion: Sludge type, biogas generation, methanogenic pathway and resistance mechanism. Bioresour. Technol. 2021, 329, 124913. [Google Scholar] [CrossRef]
  25. Roberts, M.C. Update on macrolide–lincosamide–streptogramin, ketolide, and oxazolidinone resistance genes. FEMS Microbiol. Lett. 2008, 282, 147–159. [Google Scholar] [CrossRef] [PubMed]
  26. Mustapha, N.A.; Hu, A.; Yu, C.P.; Sharuddin, S.S.; Ramli, N.; Shirai, Y.; Maeda, T. Seeking key microorganisms for enhancing methane production in anaerobic digestion of waste sewage sludge. Appl. Microbiol. Biotechnol. 2018, 102, 5323–5334. [Google Scholar] [CrossRef] [PubMed]
  27. Sandegren, L. Selection of antibiotic resistance at very low antibiotic concentrations. Ups. J. Med. 2014, 119, 103–107. [Google Scholar] [CrossRef] [PubMed]
  28. Fan, C.; He, J. Proliferation of antibiotic resistance genes in microbial consortia of sequencing batch reactors (SBRs) upon exposure to trace erythromycin or erythromycin-H2O. Water Res. 2011, 45, 3098–3106. [Google Scholar] [CrossRef] [PubMed]
  29. Feng, L.; Casas, M.E.; Ottosen, L.D.M.; Moller, H.B.; Bester, K. Removal of antibiotics during the anaerobic digestion of pig manure. Sci. Total Environ. 2017, 603–604, 219–225. [Google Scholar] [CrossRef] [PubMed]
  30. Yang, L.; Wen, Q.; Zhao, Y.; Chen, Z.; Wang, Q.; Burgmann, H. New insight into effect of antibiotics concentration and process configuration on the removal of antibiotics and relevant antibiotic resistance genes. J. Hazard. Mater. 2019, 373, 60–66. [Google Scholar] [CrossRef] [PubMed]
  31. Ni, B.J.; Zeng, S.; Wei, W.; Dai, X.; Sun, J. Impact of roxithromycin on waste activated sludge anaerobic digestion: Methane production, carbon transformation and antibiotic resistance genes. Sci. Total Environ. 2020, 703, 134899. [Google Scholar] [CrossRef] [PubMed]
  32. Johnson, D.R.; Helbling, D.E.; Lee, T.K.; Park, J.; Fenner, K.; Kohler, H.P.; Ackermann, M. Association of biodiversity with the rates of micropollutant biotransformations among full-scale wastewater treatment plant communities. Appl. Environ. Microbiol. 2015, 81, 666–675. [Google Scholar] [CrossRef] [PubMed]
  33. Chen, Y.; Zhao, Z.; Zou, H.; Yang, H.; Sun, T.; Li, M.; Chai, H.; Li, L.; Ai, H.; Shi, D.; et al. Digestive performance of sludge with different crop straws in mesophilic anaerobic digestion. Bioresour. Technol. 2019, 289, 121595. [Google Scholar] [CrossRef] [PubMed]
  34. Xiaomei, Z.; Rujing, L.; Jun, X.; Yingying, H.; Xinying, Z.; Li, X. Enhanced methane production by bimetallic metal–organic frameworks (MOFs) as cathode in an anaerobic digestion microbial electrolysis cell. Chem. Eng. J. 2022, 440, 135799. [Google Scholar] [CrossRef]
  35. Mai, D.T.; Stuckey, D.C.; Oh, S. Effect of ciprofloxacin on methane production and anaerobic microbial community. Bioresour. Technol. 2018, 261, 240–248. [Google Scholar] [CrossRef] [PubMed]
  36. Wang, Y.; Wei, W.; Huang, Q.-S.; Ni, B.-J. Methane production from algae in anaerobic digestion: Role of corncob ash supplementation. J. Clean. Prod. 2021, 327, 129485. [Google Scholar] [CrossRef]
  37. Fan, C.; Zhou, M.; Tang, X.; Zeng, G.; Xu, Q.; Song, B.; Gong, R.; Zhang, B.; Xiong, W.; Lu, Y.; et al. Triclosan enhances short-chain fatty acid production from sludge fermentation by elevating transcriptional activity of acidogenesis bacteria. Chem. Eng. J. 2020, 384, 123285. [Google Scholar] [CrossRef]
  38. Wang, H.; Li, J.; Zhao, Y.; Xu, C.; Zhang, K.; Li, J.; Yan, L.; Gu, J.D.; Wei, D.; Wang, W. Establishing practical strategies to run high loading corn stover anaerobic digestion: Methane production performance and microbial responses. Bioresour. Technol. 2020, 310, 123364. [Google Scholar] [CrossRef] [PubMed]
  39. Wu, B.; Lin, R.; Kang, X.; Deng, C.; Xia, A.; Dobson, A.D.W.; Murphy, J.D. Graphene Addition to Digestion of Thin Stillage Can Alleviate Acidic Shock and Improve Biomethane Production. ACS Sustain. Chem. Eng. 2020, 8, 13248–13260. [Google Scholar] [CrossRef]
  40. Xie, J.; Chen, Y.; Duan, X.; Feng, L.; Yan, Y.; Wang, F.; Zhang, X.; Zhang, Z.; Zhou, Q. Activated carbon promotes short-chain fatty acids production from algae during anaerobic fermentation. Sci. Total Environ. 2019, 658, 1131–1138. [Google Scholar] [CrossRef] [PubMed]
  41. Gonzalez-Gil, L.; Krah, D.; Ghattas, A.K.; Carballa, M.; Wick, A.; Helmholz, L.; Lema, J.M.; Ternes, T.A. Biotransformation of organic micropollutants by anaerobic sludge enzymes. Water Res. 2019, 152, 202–214. [Google Scholar] [CrossRef]
  42. Dyksma, S.; Jansen, L.; Gallert, C. Syntrophic acetate oxidation replaces acetoclastic methanogenesis during thermophilic digestion of biowaste. Microbiome 2020, 8, 105. [Google Scholar] [CrossRef] [PubMed]
  43. Wang, J.; Feng, K.; Lou, Y.; Lu, B.; Liu, B.; Xie, G.; Ren, N.; Xing, D. The synergistic effect of potassium ferrate and peroxymonosulfate application on biogas production and shaping microbial community during anaerobic co-digestion of a cow manure-cotton straw mixture. Bioresour. Technol. 2021, 333, 125166. [Google Scholar] [CrossRef] [PubMed]
  44. Li, K.; Wang, K.; Wang, J.; Yuan, Q.; Shi, C.; Wu, J.; Zuo, J. Performance assessment and metagenomic analysis of full-scale innovative two-stage anaerobic digestion biogas plant for food wastes treatment. J. Clean. Prod. 2020, 264, 121646. [Google Scholar] [CrossRef]
  45. Tong, J.; Fang, P.; Zhang, J.; Wei, Y.; Su, Y.; Zhang, Y. Microbial community evolution and fate of antibiotic resistance genes during sludge treatment in two full-scale anaerobic digestion plants with thermal hydrolysis pretreatment. Bioresour. Technol. 2019, 288, 121575. [Google Scholar] [CrossRef] [PubMed]
  46. Sun, H.; Zhang, H.; Wu, D.; Ding, J.; Niu, Y.; Jiang, T.; Yang, X.; Liu, Y. Deciphering the antibiotic resistome and microbial community in municipal wastewater treatment plants at different elevations in eastern and western China. Water Res. 2023, 229, 119461. [Google Scholar] [CrossRef] [PubMed]
  47. Hui, X.; Fang, W.; Wang, G.; Liu, H.; Dai, X. Waste recycling of antibiotic mycelial residue: The feasible harmless treatment and source control of antibiotic resistance. J. Clean. Prod. 2023, 401, 136789. [Google Scholar] [CrossRef]
  48. Ma, J.; Wang, P.; Gu, W.; Su, Y.; Wei, H.; Xie, B. Does lipid stress affect performance, fate of antibiotic resistance genes and microbial dynamics during anaerobic digestion of food waste? Sci. Total Environ. 2021, 756, 143846. [Google Scholar] [CrossRef] [PubMed]
  49. Su, Y.; Wang, J.; Xia, H.; Xie, B. Comparative network analysis revealing the mechanisms of antibiotic resistance genes removal by leachate recirculation under different hydraulic loadings. Sci. Total Environ. 2019, 649, 318–326. [Google Scholar] [CrossRef] [PubMed]
  50. Syafiuddin, A.; Boopathy, R. Role of anaerobic sludge digestion in handling antibiotic resistant bacteria and antibiotic resistance genes—A review. Bioresour. Technol. 2021, 330, 124970. [Google Scholar] [CrossRef] [PubMed]
  51. Li, C.; Wang, Y.; Xie, S.; Wang, R.; Sheng, H.; Yang, H.; Yuan, Z. Synergistic treatment of sewage sludge and food waste digestate residues for efficient energy recovery and biochar preparation by hydrothermal pretreatment, anaerobic digestion, and pyrolysis. Appl. Energ. 2024, 364, 123203. [Google Scholar] [CrossRef]
  52. Zhang, Y.; Liu, H.; Dai, X.; Cai, C.; Wang, J.; Wang, M.; Shen, Y.; Wang, P. Impact of application of heat-activated persulfate oxidation treated erythromycin fermentation residue as a soil amendment: Soil chemical properties and antibiotic resistance. Sci. Total Environ. 2020, 736, 139668. [Google Scholar] [CrossRef] [PubMed]
  53. Federation, W.E.; Association, A. Standard Methods for the Examination of Water and Wastewater; American Public Health Association (APHA): Washington, DC, USA, 2005. [Google Scholar]
  54. Guo, M.T.; Yuan, Q.B.; Yang, J. Ultraviolet reduction of erythromycin and tetracycline resistant heterotrophic bacteria and their resistance genes in municipal wastewater. Chemosphere 2013, 93, 2864–2868. [Google Scholar] [CrossRef] [PubMed]
  55. Liu, M.; Ding, R.; Zhang, Y.; Gao, Y.; Tian, Z.; Zhang, T.; Yang, M. Abundance and distribution of Macrolide-Lincosamide-Streptogramin resistance genes in an anaerobic-aerobic system treating spiramycin production wastewater. Water Res. 2014, 63, 33–41. [Google Scholar] [CrossRef] [PubMed]
  56. Chen, Y.; Li, P.; Huang, Y.; Yu, K.; Chen, H.; Cui, K.; Huang, Q.; Zhang, J.; Gin, K.Y.-H.; He, Y. Environmental media exert a bottleneck in driving the dynamics of antibiotic resistance genes in modern aquatic environment. Water Res. 2019, 162, 127–138. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Effects of ERY on methane production during anaerobic digestion. (a) Phase I. (b) Phase II.
Figure 1. Effects of ERY on methane production during anaerobic digestion. (a) Phase I. (b) Phase II.
Molecules 29 03489 g001
Figure 2. Variation in ERY concentration during anaerobic digestion (Cs and Ce represent initial and final concentrations of ERY, respectively).
Figure 2. Variation in ERY concentration during anaerobic digestion (Cs and Ce represent initial and final concentrations of ERY, respectively).
Molecules 29 03489 g002
Figure 3. Chord diagram showing relative abundance of ARGs in different anaerobic digestion treatments.
Figure 3. Chord diagram showing relative abundance of ARGs in different anaerobic digestion treatments.
Molecules 29 03489 g003
Figure 4. Variation in microbial community diversity and structure. (a) Alpha diversity. (b) Principal coordinate analysis (PCoA) based on Bray–Curtis distance showing distribution of microbial community.
Figure 4. Variation in microbial community diversity and structure. (a) Alpha diversity. (b) Principal coordinate analysis (PCoA) based on Bray–Curtis distance showing distribution of microbial community.
Molecules 29 03489 g004aMolecules 29 03489 g004b
Figure 5. Changes in microbial community composition in different anaerobic digestion treatments. (a) Microbial community composition at phylum level. (b) Heatmap showing 50 genera with highest relative abundance in different anaerobic digestion treatments.
Figure 5. Changes in microbial community composition in different anaerobic digestion treatments. (a) Microbial community composition at phylum level. (b) Heatmap showing 50 genera with highest relative abundance in different anaerobic digestion treatments.
Molecules 29 03489 g005aMolecules 29 03489 g005b
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

Zhang, Y.; Li, C.; Zhu, X.; Angelidaki, I. Effect of Stepwise Exposure to High-Level Erythromycin on Anaerobic Digestion. Molecules 2024, 29, 3489. https://doi.org/10.3390/molecules29153489

AMA Style

Zhang Y, Li C, Zhu X, Angelidaki I. Effect of Stepwise Exposure to High-Level Erythromycin on Anaerobic Digestion. Molecules. 2024; 29(15):3489. https://doi.org/10.3390/molecules29153489

Chicago/Turabian Style

Zhang, Yanxiang, Chunxing Li, Xinyu Zhu, and Irini Angelidaki. 2024. "Effect of Stepwise Exposure to High-Level Erythromycin on Anaerobic Digestion" Molecules 29, no. 15: 3489. https://doi.org/10.3390/molecules29153489

APA Style

Zhang, Y., Li, C., Zhu, X., & Angelidaki, I. (2024). Effect of Stepwise Exposure to High-Level Erythromycin on Anaerobic Digestion. Molecules, 29(15), 3489. https://doi.org/10.3390/molecules29153489

Article Metrics

Back to TopTop