Search

Abstract In the face of climate change, the accelerating nature crisis, and other anthropogenic impacts, conserving biodiversity demands effective monitoring. While innovative approaches have emerged to enhance biodiversity assessment, significant data gaps persist, particularly within marine ecosystems. In this study, we assess the utility of combining citizen (community) science with environmental DNA (eDNA) metabarcoding for characterizing and quantifying marine fish biodiversity. Over the summer of 2022, 32 volunteers conducted extensive water sampling in a large Norwegian fjord, yielding 96 samples. Contrasting eDNA findings with conventional observational surveys (such as a national species registration database and beach seine surveys) unveiled a substantial overlap in recorded species inventories. The eDNA citizen science initiative identified previously undocumented species and rediscovered others emphasizing an increase in warm‐water species within the study area. Additionally, eDNA data unveiled reduced diversity within the inner fjord relative to the outer fjord. To conclude, our study demonstrates the successful integration of eDNA within a citizen science framework, facilitating comprehensive biodiversity tracking across coastal marine regions. These findings hold promise for advancing marine conservation efforts by providing valuable data to inform critical decision‐making processes.

Abstract Despite the importance of bacteria in aquatic ecosystems and their predictable diversity patterns across space and time, biomonitoring tools for status assessment relying on these organisms are widely lacking. This is partly due to insufficient data and models to identify reliable microbial predictors. Here, we show metabarcoding in combination with multivariate statistics and machine learning allows to identify bacterial bioindicators for existing biological status classification systems. Bacterial beta-diversity dynamics follow environmental gradients and the observed associations highlight potential bioindicators for ecological outcomes. Spatio-temporal links spanning the microbial communities along the river allow accurate prediction of downstream biological status from upstream information. Network analysis on amplicon sequence veariants identify as good indicators genera Fluviicola, Acinetobacter, Flavobacterium, and Rhodoluna, and reveal informational redundancy among taxa, which coincides with taxonomic relatedness. The redundancy among bacterial bioindicators reveals mutually exclusive taxa, which allow accurate biological status modeling using as few as 2–3 amplicon sequence variants. As such our models show that using a few bacterial amplicon sequence variants from globally distributed genera allows for biological status assessment along river systems.

The leaching of additives from plastics and elastomers (rubbers) has raised concerns due to their potential negative impacts on the environment and the development of antibiotic resistance. In this study, we investigated the effects of chemicals extracted from two types of rubber on microbiomes derived from a benthic sea urchin and two pelagic fish species. Additionally, we examined whether bacterial communities preconditioned with rubber-associated chemicals displayed adaptations to antibiotics. At the highest tested concentrations of chemicals, we observed reduced maximum growth rates and yields, prolonged lag phases, and increased alpha diversity. While the effects on alpha and beta diversity were not always conclusive, several bacterial genera were significantly influenced by chemicals from the two rubber sources. Subsequent exposure of sea urchin microbiomes preconditioned with rubber chemicals to the antibiotic ciprofloxacin resulted in decreased maximum growth rates. This indicates a more sensitive microbiome to ciprofloxacin when preconditioned with rubber chemicals. Although no significant interaction effects between rubber chemicals and ciprofloxacin exposure were observed in bacterial alpha and beta diversity, we observed log-fold changes in two bacterial genera in response to ciprofloxacin exposure. These findings highlight the structural and functional alterations in microbiomes originating from various marine species when exposed to rubber-associated chemicals and underscore the potential risks posed to marine life.

