Next Article in Journal
Reconstitution in Proteoliposomes of the Recombinant Human Riboflavin Transporter 2 (SLC52A2) Overexpressed in E. coli
Previous Article in Journal
Independent Component Analysis for Unraveling the Complexity of Cancer Omics Datasets
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

1
Department of Biochemistry and Environmental Chemistry, Institute of Biotechnology, The John Paul II Catholic University of Lublin, Konstantynów Street 1 I, 20-708 Lublin, Poland
2
Department of Agricultural Microbiology, Institute of Soil Science and Plant Cultivation–State Research Institute (IUNG-PIB), Czartoryskich Street 8, 24-100 Puławy, Poland
3
Institute of Agrophysics, Polish Academy of Sciences, Doświadczalna Street 4, 20-290 Lublin, Poland
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2019, 20(18), 4415; https://doi.org/10.3390/ijms20184415
Submission received: 15 July 2019 / Revised: 5 September 2019 / Accepted: 6 September 2019 / Published: 8 September 2019
(This article belongs to the Section Molecular Microbiology)

Abstract

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

1. Introduction

Microbial diversity and functioning in the environment are one of the most intriguing issues in science, given that they are related to the foundations of the knowledge regarding biogeochemical cycles. Recent methodical advances, in particular the introduction of next generation sequencing, facilitate the description of microbiomes with unprecedented depth [1,2]. It is very important to recognize microbial communities in the natural environment in terms of their composition and function, as achievements in this field not only contribute to overall knowledge but also provide a basis for biotechnological progress [3,4]. One of the most common objects of this type of research are anoxic environments, which are a source of microbiota that can be employed in biogas production due to their physiological capabilities [5,6]. Methanogenesis is a multistage process in which degradation of complex organic molecules into simple one-carbon compounds is performed by a consortium of microorganisms [7,8]. Their existence in the natural environment is an effect of coevolution leading to the development of metabolic interactions that allow the flow of carbon, energy, and other intermediates for mutual benefit [9,10]. The best recognized methanogenic Archaea belong to the phylum Euryarchaeota and are classified into seven orders, namely Methanococcales, Methanobacteriales, Methanosarcinales, Methanomicrobiales, Methanopyrales, Methanocellales, and Methanomassiliicoccales [11]; yet, additional new phyla, namely Bathyarchaeota [12] and Verstraetearchaeota [13], have been postulated. Methanogens use a narrow range of substrates such as H2/CO2, formate, acetate, methanol, methylated compounds, and CO. Ethanol and 2-propanol may be used by some methanogens instead of H2 as electron donors [11]. This narrow range of substrates is delivered by bacteria, mostly from the genera Clostridium, Pseudomonas, Desulfovibrio, and Cellulomonas [14]. Organic matter, in freshwater lakes, is primarily decomposed by bacteria to acetate, CO2, and H2, which are the main substrates used by methanogens in two major pathways: hydrogenotrophic and acetotrophic [15,16].
The biotechnological potential of anaerobic microbial communities is closely connected with their biodiversity [17]. The highest biodiversity has been found in transitional zones called ecotones [18]. The observations that were first made on the macroscale (for higher plants and animals) [19] revealed the highest biodiversity at the boundaries of forests and meadows [20], or lakes and terrestrial ecosystems [21], and have been confirmed on the microscale as well [22]. Microbial biodiversity (in terms of the structure and function) in zones located between forests and grasslands [22], lakes and rivers [23,24], was found to be higher than in adjacent ecosystems. Such locations may thus have a hidden potential for biotechnology. The present study was focused on a microbial community found in the sediment of a shallow reservoir, Szczecin, that has developed over the last three decades in subsidence resulting from mining activity. Due to the continuous exploitation of coal, the area and depth of the reservoir are increasing continuously. It is also exposed to huge temperature and insolation fluctuations resulting from the climatic conditions of the mid latitudes. Consequently, the sediment is subjected to continuous changes in temperature and aeration. Furthermore, the submergence of subsequent fragments of agricultural soils leads to the release of huge amounts of biogenic compounds, and an additional portion thereof is transported from adjacent fields by surface runoff. The imbalance in the ecosystem is manifested by algal blooms occurring in the reservoir each summer. Oxygen deficiencies and organic matter supply from decaying algae create conditions that support biogenic methane (CH4) formation, which is a result of the cooperation of methanogens with fermentative microorganisms.
Generally, the biodiversity and metabolic activities of bacterial communities in lake sediments decrease with the depth in the profile [25]. Hence, the surface layer of the bottom sediments of the Szczecin reservoir was chosen to study microbial ecology and a possible source of microorganisms that can potentially be used in biogas production.
Sediment-derived enrichment cultures grown on various carbon substrates were used as a tool to elucidate the relationship between microbial community composition and methane production rate. Cultivation largely reduces the initial biodiversity. On the other hand, targeted selection of culture medium components may lead to the enrichment of microorganisms exhibiting the desired abilities and provide new insight into their ecophysiology. The complexity and multistage nature of the biological CH4 formation process suggest that investigations are only valid when the community structure is studied using high throughput methods.
The specific aim of this study was to determine the link between the microbiome composition and methane production and to describe changes occurring in the methanogenic community structure and its metabolic capabilities as a function of the enrichment stage and medium composition.

2. Results

The investigated sediments were characterized by slightly acidic reaction (5.46), redox potential (Eh) below −265 mV, and organic carbon content equal to 48 g kg d.w.−1. The content of N-NH4 forms was almost three times higher than the concentration of bioavailable phosphorus (P-PO4) (Table 1).

2.1. Methane Production

CH4 production in the sediment samples (tested in laboratory conditions at temperatures 10–40 °C) revealed a maximum of over 1.67 mg CH4 L−1 d−1, at 30 °C. The enrichment cultures were prepared in the same temperature range and in mineral media supplemented with different carbon sources. In the medium H(+), the highest methane production rate (MP, above 2.99 mg CH4 L−1 d−1, at 20 and 30 °C) was determined during stage S1. In the subsequent stages S2 (20, 30 °C) and S3 (30 °C), the activities decreased and were significantly lower, ca. 1.28 (20, 30 °C) and 1.12 mg CH4 L−1 d−1 (30 °C), respectively. The supplementation of the medium with sodium acetate (H(+)acet) resulted in MP of approximately 5.54 mg CH4 L−1 d−1 in the S1 stage, whereas one third of that value was noted in S2 and S3. Substantially higher MP was revealed when the headspace gas in the serum vials was replaced by methanogenic substrates: CO2 and H2 (20:80 v/v). In this case, the optimal temperature was 30 °C. MP increased with each subsequent stage of incubation. It was the highest in S3 (exceeding 12.8 mg CH4 L−1 d−1). This value was also the highest among all the experimental treatments. When the medium was free of tryptone and yeast extract but supplemented with sodium acetate (H(−)acet), methane production above 3.8 mg CH4 L−1 d−1 was observed at 20–40 °C in S1. In S2, MP decreased significantly, maximally by 83% with optimal thermal conditions noted at 20 °C. At this temperature, the CH4 production rate reached 1.74 mg CH4 L−1 d−1, and the reduction of MP was the lowest but still substantial, i.e., 57%, compared to S1. In S3, the value of MP was similar to that in S2. In the H(−) medium and the CO2 and H2 atmosphere (20:80 v/v), MP was the highest at 30 and 40 °C in S1. However, in the next stage (S2), 30 °C appeared to be more suitable for methanogenesis, with MP reaching as much as 10.23 mg CH4 L−1 d−1. In S3, MP was still high, i.e., 8.85 mg CH4 L−1 d−1 (Figure 1).
The changes in the duration of the lag phase of MP indicate that the microbial consortia that developed in all of the treatments are able to metabolize the substrates added and create suitable conditions for methanogenesis. The lag phase was the longest in the initial phase of the experiment (S1) and shortened after the subsequent transfers (Table 2). The rise in temperature also caused a reduction in the lag phase. Starting from S2, in almost all experimental treatments, at 30 and 40 °C, methane production began within one day of inoculation. The only exception was the community grown with acetate as the sole carbon source (H(−)acet). This community was not only characterized by the longest lag time but also saw optimal growth at 20 °C.

2.2. Physicochemical Conditions

The pH value at the beginning of the incubation in S1 was generally slightly higher than at the end (Table 3). In the subsequent transitions, changes in the reaction depended on the treatment applied. Generally, H2 and CO2 resulted in slight acidification of the medium, while addition of organic carbon substrates promoted alkalization (Table 3). Based on the oxidation-reduction potential, it was found that the microbial communities were able to reduce the growth medium, with efficiency increasing from S1 to S3.
In the first stage of the experiment (S1), optical density (OD600) was not measured due to the addition of a source material containing insoluble particles. In stages 2 (S2) and 3 (S3), the increase in turbidity was confirmed in all treatments, and significantly larger changes were revealed in the combinations with medium supplemented with yeast extract and tryptone (H+) and acetate addition than in those containing solely acetate or CO2/H2 (Table 3).

2.3. Microbial Community Structure Across Sediment and Media Combinations

