Bernd Wemheuer

Tax4Fun: predicting functional profiles from metagenomic 16S rRNA data

Motivation: The characterization of phylogenetic and functional diversity is a key element in the analysis of microbial communities. Amplicon-based sequencing of marker genes, such as 16S rRNA, is a powerful tool for assessing and comparing the structure of microbial communities at a high phylogenetic resolution. Because 16S rRNA sequencing is more cost-effective than whole metagenome shotgun sequencing, marker gene analysis is frequently used for broad studies that involve a large number of different samples. However, in comparison to shotgun sequencing approaches, insights into the functional capabilities of the community get lost when restricting the analysis to taxonomic assignment of 16S rRNA data. Results: Tax4Fun is a software package that predicts the functional capabilities of microbial communities based on 16S rRNA datasets. We evaluated Tax4Fun on a range of paired metagenome/16S rRNA datasets to assess its performance. Our results indicate that Tax4Fun provides a good approximation to functional profiles obtained from metagenomic shotgun sequencing approaches. Availability and implementation: Tax4Fun is an open-source R package and applicable to output as obtained from the SILVAngs web server or the application of QIIME with a SILVA database extension. Tax4Fun is freely available for download at http://tax4fun.gobics.de/. Contact: kasshau@gwdg.de Supplementary information: Supplementary data are available at Bioinformatics online.

Adaptation of an abundant Roseobacter RCA organism to pelagic systems revealed by genomic and transcriptomic analyses.

The RCA (Roseobacter clade affiliated) cluster, with an internal 16S rRNA gene sequence similarity of >98%, is the largest cluster of the marine Roseobacter clade and most abundant in temperate to (sub)polar oceans, constituting up to 35% of total bacterioplankton. The genome analysis of the first described species of the RCA cluster, Planktomarina temperata RCA23, revealed that this phylogenetic lineage is deeply branching within the Roseobacter clade. It shares not >65.7% of homologous genes with any other organism of this clade. The genome is the smallest of all closed genomes of the Roseobacter clade, exhibits various features of genome streamlining and encompasses genes for aerobic anoxygenic photosynthesis (AAP) and CO oxidation. In order to assess the biogeochemical significance of the RCA cluster we investigated a phytoplankton spring bloom in the North Sea. This cluster constituted 5.1% of the total, but 10–31% (mean 18.5%) of the active bacterioplankton. A metatranscriptomic analysis showed that the genome of P. temperata RCA23 was transcribed to 94% in the bloom with some variations during day and night. The genome of P. temperata RCA23 was also retrieved to 84% from metagenomic data sets from a Norwegian fjord and to 82% from stations of the Global Ocean Sampling expedition in the northwestern Atlantic. In this region, up to 6.5% of the total reads mapped on the genome of P. temperata RCA23. This abundant taxon appears to be a major player in ocean biogeochemistry.

Metagenome Survey of a Multispecies and Alga-Associated Biofilm Revealed Key Elements of Bacterial-Algal Interactions in Photobioreactors

ABSTRACT Photobioreactors (PBRs) are very attractive for sunlight-driven production of biofuels and capturing of anthropogenic CO2. One major problem associated with PBRs however, is that the bacteria usually associated with microalgae in nonaxenic cultures can lead to biofouling and thereby affect algal productivity. Here, we report on a phylogenetic, metagenome, and functional analysis of a mixed-species bacterial biofilm associated with the microalgae Chlorella vulgaris and Scenedesmus obliquus in a PBR. The biofilm diversity and population dynamics were examined through 16S rRNA phylogeny. Overall, the diversity was rather limited, with approximately 30 bacterial species associated with the algae. The majority of the observed microorganisms were affiliated with Alphaproteobacteria, Betaproteobacteria, and Bacteroidetes. A combined approach of sequencing via GS FLX Titanium from Roche and HiSeq 2000 from Illumina resulted in the overall production of 350 Mbp of sequenced DNA, 165 Mbp of which was assembled in larger contigs with a maximum size of 0.2 Mbp. A KEGG pathway analysis suggested high metabolic diversity with respect to the use of polymers and aromatic and nonaromatic compounds. Genes associated with the biosynthesis of essential B vitamins were highly redundant and functional. Moreover, a relatively high number of predicted and functional lipase and esterase genes indicated that the alga-associated bacteria are possibly a major sink for lipids and fatty acids produced by the microalgae. This is the first metagenome study of microalga- and PBR-associated biofilm bacteria, and it gives new clues for improved biofuel production in PBRs.

Microbial diversity and biochemical potential encoded by thermal spring metagenomes derived from the Kamchatka Peninsula.

