Niche breadth affects bacterial transcription patterns along a salinity gradient.

Understanding the molecular mechanisms that determine a species' life history is important for predicting their susceptibility to environmental change. While specialist species with a narrow niche breadth (NB) maximize their fitness in their optimum habitat, generalists with broad NB adapt to multiple environments. The main objective of this study was to identify general transcriptional patterns that would distinguish bacterial strains characterized by contrasted NBs along a salinity gradient. More specifically, we hypothesized that genes encoding fitness-related traits, such as biomass production, have a higher degree of transcriptional regulation in specialists than in generalists, because the fitness of specialists is more variable under environmental change. By contrast, we expected that generalists would exhibit enhanced transcriptional regulation of genes encoding traits that protect them against cellular damage. To test these hypotheses, we assessed the transcriptional regulation of fitness-related and adaptation-related genes of 11 bacterial strains in relation to their NB and stress exposure under changing salinity conditions. The results suggested that transcriptional regulation levels of fitness- and adaptation-related genes correlated with the NB and/or the stress exposure of the inspected strains. We further identified a shortlist of candidate stress marker genes that could be used in future studies to monitor the susceptibility of bacterial populations or communities to environmental changes.

The environment drives microbial trait variability in aquatic habitats.

A prerequisite to improve the predictability of microbial community dynamics is to understand the mechanisms of microbial assembly. To study factors that contribute to microbial community assembly, we examined the temporal dynamics of genes in five aquatic metagenome time-series, originating from marine offshore or coastal sites and one lake. With this trait-based approach we expected to find gene-specific patterns of temporal allele variability that depended on the seasonal metacommunity size of carrier-taxa and the variability of the milieu and the substrates to which the resulting proteins were exposed. In more detail, we hypothesized that a larger seasonal metacommunity size would result in increased temporal variability of functional units (i.e., gene alleles), as shown previously for taxonomic units. We further hypothesized that multicopy genes would feature higher temporal variability than single-copy genes, as gene multiplication can result from high variability in substrate quality and quantity. Finally, we hypothesized that direct exposure of proteins to the extracellular environment would result in increased temporal variability of the respective gene compared to intracellular proteins that are less exposed to environmental fluctuations. The first two hypotheses were confirmed in all data sets, while significant effects of the subcellular location of gene products was only seen in three of the five time-series. The gene with the highest allele variability throughout all data sets was an iron transporter, also representing a target for phage infection. Previous work has emphasized the role of phage-prokaryote interactions as a major driver of microbial diversity. Our finding therefore points to a potentially important role of iron transporter-mediated phage infections for the assembly and maintenance of diversity in aquatic prokaryotes.

Microbial iron metabolism as revealed by gene expression profiles in contrasted Southern Ocean regimes.

Iron (Fe) is a limiting nutrient in large regions of the ocean, but the strategies of prokaryotes to cope with this micronutrient are poorly known. Using a gene-specific approach from metatranscriptomics data, we investigated seven Fe-related metabolic pathways in microbial communities from high nutrient low chlorophyll and naturally Fe-fertilized waters in the Southern Ocean. We observed major differences in the contribution of prokaryotic groups at different taxonomic levels to transcripts encoding Fe-uptake mechanisms, intracellular Fe storage and replacement and Fe-related pathways in the tricarboxylic acid (TCA) cycle. The composition of the prokaryotic communities contributing to the transcripts of a given Fe-related pathway was overall independent of the in situ Fe supply, indicating that microbial taxa utilize distinct Fe-related metabolic processes. Only a few prokaryotic groups contributed to the transcripts of more than one Fe-uptake mechanism, suggesting limited metabolic versatility. Taxa-specific expression of individual genes varied among prokaryotic groups and was substantially higher for all inspected genes in Fe-limited as compared to naturally fertilized waters, indicating the link between transcriptional state and Fe regime. Different metabolic strategies regarding low Fe concentrations in the Southern Ocean are discussed for two abundant prokaryotic groups, Pelagibacteraceae and Flavobacteriaceae.

A metatranscriptomics-based assessment of small-scale mixing of sulfidic and oxic waters on redoxcline prokaryotic communities.

