Microbial core communities in activated sludge plants are strongly affected by immigration and geography.

BACKGROUND: The microbiota in wastewater treatment plants (WWTPs) and incoming wastewater is critical for the treatment process, the preservation of natural ecosystems and human health, and for the recovery of resources and achievement of sustainability goals. Both core species and conditionally rare and abundant taxa (CRAT) are considered process-critical but little is known about identity as well as true functional and ecological importance. Here, we present a comprehensive investigation of the microbiota of 84 municipal activated sludge (AS) plants with nutrient removal treating ~ 70% of all wastewater within a confined geographical area, Denmark (43,000 km2). With the use of an ecosystem-specific database (MiDAS 5.2), species-level classification allowed us to investigate the core and CRAT species, whether they were active, and important factors determining their presence. RESULTS: We established a comprehensive catalog of species with names or placeholder names showing each plant contained approx. 2,500 different species. Core and CRAT represented in total 258 species, constituting around 50% of all reads in every plant. However, not all core and CRAT could be regarded as process-critical as growth rate calculations revealed that 43% did not grow in the AS plants and were present only because of continuous immigration from the influent. Analyses of regional microbiota differences and distance decay patterns revealed a stronger effect for species than genera, demonstrating that geography had a clear effect on the AS microbiota, even across a limited geographical area such as Denmark (43,000 km2). CONCLUSIONS: The study is the first comprehensive investigation of WWTPs in a confined geographical area providing new insights in our understanding of activated sludge microbiology by introducing a concept of combining immigration and growth calculation with identifying core and CRAT to reveal the true ecosystem-critical organisms. Additionally, the clear biogeographical pattern on this scale highlights the need for more region-level studies to find regional process-critical taxa (core and CRAT), especially at species and amplicon sequence variant (ASV) level.

Effect of feeding frequency on the anaerobic digestion of berry fruit waste.

On-site anaerobic digesters for small agricultural farms typically have feeding schedules that fluctuate according to farm operations. Shocks in feeding, particularly for putrescible waste can disrupt the stable operation of a digester. The effect of intermittent feeding on the anaerobic digestion of rejected raspberries was investigated in four 3L reactors operated in semicontinuous mode for 350 days at 38 °C with a hydraulic retention time of 25 days and an organic loading rate (OLR) of 1gVS/L/d. During the acclimatisation period (147 days) the organic loading was 5 feeds per week. The feeding regime of two reactors was then changed while maintaining the same OLR and HRT to one weekly feed event in one reactor and 3 equal feeds per week in another. The feeding regime did not significantly affect specific methane yield (369 ± 47 L/kgVS on average) despite very different weekly patterns in methane production. Volatile fatty acids (VFA) comprised >83 % of the organics in the effluent, while the rest included non-inhibitory concentrations of phenolic compounds (515-556 mg gallic acid/L). The microbial composition and relative abundance of predominant groups in all reactors were the archaeal genera Methanobacterium and Methanolinea and the bacterial phyla Bacteridota and Firmicutes. Increasing the OLR to 2gVS/L/d on day 238 resulted in failure of all reactors, attributed to the insufficient alkalinity to counterbalance the VFA produced, and the pH decrease below 6. Overall results suggests that optimal digestion of raspberry waste is maintained despite variations in feeding frequency, but acidification can occur with OLR changes.

Inhibitory mechanisms on dry anaerobic digestion: Ammonia, hydrogen and propionic acid relationship.

Inhibitory pathways in dry anaerobic digestion are still understudied and current knowledge on wet processes cannot be easily transferred. This study forced instability in pilot-scale digesters by operating at short retention times (40 and 33 days) in order to understand inhibition pathways over long term operation (145 days). The first sign of inhibition at elevated total ammonia concentrations (8 g/l) was a headspace hydrogen level over the thermodynamic limit for propionic degradation, causing propionic accumulation. The combined inhibitory effect of propionic and ammonia accumulation resulted in further increased hydrogen partial pressures and n-butyric accumulation. The relative abundance of Methanosarcina increased while that of Methanoculleus decreased as digestion deteriorated. It was hypothesized that high ammonia, total solids and organic loading rate inhibited syntrophic acetate oxidisers, increasing their doubling time and resulting in its wash out, which in turn inhibited hydrogenotrophic methanogenesis and shifted the predominant methanogenic pathway towards acetoclastic methanogenesis at free ammonia over 1.5 g/l. C/N increases to 25 and 29 reduced inhibitors accumulation but did not avoid inhibition or the washout of syntrophic acetate oxidising bacteria.