Volcanic regions contain a variety of environments suitable for extremophiles. This study was focused on assessing and exploiting the prokaryotic diversity of two microbial communities derived from different Kamchatkian thermal springs by metagenomic approaches. Samples were taken from a thermoacidophilic spring near the Mutnovsky Volcano and from a thermophilic spring in the Uzon Caldera. Environmental DNA for metagenomic analysis was isolated from collected sediment samples by direct cell lysis. The prokaryotic community composition was examined by analysis of archaeal and bacterial 16S rRNA genes. A total number of 1235 16S rRNA gene sequences were obtained and used for taxonomic classification. Most abundant in the samples were members of Thaumarchaeota, Thermotogae, and Proteobacteria. The Mutnovsky hot spring was dominated by the Terrestrial Hot Spring Group, Kosmotoga, and Acidithiobacillus. The Uzon Caldera was dominated by uncultured members of the Miscellaneous Crenarchaeotic Group and Enterobacteriaceae. The remaining 16S rRNA gene sequences belonged to the Aquificae, Dictyoglomi, Euryarchaeota, Korarchaeota, Thermodesulfobacteria, Firmicutes, and some potential new phyla. In addition, the recovered DNA was used for generation of metagenomic libraries, which were subsequently mined for genes encoding lipolytic and proteolytic enzymes. Three novel genes conferring lipolytic and one gene conferring proteolytic activity were identified.

Driving forces of soil bacterial community structure, diversity, and function in temperate grasslands and forests

Soil bacteria provide a large range of ecosystem services such as nutrient cycling. Despite their important role in soil systems, compositional and functional responses of bacterial communities to different land use and management regimes are not fully understood. Here, we assessed soil bacterial communities in 150 forest and 150 grassland soils derived from three German regions by pyrotag sequencing of 16S rRNA genes. Land use type (forest and grassland) and soil edaphic properties strongly affected bacterial community structure and function, whereas management regime had a minor effect. In addition, a separation of soil bacterial communities by sampling region was encountered. Soil pH was the best predictor for bacterial community structure, diversity and function. The application of multinomial log-linear models revealed distinct responses of abundant bacterial groups towards pH. Predicted functional profiles revealed that differences in land use not only select for distinct bacterial populations but also for specific functional traits. The combination of 16S rRNA data and corresponding functional profiles provided comprehensive insights into compositional and functional adaptations to changing environmental conditions associated with differences in land use and management.

Bacterial community dynamics during polysaccharide degradation at contrasting sites in the Southern and Atlantic Oceans.

The bacterial degradation of polysaccharides is central to marine carbon cycling, but little is known about the bacterial taxa that degrade specific marine polysaccharides. Here, bacterial growth and community dynamics were studied during the degradation of the polysaccharides chitin, alginate and agarose in microcosm experiments at four contrasting locations in the Southern and Atlantic Oceans. At the Southern polar front, chitin-supplemented microcosms were characterized by higher fractions of actively growing cells and a community shift from Alphaproteobacteria to Gammaproteobacteria and Bacteroidetes. At the Antarctic ice shelf, chitin degradation was associated with growth of Bacteroidetes, with 24% higher cell numbers compared with the control. At the Patagonian continental shelf, alginate and agarose degradation covaried with growth of different Alteromonadaceae populations, each with specific temporal growth patterns. At the Mauritanian upwelling, only the alginate hydrolysis product guluronate was consumed, coincident with increasing abundances of Alteromonadaceae and possibly cross-feeding SAR11. 16S rRNA gene amplicon libraries indicated that growth of the Bacteroidetes-affiliated genus Reichenbachiella was stimulated by chitin at all cold and temperate water stations, suggesting comparable ecological roles over wide geographical scales. Overall, the predominance of location-specific patterns showed that bacterial communities from contrasting oceanic biomes have members with different potentials to hydrolyse polysaccharides.

Deciphering associations between dissolved organic molecules and bacterial communities in a pelagic marine system.

Dissolved organic matter (DOM) is the main substrate and energy source for heterotrophic bacterioplankton. To understand the interactions between DOM and the bacterial community (BC), it is important to identify the key factors on both sides in detail, chemically distinct moieties in DOM and the various bacterial taxa. Next-generation sequencing facilitates the classification of millions of reads of environmental DNA and RNA amplicons and ultrahigh-resolution mass spectrometry yields up to 10 000 DOM molecular formulae in a marine water sample. Linking this detailed biological and chemical information is a crucial first step toward a mechanistic understanding of the role of microorganisms in the marine carbon cycle. In this study, we interpreted the complex microbiological and molecular information via a novel combination of multivariate statistics. We were able to reveal distinct relationships between the key factors of organic matter cycling along a latitudinal transect across the North Sea. Total BC and DOM composition were mainly driven by mixing of distinct water masses and presumably retain their respective terrigenous imprint on similar timescales on their way through the North Sea. The active microbial community, however, was rather influenced by local events and correlated with specific DOM molecular formulae indicative of compounds that are easily degradable. These trends were most pronounced on the highest resolved level, that is, operationally defined ‘species’, reflecting the functional diversity of microorganisms at high taxonomic resolution.