Abstract Due to climate warming, ice sheets around the world are losing mass, contributing to changes across terrestrial landscapes on decadal time spans. However, landscape repercussions on climate are poorly constrained mostly due to limited knowledge on microbial responses to deglaciation. Here, we reveal the genomic succession from chemolithotrophy to photo- and heterotrophy and increases in methane supersaturation in freshwater lakes upon glacial retreat. Arctic lakes at Svalbard also revealed strong microbial signatures form nutrient fertilization by birds. Although methanotrophs were present and increased along lake chronosequences, methane consumption rates were low even in supersaturated systems. Nitrous oxide oversaturation and genomic information suggest active nitrogen cycling across the entire deglaciated landscape, and in the high Arctic, increasing bird populations serve as major modulators at many sites. Our findings show diverse microbial succession patterns, and trajectories in carbon and nitrogen cycle processes representing a positive feedback loop of deglaciation on climate warming.

Abstract Lakes are significant players for the global climate since they sequester terrestrially derived dissolved organic carbon (DOC), and emit greenhouse gases like CO2 to the atmosphere. However, the differences in environmental drivers of CO2 concentrations are not well constrained along latitudinal and thus climate gradients. Our aim here is to provide a better understanding of net heterotrophy and gas balance at the catchment scale in a set of boreal, sub-Arctic and high-Arctic lakes. We assessed water chemistry and concentrations of dissolved O2 and CO2, as well as the CO2:O2 ratio in three groups of lakes separated by steps of approximately 10 degrees latitude in South-Eastern Norway (near 60° N), sub-Arctic lakes in the northernmost part of the Norwegian mainland (near 70° N) and high-Arctic lakes on Svalbard (near 80° N). Across all regions, CO2 saturation levels varied more (6–1374%) than O2 saturation levels (85–148%) and hence CO2 saturation governed the CO2:O2 ratio. The boreal lakes were generally undersaturated with O2, while the sub-Arctic and high-Arctic lakes ranged from O2 saturated to oversaturated. Regardless of location, the majority of the lakes were CO2 supersaturated. In the boreal lakes the CO2:O2 ratio was mainly related to DOC concentration, in contrast to the sub-Arctic and high-Arctic localities, where conductivity was the major statistical determinant. While the southern part is dominated by granitic and metamorphic bedrock, the sub-Arctic sites are scattered across a range of granitic to sedimentary bed rocks, and the majority of the high-Arctic lakes are situated on limestone, resulting in contrasting lake alkalinities between the regions. DOC dependency of the CO2:O2 ratio in the boreal region together with low alkalinity suggests that in-lake heterotrophic respiration was a major source of lake CO2. Contrastingly, the conductivity dependency indicates that CO2 saturation in the sub-Arctic and high-Arctic lakes was to a large part explained by DIC input from catchment respiration and carbonate weathering.

Freshwaters represent the most threatened environments with regard to biodiversity loss and, therefore, there is a need for national monitoring programs to effectively document species distribution and evaluate potential risks for vulnerable species. The monitoring of species for effective management practices is, however, challenged by insufficient data acquisition when using traditional methods. Here we present the application of environmental DNA (eDNA) metabarcoding of amphibians in combination with quantitative PCR (qPCR) assays for an invasive pathogenic chytrid species (Batrachochytrium dendrobatidis -Bd), a potential threat to endemic and endangered amphibian species. Statistical comparison of amphibian species detection using either traditional or eDNA-based approaches showed weak correspondence. By tracking the distribution of Bd over three years, we concluded that the risk for amphibian extinction is low since Bd was only detected at five sites where multiple amphibians were present over the sampled years. Our results show that eDNA-based detection can be used for simultaneous monitoring of amphibian diversity and the presence of amphibian pathogens at the national level in order to assess potential species extinction risks and establish effective management practices. As such our study represents suggestions for a national monitoring program based on eDNA.

Abstract The resistance and resilience provided by functional redundancy, a common feature of microbial communities, is not always advantageous. An example is nitrite oxidation in partial nitritation-anammox (PNA) reactors designed for nitrogen removal in wastewater treatment, where suppression of nitrite oxidizers like Nitrospira is sought. In these ecosystems, biofilms provide microhabitats with oxygen gradients, allowing the coexistence of aerobic and anaerobic bacteria. We designed a disturbance experiment where PNA biofilms, treating water from a high-rate activated sludge process, were constantly or intermittently exposed to anaerobic sidestream wastewater, which has been proposed to inhibit nitrite oxidizers. With increasing sidestream exposure we observed decreased abundance, alpha-diversity, functional versatility, and hence functional redundancy, among Nitrospira in the PNA biofilms, while the opposite patterns were observed for anammox bacteria within Brocadia. At the same time, species turnover was observed for aerobic ammonia-oxidizing Nitrosomonas populations. The different exposure regimens were associated with metagenomic assembled genomes of Nitrosomonas, Nitrospira, and Brocadia, encoding genes related to N-cycling, substrate usage, and osmotic stress response, possibly explaining the three different patterns by niche differentiation. These findings imply that disturbances can be used to manage the functional redundancy of biofilm microbiomes in a desirable direction, which should be considered when designing operational strategies for wastewater treatment.

Boreal lakes are the most abundant lakes on Earth. Changes in acid rain deposition, climate, and catchment land use have increased lateral fluxes of terrestrial dissolved organic matter (DOM), resulting in a widespread browning of boreal freshwaters. This browning affects the aqueous communities and ecosystem processes, and boost emissions of the greenhouse gases (GHG) CH4, CO2, and N2O. In this study, we predicted biotic saturation of GHGs in boreal lakes by using a set of chemical, hydrological, climate, and land use parameters. For this purpose, concentrations of GHGs and nutrients (organic C, -P, and -N) were determined in surface water samples from 73 lakes in south-eastern Norway covering wide ranges in DOM and nutrient concentrations, as well as catchment properties and land use. The spatial variation in saturation of each GHG is related to explanatory variables. Catchment characteristics (hydrological and climate parameters) such as lake size and summer precipitation, as well as NDVI, were key determinants when fitting GAM models for CH4 and CO2 saturation (explaining 71 and 54%, respectively), while summer precipitation and land use data were the best predictors for the N2O saturation, explaining almost 50% of deviance. Our results suggest that lake size, precipitation, and terrestrial primary production in the watershed control the saturation of GHG in boreal lakes. These predictions based on the 73-lake dataset was validated against an independent dataset from 46 lakes in the same region. Together, this provides an improved understanding of drivers and spatial variation in GHG saturation in boreal lakes across wide gradients of lake and catchment properties. The assessment highlights the need to incorporate multiple explanatory parameters in prediction models of GHGs for extrapolation across the boreal biome.

The current trend of increasing input of terrestrially derived dissolved organic carbon (DOC) to boreal freshwater systems is causing increased levels of carbon dioxide (CO2) supersaturation and degassing. Phosphorus (P) is often the most limiting nutrient for bacterial growth and would thus be expected to increase overall mineralization rates and CO2 production. However, high carbon (C) to P ratios of terrestrially derived DOC could also cause elevated cell-specific respiration of the excess C in heterotrophic bacteria. Using data from a survey of 75 Scandinavian lakes along an ecosystem gradient of DOC, we estimated in situ CO2 production rates. These rates showed a unimodal response with DOC-specific CO2 production negatively related to DOC:total phosphorus (TP) ratio, and a turning point at 5 mg C L−1, indicating higher DOC turnover rates in productive than in unproductive lakes. To further assess the dependency of bacterial respiration (BR) on DOC and P, we monitored CO2 production in incubations of water with a gradient of DOC crossed with two levels of inorganic P. Finally, we crossed DOC and P with a temperature gradient to test the temperature dependency of respiration rates [as oxygen (O2) consumption]. While total CO2 production seemed to be unaffected by P additions, respiration rates, and growth yields, as estimated by ribosomal gene copy numbers, suggest increased bacterial growth and decreased cell-specific respiration under non-limited P conditions. Respiration rates showed a sigmoid response to increasing DOC availability reaching a plateau at about 20 mg C L−1 of initial DOC concentrations. In addition to these P and DOC level effects, respiration rates responded in a non-monotonic fashion to temperature with an increase in respiration rates by a factor of 2.6 (±0.2) from 15 to 25°C and a decrease above 30°C. The combined results from the survey and experiments highlight DOC as the major determinant of CO2 production in boreal lakes, with P and temperature as significant modulators of respiration kinetics.

