Adaptation of Methanogenic Inocula to Anaerobic Digestion of Maize Silage

A comparative study on polyaluminum chloride (PACl) and Moringa oleifera (MO) chemically enhanced primary treatment (CEPT) in enhanced biogas production: anaerobic digestion performance and the Gompertz model

A comparative study was performed to estimate biogas production from sludge produced by organic and inorganic chemically enhanced primary treatments (CEPTs). To this end, the effects of two coagulants, polyaluminum chloride (PACl) and Moringa oleifera (MO), on CEPT and biogas production in anaerobic digestion were surveyed within an incubation period of 24 days. The optimal dosage and pH of PACl and MO were optimized in terms of sCOD, TSS and VS parameters in the CEPT process. Next, the digestion performance of anaerobic digestion reactors fed with sludge obtained from PACl and MO coagulants at a batch mesophilic reactor (37 ± 1 °C) was surveyed from the biogas production, volatile solid reduction (VSR) and Gompertz model. At the optimal conditions (pH = 7 and dosage = 5 mg L^-1), the removal efficiency of COD, TSS and VS in CEPT assisted with PACL was 63, 81 and 56%, respectively. Moreover, CEPT assisted with MO led to the removal efficiency of COD, TSS and VS until 55, 68 and 25%, respectively. The highest methane yield (0.598 L g_{VS removed}^-1) was obtained in an anaerobic digestion reactor with sludge from the MO coagulant. The anaerobic digestion of CEPT sludge instead of primary sludge resulted in higher sCOD removal efficiency, and 43-50% of sCOD was observed compared with the removal of 32% for the primary sludge. Furthermore, the high coefficient of determination (R²) demonstrated the trustworthy predictive precision of the modified Gompertz model with actual data. The combination of CEPT and anaerobic digestion, especially using natural coagulants, provides a cost-effective and practical way to increase BMP from primary sludge.

Culture adaptation for enhanced biogas production from birch wood applying stable carbon isotope analysis to monitor changes in the microbial community

Birch wood is a potential feedstock for biogas production in Northern Europe; however, the lignocellulosic matrix is recalcitrant preventing efficient conversion to methane. To improve digestibility, birch wood was thermally pre-treated using steam explosion at 220 °C for 10 min. The steam-exploded birch wood (SEBW) was co-digested with cow manure for a period of 120 days in continuously fed CSTRs where the microbial community adapted to the SEBW feedstock. Changes in the microbial community were tracked by stable carbon isotopes- and 16S r RNA analyses. The results showed that the adapted microbial culture could increase methane production up to 365 mL/g VS day, which is higher than previously reported methane production from pre-treated SEBW. This study also revealed that the microbial adaptation significantly increased the tolerance of the microbial community against the inhibitors furfural and HMF which were formed during pre-treatment of birch. The results of the microbial analysis indicated that the relative amount of cellulosic hydrolytic microorganisms (e.g. Actinobacteriota and Fibrobacterota) increased and replaced syntrophic acetate bacteria (e.g. Cloacimonadota, Dethiobacteraceae, and Syntrophomonadaceae) as a function of time. Moreover, the stable carbon isotope analysis indicated that the acetoclastic pathway became the main route for methane production after long-term adaptation. The shift in methane production pathway and change in microbial community shows that for anaerobic digestion of SEBW, the hydrolysis step is important. Although acetoclastic methanogens became dominant after 120 days, a potential route for methane production could also be a direct electron transfer among Sedimentibacter and methanogen archaea.

Hydrodynamic disintegration of thickened excess sludge and maize silage to intensify methane production: Energy effect and impact on microbial communities

