Skip to content
BY-NC-ND 3.0 license Open Access Published by De Gruyter Open Access April 20, 2017

Production of biogas: relationship between methanogenic and sulfate-reducing microorganisms

Ivan Kushkevych EMAIL logo , Monika Vítězová , Tomáš Vítěz and Milan Bartoš
From the journal Open Life Sciences


The production of high-quality methane depends on many factors, including temperature, pH, substrate, composition and relationship of the microorganisms. The qualitative and quantitative composition of methanogenic and sulfate-reducing microorganisms and their relationship in the experimental bioreactors has never been studied. The aim of this research was to characterize, for the first time, the diversity of the methanogenic microorganisms and sulfate-reducing bacteria, and study their relationship and biogas production in experimental bioreactors. Amplification of 16S rRNA gene fragments was carried out. Purified amplicons were paired-end sequenced on an Illumina Mi-Seq platform. The dominant morphotypes of these microorganisms in the bioreactor were homologous (99%) by the sequences of 16S rRNA gene to the Methanosarcina, Thermogymnomonas, Methanoculleus genera and Archaeon deposited in GenBank. Three dominant genera of sulfate-reducing bacteria, Desulfomicrobium, Desulfobulbus and Desulfovibrio, were detected in the bioreactor. The phylogenetic trees showing their genetic relationship were constructed. The diversity and number of the genera, production of methane, hydrogen sulfide and hydrogen in the bioreactor was investigated. This research is important for understanding the relationship between methanogenic microbial populations and other bacterial physiological groups, their substrate competition and, in turn, can be helpful for controlling methanogenesis in bioreactors.

1 Introduction

Bioenergy production from agricultural, municipal, and industrial waste is efficiently accomplished through anaerobic digestion to biogas. Biogas produced usually contains 55 – 70% of methane (CH4) and 30 – 45% of carbon dioxide (CO2) and it is stored in gas holders and subsequently used as a potential source of energy [1-4].

Anaerobic fermentation is a natural process in which microorganisms convert biodegradable substrate into biogas. It occurs in marshes, wetlands, river sediments, and also in the digestive tract of ruminants [5,6]. The microorganisms are also active in landfills where they are responsible for biodegradable waste decomposition [2,5,7]. Biogas is the final product of anaerobic metabolism. The process occurs in an anaerobic environment through the consecutive biochemical breakdown of polymers to methane and carbon dioxide [8, 9]. This is a result of the metabolism of different microorganisms which include fermentative microbes (acidogens); hydrogen-producing, acetate-forming microbes (acetogens); and methane-producing microbes (methanogens) [10-13].

Recently, the interest in anaerobic fermentation is mainly focused on its use in the economic recovery of biogas from industrial and agricultural surpluses [14-16].

On the basis of homologous sequence analysis of 16S rRNAs, methanogens have been classified into one of the three primary kingdoms of living organisms: the Archaea [10, 14, 17]. Recombinant DNA technology is one of the most powerful techniques for characterizing the biochemical and genetic regulation of methanogenesis [18]. This necessitates the selection of genetic markers, an efficient genetic transformation system, and a vector system for genetic recombination as prerequisites [19]. Anaerobic methane generation systems are known as methane bioreactors. The production of methane and growth of methanogenic microorganisms in the bioreactors depend on many factors, including temperature, pH, substrate quality, composition of specific groups of microorganisms and their accumulation of the toxic metabolic products. One of such final products of sulfate-reducing bacteria metabolism is hydrogen sulfide produced in the process of dissimilatory sulfate reduction [6, 7, 20-22].

It is known that hydrogen sulfide is toxic for living organisms and can inhibit enzymes of different groups of microorganisms [6, 20]. The sulfate-reducing bacteria can compete with the methanogens for substrate components, in this case for molecular hydrogen, and produce hydrogen sulfide in high concentration in bioreactors. It, in turn, can inhibit the growth of methanogenic microorganisms and their process of methanogenesis.

Recent progress in the molecular biology of methanogens is reviewed, new digesters are described and improvements in the operation of various types of bioreactors are also discussed. However, the prevalence of the methanogenic populations and sulfate-reducing bacteria in methane bioreactors, their diversity and relationship has never been studied properly.

The aim of this research was to characterize the diversity of the methanogenic microorganisms and sulfate-reducing bacteria using amplification of gene fragments and Illumina sequencing, and to study their relationship and production biogas in the experimental bioreactors.

2 Materials and methods

2.1 Batch anaerobic fermentation tests

A total of six laboratory batch fermenters of volume 120 liters each were used. Fermenters were equipped with a low speed agitator (operated at 60 rpm, 1 min agitation time, 60 min rest), temperature control, monitoring of pH, and sampling valve. On the first day of experiments, all fermenters were filled up with inoculum substrate and maize silage, in ratio 35:65% (v/v). Inoculum and maize silage was collected at a full scale biogas plant in Čejč, Czech Republic. This biogas plant processes as feedstock maize silage and liquid pig manure and it is operated in the mesophilic temperature regime (39°C). Inoculum characteristics were pH=7.3±0.1; DM 3.84%±0.3%; ODM 72.71%±0.4%. For maize silage, dry matter (DM) was 39.15%±0.5 and organic dry matter (ODM) was 99.21%±0.6%. Fermenters were set to 39°C ± 0.1°C. The hydraulic retention time was 26 days. During the fermentation, the quality of biogas was determined daily using X-am® 7000 gas analyzer (Drägerwerk AG & Co. KGaA Germany). Concentration of CH4, CO2, H2S and H2 in the biogas produced was determined. Quantity of biogas produced was measured by the gasometer PREMGAS BK G4 (Elster, Germany) and converted to standard conditions (T0 = 273 K, p0 = 101325 Pa).