Full-scale activated sludge transplantation reveals a highly resilient community structure.

Well-functioning and stable microbial communities are critical for the operation of activated sludge (AS) wastewater treatment plants (WWTPs). Bioaugmentation represents a potentially useful approach to recover deteriorated systems or to support specific AS processes, but its application in full-scale WWTPs is generally problematic. We conducted a massive transplantation (in one day) exchanging AS from a donor to a recipient full-scale WWTP with similar process type (biological removal of nitrogen and phosphorus) and performance, but with differences in microbial community structure. The treatment performance in the recipient plant was not compromised and the effluent quality remained stable. The AS community structure of the recipient plant was initially very similar to the donor AS, but it almost completely restored the pre-transplantation structure approximately 40 days after transplantation, corresponding to 3 times the solid retention time. Most of the unique species of donor AS added to recipient AS disappeared quickly, although some disappeared more slowly the following months, indicating some survival and potentially a time limited function in the recipient plant. Moreover, the addition in higher abundance of most species already present in the recipient AS (e.g., the polyphosphate accumulating organisms) or the reduction of the abundance of unwanted bacteria (e.g., filamentous bacteria) in the recipient plant was not successful. Moreover, we observed similar abundance patterns after transplantation for species belonging to different functional guilds, so we did not observe an increase of the functional redundancy. Investigations of the microbial community structure in influent wastewater revealed that for some species the abundance trends in the recipient plant were closely correlated to their abundance in the influent. We showed that a very resilient microbial community was responsible for the outcome of the transplantation of AS at full-scale WWTP, potentially as a consequence of mass-immigration from influent wastewater. The overall results imply that massive transplantation of AS across different WWTPs is not a promising strategy to permanently solve operational problems. However, by choosing a compatible AS donor, short term mitigation of serious operational problems may be possible.

Exploring the microbial influence on seasonal nitrous oxide concentration in a full-scale wastewater treatment plant using metagenome assembled genomes.

Nitrous oxide is a highly potent greenhouse gas and one of the main contributors to the greenhouse gas footprint of wastewater treatment plants (WWTP). Although nitrous oxide can be produced by abiotic reactions in these systems, biological N2O production resulting from the imbalance of nitrous oxide production and reduction by microbial populations is the dominant cause. The microbial populations responsible for the imbalance have not been clearly identified, yet they are likely responsible for strong seasonal nitrous oxide patterns. Here, we examined the seasonal nitrous oxide concentration pattern in Avedøre WWTP alongside abiotic parameters, the microbial community composition based on 16S rRNA gene sequencing and already available metagenome-assembled genomes (MAGs). We found that the WWTP parameters could not explain the observed pattern. While no distinct community changes between periods of high and low dissolved nitrous oxide concentrations were determined, we found 26 and 28 species with positive and negative correlations to the seasonal N2O concentrations, respectively. MAGs were identified for 124 species (approximately 31% mean relative abundance of the community), and analysis of their genomic nitrogen transformation potential could explain this correlation for four of the negatively correlated species. Other abundant species were also analysed for their nitrogen transformation potential. Interestingly, only one full-denitrifier (Candidatus Dechloromonas phosphorivorans) was identified. 59 species had a nosZ gene predicted, with the majority identified as a clade II nosZ gene, mainly from the phylum Bacteroidota. A correlation of MAG-derived functional guilds with the N2O concentration pattern showed that there was a small but significant negative correlation with nitrite oxidizing bacteria and species with a nosZ gene (N2O reducers (DEN)). More research is required, specifically long-term activity measurements in relation to the N2O concentration to increase the resolution of these findings.

Seasonal microbial community dynamics complicates the evaluation of filamentous bulking mitigation strategies in full-scale WRRFs.

