1. Introduction
Nowadays, a wide variety of antibiotics are being used in animal farms to cure, prevent and also to improve the growth of animals, accounting for more than 52% of total antibiotics consumption in the world [
1,
2,
3]. Due to the rapid effect of antibiotics and low cost of them, daily use of it rocketed during the last two decades) [
4]. Although in 2017, attempts to restrict using antibiotics took place, from 2009 to 2016 use of several veterinary antibiotics (VAs) was raised by 36.8% on average [
4]. Many of these compounds have weak absorption within the animal gut and intestine during digestion, resulting in the excretion of potent parent and daughter products [
5]. A high percentage of the antibiotics (60–90%) is excreted without metabolism in urine and feces, leading to potential human and ecological health risks for soil and water [
6,
7,
8,
9]. In addition, based on a study of Alexy et al. [
10], most of the antibiotics being used are not biodegradable, with degradation extents varying between 4 and 27%.
Since a vast volume of manure is being produced each year due to concentrated animal feeding operations (CAFOs) and is mostly being applied to solid materials as a fertilizer, the long-term presence of such antibiotics in manure with even trace concentrations (i.e., ng/L) could lead to the formation of antibiotic resistance genes (ARGs) [
11]. The microbacterial resistance will result in higher medical costs, longer treatment periods, and increased mortality [
12,
13]. Animal farms typically utilize simple treatment systems, which are mainly AD, stockpiling, composting, wetlands, or lagoons ([
3,
14]). These simple treatments might not be sufficient to prevent the appearance of ARGs [
15]. There are new approaches to remove or reduce antibiotics, including activated carbon adsorption, membrane filtration, advanced oxidation processes ([
14,
16,
17,
18,
19]), however, they require advanced technical supervision as well as extreme expenses. Furthermore, the VAs could pollute the soil and water, then the human food chain through crops and animal-derived foods [
20,
21,
22,
23,
24]. Moreover, the residue of VAs in the AD process could sustain microbes under the minimum inhibitory concentration (MIC), fostering the selection for ARGs by microbes [
25,
26].
Anaerobic digesters can produce biogas as well as removing VAs and ARGs [
27,
28,
29]. In contrast, composting [
30] requires less monitoring, has a more stable working system and a generally higher removal rate for certain VAs, but consumes energy and occasionally has less capability to remove ARGs than AD. Xie et al. [
31] concluded that thermophilic composting of cow manure would result in ARG mitigation, lowering 16S rRNA with tetracycline, sulfonamide and fluoroquinolone resistance genes, however, not effective with aadA, aadA2, qacED1, tetL, cintI1, intI1, and tnpA04. A similar study of dairy manure composts showed satisfactory treatment of antibiotic-resistant E. coli and Salmonella, yet some antibiotic-resistant Enterobacter spp. and multidrug-resistant Pseudomonas spp. population raised after application of these composts to rangeland soils in Texas [
29]. Prado et al. [
32] aimed to use an aerobic reactor with activated sludge to track the fate of tetracycline (TC) and Tylosin as antibiotics. Both TC and Tylosin were not biodegradable in this type of reactor. Research also determined that the biosorption of both antibiotics appeared to be most favorable for TC.
Joy et al. [
33] investigated the behavior of three antibiotics (bacitracin, chlortetracycline, and tylosin) and two classes of ARGs (Tet and Erm), which were monitored in swine manure slurry under anaerobic conditions. First-order decay rates were determined for each antibiotic with half-lives ranging from 1 day (chlortetracycline) to 10 days (tylosin). Angenent and Wrenn [
34] examined the effects of an anaerobic sequencing batch reactor (ASBR) on the removal of antibiotic tylosin. They observed no inhibitory effect on biogas production, but some macrolide–lincosamide– streptogramin B (MLSB)-resistant bacteria appeared. Shi et al. [
35], discovered that a certain dosage of tetracycline (TC) and sulfamethoxydiazine (SMD) could reduce biogas production. They also noticed the rapid disappearance of antibiotics (more than 50%) in the first 12 h. However, they were not sure about whether it was being degraded or just absorbed into solid materials. A similar study was conducted by Beneragama et al. [
36], who confirmed the efficiency of AD of antibiotics in dairy manure. They also utilized thermophilic microorganisms (working in 55 °C or 131 °F). Results showed no inhibition in gas production and the efficiency of the reactor.
Approximately 80% of the 16,000 metric tons of antibiotics sold annually in the U.S. are used in animal husbandry [
37]. These antibiotics can be transported to run-off water, groundwater, soil, and finally, plants [
38,
39,
40,
41,
42,
43]. In 2030, antibiotic resistance would cost USD 3.4 trillion due to subsequent mortality and substitute treatment [
44]. Therefore, according to this extensive use of antibiotics and their stability in the environment, seeking an efficient and financially feasible method is vital. Besides, limited findings are available about the threshold concentrations of different classes of antibiotics in manure that can be removed during the AD process and the interaction between anaerobic digesters and antibiotics.
The objectives of this study were to evaluate (i) anaerobic digestion efficiency on the removal of chlortetracycline (CTC) and Tylosin, (ii) inhibitory behavior of VAs on the reactors, and (iii) the effect of these antibiotics on microbial dynamics in anaerobic digesters with swine manure. The novelty of this research is the imitation of on-farm mesophilic anaerobic digesters that were loaded frequently with manure from a commercial pig farm, while operating over several months. The research focused on widely used antibiotics and emphasized proper concentration and duration at which the antibiotics were administered (following manufacturer’s recommendation and average pig weights), and amount of antibiotics excreted by the animals based on a literature review. Critical operating variables of the digesters, including pH and biogas productions were monitored closely, similar to what most on-farm AD technicians are employing to monitor AD performance, which require no sophisticated analytical expertise. The bench-scale anaerobic digesters were relatively larger and semi-continuously loaded for over several months, while many of the previous studies only focused on inhibition effect and usually using batch reactors, which are different from actual on-farm AD conditions.
2. Methods
The current study focuses on utilizing an AD process to assist in removing VAs and finding the efficiency and practicality of the reactors. Antibiotics are chosen to be spiked meticulously, based on their importance, usage in feedstock and their danger to the environment. Besides, the dosage of antibiotics was close to concentrations administrated for animals, absorbed, and then excreted, to imitate the real condition. Our reactors are fed with swine manure, which has been tested for background concentration of VAs, to diminish the chance of interference. CTC was injected into anaerobic sequencing batch reactors (ASBR) with doses of 263, 280, and 298 mg/L for each spike and a total of three injections every two days while Tylosin doses were 88, 175, and 263 mg/L for each injection and a total of 5 injections every two days.
Several factors are contributing to reactor performance and biogas production. Temperature is one of them. Different types of bacteria work on various temperature ranges and some of them are highly susceptible to temperature fluctuation. Besides, the pH and alkalinity of the environment in which bacteria are growing should be near neutral and consistent. Therefore, the temperature was kept at around 39 °C (102 °F). In addition, we were recording incubator temperature and humidity for tracking the performance of our incubator. The methane forming bacteria are very sensitive to slight changes in organic loading, pH, and temperature (a temperature change greater than 2 degrees of Fahrenheit per day will affect the methane formers).
2.1. Feedstock for Anaerobic Digestion
Manure samples were collected from a mid-central Missouri commercial swine farm. The farm was VA free for the finishing pigs, located in Versailles, Missouri, USA. Furthermore, to make sure that no antibiotics existed in the solid manure used, it was analyzed to eradicate any interference or error. Because the farm has shallow pits, the manure would be less than one month old. After collecting manure, buckets full of manure were kept frozen at −20 °C (−4 °F) until they were used as feedstock for the reactors. Once manure was needed, one of these big buckets was thawed down and separated into a small bucket (usually 4 L (L) in volume). Just one of these small buckets was in the refrigerator for feeding; the rest were kept in a freezer to keep it unchanged as much as possible. Total and volatile solids (TS and VS) of each big bucket were tested to evaluate the proper feeding ratio. The total solid (TS) of the solid manure was 25.89%, and the volatile solid (VS) was 82.02% of the TS. There was no test conducted to verify the potential effect of the freezing, although there were little observed changes in biogas production between refrigerated and frozen manure in the last year of AD tests.
The inoculum was collected from semi-continuous AD jars of previous tests (Wang et al. [
45]), which were steadily producing biogas for over three months, and the feedstock was swine manure with organic loading rate (OLR) of 1 g-VS L
−1 d
−1 only. The total solid (TS) of inoculum was 2.20%, and the volatile solid (VS) was 64.92% of TS.
2.2. Experimental Design
Tests were carried out with laboratory size jars as reactors (adjusted for AD). Antibiotics were added to reactors with different concentrations of CTC and Tylosin to monitor antibiotic removal and gas production variations in those mesophilic reactors. Nine laboratory-scale jars as anaerobic bioreactors with the working volume of 1.375 L were kept at 39 ± 2 °C (102 °F to 105 °F) in the incubator. The jars are being fed with VA-free swine manure at 1g-VS per L-day, with 21 days hydraulic retention time (HRT). The volume of the feed given every two days is measured based on HRT and our reactor volume. Because our HRT is 21 days and the reactor volume is 1.375 L, so 0.131 L of our reactor liquid were removed and replace by feedstock (the digesters were fed every two days) [
46].
Each jar was connected to 10-L Tedlar bags to collect produced biogas, and the volume was measured every four days [
47]. A custom-built device was used to help distribute the biogas evenly in the bag, so the height of the bag could be measured more accurately. By utilizing a predetermined model, the volume of each bag was then estimated by bag height. Besides, to prevent any leakage of the Tedlar bags, each time two of the bags were randomly tested for possible leakage before emptying. Additionally, tubes, caps, and any connective parts were tested for leakage. After biogas measurements, bags were emptied safely and burned.
The experiment consisted of nine jars; three of them were spiked with CTC, three with Tylosin and the last three with both CTC and Tylosin, to observe the combined effect or any interaction between two types of antibiotics (
Figure 1). Furthermore, to investigate the efficiency of the AD process, we added six more jars, filled with distilled water and the headspace with N
2. Three of these jars were being kept in the incubator at the same temperature of the digester jars (39 °C), the rest were being kept in the room temperature to monitor the effect of the temperature. The same pattern of antibiotics concentration was conducted for control jars, two groups of three jars. Retention time and sampling procedures were identical.
2.3. Antibiotics
The two most widely used antibiotics were selected based on consumption rate and market share of different antibiotic classes in the United States [
4]. Antibiotics used in this experiment were CTC as chlortetracycline HCI and tylosin as tylosin tartrate. Commercial grade CTC was bought from “PharmGate Animal Health; Omaha, Nebraska” with the brand of “Pennchlor 64”. Commercial grade tylosin used was from “Elanco Animal Health; Indianapolis, Indiana” with the brand of “Elanco”. Moreover, to prepare standard samples for LCMS/MS, both antibiotics were ordered as the analytical grade from “Sigmaaldrich”, St. Louis, Missouri. Chlortetracycline hydrochloride, VETRANAL™, analytical standard, with CAS number of 64-72-2 and tylosin, United States Pharmacopeia (USP) Reference Standard, with CAS number of 1401-69-0 were used for standard solutions.
We followed the prescription on the labels to imitate the real condition in a barn. Consequently, the recommended dosage for tylosin was 66 ppm in drinking water. For Swine Dysentery, adding tylosin to drinking water should be continued for 3 to 10 days, depending upon the severity of the infection. For CTC, the recommended dosage was 22 mg/kg body weights per day. The duration of treatment is 3 to 5 days depending on the infection. Pigs are generating 4.28 L of manure per day on average. Additionally, it is assumed that the average body weight of a pig is around 68 kg (finishing pigs weigh around 45–113 kg). These assumptions would help us estimate the concentration we should inject in our reactors, by considering the excretion rate and metabolism percentage.
The stabilization time for the reactors and the microbial community was expected to be one to three months, until the biogas production, digestate pH, and alkalinity trend became flat. For the current research, the digester was fed for eight weeks or 2.5 times the retention time. Important operating variables, including organic loading rate based on total volatile solids (TVS), solid content, temperature, mixing (swirling the jar daily), and foaming (if any) were recorded. Digester alkalinity and pH were monitored weekly by measuring the digestate.
For Lower range concentration, the lowest factor in each section was used. For instance, to calculate the lower band of tylosin, 11.35 L per day as pig’s drinking volume, 50% excretion level applied and 4.28-L excretion per day was selected. For upper range concentration, the highest factor in each section was used. For instance, to calculate the upper band of tylosin, 18.93 L per day as pig’s drinking volume, 90% excretion level and 4.28-L excretion per day was selected. The average concentration is the average of the lower and upper concentration. Recalling that jars were loaded every two days with a mixture of solid manure and water, VAs added with feed had a concentration of day 1 plus day 2.
Since antibiotics are being added to the water part of the feeding (not to solid part), therefore solubility of the VAs should be checked.
Table 1 is a summary of the solubility of CTC and tylosin in the water at 20 °C:
Considering the solubility of CTC and tylosin in the water at 20 °C, there was no problem with CTC and tylosin solving limit individually. However, one set of three jars was used, which we decided to use to test the combined effect of antibiotics, so we had to mix two antibiotics in the same volume of water (0.103 L). There is always a chance of interference between two types of chemicals, especially when they are being added near their solubility limit. Thus, the decision was made to add CTC directly to the water, transfer it to the reactor and then add tylosin powder separately to the reactor. Other solvents such as methanol or ACN were dismissed because of their adverse effect and interference with the reactor’s performance (an independent test was conducted to evaluate the impact of adding methanol onto AD performance; details are not included in this paper).
Table 2 summarizes the recommended VA concentrations based on the manufacturer’s recommendation and corresponding dosages considering the ranges of dosage, water consumption, and excretion rate.
As previously mentioned, we recorded biogas production for at least two HRTs, before and after introducing the antibiotics.
Table 3 illustrates the added VAs concentration in each reactor. As shown below, the first group is being administrated only with CTC, the second group with both CTC and tylosin and the last, with only tylosin.
The pH of the digestate was measured every two days while adding antibiotic, with pH meter (PINPOINT, American Marine Inc., Ridgefield, CT, USA). Using pH data, microbial activity of the digester and the reactor performance is projected. However, pH can also be affected by alkalinity. For quality assurance, alkalinity tests were also carried out.
The CO
2 concentration of the biogas was measured with a standard combustion analyzer (Bacharach Fyrite Classic Combustion Analyzer, New Kensington, PA, USA) every eight days. The concentration of CO
2 was measured every four days during the antibiotic addition period. Comparative tests using a gas spectrometry were used to check how accurate our measurements were. Below is a comparative table that illustrates accuracy control values (
Table 4). The gas chromatograph device was (GC-2014, Shimadzu, US) with a thermal conductivity detector (TCD) using a ShinCarbon ST 80/100 Column (Restek, US) [
48].
Because the administration times for CTC and tylosin were different (6 days for CTC and 10 days for tylosin), jars with CTC spikes were sampled at the end of day 6, while tylosin-spiked jars were sampled at day 10. For jars with the combined CTC and tylosin, samples were taken at both day 6 and day 10. Since VAs were added every two days with feed, the sampling would occur two days after the last spike. Samples were frozen at −20 °C immediately. Gas production, pH and CO2 level were considered the vital data, which were recorded before, during and past spikes.
4. Results and Discussion
The presence of VAs in anaerobic digesters could have an inhibitory effect on biogas production [
36,
59,
60], because VAs could disrupt microorganisms’ dynamics, especially when the concentration is high. Because AD is not efficient in degrading VAs completely, in the long term, AD reactors can also become a fostering environment for VAs that would help the development of new ARGs [
61]. By scrutinizing the figures derived, some abnormalities were visible one week after the last spike, recalling that October 18 was the start date of the spiking antibiotics and final day was October 28 (
Figure 2a). This biogas fluctuation started with a decline in samples spiked with tylosin and also a mixture of tylosin and CTC, immediately after the first spike. For CTC samples, this drop was delayed until early November. On November 11, it grew again and then reached its lowest point on November 23.
For the CTC plus tylosin, after a drop on October 22, and again on October 30, we witness a surge after that. Reactor #4 peaks on November 7 and reactor #6 peaks on November 19. Reactor #5 climbs steadily during this period. They start to drop in mid-November and reach their low at the end of November. Similarly,
Figure 2b shows the same behavior, declining after the first spike until the end of October (last spike), followed by an upward trajectory. Likewise, this trend hits its bottom in early December.
Running a T-Test on biogas data implies that AD bacterial activity was immediately inhibited for samples that have tylosin in them (
Figure 2b,c) (
p-value = 0.005). Still, the bacteria either adapted or the inhibiting compound was removed from the system after a few weeks. Biogas production was untouched for CTC samples, yet for the mixture of CTC and tylosin, and tylosin alone, it was significantly lower, immediately after VA addition. The tylosin concentration in this experiment was 92 mg/L and less, complying with the findings of Mitchel et al. [
62]. They concluded that the bioreactor containing 92 mg/L tylosin had less biogas for nearly 30 d until the system recovered. The biogas reduction for samples with tylosin and CTC was close to 14%, and for tylosin, samples was between 8 and 19%, with no dose-dependent relationship. On the other hand, Chelliapan et al. [
63] found no biogas inhibition in an up-flow anaerobic stage reactor (UASR) containing 100–800 mg/L tylosin.
Erythromycin, another macrolide antibiotic caused 6–24% biogas reduction with 6–100 mg/L, and no dose-dependent relationship [
64].
But CTC did not disturb the bacterial activity, substantiate the evidence that CTC antibiotic may present minimal AD biogas inhibition at concentrations less than approximately 70 mg/L occurring in the current study. Yin et al. [
65] observed similar results; for a mesophilic anaerobic digester with the manure and CTC concentrations of 0, 20, 40, and 60 mg/kg. TS, no significant inhibition in biogas production occurred. Dreher et al. [
66] showed that no inhibition of biogas production happened in anaerobic sequencing batch reactor with 28 mg/L CTC but that the volumetric composition of methane decreased by about 13–15%. Mixed results of the inhibition in the literature could be due to various reactor types, inoculum/manure ratio, inoculum and manure age and source, reactor size, and batch or continuous operation [
67]. In this experiment, CTC concentration was probably lower than its required inhibitory level.
4.1. pH and CO2
pH value can demonstrate how well Acetogenesis and Methanogenesis bacteria are working. At the beginning of AD performance, Acetogenesis bacteria start to produce volatile acids that cause the pH to decrease. Subsequently, Methanogenesis bacteria convert the volatile acids to methane and CO2, and cause pH to increase. At HRTs with more than five days, the methane-forming bacteria begin to consume the volatile acids.
By comparing before and after the addition of VAs, it is evident that reactors are experiencing a fluctuating pH status (
Figure 3a–c). The graph shows that variations immediately after antibiotic spike have increased intensively, with a rising trend. Following up, in the first week of November, almost all reactors reach their plateau. From then on, the gradual decline continued until November 17. Subsequently, reactors seemed to recover themselves with an increase in pH. By the end of November, pH returns to its average level of around 7.8.
The T-test on the pH data shows that pH values were significantly lower for the samples with tylosin in them (p-value = 0.05). However, CTC did not affect the pH significantly. Since fluctuations in pH level are not sharp, this indicates that VAs did not disturb the bacterial community substantially.
Nevertheless, none of the reactors became upset or affected intensively by the VA addition. pH fluctuation was ±0.16 maximum and it never dropped under 7.60. Similarly, the CO
2 level has detectable alteration around the antibiotic spike date (
Figure 4a–c).
Biogas produced is consisting of almost 50–75% of methane, 25–40% of carbon dioxide and other gases, depending on organic material [
68]. By comparing CO
2 data and performing a T-test, results imply that CTC had a significant effect on the biogas methane content (
p-value = 0.05). At the same time, samples with tylosin only were not affected considerably. The reason could simply be that the CTC is active primarily against Gram-negative organisms by blending with the A location of the 30S subunit of bacterial ribosomes. So, they prevent peptide growth and the protein synthesis effect, which finally leads to bacteria death [
65,
69]. Methanogen bacteria are Gram-negative bacteria [
70]. Thus, at a certain level of CTC, significant biogas inhibition should be imposed to the bioreactor.
Values of VS in the digestate before and after VA addition are shown in
Figure 5. In general, every treatment sample except for the medium tylosin concentration, showed an increase in the VS percentage after VA addition. The VS values agree with the slight fluctuations observed in
Figure 2,
Figure 3 and
Figure 4, that the microbial communities were slightly affected by the VA addition and the biogas production was not halted. To recall, the initial manure VS loading was 2 g-VS/L/day and the sampling for VS were conducted between 6 and 10 days after the first injection and two days after the last injection of VA. Thus, it complies with reduction in methane production which was around 13–15% studied by Dreher et al. [
66]. Angenent et al. [
34] also reported a temporary decrease in VS removal which recovered quickly. The average VS level before and after VA addition is 0.56% and 0.64%; and VS removal is 1.44% and 1.36%, respectively.
4.2. LC-MS/MS Results and Adjustments
The plan was to try using the measured Enrofloxacin concentrations to calculate adjustment factors for the other VAs. With these factors, sample concentrations after the dilution and short-term loading in the lab digesters were recalculated, assuming there was no degradation or absorption. Should there be significant differences, these would then be caused by sampling error, degradation due to AD, or the error of the LC-MS/MS measurement, including the SPE.
Table 6 shows all of the samples, their added VAs and a comparison between spiked concentration, detected concentration by LC-MS/MS and recalculated concentration using adjustment factors.
4.3. Relatively High Recovery of Enrofloxacin in the Water-Only Samples
Although the spiked Enrofloxacin in digestate and manure samples had a very low recovery rate (226 ppb to 1433 ppb vs. 4444 ppb spiked values) (
Table 7), all but one water sample detected relatively higher Enrofloxacin concentrations (1667 ppb to 3197 ppb,
Table 1). The Enrofloxacin concentration in the first water sample (236 ppb, samples #1) was only a fraction of water samples. The water samples were made with distilled water and VAs, no solid manure. Relatively higher recovery rates suggest that there is a systematic bias in measuring the Enrofloxacin in the samples that have solids (manure and digestate).
Therefore, when sample 13 was excluded, the average of the water sample group was 2481 ppb, while the digestate samples averaged 825 ppb. On the other hand, if we disregard the Enrofloxacin concentrations as an adjustment factor and just compare the LC-MS/MS values with our calculated concentrations (assuming no degradation), provides a better outcome. In this way, external standards are utilized to evaluate samples with only water and VAs, to monitor whether the removal of VAs is due to AD or not.
4.4. Very High Recovery of Tylosin in the Water-Only Samples
As shown in
Table 8, concentrations for tylosin are very close to and sometimes higher than what we were expecting (LC-MS/MS measured 106 ppm, we expected 92 ppm for sample 2). Furthermore, the water samples that were not spiked with tylosin did yield very low tylosin concentrations (0.86 ppm and 1.79 ppm, samples 13 and 16).
It is a different case for CTC; yet the reasons for low CTC detection are still unknown. The trend for CTC concentration shows they are disappearing so fast, which may be due to its half-life degradation or anaerobic reactor removal; alternatively, this may simply be because the CTC we used was already degraded, see the discussion below.
4.5. Consistent and Proportional LC-MS/MS Tylosin Results in the Digestate Samples
Based on the tylosin results being more consistent than CTC and Enrofloxacin results, the LC-MS/MS results of the digestate samples were meticulously scrutinized. Even though the LC-MS/MS detected concentration values of tylosin that were lower than expected, they were consistently proportional to the concentrations. For example, the expected concentration of tylosin was 20, 40, and 60 ppm for samples 4, 5, and 6, with zero degradation assumption, the LC-MS/MS values were 0.5, 0.9, and 1.2 ppm. No significant correlation for CTC was found.
Figure 6 shows the LC-MS/MS results (
Y-axis, ppm) vs. spiked values (
X-axis, ppm). The consistently lower measured concentrations in the digestate samples and the high recovery rates in the water samples suggest that there was significant tylosin degradation due to the AD process.
4.6. Consistent LC-MS/MS Tylosin Measurements in the Manure External Standards
The detected tylosin concentrations of the three external standard samples were similar and had low deviation,
Table 9. For tylosin, the detected levels ranged from 39.8 to 44.3 ppm and averaged 41.5 ppm, while the expected concentration was 77.9 ppm. For CTC, the measured concentrations were again a small fraction of the expected level.
4.7. Relatively Higher Recovery of CTC in the External Standard Samples, and Two Additional CTC Standards
Studying the results of the three external standard samples showed that there was a relatively higher recovery rate for CTC. As an instance, compared with the concentration of 233.78 ppm, LC-MS/MS detected 14.96, 36.16 and 28.52 ppm. Compared with previous CTC samples, this group had a much higher recovery rate, which was fresh samples made with diluted manure and VAs. Besides, the external standards were prepared with commercial-grade antibiotics instead of analytical grade.
Table 10 shows results for freshly prepared samples with diluted manure and CTC antibiotic, at concentrations of 4ppm and 40 ppm, and the LC-MS/MS measured concentrations were 1.63 ppm and 11.7 ppm, respectively. The recovery rates of CTC were 29% and 41%. Since the samples were freshly prepared, the probability of degradation due to AD or half-life degradation was eliminated. Other possibilities are absorption to organic matter, and that inconsistent purity or degradation had already happened before application. For tylosin, LC-MS/MS detected higher concentrations (12 ppm and 67.6 ppm detected for 4 ppm and 40 ppm samples, respectively), which gives a detection rate of 300% and 169%.
4.8. Applying External Standard Correction Factors
Because the internal standards (Enrofloxacin) did not yield consistent measurement, the correction was made based on external standards instead. By applying the external standard adjustment factor, the VA concentrations were corrected accordingly. The adjustment factor was obtained from samples with manure and spiked antibiotics, without retention time for AD. In other words, we just spiked different concentrations of antibiotics in samples made with manure, then prepared those for LC-MS/MS, immediately. In this way, we may be able to track other important factors contributing to our results, such as absorption, ion suppression or enhancement and recovery rate.
Figure 7 presents the measured and corrected VA concentrations.
By using the external standard correction factor instead of the Enrofloxacin correction factor, data are more consistent, especially for tylosin (less than 6% error). It suggests that Enrofloxacin failed to act as an ideal internal standard. The inconsistency could be due to Enrofloxacin binding to the abundant organic materials or the presence of Ca2+, Mg2+ Ions.
4.9. Degradation of CTC
For CTC, results showed a high degradation rate for both the samples in water and AD (
Figure 8). For instance, almost all CTC injections with various concentrations have close to zero concentration. The low concentrations were measured for AD-treated samples, and also for CTC dissolved in water stored at room temperature and 40 °C, the temperature of the AD. In addition, the concentrations of the external standards were 234 ppm. The results suggest that the CTC degrade much faster than the tylosin, which might be due to the shorter half-life (8 days) as reported by the manufacturer. It is also possible that the CTC powder we used had already degraded. CTC concentration in external standard samples was reduced to 131, 320 and 251 ppm from its original 234 ppm. Because the purity correction was already applied and recovery rate adjustment was made, also only low CTC concentration was detected for the water and digestate samples, the CTC probably just degraded itself over a short time. A study by Winckler and Grafe [
71] showed that tetracycline in liquid manure was degraded by 50% in 82 days. Arikan [
72] reported a 75% reduction in CTC concentration with AD after 33 days, with a half-life of 18 days. Cheng et al. [
73] reported a high affinity between tetracycline and solid manure during AD. For future research, additional testing to examine the possibility of the adsorption by the glass jar used in this study as the AD reactors should be conducted, since there are very few investigations on this subject.
4.10. Degradation of Tylosin
Tylosin degraded very well with the ASBR reactor working at 39 °C and loading with swine manure every two days. By comparing the degradation rate of tylosin in ASBRs with jars filled only with water, it shows that AD is effective in reducing tylosin (
Figure 9). The degradation rate of tylosin in water averaged 33.5; however, the degradation was 100 percent with AD. A study by Kolz et al. [
74] concluded no effective degradation for tylosin B and D in anaerobic conditions up to eight months. tylosin A was degraded under the aerobic conditions with a half-life of 2 to 40 days [
5,
75]. Stone et al. [
59] also reported no significant degradation for tylosin.
4.11. Effects of Having Two Types of VA in the Digestion and Water
In these experiments, we planned to compare the effect of having two different antibiotics on the reduction efficiency of antibiotics. Identical injections and concentrations were applied in the mixture treatment, except for a series of samples, both antibiotics were spiked.
For tylosin, the results showed that there was little difference between samples. The tylosin removal was similar with or without CTC mixture, suggesting that the CTC addition had no adverse effect on tylosin degradation.
However, because CTC degradation was much faster than tylosin, and because the samples collected for CTC concentration measurements were not resolute enough (shorter time than the five-day sampling), the speed of the CTC degradation and effects of the tylosin addition could not be determined based on this dataset (
Table 9).
4.12. Contrast Tylosin Reduction in Water and AD Reactors
Tylosin tartrate showed a relatively higher removal in AD treatment when compared to samples that were dissolved in water.
Figure 10 depicts that the tylosin samples treated by AD for twelve days were removed entirely (100% removal). However, the removal of tylosin tartrate in water (dotted line) was 40% or less during the same period. The lower degradation in the water samples suggests that AD is effective in enhancing tylosin degradation in the animal manure, and that this could be the essential effect of the AD.
4.13. Bacterial Community Dynamics
Phyla Firmicutes, Bacteroidetes, Proteobacteria, and Synergistetes were dominant or co-dominant in bacteria. Different types of Clostridia consisted mostly of the Firmicutes. Methanomicrobia was the dominant Archaea among our samples.
The first group was treated with CTC and showed a slight fluctuation in
Archaea abundance, 4.37%, 5.19%, and 5.06% related to low, medium, and high concentrations of CTC, respectively. As shown in the
Table 11,
Firmicutes almost remained constant, and Bacteriodetes increased from 14.03 to 15.30% and then decreased to 12.50% for low, medium and high concentrations of CTC, respectively. The same pattern occurred for Synergistetes, going up from 1.67 to 1.75% and then dropping to 1.44%, with the mentioned level of CTC concentrations. However, the reverse happened for Proteobacteria: the abundance level reduced from 2.46 to 1.35% and then rose to 1.87%.
The second group, consisting of four samples, all are added with the mixture of CTC and tylosin with low, medium and high concentrations. Two samples were taken from medium range concentration because the administration duration of CTC was 6 days and tylosin was 10 days. Therefore, samples were taken at the end of the administration of each VA, at day 6 and day 10. There is a reverse relation between Archaea abundance and VAs concentration as well as administration duration. Archaea level is dropping with a higher concentration of VAs and longer retention time. The effect of the administration period is stronger than the dosage on the Archaea population. Firmicutes, on the other hand, have increased from 74.9 to 79.70% by increasing the dosage of VAs from low to medium. The trend is not consistent with shifting from medium to high concentration of VA; it would decrease the abundance of Firmicutes. The administration period has the same effect, but not as much as dosage. Phylum Bacteroides population increased with the rise of VAs concentration and doubled with an increase in retention time, from 6 to 10 days with a medium level of VAs. Proteobacteria and Synergistetes abundance both have dropped by increasing dosage for low to medium, but the fall is drastic for Proteobacteria, changing from 7 to 1.8%. By the end of CTC administration, when adding tylosin, their abundance recovered slightly. Surprisingly, in high concentrations of CTC and tylosin, both of these bacteria showed growth in their population. The proteobacteria population is almost doubled by having high dosage if VAs instead of low dosage.
The last group, which is being medicated with tylosin only, Archaea, decreased slightly and then almost doubled when moving from low concentration to high. For Bacteroides and Firmicutes, it is exactly the reverse, with maximum abundance around medium range concentration and the nearly same number for low and high concentrations. For proteobacteria, results showed a sharp drop with shifting from low to medium concentration, 10.2% changed to 4.03%. However, by raising the dosage, the abundance of proteobacteria returned to 15.97%.
Bacteroidia are the major classes found within the phylum of
Bacteroidetes; and are abundant in digesters that use cow manure as feedstock [
76].
Firmicutes phylum is mostly syntrophic bacteria that can decompose a variety of fatty acids, and exists in both activated sludge systems and anaerobic digesters [
77]. Within the species of
Firmicutes,
Clostridia is the dominant class. The predominance of
Clostridia in the AD sludge was related to the comparably fast hydrolysis and VFA (volatile fatty acids) fermentation happening in the digesters [
78].
Fatty acid-oxidizing bacteria, including Synergistales group which have syntrophs are connecting bonds of the chain between the primary fermenters and methanogens [
79], are abundant in thermophilic digesters [
80,
81]. The presence of
Synergistetes (syntrophic acetate oxidizers) might be an indicator of a decent acetotrophic activity in the bioreactor [
82].
There are two major categories of methanogens; acetoclastic which consumes acetate to produce methane or hydrogenotrophic that are converting CO2 and H2 to methane.
The
acetoclastic methanogenesis is linked with the
Methanosarcinales and the
hydrogenotrophic methanogenesis is linked with the
Methanomicrobiaceae family. In the current study, the
hydrogenotrophic pathway with
Methanosarcinaceae was dominant. Kim et al. [
42], Nogueira et al. [
78] and Padmasiri et al. [
83] also detected a dominant
Methanomicrobiales order on AD.
4.14. Statistical Analysis of the Effect of Different Treatments on Samples
There is a notable difference between the treated and control group, which shows antibiotics had a significant influence on altering the bacterial community in our digesters.
Figure 11 shows the samples arranged using the same two similarity measures used to generate the PCoA plots but in the form of a dendrogram. Bray–Curtis similarities and Jaccard similarities method are used for this comparison.
In both of the methods, the differences between all three treatment groups are modest and likely obscured by the variability introduced by the control samples.
Figure 12 shows a stacked bar chart at the best taxonomic resolution afforded by our primers. Again, the differences between the two datasets are stark, while the differences between treatment groups are more subtle (but present).
Figure 13 shows a heat map in which samples (columns) are ordered according to similarity using a hierarchical method (UPGMA) based on the 50 OTUs (operational taxonomic unit) (rows) with the lowest p values following serial ANOVA testing of all 629 OTUs. In short, it shows perfect discrimination between groups when the samples are clustered using only the 50 most variable OTUs. Taxonomic identity of the microbes is listed on the right-hand side of the heat map.
Figure 14 shows box plots representing the relative abundance of eight of the OTUs with the lowest p values. One can easily note the very clear pattern of microbes with sensitivity to one or the other drugs, with Proteiniphilum (lower left) being the anomaly. Some other microbes with low p values mostly had these types of patterns, either down in CTC and CTC + tylosin, or down in tylosin and CTC + tylosin.
Finally,
Figure 15 is a random forest analysis looking for “biomarkers” of each treatment group. The greater the “MeanDecreaseAccuracy”, the better that OTU is as a biomarker of the rankings shown to the right of the Figure. For example, Methanoculleus is apparently an excellent predictor of these groups by having a high relative abundance in the CTC samples and low abundance in the CTC + tylosin samples. Likewise, the Anaerorhabdus furcosa group and Flexilinea sp. can be found more on samples with high tylosin concentration and low CTC. Ruminiclostridium sp. can be abundant in conditions with high CTC levels and low tylosin.