As shown by the 16S rRNA gene analysis, Bacteria was the predominant domain in the sediment and enrichments in all stages of the experiment (Figure 2). The relative contribution of Archaea tended to increase with the successive stages of culture in media H(−)acet, H(+), and H(+)acet. Archaea species were the most prominent in S3 of the incubation in media H(+) and H(−)acet and accounted for 43 and 38%, respectively, of the total microbial community. In media H(−)CO2/H2 and H(+)CO2/H2, the Archaea contribution decreased in the subsequent stages (Figure 2). The Archaea domain was dominated by methanogens, which accounted for almost 99% of the total Archaea.
Representatives of the genus Methanothrix were the dominant methanogens in the bottom sediment; additionally, Methanocella, Methanobacterium, Methanoregula, Methanomassiliicoccus, Methanosarcina, and Methanolinea species and some methanogens that were not classified to any genus were identified (Figure 3). The culture in each medium resulted in differentiation of the methanogen composition. In the last step of the culture (S3), methanogenic communities grown on media H(−)CO2/H2, H(+), and H(+)CO2/H2 were dominated by Methanobacterium, which accounted for 100%, 57%, and 85% of the identified methanogens, respectively. In media H(−)acet and H(+)acet, the main methanogen represented Methanosarcina (97% and 41% of all the methanogens, respectively). The second important genus identified in consortia growing in H(+)acet was Methanospirillum (35%). The consortium from the medium combination H(+) additionally contained species from Methanothrix (13%), while Methanomassiliicoccus identified in H(+)acet and H(+)CO2/H2 accounted for 12% and 2% of the cultured methanogens, respectively.
The presence of non-methanogenic Archaea (Crenarchaeota (1.55%) Woesearchaeota (0.55%), Pacearchaeota (0.26%) Diapherotrites (0.015%), and Thaumarchaeota (0.015%)) was confirmed only in the native bottom sediment.
The structure of the microbial communities was clearly differentiated by available carbon sources. At the class level (Figure 4), Deltaproteobacteria (above 12%) dominated in the sediment (SZN) followed by Actinobacteria (6.2%), Betaproteobacteria (5.3%), Bacteroidia (5.3%), and Acidobacteria (5.3%), which had smaller but significant contribution in the consortium. In the first stage (S1), Deltaproteobacteria dominated in H(+)acet, H(+)CO2/H2, and H(+), Candidatus Cloacamonas in H(−)CO2/H2, and Bacteroidia in H(−)acet. The proportion of the identified microorganisms changed during the next stages of the culture. In medium H(−)CO2/H2 in the subsequent stages, an increasing proportion of Betaproteobacteria was noted (from 1.25% at the beginning to over 71% in the last stage of the culture). Increased contribution was also found in the case of Methanobacteria (from 0.4% to 26% and 18% in stages 2 and 3, respectively). The contribution of Betaproteobacteria increased with successive stages of the experiment when acetate was used (H(−)acet, from 1.6 to 25.6%) as a sole carbon source, as in the case of Methanomicrobia, (from 2.1% to 37.7%). An increase in the contribution of Methanomicrobia was also noted in the subsequent stages of H(+) media enrichment (from 3.8% to 18.2%). In the medium supplemented with multiple carbon sources H(+)CO2/H2 (tryptone, yeast extract, and CO2), there was a significant increase in the contribution of Caldisericum (to over 30%) between the first and last stages. Media based on tryptone, yeast extract, and acetate (H(+)acet) caused an increase in the contribution of Clostridia as well as Methanomicrobia in the microbial community structure (Figure 4).

2.4. Microbial Biodiversity

The microbial taxonomic diversity in the SZN sediment was high, which was reflected in the values of the biodiversity indices (Table 4). Analysis of the enrichment cultures revealed that the composition of the medium had a great impact on biodiversity. All indices used to describe the community composition as a function of the substrates added and the transition stage have shown that the H(−)CO2/H2 and H(−)acet treatments in the first cultivation stage were characterized by the greatest biodiversity and harbored nearly half of the OTUs found in SZN. A great majority of OTUs found in the initial stages of the H(−)CO2/H2 and H(−)acet treatments was lost during the subsequent transitions. The decrease was especially visible in the community that developed in H(−)CO2/H2, where only 96 OTUs were detected in the final stage of the experiment. Surprisingly, similarly low final biodiversity (97 OTUs) was found in S3 of the H(+)CO2/H2 variant. These two treatments also revealed the highest methanogenic activity which increased with subsequent transitions. The selection is confirmed by the values of the inverse Simpson index, which is strongly affected by the evenness of abundant species. Communities that developed in H(+)CO2/H2 and H(−)CO2/H2 in S3 were characterized by very low values of this index (7.6 and 2.1, respectively), which confirms that they were dominated by only a few microbial taxa (Table 4). The microbial consortia found in the other treatments were much more resistant to biodiversity loss. Interestingly, nearly 240 OTUs were still present in the communities grown on H(+) and H(+)acet in the final stage of cultivation (S3).
The beta-diversity analysis confirmed the significance of differences between the communities. The reconstructed tree describing divergence has shown clearly that the substrates added to the culturing medium were the main factors affecting the composition of the SZN-derived communities (Figure 5). CO2/H2 seems to be the primary differentiating factor. Communities that developed in H(+)CO2/H2 and H(−)CO2/H2 formed a separate branch and exhibited only a 15% similarity to those grown on solely organic compounds. The differences between the communities grown on organic substrates were smaller.

3. Discussion

3.1. Lake Sediment

Methane production in the sediment (SZN) ranged from 0.58 to over 1.67 mg CH4 L−1 d−1 depending on the incubation temperature and with a maximum at 30 °C. These values were in the range of the methane production determined for river and pond sediments investigated by Yang [26] in northern Taiwan (0.24–18.0 mg CH4 L−1 d−1) and for Lake Kinneret in Israel (0.38–7.68 mg CH4 L−1 d−1) investigated by Schwarz and co-workers [27]. Methane formation was also higher than that found in sludges from acidic bog lakes Grosse Fuchskuhle (Germany) (0.032–0.368 mg CH4 L−1 d−1) [28].
Methanogens, with the dominant Methanosarcinales (Methanothrix and Methanosarcina) able to utilize a number of different substrates, e.g., CO2, H2, and acetate (Figure 3), were responsible for the methane production in the Szczecin reservoir sediment. The contribution of representatives of Methanobacteriales, Methanomicrobiales, Methanocellales and Methanomassiliicoccales was much lower. This situation is different from that in eutrophic and mesotrophic lakes (Dianchi and Erhai, China) where the dominance of Methanomicrobiales and Methanobacteriales was confirmed, and, despite the slightly acidic reaction of the Szczecin reservoir sediment, was similar to the alkaline soda lake (Mono Lake, CA, USA) dominated by Methanosarcinales and Methanomicrobiales [29,30]. A similar methanogenic community structure was also found in several Amazonian lake sediments with variable pH (5.7 to 8.1) [31]. Deltaproteobacteria, Actinobacteria, Betaproteobacteria, Bacteroidia, and Acidobacteria (more than 5% contribution) dominated among the accompanying bacteria in the Szczecin reservoir sediment. The microbial community was slightly similar to that detected in eutrophic Chinese lakes, where Deltaproteobacteria, Betaproteobacteria, Gammaproteobacteria, Alphaproteobacteria, Anaerolineae, and Nitrospira were the dominant classes [32].

3.2. Enrichment Cultures

In each treatment, MP was observed after the lag phase. Its occurrence at the beginning of the culture-based experiment results from adaptation to changed conditions [33]. In our experiment, where the CH4 dynamics was of central importance, the lag phase duration reflects the time necessary for reduction of the culture medium to sufficiently low redox conditions. It could be expected that biodiversity loss would impair self-reduction potential of the microbial community. Our experiment shows that sediment derived from an instable subsidence reservoir (Szczecin) harbors a community that is able to maintain (or even increase) the self-reduction potential in spite of the substantial reduction of biodiversity. Furthermore, although the lag phase was the shortest in the cultures grown with H2 (which is a reducing agent itself) other substrates also supported the fast creation of conditions suitable for methanogenesis, which is important when considering further research of the biotechnological potential of the community.

3.3. Microbial Biodiversity