The activated sludge wastewater treatment process has been thoroughly researched in more than 100 years, yet there are still operational challenges that have not been fully resolved. Such a challenge is the control of filamentous bulking caused by the overgrowth of certain filamentous bacteria. In this study, we tested different mitigation strategies to reduce filamentous bulking, caused by two common filamentous genera found in full-scale water resource recovery facilities (WRRF), Candidatus Microthrix and Candidatus Amarolinea. PAX dosing, ozone addition, hydrocyclone implementation, and the addition of nano-coagulants were tested as mitigation strategies in four parallel treatment lines in a full-scale WRRF over three consecutive years. Unexpectedly, the activated sludge settleability was not affected by any of the mitigation strategies. Some of the strategies appeared to have a strong mitigating effect on the two filamentous species. However, detailed analyses of the microbial communities revealed strong recurrent seasonal variations in all four lines, including the control line which masked the real effect. After removing the effect of the seasonal variation by using a time-series decomposition approach, it was clear that the filamentous bacteria were mostly unaffected by the mitigation strategies. Only PAX dosing had some effect on Ca. Microthrix, but only on one species, Ca. Microthrix subdominans, and not on the most common Ca. Microthrix parvicella. Overall, our study shows the importance of long-term monitoring of microbial communities at species level to understand the normal seasonal pattern to effectively plan and execute full-scale experiments. Moreover, the results highlight the importance of using parallel reference treatment lines when evaluating the effect of mitigation strategies in full-scale treatment plants.

Quantification of Biologically and Chemically Bound Phosphorus in Activated Sludge from Full-Scale Plants with Biological P-Removal.

Phosphorus (P) is present in activated sludge from wastewater treatment plants in the form of metal salt precipitates, extracellular polymeric substances, or bound into the biomass, for example, as intracellular polyphosphate (poly-P). Several methods for a reliable quantification of the different P-fractions have recently been developed, and this study combines them to obtain a comprehensive P mass-balance of activated sludge from four enhanced biological phosphate removal (EBPR) plants. Chemical characterization by ICP-OES and sequential P fractionation showed that chemically bound P constituted 38-69% of total P, most likely in the form of Fe, Mg, or Al minerals. Raman microspectroscopy, solution state 31P NMR, and 31P MAS NMR spectroscopy applied before and after anaerobic P-release experiments, were used to quantify poly-P, which constituted 22-54% of total P and was found in approximately 25% of all bacterial cells. Raman microspectroscopy in combination with fluorescence in situ hybridization was used to quantify poly-P in known polyphosphate-accumulating organisms (PAO) (Tetrasphaera, Candidatus Accumulibacter, and Dechloromonas) and other microorganisms known to possess high level of poly-P, such as the filamentous Ca. Microthrix. Interestingly, only 1-13% of total P was stored by unidentified PAO, highlighting that most PAOs in the full-scale EBPR plants investigated are known.

Microbial communities across activated sludge plants show recurring species-level seasonal patterns.

Microbial communities in activated sludge (AS) are the core of sanitation in wastewater treatment plants (WWTPs). Microbial communities in AS have shown seasonal changes, however, long-term experiments (>2 years) are rarely conducted, limiting our understanding of the true seasonal dynamics in WWTPs. In this study, we resolved the microbial seasonal dynamics at the species level in four municipal full-scale WWTPs, sampled every 7-10 days, during 3-5 consecutive years. By applying a new time-series analysis approach, we revealed that the seasonal pattern was species-specific, where species belonging to the same functional guild or genus may show different seasonal dynamics. Species could be grouped into cohorts according to their seasonal patterns, where seasonal cohorts showed repeatable annual dynamics across years and plants. Species were also grouped according to their net growth rate in the AS (i.e., growing species and disappearing species). Growing species were more prevailing in spring and autumn cohorts, while disappearing species, which were only present due to the continuous immigration from influent wastewater, were mostly associated with winter and spring cohorts. Most known process-critical species, such as nitrifiers, polyphosphate accumulating organisms and filamentous organisms, showed distinct species-specific patterns. Overall, our study showed that overarching seasonal patterns affected microbial species in full-scale AS plants, with similar seasonal patterns across plants for many dominant species. These recurrent seasonal variations should be taken into account in the operation, understanding and management of the WWTPs.

Low Global Diversity of Candidatus Microthrix, a Troublesome Filamentous Organism in Full-Scale WWTPs.

