High level of intergenera gene exchange shapes the evolution of haloarchaea in an isolated Antarctic lake.

Deep Lake in Antarctica is a globally isolated, hypersaline system that remains liquid at temperatures down to -20 °C. By analyzing metagenome data and genomes of four isolates we assessed genome variation and patterns of gene exchange to learn how the lake community evolved. The lake is completely and uniformly dominated by haloarchaea, comprising a hierarchically structured, low-complexity community that differs greatly to temperate and tropical hypersaline environments. The four Deep Lake isolates represent distinct genera (∼85% 16S rRNA gene similarity and ∼73% genome average nucleotide identity) with genomic characteristics indicative of niche adaptation, and collectively account for ∼72% of the cellular community. Network analysis revealed a remarkable level of intergenera gene exchange, including the sharing of long contiguous regions (up to 35 kb) of high identity (∼100%). Although the genomes of closely related Halobacterium, Haloquadratum, and Haloarcula (>90% average nucleotide identity) shared regions of high identity between species or strains, the four Deep Lake isolates were the only distantly related haloarchaea to share long high-identity regions. Moreover, the Deep Lake high-identity regions did not match to any other hypersaline environment metagenome data. The most abundant species, tADL, appears to play a central role in the exchange of insertion sequences, but not the exchange of high-identity regions. The genomic characteristics of the four haloarchaea are consistent with a lake ecosystem that sustains a high level of intergenera gene exchange while selecting for ecotypes that maintain sympatric speciation. The peculiarities of this polar system restrict which species can grow and provide a tempo and mode for accentuating gene exchange.

Phage Genome Diversity in a Biogas-Producing Microbiome Analyzed by Illumina and Nanopore GridION Sequencing.

The microbial biogas network is complex and intertwined, and therefore relatively stable in its overall functionality. However, if key functional groups of microorganisms are affected by biotic or abiotic factors, the entire efficacy may be impaired. Bacteriophages are hypothesized to alter the steering process of the microbial network. In this study, an enriched fraction of virus-like particles was extracted from a mesophilic biogas reactor and sequenced on the Illumina MiSeq and Nanopore GridION sequencing platforms. Metagenome data analysis resulted in identifying 375 metagenome-assembled viral genomes (MAVGs). Two-thirds of the classified sequences were only assigned to the superkingdom Viruses and the remaining third to the family Siphoviridae, followed by Myoviridae, Podoviridae, Tectiviridae, and Inoviridae. The metavirome showed a close relationship to the phage genomes that infect members of the classes Clostridia and Bacilli. Using publicly available biogas metagenomic data, a fragment recruitment approach showed the widespread distribution of the MAVGs studied in other biogas microbiomes. In particular, phage sequences from mesophilic microbiomes were highly similar to the phage sequences of this study. Accordingly, the virus particle enrichment approach and metavirome sequencing provided additional genome sequence information for novel virome members, thus expanding the current knowledge of viral genetic diversity in biogas reactors.

Genome Analyses and Genome-Centered Metatranscriptomics of Methanothermobacter wolfeii Strain SIV6, Isolated from a Thermophilic Production-Scale Biogas Fermenter.

In the thermophilic biogas-producing microbial community, the genus Methanothermobacter was previously described to be frequently abundant. The aim of this study was to establish and analyze the genome sequence of the archaeal strain Methanothermobacter wolfeii SIV6 originating from a thermophilic industrial-scale biogas fermenter and compare it to related reference genomes. The circular chromosome has a size of 1,686,891 bases, featuring a GC content of 48.89%. Comparative analyses considering three completely sequenced Methanothermobacter strains revealed a core genome of 1494 coding sequences and 16 strain specific genes for M. wolfeii SIV6, which include glycosyltransferases and CRISPR/cas associated genes. Moreover, M. wolfeii SIV6 harbors all genes for the hydrogenotrophic methanogenesis pathway and genome-centered metatranscriptomics indicates the high metabolic activity of this strain, with 25.18% of all transcripts per million (TPM) belong to the hydrogenotrophic methanogenesis pathway and 18.02% of these TPM exclusively belonging to the mcr operon. This operon encodes the different subunits of the enzyme methyl-coenzyme M reductase (EC: 2.8.4.1), which catalyzes the final and rate-limiting step during methanogenesis. Finally, fragment recruitment of metagenomic reads from the thermophilic biogas fermenter on the SIV6 genome showed that the strain is abundant (1.2%) within the indigenous microbial community. Detailed analysis of the archaeal isolate M. wolfeii SIV6 indicates its role and function within the microbial community of the thermophilic biogas fermenter, towards a better understanding of the biogas production process and a microbial-based management of this complex process.

Genomics and prevalence of bacterial and archaeal isolates from biogas-producing microbiomes.

