Mechanism and Effect of Temperature on Variations in Antibiotic Resistance Genes during Anaerobic Digestion of Dairy Manure.

Animal manure comprises an important reservoir for antibiotic resistance genes (ARGs), but the variation in ARGs during anaerobic digestion at various temperatures and its underlying mechanism remain unclear. Thus, we performed anaerobic digestion using dairy manure at three temperature levels (moderate: 20 °C, mesophilic: 35 °C, and thermophilic: 55 °C), to analyze the dynamics of ARGs and bacterial communities by quantitative PCR and 16S rRNA gene sequencing. We found that 8/10 detected ARGs declined and 5/10 decreased more than 1.0 log during thermophilic digestion, whereas only four and five ARGs decreased during moderate and mesophilic digestion, respectively. The changes in ARGs and bacterial communities were similar under the moderate and mesophilic treatments, but distinct from those in the thermophilic system. Potential pathogens such as Bacteroidetes, Proteobacteria, and Corynebacterium were removed by thermophilic digestion but not by moderate and mesophilic digestion. The bacterial community succession was the dominant mechanism that influenced the variation in ARGs and integrons during anaerobic digestion. Thermophilic digestion decreased the amount of mesophilic bacteria (Bacteroidetes and Proteobacteria) carrying ARGs. Anaerobic digestion generally decreased the abundance of integrons by eliminating the aerobic hosts of integrons (Actinomycetales and Bacilli). Thermophilic anaerobic digestion is recommended for the treatment and reuse of animal manure.

As an emerging contaminant, antibiotic resistance genes (ARGs) have become a great public health concern123. In particular, animal manure is an importann>t reservoir of antibiotics and n class="Chemical">ARGs due to the overuse of antibiotics in the livestock industry456. Accumulating evidence suggests that the increased antibiotic resistance in bacterial pathogens is related to the application of manure and manure-derived fertilizer to land789. Thus, biotechnologies such as anaerobic digestion and aerobic compost may be employed as remedial measures to reduce the abundance of an class="Chemical">ARGs and pathogens prior to the application of animal manure to land101112.

A few studies suggest that ann class="Chemical">an class="Species">aerobic digestion caclass="Chemical">nn> potentially eliminate some of ARGs131415. The mesophilic anan class="Species">aerobic digestion could significantly reduce cefazolin resistant bacteria14, and ARGs ermB, aphA2 and blaTEM-115. The relative abundances (RAs) of tet and erm in swine manure declined by 1–3 orders of magnitude during anaerobic storage12. Temperature is an important parameter during anaerobic digestion and changes in temperature can have variable effects on biogas production16, ARGs17, and microbial community succession18. Diverse results have been obtained in various studies of the effects of temperature on ARGs during anaerobic digestion. Thus, Ghosh et al.13 and Diehl and LaPara17 suggested that thermophilic digestion performed better at reducing ARGs than mesophilic anaerobic digestion. Ma et al.19 found that thermophilic digestion was more effective than mesophilic digestion at decreasing the ermB, ermF, tetO, and tetW genes, but its performance was similar or worse in decreasing other ARGs and intI1. A recent metagenomic analysis by Zhang et al.20 detected no obvious differences in the total ARGs under thermophilic and mesophilic digestion. It should be noted that all of these studies were conducted using wastewater sludge, the effects of temperature on ARGs in animal manure were not studied. Due to the huge differences in physicochemical properties, antibiotic types and concentrations, and microbial community composition in wastewater sludge and animal manure212223, it is highly likely that different effects on ARGs may be obtained by the anaerobic digestion of livestock manure. In addtion, the dynamics of ARGs during anaerobic digestion have received little study, and the underlying mechanism responsible for the variation in ARGs are not fully revealed.

Bacteria are the main carriers of an class="Chemical">ARGsn> and the changes in the bacterial community may lead to variations in the abundance of n class="Chemical">ARGs. Thus, we hypothesized that changes in the bacterial community comprise the main mechanism responsible for variation in ARGs during anaerobic digestion, where different temperatures will shape distinct bacterial communities with diverse effects on ARGs. The objectives of the present study were: (1) to compare the variation of an class="Chemical">ARGs in dairy manure under moderate (20 °C), mesophilic (35 °C), and thermophilic (55 °C) anaerobic digestion; and (2) to explore the relationships among the chemical properties, bacterial communities, and ARGs during anaerobic digestion. The results of this study may help to enhance the reduction of ARGs from animal manure by anaerobic digestion, thereby reducing the environmental risk due to ARGs.