Candidatus Microthrix is one of the most common bulking filamentous microorganisms found in activated sludge wastewater treatment plants (WWTPs) across the globe. One species, Ca. M. parvicella, is frequently observed, but global genus diversity, as well as important aspects of its ecology and physiology, are still unknown. Here, we use the MiDAS ecosystem-specific 16S rRNA gene database in combination with amplicon sequencing of Danish and global WWTPs to investigate Ca. Microthrix spp. diversity, distribution, and factors affecting their global presence. Only two species were abundant across the world confirming low diversity of the genus: the dominant Ca. M. parvicella and an unknown species typically present along with Ca. M. parvicella, although usually in lower abundances. Both species were mostly found in Europe at low-to-moderate temperatures and their growth was favored in municipal WWTPs with advanced process designs. As no isolate is available for the novel species, we propose the name "Candidatus Microthrix subdominans." Ten high-quality metagenome-assembled genomes recovered from Danish WWTPs, including 6 representing the novel Ca. M. subdominans, demonstrated high genetic similarity between the two species with a likely preference for lipids, a putative capability to reduce nitrate and nitrite, and the potential to store lipids and poly-P. Ca. M. subdominans had a potentially more versatile metabolism including additional sugar transporters, higher oxygen tolerance, and the potential to use carbon monoxide as energy source. Newly designed fluorescence in situ hybridization probes revealed similar filamentous morphology for both species. Raman microspectroscopy was used to quantify the in situ levels of intracellular poly-P. Despite the observed similarities in their physiology (both by genomes and in situ), the two species showed different seasonal dynamics in Danish WWTPs through a 13-years survey, possibly indicating occupation of slightly different niches. The genomic information provides the basis for future research into in situ gene expression and regulation, while the new FISH probes provide a useful tool for further characterization in situ. This study is an important step toward understanding the ecology of Ca. Microthrix in WWTPs, which may eventually lead to optimization of control strategies for its growth in this ecosystem.

Mass-immigration determines the assembly of activated sludge microbial communities.

The assembly of bacterial communities in wastewater treatment plants (WWTPs) is affected by immigration via wastewater streams, but the impact and extent of bacterial immigrants are still unknown. Here, we quantify the effect of immigration at the species level in 11 Danish full-scale activated sludge (AS) plants. All plants have different source communities but have very similar process design, defining the same overall environmental growth conditions. The AS community composition in each plant was strongly reflected by the corresponding influent wastewater (IWW) microbial composition. Most species in AS across the plants were detected and quantified in the corresponding IWW, allowing us to identify their fate in the AS: growing, disappearing, or surviving. Most of the abundant species in IWW disappeared in AS, so their presence in the AS biomass was only due to continuous mass-immigration. In AS, most of the abundant growing species were present in the IWW at very low abundances. We predicted the AS species abundances from their abundance in IWW by using a partial least square regression model. Some species in AS were predicted by their own abundance in IWW, while others by multiple species abundances. Detailed analyses of functional guilds revealed different prediction patterns for different species. We show, in contrast to the present understanding, that the AS microbial communities were strongly controlled by the IWW source community and could be quantitatively predicted by taking into account immigration. This highlights a need to revise the way we understand, design, and manage the microbial communities in WWTPs.

"Candidatus Dechloromonas phosphoritropha" and "Ca. D. phosphorivorans", novel polyphosphate accumulating organisms abundant in wastewater treatment systems.

Members of the genus Dechloromonas are often abundant in enhanced biological phosphorus removal (EBPR) systems and are recognized putative polyphosphate accumulating organisms (PAOs), but their role in phosphate removal is still unclear. Here, we used 16S rRNA gene sequencing and fluorescence in situ hybridization (FISH) to investigate the abundance and distribution of Dechloromonas spp. in Danish and global wastewater treatment plants. The two most abundant species worldwide revealed in situ dynamics of important intracellular storage polymers, measured by FISH-Raman in activated sludge from four full-scale EBPR plants and from a lab-scale reactor fed with different substrates. Moreover, seven distinct Dechloromonas species were determined from a set of ten high-quality metagenome-assembled genomes (MAGs) from Danish EBPR plants, each encoding the potential for polyphosphate (poly-P), glycogen, and polyhydroxyalkanoates (PHA) accumulation. The two species exhibited an in situ phenotype in complete accordance with the metabolic information retrieved by the MAGs, with dynamic levels of poly-P, glycogen, and PHA during feast-famine anaerobic-aerobic cycling, legitimately placing these microorganisms among the important PAOs. They are potentially involved in denitrification showing niche partitioning within the genus and with other important PAOs. As no isolates are available for the two species, we propose the names Candidatus Dechloromonas phosphoritropha and Candidatus Dechloromonas phosphorivorans.