BACKGROUND: To elucidate biogas microbial communities and processes, the application of high-throughput DNA analysis approaches is becoming increasingly important. Unfortunately, generated data can only partialy be interpreted rudimentary since databases lack reference sequences. RESULTS: Novel cellulolytic, hydrolytic, and acidogenic/acetogenic Bacteria as well as methanogenic Archaea originating from different anaerobic digestion communities were analyzed on the genomic level to assess their role in biomass decomposition and biogas production. Some of the analyzed bacterial strains were recently described as new species and even genera, namely Herbinix hemicellulosilytica T3/55T, Herbinix luporum SD1DT, Clostridium bornimense M2/40T, Proteiniphilum saccharofermentans M3/6T, Fermentimonas caenicola ING2-E5BT, and Petrimonas mucosa ING2-E5AT. High-throughput genome sequencing of 22 anaerobic digestion isolates enabled functional genome interpretation, metabolic reconstruction, and prediction of microbial traits regarding their abilities to utilize complex bio-polymers and to perform specific fermentation pathways. To determine the prevalence of the isolates included in this study in different biogas systems, corresponding metagenome fragment mappings were done. Methanoculleus bourgensis was found to be abundant in three mesophilic biogas plants studied and slightly less abundant in a thermophilic biogas plant, whereas Defluviitoga tunisiensis was only prominent in the thermophilic system. Moreover, several of the analyzed species were clearly detectable in the mesophilic biogas plants, but appeared to be only moderately abundant. Among the species for which genome sequence information was publicly available prior to this study, only the species Amphibacillus xylanus, Clostridium clariflavum, and Lactobacillus acidophilus are of importance for the biogas microbiomes analyzed, but did not reach the level of abundance as determined for M. bourgensis and D. tunisiensis. CONCLUSIONS: Isolation of key anaerobic digestion microorganisms and their functional interpretation was achieved by application of elaborated cultivation techniques and subsequent genome analyses. New isolates and their genome information extend the repository covering anaerobic digestion community members.

Unraveling the microbiome of a thermophilic biogas plant by metagenome and metatranscriptome analysis complemented by characterization of bacterial and archaeal isolates.

BACKGROUND: One of the most promising technologies to sustainably produce energy and to mitigate greenhouse gas emissions from combustion of fossil energy carriers is the anaerobic digestion and biomethanation of organic raw material and waste towards biogas by highly diverse microbial consortia. In this context, the microbial systems ecology of thermophilic industrial-scale biogas plants is poorly understood. RESULTS: The microbial community structure of an exemplary thermophilic biogas plant was analyzed by a comprehensive approach comprising the analysis of the microbial metagenome and metatranscriptome complemented by the cultivation of hydrolytic and acido-/acetogenic Bacteria as well as methanogenic Archaea. Analysis of metagenome-derived 16S rRNA gene sequences revealed that the bacterial genera Defluviitoga (5.5 %), Halocella (3.5 %), Clostridium sensu stricto (1.9 %), Clostridium cluster III (1.5 %), and Tepidimicrobium (0.7 %) were most abundant. Among the Archaea, Methanoculleus (2.8 %) and Methanothermobacter (0.8 %) were predominant. As revealed by a metatranscriptomic 16S rRNA analysis, Defluviitoga (9.2 %), Clostridium cluster III (4.8 %), and Tepidanaerobacter (1.1 %) as well as Methanoculleus (5.7 %) mainly contributed to these sequence tags indicating their metabolic activity, whereas Hallocella (1.8 %), Tepidimicrobium (0.5 %), and Methanothermobacter (<0.1 %) were transcriptionally less active. By applying 11 different cultivation strategies, 52 taxonomically different microbial isolates representing the classes Clostridia, Bacilli, Thermotogae, Methanomicrobia and Methanobacteria were obtained. Genome analyses of isolates support the finding that, besides Clostridium thermocellum and Clostridium stercorarium, Defluviitoga tunisiensis participated in the hydrolysis of hemicellulose producing ethanol, acetate, and H2/CO2. The latter three metabolites are substrates for hydrogentrophic and acetoclastic archaeal methanogenesis. CONCLUSIONS: Obtained results showed that high abundance of microorganisms as deduced from metagenome analysis does not necessarily indicate high transcriptional or metabolic activity, and vice versa. Additionally, it appeared that the microbiome of the investigated thermophilic biogas plant comprised a huge number of up to now unknown and insufficiently characterized species.

Comparative metagenomics of biogas-producing microbial communities from production-scale biogas plants operating under wet or dry fermentation conditions.