Lateral intrusions of oxygen caused by small-scale mixing are thought to shape microbial activity in marine redoxclines. To examine the response of prokaryotes to such mixing events we employed a shipboard mixing experiment in the euxinic central Baltic Sea: oxic, nitrate containing and sulfidic water samples without detectable oxygenized substances were incubated directly or after mixing. While nitrate, nitrite and ammonium concentrations stayed approximately constant in all incubations, we observed a decrease of sulfide after the contact with oxygen in the sulfide containing incubations. The transcription of marker genes from chemolithoauthotrophic key players including archaeal nitrifiers as well as gammaproteobacterial and campylobacterial autotrophic organisms that couple denitrification with sulfur-oxidation were followed at four time points within 8.5 h. The temporally contrasting transcriptional profiles of gammaproteobacterial and campylobacterial denitrifiers that depend on the same inorganic substrates pointed to a niche separation. Particular archaeal and campylobacterial marker genes involved in nitrification, denitrification and sulfur oxidation, which depend on oxidized substrates, were highly upregulated in the anaerobic sulfidic samples. We suggest that, despite the absence of measurable oxygenated compounds in the sulfidic water, frequent intermittent small-scale intrusions stimulate the permanent upregulation of genes involved in nitrification, denitrification and sulfur oxidation.

Experimental insights into the importance of ecologically dissimilar bacteria to community assembly along a salinity gradient.

The response of local communities to marine-freshwater transitions and the processes that underlie community assembly are unclear, particularly with respect to bacteria that differ in their life strategies. Here, we implemented a transplant experiment where bacterioplankton from three regions of the Baltic Sea with differing salinities (∼3, 7 and 28 psu) were exposed to each other's environmental conditions. We found that habitat specialists were more abundant than generalists after exposure to salinity changes, irrespective of their origins. Most specialists that were selected following a salinity change were rare in the starting communities. Selection for generalists, however, was not specifically driven by the recruitment of either rare or abundant members, suggesting that taxon's initial abundance is minor relevant to the growth of generalists. Patterns in phylogenetic relatedness indicated that environmental filtering was the most influential assembly mechanism for specialists, whereas competitive interaction was more important for the assembly of generalists. Altogether, this study shows that large salinity changes promote the establishment of habitat specialists and that deterministic processes vary during community assembly for ecologically dissimilar taxa. We, therefore, propose that distinguishing assembly mechanisms of different community members helps understand and predict community dynamics in response to environmental change.

Interactions of Freshwater Cyanobacteria with Bacterial Antagonists.

Cyanobacterial and algal mass development, or blooms, have severe effects on freshwater and marine systems around the world. Many of these phototrophs produce a variety of potent toxins, contribute to oxygen depletion, and affect water quality in several ways. Coexisting antagonists, such as cyanolytic bacteria, hold the potential to suppress, or even terminate, such blooms, yet the nature of this interaction is not well studied. We isolated 31 cyanolytic bacteria affiliated with the genera Pseudomonas, Stenotrophomonas, Acinetobacter, and Delftia from three eutrophic freshwater lakes in Sweden and selected four phylogenetically diverse bacterial strains with strong-to-moderate lytic activity. To characterize their functional responses to the presence of cyanobacteria, we performed RNA sequencing (RNA-Seq) experiments on coculture incubations, with an initial predator-prey ratio of 1:1. Genes involved in central cellular pathways, stress-related heat or cold shock proteins, and antitoxin genes were highly expressed in both heterotrophs and cyanobacteria. Heterotrophs in coculture expressed genes involved in cell motility, signal transduction, and putative lytic activity. l,d-Transpeptidase was the only significantly upregulated lytic gene in Stenotrophomonas rhizophila EK20. Heterotrophs also shifted their central metabolism from the tricarboxylic acid cycle to the glyoxylate shunt. Concurrently, cyanobacteria clearly show contrasting antagonistic interactions with the four tested heterotrophic strains, which is also reflected in the physical attachment to their cells. In conclusion, antagonistic interactions with cyanobacteria were initiated within 24 h, and expression profiles suggest varied responses for the different cyanobacteria and studied cyanolytes.IMPORTANCE Here, we present how gene expression profiles can be used to reveal interactions between bloom-forming freshwater cyanobacteria and antagonistic heterotrophic bacteria. Species-specific responses in both heterotrophs and cyanobacteria were identified. The study contributes to a better understanding of the interspecies cellular interactions underpinning the persistence and collapse of cyanobacterial blooms.