The biodiversity in the SZN sediment (Table 4) was similar to that in other lake sediments worldwide (analyzed using the analogous sequencing depth) [34,35] and higher than values typically calculated for soil [36,37] or freshwater [38,39]. This confirms that the transitional environments of the sediments harbor much higher microbial diversity that the adjacent environments.
Microbial diversity was largely lost in the enrichments. The highest decline was noted between SZN and S1 of the cultures for all calculated indexes (Table 4). It should be noted that the trophic structure of the microbial community in the S1 cultures was still influenced by the availability of organic matter originating from the basal sediment (dilution 2 × 10−2). The second stage S2 was transitional, while at the third stage of cultivation when the availability of the basal material was 8 × 10−6, the community was in fact dependent only on substrates delivered by the experimental medium. The decrease in biodiversity between the second and third stage of culturing was low, which suggests that the S3 community structures were determined mainly by the experimental conditions.
Subsequent transitions resulted in a decrease in unidentified sequences, which probably points to the disappearance of species that are not able to grow when deprived of the natural sediment and are unculturable to date. Furthermore, it was found that a majority of the genera that were lost during the successive stages of the experiment were aerobic (Gaiella, Gemmatimonas, Nitrospira, Conexibacter, Methylobacter, Thermoleophilum, Kofleria) or microaerobic (Sideroxydans, Magnetococcus, Anaeromyxobacter, and Sterolibacterium—threshold of 1% of the reads in SZN. The presence of aerobic microorganisms confirms the transitional character of the SZN sediment and the dualistic character of the microbial community (aerobic/anaerobic). The loss of aerobes was a consequence of the application of the culturing conditions.
The structure of the microbial communities was shaped by the available carbon sources (Figure 4). Adaptation to the culturing conditions (available carbon sources) was also clearly visible for the methanogens. Methanothrix, which is the most frequent Archaea in the SZN sediment, disappeared gradually in a majority of the treatments. In the third stage of the experiment, it was outcompeted by other methanogenic genera (Figure 3). The differentiation reflected the metabolic capabilities of Archaea. The availability of the thermodynamically beneficial energy source (H2) led to increased contribution of hydrogenotrophic Methanobacterium (in H(+), H(−)CO2/H2, and H(+)CO2/H2) while the addition of acetate to the culture medium (H(+)acet and H(−)acet promoted Methanosarcina. In this case, our experiment confirms the previously described fact that although both are acetotrophic, Methanothrix and Methanosarcina differ in the kinetics of enzymes involved in acetate assimilation. Hence, low acetate concentrations favour Methanothrix, which is outcompeted by Methanosarcina at high acetate levels [40].
Methanothrix almost completely disappeared at S3 of the experiment. Only in the H(+) culture, at S3, did it still account for 13% of the methanogens (Figure 3). If the outcompetition of Methanothrix by Methanosarcina can indeed be explained by competition for the substrate, the disappearance of Methanothrix in H(+)acet and H(−)acet requires additional explanation. The composition of the basal medium (vitamins and trace element content) was the same in all cultures. Therefore, it seems that the primary cause of the changes in methanogen contribution would be the associated bacterial community (differentiated by the available carbon source).
Methanothrix may lack some key compounds, e.g., auxotrophic vitamins or amino acids. Recently, Hubalek et al. [41] presented a comprehensive study showing that the proportion of genes encoding auxotrophy for vitamins and amino acids in the metagenome-assembled genomes of anaerobic, hydrocarbon-degrading communities is surprisingly high compared to those linked with energy conservation. This led them to the conclusion that metabolic interactions between obligate mutualistic microbial partners should be of central importance because beyond the canonical H2-producing and syntrophic bacteria - methanogen partnership, a complex (although not fully defined) interactions play an important role in determination of the metabolism of the entire community.
Our study shows that the cooperation between acetotrophic methanogens (Methanothrix) and acetogenic bacteria may be one of those relationships. It has already been proven that Methanosarcina owe its physiological flexibility to Clostridia—the most probable source of a unique (as for Archaea) enzymatic system employing acetate kinase (AckA) and phosphoacetyl transferase (Pta) [42]. Enzymes involved in methane production from acetate may not be the only ones “imported” via gene transfer. It cannot be excluded that other genes, not yet identified but enhancing survival in environmental conditions, were also incorporated by Methanosarcina. Methanothrix seems to lack such benefits. In this work, we found a relationship between the contribution of Methanothrix and acetogenic bacteria (the latter presented as KEGG-revealed expected abundances of genes responsible for the synthesis of enzymes involved in the acetate-generating Wood–Ljungdahl pathway - carbon monoxide dehydrogenase [EC 1.2.7.4]/acetyl-CoA synthase [EC 2.3.1.169]). We hypothesize that Methanothrix gains more benefits from cooperation with acetogenic bacteria than from substrate delivery and interspecies electron transfer. Examples of such a relationship have already been demonstrated in co-cultures (alanine transfer between Methanococcus maripaludis and Desulfovibrio vulgaris [43]).
A support for the deduced relationships seems to be the disappearance of Methanothrix in treatments containing acetate, where the growth of acetogenic bacteria was inhibited by excess product concentration (via a mechanism described previously by Wang and Wang [44] (Figure 6). Full confirmation of the necessity of Methanothrix – acetogen cooperation requires comprehensive studies. The use of model co-cultures subjected to transcriptional, proteomic and metabolic analyses or shotgun metagenome/metatranscriptome analyses of environmental samples would explain the exact nature of the deduced cooperation between Methanothrix and acetogenic bacteria.
Interestingly, the proportion of methanogens was not in line with the methane production rate detected in vivo (Figure 1). The treatment that exhibited the highest methane production (growing with each transfer and characterized by the shortest lag phase) was also characterized by a very low proportion of the identified methanogens (4.3% of the reads). Instead, unexpectedly, high contribution of Caldiserica was found in H(+)CO2/H2. These bacteria are one of the most intriguing elements of the SZN consortium. Caldiserica, formerly known as OP5, was first described based on environmental 16S rRNA fragments isolated from Obsidian Pool (Yellowstone) [45]. The first culturable species of the phylum, i.e., Caldisericum exile, was isolated by Mori and co-workers a decade later [46] from a hot spring in Japan. In the present experiment, Caldiserica was found in almost all of the treatments and was particularly abundant in H(+)CO2/H2. The culturing conditions in which CO2/H2, yeast extract, and tryptone were used were highly suitable for these bacteria, whose participation in the community structure grew successfully with the consecutive culture stages reaching over 30% of the sequence reads in S3. In the other variants containing H2/CO2, yeast extract, and tryptone separately, they accounted for 0.17% and 3.2% of the sequences in the final stage of the experiment. Almost no Caldiserica representatives were detected in the treatments containing acetate (even in the presence of yeast extract (H(+)acet)).
C. exile, the only known culturable representative of the phylum to date, was described as anaerobic, thermophilic, and thiosulfate-reducing bacterium. Therefore, it could be expected that the presence of Caldiserica would hinder methane production by competition for hydrogen with methanogens. In this work, we found that the presence of huge numbers of Caldiserica did not reduce methanogenesis but seemed to even stimulate it. Similar observations have been reported by Ma and co-workers [47], who investigated degradation of hexadecane to methane as a function of sulphate concentration. In that research, the most effective culture (containing 0.5 mM sulphate) contained a high proportion of Caldiserica. These authors did not emphasize the role of Caldiserica but their analysis of the whole community suggested that there is a possibility of cooperation between incomplete-oxidizing sulphate reducers and methanogens, as incomplete oxidation of organic intermediates may generate H2 through sulphate reduction. In fact, sulphur disproportionation may be carried out with protons being either a substrate or a product of the reaction, as described below (Equations (1), (2)) [48,49]:
S2O32− + 2H+ + 2e HS + HSO3
S2O32−+H2O SO42−+HS+H+
The ecology of Caldiserica is currently being discovered and described. The first reports on representatives of the phylum were associated with extreme thermophilic environments, e.g., the hot springs mentioned above [50] and hydrothermal vents [51]. Further studies provided growing evidence that representatives of the phylum can occupy other environmental niches as well. Interestingly, their high contribution (reaching as much as 60% of the total community) has been confirmed in permafrost [52,53], which denies their exclusively thermophilic character. Additionally, a negative effect of the increased temperature has been observed in this specific location [52]. Caldiserica was also identified in lake waters both in deep anoxic parts [54] and, surprisingly, the upper layers [55]. The results presented in this study pointing to the presence of these bacteria in the sediments of the shallow subsidence reservoir support the recent discoveries of Caldiserica capability to live in mesophilic conditions or even temperatures close to zero (such as those occurring in shallow lake sediments in winter) and to cope with oxidative stress. Bearing in mind the variety of environments occupied by Caldiserica (evidence of genomic diversity), it may be expected that the methodological progress in the field of environmental genomics will soon facilitate description of other species belonging to this phylum colonising various environments and eluding culturing attempts.
Other representatives of microbial “dark matter” that have recently come into the limelight and were found in the Szczecin reservoir sediments were those of Methanomassiliicoccus. This methanogenic methylotrophic Archaea was first described in human faeces [56]. Further, its relatives were found in the intestinal tracts of other organisms [57] or faeces-affected sludges (e.g., from wastewater treatment plants) [58,59]. Most research of these microorganisms focuses on their interaction with human health [60,61,62] or the unique methylotrophic but H2-dependent metabolism [63]. Methanogenic Thermoplasmata (including Methanomassilicoccales) use a reduced methanogenic pathway, in which methanol and other methylated compounds are reduced to methane in the presence of H2 [56]. This metabolic pathway has long been considered to have minor environmental importance, as it was reported to be used by only two methanogenic species. Recently, the visibility of Thermoplasmata-related sequences has been enhanced by the description of culturable species and M. luminensis genome sequence deposition in public databases [64]. The ecology of this newly described archaeal phylum is currently being recognized. Comparative phylogenetic studies performed by Paul et al. [65] have implied that Methanomassiliicoccales may be a part of the microbiome occurring in various environments. To date, these assumptions have been confirmed for extreme environments such as hot springs [66], formation waters connected with oil reservoirs [67], wetland soils [68], lake sediments [69], and deep subsurface (coal) [70]. The present study, indicating that the Szczecin reservoir sediment (SZN) is occupied by Methanomassiliicoccaceae, is in line with the aforementioned discoveries and is the second report, after Fan and Xing [69], on their presence in lake sediments. Surprisingly, the enrichment cultures in the presented experiment revealed that this group of methanogens is especially enriched in the presence of acetate. All currently published enrichment cultures and a sole M. luminyensis isolate were obtained on methanol or methylamines as a carbon source and H2 [68]. Hence, it could be expected that Methanomassiliicoccus would find the best growth conditions in medium H(+)CO2/H2, where organic substrates and hydrogen were added. Surprisingly, only 0.17% of the sequences in this treatment were affiliated to this genus, vs. nearly 2.5% in H(+)acet, which means that it was 27 times more abundant in these conditions than in the original sediment. These results contrast with previous studies of Methanomassiliicoccus-containing enrichment cultures. In experiments presented by Lv et al. [67], addition of acetate to the culture medium resulted in replacement of Methanomassiliicococcacae by Methanosaeta and Methanosarcina (both known for acetotrophic metabolism). In the present study, Methanosarcina was also dominant in all acetate-amended cultures (which is not a surprise) but in H(+)acet Methanomassiliicoccus accounted for nearly 12% of all methanogens. The difference between these two experiments may result from the different origin of the inoculates. Lv and co-workers [67] investigated communities derived from oil production waters, while our study was developed based on a community retrieved from a shallow lake sediment; therefore, they may represent distinct species. The utilization of acetate by Methanomassiliicoccaceae in lake sediments may be a result of adaptation to in situ conditions. Shallow lake littoral zones are often overgrown by aquatic macrophytes. The roots of these plants are known to exude organic acids. Acetate is thus an abundant substrate in the sediment and, provided anaerobiosis is maintained, can be used by methanogens. Our hypothesis pointing to stimulation of Methanomassiliicoccaceae by acetate exudates seems to be confirmed also by Fan and Xing [69], who reported that representatives of the genus are more abundant in littoral- macrophyte overgrown sediments than in deeper parts of the lake (dominated by algae). Also, the sequence retrieved from the SZN metagenome exhibited high similarity with 16SrDNA fragments isolated from environmental littoral samples overgrown by reed (AB896665.1) and rice (KU522088.1; GU134476.1).

4. Materials and Methods

The research material was bottom sediment (0–5 cm) of an endorheic reservoir named Szczecin (51°20’18″N 22°59’45″E, Lublin region, Poland). The reservoir developed in 1995 in the post-mining subsidence as result of underground operations carried out in the nearby “Bogdanka” coal mine (Lublin Coal Basin). Its area is about 100 ha and the maximum depth is approximately 2.5 m [71]. The reservoir basement is formed of sandy-clay soil, which is characteristic of the area [72]. The water in the reservoir has become eutrophic since nutrients are likely to get into the water from nearby fields [73]. Sediments were taken in triplicate in the summer of 2016 (Piston sampler, Eijkelkamp, Nederland). They were transferred into tightly closed sterile jars, transported to the laboratory, and kept at 5 °C for a few days until the experiment started and at −20 °C for the molecular analyses.

4.1. Physical and Chemical Analysis

The sediment moisture was determined gravimetrically by oven-drying to a constant weight at 105 °C, immediately after collecting the samples. Reaction (pH) and redox potential (Eh) were determined using a multifunctional potential meter Sension+ MM150 and a multi-sensor for pH and Eh (HACH, USA). The carbon content in dry samples was determined by means of TOC-VCSH with an SSM-5000A module autoanalyzer (Shimadzu, Japan). The total organic carbon (TOC) amount was calculated from the difference between total carbon (TC) and inorganic carbon (IC) [74]. Bioavailable forms of nitrogen (nitrite, nitrate, ammonium) and phosphorus were measured in bottom sediment extracts (sediments with deionized water and with 0.5 M NaHCO3 for nitrogen and phosphorus forms respectively) using an AA3 autoanalyzer (Braun & Luebbe, Germany) after filtering through filter paper (Munktell, grade 390, Germany) according to the method described by Banach [75].

4.2. Incubation Experiments

4.2.1. Assessment of the Methane Production Rate in Sediment

Each stage of methanogenic incubation was prepared aseptically in anoxic conditions with a nitrogen atmosphere (glove box chamber, Labconco, Kansas City, MO, USA). At the beginning, the natural capability of methane production in the sediments (SZN) was tested. Ten millilitres of bottom sediments were placed in 60-mL serum vials, closed with butyl rubber septa, and capped. Samples were incubated at temperatures of 10, 20, 30, and 40 °C for up to 90 days, always in triplicate.

4.2.2. Cultivation Conditions of the Enrichment Cultures

The media were prepared according to Horn [76] with slight modifications. Sterile (autoclaved) medium containing (milligrams per litre) (NH4)2SO4—25; CaCl2 ·x 2H2O—10; MgCl2 x 6H2O—5; NaCl—200; NH4Cl—200; and KH2PO4—200, supplemented with a vitamin solution (10 mL) and a trace element solution (10 mL), was used as a basic mineral solution (H(−)) [77]. The medium was prepared without addition of cysteine or other reducing agents, to facilitate assessment of the natural reduction potential of the investigated microbial communities. The hermeticity of the bottles (oxygen intrusion) was monitored chromatographically. No resazurin was added. The pH of the medium was adjusted to 7 with KOH. The experimental treatments were prepared with tryptone, yeast extract, acetate, and CO2/H2 (20:80 v/v) as shown in Table 5. The initial headspace pressure in N2-filled containers was close to atmospheric. In treatments with CO2/H2, nitrogen was replaced with a mixture of the given gases with overpressure of approximately 150 kPa [78].
The enrichments were carried out in 120 mL serum vials in a 1:50 proportion of the inoculating material to the fresh medium in the dark and without shaking for up to 90 days. The bottom sediments were used as an inoculum in the first stage of the enrichment (S1), whilst the liquid enrichment from the previous stages (S1 and S2, respectively) were used in the second (S2) and third (S3) stages. Incubations of the SZN and S1 enrichment cultures were carried out at temperatures between 10 and 40 °C. The number of temperature variants at subsequent stages was reduced by those characterized by the lowest methanogenic activities.
The optical density of the enrichment (OD660) was measured using a UV-1800 (Shimadzu, Japan) spectrophotometer [79].

4.2.3. Chromatographic Analyses and Calculation of the Methane Production

The rate of CH4 production and, also loss of CO2 and H2 in some research enrichment variants, was determined using a gas chromatograph (GC 3800, Varian, USA) equipped with flame ionization (FID, 200 °C) and thermal conductivity (TCD, 120 °C) detectors in series and with the use of two types of columns: a Poraplot Q 0.53 mm ID (25 m) and a Molecular Sieve 5A 0.53 mm ID (30 m) connected together, were used with helium as the carrier gas [74,80]. Methane production rate (MP) was determined on the basis of the linear increase in the methane concentration in time and expressed as a mg of produced CH4 per litre of bottom sediments or enrichment medium per day (mg CH4 L−1 d−1).

4.3. DNA Extraction and NGS Procedure

Microbial genomic DNA from the bottom sediment and cultures was extracted using a PowerLyzer PowerSoil DNA Isolation Kit (Quiagen, Hilden, Germany) according to the manufacturer’s instructions. 0.25 g of fresh material from the bottom sediments and cell pellets harvested from 4 mL of the enrichments (after centrifugation, 13000× RPM, 5 min) in the last incubation step (S3) were used for isolation. The presence of DNA was confirmed by electrophoresis in a 1% agarose gel with 1×TBE buffer and a SimplySafe™ (EURX) stain for detection of nucleic acid. The V3–V4 region of 16S rRNA gene amplicons were sequenced with the MiSeq Illumina technology (Genomed Inc, Poland).

4.4. Bioinformatic Analysis

Amplicon sequence variants (ASVs) were resolved with DADA2 version 1.8 package [81] in R version 3.5.1 [82]. Based on the sequence quality plots, forward and reverse reads were trimmed, respectively, to 250 and 240 bp, and primer sequences were removed from all reads. The following filtering parameters were used: maxN = 0, maxEE for the forwards reads = 3 and for the reverse reads = 5, truncQ = 2. Other parameters were set to default. The error rates were estimated by learnErrors using one million reads. Sequences were dereplicated using derepFastq with default parameters and exact sequence variants were resolved using dada. Next, removeBimeraDenovo was used to remove chimeric sequences.
Taxonomy was assigned against the latest version of the RDP database (11 version) using a Naïve Bayesian Classifier [83] with the minboot parameter set to 80. The resulting taxonomy and read-count tables constructed in DADA2 were appropriately converted and imported into the phyloseq (1.22.3) package [84]. Sequences identified as chloroplast and mitochondria were removed.
16S rRNA amplicon sequencing data generated in this study were deposited in the NCBI Sequence Read Archive (SRA) under the BioProject number PRJNA514232.

4.5. Metabolic Pathways Prediction Using PICRUSt Software

Using PICRUSt (Phylogenetic Investigation of Communities by Reconstruction of Unobserved States, software) [85], the functional genes were reconstructed based on bacterial 16S rRNA gene sequences. The input ASVs from the dada2 workflow were clustered into OTUs (operational taxonomic units) using the cluster-features-closed-reference algorithm implemented in qiime2 software [86] with the 97% similarity threshold and using the GreenGenes (13.8) database as a reference. The number of 16S rRNA gene copies was normalized using the normalize_by_copy_number.py algorithm implemented in PICRUSt software. Metagenomes were predicted using the predict_metagenomes.py algorithm, and the functional genes were annotated in the KEGG database using the categorize_by_function.py algorithm.

5. Conclusions

The results obtained in this study show that our current understanding of microbial processes leading to degradation of organic matter to methane is still far from complete. We have shown that in addition to the already identified relations in the microbial world leading to methane production, there is an enormous scope of biotechnological potential hidden in the microbial “dark matter”. In light of the great diversity and hardly graspable dynamics of microbial communities in the natural environment, our knowledge is still insufficient, although it is expanding year on year. In the present work, we have shown that the enrichment cultures diversified by the available carbon substrate may serve as a useful tool for the study of microbial ecology and physiology, giving a chance to highlight the activity of microorganisms (e.g., Caldisericum, Methanomassiliicoccus, Methanothrix) that escape recognition using other methods (e.g., based on pure culture analysis). Biotechnology can gain from studies of interspecies relations. Our research suggest that Caldisericum ssp., which are likely to cooperate with methanogens, may have an unexploited potential for the biogas production industry.
Furthermore, their high methane production rate and self-reduction capacity imply that further studies should be undertaken to test the potential of the cultures isolated from lake sediments for biogasification of various organic material. Furthermore, we suggested that more efforts should be made to investigate the role of Caldiserica in methane formation processes due to its putative positive impact that may have biotechnological application.

Author Contributions

Conceptualization, A.S.-N., A.P.; Investigation, A.S.-N., A.P., J.G., A.K. (Adam Kubaczyński), A.B., A.G. (Andrzej Górski); Formal analysis, A.S.-N., A.P., J.G.; Data curation, J.G., W.G., A.K. (Agnieszka Kuźniar); Methodology, A.S.-N., A.P., J.G.; Visualization, A.S.-N., A.P., J.G.; Writing—original draft preparation, A.S.-N., A.P.; Resources, A.G. (Anna Gałązka); Supervision, Z.S.; Funding acquisition, Z.S. Article revision, A.S.-N., A.P.

Funding

This research was funded by the National Science Centre (Poland), grant number 2015/17/B/NZ9/01662.

Acknowledgments

We thank Magdalena Hunicz for her help in the laboratory work.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Abbreviations

KEGGKyoto Encyclopedia of Genes and Genomes
PICRUStPhylogenetic Investigation of Communities by Reconstruction of Unobserved States, software

References

  1. Boughner, L.A.; Singh, P. Microbial Ecology: Where are we now? Postdoc. J. 2016, 4, 3–17. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Eckert, I.M.K.; Littlefair, J.E.; Zhang, G.K.; Chain, F.J.J.; Crease, T.J.; Cristescu, M.E. Bioinformatics for biomonitoring: Species detection and diversity estimates across Next-Generation Sequencing Platforms. Adv. Ecol. Res. 2018, 59, 1–32. [Google Scholar] [CrossRef]
  3. Costessi, A.; van den Bogert, B.; May, A.; Ver Loren van Themaat, E.; Roubos, J.A.; Kolkman, M.A.B.; Butler, D.; Pirovano, W. Novel sequencing technologies to support industrial biotechnology. FEMS Microbiol. Lett. 2018, 365, fny103. [Google Scholar] [CrossRef] [PubMed]
  4. Detman, A.; Mielecki, D.; Pleśniak, Ł.; Bucha, M.; Janiga, M.; Matyasik, I.; Chojnacka, A.; Jędrysek, M.-O.; Błaszczyk, B.K.; Sikora, A. Methane-yielding microbial communities processing lactate-rich substrates: A piece of the anaerobic digestion puzzle. Biotechnol. Biofuels 2018, 11, 116. [Google Scholar] [CrossRef] [PubMed]
  5. Wang, S.; Hou, X.; Su, H. Exploration of the relationship between biogas production and microbial community under high salinity conditions. Sci. Rep. 2017, 7, 1149. [Google Scholar] [CrossRef] [PubMed]
  6. Wojcieszak, M.; Pyzik, A.; Poszytek, K.; Krawczyk, P.S.; Sobczak, A.; Lipinski, L.; Roubinek, O.; Palige, J.; Sklodowska, A.; Drewniak, L. Adaptation of methanogenic inocula to anaerobic digestion of maize silage. Front. Microbiol. 2017, 8, 1881. [Google Scholar] [CrossRef] [PubMed]
  7. Fuertez, J.; Boakye, R.; McLennan, J.; Adamsc, D.J.; Sparks, T.D.; Gottschalke, A. Developing methanogenic microbial consortia from diverse coal sources and environments. J. Nat. Gas Sci. Eng. 2017, 46, 637–650. [Google Scholar] [CrossRef]
  8. Green, M.S.; Flanegan, K.C.; Gilcrease, P.C. Characterization of a methanogenic consortium enriched from a coalbed methane well in the Powder River Basin, U.S.A. Int. J. Coal Geol. 2008, 76, 34–45. [Google Scholar] [CrossRef]
  9. Plugge, C.M. Biogas. Microb. Biotechnol. 2017, 10, 1128–1130. [Google Scholar] [CrossRef] [Green Version]
  10. Rapheal, S.V.; Swaminathan, K.R.; Lalitha, K. Metabolic characteristics of an aerobe isolated from a methylotrophic methanogenic enrichment culture. J. Biosci. 2003, 28, 235–242. [Google Scholar] [CrossRef]
  11. Enzmann, F.; Mayer, F.; Rother, M.; Holtmann, D. Methanogens: Biochemical background and biotechnological applications. AMB Expr. 2018, 8, 1. [Google Scholar] [CrossRef] [PubMed]
  12. Evans, P.N.; Parks, D.H.; Chadwick, G.L.; Robbins, S.J.; Orphan, W.J.; Golding, S.D.; Tyson, G.W. Methane metabolism in the archaeal phylum Bathyarchaeota revealed by genome-centric metagenomics. Science 2015, 350, 434–438. [Google Scholar] [CrossRef] [PubMed]
  13. Vanwonterghem, I.; Evans, P.N.; Parks, D.H.; Jensen, P.D.; Woodcroft, B.J.; Hugenholtz, P.; Tyson, G.W. Methylotrophic methanogenesis discovered in the archaeal phylum Verstraetearchaeota. Nat. Microbiol. 2016, 1, 16170. [Google Scholar] [CrossRef] [PubMed]
  14. Purwantini, E.; Torto-Alalibo, T.; Lomax, J.; Setubal, J.C.; Tyler, B.M.; Mukhopadhyay, B. Genetic resources for methane production from biomass described with the Gene Ontology. Front. Microbiol. 2014, 5, 634. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Conrad, R.; Klose, M.; Claus, P.; Enrich-Prast, A. Methanogenic pathway, C-13 isotope fractionation, and archaeal community composition in the sediment of two clear-water lakes of Amazonia. Limnol. Oceanogr. 2010, 55, 689–702. [Google Scholar] [CrossRef]
  16. Lovley, D.R.; Klug, M.J. Model for the distribution of sulfate reduction and methanogenesis in freshwater sediments. Geochim. Cosmochim. Acta 1986, 50, 11–18. [Google Scholar] [CrossRef] [Green Version]
  17. Esteves, A.I.S.; Hardoim, C.C.P.; Xavier, J.R.; Gonçalves, J.M.; Costa, R. Molecular richness and biotechnological potential of bacteria cultured from Irciniidae sponges in the north-east Atlantic. FEMS Microbiol. Ecol. 2013, 85, 519–536. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Zurita, G.; Pe’er, G.; Bellocq, M.I.; Hansbauer, M.M. Edge effects and their influence on habitat suitability calculations: A continuous approach applied to birds of the Atlantic forest. J Appl. Ecol. 2012, 49, 503–512. [Google Scholar] [CrossRef]
  19. Kark, S. Effects of ecotones on biodiversity. In Encyclopedia of Biodiversity; Levin, S.A., Ed.; Elsevier Inc.: Amsterdam, The Netherlands, 2007; pp. 1–10. [Google Scholar]
  20. Dąbrowska-Prot, E.; Wasiłowska, A. Ecological importance of meadow patches in protected forest area: Floristic diversity and the dynamics of insect communities. Pol. J. Ecol. 2010, 58, 741–758. [Google Scholar]
  21. Lachavanne, J.-B.; Juge, R. Biodiversity in Land–Inland Water Ecotones; Taylor & Francis: Abingdon, UK, 1997; p. 326. [Google Scholar]
  22. Shen, C.; Shi, Y.; Ni, Y.; Deng, Y.; Van Nostrand, J.D.; He, Z.; Zhou, J.; Chu, H. Dramatic increases of soil microbial functional gene diversity at the treeline ecotone of Changbai Mountain. Front. Microbiol. 2016, 7, 1184. [Google Scholar] [CrossRef]
  23. Mieczan, T.; Adamczuk, M.; Tarkowska-Kukuryk, M.; Nawrot, D. Effect of water chemistry on zooplanktonic and microbial communities across freshwater ecotones in different macrophyte-dominated shallow lakes. J. Limnol. 2016, 75, 262–274. [Google Scholar] [CrossRef]
  24. Jiang, H.; Dong, H.; Zhang, G.; Yu, B.; Chapman, L.R.; Fields, M.W. Microbial diversity in water and sediment of Lake Chaka, an athalassohaline lake in northwestern China. Appl. Environ. Microbiol. 2006, 72, 3832–3845. [Google Scholar] [CrossRef] [PubMed]
  25. Laskar, F.; Das Purkayastha, S.; Sen, A.; Bhattacharya, M.K.; Misra, B.B. Diversity of methanogenic Archaea in freshwater sediments of lacustrine ecosystems. J. Basic Microbiol. 2018, 58, 101–119. [Google Scholar] [CrossRef] [PubMed]
  26. Yang, S.S. Methane production in river and lake sediments in Taiwan. Environ. Geochem. Health 1998, 20, 245–249. [Google Scholar] [CrossRef]
  27. Schwarz, J.I.K.; Eckert, W.; Conrad, R. Response of the methanogenic microbial community of a profundal lake sediment (Lake Kinneret, Israel) to algal deposition. Limnol. Oceanogr. 2008, 53, 113–121. [Google Scholar] [CrossRef] [Green Version]
  28. Chan, O.C.; Wolf, M.; Hepperle, D.; Casper, P. Methanogenic archaeal community in the sediment of an artificially partitioned acidic bog lake. FEMS Microbiol. Ecol. 2002, 42, 119–129. [Google Scholar] [CrossRef] [PubMed]
  29. Yang, Y.; Dai, Y.; Li, N.; Li, B.; Xie, S.; Liu, Y. Temporal and spatial dynamics of sediment anaerobic ammonium oxidation (Anammox) bacteria in freshwater lakes. Microb. Ecol. 2017, 73, 285–295. [Google Scholar] [CrossRef]
  30. Rojas, P.; Rodríguez, N.; de la Fuente, V.; Sánchez-Mata, D.; Amils, R.; Sanz, J.L. Microbial diversity associated with the anaerobic sediments of a soda lake (Mono Lake, California, USA). Can. J. Microbiol. 2018, 64, 385–392. [Google Scholar] [CrossRef] [Green Version]
  31. Ji, Y.; Angel, R.; Klose, M.; Claus, P.; Marotta, H.; Pinho, L.; Enrich-Prast, A.; Conrad, R. Structure and function of methanogenic microbial communities in sediments of Amazonian lakes with different water types. Environ. Microbiol. 2016, 18, 5082–5100. [Google Scholar] [CrossRef]
  32. Huang, W.; Chen, X.; Jiang, X.; Zheng, B. Characterization of sediment bacterial communities in plain lakes with different trophic statuses. MicrobiologyOpen 2017, 6, e00503. [Google Scholar] [CrossRef]
  33. Xing, L.; Yang, S.; Yin, Q.; Xie, S.; Strong, P.J.; Wu, G. Effects of carbon source on methanogenic activities and pathways incorporating metagenomic analysis of microbial community. Bioresour. Technol. 2017, 244, 982–988. [Google Scholar] [CrossRef] [PubMed]
  34. Sinclair, L.; Osman, O.A.; Bertilsson, S.; Eiler, A. Microbial community composition and diversity via 16S rRNA gene amplicons: Evaluating the illumina platform. PLoS ONE 2015, 10, e0116955. [Google Scholar] [CrossRef] [PubMed]
  35. Wan, Y.; Ruan, X.; Zhang, Y.; Li, R. Illumina sequencing-based analysis of sediment bacteria community in different trophic status freshwater lakes. MicrobiologyOpen 2017, 6, e00450. [Google Scholar] [CrossRef] [PubMed]
  36. Lan, G.; Li, Y.; Jatoi, M.T.; Tan, Z.; Wu, Z.; Xie, G. Change in soil microbial community compositions and diversity following the conversion of tropical forest to rubber plantations in Xishuangbanan, Southwest China. Trop. Conserv. Sci. 2017, 10, 1–14. [Google Scholar] [CrossRef]
  37. Wu, X.; Xu, H.; Liu, G.; Ma, X.; Mu, C.; Zhao, L. Bacterial communities in the upper soil layers in the permafrost regions on the Qinghai-Tibetan plateau. Appl. Soil Ecol. 2017, 120, 81–88. [Google Scholar] [CrossRef]
  38. Chen, Y.; Dai, Y.; Wang, Y.; Wu, Z.; Xie, S.; Liu, Y. Distribution of bacterial communities across plateau freshwater lake and upslope soils. J. Environ. Sci. 2016, 43, 61–69. [Google Scholar] [CrossRef] [PubMed]
  39. Parulekar, N.N.; Kolekar, P.; Jenkins, A.; Kleiven, S.; Utkilen, H.; Johansen, A.; Sawant, S.; Kulkarni-Kale, U.; Kale, M.; Sæbø, M. Characterization of bacterial community associated with phytoplankton bloom in a eutrophic lake in South Norway using 16S rRNA gene amplicon sequence analysis. PLoS ONE 2017, 12, e0173408. [Google Scholar] [CrossRef] [PubMed]
  40. Jetten, M.S.M.; Stams, A.J.M.; Zehnder, A.J.B. Methanogenesis from acetate: A comparison of the acetate metabolism in Methanothrix soehngenii and Methanosarcina spp. FEMS Microbiol. Lett. 1992, 88, 181–197. [Google Scholar] [CrossRef]
  41. Hubalek, V.; Buck, M.; Tan, B.; Foght, J.; Wendeberg, A.; Berry, D.; Bertilsson, S.; Eiler, A. Vitamin and Amino Acid Auxotrophy in Anaerobic Consortia Operating under Methanogenic Conditions. Msystems 2017, 2, e00038-17. [Google Scholar] [CrossRef] [Green Version]
  42. Fournier, G.P.; Gogarten, J.P. Evolution of acetoclastic methanogenesis in Methanosarcina via horizontal gene transfer from cellulolytic Clostridia. J. Bacteriol. 2008, 190, 1124–1127. [Google Scholar] [CrossRef]
  43. Walker, C.B.; Redding-Johanson, A.M.; Baidoo, E.E.; Rajeev, L.; Hendrickson, E.L.; Joachimiak, M.P.; Stolyar, S.; Arkin, A.P.; Leigh, J.A.; Zhou, J.A.; et al. Functional responses of methanogenic archaea to syntrophic growth. ISME J. 2012, 6, 2045–2055. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Wang, G.; Wang, D.I. Elucidation of growth inhibition and acetic acid production by Clostridium thermoaceticum. Appl. Environ. Microbiol. 1984, 47, 294–298. [Google Scholar] [PubMed]
  45. Hugenholtz, P.; Pitulle, C.; Hershberger, K.L.; Pace, N.R. Novel division level bacterial diversity in a Yellowstone hot spring. J. Bacteriol. 1998, 180, 366–376. [Google Scholar] [PubMed]
  46. Mori, K.; Yamaguchi, K.; Sakiyama, Y.; Urabe, T.; Suzuki, K. Caldisericum exile gen. nov., sp. nov., an anaerobic, thermophilic, filamentous bacterium of a novel bacterial phylum, Caldiserica phyl. nov., originally called the candidate phylum OP5, and description of Caldisericaceae fam. nov., Caldisericales ord. no. Int. J. Syst. Evol. Microbiol. 2009, 59, 2894–2898. [Google Scholar] [CrossRef] [PubMed]
  47. Ma, T.-T.; Liu, L.-Y.; Rui, J.-P.; Yuan, Q.; Feng, D.-S.; Zhou, Z.; Dai, L.-R.; Zena, W.-Q.; Zhang, H.; Cheng, L. Coexistence and competition of sulfate-reducing and methanogenic populations in an anaerobic hexadecane-degrading culture. Biotechnol. Biofuels. 2017, 10, 207. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Finster, K. Microbiological disproportionation of inorganic sulfur compounds. J. Sulfur. Chem. 2008, 29, 281–292. [Google Scholar] [CrossRef]
  49. Loka Bharathi, P.A. Sulfur cycle. In Encyclopedia of Ecology; Jørgensen, S.E., Fath, B.D., Eds.; Academic Press: Oxford, UK, 2008; pp. 3424–3431. [Google Scholar]
  50. Mori, K.; Sunamura, M.; Yanagawa, K.; Ishibashi, J.; Miyoshi, Y.; Iino, T.; Suzuki, K.; Urabe, T. First Cultivation and ecological investigation of a bacterium affiliated with the candidate phylum OP5 from hot springs. Appl. Environ. Microbiol. 2008, 74, 6223–6229. [Google Scholar] [CrossRef] [PubMed]
  51. Stott, M.B.; Saito, J.A.; Crowe, M.A.; Dunfield, P.F.; Hou, S.; Nakasone, E.; Daughney, C.J.; Smirnova, A.V.; Mountain, B.W.; Takai, K.; et al. Culture-independent characterization of a novel microbial community at a hydrothermal vent at Brothers volcano, Kermadec arc, New Zealand. J. Geophys. Res. Solid Earth 2008, 113, B08S06. [Google Scholar] [CrossRef]
  52. Monteux, S.; Weedon, J.T.; Blume-Werry, G.; Gavazov, K.; Jassey, E.J.; Johansson, M.; Keuper, F.; Olid, C.; Dorrepaal, E. Long-term in situ permafrost thaw effects on bacterial communities and potential aerobic respiration. ISME J. 2018, 12, 2129–2141. [Google Scholar] [CrossRef] [Green Version]
  53. Martinez, M.A.; Woodcroft, B.J.; Ignacio Espinoza, J.C.; Zayed, A.A.; Singleton, C.M.; Boyd, J.A.; Li, Y.F.; Purvine, S.; Maughan, H.; Hodgkins, S.B.; et al. Discovery and ecogenomic context of a global Caldiserica-related phylum active in thawing permafrost, Candidatus Cryosericota phylum nov., Ca. Cryosericia class nov., Ca. Cryosericales ord. nov., Ca. Cryosericaceae fam. nov., comprising the four species Cryosericum septentrionale gen. nov. sp. nov., Ca. C. hinesii sp. nov., Ca. C. odellii sp. nov., Ca. C. terrychapinii sp. nov. Syst. Appl. Microbiol. 2019, 42, 54–56. [Google Scholar] [CrossRef]
  54. Tiodjio, R.E.; Sakatoku, A.; Nakamura, A.; Tanaka, D.; Fantong, W.Y.; Tchakam, K.B.; Tanyileke, G.; Ohba, T.; Hell, V.J.; Kusakaba, M.; et al. Bacterial and archaeal communities in Lake Nyos (Cameroon, Central Africa). Sci. Rep. 2014, 4, 6151. [Google Scholar] [CrossRef] [PubMed]
  55. Kurilkina, M.I.; Zakharova, Y.R.; Galachyants, Y.P.; Petrova, D.P.; Bukin, Y.S.; Domysheva, V.M.; Blinov, W.; Likhoshway, Y.V. Bacterial community composition in the water column of the deepest freshwater Lake Baikal as determined by next-generation sequencing. FEMS Microbiol. Ecol. 2016, 92, fiw094. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Dridi, B.; Fardeau, M.-L.; Ollivier, B.; Raoult, D.; Drancourt, M. Methanomassiliicoccus luminyensis gen. nov., sp. nov., a methanogenic archaeon isolated from human faeces. Int. J. Syst. Evol. Microbiol. 2012, 62, 1902–1907. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Borrel, G.; Parisot, N.; Harris, H.M.B.; Peyretaillade, E.; Gaci, N.; Tottey, W.; Bardot, O.; Raymann, K.; Gribaldo, S.; Peyret, P.; et al. Comparative genomics highlights the unique biology of Methanomassiliicoccales, a Thermoplasmatales-related seventh order of methanogenic archaea that encodes pyrrolysine. BMC Genomics 2014, 15, 679. [Google Scholar] [CrossRef] [PubMed]
  58. Chojnacka, A.; Szczęsny, P.; Błaszczyk, M.K.; Zielenkiewicz, U.; Detman, A.; Salamon, A.; Sikora, A. Noteworthy facts about a methane-producing microbial community processing acidic effluent from sugar beet molasses fermentation. PLoS ONE 2015, 10, e0128008. [Google Scholar] [CrossRef] [PubMed]
  59. Liu, C.; Li, H.; Zhang, Y.; Si, D.; Chen, Q. Evolution of microbial community along with increasing solid concentration during high-solids anaerobic digestion of sewage sludge. Bioresour. Technol. 2016, 216, 87–94. [Google Scholar] [CrossRef] [Green Version]
  60. Brugère, J.-F.; Borrel, G.; Gaci, N.; Tottey, W.; O’Toole, P.W.; Malpuech-Brugère, C. Archaebiotics: Proposed therapeutic use of archaea to prevent trimethylaminuria and cardiovascular disease. Gut Microbes 2014, 5, 5–10. [Google Scholar] [CrossRef]
  61. Dridi, B.; Henry, M.; Richet, H.; Raoult, H.; Drancourt, M. Age-related prevalence of Methanomassiliicoccus luminyensis in the human gut microbiome. Apmis 2012, 120, 773–777. [Google Scholar] [CrossRef]
  62. Bang, C.; Vierbuchen, T.; Gutsmann, T.; Heine, H.; Schmitz, R.A. Immunogenic properties of the human gut-associated archaeon Methanomassiliicoccus luminyensis and its susceptibility to antimicrobial peptides. PLoS ONE 2017, 12, e0185919. [Google Scholar] [CrossRef]
  63. Kröninger, L.; Gottschling, J.; Deppenmeier, U. Growth characteristics of Methanomassiliicoccus luminyensis and expression of methyltransferase encoding genes. Archaea 2017, 2017, 2756573. [Google Scholar] [CrossRef]
  64. Gorlas, A.; Robert, C.; Gimenez, G.; Drancourt, M.; Raoult, D. Complete genome sequence of Methanomassiliicoccus luminyensis, the largest genome of a human-associated Archaea species. J. Bacteriol. 2012, 194, 4745. [Google Scholar] [CrossRef] [PubMed]
  65. Paul, K.; Nonoh, J.O.; Mikulski, L.; Brune, A. ‘Methanoplasmatales,’ Thermoplasmatales-related archaea in termite guts and other environments, are the seventh order of methanogens. Appl. Environ. Microbiol. 2012, 78, 8245–8253. [Google Scholar] [CrossRef] [PubMed]
  66. Bowen De León, K.; Gerlach, R.; Peyton, B.; Fields, M. Archaeal and bacterial communities in three alkaline hot springs in Heart Lake Geyser Basin, Yellowstone National Park. Front. Microbiol. 2013, 4, 330. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Lv, L.; Zhou, L.; Wang, L.-Y.; Liu, J.-F.; Gu, J.-D.; Mu, B.-Z.; Yang, S.-Z. Selective inhibition of methanogenesis by sulfate in enrichment culture with production water from low-temperature oil reservoir. Int. Biodeterior. Biodegrad. 2016, 108, 133–141. [Google Scholar] [CrossRef]
  68. Söllinger, A.; Schwab, C.; Weinmaier, T.; Loy, A.; Tveit, A.T.; Schleper, C.; Urich, T. Phylogenetic and genomic analysis of Methanomassiliicoccales in wetlands and animal intestinal tracts reveals clade-specific habitat preferences. FEMS Microbiol. Ecol. 2016, 92, fiv149. [Google Scholar] [CrossRef] [PubMed]
  69. Fan, X.; Xing, P. Differences in the composition of archaeal communities in sediments from sontrasting zones of Lake Taihu. Front. Microbiol. 2016, 7, 1510. [Google Scholar] [CrossRef] [PubMed]
  70. Su, X.; Zhao, W.; Xia, D. The diversity of hydrogen-producing bacteria and methanogens within an in situ coal seam. Biotechnol. Biofuels 2018, 11, 245. [Google Scholar] [CrossRef] [PubMed]
  71. LWB. Integrated Report of the Lubelski Węgiel BOGDANKA Capital Group for the Year 2016; Lubelski Węgiel Bogdanka SA: Lublin Voivodeship, Poland, 2017; pp. 122–134. [Google Scholar]
  72. Ciosmak, M. The anthropogenic inundated area ”Szczecin” in the zone of active underground exploitation of hard coal in the Lublin coal basin. J. Ecol. Eng. 2016, 17, 70–80. [Google Scholar] [CrossRef]
  73. Sawicki, B. Geological and mining conditions of sublevel operation. In Proceedings of the 12th International Scientific and Technical Conference—Mining Natural Hazards, Katowice, Poland, 21–24 November 2005; The Central Mining Institute: Katowice, Poland, 2005; pp. 75–79. (In Polish). [Google Scholar]
  74. Szafranek-Nakonieczna, A.; Stępniewska, Z. The influence of the aeration status (ODR, Eh) of peat soils on their ability to produce methane. Wetl. Ecol. Manag. 2015, 23, 665–676. [Google Scholar] [CrossRef] [Green Version]
  75. Banach, A.M.; Banach, K.; Visser, E.J.W.; Stępniewska, Z.; Smits, A.J.M.; Roelofs, J.G.M.; Lamers, L.P.M. Effects of summer flooding on floodplain biogeochemistry in Poland; implications for increased flooding frequency. Biogeochemistry 2009, 92, 247–262. [Google Scholar] [CrossRef] [Green Version]
  76. Horn, M.A.; Matthies, C.; Küsel, K.; Schramm, A.; Drake, H.L. Hydrogenotrophic methanogenesis by moderately acid-tolerant methanogens of a methane-emitting acidic peat. Appl. Environ. Microbiol. 2003, 69, 74–83. [Google Scholar] [CrossRef] [PubMed]
  77. Balch, W.E.; Fox, G.E.; Magrum, L.J.; Woese, C.R.; Wolfe, R.S. Methanogens: Reevaluation of a unique biological group. Microbiol. Rev. 1979, 43, 260–296. [Google Scholar] [PubMed]
  78. Wolfe, R.S. Techniques for cultivating methanogens. Methods Enzymol. 2011, 494, 1–22. [Google Scholar] [CrossRef] [PubMed]
  79. Stępniewska, Z.; Goraj, W.; Kuźniar, A.; Łopacka, N.; Małysza, M. Enrichment culture and identification of endophytic methanotrophs isolated from peatland plants. Folia Microbiol. (Praha) 2017, 62, 381–391. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  80. Szafranek-Nakonieczna, A.; Bennicelli, R.P. Ability of peat soil to oxidize methane and affect of temperature and layer deposition. Pol. J. Environ. Stud. 2010, 19, 805–810. [Google Scholar]
  81. Callahan, B.J.; McMurdie, P.J.; Rosen, M.J.; Han, A.W.; Johnson, A.J.A.; Holmes, S.P. DADA2: High-resolution sample inference from Illumina amplicon data. Nat. Methods 2016, 13, 581–583. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  82. Development Core Team. A Language and Environment for Statistical Computing In R Foundation for Statistical Computing; R Development Core Team: Vienna, Austria, 2018; Available online: http: //www.R-project.org (accessed on 10 June 2018).
  83. Wang, Q.; Garrity, G.M.; Tiedje, J.M.; Cole, J.R. Naïve Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl. Environ. Microbiol. 2007, 73, 5261–5267. [Google Scholar] [CrossRef]
  84. McMurdie, P.J.; Holmes, S. phyloseq: An R Package for Reproducible Interactive Analysis and Graphics of Microbiome Census Data. PLoS ONE 2013, 8, e61217. [Google Scholar] [CrossRef]
  85. Langille, M.G.I.; Zaneveld, J.; Caporaso, J.G.; McDonald, D.; Knights, D.; Reyes, J.A.; Clemente, J.C.; Burkepile, C.E.; Thurber, R.L.V.; Knight, R.; et al. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat. Biotechnol. 2013, 31, 814–821. [Google Scholar] [CrossRef]
  86. Bolyen, E.; Rideout, J.R.; Dillon, M.R.; Bokulich, N.A.; Abnet, C.; Al-Ghalith, G.A.; Alexander, H.; Alm, E.J.; Arumugam, M.; Asnicar, F.; et al. QIIME 2: Reproducible, interactive, scalable, and extensible microbiome data science. PeerJ Preprints 2018. [Google Scholar] [CrossRef]
Figure 1. Methane production rate in particular media and stages of enrichment. Mean values with standard deviation (SD) are presented.
Figure 1. Methane production rate in particular media and stages of enrichment. Mean values with standard deviation (SD) are presented.
Ijms 20 04415 g001
Figure 2. Archaea and Bacteria contribution in the particular stages of the experiment.
Figure 2. Archaea and Bacteria contribution in the particular stages of the experiment.
Ijms 20 04415 g002
Figure 3. Changes in the methanogen community during incubation in different media.
Figure 3. Changes in the methanogen community during incubation in different media.
Ijms 20 04415 g003
Figure 4. Diversity in the microbial communities at the class level (microbial classes represented by >2% of total sequences, classes making up less than 2% were classified as others) in response to the different media composition and in the particular stage of enrichment identified by 16S rRNA gene sequences.
Figure 4. Diversity in the microbial communities at the class level (microbial classes represented by >2% of total sequences, classes making up less than 2% were classified as others) in response to the different media composition and in the particular stage of enrichment identified by 16S rRNA gene sequences.
Ijms 20 04415 g004
Figure 5. Beta-diversity of the bacterial communities in the experimental treatments. The similarity presented on the scale was calculated based on the Bray-Curtis distances, and the tree was generated by an UPGMA algorithm. The bootstrap number was set to 999.
Figure 5. Beta-diversity of the bacterial communities in the experimental treatments. The similarity presented on the scale was calculated based on the Bray-Curtis distances, and the tree was generated by an UPGMA algorithm. The bootstrap number was set to 999.
Ijms 20 04415 g005
Figure 6. Relative contribution of Methanothrix and acetogenic bacteria (the latter presented as KEGG-revealed expected abundances of genes responsible for the synthesis of enzymes involved in the acetate-generating Wood–Ljungdahl pathway - carbon monoxide dehydrogenase [EC 1.2.7.4]/acetyl-CoA synthase [EC 2.3.1.169].
Figure 6. Relative contribution of Methanothrix and acetogenic bacteria (the latter presented as KEGG-revealed expected abundances of genes responsible for the synthesis of enzymes involved in the acetate-generating Wood–Ljungdahl pathway - carbon monoxide dehydrogenase [EC 1.2.7.4]/acetyl-CoA synthase [EC 2.3.1.169].
Ijms 20 04415 g006
Table 1. Characteristics of investigated sediments (± SD).
Table 1. Characteristics of investigated sediments (± SD).
N-NO2N-NO3N-NH4P-PO4TCICTOCMoisturepHEh
mg kg d.w.−1g kg d.w.−1%mV
0.11 ± 0.010.30 ± 0.0428.21 ± 4.8810.05 ± 0.2548.0 ± 0.80.0048.0 ± 0.858.71 ± 2.285.46 ± 0.02−265.63 ± 2.89
Table 2. Lag phase duration in particular media and stages of enrichment.
Table 2. Lag phase duration in particular media and stages of enrichment.
Stages of Bacterial CultivationTemperature (°C)Lag Phase (d)
H(−)CO2/H2 H(−)acetH(+)H(+)CO2/H2H(+)acet
S1101431141431
201021101010
30110111
4017111
S22066616
3017111
4016111
S320-4---
301-111
Table 3. Changes in physicochemical parameters between the start and the end of the particular stages of incubation (± SD).
Table 3. Changes in physicochemical parameters between the start and the end of the particular stages of incubation (± SD).
Stages of Bacterial Cultivation/ParametersMedium Variants
H(−)CO2/H2 H(−)acetH(+)H(+)CO2/H2H(+)acet
Δ (start to end)
S1pH−0.5 ± 0.01−0.24 ± 0.03−0.13 ± 0.02−0.21 ± 0.00−0.14 ± 0.01
Eh (mV)−383.20 ± 0.32−264.44 ± 0.06−337.27 ± 0.40−368.20 ± 0.10−303.87 ± 0.25
ODnanananana
S2pH−0.24 ± 0.011.07 ± 0.000.37 ± 0.00−0.36 ± 0.010.79 ± 0.02
Eh (mV)−330.60 ± 0.10−463.17 ± 0.30−347.77 ± 0.49−358.20 ± 0.10−467.53 ± 0.84
OD0.09 ± 0.020.08 ± 0.020.15 ± 0.010.09 ± 0.020.14 ± 0.01
S3pH−1.18 ± 0.021.49 ± 0.000.16 ± 0.04−1.56 ± 0.040.76 ± 0.01
Eh (mV)−405.03 ± 0.49−466.33 ± 0.15−400.50 ± 1.42−465.80 ± 0.21−433.16 ± 1.05
OD0.09 ± 0.010.07 ± 0.020.16 ± 0.010.16 ± 0.020.16 ± 0.01
Table 4. Biodiversity indices of the experimental treatments.
Table 4. Biodiversity indices of the experimental treatments.
TreatmentStageObservedChao1ACEShannonSimpsonInvSimpsonFisher
SZN-25362585.962565.617.061.00422.83559.06
H(−)CO2/H2S112271232.691229.265.450.9848.19207.52
S2230230.00230.003.250.9111.2728.11
S39696.0096.001.400.532.1110.55
H(−)acetS112061206.561206.805.410.9864.49202.84
S2317317.00317.143.640.9415.3941.54
S3184184.00184.202.590.805.0022.77
H(+)S1996996.30996.435.100.9733.28158.44
S2353353.00353.144.070.9625.1748.10
S3238238.00238.003.380.9211.9630.61
H(+)CO2/H2S1574574.00574.144.930.9866.4085.08
S2121121.00121.001.890.693.2713.60
S39797.0097.002.670.877.6310.80
H(+)acetS1530530.00530.004.720.9739.5769.79
S2335335.00335.143.540.9110.8344.04
S3241241.00241.283.830.9520.5331.56
Table 5. Variants of the media.
Table 5. Variants of the media.
Media codeTryptone (g L−1)Yeast Extract (g L−1)CH3COONa (g L−1)Under Gases
SZNNatural sediment with no additivesN2
H(−)CO2/H2000CO2/H2 (20:80 v/v)
H(−)acet002N2
H(+)0.50.50N2
H(+)CO2/H20.50.50CO2/H2 (20:80 v/v)
H(+)acet0.50.52N2
H(−)—basic mineral solution, H(+)—basic mineral solution supplemented with tryptone and yeast extract.

Share and Cite

MDPI and ACS Style

Szafranek-Nakonieczna, A.; Pytlak, A.; Grządziel, J.; Kubaczyński, A.; Banach, A.; Górski, A.; Goraj, W.; Kuźniar, A.; Gałązka, A.; Stępniewska, Z. Changes in the Substrate Source Reveal Novel Interactions in the Sediment-Derived Methanogenic Microbial Community. Int. J. Mol. Sci. 2019, 20, 4415. https://doi.org/10.3390/ijms20184415

AMA Style

Szafranek-Nakonieczna A, Pytlak A, Grządziel J, Kubaczyński A, Banach A, Górski A, Goraj W, Kuźniar A, Gałązka A, Stępniewska Z. Changes in the Substrate Source Reveal Novel Interactions in the Sediment-Derived Methanogenic Microbial Community. International Journal of Molecular Sciences. 2019; 20(18):4415. https://doi.org/10.3390/ijms20184415

Chicago/Turabian Style

Szafranek-Nakonieczna, Anna, Anna Pytlak, Jarosław Grządziel, Adam Kubaczyński, Artur Banach, Andrzej Górski, Weronika Goraj, Agnieszka Kuźniar, Anna Gałązka, and Zofia Stępniewska. 2019. "Changes in the Substrate Source Reveal Novel Interactions in the Sediment-Derived Methanogenic Microbial Community" International Journal of Molecular Sciences 20, no. 18: 4415. https://doi.org/10.3390/ijms20184415

APA Style

Szafranek-Nakonieczna, A., Pytlak, A., Grządziel, J., Kubaczyński, A., Banach, A., Górski, A., Goraj, W., Kuźniar, A., Gałązka, A., & Stępniewska, Z. (2019). Changes in the Substrate Source Reveal Novel Interactions in the Sediment-Derived Methanogenic Microbial Community. International Journal of Molecular Sciences, 20(18), 4415. https://doi.org/10.3390/ijms20184415

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop