Metatranscriptomics captures dynamic shifts in mycorrhizal coordination in boreal forests.

Carbon storage and cycling in boreal forests-the largest terrestrial carbon store-is moderated by complex interactions between trees and soil microorganisms. However, existing methods limit our ability to predict how changes in environmental conditions will alter these associations and the essential ecosystem services they provide. To address this, we developed a metatranscriptomic approach to analyze the impact of nutrient enrichment on Norway spruce fine roots and the community structure, function, and tree-microbe coordination of over 350 root-associated fungal species. In response to altered nutrient status, host trees redefined their relationship with the fungal community by reducing sugar efflux carriers and enhancing defense processes. This resulted in a profound restructuring of the fungal community and a collapse in functional coordination between the tree and the dominant Basidiomycete species, and an increase in functional coordination with versatile Ascomycete species. As such, there was a functional shift in community dominance from Basidiomycetes species, with important roles in enzymatically cycling recalcitrant carbon, to Ascomycete species that have melanized cell walls that are highly resistant to degradation. These changes were accompanied by prominent shifts in transcriptional coordination between over 60 predicted fungal effectors, with more than 5,000 Norway spruce transcripts, providing mechanistic insight into the complex molecular dialogue coordinating host trees and their fungal partners. The host-microbe dynamics captured by this study functionally inform how these complex and sensitive biological relationships may mediate the carbon storage potential of boreal soils under changing nutrient conditions.

Prevalent reliance of bacterioplankton on exogenous vitamin B1 and precursor availability.

Vitamin B1 (B1 herein) is a vital enzyme cofactor required by virtually all cells, including bacterioplankton, which strongly influence aquatic biogeochemistry and productivity and modulate climate on Earth. Intriguingly, bacterioplankton can be de novo B1 synthesizers or B1 auxotrophs, which cannot synthesize B1 de novo and require exogenous B1 or B1 precursors to survive. Recent isolate-based work suggests select abundant bacterioplankton are B1 auxotrophs, but direct evidence of B1 auxotrophy among natural communities is scant. In addition, it is entirely unknown if bulk bacterioplankton growth is ever B1-limited. We show by surveying for B1-related genes in estuarine, marine, and freshwater metagenomes and metagenome-assembled genomes (MAGs) that most naturally occurring bacterioplankton are B1 auxotrophs. Pyrimidine B1-auxotrophic bacterioplankton numerically dominated metagenomes, but multiple other B1-auxotrophic types and distinct uptake and B1-salvaging strategies were also identified, including dual (pyrimidine and thiazole) and intact B1 auxotrophs that have received little prior consideration. Time-series metagenomes from the Baltic Sea revealed pronounced shifts in the prevalence of multiple B1-auxotrophic types and in the B1-uptake and B1-salvaging strategies over time. Complementarily, we documented B1/precursor limitation of bacterioplankton production in three of five nutrient-amendment experiments at the same time-series station, specifically when intact B1 concentrations were ≤3.7 pM, based on bioassays with a genetically engineered Vibrio anguillarum B1-auxotrophic strain. Collectively, the data presented highlight the prevalent reliance of bacterioplankton on exogenous B1/precursors and on the bioavailability of the micronutrients as an overlooked factor that could influence bacterioplankton growth and succession and thereby the cycling of nutrients and energy in aquatic systems.

Meta-omic analyses of Baltic Sea cyanobacteria: diversity, community structure and salt acclimation.

Cyanobacteria are important phytoplankton in the Baltic Sea, an estuarine-like environment with pronounced north to south gradients in salinity and nutrient concentrations. Here, we present a metagenomic and -transcriptomic survey, with subsequent analyses targeting the genetic identity, phylogenetic diversity, and spatial distribution of Baltic Sea cyanobacteria. The cyanobacterial community constituted close to 12% of the microbial population sampled during a pre-bloom period (June-July 2009). The community was dominated by unicellular picocyanobacteria, specifically a few highly abundant taxa (Synechococcus and Cyanobium) with a long tail of low abundance representatives, and local peaks of bloom-forming heterocystous taxa. Cyanobacteria in the Baltic Sea differed genetically from those in adjacent limnic and marine waters as well as from cultivated and sequenced picocyanobacterial strains. Diversity peaked at brackish salinities 3.5-16 psu, with low N:P ratios. A shift in community composition from brackish to marine strains was accompanied by a change in the repertoire and expression of genes involved in salt acclimation. Overall, the pre-bloom cyanobacterial population was more genetically diverse, widespread and abundant than previously documented, with unicellular picocyanobacteria being the most abundant clade along the entire Baltic Sea salinity gradient.

Towards high-throughput parallel imaging and single-cell transcriptomics of microbial eukaryotic plankton.

Single-cell transcriptomics has the potential to provide novel insights into poorly studied microbial eukaryotes. Although several such technologies are available and benchmarked on mammalian cells, few have been tested on protists. Here, we applied a microarray single-cell sequencing (MASC-seq) technology, that generates microscope images of cells in parallel with capturing their transcriptomes, on three species representing important plankton groups with different cell structures; the ciliate Tetrahymena thermophila, the diatom Phaeodactylum tricornutum, and the dinoflagellate Heterocapsa sp. Both the cell fixation and permeabilization steps were adjusted. For the ciliate and dinoflagellate, the number of transcripts of microarray spots with single cells were significantly higher than for background spots, and the overall expression patterns were correlated with that of bulk RNA, while for the much smaller diatom cells, it was not possible to separate single-cell transcripts from background. The MASC-seq method holds promise for investigating "microbial dark matter", although further optimizations are necessary to increase the signal-to-noise ratio.

New chemical and microbial perspectives on vitamin B1 and vitamer dynamics of a coastal system.

Vitamin B1 (thiamin, B1) is an essential micronutrient for cells, yet intriguingly in aquatic systems most bacterioplankton are unable to synthesize it de novo (auxotrophy), requiring an exogenous source. Cycling of this valuable metabolite in aquatic systems has not been fully investigated and vitamers (B1-related compounds) have only begun to be measured and incorporated into the B1 cycle. Here, we identify potential key producers and consumers of B1 and gain new insights into the dynamics of B1 cycling through measurements of B1 and vitamers (HMP: 4-amino-5-hydroxymethyl-2-methylpyrimidine, HET: 4-methyl-5-thiazoleethanol, FAMP: N-formyl-4-amino-5-aminomethyl-2-methylpyrimidine) in the particulate and dissolved pool in a temperate coastal system. Dissolved B1 was not the primary limiting nutrient for bacterial production and was relatively stable across seasons with concentrations ranging from 74-117 pM, indicating a balance of supply and demand. However, vitamer concentration changed markedly with season as did transcripts related to vitamer salvage and transport suggesting use of vitamers by certain bacterioplankton, e.g. Pelagibacterales. Genomic and transcriptomic analyses showed that up to 78% of the bacterioplankton taxa were B1 auxotrophs. Notably, de novo B1 production was restricted to a few abundant bacterioplankton (e.g. Vulcanococcus, BACL14 (Burkholderiales), Verrucomicrobiales) across seasons. In summer, abundant picocyanobacteria were important putative B1 sources, based on transcriptional activity, leading to an increase in the B1 pool. Our results provide a new dynamic view of the players and processes involved in B1 cycling over time in coastal waters, and identify specific priority populations and processes for future study.

Whole genome sequence of the deep-sea sponge Geodia barretti (Metazoa, Porifera, Demospongiae).

Sponges are among the earliest branching extant animals. As such, genetic data from this group are valuable for understanding the evolution of various traits and processes in other animals. However, like many marine organisms, they are notoriously difficult to sequence, and hence, genomic data are scarce. Here, we present the draft genome assembly for the North Atlantic deep-sea high microbial abundance species Geodia barretti Bowerbank 1858, from a single individual collected on the West Coast of Sweden. The nuclear genome assembly has 4,535 scaffolds, an N50 of 48,447 bp and a total length of 144 Mb; the mitochondrial genome is 17,996 bp long. BUSCO completeness was 71.5%. The genome was annotated using a combination of ab initio and evidence-based methods finding 31,884 protein-coding genes.

Composition and Seasonality of Membrane Transporters in Marine Picoplankton.

In this study, we examined transporter genes in metagenomic and metatranscriptomic data from a time-series survey in the temperate marine environment of the Baltic Sea. We analyzed the abundance and taxonomic distribution of transporters in the 3µm-0.2µm size fraction comprising prokaryotes and some picoeukaryotes. The presence of specific transporter traits was shown to be guiding the succession of these microorganisms. A limited number of taxa were associated with the dominant transporter proteins that were identified for the nine key substrate categories for microbial growth. Throughout the year, the microbial taxa at the level of order showed highly similar patterns in terms of transporter traits. The distribution of transporters stayed the same, irrespective of the abundance of each taxon. This would suggest that the distribution pattern of transporters depends on the bacterial groups being dominant at a given time of the year. Also, we find notable numbers of secretion proteins that may allow marine bacteria to infect and kill prey organisms thus releasing nutrients. Finally, we demonstrate that transporter proteins may provide clues to the relative importance of biogeochemical processes, and we suggest that virtual transporter functionalities may become important components in future population dynamics models.