Aquatic N-fixation is generally associated with the growth and mass development of Cyanobacteria in nitrogen-deprived photic zones. However, sequenced genomes and environmental surveys suggest active aquatic N-fixation also by many non-cyanobacterial groups. Here, we revealed the seasonal variation and genomic diversity of potential N-fixers in a humic bog lake using metagenomic data and nif gene clusters analysis. Groups with diazotrophic operons were functionally divergent and included Cholorobi, Geobacter, Desulfobacterales, Methylococcales, and Acidobacteria. In addition to nifH (a gene that encodes the dinitrogenase reductase component of the molybdenum nitrogenase), we also identified sequences corresponding to vanadium and iron-only nitrogenase genes. Within the Chlorobi population, the nitrogenase (nifH) cluster was included in a well-structured retrotransposon. Furthermore, the presence of light-harvesting photosynthesis genes implies that anoxygenic photosynthesis may fuel nitrogen fixation under the prevailing low-irradiance conditions. The presence of rnf genes (related to the expression of H+/Na+-translocating ferredoxin: NAD+ oxidoreductase) in Methylococcales and Desulfobacterales suggests that other energy-generating processes may drive the costly N-fixation in the absence of photosynthesis. The highly reducing environment of the anoxic bottom layer of Trout Bog Lake may thus also provide a suitable niche for active N-fixers and primary producers. While future studies on the activity of these potential N-fixers are needed to clarify their role in freshwater nitrogen cycling, the metagenomic data presented here enabled an initial characterization of previously overlooked diazotrophs in freshwater biomes.

The composition of intestinal microbiota commonly varies among animal hosts and may affect host health. However, we have limited knowledge about the different relative roles of assembly processes, such as drift, dispersal and environmental selection, for the composition of gut microbiota. Here, we conducted a field study analyzing intestinal microbial communities of two fish species that either have (perch) or lack (roach) a stomach. We used a suite of statistical tools to evaluate the role of different assembly processes for intestine microbiota, including null model analysis (Chase et al., 2011; Fine and Kembel, 2011; Stegen et al., 2013), SourceTracker analysis (Knights et al., 2011) and several multivariate analyses, such as pRDA and PLS analysis. Drift, dispersal (i.e., microbes associated with food sources) and environmental factors (i.e., diet, host habitats), appeared to be of equal importance for the assembly of intestinal microbial communities in roach, while drift appeared most important in perch, followed by dispersal and environmental selection. Furthermore, we found that microbes associated with macroinvertebrates had a positive association to fish body condition (weight/length3) whereas microbes associated with zooplankton had a negative association to fish body condition. These results emphasize the important combined roles of drift, dispersal and environmental selection in shaping the host-associated microbial communities. We conclude that general conclusions about fish as a whole are not justified since different species differ in the relative roles of these important drivers of community assembly.

