Understanding the effects of pretreatment on anaerobic digestion of sludge waste from wastewater treatment plants is becoming increasingly important, as impetus moves towards the utilization of sludge for renewable energy production. Although the field of sludge pretreatment has progressed significantly over the past decade, critical questions concerning the underlying microbial interactions remain unanswered. In this study, a metagenomic approach was adopted to investigate the microbial composition and gene content contributing to enhanced biogas production from sludge subjected to a novel pretreatment method (maintaining pH at 10 for 8 days) compared to other documented methods (ultrasonic, thermal and thermal-alkaline).
Our results showed that pretreated sludge attained a maximum methane yield approximately 4-fold higher than that of the blank un-pretreated sludge set-up at day 17. Both the microbial and metabolic consortium shifted extensively towards enhanced biodegradation subsequent to pretreatment, providing insight for the enhanced methane yield. The prevalence of Methanosaeta thermophila and Methanothermobacter thermautotrophicus, together with the functional affiliation of enzymes-encoding genes suggested an acetoclastic and hydrogenotrophic methanogenesis pathway. Additionally, an alternative enzymology in Methanosaeta was observed.
This study is the first to provide a microbiological understanding of improved biogas production subsequent to a novel waste sludge pretreatment method. The knowledge garnered will assist the design of more efficient pretreatment methods for biogas production in the future.
Activated sludge technology is currently the most broadly-implemented biological method for biomass conversion in wastewater treatment plants (WWTPs) . However, vast quantities of highly organic waste activated sludge (WAS) is produced during this process, and this by-product mass continues to increase with the expansion of population and industry [1–4]. Sludge disposal by landfill or incineration may no longer be feasible in the near future due to land scarcity, high waste charge and increasingly stringent environmental control regulations [4, 5]. As a result, the strategy for sludge management is shifting towards its re-utilization as a potential source for renewable energy [6–8]. In this regard, the anaerobic digestion process represents an attractive means of sludge reduction while producing renewable energy in the form of biogas [8, 9]. Identifying efficient ways to improve methane production, a major biogas product from anaerobic digestion, has now become a topic of interest for numerous researchers [3, 10, 11].
The performance of an anaerobic digestion system has been shown to be tied closely to its microbial community structure . Methane production from WAS is a complex, multi-step process which involves multiple syntrophic interactions within the microbial consortium . Complex compounds (polysaccharides, proteins, nucleic acids, and lipids) are first converted to oligomers and monomers through the action of extracellular hydrolytic enzymes produced by the primary fermenting bacteria . Subsequently, the intermediate products are further transformed into acetate, carbon dioxide, hydrogen and formate by secondary fermenters . The final methanogenesis step is then conducted by methanogenic archaea, whose energy substrates are highly restricted to acetate, H2, CO2, formate or certain C1 compounds . To improve methane yield from sludge, enormous research efforts have been devoted to the development of pretreatment methods to accelerate sludge hydrolysis, including thermal , thermal-alkaline , ultrasonic , mechanical and thermo-chemical methods . Although the field of pretreatment research has progressed significantly in the past decade, many significant questions related to their effects on the underlying microbial interactions remain unanswered.
As the field of pretreatment method research is nearing a threshold, the accomplishments of the past are pushing on the door of microbiology to provide new insights [20, 21]. In this study, we aimed to investigate the microbial composition and gene content of sludge subjected to our novel pretreatment method (maintaining pH at 10 for 8 days) which leads to significantly enhanced methane generation compared to other documented methods (ultrasonic, thermal and thermal-alkaline) . It was reported in our previous study that both sludge hydrolysis and short-chain fatty acids (SCFAs; eg. acetic, butyric and propionic acid) accumulation were significantly enhanced when WAS was anaerobically fermented under the condition of pH 10 for 8 days . This phenomenon was suspected to be due to biotic factors rather than abiotic ones (e.g. alkaline hydrolysis) since much higher SCFAs accumulation and enzyme activities were observed in un-autoclaved sludge compared to autoclaved sludge . Further investigation by our group suggested that the solubilization of the sludge matrix, usually a hydrolysis event by the embedded extracellular enzymes, may contribute to the significant SCFAs improvement .
In this study, a shotgun metagenomic approach was chosen to study potential shifts in microbial communities and/or gene contents that could help explain elevated productions of methane under our novel pretreatment method [24, 25]. The latest advances in pyrosequencing technology afford new opportunities to undertake such metagenomic studies to explore the dynamics of microbial communities in time, space or under fluctuating environmental conditions with un-precedented levels of microbial diversity coverage and depth [26–28]. In addition to elucidating the microbiology underpinning the sludge pretreatment process , our study sought to improve knowledge of the diversity and physiology of participating syntrophs and methanogens, as well as the mechanism behind the establishment and maintenance of mutualistic cooperation. This knowledge will help to establish a better control over the hydrolysis and methanogenic processes, and promote pretreatment as part of a pertinent strategy for sludge management in WWTP. To the best of our knowledge, this study is the first to adopt a whole genome shotgun approach to link the knowledge between microbiology and engineering of sludge pretreatment methods. Additional statistical soundness was conferred to the surmised conclusion, as the comparative analysis was conducted between meta-datasets generated by the same sequencing method, from well-characterized experimental designs which only differ in known parameters by manipulation.
Results and discussion
Enhanced methane production by using a newly devised pretreatment strategy
Collected WAS was subjected to a novel pretreatment method (maintaining pH 10 for 8 days) , and its effects on methane production enhancement was characterized in this study. Figure 1 illustrates the experimental flowchart of obtaining end samples from the biogas-producing anaerobic digesters. Pretreated sludge for 30 days (P30), un-pretreated sludge for 30 days (UP30), or 40 days (UP40) were used as substrate. Pretreated sludge was significantly more effective in producing biogas, as total gas volume produced on day 8 (604.32 ml/g VSS-added) is 7.71 and 1.67 times of that in the un-pretreated sludge bioreactor on day 8 (78.43 ml/g VSS-added) and day 16 (360.99 ml/g VSS-added), respectively (Figure 2A-B). Methane production was 389.8 mL-CH4/g VSS-added using pretreated sludge as a substrate, which was 3.78 times that of un-pretreated sludge (103.2 mL-CH4/g VSS-added) (see Figure 2C). The generated methane on average constituted 70.5% and 59.1% of the total generated gas composition in the two respective bioreactors, while the rest was primarily CO2 (pretreated sludge, 3.6%; un-pretreated sludge, 9.2%) and other small amounts of N2, H2, NH3 and H2S (data not shown) (Figure 2D-E). It is worthwhile to note that the methane content was increased by 19.3%, and carbon dioxide content was decreased by 60.9% after the pretreatment. The reproducibility of this new pretreatment strategy validated the experimental success seen before and its technical robustness. This warranted investigation of biological mechanisms behind this significant bioreactor improvement in performance.
Biogas producing microbial community residing in anaerobic digester inferred by metagenome sequencing
To analyze biogas-producing microbes residing in our studied bioreactors in terms of community structure, gene content, metabolic capabilities and the role of specific organisms in biogas formation, a metagenomic approach using the 454 pyrosequencing was used. Statistical data summarizing the output sequencing quantity of the three independent runs of P30, UP30 and UP40 is given in Table 1. To detect differentially abundant features between the microbial communities, we used the Meta Genome Rapid Annotation using Subsystem Technology (MG-RAST) analysis pipeline, which has shown to be applicable to various metagenomic data performing taxonomic classification as well as functional annotation .
Statistics of 454 GS Junior pyrosequence datasets presented in this study
pH10 for 8d
Mean sequence length
Legends: By taxonomic classification based on M5NR database, Bacteria (82.13% on average) and Archaea (11.86% on average) are the dominant domains, with variable proportions in the bioreactor samples. UP 30 and UP40 are collected from same reactor at different time since semi-continuous culture.
Composition of the biogas-producing microbial community was obtained by taxonomic classification of reads base on the M5NR database on the MG-RAST platform. In all three reactor samples, Bacteria were the dominant superkingdom (Table 1). Domain-based allocation was primarily assigned to Bacteria (82.13% on average) and Archaea (11.86% on average), whilst the assigned reads to Eukaryota and Viruses altogether accounted for less than half a percent (Table 1). The bacterial proportion of the microbial community in P30 was the highest (85.46%) followed by UP30 (81.82%) and UP40 (79.10%). As for the Archaea domain, it ranged from 8.32% in P30, to 12.13% in UP30, and to 15.13% in UP40. The clustering pattern shown in Additional file 1 provides a graphical representation of the overall taxonomic similarity—the microbiomes of the unpretreated sludge bioreactors (UP30 and UP40) displayed high resemblance, in which the mature unpretreated sludge bioreactor (UP40, operated for longer time) is comparatively closer to the significantly diverse pretreated sludge bioreactor. The result provided the first insight into both the temporal dynamic and the enormous impact of pretreatment on microbiology. Downstream analyses focused on the Bacteria and Archaea domains (unless otherwise specified, the percentages below are representative of the identified reads within each domain per independent run).
Within the Bacteria domain, the top abundant phyla were Proteobacteria (UP30, 32.17%; UP40, 30.45%; P30, 42.28%), Bacteroidetes (UP30, 22.91%; UP40, 22.25%; P30, 23.38%), Firmicutes (UP30, 16.54%; UP40, 17.65%; P30, 10.97%), Actinobacteria (UP30, 7.13%; UP40, 6.94%; P30, 6.83%) and Chloroflexi (UP30, 5.00%; UP40, 5.21%; P30, 3.37%), they collectively account for over 0.7 of the bacterial reads (normalized between 0 and 1) for each of the three bioreactors (Figure 3A). Proteobacteria, Bacteroidetes and Firmicutes have been reported to be dominant phyla as well in similar analysis of anaerobic digestion of sludge . Further resolution at the class level revealed that the microbial compositions overlapped between the bioreactor samples. To be more precise, the dominant bacterial lineages presented in this study were related to anaerobic digestion, including Bacteroidia (UP30, 13.50%; UP40, 12.69%; P30, 8.90%), δ-proteobacteria (UP30, 13.08%; UP40, 12.71%; P30, 8.93%), Clostridia (UP30, 11.42%; UP40, 12.21%; P30, 7.31%), γ-proteobacteria (UP30, 8.12%; UP40, 7.90%; P30, 11.18%), Actinobacteria (class) (UP30, 7.37%; UP40, 7.20%; P30, 7.06%), α-proteobacteria (UP30, 5.78%; UP40, 5.28%; P30, 9.33%), β-proteobacteria (UP30, 5.73%; UP40, 4.86%; P30, 13.58%), Bacilli (UP30, 5.16%; UP40, 5.25%; P30, 3.53%), Flavobacteria (UP30, 4.70%; UP40, 4.76%; P30, 6.98%) and Chloroflexi (class) (UP30, 2.65%; UP40, 2.56%; P30, 1.84%) (Figure 3B).
With time, there was moderate temporal variation in the relative abundances of these bacterial members in the anaerobic digester fed with un-pretreated sludge (UP30 and UP40). While the relative proportion of these communities remained comparable between UP30 and UP40, a larger extent of differences could be observed when P30 was taken into comparison for certain lineages. The over-represented bacterial members in P30 included α-proteobacteria (9.33%), β-proteobacteria (13.58%) and γ-proteobacteria (11.18%); which are respectively 69.06%, 158.21% and 39.60% more than UP30 and UP40 on average. Overall, these sequences were originated from the WAS substrate, AGS inoculum and the biogas-producing symbiont within the system, hence the differences observed in their relative abundance were attributable to the sole variable of either using pretreated or un-pretreated sludge in the bioreactors. These results suggested that the bacterial communities that underlie the anaerobic digesters were dynamic, and they responded rapidly to the pretreated sludge substrate and change substantially over time.
Based on the phylogenetic affiliation of the metagenomic sequences, it is possible to form hypotheses regarding the metabolic functions of the groups . Proteobacteria are microorganisms involved in the initial steps of degradation, studies have shown that they are the main consumers of propionate, butyrate and acetate . Bacteroidetes are known proteolytic bacteria, responsible for the degradation of protein and subsequent fermentation of amino acids into acetate for acetoclastic methanogens [31, 32]. Concerning Firmicutes, they are syntrophic bacteria which degrade volatile fatty acids such as butyrate . The H2 generated from this process could be then uptaken by the hydrogenotrophic methanogens . The bacterial class Clostridia is frequently found in co-culture with other species in biomasss conversion systems and it includes a number of anaerobic species that are commonly associated with the decomposition of lignocelluloses and municipal solid waste [34, 35]. As for Chloroflexi, this group is often found in various wasterwater treatments such as anaerobic digesters and biological nutrient removal processes; their potential role in carbohydrate degradation has been reported in several studies [31, 36–38]. Overall, the resident bacteria manifested in the bioreactors represented a biomass decomposing community.
Archaeal representatives were less diverse and consisted of only three major families of methanogens, belonging to Methanosarcinaceae, Methanobacteriaceae and Methanosaetaceae (Table 2). The predominance of selected methanogen lineages in the Archaeal domain (represented as percentage of total identified reads in Archaea) has been observed in both production and laboratory-scale biogas reactors [23, 39, 40], and concluded in a recent meta-analysis of the collective microbial diversity [41, 42]. At the genus level, Methanosarcina was determined to be the most abundant methanogen in all three bioreactors (13.51% on average—UP30, 15.10%; UP40, 13.32%; P30, 12.10%), while Methanosaeta (10.67% on average—UP30, 13.13%; UP40, 10.31%; P30, 8.56%) and Methanothermobacter (10.35% on average— UP30, 10.73%; UP40, 9.41%; P30, 10.91%) contributed to slightly lesser portions. Anaerobic digesters are typical habitats to these three genera of methanogens , and the enrichment of Methanosarcina species was in congruence with other evaluation studies on primary sludge and WAS anaerobic digesters [44, 45]. The observation of a lesser portion of Methanothermobacter was expected, as the operating condition of the bioreactor (approximate 35°C) was not in favour of this thermophilic Archeon’s proliferation (55-65°C) .
Relative abundance of Archaea in pretreated and un-pretreated sludge bioreactors at day 30 and 40
Archaeal taxonomic affiliation
Percentage of all Archaeal reads
(phylum/ family/ genus/ species)
Legends: Archaeal domain consists three major families belonging to Methanosarcinaceae, Methanobacteriaceae and Methanosaetacea, while Methanosarcina, Methanosaeta and Methanothermobacter represent top abundant methanogens on the genus level. The decrease in abundance of Methanosarcina-related species in P30 was correlated with an enhanced bioproduction of acetate in the pretreated sludge. Except for phylum level, only candidates of abundance >8% are presented. Species are represented as percentage of all methanogens.
Methanosaeta species were observed to be more abundant in samples of bioreactor digesting un-pretreated sludge (UP30; 13.13%) than that with pretreated sludge (P30; 8.56%), but the population decreases with time (UP40; 8.56%). A similar observation was found in an earlier surveillance of the methanogenic population dynamics in anaerobic digesters, where the hitherto abundant Methanosaeta population decreased rapidly as the acetate concentration increased . As acetoclastic methanogens, both Methanosarcina and Methanosaeta are able to split acetate, oxidize the carboxyl-group to CO2 and reduce the methyl group to CH4. At the same time, while Methanosarcina thrive in environments with high acetate concentration; lower acetate concentrations benefit the dominance of Methanosaeta owing to their high affinity for acetate [12, 46–48]. As aforementioned, an enhanced bioproduction of acetic acid was reported in the pretreated sludge , the revealed shift in the methanogenic community could therefore be interpreted as an ecological consequence [30, 49, 50] of the sludge pretreatment.
Methane production via acetoclastic and hydrogenotrophic pathways
To identify the methane-producing organisms in the bioreactor samples, taxonomic and functional affiliation of the metagenomic reads were evaluated in parallel (represented as percentage of total identified methanogens). On the species level, it is common to all three bioreactor samples that Methanosaeta thermophila (12.36% on average—UP30, 14.48%; UP40, 12.36%; P30, 10.25%) and Methanothermobacter thermautotrophicus (9.34% on average—UP30, 9.50%; UP40, 8.59%; P30, 9.94%) are the over-represented methanogens (Table 2). Here, Methanosarcina species were not represented as one of the top dominant methanogens (13.51% on average—UP30, 15.10%; UP40, 13.32%; P30, 12.10%), this observation was linked to the intra-divergence of each family . M. thermophila is an obligate methanogen which consumes acetate only , whereas M. thermautotrophicus conserves energy by using H2 to reduce CO2 to CH4[33, 52]. While the pH was adjusted to be within the optimum range (7.0) for the two dominate methanogens after pretreatment (Figure 1), it was interesting to detect the prevalence of the thermophilic M. thermautotrophicus in the bioreactors that operated at mesophilic temperature . At the same time, functional enzyme-encoding genes for the two methanogenesis pathways were identified with reference to KEGG and Metacyc pathway database entries (Figure 4, Additional file 2) [54, 55]. Nonetheless, based on these results, it was evinced that the production of methane in the studied bioreactors was performed by methanogenic Archaea via acetoclastic and hydrogenotrophic pathways.
Acetoclastic and hydrogenotrohpic pathways are indeed common methanogenesis pathways reported in various reports [42, 56]; the affirmation of these pathways permitted subsequent in-depth investigation of the encompassed enzymology. The first step in acetoclastic methanogenesis is the formation of acetyl-CoA from acetate . Analysis of the recently completed genome of Methanosaeta thermophilia confirmed that the majority of the acetoclastic pathway are similar for Methanosaeta and Methanosarcina, except the enzymes employed for the first step of catalyzation . It was proposed that acetoclastic methanogenesis in Methanosaeta proceeded with a modified version of the pathway compared with Methanosarcina, which utilizes the acetate kinase/phosphotransacetylase pathway to convert acetate to acetyl-CoA [57, 58]. Taking M. thermophila as an example, its genome does not include a readily identifiable acetate kinase, and it has been postulated that this archeon utilizes an acetate transporter coupled with acetyl-CoA synthetases to convert acetate to acetyl-CoA, and the hydrolysis of pyrophosphate by inorganic pyrophosphatase (PPase) drives this reaction forward . In this study, the analysis of metagenomic datasets using SEED annotation indicated the presence of acetyl-CoA synthetases (EC 22.214.171.124) and inorganic pyrophosphatase (EC 126.96.36.199) in the compiled Methanosaeta bin. While the sequencing depth might account for the absence of genes for the acetate kinase/phosphotransacetylase pathway, the discovery of these two enzymes advocated an alternative pathway in acetotrohpic methanogenesis for Methanosaeta. Results from a recent study on a terephthalate-degrading bioreactor supported this hypothesis .
Functional affiliation related to higher methane production after sludge pretreatment
Degradation of the highly organic polymers represents the first and overall rate-limiting step for the mineralization of organic matter in activated sludge and anaerobic digested sludge treatment systems [59–61]. In our work, the studied pretreatment involved alteration of pH, which is an important parameter affecting both bacterial activity and metabolite pathways. KO-based annotations were used to understand how these phylogenetic trends could be used to predict the metabolic potential of these microbes. Figure 5 shows the subsystems that are related to higher methane production, including metabolism of amino acids, energy, carbohydrates, nucleotides, lipids, cofactors and vitamins, xenobiotics, as well as the fermentation of different substrates. These results revealed a general elevated expression of these faculties in the bioreactors fed with un-pretreated sludge over time (compare UP30 and UP40), whereas these levels were highest in the bioreactor digesting pretreated sludge (P30). Distribution of the functional systems was most divergent in P30, which showed predominance in metabolism consistent with a community shifted towards an enhanced biomass degradation metabolism. Herewith, the downstream analysis focused on the degradation of carbohydrates for a number of reasons. Firstly, the matrix of extra-cellular polymeric substances (EPS) that crosslink cells together remains the primary solid part of biological sludge, of which the polysaccharide is the predominant component [62–64]. Secondly, complex carbohydrate is the commonplace recalcitrant in the hydrolysis process, and bioprospecting for carbohydrate active enzymes is a topic of interest in metagenomics . As a result, the focus of our study was hence placed on the degradation of this dominant recalcitrant component. (Please see Additional file 3 for further insights towards protein and lipid degradation.)
As sequencing technology provides access to a remarkable array of microbial functional capacity, sequence-based data mining is an important prospect in metagenomic projects [66, 67]. To understand how the microbial community mediates the solibulization of the sludge matrix, a carbohydrate-active enzyme (CAZy) characterization of the metagenomic datasets was performed after ab initio gene prediction [68, 69]. Conventional sequence homology-based enzyme discovery introduces a bias towards the identification of candidates similar to known enzymes, rather than enzymes with low sequence identity and potentially divergent biochemical properties . The entries in the CAZy database contains both experimentally verified and putative carbohydrate-active enzyme domains, hence this search strategy would be able to provide a better insight into the catalysis of biochemical reactions [68, 70].
In this study, a total of 1917 and 107 gene modules were recognized across 52 glycosyl hydrolases (GH) and 9 carbohydrate binding modules (CBM) respectively (Tables 3 and 4). GH is defined as a widespread group of enzymes which hydrolyzes the glycosidic bond between carbohydrates or between that and a non-carbohydrate moiety . On the other hand, CBM are contiguous amino acids within a carbohydrate-active enzyme with a discreet fold which bears the carbohydrate-binding activity . In contrast to the small number of enzymes devoted to the hydrolysis of the main chain of cellulose, hemicelluloses and pectins (GH 5, 6, 7, 8, 9, 12, 44, 45, 48 and 74), all three metagenomes displayed a variety of enzymes that digest the side chains of these polymers and oligosaccharides. Families of GH 2 and 3, which contain a large range of glycosidases cleaving non-reducing carbohydrates in oligosaccharides and the side chains of hemicelluloses and pectins, were particularly abundant in the bioreactor samples (29.10% of all GH—UP30, 27.56%; UP40, 29.55%; P30, 24.20%). CBM 48 was the most dominant in the three metagenomes, it binds glycogen and is commonly found in a range of enzymes that act on branched substrates, such as the hydrolysis of glycogen, amylopectin and dextrin by isoamylase, pulluanase and branching enzyme. Overall, the number of GH hit returns appeared to increase with time in the bioreactor fed with un-pretreated sludge (468 counts in UP30; 792 counts in UP40); while the bioreactor digesting the pretreated sludge harboured 1.40 times the GH count in the control counterpart at the same time. As for CBM, P30 (51 counts) had a higher abundance than UP40 (38 counts), and made up to approximately 2.83 times that of UP30. Generally, there appeared to be a dynamic process in enzymology within the anaerobic digesters with time and substrate nature The heightened GH and CBM abundance in the bioreactor with pretreated sludge may possibly correlate with the higher sludge solubilization proposed in previous studies .
Identification of glycoside hydrolase (GH) family classified by CAZy database
Number of counts
Number of counts
CAZy GH family
CAZy GH family
GH 42 M
GH 65 m
GH 65 N
GH 67 M
GH 67 N
GH 32 N
Total number of counts
Legend: Dynamic process in enzymology with time and substrate nature was observed in the bioreactors, the heightened GH and abundance (in bold) in the pretreated-sludge bioreactor were correlated with a higher sludge solubilization, and in turns a higher SCFA production.
Identification of carbohydrate binding module (CBM) family classified by CAZy database
Number of counts
CAZy CBM family
Total number of counts
Legend: Dynamic process in enzymology with time and substrate nature was observed in the bioreactors, the heightened CBM abundance (in bold) in the pretreated-sludge bioreactor were correlated with a higher sludge solubilization, and in turns a higher SCFA production.
The treatment and disposal of excess sludge accumulation represents a bottleneck in wastewater treatment plants worldwide due to environmental, economic, social and legal factors. Reduction of excess sludge is becoming one of the biggest challenges in biological wastewater treatment [71, 72]. The anaerobic digestion process is a promising technique, along with which waste is eliminated and methane is produced as a valuable renewable energy source. Understanding the impact of pretreatment on the microbiology at different stages of anaerobic digestion is important; it ultimately impacts the performance of a bioreactor. This study represents the first attempt to gain an in-depth metagenomic perspective, with regard to taxonomic and functional aspects, on increased biogas (methane) production when a newly devised pretreatment method is used. Dual taxonomic and functional analysis indicated the microbial and metabolic consortium shifted extensively towards enhanced biodegradation, also seen in the methanogenesis pathways. Altogether, these results presented here help further microbiological understanding in sludge pretreatment research for anaerobic digestion. As microbiology is the ultimate driver for the anaerobic digestion process, this new perspective encourages a closer engagement between the engineering and microbiology knowledge pools. Further studies involving deeper sequencing of the metagenomes to characterize the decomposing and methanogenic microbiome at higher frequency and replicate number parallel to the physical-chemical characterization are required for each WAS pretreatment intervention.
Collection of waste activated sludge and anaerobic granular sludge
The waste activated sludge (WAS) used as substrate for methane production was collected from the secondary sedimentation tank of a municipal wastewater treatment plant in Shanghai, China. The anaerobic granular sludge (AGS) used as the inoculum for methane production was obtained from the upward-flow anaerobic sludge blanket (USAB) reactor of a food wastewater treatment plant in Yixing, China. The main features of WAS and AGS can be found in Additional file 4. The AGS was cultured in a laboratory up-flow anaerobic sludge bed (UASB) by synthetic wastewater (Additional file 5) prior to its use as the inoculum for methane production, where the use of synthetic wastewater for AGS culturing is a common practice in anaerobic digestion studies investigating the microbiology and biodegradability [4, 73, 74]. The hydrolytic retention time of UASB was 6 h, and the AGS concentration in UASB was approximately 29165 mg/L.
Pretreatment of waste activated sludge
Biogas-producing, anaerobic digester samples were obtained from a reproduction of the experiment described in our previous work , which introduced a new sludge pretreatment strategy (pH 10 for 8 d) for enhanced methane production from WAS. During the pretreatment, the pH of the collected WAS was adjusted to 10.0 by addition of 5 M NaOH or 4 M HCl with an automatic titrator, in a tightly-sealed reactor which was mechanically stirred (80 rpm) at 35 ± 1°C for 8 d. A blank test was also conducted, in which the sludge was not pretreated but mechanically stirred for 8 d. Then the pH in all reactors was adjusted to pH 7.0, and 30 mL of AGS was added to each reactor for methane production. All reactors were sealed with rubber stoppers and mechanically stirred at 80 rpm without further pH control. The retention time of all reactors was 17 d and the organic loading rate was 1.88 kg/COD/m3/d. The gas generation was detected by water displacement and the methane concentration was measured by a gas chromatograph (GC-14B, Shimadzu, Japan). The analyses of VSS were conducted in accordance with APHA Standard Methods . In this study, the methane yield was reported as the amount of methane generated per gram of volatile suspended solid added (mL-CH4/gVSS-added) unless otherwise stated. The full unit of methane production is ml-CH4/gVSS-added, in which the VSS-added refers to the initial VSS to facilitate comparison of methane production. Moreover, there was no SCFA in the original sludge before anaerobic fermentation, meaning that the initial VSS did not contain volatile fatty acids—hence the risk of overestimation of specific methane yield is improbable .
Preparation of anaerobic digester samples
The granular samples cultured with pretreated (pH 10 for 8 d) and un-pretreated WAS were collected when the methane production rate reached the maximum on day 17, and then further cultured by a semi-continuous method to maintain the methane production rate. In replicates of three, the semi-continuous flow reactors with working volume of 1.0 L received 470 mL of pH 10 pretreated sludge and un-pretreated sludge, respectively. Everyday, 52 mL sludge mixture was manually wasted from each reactor, and 52 mL pretreated and un-pretreated sludge was added to the respective reactors; AGS was fed back to reactors by sifting the wasted sludge through a sifter with aperture of 0.2 mm. Further procedure details were described in .
As a stable methane generation continued to be maintained by the semi-continuous culturing method, 50 mL anaerobic digester samples were collected from the pretreatment set-ups (as pooled sample) at day 30 (P30), control set-up at day 30 (UP30) and control set-up at day 40 (UP40) for microbial studies. The two samples collected at day 30 (P30, UP30) served to characterize taxonomic and functional difference between pretreated and unpretreated samples. Both sludge bioreactor were seeded with the same synthetic wastewater cultured AGS substrate, and later gradually replaced by the experimental sludge (pretreated or unpretreated WAS) completely as the batch experiment continued. Since both bioreactors received the same inoculums, therefore the only existing variable was the pretreatment, ascribing to the observed difference between the microbiomes. As for UP40, this additional un-pretreated sample at day 40 was prepared to provide further insight of the microbial dynamic in response to time. Figure 1 shows the experimental scheme for the preparation of the anaerobic digester samples.
DNA extraction and 454 sequencing
Samples (P30, UP30 and UP40) were kept at −20°C and thawed once only for nucleic acid extraction using PowerSoil® DNA Isolation Kit (Mobio, Carlsbad, USA), conducted within 24 hours of sample collection. DNA quantification was performed using Nanodrop 2000 spectrophotometer (Nanodrop Technologies, Wilmington, USA). Sequencing of total DNA derived from each was done by the shotgun sequencing approach on the 454 GS Junior system, and the sequencing quality was maintained by 454 propriety software. The raw sequence reads are submitted under the Sequence Read Archive (SRA) accession number SAMN01924662(P30), SAMN01924663 (UP30), SAMN01924664 (UP40).
To obtain a quantitative picture of the taxonomic composition, sequence datasets were characterized without a prior assembly step . The unassembled sequence reads were submitted to the online metagenome analysis tool Metagenome Rapid Annotation using Subsystem Technology (MG-RAST analysis pipeline 3.0) . Following the quality trimming with default setting, taxonomic and functional profiling was performed using the M5 non-redundant protein database (M5NR) and KEGG Orthology (KO) reference database  Parameters used were e-value cutoff of 1 × 10-5, minimum alignment length of 50 base pairs, and a minimum percentage identity at 50%. Relative abundance of the Archaeal and Bacterial groups were determined by the percentages of respective reads over total assigned reads. Parallel processing was performed similarly on Megan4 to confirm the clustering pattern of taxonomic structure . Only the taxonomic and functional annotation provided by MG-RAST was relied upon for all subsequent analysis.
To identify important carbohydrate- metabolism related gene candidates, pyrosequence reads were first subjected to gene prediction by a heuristic approach using MetaGeneMark 1.0 . Next, the amino acid sequences were screened against catalytic domains (CD) characteristic for particular glycoside hydrolase (GH) families or for a particular carbohydrate binding module (CBM) as classified by the CAZy database . HMM models were downloaded from Pfam database version 25.0 and used as a database for pfam_scan  (e-value cutoff of 1x10-4). Methanogenesis related gene candidates were identified with reference to KEGG and Metacyc pathway database [54, 55].
FCL is a Full Professor at School of Biological Sciences, former Dean of Faculty of Science, Board member of Initiative on Clean Energy and Environment, The University of Hong Kong, and Fellow of American Association for the Advancement of Science (AAAS). RKHH is a post-doc, JL, HMT, MSB and TJP are PhD students, MTW is an MPhil student in the same institution. YC is a Full Professor and DZ is a post-doc at State Key Laboratory of Pollution Control and Resources Reuse, School of Environmental Science and Engineering, Tongji University.
Anaerobic granular sludge
Carbohydrate binding modules
Sequence Read Archive
Short-chain fatty acids
Upward-flow anaerobic sludge blanket
Volatile suspended solid
Waste activated sludge
Wastewater treatment plants
This study was partially supported by Initiative on Clean Energy and Environment, The University of Hong Kong (HKU-ICEE) and a Bioinformatics Center fund from Nanjing Agricultural University. We are grateful to P. Benzie, C. Bush and I. Tasovski for editing the manuscript.
School of Biological Sciences, The University of Hong Kong
State Key Laboratory of Pollution Control and Resources Reuse, School of Environmental Science and Engineering, Tongji University
Bioinformatics Center, Nanjing Agricultural University
Kim S, Kim J, Park C, Kim TH, Lee M, Kim SW, Lee J: Effects of various pretreatments for enhanced anaerobic digestion with waste activated sludge.J Biosci Bioeng 2003, 95:271–275.
Reynolds TDR, P A: Unit operations and process in environmental engineering. 2nd edition. Boston, MA, USA: PWS Publishing; 1996.
Jahng D, Cui R: Enhanced methane production from anaerobic digestion of disintegrated and deproteinized excess sludge.Biotechnol Lett 2006, 28:531–538.View Article
Zhang D, Chen Y, Zhao Y, Zhu X: New sludge pretreatment method to improve methane production in waste activated sludge digestion.Environ Sci Technol 2010, 44:4802–4808.View Article
Weemaes M, Grootaerd H, Simoens F, Verstraete W: Anaerobic digestion of ozonized biosolids.Water Res 2000, 34:2330–2336.View Article
Shin SG, Lee S, Lee C, Hwang K, Hwang S: Qualitative and quantitative assessment of microbial community in batch anaerobic digestion of secondary sludge.Bioresource Technol 2010, 101:9461–9470.View Article
Mossakowska A, Hellstrom BG, Hultman B: Strategies for sludge handling in the Stockholm region.Water Sci Technol 1998, 38:111–118.
Odegaard H, Paulsrud B, Karlsson I: Wastewater sludge as a resource: sludge disposal strategies and corresponding treatment technologies aimed at sustainable handling of wastewater sludge.Water Sci Technol 2002, 46:295–303.
Riviere D, Desvignes V, Pelletier E, Chaussonnerie S, Guermazi S, Weissenbach J, Li T, Camacho P, Sghir A: Towards the definition of a core of microorganisms involved in anaerobic digestion of sludge.ISME J 2009, 3:700–714.View Article
Supaphol S, Jenkins SN, Intomo P, Waite IS, O’Donnell AG: Microbial community dynamics in mesophilic anaerobic co-digestion of mixed waste.Bioresource Technol 2011, 102:4021–4027.View Article
Ting CH, Lin KR, Lee DJ, Tay JH: Production of hydrogen and methane from wastewater sludge using anaerobic fermentation.Water Sci Technol 2004, 50:223–228.
Raskin L, Griffin ME, McMahon KD, Mackie RI: Methanogenic population dynamics during start-up of anaerobic digesters treating municipal solid waste and biosolids.Biotechnol Bioeng 1998, 57:342–355.View Article
Stams AJ: Metabolic interactions between anaerobic bacteria in methanogenic environments.Antonie Van Leeuwenhoek 1994, 66:271–294.View Article
Thauer RK: Biochemistry of methanogenesis: a tribute to Marjory Stephenson. 1998 Marjory Stephenson Prize Lecture.Microbiology 1998,144(Pt 9):2377–2406.View Article
Li YY, Noike T: Upgrading of anaerobic-digestion of waste activated-sludge by thermal pretreatment.Water Sci Technol 1992, 26:857–866.
Vlyssides AG, Karlis PK: Thermal-alkaline solubilization of waste activated sludge as a pre-treatment stage for anaerobic digestion.Bioresource Technol 2004, 91:201–206.View Article
Neis U, Tiehm A, Nickel K, Zellhorn M: Ultrasonic waste activated sludge disintegration for improving anaerobic stabilization.Water Res 2001, 35:2003–2009.View Article
Tanaka S, Kobayashi T, Kamiyama K, Bildan MLNS: Effects of thermochemical pretreatment on the anaerobic digestion of waste activated sludge.Water Sci Technol 1997, 35:209–215.
Wei H, Tucker MP, Baker JO, Harris M, Luo Y, Xu Q, Himmel ME, Ding SY: Tracking dynamics of plant biomass composting by changes in substrate structure, microbial community, and enzyme activity.Biotechnol Biofuels 2012, 5:20.View Article
Lau MW, Gunawan C, Dale BE: The impacts of pretreatment on the fermentability of pretreated lignocellulosic biomass: a comparative evaluation between ammonia fiber expansion and dilute acid pretreatment.Biotechnol Biofuels 2009, 2:30.View Article
Yuan HY, Chen YG, Zhang HX, Jiang S, Zhou Q, Gu GW: Improved bioproduction of short-chain fatty acids (SCFAs) from excess sludge under alkaline conditions.Environ Sci Technol 2006, 40:2025–2029.View Article
He PJ, Yu GH, Shao LM, He PP: Toward understanding the mechanism of improving the production of volatile fatty acids from activated sludge at pH 10.0.Water Res 2008, 42:4637–4644.View Article
Venter JC, Remington K, Heidelberg JF, Halpern AL, Rusch D, Eisen JA, Wu DY, Paulsen I, Nelson KE, Nelson W: Environmental genome shotgun sequencing of the Sargasso Sea.Science 2004, 304:66–74.View Article
Sogin ML, Morrison HG, Huber JA, Mark Welch D, Huse SM, Neal PR, Arrieta JM, Herndl GJ: Microbial diversity in the deep sea and the underexplored “rare biosphere”.P Natl Acad Sci USA 2006, 103:12115–12120.View Article
Krober M, Bekel T, Diaz NN, Goesmann A, Jaenicke S, Krause L, Miller D, Runte KJ, Viehover P, Puhler A, Schluter A: Phylogenetic characterization of a biogas plant microbial community integrating clone library 16S-rDNA sequences and metagenome sequence data obtained by 454-pyrosequencing.J Biotechnol 2009, 142:38–49.View Article
Hollister EB, Engledow AS, Hammett AJM, Provin TL, Wilkinson HH, Gentry TJ: Shifts in microbial community structure along an ecological gradient of hypersaline soils and sediments.ISME J 2010, 4:829–838.View Article
Wirth R, Kovács E, Maróti G, Bagi Z, Rákhely G, Kovács KL: Characterization of a biogas-producing microbial community by short-read next generation DNA sequencing.Biotechnol Biofuels 2012, 5:4.View Article
Zhang H, Banaszak JE, Parameswaran P, Alder J, Krajmalnik-Brown R, Rittmann BE: Focused-Pulsed sludge pre-treatment increases the bacterial diversity and relative abundance of acetoclastic methanogens in a full-scale anaerobic digester.Water Res 2009, 43:4517–4526.View Article
Dinsdale EA, Edwards RA, Hall D, Angly F, Breitbart M, Brulc JM, Furlan M, Desnues C, Haynes M, Li LL: Functional metagenomic profiling of nine biomes.Nature 2008, 452:629-U628.View Article
Okabe S, Ariesyady HD, Ito T: Functional bacterial and archaeal community structures of major trophic groups in a full-scale anaerobic sludge digester.Water Res 2007, 41:1554–1568.View Article
Deppenmeier U: The unique biochemistry of methanogenesis.Prog Nucleic Acid Re 2002, 71:223–283.View Article
Blackall LL, Burrell PC, O’Sullivan C, Song H, Clarke WP: Identification, detection, and spatial resolution of Clostridium populations responsible for cellulose degradation in a methanogenic landfill leachate bioreactor.Appl Environ Microb 2004, 70:2414–2419.View Article
Klocke M, Mahnert P, Mundt K, Souidi K, Linke B: Microbial community analysis of a biogas-producing completely stirred tank reactor fed continuously with fodder beet silage as mono-substrate.Syst Appl Microbiol 2007, 30:139–151.View Article
Sekiguchi Y, Takahashi H, Kamagata Y, Ohashi A, Harada H: In situ detection, isolation, and physiological properties of a thin filamentous microorganism abundant in methanogenic granular sludges: a novel isolate affiliated with a clone cluster, the green non-sulfur bacteria, subdivision I.Appl Environ Microb 2001, 67:5740–5749.View Article
Blackall LL, Bjornsson L, Hugenholtz P, Tyson GW: Filamentous Chloroflexi (green non-sulfur bacteria) are abundant in wastewater treatment processes with biological nutrient removal.Microbiol-Sgm 2002, 148:2309–2318.
Sekiguchi Y, Yamada T, Hanada S, Ohashi A, Harada H, Kamagata Y: Anaerolinea thermophila gen. nov., sp nov and Caldilinea aerophila gen. nov., sp nov., novel filamentous thermophiles that represent a previously uncultured lineage of the domain Bacteria at the subphylum level.Int J Syst Evol Micr 2003, 53:1843–1851.View Article
Zakrzewski M, Goesmann A, Jaenicke S, Junemann S, Eikmeyer F, Szczepanowski R, Al-Soud WA, Sorensen S, Puhler A, Schluter A: Profiling of the metabolically active community from a production-scale biogas plant by means of high-throughput metatranscriptome sequencing.J Biotechnol 2012, 158:248–258.View Article
Jaenicke S, Ander C, Bekel T, Bisdorf R, Droge M, Gartemann KH, Junemann S, Kaiser O, Krause L, Tille F: Comparative and joint analysis of two metagenomic datasets from a biogas fermenter obtained by 454-pyrosequencing.PLoS One 2011, 6:e14519.View Article
Liu WT, Lykidis A, Chen CL, Tringe SG, McHardy AC, Copeland A, Kyrpides NC, Hugenholtz P, Macarie H, Olmos A, Monroy O: Multiple syntrophic interactions in a terephthalate-degrading methanogenic consortium.ISME J 2011, 5:122–130.View Article
Nelson MC, Morrison M, Yu Z: A meta-analysis of the microbial diversity observed in anaerobic digesters.Bioresour Technol 2011, 102:3730–3739.View Article
Whitman WB, Liu YC: Metabolic, phylogenetic, and ecological diversity of the methanogenic archaea.Ann Ny Acad Sci 2008, 1125:171–189.View Article
Raskin L, McMahon KD, Zheng DD, Stams AJM, Mackie RI: Microbial population dynamics during start-up and overload conditions of anaerobic digesters treating municipal solid waste and sewage sludge.Biotechnol Bioeng 2004, 87:823–834.View Article
Raskin L, Stroot PG, McMahon KD, Mackie RI: Anaerobic codigestion of municipal solid waste and biosolids under various mixing conditions - I. Digester performance.Water Res 2001, 35:1804–1816.View Article
Liesack W, Grosskopf R, Janssen PH: Diversity and structure of the methanogenic community in anoxic rice paddy soil microcosms as examined by cultivation and direct 16S rRNA gene sequence retrieval.Appl Environ Microb 1998, 64:960–969.
Shigematsu T, Tang Y, Kobayashi T, Kawaguchi H, Morimura S, Kida K: Effect of dilution rate on metabolic pathway shift between aceticlastic and nonaceticlastic methanogenesis in chemostat cultivation.Appl Environ Microbiol 2004, 70:4048–4052.View Article
Nettmann E, Bergmann I, Pramschufer S, Mundt K, Plogsties V, Herrmann C, Klocke M: Polyphasic analyses of methanogenic archaeal communities in agricultural biogas plants (vol 76, pg 2540, 2010).Appl Environ Microb 2011, 77:394–394.View Article
Lee C, Kim J, Hwang K, O’Flaherty V, Hwang S: Quantitative analysis of methanogenic community dynamics in three anaerobic batch digesters treating different wastewaters.Water Res 2009, 43:157–165.View Article
Smith KS, Ingram-Smith C: Methanosaeta, the forgotten methanogen?Trends Microbiol 2007, 15:150–155.View Article
Smith DR, DoucetteStamm LA, Deloughery C, Lee HM, Dubois J, Aldredge T, Bashirzadeh R, Blakely D, Cook R, Gilbert K: Complete genome sequence of Methanobacterium thermoautotrophicum Delta H: Functional analysis and comparative genomics.J Bacteriol 1997, 179:7135–7155.
Boone DW WB, Rouviere P: Diversity and taxonomy of methanogens. In Methanogenesis ecology, physiology, biochemistry & genetics. Edited by: JG F. New York: Chapman & Hall; 1993:35–80.
Kanehisa M, Goto S: KEGG: Kyoto Encyclopedia of Genes and Genomes.Nucleic Acids Res 2000, 28:27–30.View Article
Karp PD, Caspi R, Foerster H, Fulcher CA, Kaipa P, Krummenacker M, Latendresse M, Paley S, Rhee SY, Shearer AG: The MetaCyc Database of metabolic pathways and enzymes and the BioCyc collection of Pathway/Genome Databases.Nucleic Acids Res 2008, 36:D623-D631.
Krause L, Diaz NN, Edwards RA, Gartemann KH, Kromeke H, Neuweger H, Puhler A, Runte KJ, Schluter A, Stoye J: Taxonomic composition and gene content of a methane-producing microbial community isolated from a biogas reactor.J Biotechnol 2008, 136:91–101.View Article
Metcalf WW, Rother M: Anaerobic growth of Methanosarcina acetivorans C2A on carbon monoxide: An unusual way of life for a methanogenic archaeon.P Natl Acad Sci USA 2004, 101:16929–16934.View Article
Singhwissmann K, Ferry JG: Transcriptional regulation of the phosphotransacetylase-encoding and acetate kinase-encoding genes (Pta and Ack) from methanosarcina-thermophila.J Bacteriol 1995, 177:1699–1702.
Frolund B, Griebe T, Nielsen PH: Enzymatic-activity in the activated-sludge floc matrix.Appl Microbiol Biot 1995, 43:755–761.View Article
Gessesse A, Dueholm T, Petersen SB, Nielsen PH: Lipase and protease extraction from activated sludge.Water Res 2003, 37:3652–3657.View Article
Pletschke BI, Burgess JE: Hydrolytic enzymes in sewage sludge treatment: A mini-review.Water Sa 2008, 34:343–349.
Carrere H, Dumas C, Battimelli A, Batstone DJ, Delgenes JP, Steyer JP, Ferrer I: Pretreatment methods to improve sludge anaerobic degradability: a review.J Hazard Mater 2010, 183:1–15.View Article
Cescutti P, Toffanin R, Pollesello P, Sutherland IW: Structural determination of the acidic exopolysaccharide produced by a Pseudomonas sp. strain 1.15.Carbohydr Res 1999, 315:159–168.View Article
Li LL, McCorkle SR, Monchy S, Taghavi S, van der Lelie D: Bioprospecting metagenomes: glycosyl hydrolases for converting biomass.Biotechnol Biofuels 2009, 2:10.View Article
Rubin EM, Hess M, Sczyrba A, Egan R, Kim TW, Chokhawala H, Schroth G, Luo SJ, Clark DS, Chen F: Metagenomic discovery of biomass-degrading genes and genomes from Cow rumen.Science 2011, 331:463–467.View Article
Suen G, Scott JJ, Aylward FO, Adams SM, Tringe SG, Pinto-Tomas AA, Foster CE, Pauly M, Weimer PJ, Barry KW: An insect herbivore microbiome with high plant biomass-degrading capacity.PLoS Genet 2010,6(9):e1001129.View Article
Cantarel BL, Coutinho PM, Rancurel C, Bernard T, Lombard V, Henrissat B: The Carbohydrate-Active EnZymes database (CAZy): an expert resource for Glycogenomics.Nucleic Acids Res 2009, 37:D233–238.View Article
Zhu W, Lomsadze A, Borodovsky M: Ab initio gene identification in metagenomic sequences.Nucleic Acids Res 2010, 38:e132.View Article
Li LL, Taghavi S, McCorkle SM, Zhang YB, Blewitt MG, Brunecky R, Adney WS, Himmel ME, Brumm P, Drinkwater C: Bioprospecting metagenomics of decaying wood: mining for new glycoside hydrolases.Biotechnol Biofuels 2011, 4:23.View Article
Chu LB, Wang JL, Wang B, Xing XH, Yan ST, Sun XL, Jurcik B: Changes in biomass activity and characteristics of activated sludge exposed to low ozone dose.Chemosphere 2009, 77:269–272.View Article
Wei YS, Van Houten RT, Borger AR, Eikelboom DH, Fan YB: Minimization of excess sludge production for biological wastewater treatment.Water Res 2003, 37:4453–4467.View Article
Bassin JP, Pronk M, Muyzer G, Kleerebezem R, Dezotti M, van Loosdrecht MCM: Effect of elevated salt concentrations on the aerobic granular sludge process: linking microbial activity with microbial community structure.Appl Environ Microb 2011, 77:7942–7953.View Article
del Rio AV, Morales N, Isanta E, Mosquera-Corral A, Campos JL, Steyer JP, Carrere H: Thermal pre-treatment of aerobic granular sludge: Impact on anaerobic biodegradability.Water Res 2011, 45:6011–6020.View Article
Eaton AD, Franson MAH, Association APH: Association AWW. Standard Methods for the Examination of Water & Wastewater. Federation WE: American Public Health Association; 1998.
Meyer F, Paarmann D, D’Souza M, Olson R, Glass EM, Kubal M, Paczian T, Rodriguez A, Stevens R, Wilke A: The metagenomics RAST server - a public resource for the automatic phylogenetic and functional analysis of metagenomes.BMC Bioinforma 2008, 9:386.View Article
Huson DH, Mitra S, Ruscheweyh HJ, Weber N, Schuster SC: Integrative analysis of environmental sequences using MEGAN4.Genome Res 2011, 21:1552–1560.View Article
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.