Characterizing the growing microorganisms at species level in 46 anaerobic digesters at Danish wastewater treatment plants: A six-year survey on microbial community structure and key drivers.

Anaerobic digestion (AD) is a key technology at many wastewater treatment plants (WWTPs) for converting primary and surplus activated sludge to methane-rich biogas. However, the limited number of surveys and the lack of comprehensive datasets have hindered a deeper understanding of the characteristics and associations between key variables and the microbial community composition. Here, we present a six-year survey of 46 anaerobic digesters, located at 22 WWTPs in Denmark, which is the first and largest known study of the microbial ecology of AD at WWTPs at a regional scale. For three types of AD (mesophilic, mesophilic with thermal hydrolysis pretreatment, and thermophilic), we present the typical value range of 12 key parameters including operational variables and performance parameters. High-resolution bacterial and archaeal community analyses were carried out at species level using amplicon sequencing of >1,000 samples and the new ecosystem-specific MiDAS 3 reference database. We detected 42 phyla, 1,600 genera, and 3,584 species in the bacterial community, where 70% of the genera and 93% of the species represented environmental taxa that were only classified based on MiDAS 3 de novo placeholder taxonomy. More than 40% of the bacterial species were found not to grow in the mesophilic and thermophilic digesters and were only present due to immigration with the feed sludge. Ammonium concentration was the main driver shaping the bacterial community while temperature and pH were main drivers for the archaea in the three types of ADs. Sub-setting for the growing microbes improved significantly the correlation analyses and revealed the main drivers for the presence of specific species. Within mesophilic digesters, feed sludge composition and other key parameters (organic loading rate, biogas yield, and ammonium concentration) correlated with specific growing species. This survey provides a comprehensive insight into community structure at species level, providing a foundation for future studies of the ecological significance/characteristics and function of the many novel or poorly described taxa.

Transition of microbial communities and degradation pathways in anaerobic digestion at decreasing retention time.

Tuning of operational variables is a common practice to control the anaerobic digestion process and, in advanced applications, to promote the accumulation of fermentation products. However, process variables are interrelated. In this study, the hydraulic retention time (HRT) was decoupled from the organic loading rate (OLR) in order to isolate the effect of HRT as a selective pressure on: process performance, metabolic rates (hydrolytic, acetogenic, and methanogenic) and the microbial community. Four mesophilic anaerobic digesters were subjected to a sequential decrease in HRT from 15 to 8, 4 and 2 days while keeping the OLR constant at chemical oxygen demand of 1 gCOD L r-1 d-1. The results showed that HRT alone was insufficient to washout methanogens from the digesters, which in turn prevented the accumulation of volatile fatty acids (VFA). Methanosaeta was the dominant genus in the four digesters at all HRTs. Metabolic rates showed that process performance was controlled by hydrolysis, with a clear shift in acetogenic rates, from butyrate and propionate degradation to ethanol degradation at 4 and 2d HRT. The change in acetogenic pathways was attributed to a shift in the fermentation pathways co-current with changes in fermentative bacteria. At 2d HRT, biofilm was formed on the walls and paddles of the digesters, probably as a survival strategy. Most of the taxa in the biofilm were also present in the digester media. Overall, it is the combination of HRT with other operational parameters which promotes the washout of methanogens and the accumulation of VFA.

Volatile fatty acids production from biowaste at mechanical-biological treatment plants: Focusing on fermentation temperature.

The impact of temperature (20, 35, 45, 55, 70 °C) on volatile fatty acid (VFA) production from biowaste collected at a mechanical-biological treatment plant was analysed. Additionally, relevant streams of the treatment plant were characterised to assess seasonality effects and conceive the integration of a fermentation unit. Batch fermentation tests at 35 °C showed the highest VFA yields (0.49-0.59 gCODVFA/gVS). The VFA yield at 35 °C was 2%, 6%, 10% and 14% higher than at 55, 45, 20 and 70 °C, respectively. The VFA profile was not affected by the fermentation temperature nor seasonality and was dominated by acetic, propionic and butyric acid (75-86% CODVFA). The concentration of non-VFA soluble COD and ammoniacal nitrogen in the fermentation liquor increased with temperature. The fermentation unit in the treatment plant was conceived after the pulper and hydrocyclones and before the anaerobic digester, while the fermenter temperature depends on the VFA application.