Abstract Background Gut microbiota provide functions of importance to influence hosts’ food digestion, metabolism, and protection against pathogens. Factors that affect the composition and functions of gut microbial communities are well studied in humans and other animals; however, we have limited knowledge of how natural food web factors such as stress from predators and food resource rations could affect hosts’ gut microbiota and how it interacts with host sex. In this study, we designed a two-factorial experiment exposing perch (Perca fluviatilis) to a predator (pike, Esox lucius), and different food ratios, to examine the compositional and functional changes of perch gut microbiota based on 16S rRNA amplicon sequencing. We also investigated if those changes are host sex dependent. Results We showed that overall gut microbiota composition among individual perch significantly responded to food ration and predator presence. We found that species richness decreased with predator presence, and we identified 23 taxa from a diverse set of phyla that were over-represented when a predator was present. For example, Fusobacteria increased both at the lowest food ration and at predation stress conditions, suggesting that Fusobacteria are favored by stressful situations for the host. In concordance, both food ration and predation stress seemed to influence the metabolic repertoire of the gut microbiota, such as biosynthesis of other secondary metabolites, metabolism of cofactors, and vitamins. In addition, the identified interaction between food ration and sex emphasizes sex-specific responses to diet quantity in gut microbiota. Conclusions Collectively, our findings emphasize an alternative state in gut microbiota with responses to changes in natural food webs depending on host sex. The obtained knowledge from this study provided us with an important perspective on gut microbiota in a food web context.

ABSTRACT Oxygen-stratified lakes are typical for the boreal zone and also a major source of greenhouse gas emissions in the region. Due to shallow light penetration, restricting the growth of phototrophic organisms, and large allochthonous organic carbon inputs from the catchment area, the lake metabolism is expected to be dominated by heterotrophic organisms. In this study, we test this assumption and show that the potential for autotrophic carbon fixation and internal carbon cycling is high throughout the water column. Further, we show that during the summer stratification carbon fixation can exceed respiration in a boreal lake even below the euphotic zone. Metagenome-assembled genomes and 16S profiling of a vertical transect of the lake revealed multiple organisms in an oxygen-depleted compartment belonging to novel or poorly characterized phyla. Many of these organisms were chemolithotrophic, potentially deriving their energy from reactions related to sulfur, iron, and nitrogen transformations. The community, as well as the functions, was stratified along the redox gradient. The autotrophic potential in the lake metagenome below the oxygenic zone was high, pointing toward a need for revising our concepts of internal carbon cycling in boreal lakes. Further, the importance of chemolithoautotrophy for the internal carbon cycling suggests that many predicted climate change-associated fluctuations in the physical properties of the lake, such as altered mixing patterns, likely have consequences for the whole-lake metabolism even beyond the impact to the phototrophic community. IMPORTANCE Autotrophic organisms at the base of the food web are the only life form capable of turning inorganic carbon into the organic form, facilitating the survival of all other organisms. In certain environments, the autotrophic production is limited by environmental conditions and the food web is supported by external carbon inputs. One such environment is stratified boreal lakes, which are one of the biggest natural sources of greenhouse gas emissions in the boreal region. Thus, carbon cycling in these habitats is of utmost importance for the future climate. Here, we demonstrate a high potential for internal carbon cycling via phototrophic and novel chemolithotrophic organisms in the anoxic, poorly illuminated layers of a boreal lake. Our results significantly increase our knowledge on the microbial communities and their metabolic potential in oxygen-depleted freshwaters and help to understand and predict how climate change-induced alterations could impact the lake carbon dynamics.

ABSTRACT Genome streamlining is frequently observed in free-living aquatic microorganisms and results in physiological dependencies between microorganisms. However, we know little about the specificity of these microbial associations. In order to examine the specificity and extent of these associations, we established mixed cultures from three different freshwater environments and analyzed the cooccurrence of organisms using a metagenomic time series. Free-living microorganisms with streamlined genomes lacking multiple biosynthetic pathways showed no clear recurring pattern in their interaction partners. Free-living freshwater bacteria form promiscuous cooperative associations. This notion contrasts with the well-documented high specificities of interaction partners in host-associated bacteria. Considering all data together, we suggest that highly abundant free-living bacterial lineages are functionally versatile in their interactions despite their distinct streamlining tendencies at the single-cell level. This metabolic versatility facilitates interactions with a variable set of community members.