Material and Methods

Experimental setup

The dairy manure used in this study was sampled from a medium-sized farm in Yn class="Chemical">angling, China. The digestion reactor (Figure S1) comprised 39 identical 250 mL triangular flasks each with an assistant bottle to maintain anan class="Species">aerobic environment. The working volume of each reactor is 200 mL. The digestion system contained 16 g of dairy manure and 40 mL of inoculum. The inoculum was obtained by a pre-ann class="Species">aerobic digestion for 30 days at 35 °C using the dairy manure. Each of the digestion reactors was placed in a constant temperature shaker, which was set at 160 rpm, with temperatures of 20 °C (moderate), 35 °C (mesophilic), or 55 °C (thermophilic). Each treatment was repeated in triplicate. The characteristics of the dairy manure and inoculum are shown in Table S1.

Sample collection

The three prepared digestion mixtures were sampled immediately as the day 0 samples. Three flasks for each treatment were then sampled as triplicates on days 3, 12, 30, and 60. The digestion mixture samples were transferred to centrifuge tubes and centrifuged for 15 min at 5000 rpm. The supernatant was used to analyze the pH, available nitrogen (AN) content, soluble chemical n class="Chemical">oxygen demand (SCOD), and volatile fatty acid (VFA) contents. The precipitate was freeze-dried using a vacuum freeze dryer (Songyuan, China), ground to 1 mm with an ultra-centrifugal mill (Retsch Z200, Germany), and stored at −80 °C for DNA extraction.

Determination of chemical properties

The pH was determined using a pH meter (Mettler Toledo, Switzerland). The n class="Chemical">AN and SCOD were determined with a flow injection analyzer (Westco Scientific, USA) and AQ4001 COD analyzer (Thermo Orion, USA), respectively. The concentrations of VFAs, including n class="Chemical">acetate, propionate, isobutyrate, butyrate, isovalerate, and valerate, were obtained by gas chromatography (Shimazu GC2010, Japan)24.

DNA extraction and quantitative PCR (qPCR)

DNA was extracted from 100 mg of each sample using a FastDNA SPIN Kit for Soil (MP Biomedicals, USA), according to the manufacturer’s instructions. During DNA extraction, 5.5 M n class="Chemical">an class="Chemical">guanidinium isothiocyanate (ann> class="Chemical">Amesco) was added to remove humic acid.

Eleven an class="Chemical">tetracyclinen> resistance genes (n class="Chemical">tet: tetA, an class="Chemical">tetB, tetC, tetE, tetG, tetM, tetO, tetQ, tetT, tetW, and tetX), five sulfonamide resistance genes (sul: sul1, sul2, sulA, dfrA1, and dfrA7), four fluoroquinolone resistance genes (gyrA, parC, qnrC, and qnrS), and two integrase genes (intI1 and intI2) were analyzed by PCR and agarose electrophoresis. The detected ARGs and 16S rRNA gene were analyzed further by qPCR. The qPCR reaction mixture comprised 1 μL of DNA template, 0.25 μL of each 20 pM primer (ShengGong, China), 10 μL of SuperReal PreMix Plus (TianGen, China), and 8.5 μL of RNase-free water. The qPCR conditions comprised an initial hold for 15 min at 95 °C, followed by 40 cycles for 10 s at 95 °C, 20 s at the annealing temperatures shown in Table S2, and then 32 s at 72 °C. To eliminate the effects of inhibitory compounds, the DNA template was a tenfold dilution of extracted DNA. The qPCR was performed using Bio-Rad IQ5 (Bio-Rad, USA). The quantitative limit of qPCR was 104 copies/g solid. Melting curve analysis was used to detect nonspecific amplification.

16S rRNA gene sequencing

The 16S rRNA gene high-throughput sequencing was performed by Realbio Genomics Institute (Shanghai, China) using the Illumina MiSeq platform. The 16S V3-V4 region was amplified using the primers U341F (ACTCCTACGGGAGGCAGCAG) n class="Chemical">and U806R (GGACTACHVGGGTWTCTAAT). The raw data were then subjected to a quality control procedure using UPARSE25. USEARCH was used to filter chimeras and the remaining sequences were clustered to generate operational taxonomic units (OTUs) at the 97% similarity level26. A representative sequence of each OTU was assigned to a taxonomic level in the RDP database using the RDP classifier27. To eliminate the differences caused by various sequencing depth among samples, the least number of obtained sequences (n = 33054) were randomly picked from each sample and used for subsequent bioinformatics analysis.