The aim of this project was to study the combination of two methods to increase methane production: feedstock pretreatment by hydrodynamic disintegration and co-digestion of maize silage (MS) with thickened excess sludge (TES). Disintegration of TES alone resulted in a 15% increase in specific methane production from 0.192 Nml/gVS (TES + MS) to 0.220 Nml/gVS (pretreated TES + MS). The energy balance revealed additional energy (0.14 Wh) would cover only the energy expenditure for the mechanical pretreatment and would not allow for net energy profit. Identification of the methanogenic consortia by 16S rRNA gene amplicon sequencing revealed that Chloroflexi, Bacteroidota, Firmicutes, Proteobacteria and Actinobacteriota were five most abundant bacteria phyla, with Methanothrix and Methanolinea as the dominant methanogens. Principal component analysis did not show any effect of feedstock pretreatment on methanogenic consortia. Instead, the composition of inoculum was the decisive factor in shaping the microbial community structure.

PLSDA-batch: a multivariate framework to correct for batch effects in microbiome data

Microbial communities are highly dynamic and sensitive to changes in the environment. Thus, microbiome data are highly susceptible to batch effects, defined as sources of unwanted variation that are not related to and obscure any factors of interest. Existing batch effect correction methods have been primarily developed for gene expression data. As such, they do not consider the inherent characteristics of microbiome data, including zero inflation, overdispersion and correlation between variables. We introduce new multivariate and non-parametric batch effect correction methods based on Partial Least Squares Discriminant Analysis (PLSDA). PLSDA-batch first estimates treatment and batch variation with latent components, then subtracts batch-associated components from the data. The resulting batch-effect-corrected data can then be input in any downstream statistical analysis. Two variants are proposed to handle unbalanced batch x treatment designs and to avoid overfitting when estimating the components via variable selection. We compare our approaches with popular methods managing batch effects, namely, removeBatchEffect, ComBat and Surrogate Variable Analysis, in simulated and three case studies using various visual and numerical assessments. We show that our three methods lead to competitive performance in removing batch variation while preserving treatment variation, especially for unbalanced batch $\times $ treatment designs. Our downstream analyses show selections of biologically relevant taxa. This work demonstrates that batch effect correction methods can improve microbiome research outputs. Reproducible code and vignettes are available on GitHub.

Magnetite-based nanoparticles and nanocomposites for recovery of overloaded anaerobic digesters

The effect of magnetite nanoparticles and nanocomposites (magnetite nanoparticles impregnated into graphene oxide) supplement on the recovery of overloaded laboratory batch anaerobic reactors was assessed using two types of starting inoculum: anaerobic granular sludge (GS) and flocculent sludge (FS). Both nanomaterials recovered methane production at a dose of 0.27 g/L within 40 days in GS. Four doses of magnetite nanoparticles from 0.075 to 1 g/L recovered the process in FS systems between 30 and 50 days relaying on the dose. The presence of nanomaterials helped to reverse the effect of volatile fatty acids inhibition and enabled microbial communities to recover but also favoured the development of certain microorganisms over others. In GS reactors, the methanogenic population changed from being mostly acetoclastic (Methanothrix soehngenii) to being dominated by hydrogenotrophic species (Methanobacterium beijingense). Nanomaterial amendment may serve as a preventative measure or provide an effective remedial solution for system recovery following overloading.

Comprehensive ADM1 Extensions to Tackle Some Operational and Metabolic Aspects in Anaerobic Digestion

Modelling in anaerobic digestion will play a crucial role as a tool for smart monitoring and supervision of the process performance and stability. By far, the Anaerobic Digestion Model No. 1 (ADM1) has been the most recognized and exploited model to represent this process. This study aims to propose simple extensions for the ADM1 model to tackle some overlooked operational and metabolic aspects. Extensions for the discontinuous feeding process, the reduction of the active working volume, the transport of the soluble compound from the bulk to the cell interior, and biomass acclimation are presented in this study. The model extensions are included by a change in the mass balance of the process in batch and continuous operation, the incorporation of a transfer equation governed by the gradient between the extra- and intra- cellular concentration, and a saturation-type function where the time has an explicit influence on the kinetic parameters, respectively. By adding minimal complexity to the existing ADM1, the incorporation of these phenomena may help to understand some underlying process issues that remain unexplained by the current model structure, broadening the scope of the model for control and monitoring industrial applications.