BACKGROUND: Decomposition of biomass for biogas production can be practiced under wet and dry fermentation conditions. In contrast to the dry fermentation technology, wet fermentation is characterized by a high liquid content and a relatively low total solid content. In this study, the composition and functional potential of a biogas-producing microbial community in an agricultural biogas reactor operating under wet fermentation conditions was analyzed by a metagenomic approach applying 454-pyrosequencing. The obtained metagenomic dataset and corresponding 16S rRNA gene amplicon sequences were compared to the previously sequenced comparable metagenome from a dry fermentation process, meeting explicitly identical boundary conditions regarding sample and community DNA preparation, sequencing technology, processing of sequence reads and data analyses by bioinformatics tools. RESULTS: High-throughput metagenome sequencing of community DNA from the wet fermentation process applying the pyrosequencing approach resulted in 1,532,780 reads, with an average read length of 397 bp, accounting for approximately 594 million bases of sequence information in total. Taxonomic comparison of the communities from wet and dry fermentation revealed similar microbial profiles with Bacteria being the predominant superkingdom, while the superkingdom Archaea was less abundant. In both biogas plants, the bacterial phyla Firmicutes, Bacteroidetes, Spirochaetes and Proteobacteria were identified with descending frequencies. Within the archaeal superkingdom, the phylum Euryarchaeota was most abundant with the dominant class Methanomicrobia. Functional profiles of the communities revealed that environmental gene tags representing methanogenesis enzymes were present in both biogas plants in comparable frequencies. 16S rRNA gene amplicon high-throughput sequencing disclosed differences in the sub-communities comprising methanogenic Archaea between both processes. Fragment recruitments of metagenomic reads to the reference genome of the archaeon Methanoculleus bourgensis MS2(T) revealed that dominant methanogens within the dry fermentation process were highly related to the reference. CONCLUSIONS: Although process parameters, substrates and technology differ between the wet and dry biogas fermentations analyzed in this study, community profiles are very similar at least at higher taxonomic ranks, illustrating that core community taxa perform key functions in biomass decomposition and methane synthesis. Regarding methanogenesis, Archaea highly related to the type strain M. bourgensis MS2(T) dominate the dry fermentation process, suggesting the adaptation of members belonging to this species to specific fermentation process parameters.

Effect of the strain Bacillus amyloliquefaciens FZB42 on the microbial community in the rhizosphere of lettuce under field conditions analyzed by whole metagenome sequencing.

Application of the plant associated bacterium Bacillus amyloliquefaciens FZB42 on lettuce (Lactuca sativa) confirmed its capability to promote plant growth and health by reducing disease severity (DS) caused by the phytopathogenic fungus Rhizoctonia solani. Therefore this strain is commercially applied as an eco-friendly plant protective agent. It is able to produce cyclic lipopeptides (CLP) and polyketides featuring antifungal and antibacterial properties. Production of these secondary metabolites led to the question of a possible impact of strain FZB42 on the composition of microbial rhizosphere communities after its application. Rating of DS and lettuce growth during a field trial confirmed the positive impact of strain FZB42 on the health of the host plant. To verify B. amyloliquefaciens as an environmentally compatible plant protective agent, its effect on the indigenous rhizosphere community was analyzed by metagenome sequencing. Rhizosphere microbial communities of lettuce treated with B. amyloliquefaciens FZB42 and non-treated plants were profiled by high-throughput metagenome sequencing of whole community DNA. Fragment recruitments of metagenome sequence reads on the genome sequence of B. amyloliquefaciens FZB42 proved the presence of the strain in the rhizosphere over 5 weeks of the field trial. Comparison of taxonomic community profiles only revealed marginal changes after application of strain FZB42. The orders Burkholderiales, Actinomycetales and Rhizobiales were most abundant in all samples. Depending on plant age a general shift within the composition of the microbial communities that was independent of the application of strain FZB42 was observed. In addition to the taxonomic profiling, functional analysis of annotated sequences revealed no major differences between samples regarding application of the inoculant strain.

Metagenome analyses reveal the influence of the inoculant Lactobacillus buchneri CD034 on the microbial community involved in grass ensiling.

Silage is green fodder conserved by lactic acid fermentation performed by epiphytic lactic acid bacteria under anaerobic conditions. To improve the ensiling process and the quality of the resulting silage, starter cultures are added to the fresh forage. A detailed analysis of the microbial community playing a role in grass ensiling has been carried out by high throughput sequencing technologies. Moreover, the influence of the inoculant Lactobacillus buchneri CD034 on the microbial community composition was studied. For this purpose, grass was ensiled untreated or inoculated with L. buchneri CD034. The fresh forage as well as silages after 14 and 58 days of fermentation were characterized physico-chemically. Characteristic silage conditions such as increased titers of lactic acid bacteria and higher concentrations of acetic acid were observed in the inoculated silage in comparison to the untreated samples. Taxonomic community profiles deduced from 16S rDNA amplicon sequences indicated that the relative abundance of Lactococci diminished in the course of fermentations and that the proportion of bacteria belonging to the phyla Proteobacteria and Bacteroidetes increased during the fermentation of untreated silage. In the inoculated silage, members of these phyla were repressed due to an increased abundance of Lactobacilli. In addition, metagenome analyses of silage samples confirmed taxonomic profiles based on 16S rDNA amplicons. Moreover, Lactobacillus plantarum, Lactobacillus brevis and Lactococcus lactis were found to be dominant species within silages as analyzed by means of fragment recruitments of metagenomic sequence reads on complete reference genome sequences. Fragment recruitments also provided clear evidence for the competitiveness of the inoculant strain L. buchneri CD034 during the fermentation of the inoculated silage. The inoculation strain was able to outcompete other community members and also affected physico-chemical characteristics of the silage.