Metatranscriptomic data reveal the effect of different community properties on multifunctional redundancy.

The assessment of functional redundancy (FR) is essential to understand community structure-function relationships because FR buffers the functional performance of communities against changes in community composition. We introduce a novel metatranscriptome-based approach to quantify FR, which permits multifunctional aspects to be addressed. FR among prokaryotes was ranked in water samples after exposure to changing salinity. FR was higher for functional categories with mostly broad functions shared among many taxa than for functional categories containing many narrow functions. Furthermore, community characteristics had a higher impact on FR than environmental conditions. The metric also allows FR to be estimated between selected groups of taxa, and FR was high between more closely related organisms if communities were grown in similar environmental conditions. Overall, our data revealed a pronounced influence of functional diversity on the one hand but also the characteristics of individual community members on FR, which was specifically high in those communities whose members were more sensitive to salinity changes.

The transcriptional response of prokaryotes to phytoplankton-derived dissolved organic matter in seawater.

To better understand the functional responses in prokaryotes to dissolved organic matter (DOM), we compared the transcriptional pattern of natural prokaryotic communities grown in continuous cultures on seawater amended with phytoplankton-derived DOM. Metatranscriptomic reads were classified taxonomically (by genomic binning) and functionally (using Kyoto Encyclopedia of Genes and Genomes), and the relative gene expression of individual taxa (genome bins) was compared with the total community response. In the first experiment comparing seawater and seawater amended with diatom-derived DOM, metatranscriptomes revealed pronounced differences in pathways involved in carbohydrate and lipid metabolism. In the second experiment comparing seawater amended with cyanobacteria- and diatom-derived DOM, metatranscriptomes had similar functional profiles, likely reflecting more similar DOM regimes in this experimental setup. Among the five most abundant taxa investigated in more detail, two featured pronounced differences in transcript abundance between treatments suggesting that they were specialized in the use of only one of the two DOM regimes. However, these two taxa were less involved in carbohydrate metabolism than others and had few genes that were significantly regulated in response to the DOM source. Our results indicate that both substrate composition and the competitive interplay of community members were decisive for the functional response of a microbial system.

The cost of adaptability: resource availability constrains functional stability under pulsed disturbances.

Global change exposes ecosystems to changes in the frequency, magnitude, and concomitancy of disturbances, which impact the composition and functioning of these systems. Here, we experimentally evaluate the effects of salinity disturbances and eutrophication on bacterial communities from coastal ecosystems. The functional stability of these communities is critically important for maintaining water quality, productivity, and ecosystem services, such as fishery yields. Microbial functional stability can be maintained via resistance and resilience, which are reflected in genomic traits such as genome size and codon usage bias and may be linked to metabolic costs. However, little is known about the mechanisms that select these traits under varying nutrient regimes. To study the impact of pulsed disturbances on community assembly and functioning depending on metabolic costs, we performed a 41-day pulse disturbance experiment across two levels of resource availability. Our setup triggered stochastic community re-assembly processes in all treatments. In contrast, we observed consistent and resource availability-dependent patterns of superordinate community functioning and structural patterns, such as functional resistance in response to disturbances, genomic trait distributions, and species diversity. Predicted genomic traits reflected the selection for taxa possessing resistant- and resilience-related traits, particularly under high nutrient availability. Our findings are a step toward unraveling the compositional and genomic underpinnings of functional resistance in microbial communities after exposure to consecutive pulse disturbances. Our work demonstrates how resource availability alleviates metabolic constraints on resistance and resilience, and this has important consequences for predicting water quality and ecosystem productivity of environments exposed to global change. IMPORTANCE: Understanding the communities' responses to disturbances is a prerequisite to predicting ecosystem dynamics and, thus, highly relevant considering global change. Microbial communities play key roles in numerous ecosystem functions and services, and the large diversity, rapid growth, and phenotypic plasticity of microorganisms are thought to allow high resistance and resilience. While potential metabolic costs associated with adaptations to fluctuating environments have been debated, little evidence supports trade-offs between resource availability, resistance, and resilience. Here, we experimentally assessed the compositional and functional responses of an aquatic microbial model community to disturbances and systematically manipulated resource availability. Our results demonstrate that the capacity to tolerate environmental fluctuations is constrained by resource availability and reflected in the selection of genomic traits.