Microbial Community Successional Changes in a Full-Scale Mesophilic Anaerobic Digester from the Start-Up to the Steady-State Conditions

Anaerobic digestion is a widely used technology for sewage sludge stabilization and biogas production. Although the structure and composition of the microbial communities responsible for the process in full-scale anaerobic digesters have been investigated, little is known about the microbial successional dynamics during the start-up phase and the response to variations occurring in such systems under real operating conditions. In this study, bacterial and archaeal population dynamics of a full-scale mesophilic digester treating activated sludge were investigated for the first time from the start-up, performed without adding external inoculum, to steady-state operation. High-throughput 16S rRNA gene sequencing was used to describe the microbiome evolution. The large majority of the reads were affiliated to fermentative bacteria. Bacteroidetes increased over time, reaching 22% of the total sequences. Furthermore, Methanosaeta represented the most abundant methanogenic component. The specific quantitative data generated by real-time PCR indicated an enrichment of bacteria and methanogens once the steady state was reached. The analysis allowed evaluation of the microbial components more susceptible to the shift from aerobic to anaerobic conditions and estimation of the microbial components growing or declining in the system. Additionally, activated sludge was investigated to evaluate the microbial core selected by the WWTP operative conditions.

Archaeal community dynamics in biogas fermentation at various temperatures assessed by mcrA amplicon sequencing using different primer pairs

In this study, the taxonomic and functional diversity of methanogenic archaea in two parallel 120 l fermenters operated at different temperatures and fed with maize silage was estimated by mcrA metabarcoding analysis using two typical primer pairs (ML and MLA) amplifying part of the functional methyl coenzyme M reductase (mcrA) gene. The alpha diversity indices showed that the ML primer pair detected a higher Operational Taxonomic Unit (OTU) abundance compared to the MLA primer pair and methanogen diversity was significantly lower in the 60 °C fermenters. The beta diversity analysis showed the methanogenic community clustered together at 50 °C and 40° and was statistically different from the 60 °C community. Similar, to alpha diversity, beta diversity was also significantly different between primer pairs. At all temperatures analysed, the primer pairs showed a different abundance of the different methanogenic OTUs, e.g. more OTUs relative to Methanoculleus sp. with the ML primer pair, and more OTUs corresponding to Methanobacterium sp. with the MLA primer pair. Moreover, OTUs corresponding to Methanosphaera sp. and Methanobrevibacter sp. were found only by using ML primer pair, while the MLA primer pair detected sequences corresponding to Methanothrix sp.

Soil-Derived Inocula Enhance Methane Production and Counteract Common Process Failures During Anaerobic Digestion

Although soil-borne methanogens are known to be highly diverse and adapted to extreme environments, their application as potential (anaerobic) inocula to improve anaerobic digestion has not been investigated until now. The present study aimed at evaluating if soil-derived communities can be beneficial for biogas (methane, CH₄) production and endure unfavorable conditions commonly associated with digestion failure. Nine study sites were chosen and tested for suitability as inoculation sources to improve biogas production via in situ measurements (CH₄ fluxes, physical and chemical soil properties, and abundance of methanogens) and during a series of anaerobic digestions with (a) combinations of both sterile or unsterile soil and diluted fermenter sludge, and (b) pH-, acetate-, propionate-, and ammonium-induced disturbance. Amplicon sequencing was performed to assess key microbial communities pivotal for successful biogas production. Four out of nine tested soil inocula exerted sufficient methanogenic activity and repeatedly allowed satisfactory CH₄/biogas production even under deteriorated conditions. Remarkably, the significantly highest CH₄ production was observed using unsterile soil combined with sterile sludge, which coincided with both a higher relative abundance of methanogens and predicted genes involved in CH₄ metabolism in these variants. Different bacterial and archaeal community patterns depending on the soil/sludge combinations and disturbance variations were established and these patterns significantly impacted CH₄ production. Methanosarcina spp. seemed to play a key role in CH₄ formation and prevailed even under stressed conditions. Overall, the results provided evidence that soil-borne methanogens can be effective in enhancing digestion performance and stability and, thus, harbor vast potential for further exploitation.