Marine and brackish surface waters are highly dynamic habitats that undergo repeated seasonal variations in microbial community composition and function throughout time. While succession of the various microbial groups has been well investigated, little is known about the underlying gene-expression of the microbial community. We investigated microbial interactions via metatranscriptomics over a spring to fall seasonal cycle in the brackish Baltic Sea surface waters, a temperate brackish water ecosystem periodically promoting massive cyanobacterial blooms, which have implications for primary production, nutrient cycling, and expansion of hypoxic zones. Network analysis of the gene expression of all microbes from 0.22 to 200 μm in size and of the major taxonomic groups dissected the seasonal cycle into four components that comprised genes peaking during different periods of the bloom. Photoautotrophic nitrogen-fixing Cyanobacteria displayed the highest connectivity among the microbes, in contrast to chemoautotrophic ammonia-oxidizing Thaumarchaeota, while heterotrophs dominated connectivity among pre- and post-bloom peaking genes. The network was also composed of distinct functional connectivities, with an early season balance between carbon metabolism and ATP synthesis shifting to a dominance of ATP synthesis during the bloom, while carbon degradation, specifically through the glyoxylate shunt, characterized the post-bloom period, driven by Alphaproteobacteria as well as by Gammaproteobacteria of the SAR86 and SAR92 clusters. Our study stresses the exceptionally strong biotic driving force executed by cyanobacterial blooms on associated microbial communities in the Baltic Sea and highlights the impact cyanobacterial blooms have on functional microbial community composition.

ABSTRACT Syntrophy among Archaea and Bacteria facilitates the anaerobic degradation of organic compounds to CH4 and CO2. Particularly during aliphatic and aromatic hydrocarbon mineralization, as in the case of crude oil reservoirs and petroleum-contaminated sediments, metabolic interactions between obligate mutualistic microbial partners are of central importance. Using micromanipulation combined with shotgun metagenomic approaches, we describe the genomes of complex consortia within short-chain alkane-degrading cultures operating under methanogenic conditions. Metabolic reconstruction revealed that only a small fraction of genes in the metagenome-assembled genomes encode the capacity for fermentation of alkanes facilitated by energy conservation linked to H2 metabolism. Instead, the presence of inferred lifestyles based on scavenging anabolic products and intermediate fermentation products derived from detrital biomass was a common feature. Additionally, inferred auxotrophy for vitamins and amino acids suggests that the hydrocarbon-degrading microbial assemblages are structured and maintained by multiple interactions beyond the canonical H2-producing and syntrophic alkane degrader-methanogen partnership. Compared to previous work, our report points to a higher order of complexity in microbial consortia engaged in anaerobic hydrocarbon transformation. IMPORTANCE Microbial interactions between Archaea and Bacteria mediate many important chemical transformations in the biosphere from degrading abundant polymers to synthesis of toxic compounds. Two of the most pressing issues in microbial interactions are how consortia are established and how we can modulate these microbial communities to express desirable functions. Here, we propose that public goods (i.e., metabolites of high energy demand in biosynthesis) facilitate energy conservation for life under energy-limited conditions and determine the assembly and function of the consortia. Our report suggests that an understanding of public good dynamics could result in new ways to improve microbial pollutant degradation in anaerobic systems.

Understanding the distribution of taxa and associated traits across different environments is one of the central questions in microbial ecology. High-throughput sequencing (HTS) studies are presently generating huge volumes of data to address this biogeographical topic. However, these studies are often focused on specific environment types or processes leading to the production of individual, unconnected datasets. The large amounts of legacy sequence data with associated metadata that exist can be harnessed to better place the genetic information found in these surveys into a wider environmental context. Here we introduce a software program, seqenv, to carry out precisely such a task. It automatically performs similarity searches of short sequences against the “nt” nucleotide database provided by NCBI and, out of every hit, extracts–if it is available–the textual metadata field. After collecting all the isolation sources from all the search results, we run a text mining algorithm to identify and parse words that are associated with the Environmental Ontology (EnvO) controlled vocabulary. This, in turn, enables us to determine both in which environments individual sequences or taxa have previously been observed and, by weighted summation of those results, to summarize complete samples. We present two demonstrative applications of seqenv to a survey of ammonia oxidizing archaea as well as to a plankton paleome dataset from the Black Sea. These demonstrate the ability of the tool to reveal novel patterns in HTS and its utility in the fields of environmental source tracking, paleontology, and studies of microbial biogeography. To install seqenv, go to: https://github.com/xapple/seqenv.

