Next Article in Journal
Free-Living Amoebae in Soil Samples from Santiago Island, Cape Verde
Previous Article in Journal
Supranational Assessment of the Quality of Probiotics: Collaborative Initiative between Independent Accredited Testing Laboratories
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Indicative Marker Microbiome Structures Deduced from the Taxonomic Inventory of 67 Full-Scale Anaerobic Digesters of 49 Agricultural Biogas Plants

1
Center for Biotechnology (CeBiTec), Bielefeld University, Universitätsstrasse 27, 33615 Bielefeld, Germany
2
Department Bioengineering, Leibniz Institute for Agricultural Engineering and Bioeconomy, Max-Eyth-Allee 100, 14469 Potsdam, Germany
3
Bioprocess Engineering, Otto von Guericke University, Universitätsplatz 2, 39106 Magdeburg, Germany
4
Bioprocess Engineering, Max Planck Institute for Dynamics of Complex Technical Systems, Sandtorstraße 1, 39106 Magdeburg, Germany
5
Microbiology, Anhalt University of Applied Sciences, Bernburger Straße 55, 06366 Köthen, Germany
6
Biochemical Conversion Department, DBFZ Deutsches Biomasseforschungszentrum Gemeinnützige GmbH, Torgauer Straße 116, 04347 Leipzig, Germany
7
The State Institute of Agricultural Engineering and Bioenergy, University of Hohenheim, Garbenstraße 9, 70599 Stuttgart, Germany
8
Institute of Agricultural Engineering, Kiel University, Max-Eyth-Str. 6, 24118 Kiel, Germany
9
Institute for Agricultural Engineering and Animal Husbandry, Bavarian State Research Center for Agriculture, Vöttinger Str. 36, 85354 Freising, Germany
*
Author to whom correspondence should be addressed.
Shared senior authorship.
Microorganisms 2021, 9(7), 1457; https://doi.org/10.3390/microorganisms9071457
Submission received: 11 June 2021 / Revised: 1 July 2021 / Accepted: 3 July 2021 / Published: 7 July 2021
(This article belongs to the Section Environmental Microbiology)

Abstract

:
There are almost 9500 biogas plants in Germany, which are predominantly operated with energy crops and residues from livestock husbandry over the last two decades. In the future, biogas plants must be enabled to use a much broader range of input materials in a flexible and demand-oriented manner. Hence, the microbial communities will be exposed to frequently varying process conditions, while an overall stable process must be ensured. To accompany this transition, there is the need to better understand how biogas microbiomes respond to management measures and how these responses affect the process efficiency. Therefore, 67 microbiomes originating from 49 agricultural, full-scale biogas plants were taxonomically investigated by 16S rRNA gene amplicon sequencing. These microbiomes were separated into three distinct clusters and one group of outliers, which are characterized by a specific distribution of 253 indicative taxa and their relative abundances. These indicative taxa seem to be adapted to specific process conditions which result from a different biogas plant operation. Based on these results, it seems to be possible to deduce/assess the general process condition of a biogas digester based solely on the microbiome structure, in particular on the distribution of specific indicative taxa, and without knowing the corresponding operational and chemical process parameters. Perspectively, this could allow the development of detection systems and advanced process models considering the microbial diversity.

1. Introduction