Anaerobic digestion of bean straw applying a fungal pre-treatment and using cow manure as co-substrate

The significant amounts of agriculture residues such as bean straw (BS) in rural areas, advises its valorisation for energy recovery. The feasibility of using BS for biogas production through anaerobic digestion was assessed. Prior to this, a fungal pre-treatment to hydrolyse BS with Pleutorus ostreatus was studied at 30°C and 100 rpm in orbital incubators with 1, 10 and 30 mg fungus/g straw for 14, 21 and 28 days. Then, anaerobic digestion experiments were performed in batch with cow manure (CM) as co-substrate and pre-treated BS at ratios (g/g total solids) of 1/2, 1/3, 1/5 and 0/1. Maximum lignin (18%) and hemicellulose (44%) degradation occurred at 30 mg fungus/g straw and 28 days, along with the highest total methane yield (38 mL CH₄/g VS loaded). The total amount of methane decreased when increasing CM in the experiments (701.4-191.5 mL CH₄), suggesting inhibition owed to a component of CM. Self-sustained biogas production of BS occurred due to the presence of bacteria (i.e. Bacilli and Bacteroidia) and archea (i.e. Methanobacteria and Methanomicrobia). However, the usage of a full-active inoculum should be studied for higher biogas production rates.

Deciphering the effects of temperature on bio-methane generation through anaerobic digestion

Anaerobic digestion (AD) is a sustainable wastewater treatment technology which facilitates energy, nutrient, and water recovery from organic wastes. The agricultural and industrial wastes are suitable substrates for the AD, as they contain a high level of biodegradable compounds. The aim of this study was to examine the AD of three different concentrations of phenol (100, 200, and 300 mg/L) containing wastewater with and without co-substrate (acetate) at four different temperatures (25, 35, 45, and 55 °C) to produce methane (CH₄)-enriched biogas. It was observed that the chemical oxygen demand (COD) and phenol removal efficiencies of up to 76% and 72%, respectively, were achieved. The CH₄ generation was found higher in anaerobic batch reactors (ABRs) using acetate as co-substrate, with the highest yield of 189.1 μL CH₄ from 500 μL sample injected, obtained using 200 mg/L of phenol at 35 °C. The results revealed that the performance of ABR in terms of degradation efficiency, COD removal, and biogas generation was highest at 35 °C followed by 55, 45, and 25 °C indicating 35 °C to be the optimum temperature for AD of phenolic wastewater with maximum energy recovery. Scanning electron microscopy (SEM) revealed that the morphology of the anaerobic sludge depends greatly on the temperature at which the system is maintained which in turn affects the performance and degradation of toxic contaminants like phenol. It was observed that the anaerobic sludge maintained at 35 °C showed uniform channels leading to higher permeability through enhanced mass transfer to achieve higher degradation rates. However, the denser sludge as in the case of 55 °C showed lesser permeability leading to limited transfer and thus reduced treatment. Quantitative real-time PCR (qPCR) analysis revealed a more noteworthy change in the population of the microbial communities due to temperature than the presence of phenol with the methanogens being the dominating species at 35 °C. The findings suggest that the planned operation of the ABR could be a promising choice for CH₄-enriched biogas and COD removal from phenolic wastewater.

Impact of Organic Loading Rate in Volatile Fatty Acids Production and Population Dynamics Using Microalgae Biomass as Substrate