Pronounced seasonal dynamics of freshwater chitinase genes and chitin processing.

Seasonal variation in activity of enzymes involved in polymer degradation, including chitinases, has been observed previously in freshwater environments. However, it is not known whether the seasonal dynamics are due to shifts in the activity of bacteria already present, or shifts in community structure towards emergence or disappearance of chitinolytic organisms. We traced seasonal shifts in the chitinase gene assemblage in a temperate lake and linked these communities to variation in chitinase activity. Chitinase genes from 20 samples collected over a full yearly cycle were characterized by pyrosequencing. Pronounced temporal shifts in composition of the chitinase gene pool (beta diversity) occurred along with distinct shifts in richness (alpha diversity) as well as chitin processing. Changes in the chitinase gene pool correlated mainly with temperature, abundance of crustacean zooplankton and phytoplankton blooms. Also changes in the physical structure of the lake, e.g. stratification and mixing were associated with changes in the chitinolytic community, while differences were minor between surface and suboxic hypolimnetic water. The lake characteristics influencing the chitinolytic community are all linked to changes in organic particles and we suggest that seasonal changes in particle quality and availability foster microbial communities adapted to efficiently degrade them.

Trait-trait relationships and tradeoffs vary with genome size in prokaryotes.

We report genomic traits that have been associated with the life history of prokaryotes and highlight conflicting findings concerning earlier observed trait correlations and tradeoffs. In order to address possible explanations for these contradictions we examined trait-trait variations of 11 genomic traits from ~18,000 sequenced genomes. The studied trait-trait variations suggested: (i) the predominance of two resistance and resilience-related orthogonal axes and (ii) at least in free living species with large effective population sizes whose evolution is little affected by genetic drift an overlap between a resilience axis and an oligotrophic-copiotrophic axis. These findings imply that resistance associated traits of prokaryotes are globally decoupled from resilience related traits and in the case of free-living communities also from traits associated with resource availability. However, further inspection of pairwise scatterplots showed that resistance and resilience traits tended to be positively related for genomes up to roughly five million base pairs and negatively for larger genomes. Genome size distributions differ across habitats and our findings therefore point to habitat dependent tradeoffs between resistance and resilience. This in turn may preclude a globally consistent assignment of prokaryote genomic traits to the competitor - stress-tolerator - ruderal (CSR) schema that sorts species depending on their location along disturbance and productivity gradients into three ecological strategies and may serve as an explanation for conflicting findings from earlier studies. All reviewed genomic traits featured significant phylogenetic signals and we propose that our trait table can be applied to extrapolate genomic traits from taxonomic marker genes. This will enable to empirically evaluate the assembly of these genomic traits in prokaryotic communities from different habitats and under different productivity and disturbance scenarios as predicted via the resistance-resilience framework formulated here.

Global phylogeography of chitinase genes in aquatic metagenomes.

Phylogeny-based analysis of chitinase and 16S rRNA genes from metagenomic data suggests that salinity is a major driver for the distribution of both chitinolytic and total bacterial communities in aquatic systems. Additionally, more acidic chitinase proteins were observed with increasing salinity. Congruent habitat separation was further observed for both genes according to latitude and proximity to the coastline. However, comparison of chitinase and 16S rRNA genes extracted from different geographic locations showed little congruence in distribution. There was no indication that dispersal limited the global distribution of either gene.

Cryopreservation and Resuscitation of Natural Aquatic Prokaryotic Communities.