RNA-based assessment of diversity and composition of active archaeal communities in the German Bight.

Archaea play an important role in various biogeochemical cycles. They are known extremophiles inhabiting environments such as thermal springs or hydrothermal vents. Recent studies have revealed a significant abundance of Archaea in moderate environments, for example, temperate sea water. Nevertheless, the composition and ecosystem function of these marine archaeal communities is largely unknown. To assess diversity and composition of active archaeal communities in the German Bight, seven marine water samples were taken and studied by RNA-based analysis of ribosomal 16S rRNA. For this purpose, total RNA was extracted from the samples and converted to cDNA. Archaeal community structures were investigated by pyrosequencing-based analysis of 16S rRNA amplicons generated from cDNA. To our knowledge, this is the first study combining next-generation sequencing and metatranscriptomics to study archaeal communities in marine habitats. The pyrosequencing-derived dataset comprised 62,045 archaeal 16S rRNA sequences. We identified Halobacteria as the predominant archaeal group across all samples with increased abundance in algal blooms. Thermoplasmatales (Euryarchaeota) and the Marine Group I (Thaumarchaeota) were identified in minor abundances. It is indicated that archaeal community patterns were influenced by environmental conditions.

A comparative metagenome survey of the fecal microbiota of a breast- and a plant-fed Asian elephant reveals an unexpectedly high diversity of glycoside hydrolase family enzymes.

A phylogenetic and metagenomic study of elephant feces samples (derived from a three-weeks-old and a six-years-old Asian elephant) was conducted in order to describe the microbiota inhabiting this large land-living animal. The microbial diversity was examined via 16S rRNA gene analysis. We generated more than 44,000 GS-FLX+454 reads for each animal. For the baby elephant, 380 operational taxonomic units (OTUs) were identified at 97% sequence identity level; in the six-years-old animal, close to 3,000 OTUs were identified, suggesting high microbial diversity in the older animal. In both animals most OTUs belonged to Bacteroidetes and Firmicutes. Additionally, for the baby elephant a high number of Proteobacteria was detected. A metagenomic sequencing approach using Illumina technology resulted in the generation of 1.1 Gbp assembled DNA in contigs with a maximum size of 0.6 Mbp. A KEGG pathway analysis suggested high metabolic diversity regarding the use of polymers and aromatic and non-aromatic compounds. In line with the high phylogenetic diversity, a surprising and not previously described biodiversity of glycoside hydrolase (GH) genes was found. Enzymes of 84 GH families were detected. Polysaccharide utilization loci (PULs), which are found in Bacteroidetes, were highly abundant in the dataset; some of these comprised cellulase genes. Furthermore the highest coverage for GH5 and GH9 family enzymes was detected for Bacteroidetes, suggesting that bacteria of this phylum are mainly responsible for the degradation of cellulose in the Asian elephant. Altogether, this study delivers insight into the biomass conversion by one of the largest plant-fed and land-living animals.

The green impact: bacterioplankton response toward a phytoplankton spring bloom in the southern North Sea assessed by comparative metagenomic and metatranscriptomic approaches

Phytoplankton blooms exhibit a severe impact on bacterioplankton communities as they change nutrient availabilities and other environmental factors. In the current study, the response of a bacterioplankton community to a Phaeocystis globosa spring bloom was investigated in the southern North Sea. For this purpose, water samples were taken inside and reference samples outside of an algal spring bloom. Structural changes of the bacterioplankton community were assessed by amplicon-based analysis of 16S rRNA genes and transcripts generated from environmental DNA and RNA, respectively. Several marine groups responded to bloom presence. The abundance of the Roseobacter RCA cluster and the SAR92 clade significantly increased in bloom presence in the total and active fraction of the bacterial community. Functional changes were investigated by direct sequencing of environmental DNA and mRNA. The corresponding datasets comprised more than 500 million sequences across all samples. Metatranscriptomic data sets were mapped on representative genomes of abundant marine groups present in the samples and on assembled metagenomic and metatranscriptomic datasets. Differences in gene expression profiles between non-bloom and bloom samples were recorded. The genome-wide gene expression level of Planktomarina temperata, an abundant member of the Roseobacter RCA cluster, was higher inside the bloom. Genes that were differently expressed included transposases, which showed increased expression levels inside the bloom. This might contribute to the adaptation of this organism toward environmental stresses through genome reorganization. In addition, several genes affiliated to the SAR92 clade were significantly upregulated inside the bloom including genes encoding for proteins involved in isoleucine and leucine incorporation. Obtained results provide novel insights into compositional and functional variations of marine bacterioplankton communities as response to a phytoplankton bloom.