Volatile fatty acids (VFAs) are regarded as building blocks with a wide range of applications, including biofuel production. The traditional anaerobic digestion used for biogas production can be alternatively employed for VFAs production. The present study aimed at maximizing VFAs productions from Chlorella vulgaris through anaerobic digestion by assessing the effect of stepwise organic loading rates (OLR) increases (3, 6, 9, 12 and 15 g COD L^-1 d^-1). The biological system was proven to be robust as organic matter conversion efficiency into VFAs increased from 0.30 ± 0.02 COD-VFAs/COD_in at 3 g COD L^-1 d^-1 to 0.37 ± 0.02 COD-VFAs/COD_in at 12 g COD L^-1d^-1. Even though, the hydrolytic step was similar for all studied scenario sCOD/tCOD = 0.52-0.58), the highest OLR (15 g COD L^-1 d^-1) did not show any further increase in VFAs conversion (0.29 ± 0.01 COD-VFAs/COD_in). This fact suggested acidogenesis inhibition at 15 g COD L^-1d^-1. Butyric (23-32%), acetic (19-26%) and propionic acids (11-17%) were the most abundant bioproducts. Population dynamics analysis revealed microbial specialization, with a high presence of Firmicutes followed by Bacteroidetes. In addition, this investigation showed the microbial adaptation of Euryarchaeota species at the highest OLR (15 g COD L^-1d^-1), evidencing one of the main challenges in VFAs production (out-competition of archaea community to avoid product consumption). Stepwise OLR increase can be regarded as a tool to promote VFAs productions. However, acidogenic inhibition was reported at the highest OLR instead of the traditional hydrolytic barriers. The operational conditions imposed together with the high VFAs and ammonium concentrations might have affected the system yields. The relative abundance of Firmicutes (74%) and Bacteroidetes (20%), as main phyla, together with the reduction of Euryarchaeota phylum (0.5%) were found the best combination to promote organic matter conversion into VFAs.

Full-scale anaerobic reactor samples would be more suitable than lab-scale anaerobic reactor and natural samples to inoculate the wheat straw batch anaerobic digesters

This study evaluated the effects of the inocula from natural wetland, lab-scale and full-scale anaerobic reactors on wheat straw anaerobic digestion. Three replicate batch reactors were constructed for each inoculum to investigate the reactor performances and microbial communities. Reactors seeded with full-scale reactor samples were started up most rapidly, achieved the highest methane production, and were recognized as the higher efficient reactors. The dominance of acetoclastic methanogens, including Methanosaeta and Methanoscrina, was crucial for the higher efficient reactors, whereas hydrogenotrophic methanogens were dominant in other reactors. Genus Treponema, which could enhance the cellulose degradation and conduct homoacetogenesis, was first reported to be dominant in the bacterial communities of high efficient reactors. Inoculum sources and process conditions were suggested to be the deterministic factors in shaping the microbial communities in the higher efficient reactors. These findings contribute to the startup of new anaerobic reactors.

Bacterial and Fungal Diversity Inside the Medieval Building Constructed with Sandstone Plates and Lime Mortar as an Example of the Microbial Colonization of a Nutrient-Limited Extreme Environment (Wawel Royal Castle, Krakow, Poland)

Biodeterioration is a serious threat to cultural heritage objects and buildings. The deterioration of a given material often incurs irreparable losses in terms of uniqueness and historical value. Hence preventive actions should be taken. One important challenge is to identify microbes involved in the biodeterioration process. In this study, we analyzed the microbial diversity of an ancient architectonical structure of the Rotunda of Sts. Felix and Adauctus, which is a part of the Wawel Royal Castle located in Krakow, Poland. The Rotunda is unavailable to tourists and could be treated as an extreme habitat due to the low content of nutrients coming either from sandstone plates bound with lime mortar or air movement. Microbial diversity was analyzed with the use of the high-throughput sequencing of marker genes corresponding to fragments of 16S rDNA (for Bacteria) and ITS2 (internal transcribed spacer 2) (for Fungi). The results showed that the microbial community adhered to wall surfaces is, to a large extent, endemic. Furthermore, alongside many microorganisms that could be destructive to masonry and mortar (e.g., Pseudomonas, Aspergillus), there were also bacteria, such as species of genera Bacillus, Paenisporosarcina, and Amycolatopsis, that can positively affect wall surface properties by reducing the damage caused by the presence of other microorganisms. We also showed that airborne microorganisms probably have little impact on the biodeterioration process as their abundance in the microbial community adhered to the ancient walls was very low.