Exploring the potential of co-fermenting sewage sludge and lipids in a resource recovery scenario.

In this study, co-fermentation of primary sludge (PS) or waste activated sludge (WAS) with lipids was explored to improve volatile fatty acid production. PS and WAS were used as base substrate to facilitate lipid fermentation at 20 °C under semi-aerobic conditions. Mono-fermentation tests showed higher VFA yields for PS (32-89 mgCOD gVS-1) than for WAS (20-41 mgCOD gVS-1) where propionate production was favoured. The principal component analysis showed that the base substrate had a notable influence on co-fermentation yields and profile. Co-fermentation with WAS resulted in a greater extent of oleic acid degradation (up to 4.7%) and evidence of chain elongation producing valerate. The occurrence of chain elongation suggests that co-fermentation can be engineered to favour medium-chain fatty acids without the addition of external commodity chemicals. BMP tests showed that neither mono-fermentation nor co-fermentation had an impact on downstream anaerobic digestion.

Resource Recovery from Wastewater by Biological Technologies: Opportunities, Challenges, and Prospects.

Limits in resource availability are driving a change in current societal production systems, changing the focus from residues treatment, such as wastewater treatment, toward resource recovery. Biotechnological processes offer an economic and versatile way to concentrate and transform resources from waste/wastewater into valuable products, which is a prerequisite for the technological development of a cradle-to-cradle bio-based economy. This review identifies emerging technologies that enable resource recovery across the wastewater treatment cycle. As such, bioenergy in the form of biohydrogen (by photo and dark fermentation processes) and biogas (during anaerobic digestion processes) have been classic targets, whereby, direct transformation of lipidic biomass into biodiesel also gained attention. This concept is similar to previous biofuel concepts, but more sustainable, as third generation biofuels and other resources can be produced from waste biomass. The production of high value biopolymers (e.g., for bioplastics manufacturing) from organic acids, hydrogen, and methane is another option for carbon recovery. The recovery of carbon and nutrients can be achieved by organic fertilizer production, or single cell protein generation (depending on the source) which may be utilized as feed, feed additives, next generation fertilizers, or even as probiotics. Additionlly, chemical oxidation-reduction and bioelectrochemical systems can recover inorganics or synthesize organic products beyond the natural microbial metabolism. Anticipating the next generation of wastewater treatment plants driven by biological recovery technologies, this review is focused on the generation and re-synthesis of energetic resources and key resources to be recycled as raw materials in a cradle-to-cradle economy concept.

Metagenome from a Spirulina digesting biogas reactor: analysis via binning of contigs and classification of short reads.

BACKGROUND: Anaerobic digestion is a biological process in which a consortium of microorganisms transforms a complex substrate into methane and carbon dioxide. A good understanding of the interactions between the populations that form this consortium can contribute to a successful anaerobic digestion of the substrate. In this study we combine the analysis of the biogas production in a laboratory anaerobic digester fed with the microalgae Spirulina, a protein rich substrate, with the analysis of the metagenome of the consortium responsible for digestion, obtained by high-throughput DNA sequencing. The obtained metagenome was also compared with a metagenome from a full scale biogas plant fed with cellulose rich material. RESULTS: The optimal organic loading rate for the anaerobic digestion of Spirulina was determined to be 4.0 g Spirulina L(-1) day(-1) with a specific biogas production of 350 mL biogas g Spirulina (-1) with a methane content of 68 %. Firmicutes dominated the microbial consortium at 38 % abundance followed by Bacteroidetes, Chloroflexi and Thermotogae. Euryarchaeota represented 3.5 % of the total abundance. The most abundant organism (14.9 %) was related to Tissierella, a bacterium known to use proteinaceous substrates for growth. Methanomicrobiales and Methanosarcinales dominated the archaeal community. Compared to the full scale cellulose-fed digesters, Pfam domains related to protein degradation were more frequently detected and Pfam domains related to cellulose degradation were less frequent in our sample. CONCLUSIONS: The results presented in this study suggest that Spirulina is a suitable substrate for the production of biogas. The proteinaceous substrate appeared to have a selective impact on the bacterial community that performed anaerobic digestion. A direct influence of the substrate on the selection of specific methanogenic populations was not observed.