Experimental reproducibility in aquatic microbial ecology is critical to predict the dynamics of microbial communities. However, controlling the initial composition of naturally occurring microbial communities that will be used as the inoculum in experimental setups is challenging, because a proper method for the preservation of those communities is lacking. To provide a feasible method for preservation and resuscitation of natural aquatic prokaryote assemblages, we developed a cryopreservation procedure applied to natural aquatic prokaryotic communities. We studied the impact of inoculum size, processing time, and storage time on the success of resuscitation. We further assessed the effect of different growth media supplemented with dissolved organic matter (DOM) prepared from naturally occurring microorganisms on the recovery of the initially cryopreserved communities obtained from two sites that have contrasting trophic status and environmental heterogeneity. Our results demonstrated that the variability of the resuscitation process among replicates decreased with increasing inoculum size. The degree of similarity between initial and resuscitated communities was influenced by both the growth medium and origin of the community. We further demonstrated that depending on the inoculum source, 45-72% of the abundant species in the initially natural microbial communities could be detected as viable cells after cryopreservation. Processing time and long-term storage up to 12 months did not significantly influence the community composition after resuscitation. However, based on our results, we recommend keeping handling time to a minimum and ensure identical incubation conditions for repeated resuscitations from cryo-preserved aliquots at different time points. Given our results, we recommend cryopreservation as a promising tool to advance experimental research in the field of microbial ecology.

Ecosystem-wide metagenomic binning enables prediction of ecological niches from genomes.

The genome encodes the metabolic and functional capabilities of an organism and should be a major determinant of its ecological niche. Yet, it is unknown if the niche can be predicted directly from the genome. Here, we conduct metagenomic binning on 123 water samples spanning major environmental gradients of the Baltic Sea. The resulting 1961 metagenome-assembled genomes represent 352 species-level clusters that correspond to 1/3 of the metagenome sequences of the prokaryotic size-fraction. By using machine-learning, the placement of a genome cluster along various niche gradients (salinity level, depth, size-fraction) could be predicted based solely on its functional genes. The same approach predicted the genomes' placement in a virtual niche-space that captures the highest variation in distribution patterns. The predictions generally outperformed those inferred from phylogenetic information. Our study demonstrates a strong link between genome and ecological niche and provides a conceptual framework for predictive ecology based on genomic data.

High ratio of bacteriochlorophyll biosynthesis genes to chlorophyll biosynthesis genes in bacteria of humic lakes.

Recent studies highlight the diversity and significance of marine phototrophic microorganisms such as picocyanobacteria, phototrophic picoeukaryotes, and bacteriochlorophyll- and rhodopsin-holding phototrophic bacteria. To assess if freshwater ecosystems also harbor similar phototroph diversity, genes involved in the biosynthesis of bacteriochlorophyll and chlorophyll were targeted to explore oxygenic and aerobic anoxygenic phototroph composition in a wide range of lakes. Partial dark-operative protochlorophyllide oxidoreductase (DPOR) and chlorophyllide oxidoreductase (COR) genes in bacteria of seven lakes with contrasting trophic statuses were PCR amplified, cloned, and sequenced. Out of 61 sequences encoding the L subunit of DPOR (L-DPOR), 22 clustered with aerobic anoxygenic photosynthetic bacteria, whereas 39 L-DPOR sequences related to oxygenic phototrophs, like cyanobacteria, were observed. Phylogenetic analysis revealed clear separation of these freshwater L-DPOR genes as well as 11 COR gene sequences from their marine counterparts. Terminal restriction fragment length analysis of L-DPOR genes was used to characterize oxygenic aerobic and anoxygenic photosynthesizing populations in 20 lakes differing in physical and chemical characteristics. Significant differences in L-DPOR community composition were observed between dystrophic lakes and all other systems, where a higher proportion of genes affiliated with aerobic anoxygenic photosynthetic bacteria was observed than in other systems. Our results reveal a significant diversity of phototrophic microorganisms in lakes and suggest niche partitioning of oxygenic and aerobic anoxygenic phototrophs in these systems in response to trophic status and coupled differences in light regime.

Effect of large magnetotactic bacteria with polyphosphate inclusions on the phosphate profile of the suboxic zone in the Black Sea.

The Black Sea is the world's largest anoxic basin and a model system for studying processes across redox gradients. In between the oxic surface and the deeper sulfidic waters there is an unusually broad layer of 10-40 m, where neither oxygen nor sulfide are detectable. In this suboxic zone, dissolved phosphate profiles display a pronounced minimum at the upper and a maximum at the lower boundary, with a peak of particulate phosphorus in between, which was suggested to be caused by the sorption of phosphate on sinking particles of metal oxides. Here we show that bacterial polyphosphate inclusions within large magnetotactic bacteria related to the genus Magnetococcus contribute substantially to the observed phosphorus peak, as they contain 26-34% phosphorus compared to only 1-5% in metal-rich particles. Furthermore, we found increased gene expression for polyphosphate kinases by several groups of bacteria including Magnetococcaceae at the phosphate maximum, indicating active bacterial polyphosphate degradation. We propose that large magnetotactic bacteria shuttle up and down within the suboxic zone, scavenging phosphate at the upper and releasing it at the lower boundary. In contrast to a passive transport via metal oxides, this bacterial transport can quantitatively explain the observed phosphate profiles.

Bacterial community composition in Central European running waters examined by temperature gradient gel electrophoresis and sequence analysis of 16S rRNA genes.

The bacterial community composition in small streams and a river in central Germany was examined by temperature gradient gel electrophoresis (TGGE) with PCR products of 16S rRNA gene fragments and sequence analysis. Complex TGGE band patterns suggested high levels of diversity of bacterial species in all habitats of these environments. Cluster analyses demonstrated distinct differences among the communities in stream and spring water, sandy sediments, biofilms on stones, degrading leaves, and soil. The differences between stream water and sediment were more significant than those between sites within the same habitat along the stretch from the stream source to the mouth. TGGE data from an entire stream course suggest that, in the upper reach of the stream, a special suspended bacterial community is already established and changes only slightly downstream. The bacterial communities in water and sediment in an acidic headwater with a pH below 5 were highly similar to each other but deviated distinctly from the communities at the other sites. As ascertained by nucleotide sequence analysis, stream water communities were dominated by Betaproteobacteria (one-third of the total bacteria), whereas sediment communities were composed mainly of Betaproteobacteria and members of the Fibrobacteres/Acidobacteria group (each accounting for about 25% of bacteria). Sequences obtained from bacteria from water samples indicated the presence of typical cosmopolitan freshwater organisms. TGGE bands shared between stream and soil samples, as well as sequences found in bacteria from stream samples that were related to those of soil bacteria, demonstrated the occurrence of some species in both stream and soil habitats. Changes in bacterial community composition were correlated with geographic distance along a stream, but in comparisons of different streams and rivers, community composition was correlated only with environmental conditions.

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.

Phenotypic plasticity in heterotrophic marine microbial communities in continuous cultures.

Phenotypic plasticity (PP) is the development of alternate phenotypes of a given taxon as an adaptation to environmental conditions. Methodological limitations have restricted the quantification of PP to the measurement of a few traits in single organisms. We used metatranscriptomic libraries to overcome these challenges and estimate PP using the expressed genes of multiple heterotrophic organisms as a proxy for traits in a microbial community. The metatranscriptomes captured the expression response of natural marine bacterial communities grown on differing carbon resource regimes in continuous cultures. We found that taxa with different magnitudes of PP coexisted in the same cultures, and that members of the order Rhodobacterales had the highest levels of PP. In agreement with previous studies, our results suggest that continuous culturing may have specifically selected for taxa featuring a rather high range of PP. On average, PP and abundance changes within a taxon contributed equally to the organism's change in functional gene abundance, implying that both PP and abundance mediated observed differences in community function. However, not all functional changes due to PP were directly reflected in the bulk community functional response: gene expression changes in individual taxa due to PP were partly masked by counterbalanced expression of the same gene in other taxa. This observation demonstrates that PP had a stabilizing effect on a community's functional response to environmental change.