2.2 Determination of the physical and chemical characteristics

The physical and chemical characteristics of fermenter content at each bioreactor, including temperature, pH, redox potential, dry matter content, organic dry matter content, biogas composition, sulfate content and acetate content were determined. DM content was determined by oven drying at 105±5°C followed by cooling in a desiccator and weighing until a constant weight. EcoCELL 111 (BMT Medical Technology Ltd., Czech Republic) was used in accordance with the Czech Standard Method (CSN EN 14346 2007) [23]. ODM content was determined by incineration of the samples in a muffle furnace at 550°C ± 5°C in accordance with the Czech Standard Method (CSN EN 15169 2007) [24], using a muffle furnace LMH 11/12 (LAC, Ltd., Czech Republic). The pH, and redox potential were determined by using pH/Cond meter 3320 (WTW GmbH, Germany) in accordance with the standard (CSN EN 12176 1999) [25]. Temperature of samples was determined by using the high accuracy PT100 RTD thermometer HH804U (OMEGA Engineering, INC., USA). Acetate content was determined using capillary electrophoresis Ionosep 2005 (RECMAN, Czech Republic). For determination of sulfate concentration, the spectrophotometer DR 3800 (Hach Lange GmbH, Germany) was used.

2.3 Isolation of DNA from collected samples

For DNA isolation samples from all 6 reactors were taken. These samples were combined and then examined. The QIAamp Fast DNA Stool Mini Kit (QIAGEN GmbH, Germany) is designed for rapid purification of total DNA from stool samples and was used for DNA extraction from biogas digester samples. DNA extraction was performed in accordance to the manufacturers’ instructions, with minor adjustments as described below. Briefly, 100 mg of sample was washed with 1.4 mL of ASL buffer (QIAGEN GmbH, Germany), heated at 95°C for 10 minutes. After centrifugation, an InhibitEX tablet was put in the supernatant to remove impurities and PCR inhibitors. After thorough centrifugation, 200 μL of the supernatant was added to 15 μL of proteinase K, and 200 μL of buffer AL (QIAGEN GmbH, Germany) was added. The mixture was heated to 70°C for 10 minutes and 200 μL of absolute ethanol was added. The sample mixture was then passed through the QIAamp kit column, followed by two washes with buffers AW1 and AW2 (QIAGEN GmbH, Germany). The DNA was eluted in a volume of 200 μL of elution buffer.

2.4 Amplification of gene fragments and Illumina sequencing

Amplification of 16S rRNA gene fragments was carried out using universal primers for amplification of the V3 and V4 variable regions of the 16S rRNA [26]. The primers were flanked by molecular barcodes for sample identification. The PCR reaction was prepared using Maxima™ Probe qPCR Master Mix (Thermo Fisher Scientific, USA) with cycling conditions as follows: 95°C for 10 min. followed by 30 cycles of incubation at 94°C for 30 s, 60°C for 30 s and 72°C for 120 s, and a final extension step at 72°C for 2 min. PCR products were visualized using electrophoresis on 1.5% agarose gels and purified from the gel using the QIAquick Gel Extraction Kit (QIAGEN GmbH, Germany). DNA was quantified using the Quant-iTPicoGreen dsDNA Assay (Thermo Fisher Scientific, USA) and equimolar amounts of the PCR products were pooled together.

Purified amplicons were paired-end sequenced on an Illumina Mi-Seq platform. QIIME data analysis package was used for 16S rRNA data analysis [27]. Quality filtering on raw sequences was performed in accordance to base quality score distributions, average base content per read and GC distribution in the reads. Chimeras and reads that did not cluster with other sequences were removed. The obtained sequences with qual scores higher than 20 were shortened to the same length of 350 bp and classified with RDP Seqmatch with an operational taxonomic unit (OTU) discrimination level set to 97%. The relative abundance of the taxonomic groups was calculated to the microorganisms detected in this study.

Sequences were compared using the BLAST feature of the National Center for Biotechnology Information (NCBI) [28]. The sequences were uploaded to Mega7 for comparative genomic analyses [29]. Alignments of sequences were performed in Mega7 using Clustal W and clustering was performed by the neighbor-joining method [30].