Comparative Fungal Community Analyses Using Metatranscriptomics and Internal Transcribed Spacer Amplicon Sequencing from Norway Spruce.

The health, growth, and fitness of boreal forest trees are impacted and improved by their associated microbiomes. Microbial gene expression and functional activity can be assayed with RNA sequencing (RNA-Seq) data from host samples. In contrast, phylogenetic marker gene amplicon sequencing data are used to assess taxonomic composition and community structure of the microbiome. Few studies have considered how much of this structural and taxonomic information is included in transcriptomic data from matched samples. Here, we described fungal communities using both host-derived RNA-Seq and fungal ITS1 DNA amplicon sequencing to compare the outcomes between the methods. We used a panel of root and needle samples from the coniferous tree species Picea abies (Norway spruce) growing in untreated (nutrient-deficient) and nutrient-enriched plots at the Flakaliden forest research site in boreal northern Sweden. We show that the relationship between samples and alpha and beta diversity indicated by the fungal transcriptome is in agreement with that generated by the ITS data, while also identifying a lack of taxonomic overlap due to limitations imposed by current database coverage. Furthermore, we demonstrate how metatranscriptomics data additionally provide biologically informative functional insights. At the community level, there were changes in starch and sucrose metabolism, biosynthesis of amino acids, and pentose and glucuronate interconversions, while processing of organic macromolecules, including aromatic and heterocyclic compounds, was enriched in transcripts assigned to the genus Cortinarius IMPORTANCE A deeper understanding of microbial communities associated with plants is revealing their importance for plant health and productivity. RNA extracted from plant field samples represents the host and other organisms present. Typically, gene expression studies focus on the plant component or, in a limited number of studies, expression in one or more associated organisms. However, metatranscriptomic data are rarely used for taxonomic profiling, which is currently performed using amplicon approaches. We created an assembly-based, reproducible, and hardware-agnostic workflow to taxonomically and functionally annotate fungal RNA-Seq data obtained from Norway spruce roots, which we compared to matching ITS amplicon sequencing data. While we identified some limitations and caveats, we show that functional, taxonomic, and compositional insights can all be obtained from RNA-Seq data. These findings highlight the potential of metatranscriptomics to advance our understanding of interaction, response, and effect between host plants and their associated microbial communities.

Energy efficiency and biological interactions define the core microbiome of deep oligotrophic groundwater.

While oligotrophic deep groundwaters host active microbes attuned to the low-end of the bioenergetics spectrum, the ecological constraints on microbial niches in these ecosystems and their consequences for microbiome convergence are unknown. Here, we provide a genome-resolved, integrated omics analysis comparing archaeal and bacterial communities in disconnected fracture fluids of the Fennoscandian Shield in Europe. Leveraging a dataset that combines metagenomes, single cell genomes, and metatranscriptomes, we show that groundwaters flowing in similar lithologies offer fixed niches that are occupied by a common core microbiome. Functional expression analysis highlights that these deep groundwater ecosystems foster diverse, yet cooperative communities adapted to this setting. We suggest that these communities stimulate cooperation by expression of functions related to ecological traits, such as aggregate or biofilm formation, while alleviating the burden on microorganisms producing compounds or functions that provide a collective benefit by facilitating reciprocal promiscuous metabolic partnerships with other members of the community. We hypothesize that an episodic lifestyle enabled by reversible bacteriostatic functions ensures the subsistence of the oligotrophic deep groundwater microbiome.

Genome sequence of segmented filamentous bacteria present in the human intestine.

Segmented filamentous bacteria (SFB) are unique immune modulatory bacteria colonizing the small intestine of a variety of animals in a host-specific manner. SFB exhibit filamentous growth and attach to the host's intestinal epithelium, offering a physical route of interaction. SFB affect functions of the host immune system, among them IgA production and T-cell maturation. Until now, no human-specific SFB genome has been reported. Here, we report the metagenomic reconstruction of an SFB genome from a human ileostomy sample. Phylogenomic analysis clusters the genome with SFB genomes from mouse, rat and turkey, but the genome is genetically distinct, displaying 65-71% average amino acid identity to the others. By screening human faecal metagenomic datasets, we identified individuals carrying sequences identical to the new SFB genome. We thus conclude that a unique SFB variant exists in humans and foresee a renewed interest in the elucidation of SFB functionality in this environment.