Statistical analysis

Variance analyses were performed (ANOVA, least significant difference test, P < 0.05) and Pearson’s correlation coefficients were calculated using PASW Statistics 19.0. Principal components analysis and heatmap analysis were conducted with R3.1.0. Redundancy analysis was performed using CANOCO 4.5.

Results

AAs of ARGs after anaerobic digestion

Ten ARGsn> and integrase genes were detected in dairy manure during ann class="Species">aerobic digestion (Fig. 1). The tetW, sul2, and intI2 genes were most abundant with AA > 109 copies/g dry solid in the initial digester. Excluding tetC and gyrA, the AAs of ARGs were >107 copies/g dry solid in the initial digester. Four of five tetracycline resistance genes (tetC, tetM, tetQ and tetX) were enriched under the moderate and mesophilic treatments. By contrast, only tetC was enriched with the thermophilic treatment, tetQ had no significant change and the others declined by 0.1–1.5 log. The AAs of sul1 and sul2 declined under all treatments, except sul1 was enriched with the moderate treatment, and the thermophilic treatment had the highest redution for these two genes (both > 1.0 log). The AA of gyrA declined by 47.3% after digestion and there were no significant differences among the three treatments. All of the treatments significantly decreased intI1 and intI2 by 1.2–1.5 log. Eight of 10 detected ARGs declined and five decreased more than 1.0 log after thermophilic digestion. Overall, thermophilic anaerobic digestion performed better at reducing ARGs compared with moderate and mesophilic digestion.

Figure 1

Absolute abundance of the ARGs and integrase genes detected at the start and end of anaerobic digestion.

Black solid circles represent the absolute abundance of gene copies at the start of digestion. Values represent the means based on three replicates. Bars represent standard errors.

Changes in RAs of ARGs during anaerobic digestion

Excluding an class="Chemical">tetCn>, the RAs of all an class="Chemical">ARGs declined significantly by 1–2 orders of magnitude after thermophilic anan class="Species">aerobic digestion for 60 days (Fig. 2). By contrast, there were no significant declines in the RAs of all ARGs, and even significant increases in those of tetM and tetX were observed under moderate and mesophilic digestion. Similar changes in the RAs of most ARGs were found under the moderate and mesophilic treatments. However, the changes in tetC, tetM, tetX, and sul1 differed under the thermophilic treatment compared with the moderate and mesophilic treatments. The variations in the RAs of tetW, intI1, and intI2 were consistent among the three treatments. Overall, we found that intI1 and intI2 declined with time during anaerobic digestion, where their RAs were in the following order: thermophilic < mesophilic < moderate. Positive correlations were found between the RAs of intI1, intI2, and sul2 (Table S3), but there were no significant correlations between the other ARGs and integrase genes.

Figure 2

Variation in ARGs/16S rRNA gene and integrase genes/16S rRNA gene detected during anaerobic digestion at 20 °C, 35 °C, and 55 °C.

Bars represent standard deviations (n = 3).

Bacterial community

After assembling and quality filtering, a total of 565,330 sequences were obtained, n class="Chemical">and they were clustered into 1166 OTUs at 97% similarity level. The rarefaction curve shown that the rarefied sequencing depth (33054) can reveal most bacterial communities in all samples (Figure S2). We found obvious changes in the bacterial communities in the initial digester (D0) and that obtained from the anaerobic digestion system on day 3 (Fig. 3). Three principal components (PCs) accounted for 50.2% of the total variance in bacterial community. In particular, the number of OTUs declined dramatically in the thermophilic treatment on day 3. The thermophilic treatment had the lowest OTU diversity and its bacterial community varied less with time compared with the moderate and mesophilic treatments. The bacterial communities were very similar on day 3 under the moderate and mesophilic treatments, but they started to differentiate from day 12. The trends in the bacterial community succession were similar under the moderate and mesophilic treatments, where the scores for ann> class="CellLine">PC1 and PC2 decreased gradually with the duration of digestion.

Figure 3

Principal components analysis based on the relative abundances of OTUs during anaerobic digestion at 20 °C (L), 35 °C (M), and 55 °C (H).