Changes in the Substrate Source Reveal Novel Interactions in the Sediment-Derived Methanogenic Microbial Community

Methanogenesis occurs in many natural environments and is used in biotechnology for biogas production. The efficiency of methane production depends on the microbiome structure that determines interspecies electron transfer. In this research, the microbial community retrieved from mining subsidence reservoir sediment was used to establish enrichment cultures on media containing different carbon sources (tryptone, yeast extract, acetate, CO₂/H₂). The microbiome composition and methane production rate of the cultures were screened as a function of the substrate and transition stage. The relationships between the microorganisms involved in methane formation were the major focus of this study. Methanogenic consortia were identified by next generation sequencing (NGS) and functional genes connected with organic matter transformation were predicted using the PICRUSt approach and annotated in the KEGG. The methane production rate (exceeding 12.8 mg CH₄ L⁻¹ d⁻¹) was highest in the culture grown with tryptone, yeast extract, and CO₂/H_2. The analysis of communities that developed on various carbon sources casts new light on the ecophysiology of the recently described bacterial phylum Caldiserica and methanogenic Archaea representing the genera Methanomassiliicoccus and Methanothrix. Furthermore, it is hypothesized that representatives of Caldiserica may support hydrogenotrophic methanogenesis.

Acclimatization contributes to stable anaerobic digestion of organic fraction of municipal solid waste under extreme ammonia levels: Focusing on microbial community dynamics

The organic fraction of municipal solid waste (OFMSW) is an abundant and sustainable substrate for the anaerobic digestion (AD) process, yet ammonia released during OFMSW hydrolysis could result in suboptimal biogas production. Acclimatized ammonia tolerant microorganisms offer an efficient way to alleviate ammonia inhibition during AD. This study aimed to achieve an efficient AD of OFMSW under extreme ammonia levels and elucidate the dynamics of the acclimatized microbial community. Thus, two mesophilic continuous stirred tank reactors (CSTR), fed only with OFMSW, were successfully acclimatized up to 8.5 g NH₄⁺-N/L, and their methane yields fluctuated <10%, compared to the methane yields without ammonia addition. Microbiological analyses showed that Methanosaeta concilii and Methanosarcina soligelidi were the dominant methanogens at low and high ammonia levels, respectively. Whilst, a unique metabolic pathway shift, from aceticlastic to hydrogenotrophic methanogenesis, of M. soligelidi was identified during the acclimatization process.

Assessment of the start-up process of anaerobic digestion utilizing swine manure: 13C fractionation of biogas and microbial dynamics

The aim of this study was to investigate how the microbial community structure adapts during the start-up phase and how the ¹³C fractionation of biogas reflects the microbial population dynamics in two parallel swine manure-fed anaerobic digesters. Two swine manure-fed reactors for the start-up of continuously stirred tank reactors at mesophilic condition were evaluated. Changes in community structure were monitored using 16S rRNA high-throughput sequencing to measure the abundance of fermenting bacteria and methanogens. Digesters with relatively stable Methanosarcinaceae started up successfully and contained high gas production and low levels of propionate. In contrast, the digester that experienced a difficult start-up period had reduced Methanosarcinaceae along with accumulated propionate and low gas production. Specific gas production, specific methane production, and ¹³C fractionation of biogas were influenced significantly by Methanosarcinaceae, Methanobacteriaceae, and Clostridiaceae, indicating that the ¹³C fractionation of biogas had significant potential to reflect microbial population changes and digester performance during the start-up period.