Identification of Novel Putative Bacterial Feruloyl Esterases From Anaerobic Ecosystems by Use of Whole-Genome Shotgun Metagenomics and Genome Binning.

Feruloyl esterases (FAEs) can reduce the recalcitrance of lignocellulosic biomass to enzymatic hydrolysis, thereby enhancing biorefinery potentials or animal feeding values of the biomass. In addition, ferulic acid, a product of FAE activity, has applications in pharmaceutical and food/beverage industries. It is therefore of great interest to identify new FAEs to enhance understanding about this enzyme family. For this purpose, we used whole-genome shotgun metagenomics and genome binning to explore rumens of dairy cows, large intestines of horses, sediments of freshwater and forest topsoils to identify novel prokaryotic FAEs and trace the responsible microorganisms. A number of prokaryotic genomes were recovered of which, genomes of Clostridiales order and Candidatus Rhabdochlamydia genus showed FAE coding capacities. In total, five sequences were deemed as putative FAE. The BLASTP search against non-redundant protein database of NCBI indicated that these putative FAEs represented novel sequences within this enzyme family. The phylogenetic analysis showed that at least three putative sequences shared evolutionary lineage with FAEs of type A and thus could possess specific activities similar to this type of FAEs, something that is not previously found outside fungal kingdom. We nominate Candidatus Rhabdochlamydia genus as a novel FAE producing taxonomic unit.

BARM and BalticMicrobeDB, a reference metagenome and interface to meta-omic data for the Baltic Sea.

The Baltic Sea is one of the world's largest brackish water bodies and is characterised by pronounced physicochemical gradients where microbes are the main biogeochemical catalysts. Meta-omic methods provide rich information on the composition of, and activities within, microbial ecosystems, but are computationally heavy to perform. We here present the Baltic Sea Reference Metagenome (BARM), complete with annotated genes to facilitate further studies with much less computational effort. The assembly is constructed using 2.6 billion metagenomic reads from 81 water samples, spanning both spatial and temporal dimensions, and contains 6.8 million genes that have been annotated for function and taxonomy. The assembly is useful as a reference, facilitating taxonomic and functional annotation of additional samples by simply mapping their reads against the assembly. This capability is demonstrated by the successful mapping and annotation of 24 external samples. In addition, we present a public web interface, BalticMicrobeDB, for interactive exploratory analysis of the dataset.

Diversity and Expression of Bacterial Metacaspases in an Aquatic Ecosystem.

Metacaspases are distant homologs of metazoan caspase proteases, implicated in stress response, and programmed cell death (PCD) in bacteria and phytoplankton. While the few previous studies on metacaspases have relied on cultured organisms and sequenced genomes, no studies have focused on metacaspases in a natural setting. We here present data from the first microbial community-wide metacaspase survey; performed by querying metagenomic and metatranscriptomic datasets from the brackish Baltic Sea, a water body characterized by pronounced environmental gradients and periods of massive cyanobacterial blooms. Metacaspase genes were restricted to ~4% of the bacteria, taxonomically affiliated mainly to Bacteroidetes, Alpha- and Betaproteobacteria and Cyanobacteria. The gene abundance was significantly higher in larger or particle-associated bacteria (>0.8 µm), and filamentous Cyanobacteria dominated metacaspase gene expression throughout the bloom season. Distinct seasonal expression patterns were detected for the three metacaspase genes in Nodularia spumigena, one of the main bloom-formers. Clustering of normalized gene expression in combination with analyses of genomic and assembly data suggest functional diversification of these genes, and possible roles of the metacaspase genes related to stress responses, i.e., sulfur metabolism in connection to oxidative stress, and nutrient stress induced cellular differentiation. Co-expression of genes encoding metacaspases and nodularin toxin synthesis enzymes was also observed in Nodularia spumigena. The study shows that metacaspases represent an adaptation of potentially high importance for several key organisms in the Baltic Sea, most prominently Cyanobacteria, and open up for further exploration of their physiological roles in microbes and assessment of their ecological impact in aquatic habitats.