As new sequencing technologies become cheaper and older ones disappear, laboratories switch vendors and platforms. Validating the new setups is a crucial part of conducting rigorous scientific research. Here we report on the reliability and biases of performing bacterial 16S rRNA gene amplicon paired-end sequencing on the MiSeq Illumina platform. We designed a protocol using 50 barcode pairs to run samples in parallel and coded a pipeline to process the data. Sequencing the same sediment sample in 248 replicates as well as 70 samples from alkaline soda lakes, we evaluated the performance of the method with regards to estimates of alpha and beta diversity. Using different purification and DNA quantification procedures we always found up to 5-fold differences in the yield of sequences between individually barcodes samples. Using either a one-step or a two-step PCR preparation resulted in significantly different estimates in both alpha and beta diversity. Comparing with a previous method based on 454 pyrosequencing, we found that our Illumina protocol performed in a similar manner - with the exception for evenness estimates where correspondence between the methods was low. We further quantified the data loss at every processing step eventually accumulating to 50% of the raw reads. When evaluating different OTU clustering methods, we observed a stark contrast between the results of QIIME with default settings and the more recent UPARSE algorithm when it comes to the number of OTUs generated. Still, overall trends in alpha and beta diversity corresponded highly using both clustering methods. Our procedure performed well considering the precisions of alpha and beta diversity estimates, with insignificant effects of individual barcodes. Comparative analyses suggest that 454 and Illumina sequence data can be combined if the same PCR protocol and bioinformatic workflows are used for describing patterns in richness, beta-diversity and taxonomic composition.

Many freshwater phytoplankton species have the potential to form transient nuisance blooms that affect water quality and other aquatic biota. Heterotrophic bacteria can influence such blooms via nutrient regeneration but also via antagonism and other biotic interactions. We studied the composition of bacterial communities associated with three bloom-forming freshwater phytoplankton species, the diatom Aulacoseira granulata and the cyanobacteria Microcystis aeruginosa and Cylindrospermopsis raciborskii. Experimental cultures incubated with and without lake bacteria were sampled in three different growth phases and bacterial community composition was assessed by 454-Pyrosequencing of 16S rRNA gene amplicons. Betaproteobacteria were dominant in all cultures inoculated with lake bacteria, but decreased during the experiment. In contrast, Alphaproteobacteria, which made up the second most abundant class of bacteria, increased overall during the course of the experiment. Other bacterial classes responded in contrasting ways to the experimental incubations causing significantly different bacterial communities to develop in response to host phytoplankton species, growth phase and between attached and free-living fractions. Differences in bacterial community composition between cyanobacteria and diatom cultures were greater than between the two cyanobacteria. Despite the significance, major differences between phytoplankton cultures were in the proportion of the OTUs rather than in the absence or presence of specific taxa. Different phytoplankton species favoring different bacterial communities may have important consequences for the fate of organic matter in systems where these bloom forming species occur. The dynamics and development of transient blooms may also be affected as bacterial communities seem to influence phytoplankton species growth in contrasting ways.