Comparative analysis of deep sequenced methanogenic communities: identification of microorganisms responsible for methane production

BACKGROUND

Although interactions between microorganisms involved in biogas production are largely uncharted, it is commonly accepted that methanogenic Archaea are essential for the process. Methanogens thrive in various environments, but the most extensively studied communities come from biogas plants. In this study, we employed a metagenomic analysis of deeply sequenced methanogenic communities, which allowed for comparison of taxonomic and functional diversity as well as identification of microorganisms directly involved in various stages of methanogenesis pathways.

RESULTS

A comprehensive metagenomic approach was used to compare seven environmental communities, originating from an agricultural biogas plant, cattle-associated samples, a lowland bog, sewage sludge from a wastewater treatment plant and sediments from an ancient gold mine. In addition to the native consortia, two laboratory communities cultivated on maize silage as the sole substrate were also analyzed. Results showed that all anaerobic communities harbored genes of all known methanogenesis pathways, but their abundance varied greatly between environments and that genes were encoded by different methanogens. Identification of microorganisms directly involved in different stages of methane production revealed that hydrogenotrophic methanogens, such as Methanoculleus, Methanobacterium, Methanobrevibacter, Methanocorpusculum or Methanoregula, predominated in most native communities, whereas acetoclastic Methanosaeta seemed to be the key methanogen in the wastewater treatment plant. Furthermore, in many environments, the methylotrophic pathway carried out by representatives of Methanomassiliicoccales, such as Candidatus Methanomethylophilus and Candidatus Methanoplasma, seemed to play an important role in methane production. In contrast, in stable laboratory reactors substrate versatile Methanosarcina predominated.

CONCLUSIONS

The metagenomic approach presented in this study allowed for deep exploration and comparison of nine environments in which methane production occurs. Different abundance of methanogenesis-related functions was observed and the functions were analyzed in the phylogenetic context in order to identify microbes directly involved in methane production. In addition, a comparison of two metagenomic analytical tools, MG-RAST and MetAnnotate, revealed that combination of both allows for a precise characterization of methanogenic communities.

Responses and successions of sulfonamides, tetracyclines and fluoroquinolones resistance genes and bacterial community during the short-term storage of biogas residue and organic manure under the incubator and natural conditions

Biogas residue and organic manure are frequently used for crop planting. However, the evaluation of antibiotic resistant bacteria (ARB), antibiotic resistance genes (ARGs) and bacterial community before their applications to fields is still lacking. This study monitored the variations of bacteria resistant to sulfadiazine, tetracycline and norfloxacin, 57 resistance genes for sulfonamides, tetracyclines and fluoroquinolones as well as the bacterial community during the 28-day aerobic storage of biogas residue and organic manure by using viable plate counts, high-throughput qPCR and Illumina MiSeq sequencing methods. Then two storage conditions, incubator (25 °C) and natural environment, were used to assess the responses of ARB and ARGs to the environmental factors. Results showed that a total of 35 and 21 ARGs were detected in biogas residue and organic manure, respectively. ARB and ARGs were enriched up to 8.01-fold in biogas residue after the 28-day storage, but varied in a narrow range during the storage of organic manure. Compared with the incubator condition, the proliferation of ARB and ARGs in biogas residue under the natural condition was relatively inhibited by the varied and complicated environmental factors. However, we found that there was no significant difference of ARB and ARGs in organic manure between the incubator and natural conditions. Bacterial community was also shifted during the storage of biogas residue, especially Bacteroidetes_VC2.1_Bac22, Aequorivita, Luteimonas and Arenimonas. Network analysis revealed that the relationship in biogas residue was much more complicated than that in organic manure, which ultimately resulted in large successions of ARB and ARGs during the short-term storage of biogas residue. Therefore, we suggest that further measures should be taken before the application of biogas residue to fields.