The research results were analyzed using software packages Statistica12 ( and Origin7.0 ( The basic statistic parameters (arithmetic average (M) and standard error (m) of the arithmetic average, M ± m) were calculated using the experimental data. For estimation of the validity of difference between the statistical characteristics of the data, Student’s index was calculated. The difference was valid when P<0.05 [31, 32].

3 Results

The physical and chemical characteristics of samples taken on the 14th day from the experimental bioreactors, operated at temperature 39°C±0.5, were pH (7.2±0.1), redox (-5.6±0.2), total solids (3.76%±0.3), volatile solids (70.85%±0.9). Biogas composition during the tests is shown in figure 1. Concentration of acetic acid ranged from 460 to 490 mg/L and concentration of sulfate ranged from 305 to 333 mg/L, while the sulfate concentration in fermenters on the first day of experiment ranged from 380 to 450 mg/L. The results of our research show that methane was intensively produced for the first 12 days, after this time, the methane production achieved plateau and was almost unchanged (Fig. 1). This can be explained by the fact that the methanogenic microorganisms can achieve the stationary growth phase after 12 days. As growth slowed, methane production (concentration) achieved a stable level (46.2–49.3%vol) during the 15–26 day period. An interesting pattern was observed in the production of hydrogen and hydrogen sulfide. The level of hydrogen rapidly grew the first 5 days and after that it was reduced. This can be explained by the fact that both hydrogenotrophic methanogens and sulfate-reducing bacteria consume this simple molecule. This is also evident when isolated microorganisms are examined and compared. Both acetotrophic and hydrogenotrophic metanogens were isolated. Sulfate-reducing bacteria can use hydrogen as an electron donor in the process of dissimilatory sulfate reduction. The final product of this process is hydrogen sulfide. The highest concentration of hydrogen sulfide was achieved on the 14th day, which can confirm the high number of sulfate-reducing bacteria (1050 OUT·mL-1) in the bioreactor. Different studies have found both unionized and total hydrogen sulfide concentrations important in inhibition of the SRB bacteria and methanogens. However, our results do not confirm the inhibitory effect of H2S at maximal concentration 390 ppm, as evident by the biogas production trend. The distribution of main genera in the bioreactor was investigated using amplification of 16S rRNA gene and Illumina sequencing. To clarify the genetic relationship of the methanogenic and sulfate-reducing populations of microorganisms in the bioreactor, sample of 16S rRNA gene sequences were compared with sequences of different strains from GenBank. The genomic sequences of the microorganisms are available in GenBank under accession no. KY172649, KY172662, KY172650, KY194790, KY172816, KY172822, KY172821, KY172646, KY172643, KY172641, KY172640, KY172645, KY172824, KY172823, KY126837. The obtained sequences were compared with reference strains using nucleotide Blast:Search. The sequences of the 16S rRNA gene of the methanogens were homologous (98– 99%) to genera of Methanosarcina, Thermogymnomonas, Methanoculleus, and Archaeon (Table 1). The sequences of sulfate-reducing bacteria were homologous (98–99%) to the genera of Desulfomicrobium, Desulfobulbus and Desulfovibrio (Table 2). It should be noted that most described sequences of the methanogenic strains and sulfate-reducing bacteria in GenBank are identified only to the domain, kingdom, family or genus and, in some cases, to species. Most of them are uncultured (

Fig. 1 Methane, hydrogen and hydrogen sulfide production in the biogas generated during fermentation test
Fig. 1

Methane, hydrogen and hydrogen sulfide production in the biogas generated during fermentation test

Table 1

The results of sequence analysis of the 16S rRNA gene of methanogenic microorganisms

Detected sequences and their accession number in GenBank (length of the gene fragment)Reference strains in GenBankAccession numberIdentity (%)
Sequence 1Methanosarcina mazei strain GS14-2aM 16S rRNA gene, partial sequenceKX826992.199
KY172649Methanosarcina sp. 1H1 gene for 16S rRNA, partial sequenceLC170394.199
(424 bp)Uncultured Methanosarcina sp. clone T6190SA18-18 16S rRNA gene, partial sequenceKU355742.199
Sequence 2Uncultured archaeon clone g10-56 16S rRNA gene, partial sequenceJX576125.198
KY172662Uncultured Thermoplasmata archaeon clone g8-4 16S rRNA gene, partial sequenceJX576112.198
(422 bp)Uncultured Methanomassiliicoccus sp. clone LZNG25 16S rRNA gene, partial sequenceJX456453.198
Sequence 3Uncultured Thermogymnomonas sp. partial 16S rRNA gene, isolate OTU_11LT624815.199
KY172650Uncultured Thermogymnomonas sp. partial 16S rRNA gene, isolate OTU_4LT624808.199
(422 bp)Uncultured archaeon partial 16S rRNA gene, clone AKA055LN874207.199
Sequence 4Uncultured Thermoplasmatales archaeon partial 16S rRNA gene, isolate OTU_9LT624813.199
KY194790Uncultured Thermoplasmatales archaeon partial 16S rRNA gene, isolate OTU_6LT624810.199
(423 bp)Uncultured archaeon clone WWA-D10 16S rRNA gene, partial sequenceKM870439.199
Sequence 5Methanoculleus sp. strain Biowerk_c-HAW 16S rRNA gene, partial sequenceKX619406.199
KY172816Methanoculleus bourgensis isolate BA1 genome assembly, chromosome: ILT549891.199
(420 bp)Methanoculleus sp. MAB1 isolate MAB1 genome assembly, chromosome: chrILT158599.199
Table 2

The results of sequence analysis of the 16S rRNA gene of sulfate-reducing bacteria

Detected sequences and their accession number in GenBank (length of the gene fragment)Reference strains in GenBankAccession numberIdentity (%)
Sequence 1Desulfovibrio desulfuricans strain E4 16S rRNA, partial sequenceKJ459863.199
KY172822Desulfovibrio desulfuricans strain E2 16S rRNA, partial sequenceKJ459861.199
(465 bp)Uncultured bacterium partial 16S rRNA gene, clone 48h8HG531899.199
Sequence 2Uncultured Desulfovibrio sp. clone MFC-2-L19 16S rRNA gene, partial sequenceJX944554.199
KY172821Uncultured bacterium clone MFC4P_127 16S rRNA gene, partial sequenceJF309175.199
(465 bp)Desulfovibrio simplex strain JCM 16812 16S rRNA gene, partial sequenceNR_113296.199
Sequence 3Desulfobulbus propionicus strain DSM 2032 16S rRNA gene, complete sequenceNR_074930.199
KY172646Desulfobulbus propionicusDSM 2032, complete genomeCP002364.199
(466 bp)Uncultured bacterium clone B01 16S rRNA gene, partial sequenceEU136253.199
Sequence 4Desulfobulbus sp. canine oral taxon 078 clone OC011 16S rRNA gene, partial sequenceJN713241.196
KY172643Uncultured Desulfobulbus sp. partial 16S rRNA gene, isolate OTU 265LT625152.196
(466 bp)Uncultured Desulfobulbus sp. partial 16S rRNA gene, isolate OTU 199LT625082.196
Sequence 5Uncultured Desulfomicrobium sp. gene for 16S rRNA, partial sequence, clone:3CP(-)_3B908615.196
KY172641Uncultured Desulfomicrobium sp. gene for 16S rRNA, partial sequence, clone:3CP(-)_2AB908614.196
(466 bp)Uncultured Desulfomicrobium sp. gene for 16S rRNA, partial sequence,clone: 3CP(+)_10AB908535.196
Sequence 6Desulfobulbus sp. Prop6 16S rRNA gene, partial sequenceKU845305.199
KY172640Uncultured bacterium clone 02d01 16S rRNA gene, partial sequenceGQ138500.199
(466 bp)Uncultured bacterium clone 04c04 16S rRNA gene, partial sequenceGQ134634.199
Sequence 7Uncultured Desulfobulbus sp. partial 16S rRNA gene, isolate OTU 265LT625152.194
KY172645Uncultured Desulfobulbus sp. partial 16S rRNA gene, isolate OTU 199LT625082.194
(466 bp)Desulfobulbus sp. oral taxon 041 clone WWP_SS1_G06 16S rRNA gene, partial sequenceEU398208.194
Sequence 8Uncultured sulfate-reducing bacterium clone 2R1V11 16S rRNA gene, partial sequenceEF592786.193
KY172824Uncultured delta proteobacterium clone 2R1U31 16S rRNA gene, partial sequenceEU104788.193
(466 bp)Uncultured sulfate-reducing bacterium clone 2R1V05 16S rRNA gene, partial sequenceEF592783.192
Sequence 9Uncultured Desulfovibrionales bacterium clone Flu2 26 16S rRNA gene, partial sequenceJQ701289.184
KY172823Desulfovibrio sp. S4 gene for 16S rRNA, partial sequenceLC186051.183
(443 bp)Uncultured Desulfovibrio sp. gene for 16S rRNA, partial sequence, clone: LR333B-40LC001349.183
Sequence 10Uncultured Desulfovibrio sp. clone MFC-2-L19 16S rRNA gene, partialJX944554.195
KY126837Desulfovibrio simplex strain JCM 16812 16S rRNA gene, partial sequenceNR_113296.195
(465 bp)Desulfovibrio intestinalis partial 16S rRNA gene, strain JG-G12AJ295680.195

Based on all of the 16S rRNA gene sequences of methanogens and sulfate-reducing bacteria from the bioreactor, a phylogenetic tree demonstrating their genetic relationship was built (Fig. 2, 3). The detected genera of methanogens were homologous with Methanosarcina mazei strain GS14-2aM, uncultured archaeon clone g10-56, uncultured hermogymnomonas p. isolate OTU_11, ncultured hermoplasmatales rchaeon isolate OTU_9, Methanoculleus sp. strain Biowerk_c-HAW (Fig. 2). Another phylogenetic tree was built for 16S rRNA gene sequences of sulfate-reducing bacteria which were homologous with Desulfovibrio desulfuricans strain E4 16S, uncultured Desulfovibrio sp. clone MFC-2-L19, Desulfobulbus propionicus strain DSM 2032, Desulfobulbus sp. clone OC011, uncultured esulfomicrobium p. clone: 3CP(-)_3, Desulfobulbus sp. Prop6, uncultured esulfobulbus p. isolate OTU 265, uncultured sulfate-reducing bacterium clone 2R1V11, uncultured Desulfovibrionales bacterium clone Flu2_26 (Fig. 3).

Fig. 2 Phylogenetic tree of relationship sequences of 16S rRNA gene of the methanogenic populations in methane bioreactor
Fig. 2

Phylogenetic tree of relationship sequences of 16S rRNA gene of the methanogenic populations in methane bioreactor

Fig. 3 Phylogenetic tree of relationship sequences of 16S rRNA gene of the sulfate-reducing bacteria in methane bioreactor
Fig. 3

Phylogenetic tree of relationship sequences of 16S rRNA gene of the sulfate-reducing bacteria in methane bioreactor

Percentage ratio of methanogenic and sulfate-reducing microorganisms was calculated by OUT·mL-1 determined from Illumina sequencing (Fig. 4). Results indicate that three genera of both physiological groups of microorganisms were detected in the experimental bioreactors. The dominant genus of methanogens was Methanosarcina, which was 62% of all detected methanogens. Two other genera, Thermogymnomonas and Methanoculleus, were 16 and 2%, respectively, as well as other Archaeons which were not identified to the genus (20%). For sulfate-reducing bacteria, three genera were identified: Desulfomicrobium, Desulfobulbus and Desulfovibrio in percentage ratios of 48, 39.8, and 12%, respectively. Other sulfate-reducing bacteria (0.19%) were not identified to the genus.

Fig. 4 Qualitative and percentage composition of methanogenic and sulfate-reducing microorganisms
Fig. 4

Qualitative and percentage composition of methanogenic and sulfate-reducing microorganisms

4 Discussion

Methane is the final product of anaerobic metabolism carried out by communities of hydrolytic bacteria and fungi, acid-producing intermediary organisms, and finally, methanogenic microorganisms [17]. Methane-producing communities are very stable and resilient, but they are also complex and largely undefined. The results of our studies are consistent to other research described in literature [5, 14, 33, 34]. Production of methane depends on many factors, including physical (temperature) and chemical (pH), type of substrate, its concentration and accessibility for microorganisms, composition and the ratio of the microorganisms and their metabolic compounds in the bioreactors. Therefore, bioreactors are a complex fermentative system which includes various oxidations and reduction processes with changes in redox potential. Our research demonstrates that the intensive production of methane in the experimental bioreactor lasts the first 12 days and correlates with the accumulation of hydrogen sulfide. The highest concentration of hydrogen sulfide (390 ppm) was detected on the 14th day and it correlates with the titer of sulfate-reducing bacteria (1050 OUT·mL-1). The sulfate-reducing bacteria (SRB) are a heterogeneous group of microorganisms which use sulfate as an electron acceptor in the process of dissimilatory sulfate reduction [35, 36]. The final product of this process is hydrogen sulfide [6, 20, 22, 37]. For sulfate reduction, SRB need exogenous electron donors, such as: organic compounds (e.g., lactate, propionate, butyrate, ethanol, etc.) and molecular hydrogen. Organic compounds for SRB can be simultaneously electron donors and carbon sources and oxidized completely (to CO2) or incompletely (to acetate and CO2) [6, 7]. SRB, oxidizing organic compounds incompletely, belong to the group called “Acetogenic sulfate-reducing bacteria” [7]. Detected SRB, Desulfovibrio(48%), Desulfomicrobium(39.8%), and Desulfobulbus (12%) genera, are acetogenic microorganisms which oxidize organic compounds, incompletely, to acetate and CO2. Produced acetate is consumed by methanogens e.g., species of the Methanosarcina genus, which are dominant microorganisms (62%) in the bioreactor (see Fig. 4). Our results demonstrated the presence of acetotrophic and hydrogenotrophic methanogenic archaea. The species of the Methanosarcina genus can form multicellular colonies and are anaerobic methanogens. They are widespread in the rumen of cows, sheep, goats, deer, and the large intestine of humans [7]. Recently, there has been a study on M. barkeri, because this species has the enzyme methylamine methyltransferase, which catabolizes methylamines leading to methane production. Methanosarcina sp. possess all three known pathways for methanogenesis, and can utilize a broad spectrum of substrates, including hydrogen. All the other methanogens can utilize no more than two methanogenic substrates and possess a single pathway for methanogenesis [38]. It also has a number of distinct morphological forms, including single cells with and without a cell envelope, as well as multicellular packets and lamina [39].

The other dominant genus of microbial methanogens in the studied bioreactor was Thermogymnomonas (16%) . The species of this genus were also isolated as a novel thermoacidophilic, cell wall-less archaeon from a solfataric field in Ohwaku-dani, Hakone, Japan. The cells were irregular cocci, sometimes lobed, club-shaped or catenated, and were highly variable in size, ranging from 0.8 to 8.0 µm in diameter [33]. Itoh et al. (2007) identified this strain as Thermogymnomonas acidicola. The strain grew at temperatures in the range 38–68°C (optimally at 60°C) and at pH 1.8–4.0 (optimally at around pH 3.0). Strain IC-189T was an obligate anaerobic and heterotrophic microorganism, requiring yeast extract for growth. Yeast extract, glucose and mannose served as carbon and energy sources. Therefore, strain IC-189T represents a novel genus (order Thermoplasmatales) and species [33].

The Methanoculleus genus was found in the lowest number among all detected methanogens. The species of the genus were frequently described as playing an important role in different biogas reactor systems [34, 40, 41]. Methanoculleus bourgensis was always detected as the dominant in biogas systems. The prevalence of M. bourgensis in reactors performing syntrophic acetate oxidation under high ammonium concentrations [42, 43, 44], indicates the importance of this methanogen in corresponding communities. Isolation and/or cultivation of M. bourgensis, together with acetate-oxidizing bacteria [45] such as Clostridium ultunense [46], led to the assumption that syntrophic association may play an important role for members of the Methanoculleus genus [42, 47]. The 16S rRNA gene sequence analysis classified the isolate as a member of the species M. bourgensis with 99% sequence identity to the 16S rRNA gene of strain MS2T [48, 49]. Genomic DNA of strain BA1 was isolated and sequenced, applying the paired-end protocol on an Illumina MiSeq system [34, 48]. In our studies, using Illumina sequencing, the sequences belonging to this genus were the most often detected in all bioreactors.

The detected genera of methanogens and SRB described can compete by molecular hydrogen (Fig. 5). Hydrogen and CO2 can be used by the methanogens for their growth and methane production and, simultaneously, SRB can also use H2 as a simple electron donor. Therefore, competition for molecular hydrogen may occur between methanogens and SRB. Accordingly, decrease of H2 was observed after the 6th day. Perhaps, the highest consumption of hydrogen by both microorganism groups occurred during this period because afterwards hydrogen levels stabilized. However, SRB produce hydrogen sulfide which can be toxic for methanogens, and may inhibit methanogenesis.

Fig. 5 The hypothetical scheme of the relationship between methanogenic and sulfate-reducing microorganisms
Fig. 5

The hypothetical scheme of the relationship between methanogenic and sulfate-reducing microorganisms

5 Conclusions

We conclude that studies of the diversity of methanogenic microorganisms under the influence of various factors in the bioreactors require further understanding of the process of biogas production. The sulfate-reducing bacteria can compete with the hydrogenotrophic methanogens for substrate components, in this case molecular hydrogen, and produce hydrogen sulfide in high concentrations in the bioreactors. This, in turn, can inhibit the growth of methanogenic microorganisms and their process of methanogenesis. It can cause unbalance of other microbial communities and instability of the fermentation process. However, this was not tested in our experiments at a maximum concentration of H2S (390 ppm). The study of new isolates is still important for selecting the optimal conditions for methane production process.


This study was supported by Masaryk University (project TAČR GAMA – internal project CTT MU “Technology for qualitative biogas treatment” Registration Project ID: 51047). The authors wish to acknowledge the institutional support of the Faculty of AgriSciences, Mendel University in Brno funded by the Ministry of Education, Youth and Sports of the Czech Republic. The authors express their gratitude to Igor Starunko, Head of Editorial Office of Studia Biologica Journal at Ivan Franko National University of Lviv for his help in the preparation of graphic material.

  1. Conflict of Interest: Authors state no conflict of interest


[1] Scherer P.A., Vollmer G.R., Fakhouri T., Martensen S., Development of methanogenic process to degrade exhaustively the organic fraction of municipal grey waste under thermophilic and hyperthermophilic conditions, Water Sci. Technol., 2000, 41, 83-9110.2166/wst.2000.0059Search in Google Scholar

[2] Schink B., Energetics of syntrophic cooperation in methanogenic degradation, Microb. Mol. Biol. Rev., 1997, 61, 262-28010.1128/mmbr.61.2.262-280.1997Search in Google Scholar

[3] Weiland P., Biogas production: current state and perspectives, Appl. Microbiol. Biotechnol., 2010, 85, 849-86010.1007/s00253-009-2246-7Search in Google Scholar

[4] Wilkie A., Biomethane from Biomass, Biowaste, and Biofuels, 195-205. In Wall J., Harwood C., Demain A. (ed.), Bioenergy. ASM Press, Washington, 200810.1128/9781555815547.ch16Search in Google Scholar

[5] Krich K., Augenstein D., Batmale J.P., Benemann J., Rutledge B., Salour D., Biomethane from Dairy Waste. A Sourcebook for the Production and Use of Renewable Natural Gas in California. USDA Rural Development, 2005Search in Google Scholar

[6] Kushkevych I., Dissimilatory sulfate reduction in the intestinal sulfate-reducing bacteria, Studia Biologica, 2016, 10(1), 197-22810.30970/sbi.1001.560Search in Google Scholar

[7] Barton L.L., Hamilton W.A., Sulphate-Reducing Bacteria. Environmental and Engineered Systems. Cambridge University Press, 2010, 55310.1017/CBO9780511541490Search in Google Scholar

[8] Ahring B., Ibrahim A.A., Mladenovska Z., Effect of temperature increase from 55 to 65°C on performance and microbial population dynamics of an anaerobic reactor treating cattle manure. Water Resour., 2001, 35, 2446-245210.1016/S0043-1354(00)00526-1Search in Google Scholar

[9] Ziemiński K., Frąc M., Methane fermentation process as anaerobic digestion of biomass: Transformations, stages and microorganisms. African J. Biotech., 2012, 11(18), 4127-413910.5897/AJBX11.054Search in Google Scholar

[10] Conrad R., Contribution of hydrogen to methane production and control of hydrogen concentration in methanogenic soils and sediments. FEMS Microbiol. Ecol., 1999, 28, 193-20210.1111/j.1574-6941.1999.tb00575.xSearch in Google Scholar

[11] Demirel B., Scherer P., The roles of acetotrophic and hydrogenotrophic methanogens during anaerobic conversion of biomass to methane: a review. Rev. Environ. Sci. Biotechnol., 2008, 7, 173-19010.1007/s11157-008-9131-1Search in Google Scholar

[12] Jackel U., Thummes K., Kampfer P., Thermophilic methane production and oxidation in compost. FEMS Microbiol. Ecol., 2005, 52, 175-18410.1016/j.femsec.2004.11.003Search in Google Scholar

[13] Santosh Y., Sreekrishnan T.R., Kohli S., Rana V., Enhancement of biogas production from solid substrates using different techniques,Bioresour. Technol., 2004, 95, 1-1010.1016/j.biortech.2004.02.010Search in Google Scholar

[14] Griffin M.E., McMahon K.D., Mackie R.I., Raskin L., Methanogenic population dynamics during start-up of anaerobic digesters treating municipal soild waste and biosolids. Biotechnol. Eng., 2000, 57, 342-35510.1002/(SICI)1097-0290(19980205)57:3<342::AID-BIT11>3.0.CO;2-ISearch in Google Scholar

[15] Grothenhuis J.T., Smith M., Plugge C.M., Yuansheng X., Lammeren A.A., Stams A.J., et al., Bacteriological composition and structure of granular sludge adapted to different substrates. Appl. Environ. Microbiol., 1991, 57, 1942-194910.1128/aem.57.7.1942-1949.1991Search in Google Scholar

[16] Ilyin V.K., Korniushenkova I.N., Starkova L.V., Lauriniavichius K.S., Study of methanogenesis during bioutilization of plant residuals. Acta Astronautica, 2005, 56, 465-47010.1016/j.actaastro.2004.05.077Search in Google Scholar

[17] Zeikus J.G., The biology of methanogenic bacteria, Bact. Rev., 1977, 41, 514-54110.1128/br.41.2.514-541.1977Search in Google Scholar

[18] Amon T., Amon B., Kryvoruchko V., Zollitsch W., Mayer K., Gruber L., Biogas production from maize and dairy cattlemanure – influence of biomass composition on the methane yield, Agric. Ecosys.Environ.,2007,118, 173-18210.1016/j.agee.2006.05.007Search in Google Scholar

[19] Chynoweth D.P., Turick C.E., Owens J.M., Jerger D.E., Peck M.W., Biochemical methane potential of biomass and waste feedstocks, Biomass Bioen.,1993, 5, 95-11110.1016/0961-9534(93)90010-2Search in Google Scholar

[20] Kushkevych I.V., Effect of hydrogen sulfide at differential concentrations on the process of dissimilatory sulfate reduction by the bacteria Desulfovibrio piger, Sci Notes of Ternopil Nat. Ped. Univ. Series Biol, 2013, 4(57), 74-80Search in Google Scholar

[21] Kushkevych I.V., Dissimilatory sulfate reduction by various Desulfovibrio sp. strains of the human intestine, Microbiol. and Biotechnol., 2013, 3(23): 50-6310.18524/2307-4663.2013.3(23).48939Search in Google Scholar

[22] Kushkevych I.V., Dissimilatory sulfate reduction in bacterium Desulfovibrio piger Vib-7 under the effect of medium with differential acidity, American J. Microbiol. & Biotechnol, 2014, 1(2), 49-55.10.2174/1874285801509010055Search in Google Scholar

[23] Czech Standards Institute. Characterization of waste – Calculation of dry matter by determination of dry residue or water content ČSN EN 14346, 2007Search in Google Scholar

[24] Czech Standards Institute. Characterization of waste – Determination of loss on ignition in waste, sludge and sediments ČSN EN 15169, 2007Search in Google Scholar

[25] Czech Standards Institute. Characterization of sludge – Determination of pH-value ČSN EN 12176, 1999Search in Google Scholar

[26] Nossa C.W., Oberdorf W.E., Yang L., Aas J.A., Paster B.J., Desantis T.Z., et al., Design of 16S rRNA gene primers for 454 pyrosequencing of the human foregut microbiome, World J. Gastroenterol., 2010, 16(33), 4135-414410.3748/wjg.v16.i33.4135Search in Google Scholar

[27] Caporaso J.G., Kuczynski J., Stombaugh J., Bittinger K., Bushman F.D., Costello E.K., et al., QIIME allows analysis of high-throughput community sequencing data. Nat. Methods, 2010, 7(5), 335-33610.1038/nmeth.f.303Search in Google Scholar

[28] Altschul S.F., Gish W., Miller W., Myers E.W., Lipman D.J., Basic local alignment search tool. J. Mol. Biol., 1990, 215(3), 403-1010.1016/S0022-2836(05)80360-2Search in Google Scholar

[29] Kearse M., Moir R., Wilson A., Stones-Havas S., Cheung M., Sturrock S., et al., Geneious Basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics, 2012, 28(12), 1647-910.1093/bioinformatics/bts199Search in Google Scholar PubMed PubMed Central

[30] Larkin M.A., Blackshields G., Brown N.P., Chenna R., McGettigan P.A., McWilliam H., et al., Clustal W. and Clustal X version 2.0. Bioinformatics, 2007, 23(21), 2947-810.1093/bioinformatics/btm404Search in Google Scholar PubMed

[31] Chen P.Y., Popovich P.M., Correlation: Parametric and Nonparametric Measures. Sage University Papers Series on Quantitative Applications in the Social Sciences, 200210.4135/9781412983808.n1Search in Google Scholar

[32] Bailey N.T.J., Statistical Methods in Biology, third ed. Cambridge University Press, Cambridge, 199510.1017/CBO9781139170840Search in Google Scholar

[33] Itoh T., Yoshikawa N., Takashina T., Thermogymnomonasacidicola gen. nov., sp. nov., a novel thermoacidophilic, cell wall-less archaeon in the order Thermoplasmatales, isolated from a solfataric soil in Hakone, Japan, Int. J. Syst. Evol. Microbiol., 2007, 57(11), 2557-6110.1099/ijs.0.65203-0Search in Google Scholar PubMed

[34] Maus I., Wibberg D., Winkler A., Pühler A., Schnürer A., Schlütera A., Complete Genome Sequence of the Methanogen Methanoculleus bourgensis BA1 Isolated from a Biogas Reactor Genome, Announcements, 2016, 4(3) e00568-1610.1128/genomeA.00568-16Search in Google Scholar PubMed PubMed Central

[35] Kushkevych I.V., Dose-dependent effect of electron acceptor and donor on dissimilatory sulfate reduction by bacteria Desulfovibrio piger Vib-7 of human intestine, Studia Biologica, 2014, 8(1), 103-11610.30970/sbi.0801.313Search in Google Scholar

[36] Kushkevych I.V., Growth of the Desulfomicrobium sp. strains, their sulfate- and lactate usage,production of sulfide and acetate by the strains isolated from the human large intestine, Microbiol. Discovery, 2014, 2, 1-810.7243/2052-6180-2-1Search in Google Scholar

[37] Kushkevych I., Bolis M., Bartoš M., Model-based characterization of the kinetic parameters of dissimilatory sulfate reduction under the effect of different initial density of Desulfovibrio piger Vib-7 bacterial cells, The Open Microbiol. J., 2015, 9, 55-6910.2174/1874285801509010055Search in Google Scholar

[38] Ferry J., Methanogenesis: Ecology, Physiology, Biochemistry & Genetics, Chapman & Hall Inc, New York, 199310.1007/978-1-4615-2391-8Search in Google Scholar

[39] Conway de Macario E., Maeder D. L., Macario A.J.L., Breaking the mould: archaea with all four chaperoning systems, Biochem. and Biophys. Res. Commun., 2003, 301, 811-81210.1016/S0006-291X(03)00047-0Search in Google Scholar

[40] Jaenicke S., Ander C., Bekel T., Bisdorf R., Dröge M., Gartemann K.H. et al., Comparative and joint analysis of two metagenomic datasets from a biogas fermenter obtained by 454-pyrosequencing, PLoS One, 2011, 6, e1451910.1371/journal.pone.0014519Search in Google Scholar PubMed PubMed Central

[41] Stolze Y., Zakrzewski M., Maus I., Eikmeyer F., Jaenicke S., Rottmann N.,et al., Comparative metagenomics of biogas-producing microbial communities from production-scale biogas plants operating under wet or dry fermentation conditions, Biotechnol. Biofuels, 2015, 8, 1410.1186/s13068-014-0193-8Search in Google Scholar PubMed PubMed Central

[42] Westerholm M., Levén L., Schnürer A., Bioaugmentation of syntrophic acetate-oxidizing culture in biogas reactors exposed to increasing levels of ammonia, Appl. Environ. Microbiol., 2012, 78, 7619-762510.1128/AEM.01637-12Search in Google Scholar PubMed PubMed Central

[43] Westerholm M., Müller B., Isaksson S., Schnürer A., Trace element and temperature effects on microbial communities and links to biogas digester performance at high ammonia levels, Biotechnol. Biofuels, 2015, 8, 15410.1186/s13068-015-0328-6Search in Google Scholar PubMed PubMed Central

[44] Moestedt J., Müller B., Westerholm M., Schnürer A., Ammoniathreshold for inhibition of anaerobic digestion of thin stillage and the importance of organic loading rate, Microb. Biotechnol., 2016, 9, 180-19410.1111/1751-7915.12330Search in Google Scholar PubMed PubMed Central

[45] Fotidis I.A., Karakashev D., Angelidaki I., Bioaugmentation with anacetate-oxidising consortium as a tool to tackle ammonia inhibition ofanaerobic digestion, Bioresour Technol., 2013, 146, 57-6210.1016/j.biortech.2013.07.041Search in Google Scholar PubMed

[46] Schnürer A., Schink B.H., Svensson B.H., Clostridium ultunense sp. nov., a mesophilic bacterium oxidizing acetate in syntrophic association with a hydrogenotrophic methanogenic bacterium, Int. J. Syst. Bacteriol., 1996, 46, 1145-115210.1099/00207713-46-4-1145Search in Google Scholar PubMed

[47] Fotidis I.A., Wang H., Fiedel N.R., Luo G., Karakashev D.B., Angelidaki I., Bioaugmentation as a solution to increase methane production from an ammonia-rich substrate, Environ. Sci. Technol., 2014, 48, 7669-767610.1021/es5017075Search in Google Scholar PubMed

[48] Maus I., Wibberg D., Stantscheff R., Stolze Y., Blom J., EikmeyerF.G.,et al., Insights into the annotated genome sequence of Methanoculleus bourgensis MS2 (T), related to dominant methanogens in biogas-producing plants, J. Biotechnol., 2014, 201, 43-5310.1016/j.jbiotec.2014.11.020Search in Google Scholar PubMed

[49] Wibberg D., Blom J., Jaenicke S., Kollin F., Rupp O., Scharf B., et al., Complete genome sequencing of Agrobacterium sp.H13-3, the former Rhizobium lupine H13-3, reveals a tripartite genome consisting of a circular and a linear chromosome and an accessory plasmid but lacking a tumor-inducing Ti-plasmid, J. Biotechnol., 2011, 155, 50-6210.1016/j.jbiotec.2011.01.010Search in Google Scholar PubMed

Received: 2016-11-11
Accepted: 2017-3-1
Published Online: 2017-4-20

© 2017 Ivan Kushkevych et al.

This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 3.0 License.

Downloaded on 26.11.2022 from
Scroll Up Arrow