In order to characterize how disturbances to microbial communities are propagated over temporal and spatial scales in aquatic environments, the dynamics of bacterial assemblages throughout a subtropical coastal embayment were investigated via SSU rRNA gene analyses over an 8-month period, which encompassed a large storm event. During non-perturbed conditions, sampling sites clustered into three groups based on their microbial community composition: an offshore oceanic group, a freshwater group, and a distinct and persistent coastal group. Significant differences in measured environmental parameters or in the bacterial community due to the storm event were found only within the coastal cluster of sampling sites, and only at 5 of 12 locations; three of these sites showed a significant response in both environmental and bacterial community characteristics. These responses were most pronounced at sites close to the shoreline. During the storm event, otherwise common bacterioplankton community members such as marine Synechococcus sp. and members of the SAR11 clade of Alphaproteobacteria decreased in relative abundance in the affected coastal zone, whereas several lineages of Gammaproteobacteria, Betaproteobacteria, and members of the Roseobacter clade of Alphaproteobacteria increased. The complex spatial patterns in both environmental conditions and microbial community structure related to freshwater runoff and wind convection during the perturbation event leads us to conclude that spatial heterogeneity was an important factor influencing both the dynamics and the resistance of the bacterioplankton communities to disturbances throughout this complex subtropical coastal system. This heterogeneity may play a role in facilitating a rapid rebound of regions harboring distinctly coastal bacterioplankton communities to their pre-disturbed taxonomic composition.

The recognition and discrimination of phytoplankton species is one of the foundations of freshwater biodiversity research and environmental monitoring. This step is frequently a bottleneck in the analytical chain from sampling to data analysis and subsequent environmental status evaluation. Here we present phytoplankton diversity data from 49 lakes including three seasonal surveys assessed by next generation sequencing (NGS) of 16S ribosomal RNA chloroplast and cyanobacterial gene amplicons and also compare part of these datasets with identification based on morphology. Direct comparison of NGS to microscopic data from three time-series showed that NGS was able to capture the seasonality in phytoplankton succession as observed by microscopy. Still, the PCR-based approach was only semi-quantitative, and detailed NGS and microscopy taxa lists had only low taxonomic correspondence. This is probably due to, both, methodological constraints and current discrepancies in taxonomic frameworks. Discrepancies included Euglenophyta and Heterokonta that were scarce in the NGS but frequently detected by microscopy and Cyanobacteria that were in general more abundant and classified with high resolution by NGS. A deep-branching taxonomically unclassified cluster was frequently detected by NGS but could not be linked to any group identified by microscopy. NGS derived phytoplankton composition differed significantly among lakes with different trophic status, showing that our approach can resolve phytoplankton communities at a level relevant for ecosystem management. The high reproducibility and potential for standardization and parallelization makes our NGS approach an excellent candidate for simultaneous monitoring of prokaryotic and eukaryotic phytoplankton in inland waters.

Nutrient limitation and resource competition in bacterial and phytoplankton communities may appear different when considering different levels of taxonomic resolution. Nutrient amendment experiments conducted in a boreal lake on three occasions during one open water season revealed complex responses in overall bacterioplankton and phytoplankton abundance and biovolume. In general, bacteria were dominant in spring, while phytoplankton was clearly the predominant group in autumn. Seasonal differences in the community composition of bacteria and phytoplankton were mainly related to changes in observed taxa, while the differences across nutrient treatments within an experiment were due to changes in relative contributions of certain higher- and lower-level phylogenetic groups. Of the main bacterioplankton phyla, only Actinobacteria had a treatment response that was visible even at the phylum level throughout the season. With increasing resolution (from 75 to 99% sequence similarity) major responses to nutrient amendments appeared using 454 pyrosequencing data of 16S rRNA amplicons. This further revealed that OTUs (defined by 97% sequence similarity) annotated to the same highly resolved freshwater groups appeared to occur during different seasons and were showing treatment-dependent differentiation, indicating that OTUs within these groups were not ecologically coherent. Similarly, phytoplankton species from the same genera responded differently to nutrient amendments even though biovolumes of the majority of taxa increased when both nitrogen and phosphorus were added simultaneously. The bacterioplankton and phytoplankton community compositions showed concurrent trajectories that could be seen in synchronous succession patterns over the season. Overall, our data revealed that the response of both communities to nutrient changes was highly dependent on season and that contradictory results may be obtained when using different taxonomic resolutions.