The production of biogas by anaerobic conversion of organic matter from various residues is a worldwide applied technology that can be integrated into sustainable bioeconomic concepts [1,2,3]. The production of biogas offers a number of benefits such as residue treatment [3], energy generation (electricity, heat and fuels) [3,4,5], mitigation of greenhouse gas emissions (in particular methane and nitrous oxide) [6,7], production of high-value fertilizers [8,9,10] as well as closing nutrient cycles [11,12].
However, in the past two decades, more than 9500 biogas plants in Germany have mainly been operated with energy crops (48.9%) and livestock manure (44.5%) [13], which are mainly co-digested to enhance the digestion process and thus increase biogas yields [14]. According to the current efforts of developing a bio-based circular economy (German National Bioeconomy Strategy, https://biooekonomie.de/nationale-biooekonomiestrategie, accessed on 5 July 2021), there will be a transition to a residue-based biogas production. This is, however, associated with high demands on the system technology and, above all, the control of the complex, microbial-mediated process [3].
The anaerobic conversion of organic matter requires the collaboration of hundreds to thousands of microbial taxa of which Bacteria and Archaea have been most intensively studied [15,16,17,18,19]. Currently, taxonomic profiling of microbial communities is mostly achieved by 16S rRNA gene amplicon sequencing (e.g., [20]), whereby most studies only considered a limited number of anaerobic digesters. Previously, the microbial communities of 21 full-scale biogas plants were taxonomically profiled whereby the community composition was mainly affected by the process temperature and the type of the supplied feedstocks [21]. Likewise, in a more recent study which investigated the microbial community structures of 20 different full-scale anaerobic digestion plants over one year, the respective compositional patterns were related to the supplied feedstocks in terms of agricultural residues, bio-waste, municipal solid waste and sewage sludge [15]. In addition to the temperature and the supplied feedstocks, also the ammonium/ammonia concentration plays an important role in shaping the microbial community structure [22,23].
Despite intensive research, most of the microorganisms that are involved in the anaerobic digestion process, their activity under specific process conditions, as well as their response to changing environmental parameters are still unknown (e.g., [15,16,24]). In this regard, some studies revealed that changes in the microbial community compositions are rather reflected in the variability of rarely existing taxa and not by changes in the relative abundances of members of the so-called core microbiome [15,25,26,27]. However, the question arises, what are specific taxa for prevalent process conditions. Moreover, several relationships of different microbial taxa have already been identified within anaerobic digestion systems, indicating the occurrence of similar niche adaptations and close collaborations (e.g., [24,28,29,30]) while the overall interrelationships of the microbial taxa are still less understood. For example, there is evidence that members of the bacterial phyla Cloacimonetes and Bacteroidetes combined with the archaeal genus Methanosaeta are supposed to be indicative taxa for well-running mesophilic biogas processes due to their sensitive response to increasing ammonium nitrogen and volatile fatty acid concentrations, but knowledge is still lacking on the microbial interconnectivity, e.g., regarding a functional compensation by other microbiome members [25,29,31]. The phylum Actinobacteria is another interesting example. Members of this phylum are well-known as human pathogens and drug producers [32,33], but their role in anaerobic digestion systems is still unclear, although they have been found in high abundances [16], particularly in manure-based small biogas plants [29]. Last but not least, members of the archaeal phylum Bathyarchaeota were also recently recorded within anaerobic digestion systems [34], whereby their system ecological role is nearly unclear as members of this phylum have a wide range of metabolic capacities, including hydrolysis, acetogenesis, methane metabolism, dissimilatory nitrogen and sulfur reduction as well as interactions with acetoclastic methanogens or heterotrophic bacteria [35]. Thus, elucidating the biogas microbiome, as described by the taxonomical, functional and ecological diversity, is still required and of great value in order to optimize management measures by means of sufficient methane yields, to derive benchmarks that indicate impending critical process conditions and in particular to accompany the transition to a residue-based biogas production by providing recommendations for sustainable biogas plant operations based on resilient microbiomes.
The objective of this study was to investigate the microbial diversity of a representative number of German biogas plants with stable process conditions which differ in terms of plant configuration, supplied feedstocks, process temperatures, and thus in the prevalent chemical parameters. We hypothesize that (i) the investigated microbial communities vary in their structural composition according to the prevalent environmental conditions whereby marker microbiome cluster can be derived for specific process conditions, (ii) the marker microbiome clusters are characterized by indicative taxa, and (iii) these indicative taxa can be used to deduce/assess the general process conditions of the corresponding biogas digesters.

2. Materials and Methods

2.1. Sample Sources, Sample Collection and Chemical Analyses

In this study the microbial diversity of 67 anaerobic digesters belonging to 49 different agricultural biogas plants was investigated. These biogas plants constitute a representative number of the biogas plants operated in Germany which were monitored over one year as part of the “Third National Biogas Measuring Programm”. The selection of the biogas plants was based in particular on the following criteria: (i) Different feedstock compositions, e.g., plants which solely converted plant biomasses or animal manure as well as co-fermentation biogas plants, (ii) different plant technologies, e.g., fermenter designs, and (iii) different process operation strategies based on, e.g., process temperature, organic loading rates and retention times.
Sampling points for the microbiological inventory were chosen after the biogas plants have shown stable process conditions for at least three months. Samples of the digester content were taken from the main digesters, reflecting the digester content at the current operational conditions [31,36]. To reduce the microbial activity, taken samples were stored on ice and directly transferred to the laboratory where aliquots were taken and stored at −20 C until subsequent chemical and microbiological analyses. The biogas plant operators provided information concerning the prevalent process parameters and, if available, the produced biogas amount and composition. For the entire sample set operational parameters were available, which include the reactor types and volumes, the supplied feedstocks (varying mixtures of energy crops and residues from livestock husbandry), organic loading rate (OLR), hydraulic retention time (HRT) and process temperature (Table S1). To compare the feedstock inputs, the absolute mass-based values were normalized to 100% for each digester. For all samples the following chemical analyses were carried out: Total solids (TS), volatile solids (VS), pH value, total ammonium nitrogen (TAN), and volatile fatty acids (VFA) according to Liebetrau et al. [37]. The free ammonia (NH 3 ) concentration was calculated as a function of the TAN concentration, the pH value and the temperature, according to Hansen et al. [38].

2.2. Taxonomic Profiling by Means of 16S rRNA Gene Amplicon Sequencing

The taxonomic profiling of the 67 biogas microbiomes was done as already described by Theuerl et al. [29]. Briefly, DNA was extracted from three replicates per sample using the FastDNA® Spin Kit for Soil (MP Biomedicals, France) in combination with the Genomic DNA Clean & Concentrator TM kit (Zymo Research, Irvine, CA, USA). The 16S rRNA gene amplicon library preparation was done applying the “16S Metagenomic Sequencing Library Preparation” protocol (Illumina Inc., San Diego, CA, USA) while using the universal primer pair Pro341F (5′-CCTACGGGNBGCASCAG-3′) and Pro805R (5′-GACTACNVGGGTATCTAATCC-3′) for amplification of the V3–V4 hypervariable regions [39]. The amplicon libraries were sequenced on the Illumina MiSeq platform applying the paired-end protocol for 300 bp reads and afterwards demultiplexed. The bioinformatic processing of the 16S rRNA gene amplicon sequences started with a quality control of the raw sequencing reads with FastQC. Afterwards the forward and reverse reads were merged with FLASH [40] and the primers were removed with cutadapt [41]. Sickle was used for trimming by quality of the merged reads [42]. The high-quality reads of each sample were subsampled to a given depth (50,000 reads) using seqtk (https://github.com/lh3/seqtk, accessed on 5 July 2021). Within the QIIME platform [43] usearch61 was used, for de novo (abundance based) and reference-based chimera detection. After this step, an open reference based operational taxonomic unit (OTU) clustering and a taxonomic assignment was accomplished while using the 16S rDNA SILVA database (release 132) as reference. Finally, triplicates of the biogas digesters were used to calculate median abundance values which were afterwards normalized to 100%.

2.3. Statistical Data Analysis

In order to evaluate correlations among and between biotic and abiotic system parameters, non-metric multidimensional scaling (NMDS) [44,45,46] was performed with the R Project for Statistical Computing [47] using the package “vegan” [48] and “gplots” [49]. The distance matrix was calculated using the McQuitty algorithm [50]. The calculated (dis)similarities were first visualized as a tree diagram to define the existing clusters of the investigated anaerobic digesters based on their taxonomic profiles. Subsequently, the calculated distance matrix was used to visualize colored NMDS ordination plots which were combined with environmental vectors based on the process parameters calculated by the function “envfit” [48]. The results were sorted according to the R 2 values while only vectors with R 2 > 0.3 and p < 0.001 were considered.
To identify taxa which are indicative for specific environmental conditions, an indicator species analysis (ISA) according to [51] was performed using the “indicspecies” package of R with a point biserial correlation coefficient as statistic value [52,53]. Considering only indicative taxa with a point biserial correlation coefficient above 0.5 and p-values below 0.0005.
Pearson correlations between the process parameters and the indicator species abundances were calculated in R using the “cor” function of the package “stats” [47]. The results were sorted according to the correlation values, while moderate correlations between ±0.50 and ±0.75 and strong correlations higher/lower than ±0.75 were considered.

3. Results and Discussion

3.1. Operational and Chemical Characteristics of the 67 Analyzed Biogas Digesters

The investigated 67 anaerobic digesters belonging to 49 agricultural full-scale biogas plants differed in their prevalent operational process conditions (Figure 1, Table S1). The main feedstocks were maize and grass silage (in 60 and 48 digesters) as well as liquid and solid cattle manure (in 36 and 20 digesters) which ranged from 2.1 to 99.3%, 0.3 to 69.0%, 2.6 to 95.5% and 4.6 to 30.6% of the total feedstock supply, respectively. Further rarely or occasionally used feedstocks were sugar beet silage in ten digesters (1.4–13.7%), whole crop rye silage in 24 digesters (0.3–33.5%), cereals in 28 digesters (0.3–19.4%), liquid and solid swine manure in eight digesters (8.5–70.5%), chicken/poultry manure in eight digesters (3.4–25.0%), sheep manure in one digester (4.2%) and horse manure in one digester (26.1%). The relative shares of the used feedstocks are in good agreement with the average feedstock compositions of German biogas plants [13].
The OLRs of the biogas digesters ranged from 0.4 to 21.2 kg VS m 3 d 1 with corresponding HRTs from 3.6 to 475.0 days. The process temperature varied between 36.0 C and 56.3 C, whereby 57 of the 67 analyzed anaerobic digesters were operated at a process temperature between 39 C and 49 C, a range between the optimal condition for mesophilic (33–38 C) or thermophilic (50–60 C) process operation, respectively [23,54]. Three biogas plants were operated at mesophilic conditions (≤38 C) and seven biogas plants at thermophilic conditions (≥50 C). Due to the broad variety of the operational process parameters, differences in the chemical process parameters were expected. Hence, the pH values ranged from 6.1 to 8.3, the TKN concentrations from 2.6 to 7.7 gL 1 , the TAN concentrations from 1.2 to 4.7 L 1 , the NH 3 concentrations from 0.5 to 157.5 mgL 1 and the VFA concentrations from below the detection limit up to 14.3 gL 1 (Figure 1, Table S2).

3.2. Taxonomic Profiling of the 67 Analyzed Anaerobic Digesters

This study investigated the taxonomic diversity of 67 anaerobic digesters originating from 49 different agricultural full-scale biogas plants. The taxonomic profiles of the biogas microbiomes were determined by 16S rRNA gene amplicon sequencing. A total of 17,556 different OTUs with an average number of 2260 ± 528 (Table S3) were identified per digester. The taxonomic profiles revealed that 89.4% ± 6.5% of the detected taxa were assigned to the domain Bacteria, 9.4% ± 6.3% to the domain Archaea and 1.2% ± 0.6% to taxonomically unclassified organisms. The dominant bacterial phyla were Firmicutes (60.4% ± 8.2%), Bacteroidetes (15.3% ± 5.1%), Cloacimonetes (3.2% ± 3.7%), Actinobacteria (1.6% ± 3.0%), Tenericutes (1.4% ± 0.6%) and Atribacteria (1.2% ± 0.8%) (Figure S1). On the order level, Clostridiales were dominantly found in all analyzed digesters with 32.5% ± 10.8%, followed by “uncultured Clostridia MBA03” with 13.4% ± 9.1%, Bacteroidales with 10.8% ± 4.4%, “uncultured Clostridia DTU014” with 5.3% ± 4.6%, Sphingobacteriales with 4.4% ± 3.4%, Lactobacillales with 3.4% ± 3.8% and Cloacimonadales with 3.2% ± 3.7% (Figure 2). Among the domain Archaea, members of the phylum Euryarcheota were predominantly found in all analyzed anaerobic digesters with 9.3% ± 6.3%. Other archaeal phyla, such as Crenarchaeota, Nanoarchaeota and Diapherotrites, were found with low relative abundances of maximal 0.2%, 0.1% and 0.03%, respectively. On the order level, Methanosarcinales dominated the analyzed digesters with 7.9% ± 5.8%, followed by Methanobacteriales with 1.3% ± 2.3% (Figure 2). The recorded taxonomic profiles, in general, correspond to the expected taxonomic composition of microbial communities in biogas digesters [15,21,22,55].

3.3. Marker Microbiome Clusters Depending on Prevalent Process Conditions

In order to identify potential marker microbiome clusters for prevalent process conditions, different statistical analyses were carried out. In a first step, a McQuitty (dis)similarity matrix was calculated and visualized in a tree diagram using the taxonomic profiles of the 67 investigated anaerobic digesters (Figure 3). Due to the taxonomic diversity of the analyzed anaerobic digesters, they were divided into two main clusters I and II and two groups of outliers. Cluster I represents digesters featuring process temperatures higher than 45 C, whereas cluster II comprises digesters with process temperatures lower than 45 C. Cluster II was further separated into two distinct subclusters IIa and IIb (Figure 3).
Subsequently, the defined microbiome clusters were used to visualize colored NMDS ordination plots which arrange objects, in this case the recorded taxonomic profiles, as close as possible according to their (dis)similarities. The identified clusters were coloured according to the defined clusters of the tree diagram: Cluster I in green, cluster IIa in red, cluster IIb in blue and the outliers in grey (Figure S2). Many of the investigated samples were assigned to one of the specific microbiome clusters I, IIa and IIb based on the calculated McQuitty (dis)similarity matrix visualized as tree diagram (Figure 3). However, the NMDS ordination plot revealed, in addition to the outliers of the tree diagram (B11, B37, B42, B43, B49), samples which also lie outside of the identified clusters (B12, B14, B19, B25, B29, B35, B38, B44) (Figure S2). These samples will be summarized as group of outliers for the following analyses.
As expected, the NMDS analysis revealed cluster I (green in Figure 3 and Figure 4) for which the process temperature had a significant impact (R 2 = 0.74, p = 0.001) on the formation of the microbial community structure. This cluster consists of 19 anaerobic digesters from 11 biogas plants (B07, B13, B15, B18, B20, B21, B24, B33, B39, B41, B45); all of them showed a process temperature higher than 45 C (Table S1). The separation of this cluster from the other anaerobic digesters is in accordance with previously published studies which identified the process temperature as one of the main factors for the shaping of microbial communities within anaerobic digesters [22,28,31,56,57].
Among the available operational and chemical process parameters, the process temperature had the strongest impact on separating the two main microbiome clusters I and II (Figure 3). Due to their lower correlation values, all other process parameters can only to a certain degree be used to determine the influence of the abiotic environment on the taxonomic diversity. Considering a correlation tendency, also the amount of liquid cattle manure (R 2 = 0.32, p = 0.001), rye and maize silage (each with a R 2 of 0.30, p = 0.001), the TAN concentration (R 2 = 0.40, p = 0.001), the TKN concentration (R 2 = 0.31, p = 0.001) as well as the NH 3 concentration (R 2 = 0.51, p = 0.001), might explain the taxonomic composition of the clusters (Figure 4). Most probably these low correlations are related to the chemical similarities of the used feedstocks in terms of energy crops and residues from livestock husbandry in combination with an almost similar process operation (Table S1). This is in accordance with previous studies which have shown stronger significant correlations of the microbial community structure with prevalent process parameters, e.g., the supplied feedstocks. Hence, differentiations occur when agricultural biogas plants are compared with anaerobic digestion plants using either chemically complex and/or heterogeneous feedstocks in function of time, (e.g., biowaste which is the most variable feedstock as its chemical composition, for example, is subjected to seasonal differences) or more homogeneous feedstocks (e.g., wastewater sludge and industrial waste) [15,22,24,28,31]. Independently of these similarities within the process parameters of the investigated biogas plants, the question arises whether indicative microbiome clusters can be derived for specific process conditions. In this regard, the NMDS analysis revealed the two subclusters, IIa and IIb, of biogas plants operated at a process temperature lower than 45 C. The subcluster Ia (marked in red in Figure 3 and Figure 4) consists of 12 anaerobic digesters from ten biogas plants (B01, B04, B09, B10, B22, B27, B28, B30, B34, B46) which were operated with a balanced mixture of energy crops (mainly maize and grass silage and occasional addition of sugar beet silage or cereals) and residues from livestock husbandry (liquid and solid cattle and swine manure) with a process temperature between 36 C and 43 C (Table S1). These plants showed none of the commonly known critical features in their chemical parameters (Table S2) such as elevated TAN or VFA concentrations [19,58,59]. Compared to this, subcluster IIb (marked in blue in Figure 3 and Figure 4) consists of 21 anaerobic digesters from 15 biogas plants (B02, B03, B05, B06, B08, B16, B17, B23, B26, B31, B32, B36, B40, B47, B48) which were operated at similar feedstock and temperature conditions as the digesters from cluster IIa, but differ by slightly higher values of the pH (7.2 ± 0.2 vs. 7.8 ± 0.7), the VFA concentration (0.07 ± 0.07 gL 1 vs. 0.54 ± 0.70 gL 1 ) and the TAN concentration (1.8 ± 0.4 gL 1 vs. 2.2 ± 0.7 gL 1 ) (Table S2).
For the previously identified group of outliers, other environmental factors than the prevalent process temperature or slight variations of certain chemical parameters led to a different composition of the corresponding microbial communities. With respect to the supplied feedstocks, it becomes obvious that the corresponding microbiomes are either affected by comparatively high amounts of solid and liquid cattle manure (samples B12, B14, B19, B29, B42, B43, B44) or by relatively high amounts of chicken and poultry manure resulting in elevated pH values as well as TAN and NH 3 concentrations (B11, B25, B38, B35, B37, B49) (Figure 4, Tables S1 and S2).
In accordance with these results, former cluster (PCoA, NMDS) and co-occurrence network analyses indicate the existence of specific microbial assemblages for different anaerobic digestion systems [15,22,28,31] while it has to be considered that these studies compared agricultural biogas plants with anaerobic digestion systems converting either biowaste or wastewater which significantly differ. In contrast, this study compared agricultural biogas plants that are more similar regarding their input feedstocks and process parameters (see Section 3.1). At this point, the question arises whether there are indicative taxa which are responsible for the identified clustering.

3.4. Indicative Taxa for Prevalent Process Conditions

In order to find indicative taxa for the three identified clusters I, IIa, and IIb as well as for the group of outliers, an Indicator Species Analysis (ISA) was calculated based on the taxonomic profiles of all digester samples. Considering only indicative taxa with a point biserial correlation coefficient above 0.5 and p-values below 0.0005, overall 253 taxa were identified as differentially abundant between the three defined clusters and cluster combinations: 23 indicative taxa were identified for cluster I, 68 taxa for cluster IIa, and four taxa for cluster IIb (Table S4). For the outlier group, 123 taxa were identified as potentially indicative, which have to be handled carefully due to the broad dispersion of the outlier samples. The remaining 35 recorded indicative taxa refer to combinations of two or three clusters, whereby 28 of them are indicative for the combination of cluster IIa and IIb.
The microbial composition of cluster I, which consists of anaerobic digesters with elevated process temperatures between 45 C and 56 C, is mainly characterized by members of uncultured genera of the class Clostridia (Figure 5, Table S4). In particular, two groups, namely “uncultured Clostridia MBA03” and “uncultured Clostridia DTU014”, showed 24.0% ± 4.9% and 11.2% ± 2.5%, respectively, the highest relative abundances of the indicator taxa. Further indicator taxa with high abundances were assigned to the genera “uncultured Lentimicrobiaceae” with 8.4% ± 2.7%, Proteiniphilum with 6.9% ± 2.0%, Defluviitoga with 2.6% ± 2.5%, “uncultured Clostridia” with 1.8% ± 0.5% and Halocella with 1.7% ± 0.9%. The archaeal community was solely represented by the genus Methanobacterium with a relative abundance of 0.9% ± 0.4% (Figure 5, Table S4). These findings are in general accordance with previously published studies which highlighted the process temperature as one of the main factors for shaping of the microbial community structure within anaerobic digesters. These microbial communities are predominated by members of the phylum Firmicutes, especially of the class Clostridia and the phylum Thermotogae as well as hydrogenotrophic methanogens such as the archaeal genus Methanobacterium [22,24,28,31,56,60]. Surprisingly, members from the phylum Bacteroidetes by means of the genera “uncultured Lentimmicrobiaceae” and Proteiniphilum were identified as indicative taxa for cluster I. It is commonly reported that the relative abundance of members of the phylum Bacteroidetes is significantly lowered at thermophilic process conditions [22,23,28,61] which is supported as the known and described species of these both genera have their optimal growth conditions at 30–37 C (Lentimicrobium saccharophilum, [62]), at 37 C (Proteiniphilum acetatigenes, [63]) or at 35–40 C (Proteiniphilum saccharofermentans, [64]). However, as a reliable assignment of the detected sequences is only possible at the genus level [20,65,66], both mentioned genera probably comprise thermophilic or at least thermo-tolerant species.
However, Pearson correlations of the indicative taxa with the process parameters showed moderate positive correlations (between 0.50 and 0.75) with the amount of whole crop rye silage and/or the NH 3 concentration in particular for the bacterial genera “uncultured Clostridia MBA03”, “uncultured Clostridia DTU014”, “uncultured Lentimicrobiaceae”, Halocella and the archaeal genus Methanobacterium (Table S6). Taking into further consideration that the bacterial genera “uncultured Clostridia MBA03” and “uncultured Clostridia DTU014” showed a moderate negative correlation (−0.54 and −0.52, respectively) with the feedstock liquid cattle manure (Table S6), it might be assumed that the members of these both genera are involved in the degradation of plant biomass. This strengthens their importance of being indicative for the process condition of the anaerobic digestion systems they are living in.
Cluster IIa consists of biogas plants characterized by very low values of the prevalent chemical process parameters (Table S2). The ISA revealed that most of the 68 recognized indicative taxa for this cluster belong to the bacterial phyla Bacteroidetes (16.4%), Firmicutes (16.4%) and Patescibacteria (11.9%), followed by the phyla Spirochaetes (9%), Chloroflexi (9%), Planctomycetes (4.5%) as well as Atribacteria, Cloacimonetes, Kiritimatiellaeota and Synergistetes (each 3%) (Figure 5, Table S4). This picture is also supported considering the 22 most abundant genera: Six belong to the phylum Firmicutes (with relative abundances between 0.2% ± 0.1% and 7.6% ± 2.8%), five to the phylum Bacteroidetes (with relative abundances between 0.6% ± 1.2% and 1.9% ± 1.1%), and two genera belong to each of the phyla Patescibacteria, Spirochaetes and Chloroflexi (with relative abundances between 0.2% ± 0.2% and 1.3% ± 1.2%) (Table S4). With one genus each, the bacterial phyla Cloacimonetes (Candidatus Cloacimonas, 1.0% ± 0.8%), and Synergistetes (“uncultured Synergistaceae”, 0.5% ± 0.1%) also belong to the most abundant indicative taxa of cluster IIa. For the archaeal community, the genus Methanosaeta was identified as an indicative taxon with a relative abundance of 6.3% ± 2.9% (Figure 5, Table S4). This is partly in accordance with previously published results by Theuerl et al. [29,31], where the occurrence of members from the bacterial phyla Bacteroidetes and Cloacimonetes and the archaeal genus Methanosaeta (syn. Methanothrix) were related to unstressed anaerobic digestion systems with low concentrations of potential process inhibitory factors. In comparison, this study showed that much more members of various phyla are potentially indicative for unstressed process conditions, besides the already known taxa. In particular, these include members of the phyla Firmicutes, Patescibacteria, Spirochaetes and Chloroflexi.
However, out of the 68 indicative taxa for cluster IIa, for ten taxa belonging to the phyla Bacteroidetes, Firmicutes, Patescibacteria, Spirochaetes, Atribacteria and Kiritimatiellaeota positive correlations (from 0.51 to 0.79) were detected, in particular for the process parameters liquid cattle manure and swine manure (Table S6). In addition, two indicative taxa assigned inter alia to the nine most abundant indicative taxa showed moderate negative correlations with the NH 3 concentration (Table S6) which supports the assumption that the identified taxa are indicative for mesophilic unstressed anaerobic digestion systems with low concentrations of potential process inhibitory factors, here in terms of the NH 3 concentration.
Cluster IIb includes digester samples from co-fermentation plants, which differ from cluster IIa by slightly elevated pH values, VFA and TAN concentrations (Table S2). The corresponding microbial communities were predominantly characterized by the indicative genus Petrimonas (0.9% ± 0.4%) of the phylum Bacteroidetes (Table S4). The remaining three indicative genera “uncultured Chthonomonadales”, Ruminiclostridium 6 and “uncultured BRC1” showed relative abundances below 0.01%.
Since the process conditions of cluster IIa and IIb are very similar and only differ slightly in values for the parameters pH, TAN and VFA, it was expected that a large group of indicative taxa were identified for the combination of these two clusters vs. cluster I. Out of the 28 indicative taxa for cluster IIa and IIb, 43% belong to the phylum Firmicutes, 29% to the phylum Bacteroidetes and 7% to the phylum Cloacimonetes. The phyla Actinobacteria, Chloroflexi, Fibrobacteres, Synergistetes, Tenericutes and Verrucomicrobia are each represented by one genus among the identified indicative taxa. The most abundant genera were Sedimentibacter (Firmicutes, 6.6% ± 2.3%), “uncultured Ruminococcaceae” (Firmicutes, 2.0% ± 0.5%), Herbinix (Firmicutes, 1.6% ± 0.6%), “uncultured Bacteroidales UCG-001” (Bacteroidetes, 2.1% ± 1.4%), Fermentimonas (Bacteroidetes, 0.9% ± 0.4%) and “uncultured Cloacimonadaceae W5” (Cloacimonetes, 3.3% ± 2.8%) (Table S4).
Compared to cluster I, correlations of the indicative taxa of cluster II (a and b) with the process parameters showed an opposite behavior by means of negative correlations with the NH 3 concentration. Hereby, 17 indicative taxa showed moderate negative correlation values between −0.50 and −0.62 for the NH 3 concentration (Table S6), indicating that these taxa from cluster II are potentially sensitive to elevated NH 3 conditions.
For the group of outliers, 123 indicator taxa were identified. These taxa have to be considered carefully due to the broad dispersion of the outlier samples. However, 40% of these taxa were assigned to the phylum Actinobacteria and 30% to the phylum Proteobacteria (Table S4). In no other cluster Proteobacteria were identified as indicator taxa and only four genera of the phylum Actinobacteria were identified in cluster II. Additionally, approximately 50% of the indicative taxa assigned to the phylum Actinobacteria showed moderate positive correlations (from 0.50 to 0.70) with the parameters solid cattle manure, chicken manure, poultry manure as well as horse manure (Table S6). Thus, it can be assumed that taxa of these phyla seem to be indicators for biogas digesters with higher amounts of animal manure. Regarding the high occurrence of members from the phylum Actinobacteria, this assumption is supported by several comparative studies including anaerobic digestion systems with higher shares of excrements (manure from livestock husbandry or wastewater sludge) [15,22,29,55,67]. However, studies considering the microbial community structure of anaerobic digestion systems compared with the gastrointestinal tract of farm animals as well as the corresponding manures, in order to elucidate the system specificity of Actinobacteria and especially Proteobacteria and their origin, are still lacking.
In conclusion, the taxonomic profiling enabled the definition of three biogas microbiome clusters. Within these clusters, taxa and their specific abundance ranges were identified which are indicative for the corresponding type of process condition.

3.5. Clustering of the Biogas Microbiomes by Indicative Taxa and Generally Occurring Taxa

In order to verify the identified indicative taxa as deduced from ISA for their usability to assess the identified clusters based on the entire taxonomic microbial diversity, further NMDS analyses were calculated (Figure 6A). These analyses revealed that it is possible to distinguish the same three main clusters and the group of outliers which could be distinctly separated from each other by just using the 253 indicative genera and their relative abundances instead of 1021 genera from the entire microbial community. Besides the mentioned indicative taxa (see Section 3.4), also 415 generally occurring genera were observed within the taxonomic profiles of the analysed digesters. They mainly show similar relative abundances within the three clusters and the group of outliers, or between cluster combinations. The most abundant (>1%) generally occurring genera were Ruminiclostridium, Clostridium sensu stricto 1, “uncultured Syntrophomonadaceae” and Ruminiclostridium 1 (Table S5). Further generally occurring taxa with partly quite different average relative abundances within the clusters were “uncultured Rikenellaceae” (cluster I: <0.1%, cluster IIa: 1.5%, cluster IIb: 0.7%, outlier: 0.9%), Methanobrevibacter (cluster I: <0.1%, cluster IIa: 0.4%, cluster IIb: 0.7%, outlier: 1.3%) and Pseudomonas (cluster I: 0.2%, cluster IIa: <0.1%, cluster IIb: <0.1%, outlier: 1.2%). These taxa were not assigned as indicator taxa, which could be explained by abundance variations in the respective clusters. As expected, it could be shown that the generally occurring taxa do not enable any clustering (Figure 6B). This indicates that certain process conditions are not reflected by members of the core-microbiome, but rather by specific taxa, which are adapted to the prevailing conditions as it was shown in wastewater treatment systems or anaerobic digestion plants (e.g., [15,27,68]).
However, from an ecological point of view according to the Pareto principle, generalists which exist in high and even abundances under various conditions are supposed to be responsible for 80% of the energy flux and hence they are important drivers for the overall process functioning [15,31,69]. Consequently the specialists, by means of the identified indicative taxa, of the biogas microbiomes occupy specific ecological niches and hence can be used to assess certain process conditions more precisely than the whole microbiome or the usually used process parameters. It is hypothesized that the reaction of the microbial communities, especially the indicative taxa, enables an earlier response to changing environmental conditions, than the measurable chemical process parameters show elevated or critical values. Thus process specific indicative taxa might be used as sensitive markers for an early insight into the process condition of agricultural biogas plants in the future.

4. Conclusions

In this study, 67 microbiomes originating from 49 agricultural full-scale biogas plants converting energy crops and/or residues from livestock husbandry were taxonomically investigated. As it was hypothesized, the investigated microbial communities vary in their structural composition according to specific process conditions whereby three marker microbiome clusters, I, IIa and IIb, and one group of outliers were identified. For the first time, this study has shown that these marker microbiome clusters are characterized by specific genera, here by 253 indicative taxa which are differentially abundant within the defined clusters. Correlation analyses of the 253 indicative taxa against the prevailing process parameters supported and specified the occurrence of microbial indicators for certain process conditions, in particular the presence of potential indicators for mesophilic unstressed conditions as well as the presence of taxa that are adapted to higher NH 3 concentrations with an overall adequate process stability. For the first time, this study showed the occurrence of potential sensitive and resistant/resilient taxa e.g., to elevated NH 3 concentrations, which leads to the assumption that such indicative taxa provide better information of the actual process conditions than the entire microbiome (with more than 1000 different genera) or the usually used chemical process parameters such as the acid and/or the ammonium/ammonia concentration. However, at this point the question arises how these indicative taxa and hence the entire system will respond to changing environmental conditions, especially the transition to a residue-based biogas production and the related frequently varying process conditions, since they are obviously highly adapted to the prevalent process conditions of their anaerobic digestion system. To follow up, there is a need to find out, whether these indicative taxa have specific tolerance ranges to the prevailing environmental conditions in order to develop easy-to-use detection systems which information can be used for the development of advanced, data-driven process models that require benchmarks for both stable/resilient and impending critical process conditions in order to recommend sustainable biogas plant operations.
From a scientific point of view, insights into the metabolic potential of entire microbial communities from various anaerobic digestion systems of the here defined microbiome clusters through metagenomic analyses as well as of genetic specifications of the occurring, in particular yet unknown microorganisms through genome reconstruction approaches combined with cultivation approaches would enable a deeper understanding of the interconnections between process conditions and the composition of the microbial community with its functional potential. Additionally, the actually realized functions as deduced from metaproteomic analyses would enable the analysis of the metabolic solutions of cluster specific microbial communities to fulfill the anaerobic digestion chain in biogas plants.
Taking into account that the currently available methods for investigating the microbial diversity are highly complex in terms of sample preparation and, in particular, data evaluation and interpretation, identifying microbial indicators (at the taxonomic and functional level) for specific process conditions is a precondition for the development of rapid and economical microbial detection systems and hence new process models. These can be used by plant operators and consultants for monitoring, assessing and managing all types of anaerobic digestion systems using e.g., sewage sludge, residues from livestock husbandry, landscape management, from food processing and consumption as well as from biorefineries, or the organic fraction of municipal solid waste.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/microorganisms9071457/s1: Figure S1: Taxonomic profiles at phylum level for the 67 investigated biogas microbiomes from 49 agricultural, full-scale biogas plants. Phyla with a maximal relative abundance less than 1% were summarized. Taxa which were not taxonomically assigned at phylum level were summarized as “unclassified”. Figure S2: Non-metric multidimensional scaling (NMDS) based on the taxonomic diversity of the 67 analyzed digester samples from the 49 agricultural, full-scale biogas plants. The color code of the clusters is deduced from the tree diagram in Figure 3: I = green, IIa = red, IIb = blue, outliers = gray. The given environmental vectors symbolize the amount of maize silage (R 2 = 0.30, p = 0.001), whole crop rye silage (R 2 = 0.30, p = 0.001), liquid cattle manure (R 2 = 0.32, p = 0.001), the process temperature (R 2 = 0.74, p = 0.001), total solids (TS) (R 2 = 0.33, p = 0.001), the total ammonium nitrogen (TAN) (R 2 = 0.40, p = 0.001), the total Kjeldahl nitrogen (TKN) (R 2 = 0.31, p = 0.001) and the free ammonia nitrogen concentration (NH 3 ) (R 2 = 0.51, p = 0.001). Figure S3: Taxonomic profiles at genus level for the 67 investigated biogas microbiomes from 49 agricultural, full-scale biogas plants sorted by their respective cluster. Genera with a maximal relative abundance <5.5% were summarized. Taxa which were not assigned at genus level were summarized as “unclassified”. Table S1: Overview of the operational parameters of the 67 analyzed biogas digesters from 49 agricultural, full-scale biogas plants. The used feedstocks were normalized to 100% for each biogas digester. Table S2: Overview of the chemical parameters of the 67 analyzed biogas digesters from 49 agricultural full-scale biogas plants. VFA = volatile fatty acids. Table S3: Overview of the detected operational taxonomic units (OTUs) as well as the calculated α -diversity indices in terms of Shannon, Simpson and Chao1 for the investigated 67 biogas digesters from 49 agricultural, full-scale biogas plants. Table S4: The 23 most abundant indicative taxa above a maximal relative abundance of 0.5%, their average abundances and standard deviations within the defined clusters based on the taxonomic profiles. Statistically significant differences were calculated using the “indicspecies” package of R with a point biserial correlation coefficient as statistic value [De Cáceres & Legendre, 2009; De Cáceres et al., 2011]. Table S5: The most abundant generally occurring taxa above a maximal relative abundance of 0.5%, their average abundances and standard deviations within the defined clusters based on the taxonomic profiles. Table S6: Correlations of 123 indicator taxa against 18 process parameters. Shown are only moderate correlations between ±0.50 and ±0.75 and strong correlations above/below ±0.75.

Author Contributions

Conceptualization: J.H., S.T., A.P., A.S.; Investigation (sampling, process engineering and chemical parameters): M.P., B.H., T.M., M.E.; Formal analysis and Data curation (process engineering and chemical parameters): M.P., B.H., T.M., M.E., J.K., J.H., S.T.; Methodology (16S rRNA gene amplicons): J.H., A.S.; Investigation, Formal analysis, Data curation and Software (16S rRNA gene amplicons): J.H.; Investigation and Formal analysis (NMDS, ISA, Correlations): J.H., J.K., S.T.; Visualization: J.H., S.T.; Writing—original draft preparation (Intro): A.S., S.T., J.H.; Writing—original draft preparation (MM): J.H., J.K., S.T.; Writing—original draft preparation (RD): J.H., S.T.; Writing—review and editing: J.H., J.K., D.B., M.P., B.H., T.M., M.E., A.P., A.S., S.T.; Resources: A.P., A.S., S.T.; Supervision: A.P., A.S., S.T.; Project administration: S.T., A.P., A.S., D.B.; Funding acquisition: A.P., A.S., S.T., D.B.; All authors have read and agreed to the published version of the manuscript.

Funding

J.H., J.K., S.T., M.P., B.H., T.M. and M.E. were funded by the German Federal Ministry of Food and Agriculture (BMEL) via the Fachagentur für Nachwachsende Rohstoffe e.V. (FNR) of the joint research project “Biogas measuring programm III” (FKZ 22403915, 22404015, 22404115, 22403515, 22403615, 22403715, 22403815). A.S. was funded by the European Union’s Horizon 2020 research and innovation program (grant agreement No. 818431, SIMBA). A.P. was funded by the BMBF funded project “Bielefeld-Gießen Center for Microbial Bioinformatics–BiGi” (Grant number 031A533, 031L0103) within the German Network for Bioinformatics Infrastructure (de.NBI).

Data Availability Statement

The raw 16S rRNA gene amplicon sequences are deposited in the European Nucleotide Archive (ENA) under the study accession PRJEB39821.

Acknowledgments

Beate-Kristin Kröck and Katharina Hanuschka are gratefully acknowledged for their valuable support with the laboratory work. Sequencing support provided by Tobias Busche and the IIT Biotech GmbH is also gratefully acknowledged. The project partners of the joint research project “Biogas Measuring Program III” are gratefully acknowledged for providing samples and data on operating conditions of the biogas digesters. We acknowledge support for the publication costs by the Open Access Publication Fund of Bielefeld University.

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.

References

  1. Hagman, L.; Blumenthal, A.; Eklund, M.; Svensson, N. The role of biogas solutions in sustainable biorefineries. J. Clean. Prod. 2018, 172, 3982–3989. [Google Scholar] [CrossRef] [Green Version]
  2. Liebetrau, J.; Sträuber, H.; Kretzschmar, J.; Denysenko, V.; Nelles, M. Anaerobic digestion. Biorefineries 2017, 166, 281–299. [Google Scholar]
  3. Theuerl, S.; Herrmann, C.; Heiermann, M.; Grundmann, P.; Landwehr, N.; Kreidenweis, U.; Prochnow, A. The future agricultural biogas plant in Germany: A vision. Energies 2019, 12, 396. [Google Scholar] [CrossRef] [Green Version]
  4. Hahn, H.; Krautkremer, B.; Hartmann, K.; Wachendorf, M. Review of concepts for a demand-driven biogas supply for flexible power generation. Renew. Sustain. Energy Rev. 2014, 29, 383–393. [Google Scholar] [CrossRef]
  5. Lauer, M.; Thrän, D. Flexible biogas in future energy systems—Sleeping beauty for a cheaper power generation. Energies 2018, 11, 761. [Google Scholar] [CrossRef] [Green Version]
  6. Baral, K.R.; Jégo, G.; Amon, B.; Bol, R.; Chantigny, M.H.; Olesen, J.E.; Petersen, S.O. Greenhouse gas emissions during storage of manure and digestates: Key role of methane for prediction and mitigation. Agric. Syst. 2018, 166, 26–35. [Google Scholar] [CrossRef]
  7. Mohankumar Sajeev, E.P.; Winiwarter, W.; Amon, B. Greenhouse gas and ammonia emissions from different stages of liquid manure management chains: Abatement options and emission interactions. J. Environ. Qual. 2018, 47, 30–41. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Möller, K.; Müller, T. Effects of anaerobic digestion on digestate nutrient availability and crop growth: A review. Eng. Life Sci. 2012, 12, 242–257. [Google Scholar] [CrossRef]
  9. Schröder, P.; Beckers, B.; Daniels, S.; Gnädinger, F.; Maestri, E.; Marmiroli, N.; Mench, M.; Millan, R.; Obermeier, M.; Oustrière, N.; et al. Intensify production, transform biomass to energy and novel goods and protect soils in Europe—A vision how to mobilize marginal lands. Sci. Total Environ. 2018, 616, 1101–1123. [Google Scholar] [CrossRef] [PubMed]
  10. Nkoa, R. Agricultural benefits and environmental risks of soil fertilization with anaerobic digestates: A review. Agron. Sustain. Dev. 2014, 34, 473–492. [Google Scholar] [CrossRef] [Green Version]
  11. Arthurson, V. Closing the global energy and nutrient cycles through application of biogas residue to agricultural land–potential benefits and drawback. Energies 2009, 2, 226–242. [Google Scholar] [CrossRef] [Green Version]
  12. Valentinuzzi, F.; Cavani, L.; Porfido, C.; Terzano, R.; Pii, Y.; Cesco, S.; Marzadori, C.; Mimmo, T. The fertilising potential of manure-based biogas fermentation residues: Pelleted vs. liquid digestate. Heliyon 2020, 6, e03325. [Google Scholar] [CrossRef] [PubMed]
  13. Daniel-Gromke, J.; Rensberg, N.; Denysenko, V.; Stinner, W.; Schmalfuß, T.; Scheftelowitz, M.; Nelles, M.; Liebetrau, J. Current developments in production and utilization of biogas and biomethane in Germany. Chem. Ing. Tech. 2018, 90, 17–35. [Google Scholar] [CrossRef]
  14. Pan, S.Y.; Tsai, C.Y.; Liu, C.W.; Wang, S.W.; Kim, H.; Fan, C. Anaerobic Co-Digestion of Agricultural Wastes Towards Circular Bioeconomy. iScience 2021, 24, 102704. [Google Scholar] [CrossRef]
  15. Calusinska, M.; Goux, X.; Fossépré, M.; Muller, E.E.; Wilmes, P.; Delfosse, P. A year of monitoring 20 mesophilic full-scale bioreactors reveals the existence of stable but different core microbiomes in bio-waste and wastewater anaerobic digestion systems. Biotechnol. Biofuels 2018, 11, 1–19. [Google Scholar] [CrossRef]
  16. Campanaro, S.; Treu, L.; Rodriguez-R, L.M.; Kovalovszki, A.; Ziels, R.M.; Maus, I.; Zhu, X.; Kougias, P.G.; Basile, A.; Luo, G.; et al. New insights from the biogas microbiome by comprehensive genome-resolved metagenomics of nearly 1600 species originating from multiple anaerobic digesters. Biotechnol. Biofuels 2020, 13, 1–18. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. De Vrieze, J.; Ijaz, U.Z.; Saunders, A.M.; Theuerl, S. Terminal restriction fragment length polymorphism is an “old school” reliable technique for swift microbial community screening in anaerobic digestion. Sci. Rep. 2018, 8, 1–12. [Google Scholar] [CrossRef]
  18. Schnürer, A. Biogas production: Microbiology and technology. Anaerobes Biotechnol. 2016, 156, 195–234. [Google Scholar]
  19. Theuerl, S.; Klang, J.; Prochnow, A. Process disturbances in agricultural biogas production—Causes, mechanisms and effects on the biogas microbiome: A review. Energies 2019, 12, 365. [Google Scholar] [CrossRef] [Green Version]
  20. Hassa, J.; Maus, I.; Off, S.; Pühler, A.; Scherer, P.; Klocke, M.; Schlüter, A. Metagenome, metatranscriptome, and metaproteome approaches unraveled compositions and functional relationships of microbial communities residing in biogas plants. Appl. Microbiol. Biotechnol. 2018, 102, 5045–5063. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Sundberg, C.; Al-Soud, W.A.; Larsson, M.; Alm, E.; Yekta, S.S.; Svensson, B.H.; Sørensen, S.J.; Karlsson, A. 454 pyrosequencing analyses of bacterial and archaeal richness in 21 full-scale biogas digesters. FEMS Microbiol. Ecol. 2013, 85, 612–626. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. De Vrieze, J.; Saunders, A.M.; He, Y.; Fang, J.; Nielsen, P.H.; Verstraete, W.; Boon, N. Ammonia and temperature determine potential clustering in the anaerobic digestion microbiome. Water Res. 2015, 75, 312–323. [Google Scholar] [CrossRef] [PubMed]
  23. Westerholm, M.; Isaksson, S.; Lindsjö, O.K.; Schnürer, A. Microbial community adaptability to altered temperature conditions determines the potential for process optimisation in biogas production. Appl. Energy 2018, 226, 838–848. [Google Scholar] [CrossRef]
  24. Mei, R.; Nobu, M.K.; Narihiro, T.; Kuroda, K.; Sierra, J.M.; Wu, Z.; Ye, L.; Lee, P.K.; Lee, P.H.; Van Lier, J.B.; et al. Operation-driven heterogeneity and overlooked feed-associated populations in global anaerobic digester microbiome. Water Res. 2017, 124, 77–84. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Klang, J.; Szewzyk, U.; Bock, D.; Theuerl, S. Nexus between the microbial diversity level and the stress tolerance within the biogas process. Anaerobe 2019, 56, 8–16. [Google Scholar] [CrossRef]
  26. Jousset, A.; Bienhold, C.; Chatzinotas, A.; Gallien, L.; Gobet, A.; Kurm, V.; Küsel, K.; Rillig, M.C.; Rivett, D.W.; Salles, J.F.; et al. Where less may be more: How the rare biosphere pulls ecosystems strings. ISME J. 2017, 11, 853–862. [Google Scholar] [CrossRef]
  27. Werner, J.J.; Knights, D.; Garcia, M.L.; Scalfone, N.B.; Smith, S.; Yarasheski, K.; Cummings, T.A.; Beers, A.R.; Knight, R.; Angenent, L.T. Bacterial community structures are unique and resilient in full-scale bioenergy systems. Proc. Natl. Acad. Sci. USA 2011, 108, 4158–4163. [Google Scholar] [CrossRef] [Green Version]
  28. Ziels, R.M.; Svensson, B.H.; Sundberg, C.; Larsson, M.; Karlsson, A.; Yekta, S.S. Microbial rRNA gene expression and co-occurrence profiles associate with biokinetics and elemental composition in full-scale anaerobic digesters. Microb. Biotechnol. 2018, 11, 694–709. [Google Scholar] [CrossRef]
  29. Theuerl, S.; Klang, J.; Prochnow, A. Microbiome Diversity and Community-Level Change Points within Manure-based small Biogas Plants. Microorganisms 2020, 8, 1169. [Google Scholar] [CrossRef]
  30. Westerholm, M.; Moestedt, J.; Schnürer, A. Biogas production through syntrophic acetate oxidation and deliberate operating strategies for improved digester performance. Appl. Energy 2016, 179, 124–135. [Google Scholar] [CrossRef] [Green Version]
  31. Theuerl, S.; Klang, J.; Heiermann, M.; De Vrieze, J. Marker microbiome clusters are determined by operational parameters and specific key taxa combinations in anaerobic digestion. Bioresour. Technol. 2018, 263, 128–135. [Google Scholar] [CrossRef]
  32. Lewin, G.R.; Carlos, C.; Chevrette, M.G.; Horn, H.A.; McDonald, B.R.; Stankey, R.J.; Fox, B.G.; Currie, C.R. Evolution and ecology of Actinobacteria and their bioenergy applications. Annu. Rev. Microbiol. 2016, 70, 235–254. [Google Scholar] [CrossRef] [Green Version]
  33. Wang, C.; Dong, D.; Wang, H.; Müller, K.; Qin, Y.; Wang, H.; Wu, W. Metagenomic analysis of microbial consortia enriched from compost: New insights into the role of Actinobacteria in lignocellulose decomposition. Biotechnol. Biofuels 2016, 9, 1–17. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Maus, I.; Rumming, M.; Bergmann, I.; Heeg, K.; Pohl, M.; Nettmann, E.; Jaenicke, S.; Blom, J.; Pühler, A.; Schlüter, A.; et al. Characterization of Bathyarchaeota genomes assembled from metagenomes of biofilms residing in mesophilic and thermophilic biogas reactors. Biotechnol. Biofuels 2018, 11, 1–13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Zhou, Z.; Pan, J.; Wang, F.; Gu, J.D.; Li, M. Bathyarchaeota: Globally distributed metabolic generalists in anoxic environments. FEMS Microbiol. Rev. 2018, 42, 639–655. [Google Scholar] [CrossRef] [PubMed]
  36. Theuerl, S.; Kohrs, F.; Benndorf, D.; Maus, I.; Wibberg, D.; Schlüter, A.; Kausmann, R.; Heiermann, M.; Rapp, E.; Reichl, U.; et al. Community shifts in a well-operating agricultural biogas plant: How process variations are handled by the microbiome. Appl. Microbiol. Biotechnol. 2015, 99, 7791–7803. [Google Scholar] [CrossRef] [PubMed]
  37. Liebetrau, J.; Pfeiffer, D.; Thrän, D. Collection of Methods for Biogas—Methods to Determine Parameters for Analysis Purposes and Parameters That Describe Processes in the Biogas Sector; DBFZ Deutsches Biomasseforschungszentrum Gemeinnützige GmbH: Leipzig, Germany; Fischer Druck: Peine, Germany, 2016. [Google Scholar]
  38. Hansen, K.H.; Angelidaki, I.; Ahring, B.K. Anaerobic digestion of swine manure: Inhibition by ammonia. Water Res. 1998, 32, 5–12. [Google Scholar] [CrossRef]
  39. Takahashi, S.; Tomita, J.; Nishioka, K.; Hisada, T.; Nishijima, M. Development of a prokaryotic universal primer for simultaneous analysis of Bacteria and Archaea using next-generation sequencing. PLoS ONE 2014, 9, e105592. [Google Scholar] [CrossRef] [Green Version]
  40. Magoč, T.; Salzberg, S.L. FLASH: Fast length adjustment of short reads to improve genome assemblies. Bioinformatics 2011, 27, 2957–2963. [Google Scholar] [CrossRef]
  41. Martin, M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet. J. 2011, 17, 10–12. [Google Scholar] [CrossRef]
  42. Joshi Nikhil, F.J. Sickle: A Sliding-Window, Adaptive, Quality-Based Trimming Tool for FastQ Files. 2011. Available online: https://github.com/najoshi/sickle (accessed on 5 July 2021).
  43. Caporaso, J.G.; Kuczynski, J.; Stombaugh, J.; Bittinger, K.; Bushman, F.D.; Costello, E.K.; Fierer, N.; Pena, A.G.; Goodrich, J.K.; Gordon, J.I.; et al. QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 2010, 7, 335–336. [Google Scholar] [CrossRef] [Green Version]
  44. Clarke, K.R. Non-parametric multivariate analyses of changes in community structure. Aust. J. Ecol. 1993, 18, 117–143. [Google Scholar] [CrossRef]
  45. Ramette, A. Multivariate analyses in microbial ecology. FEMS Microbiol. Ecol. 2007, 62, 142–160. [Google Scholar] [CrossRef] [Green Version]
  46. Paliy, O.; Shankar, V. Application of multivariate statistical techniques in microbial ecology. Mol. Ecol. 2016, 25, 1032–1057. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2020. [Google Scholar]
  48. Oksanen, J.; Blanchet, F.G.; Friendly, M.; Kindt, R.; Legendre, P.; McGlinn, D.; Minchin, P.R.; O’Hara, R.B.; Simpson, G.L.; Solymos, P.; et al. vegan: Community Ecology Package; R Package Version 2.5-7; R Foundation for Statistical Computing: Vienna, Austria, 2020. [Google Scholar]
  49. Warnes, G.R.; Bolker, B.; Bonebakker, L.; Gentleman, R.; Huber, W.; Liaw, A.; Lumley, T.; Maechler, M.; Magnusson, A.; Moeller, S.; et al. gplots: Various R Programming Tools for Plotting Data; R Package Version 3.1.1; R Foundation for Statistical Computing: Vienna, Austria, 2020. [Google Scholar]
  50. McQuitty, L.L. Similarity analysis by reciprocal pairs for discrete and continuous data. Educ. Psychol. Meas. 1966, 26, 825–831. [Google Scholar] [CrossRef]
  51. Dufrêne, M.; Legendre, P. Species assemblages and indicator species: The need for a flexible asymmetrical approach. Ecol. Monogr. 1997, 67, 345–366. [Google Scholar] [CrossRef]
  52. De Cáceres, M.; Legendre, P. Associations between species and groups of sites: Indices and statistical inference. Ecology 2009, 90, 3566–3574. [Google Scholar] [CrossRef] [PubMed]
  53. De Cáceres, M.; Sol, D.; Lapiedra, O.; Legendre, P. A framework for estimating niche metrics using the resemblance between qualitative resources. Oikos 2011, 120, 1341–1350. [Google Scholar] [CrossRef]
  54. Kim, J.; Lee, C. Response of a continuous anaerobic digester to temperature transitions: A critical range for restructuring the microbial community structure and function. Water Res. 2016, 89, 241–251. [Google Scholar] [CrossRef]
  55. Maus, I.; Kim, Y.S.; Wibberg, D.; Stolze, Y.; Off, S.; Antonczyk, S.; Pühler, A.; Scherer, P.; Schlüter, A. Biphasic study to characterize agricultural biogas plants by high-throughput 16S rRNA gene amplicon sequencing and microscopic analysis. J. Microbiol. Biotechnol. 2017, 27, 321–334. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Heyer, R.; Benndorf, D.; Kohrs, F.; De Vrieze, J.; Boon, N.; Hoffmann, M.; Rapp, E.; Schlüter, A.; Sczyrba, A.; Reichl, U. Proteotyping of biogas plant microbiomes separates biogas plants according to process temperature and reactor type. Biotechnol. Biofuels 2016, 9, 1–16. [Google Scholar] [CrossRef] [Green Version]
  57. Hupfauf, S.; Plattner, P.; Wagner, A.O.; Kaufmann, R.; Insam, H.; Podmirseg, S.M. Temperature shapes the microbiota in anaerobic digestion and drives efficiency to a maximum at 45 C. Bioresour. Technol. 2018, 269, 309–318. [Google Scholar] [CrossRef] [PubMed]
  58. Drosg, B. Process Monitoring in Biogas Plants; IEA Bioenergy Paris: Paris, France, 2013. [Google Scholar]
  59. Amha, Y.M.; Anwar, M.Z.; Brower, A.; Jacobsen, C.S.; Stadler, L.B.; Webster, T.M.; Smith, A.L. Inhibition of anaerobic digestion processes: Applications of molecular tools. Bioresour. Technol. 2018, 247, 999–1014. [Google Scholar] [CrossRef] [PubMed]
  60. Pap, B.; Györkei, Á.; Boboescu, I.Z.; Nagy, I.K.; Bíró, T.; Kondorosi, É.; Maróti, G. Temperature-dependent transformation of biogas-producing microbial communities points to the increased importance of hydrogenotrophic methanogenesis under thermophilic operation. Bioresour. Technol. 2015, 177, 375–380. [Google Scholar] [CrossRef] [PubMed]
  61. Regueiro, L.; Carballa, M.; Lema, J.M. Outlining microbial community dynamics during temperature drop and subsequent recovery period in anaerobic co-digestion systems. J. Biotechnol. 2014, 192, 179–186. [Google Scholar] [CrossRef] [PubMed]
  62. Sun, L.; Toyonaga, M.; Ohashi, A.; Tourlousse, D.M.; Matsuura, N.; Meng, X.Y.; Tamaki, H.; Hanada, S.; Cruz, R.; Yamaguchi, T.; et al. Lentimicrobium saccharophilum gen. nov., sp. nov., a strictly anaerobic bacterium representing a new family in the phylum Bacteroidetes, and proposal of Lentimicrobiaceae fam. nov. Int. J. Syst. Evol. Microbiol. 2016, 66, 2635–2642. [Google Scholar] [CrossRef] [Green Version]
  63. Chen, S.; Dong, X. Proteiniphilum acetatigenes gen. nov., sp. nov., from a UASB reactor treating brewery wastewater. Int. J. Syst. Evol. Microbiol. 2005, 55, 2257–2261. [Google Scholar] [CrossRef] [Green Version]
  64. Hahnke, S.; Langer, T.; Koeck, D.E.; Klocke, M. Description of Proteiniphilum saccharofermentans sp. nov., Petrimonas mucosa sp. nov. and Fermentimonas caenicola gen. nov., sp. nov., isolated from mesophilic laboratory-scale biogas reactors, and emended description of the genus Proteiniphilum. Int. J. Syst. Evol. Microbiol. 2016, 66, 1466–1475. [Google Scholar] [CrossRef] [PubMed]
  65. Johnson, J.S.; Spakowicz, D.J.; Hong, B.Y.; Petersen, L.M.; Demkowicz, P.; Chen, L.; Leopold, S.R.; Hanson, B.M.; Agresta, H.O.; Gerstein, M.; et al. Evaluation of 16S rRNA gene sequencing for species and strain-level microbiome analysis. Nat. Commun. 2019, 10, 1–11. [Google Scholar] [CrossRef] [Green Version]
  66. Kim, M.; Morrison, M.; Yu, Z. Evaluation of different partial 16S rRNA gene sequence regions for phylogenetic analysis of microbiomes. J. Microbiol. Methods 2011, 84, 81–87. [Google Scholar] [CrossRef] [PubMed]
  67. Hülsemann, B.; Zhou, L.; Merkle, W.; Hassa, J.; Müller, J.; Oechsner, H. Biomethane potential test: Influence of inoculum and the digestion system. Appl. Sci. 2020, 10, 2589. [Google Scholar] [CrossRef] [Green Version]
  68. Shade, A.; Jones, S.E.; Caporaso, J.G.; Handelsman, J.; Knight, R.; Fierer, N.; Gilbert, J.A. Conditionally rare taxa disproportionately contribute to temporal changes in microbial diversity. MBio 2014, 5, e01371-14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. De Vrieze, J.; Verstraete, W. Perspectives for microbial community composition in anaerobic digestion: From abundance and activity to connectivity. Environ. Microbiol. 2016, 18, 2797–2809. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Range of the operational and chemical process parameters from the 67 analyzed anaerobic digesters. The different units of the paramaters are indicated on the x-axis. Supplied feedstocks were normalized to 100% for each digester. The acid concentrations were summarized as volatile fatty acids (VFA) given as acetic acid equivalents. OLR = organic loading rate, HRT = hydraulic retention time.
Figure 1. Range of the operational and chemical process parameters from the 67 analyzed anaerobic digesters. The different units of the paramaters are indicated on the x-axis. Supplied feedstocks were normalized to 100% for each digester. The acid concentrations were summarized as volatile fatty acids (VFA) given as acetic acid equivalents. OLR = organic loading rate, HRT = hydraulic retention time.
Microorganisms 09 01457 g001
Figure 2. Taxonomic profiles on order level of the 67 analyzed anaerobic digesters from 49 agricultural, full-scale biogas plants as deduced from 16S rRNA gene amplicon data. Orders with a maximal relative abundance less than 1% were summarized. Taxa which were not taxonomically assigned at order level were summarized as “unclassified”.
Figure 2. Taxonomic profiles on order level of the 67 analyzed anaerobic digesters from 49 agricultural, full-scale biogas plants as deduced from 16S rRNA gene amplicon data. Orders with a maximal relative abundance less than 1% were summarized. Taxa which were not taxonomically assigned at order level were summarized as “unclassified”.
Microorganisms 09 01457 g002
Figure 3. Tree diagram based on a McQuitty (dis)similarities matrix calculated from the taxonomic profiles of the analyzed 67 anaerobic digesters from 49 agricultural biogas plants.
Figure 3. Tree diagram based on a McQuitty (dis)similarities matrix calculated from the taxonomic profiles of the analyzed 67 anaerobic digesters from 49 agricultural biogas plants.
Microorganisms 09 01457 g003
Figure 4. Non-metric multidimensional scaling (NMDS) based on the taxonomic diversity of the 67 analyzed digester samples from 49 agricultural, full-scale biogas plants. The color code is deduced from the tree diagram in Figure 3: Cluster I = green, cluster IIa = red, cluster IIb = blue, outlier = gray. The given environmental vectors symbolize the amount of maize silage (R 2 = 0.30, p = 0.001), whole crop rye silage (R 2 = 0.30, p = 0.001), liquid cattle manure (R 2 = 0.32, p = 0.001), the process temperature (R 2 = 0.74, p = 0.001), total solids (TS) (R 2 = 0.33, p = 0.001), the total ammonium nitrogen (TAN) (R 2 = 0.40, p = 0.001), the total Kjeldahl nitrogen (TKN) (R 2 = 0.31, p = 0.001) and the free ammonia nitrogen concentration (NH 3 ) (R 2 = 0.51, p = 0.001).
Figure 4. Non-metric multidimensional scaling (NMDS) based on the taxonomic diversity of the 67 analyzed digester samples from 49 agricultural, full-scale biogas plants. The color code is deduced from the tree diagram in Figure 3: Cluster I = green, cluster IIa = red, cluster IIb = blue, outlier = gray. The given environmental vectors symbolize the amount of maize silage (R 2 = 0.30, p = 0.001), whole crop rye silage (R 2 = 0.30, p = 0.001), liquid cattle manure (R 2 = 0.32, p = 0.001), the process temperature (R 2 = 0.74, p = 0.001), total solids (TS) (R 2 = 0.33, p = 0.001), the total ammonium nitrogen (TAN) (R 2 = 0.40, p = 0.001), the total Kjeldahl nitrogen (TKN) (R 2 = 0.31, p = 0.001) and the free ammonia nitrogen concentration (NH 3 ) (R 2 = 0.51, p = 0.001).
Microorganisms 09 01457 g004
Figure 5. Exemplary distribution of the 21 most abundant indicative taxa and their relative abundances within the defined clusters IIa, IIb, I based on the taxonomic profiles. Indicator taxa of the respective clusters are boxed in red for cluster IIa, blue for cluster IIb and green for cluster I. Statistically significant differences were calculated using the “indicspecies” package of R with a point biserial correlation coefficient as statistic value [52,53].
Figure 5. Exemplary distribution of the 21 most abundant indicative taxa and their relative abundances within the defined clusters IIa, IIb, I based on the taxonomic profiles. Indicator taxa of the respective clusters are boxed in red for cluster IIa, blue for cluster IIb and green for cluster I. Statistically significant differences were calculated using the “indicspecies” package of R with a point biserial correlation coefficient as statistic value [52,53].
Microorganisms 09 01457 g005
Figure 6. Non-metric multidimensional scaling (NMDS) (A) based on the identified indicative taxa as deduced from indicator species analyses (ISA) and (B) the generally occurring taxa of the 67 analyzed digesters from 49 agricultural, full-scale biogas plants. The color code corresponds to the clusters in Figure 4: I = green, IIa = red, IIb = blue, outlier = gray.
Figure 6. Non-metric multidimensional scaling (NMDS) (A) based on the identified indicative taxa as deduced from indicator species analyses (ISA) and (B) the generally occurring taxa of the 67 analyzed digesters from 49 agricultural, full-scale biogas plants. The color code corresponds to the clusters in Figure 4: I = green, IIa = red, IIb = blue, outlier = gray.
Microorganisms 09 01457 g006
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hassa, J.; Klang, J.; Benndorf, D.; Pohl, M.; Hülsemann, B.; Mächtig, T.; Effenberger, M.; Pühler, A.; Schlüter, A.; Theuerl, S. Indicative Marker Microbiome Structures Deduced from the Taxonomic Inventory of 67 Full-Scale Anaerobic Digesters of 49 Agricultural Biogas Plants. Microorganisms 2021, 9, 1457. https://doi.org/10.3390/microorganisms9071457

AMA Style

Hassa J, Klang J, Benndorf D, Pohl M, Hülsemann B, Mächtig T, Effenberger M, Pühler A, Schlüter A, Theuerl S. Indicative Marker Microbiome Structures Deduced from the Taxonomic Inventory of 67 Full-Scale Anaerobic Digesters of 49 Agricultural Biogas Plants. Microorganisms. 2021; 9(7):1457. https://doi.org/10.3390/microorganisms9071457

Chicago/Turabian Style

Hassa, Julia, Johanna Klang, Dirk Benndorf, Marcel Pohl, Benedikt Hülsemann, Torsten Mächtig, Mathias Effenberger, Alfred Pühler, Andreas Schlüter, and Susanne Theuerl. 2021. "Indicative Marker Microbiome Structures Deduced from the Taxonomic Inventory of 67 Full-Scale Anaerobic Digesters of 49 Agricultural Biogas Plants" Microorganisms 9, no. 7: 1457. https://doi.org/10.3390/microorganisms9071457

APA Style

Hassa, J., Klang, J., Benndorf, D., Pohl, M., Hülsemann, B., Mächtig, T., Effenberger, M., Pühler, A., Schlüter, A., & Theuerl, S. (2021). Indicative Marker Microbiome Structures Deduced from the Taxonomic Inventory of 67 Full-Scale Anaerobic Digesters of 49 Agricultural Biogas Plants. Microorganisms, 9(7), 1457. https://doi.org/10.3390/microorganisms9071457

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