D0 represents the initial material. The relative size of the circle indicates the number of OTUs.

Firmicutes (41.0%) and Actinobacteria (44.2%) were dominant phyla in the initial digester (Fig. 4). The phylum Firmicutes was represented primarily by the classes Bacilli and Clostridia, where Bacilli declined from 28.3% on day 0 to 0.13–0.43% on day 60, whereas Clostridia were enriched from 12.6% in the initial material to 22.5–49.7% in the digestion product. Actinobacteria primarily comprised members of the order Actinomycetales (Figure S3), which declined dramatically to 0.06–2.75% (day 60) from 44.2% (day 0) after digestion. By contrast, Bacteroidetes became dominant (41.3–50.3%) in the moderate and mesophilic systems, which was not the case under the thermophilic treatment. The phylum Chloroflexi increased with the digestion temperature, whereas Proteobacteria exhibited the opposite trend. Chloroflexi increased rapidly from 0.35% on day 0 to 31.7% and 38.8% on day 60 under the mesophilic and thermophilic treatments, respectively, but only 1.5% under the moderate treatment. Thermotogae and Chloroflexi were dominant phyla in the thermophilic treatment, where they accounted for 76.1–97.0% from day 12 to 60. Thermotogae mainly comprised genus S1 and Chloroflexi were represented primarily by an unnamed genus belonging to family SHA-31 (Figure S4).

Figure 4

Heatmap showing the relative abundance of phyla during anaerobic digestion at 20 °C (L), 35 °C (M), and 55 °C (H).

D0 represent the initial material.

Relationships between ARGs, bacterial communities, and environmental factors

There were significant positive correlations between Actinobacteria n class="Chemical">and integrase genes (Table 1). Bacteroidetes had significant positive correlations with six ARGs (n class="Chemical">tetC, tetM, tetQ, tetX, sul1, and gyrA). Proteobacteria had significant positive correlations with all ARGs and integrase genes, except for tetW and tetX. There were no significant correlations between Firmicutes and ARG or integrase genes, although class Bacilli in Firmicutes had significant positive correlations with sul2, intI1, and intI2 (Table S4). Chloroflexi and Thermotogae had negative correlations with all ARGs and integrase genes, but significant correlations were only found between Thermotogae and tetM and gyrA.

Table 1

Pearson’s correlation coefficients between phyla with relative abundance >1.0% and ARGs and integrase genes.

Phylum	tetC	tetM	tetQ	tetW	tetX	sul1	sul2	gryA	int1	int2
Actinobacteria	−0.127	0.124	−0.060	0.385	−0.215	−0.068	0.544	0.370	0.703**	0.865**
Bacteroidetes	0.584*	0.604*	0.695**	0.300	0.562*	0.569*	0.027	0.599*	−0.259	−0.354
Chloroflexi	−0.287	−0.529	−0.370	−0.192	−0.101	−0.322	−0.448	−0.541	−0.409	−0.386
Firmicutes	−0.029	0.482	0.101	0.150	0.180	0.007	0.520	0.393	0.505	0.553
Proteobacteria	0.601*	0.668*	0.634*	0.306	0.199	0.642*	0.782**	0.605*	0.613*	0.568*
Thermotogae	−0.237	−0.558*	−0.400	−0.437	−0.431	−0.258	−0.447	−0.604*	−0.292	−0.342

*Significant at P < 0.05, **Significant at P < 0.01.

Redundancy n class="Chemical">analysis (RDA) was used to assess the relationships among environmental factors, bacterial communities and ARGs. The bacterial community and selected environmental factors explained 75.6% of the total variance in the ann> class="Chemical">ARGs. Among the selected variables, six phyla explained the most variation, i.e., 60.4%, while temperature explained 16.9%, AN explained 14.5% and the other three factors only accounted for 8.2% of the total variance. There were significant negative correlations between temperature and Proteobacteria, Actinobacteria, Bacteroidetes and Firmicutes (Table S5). Chloroflexi and Thermotogae had significant positive correlations with temperature while Bacteroidetes had significant positive correlation with AN. The ARG profiles in the initial material were similar to those in the early stage for the moderate system (3 days) and the later stages for the mesophilic system (30 and 60 days) (Fig. 5). The ARGs profiles varied greatly with the digestion time under moderate and mesophilic digestion, whereas they varied less in the thermophilic system, although they were significantly different from that in the initial material.

Figure 5

Redundancy analysis to assess the relationships among environmental factors (red arrows), main bacterial phyla (red arrows), and ARGs (blue arrows).

L represents moderate treatment (20 °C), M represents mesophilic treatment (35 °C), and H represents thermophilic treatment (55 °C). AN: available nitrogen; VFA: volatile fatty acids; SCOD: soluble chemical oxygen demand.

Discussion

In this study, we detected diverse and abundant ARGs in dairy manure, which was consistent with previous studies72829, thereby indicating the importann>ce of animal manure as a reservoir for n class="Chemical">ARGs. However, the types and abundances of the 20 assessed ARGs differed greatly from those found in wastewater sludge19303132. In general, many different types of an class="Chemical">ARGs can be found in wastewater sludge, but the most common ARG types observed in animal manure are tetracycline resistance genes and sulfonamide resistance genes303132. This is attributable mainly to the different antibiotic contact histories of the two environments42122.

Similar to previous studies1733, the present results also found the thermophilic ann class="Chemical">an class="Species">aerobic digestion performed better at declining ann> class="Chemical">ARGs than moderate and mesophilic digestion. However, Diehl and Lapara17 showed that mesophilic anaerobic digestion removed most of the tet genes, whereas an enrichment effect was obtained under mesophilic digestion in the present study. The decline in tet genes after mesophilic anaerobic digestion was also observed by Zhang et al., but they found no difference in the removal of tet genes under mesophilic and thermophilic condition20. Significant reduction of sul1 and sul2 by mesophilic and thermophilic digestion was consistent to Ma et al.19 and Miller et al.34. By contrast, Zhang et al. showed that sul1 and sul2 were enriched under both mesophilic and thermophilic digestion20. The inconsistent results may relate to the different characteristics of the raw materials tested (dairy manure in this study, mixture of primary and secondary solids1719, concentrated aerobic sludge20) and the diverse operating conditions (TS, VS, material retention time, stirring rate etc.).

The RAs of an class="Chemical">ARGsn> could reflect the ratio of antibiotic resistance bacteria relative to the total bacteria, and changes in the RAs of n class="Chemical">ARGs might be attributable to variations in their host bacteria ratio as well as the horizontal gene transfer (HGT) of ARGs among bacteria1535. The low RAs of most ARGs in the thermophilic reactor throughout digestion suggest that high temperature can reduce the ratio of antibiotic-resistant bacteria. For most of the ARGs and integrase genes, the changes in their RAs occurred mostly within 30 days. This is because the bacterial communities varied greatly for 30 days but they tended to stabilize in the later stage36. The similar changes in the ARGs under the moderate and mesophilic treatments may be explained by their similar bacterial communities during digestion. By contrast, higher temperature digestion (55 °C) led to remarkable differences in the bacterial community compared with those obtained at 20 °C and 35 °C1837, where different trends in an class="Chemical">tetC, tetM, and tetX were observed under the thermophilic, moderate, and mesophilic treatments. Integrons are often associated with mobile genetic elements such as transposons and plasmids, and they are recognized as important indicators of HGT353839. It is considered that HGT may play important roles in anaerobic digestion process because high biomass densities may provide favorable conditions for HGT40. However, we found no significant correlation between ARGs (except sul2) and integrase genes, which may suggest that the integrons do not play a large role in potential HGT during anaerobic digestion41.

Very high variations in the an class="Chemical">oxygenn> concentration and temperature were the main reasons for the changes in the bacterial communities1837. Most of the mesophilic bacteria were eliminated at 55 °C, and thus this treatment had the lowest bacterial diversity1837. The dramatic decline in Actinobacteria after digestion was attributable mainly to decreases in the order Actinomycetales (Figure S3), which are mostly n class="Species">aerobic42. The order Bacteroidales was dominant in the phylum Bacteroidetes and they are obligate anaerobes but not thermotolerant43, and thus they increased in the 20 °C and 35 °C treatments, whereas they were not found in the 55 °C treatment. Two classes of Firmicutes exhibited differences because Clostridia are anaerobic whereas Bacilli are obligate or facultative aerobes44. Chloroflexi and Thermotogae are anaerobic thermophilic bacteria37, and they were only found in the 35 °C and 55 °C treatments. It has been reported that Chloroflexi can use halogenated organics as energy sources45 and that Thermotogae can utilize different complex carbohydrates to produce an class="Chemical">hydrogen gas during anaerobic digestion46.

This study confirmed that changes in the bacterial communities are key drivers of chn class="Chemical">anges in the variation of ARGs during anann> class="Species">aerobic digestion. Similar results were obtained based on tests of aerobic composting and water chlorination processes4147. Miller et al.48 also found that the ARG level of a digester was attributed to the survival or death of antibiotic resistant bacteria during anaerobic digestion. The variations in six phyla (Actinobacteria, Bacteroidetes, Chloroflexi, Firmicutes, Proteobacteria, and Thermotogae) were the main causes of the changes during anaerobic digestion and the differences at various temperatures. By contrast, the chemical characteristics of the digestion system were not correlated with the changes in ARGs. The results of the present study indicate that Proteobacteria and Bacteroidetes are important hosts of ARGs during anaerobic digestion (Table 1), which is consistent with previous reports that Bacteroidetes and Proteobacteria can resist a wide variety of antibiotics, including tetracycline, β-lactams, aminoglycosides, and erythromycin4950. The high temperature treatment effectively eliminated Proteobacteria and Bacteroidetes, thereby reducing tetM, tetQ, gyrA, and sul1 with high efficiency. Similarly, the increase in Bacteroidetes was responsible for the increases of tetC, tetM, tetQ, tetX, and sul1 under moderate and mesophilic temperatures during anaerobic digestion. In addition, the inhibition of HGT at high temperatures may also contribute to the lower RAs of ARGs in the thermophilic system4851.

Ann class="Chemical">an class="Species">aerobic digestion caclass="Chemical">nn> remove some pathogens from raw waste and thermophilic digestion is particularly effective5253. Several potentially opportunistic n class="Species">human pathogens such as Bacilli and Actinomycetales were removed by anaerobic digestion in this study. However, the abundance of Bacteroidetes, which includes some species that can cause infections of the peritoneal cavity, gastrointestinal surgery, and an class="Disease">appendicitis54, increased after moderate and mesophilic digestion. Actinomycetales and Corynebacterium are potential pathogens that can harbor sul and integrons755, where they comprised 29.1% of the total bacteria in the raw material, but they were removed by thermophilic digestion, whereas 0.26% and 0.17% remained under moderate and mesophilic digestion. Proteobacteria are also potentially pathogenic bacteria56, which were removed by thermophilic digestion, but 5.5% and 1.7% remained after moderate and mesophilic digestion. These results suggest that the digestion products from moderate and mesophilic treatments are potential risks for agricultural application1552.

Interestingly, in agreement with the results of most previous studies121314171934, integrase genes were significantly reduced at different temperatures, thereby demonstrating the universal reduction effect of n class="Chemical">anaerobic digestion on integrase genes. Proteobacteria, n class="Disease">order Actinomycetales in Actinobacteria, and class Bacilli in Firmicutes were the main carriers of integrons in the present study. Actinomycetales were the main producer of antibiotics, and thus they were common hosts of ARGs and integrons5758. Proteobacteria and Bacilli have also been reported as common carriers of integrons5960. The universal reduction of integrons by anaerobic digestion can be explained by the significant reduction of aerobic bacteria such as Actinomycetales and Bacilli after anaerobic digestion of the raw material.

In conclusion, the present study demonstrated that bacterial community succession was the main mechanism that affected the variation in n class="Chemical">an class="Chemical">ARGs ann>d integrase genes. Thermophilic ann class="Species">aerobic digestion reduced the abundance of mesophilic bacteria carrying ARGs and its an class="Chemical">ARG reduction was better than that of the moderate and mesophilic treatments. Anaerobic digestion universally reduced the abundance of integrase genes because it removed significant amounts of the aerobic hosts of integrons. Thermophilic digestion also has advantages in terms of methane production as well as speed and stability compared with moderate and mesophilic digestion33. The consumption of more energy to maintain a high temperature is the main disadvantage of thermophilic anaerobic digestion. However, given the reduction of ARGs, integrase genes, and pathogenic bacteria, as well as methane production, thermophilic anaerobic digestion is highly recommended for the treatment and reuse of animal manure.

Additional Information

How to cite this article: Sun, W. et al. Mechanism n class="Chemical">and Effect of Temperature on Variations in Antibiotic Resistance Genes during Anan class="Species">aerobic Digestion of Dairy Manure. Sci. Rep. 6, 30237; doi: 10.1038/n class="Chemical">srep30237 (2016).