Vertical partitioning of freshwater bacterioplankton community in a deep mesotrophic lake with a fully oxygenated hypolimnion (Lake Biwa, Japan)

Spatio-temporal changes of small protist and free-living bacterial communities in a temperate dimictic lake: insights from metabarcoding and machine learning

Microbial communities, which include prokaryotes and protists, play an important role in aquatic ecosystems and influence ecological processes. To understand these communities, metabarcoding provides a powerful tool to assess their taxonomic composition and track spatio-temporal dynamics in both marine and freshwater environments. While marine ecosystems have been extensively studied, there is a notable research gap in understanding eukaryotic microbial communities in temperate lakes. Our study addresses this gap by investigating the free-living bacteria and small protist communities in Lake Roś (Poland), a dimictic temperate lake. Metabarcoding analysis revealed that both the bacterial and protist communities exhibit distinct seasonal patterns that are not necessarily shaped by dominant taxa. Furthermore, machine learning and statistical methods identified crucial amplicon sequence variants (ASVs) specific to each season. In addition, we identified a distinct community in the anoxic hypolimnion. We have also shown that the key factors shaping the composition of analysed community are temperature, oxygen, and silicon concentration. Understanding these community structures and the underlying factors is important in the context of climate change potentially impacting mixing patterns and leading to prolonged stratification.

Out of sight, but not out of season: Nitrifier distributions and population dynamics in a large oligotrophic lake

Nitrification is an important control on the form and distribution of nitrogen in freshwater ecosystems. However, the seasonality of nitrogen pools and the diversity of organisms catalyzing this process have not been well documented in oligotrophic lakes. Here, we show that nitrogen pools and nitrifying organisms in Flathead Lake are temporally and vertically dynamic, with nitrifiers displaying specific preferences depending on the season. While the ammonia-oxidizing bacteria (AOB) Nitrosomonadaceae and nitrite-oxidizing bacteria (NOB) Nitrotoga dominate at depth in the summer, the ammonia-oxidizing archaea (AOA) Nitrososphaerota and NOB Nitrospirota become abundant in the winter. Given clear seasonality in ammonium, with higher concentrations during the summer, we hypothesize that the succession between these two nitrifying groups may be due to nitrogen affinity, with AOB more competitive when ammonia concentrations are higher and AOA when they are lower. Nitrifiers in Flathead Lake share more than 99% average nucleotide identity with those reported in other North American lakes but are distinct from those in Europe and Asia, indicating a role for geographic isolation as a factor controlling speciation among nitrifiers. Our study shows there are seasonal shifts in nitrogen pools and nitrifying populations, highlighting the dynamic spatial and temporal nature of nitrogen cycling in freshwater ecosystems.

Aqueous Geochemical Controls on the Sestonic Microbial Community in Lakes Michigan and Superior

Despite being the largest freshwater lake system in the world, relatively little is known about the sestonic microbial community structure in the Laurentian Great Lakes. The goal of this research was to better understand this ecosystem using high-throughput sequencing of microbial communities as a function of water depth at six locations in the westernmost Great Lakes of Superior and Michigan. The water column was characterized by gradients in temperature, dissolved oxygen (DO), pH, and other physicochemical parameters with depth. Mean nitrate concentrations were 32 μmol/L, with only slight variation within and between the lakes, and with depth. Mean available phosphorus was 0.07 μmol/L, resulting in relatively large N:P ratios (97:1) indicative of P limitation. Abundances of the phyla Actinobacteria, Bacteroidetes, Cyanobacteria, Thaumarchaeota, and Verrucomicrobia differed significantly among the Lakes. Candidatus Nitrosopumilus was present in greater abundance in Lake Superior compared to Lake Michigan, suggesting the importance of ammonia-oxidating archaea in water column N cycling in Lake Superior. The Shannon diversity index was negatively correlated with pH, temperature, and salinity, and positively correlated with DO, latitude, and N₂ saturation. Results of this study suggest that DO, pH, temperature, and salinity were major drivers shaping the community composition in the Great Lakes.

Long-Read-Resolved, Ecosystem-Wide Exploration of Nucleotide and Structural Microdiversity of Lake Bacterioplankton Genomes

Reconstruction of metagenome-assembled genomes (MAGs) has become a fundamental approach in microbial ecology. However, a MAG is hardly complete and overlooks genomic microdiversity because metagenomic assembly fails to resolve microvariants among closely related genotypes. Aiming at understanding the universal factors that drive or constrain prokaryotic genome diversification, we performed an ecosystem-wide high-resolution metagenomic exploration of microdiversity by combining spatiotemporal (2 depths × 12 months) sampling from a pelagic freshwater system, high-quality MAG reconstruction using long- and short-read metagenomic sequences, and profiling of single nucleotide variants (SNVs) and structural variants (SVs) through mapping of short and long reads to the MAGs, respectively. We reconstructed 575 MAGs, including 29 circular assemblies, providing high-quality reference genomes of freshwater bacterioplankton. Read mapping against these MAGs identified 100 to 101,781 SNVs/Mb and 0 to 305 insertions, 0 to 467 deletions, 0 to 41 duplications, and 0 to 6 inversions for each MAG. Nonsynonymous SNVs were accumulated in genes potentially involved in cell surface structural modification to evade phage recognition. Most (80.2%) deletions overlapped with a gene coding region, and genes of prokaryotic defense systems were most frequently (>8% of the genes) overlapped with a deletion. Some such deletions exhibited a monthly shift in their allele frequency, suggesting a rapid turnover of genotypes in response to phage predation. MAGs with extremely low microdiversity were either rare or opportunistic bloomers, suggesting that population persistency is key to their genomic diversification. The results concluded that prokaryotic genomic diversification is driven primarily by viral load and constrained by a population bottleneck. IMPORTANCE Identifying intraspecies genomic diversity (microdiversity) is crucial to understanding microbial ecology and evolution. However, microdiversity among environmental assemblages is not well investigated, because most microbes are difficult to culture. In this study, we performed cultivation-independent exploration of bacterial genomic microdiversity in a lake ecosystem using a combination of short- and long-read metagenomic analyses. The results revealed the broad spectrum of genomic microdiversity among the diverse bacterial species in the ecosystem, which has been overlooked by conventional approaches. Our ecosystem-wide exploration further allowed comparative analysis among the genomes and genes and revealed factors behind microbial genomic diversification, namely, that diversification is driven primarily by resistance against viral infection and constrained by the population size.

Quantification of archaea-driven freshwater nitrification from single cell to ecosystem levels

Deep oligotrophic lakes sustain large populations of the class Nitrososphaeria (Thaumarchaeota) in their hypolimnion. They are thought to be the key ammonia oxidizers in this habitat, but their impact on N-cycling in lakes has rarely been quantified. We followed this archaeal population in one of Europe's largest lakes, Lake Constance, for two consecutive years using metagenomics and metatranscriptomics combined with stable isotope-based activity measurements. An abundant (8-39% of picoplankton) and transcriptionally active archaeal ecotype dominated the nitrifying community. It represented a freshwater-specific species present in major inland water bodies, for which we propose the name "Candidatus Nitrosopumilus limneticus". Its biomass corresponded to 12% of carbon stored in phytoplankton over the year´s cycle. Ca. N. limneticus populations incorporated significantly more ammonium than most other microorganisms in the hypolimnion and were driving potential ammonia oxidation rates of 6.0 ± 0.9 nmol l^‒1 d^‒1, corresponding to potential cell-specific rates of 0.21 ± 0.11 fmol cell^-1 d^-1. At the ecosystem level, this translates to a maximum capacity of archaea-driven nitrification of 1.76 × 10⁹ g N-ammonia per year or 11% of N-biomass produced annually by phytoplankton. We show that ammonia-oxidizing archaea play an equally important role in the nitrogen cycle of deep oligotrophic lakes as their counterparts in marine ecosystems.

Local Habitat Filtering Shapes Microbial Community Structure in Four Closely Spaced Lakes in the High Arctic

Arctic lakes are experiencing increasingly shorter periods of ice cover due to accelerated warming at northern high latitudes. Given the control of ice cover thickness and duration over many limnological processes, these changes will have pervasive effects. However, due to their remote and extreme locations even first-order data on lake ecology is lacking for many ecosystems. The aim of this study was to characterize and compare the microbial communities of four closely spaced lakes in Stuckberry Valley (northern Ellesmere Island, Canadian Arctic Archipelago), in the coastal margin zone of the Last Ice Area, that differed in their physicochemical, morphological and catchment characteristics. We performed high-throughput amplicon sequencing of the V4 16S rRNA gene to provide inter- and intra-lake comparisons. Two deep (>25 m) and mostly oxygenated lakes showed highly similar community assemblages that were distinct from those of two shallower lakes (<10 m) with anoxic bottom waters. Proteobacteria, Verrucomicrobia, and Planctomycetes were the major phyla present in the four water bodies. One deep lake contained elevated proportions of Cyanobacteria and Thaumarchaeota that distinguished it from the others, while the shallow lakes had abundant communities of predatory bacteria, as well as microbes in their bottom waters that contribute to sulfur and methane cycles. Despite their proximity, our data suggest that local habitat filtering is the primary determinant of microbial diversity in these systems. This study provides the first detailed examination of the microbial assemblages of the Stuckberry lakes system, resulting in new insights into the microbial ecology of the High Arctic.

Species-Level Spatio-Temporal Dynamics of Cyanobacteria in a Hard-Water Temperate Lake in the Southern Baltics

Cyanobacteria are important primary producers in temperate freshwater ecosystems. However, studies on the seasonal and spatial distribution of cyanobacteria in deep lakes based on high-throughput DNA sequencing are still rare. In this study, we combined monthly water sampling and monitoring in 2019, amplicon sequence variants analysis (ASVs; a proxy for different species) and quantitative PCR targeting overall cyanobacteria abundance to describe the seasonal and spatial dynamics of cyanobacteria in the deep hard-water oligo-mesotrophic Lake Tiefer See, NE Germany. We observed significant seasonal variation in the cyanobacterial community composition (p < 0.05) in the epi- and metalimnion layers, but not in the hypolimnion. In winter-when the water column is mixed-picocyanobacteria (Synechococcus and Cyanobium) were dominant. With the onset of stratification in late spring, we observed potential niche specialization and coexistence among the cyanobacteria taxa driven mainly by light and nutrient dynamics. Specifically, ASVs assigned to picocyanobacteria and the genus Planktothrix were the main contributors to the formation of deep chlorophyll maxima along a light gradient. While Synechococcus and different Cyanobium ASVs were abundant in the epilimnion up to the base of the euphotic zone from spring to fall, Planktothrix mainly occurred in the metalimnetic layer below the euphotic zone where also overall cyanobacteria abundance was highest in summer. Our data revealed two potentially psychrotolerant (cold-adapted) Cyanobium species that appear to cope well under conditions of lower hypolimnetic water temperature and light as well as increasing sediment-released phosphate in the deeper waters in summer. The potential cold-adapted Cyanobium species were also dominant throughout the water column in fall and winter. Furthermore, Snowella and Microcystis-related ASVs were abundant in the water column during the onset of fall turnover. Altogether, these findings suggest previously unascertained and considerable spatiotemporal changes in the community of cyanobacteria on the species level especially within the genus Cyanobium in deep hard-water temperate lakes.

Microdiversity and phylogeographic diversification of bacterioplankton in pelagic freshwater systems revealed through long-read amplicon sequencing

BACKGROUND

Freshwater ecosystems are inhabited by members of cosmopolitan bacterioplankton lineages despite the disconnected nature of these habitats. The lineages are delineated based on > 97% 16S rRNA gene sequence similarity, but their intra-lineage microdiversity and phylogeography, which are key to understanding the eco-evolutional processes behind their ubiquity, remain unresolved. Here, we applied long-read amplicon sequencing targeting nearly full-length 16S rRNA genes and the adjacent ribosomal internal transcribed spacer sequences to reveal the intra-lineage diversities of pelagic bacterioplankton assemblages in 11 deep freshwater lakes in Japan and Europe.

RESULTS

Our single nucleotide-resolved analysis, which was validated using shotgun metagenomic sequencing, uncovered 7-101 amplicon sequence variants for each of the 11 predominant bacterial lineages and demonstrated sympatric, allopatric, and temporal microdiversities that could not be resolved through conventional approaches. Clusters of samples with similar intra-lineage population compositions were identified, which consistently supported genetic isolation between Japan and Europe. At a regional scale (up to hundreds of kilometers), dispersal between lakes was unlikely to be a limiting factor, and environmental factors or genetic drift were potential determinants of population composition. The extent of microdiversification varied among lineages, suggesting that highly diversified lineages (e.g., Iluma-A2 and acI-A1) achieve their ubiquity by containing a consortium of genotypes specific to each habitat, while less diversified lineages (e.g., CL500-11) may be ubiquitous due to a small number of widespread genotypes. The lowest extent of intra-lineage diversification was observed among the dominant hypolimnion-specific lineage (CL500-11), suggesting that their dispersal among lakes is not limited despite the hypolimnion being a more isolated habitat than the epilimnion.

CONCLUSIONS

Our novel approach complemented the limited resolution of short-read amplicon sequencing and limited sensitivity of the metagenome assembly-based approach, and highlighted the complex ecological processes underlying the ubiquity of freshwater bacterioplankton lineages. To fully exploit the performance of the method, its relatively low read throughput is the major bottleneck to be overcome in the future. Video abstract.

Reservoir water stratification and mixing affects microbial community structure and functional community composition in a stratified drinking reservoir

To investigate how the aquatic bacterial community of a stratified reservoir drives the evolution of water parameters, the microbial community structure and network characteristics of bacteria in a stratified reservoir were investigated using Illumina MiSeq sequencing technology. A total of 42 phyla and 689 distinct genera were identified, which showed significant seasonal variation. Additionally, stratified variations in the bacterial community strongly reflected the vertical gradient and seasonal changes in water temperature, dissolved oxygen, and nutrition concentration. Furthermore, principal coordinate analysis indicated that most microorganisms were likely influenced by changes in water stratification conditions, exhibiting significant differences during the stratification period and mixing period based on Adonis, MRPP, and Anosim. Compared to the stratification period, 123 enhanced operational taxonomic units (OTUs; 29%) and 226 depleted OTUs (52%) were identified during the mixing period. Linear discriminant analysis effect size results showed that 15 major genera were enriched in the mixing period and 10 major genera were enriched in the stratification period. Importantly, network analysis revealed that the keystone species belonged to hgcI_clade, CL500-29, Acidibacter, Paucimonas, Flavobacterium, Prochlorothrix, Xanthomonadales, Chloroflexia, Burkholderiales, OPB56, KI89A_clade, Synechococcus, Caulobacter or were unclassified. Redundancy analysis showed that temperature, dissolved oxygen, pH, chlorophyll-α, total phosphorus, nitrate, and ammonia were important factors influencing the water bacterial community and function composition, which were consistent with the results of the Mantel test analysis. Furthermore, random forest analysis showed that temperature, dissolved oxygen, ammonia, and total dissolved phosphorous were the most important variables predicting water bacterial community and function community α- and β-diversity (P < 0.05). Overall, these results provide insight into the interactions between the microbial community and water quality evolution mechanism in Zhoucun reservoir.

Planctomycetes as a Vital Constituent of the Microbial Communities Inhabiting Different Layers of the Meromictic Lake Sælenvannet (Norway)

Meromictic lakes are permanently stratified lakes that display steep gradients in salinity, oxygen and sulphur compounds tightly linked to bacterial community structure and diversity. Lake Sælenvannet is a meromictic lake located south of Bergen, Norway. The 26 m deep lake is connected to the open sea and permanently stratified into two layers separated by a chemocline. The upper water layer is brackish with major input from water runoff from the surroundings. The bottom layer consists of old saline water with low or no oxygen concentrations. Bacteria from phylum Planctomycetes are reported to be ubiquitous in lake environments. They are involved in the degradation of complex carbon sources in aquatic environments and are also linked to anaerobic processes such as fermentation and sulphur reduction. To study Planctomycete distribution along a chemical gradient, we sampled the water column throughout Lake Sælenvannet in 2012 and profiled the microbial community using 16S rRNA amplicon sequencing (metabarcoding) with 454 pyrosequencing. Planctomycetes related 16S rRNA gene sequences were found to be present both in the oxic and anoxic parts of the lake and showed an uneven distribution throughout the water column, with the highest relative abundance of 10% found in the saline anoxic layer at 15 m depth. In a follow-up study in 2014, samples from eight different depths were collected for enrichment and isolation of novel Planctomycetes. This study resulted in successful isolation in pure culture of 10 isolates affiliated to four different genera from the family Planctomycetaceae. One strain closely related to Blastopirellula cremea was isolated from 9 m depth, and two novel strains affiliated to the genera Stieleria and Gimesia were isolated at 7 and 9 m depths, respectively. Furthermore, seven isolates with identical 16S rRNA gene sequences were retrieved from seven different depths which varied greatly in salinity and chemical composition. These isolates likely represent a new species affiliated to Rubinisphaera. The adaptation of this novel Planctomycete to water depths spanning the entire chemical gradient could indicate a high phenotypic plasticity and/or a very efficient survival strategy. Overall, our results show the presence of a diverse group of Planctomycetes in Lake Sælenvannet, with a strong potential for novel adaptations to chemical stress factors.

Microbial communities of the Laurentian Great Lakes reflect connectivity and local biogeochemistry

The Laurentian Great Lakes are a vast, interconnected freshwater system spanning strong physicochemical gradients, thus constituting a powerful natural laboratory for addressing fundamental questions about microbial ecology and evolution. We present a comparative analysis of pelagic microbial communities across all five Laurentian Great Lakes, focusing on Bacterial and Archaeal picoplankton characterized via 16S rRNA amplicon sequencing. We collected samples throughout the water column from the major basins of each lake in spring and summer over 2 years. Two oligotypes, classified as LD12 (Alphaproteobacteria) and acI‐B1 (Actinobacteria), were among the most abundant in every sample. At the same time, microbial communities showed distinct patterns with depth during summer stratification. Deep hypolimnion samples were frequently dominated by a Chloroflexi oligotype that reached up to 19% relative abundance. Stratified surface communities differed between the colder, less productive upper lakes (Superior, Michigan, Huron) and warmer, more productive lower lakes (Erie, Ontario), in part due to an Actinobacteria oligotype (acI‐C2) that averaged 7.7% of sequences in the lower lakes but <0.2% in the upper lakes. Together, our findings suggest that both hydrologic connectivity and local selective pressures shape microbial communities in the Great Lakes and establish a framework for future investigations.

A single Thaumarchaeon drives nitrification in deep oligotrophic Lake Constance

Ammonia released during organic matter mineralization is converted during nitrification to nitrate. We followed spatiotemporal dynamics of the nitrifying microbial community in deep oligotrophic Lake Constance. Depth-dependent decrease of total ammonium (0.01-0.84 μM) indicated the hypolimnion as the major place of nitrification with ¹⁵ N-isotope dilution measurements indicating a threefold daily turnover of hypolimnetic total ammonium. This was mirrored by a strong increase of ammonia-oxidizing Thaumarchaeota towards the hypolimnion (13%-21% of bacterioplankton) throughout spring to autumn as revealed by amplicon sequencing and quantitative polymerase chain reaction. Ammonia-oxidizing bacteria were typically two orders of magnitude less abundant and completely ammonia-oxidizing (comammox) bacteria were not detected. Both, 16S rRNA gene and amoA (encoding ammonia monooxygenase subunit B) analyses identified only one major species-level operational taxonomic unit (OTU) of Thaumarchaeota (99% of all ammonia oxidizers in the hypolimnion), which was affiliated to Nitrosopumilus spp. The relative abundance distribution of the single Thaumarchaeon strongly correlated to an equally abundant Chloroflexi clade CL500-11 OTU and a Nitrospira OTU that was one order of magnitude less abundant. The latter dominated among recognized nitrite oxidizers. This extremely low diversity of nitrifiers shows how vulnerable the ecosystem process of nitrification may be in Lake Constance as Central Europe's third largest lake.

Genome-resolved viral and cellular metagenomes revealed potential key virus-host interactions in a deep freshwater lake

Metagenomics has dramatically expanded the known virosphere, but freshwater viral diversity and their ecological interaction with hosts remain poorly understood. Here, we conducted a metagenomic exploration of planktonic dsDNA prokaryotic viruses by sequencing both virion (<0.22 μm) and cellular (0.22-5.0 μm) fractions collected spatiotemporally from a deep freshwater lake (Lake Biwa, Japan). This simultaneously reconstructed 183 complete (i.e., circular) viral genomes and 57 bacterioplankton metagenome-assembled genomes. Analysis of metagenomic read coverage revealed vertical partitioning of the viral community analogous to the vertically stratified bacterioplankton community. The hypolimnetic community was generally stable during stratification, but occasionally shifted abruptly, presumably due to lysogenic induction. Genes involved in assimilatory sulfate reduction were encoded in 20 (10.9%) viral genomes, including those of dominant viruses, and may aid viral propagation in sulfur-limited freshwater systems. Hosts were predicted for 40 (21.9%) viral genomes, encompassing 10 phyla (or classes of Proteobacteria) including ubiquitous freshwater bacterioplankton lineages (e.g., Ca. Fonsibacter and Ca. Nitrosoarchaeum). Comparison with viral genomes derived from published metagenomes revealed viral phylogeographic connectivity in geographically isolated habitats. Notably, analogous to their hosts, actinobacterial viruses were among the most diverse, ubiquitous and abundant viral groups in freshwater systems, with potential high lytic activity in surface waters.

Evolution in action: habitat transition from sediment to the pelagial leads to genome streamlining in Methylophilaceae

The most abundant aquatic microbes are small in cell and genome size. Genome-streamlining theory predicts gene loss caused by evolutionary selection driven by environmental factors, favouring superior competitors for limiting resources. However, evolutionary histories of such abundant, genome-streamlined microbes remain largely unknown. Here we reconstruct the series of steps in the evolution of some of the most abundant genome-streamlined microbes in freshwaters ("Ca. Methylopumilus") and oceans (marine lineage OM43). A broad genomic spectrum is visible in the family Methylophilaceae (Betaproteobacteria), from sediment microbes with medium-sized genomes (2-3 Mbp genome size), an occasionally blooming pelagic intermediate (1.7 Mbp), and the most reduced pelagic forms (1.3 Mbp). We show that a habitat transition from freshwater sediment to the relatively oligotrophic pelagial was accompanied by progressive gene loss and adaptive gains. Gene loss has mainly affected functions not necessarily required or advantageous in the pelagial or is encoded by redundant pathways. Likewise, we identified genes providing adaptations to oligotrophic conditions that have been transmitted horizontally from pelagic freshwater microbes. Remarkably, the secondary transition from the pelagial of lakes to the oceans required only slight modifications, i.e., adaptations to higher salinity, gained via horizontal gene transfer from indigenous microbes. Our study provides first genomic evidence of genome reduction taking place during habitat transitions. In this regard, the family Methylophilaceae is an exceptional model for tracing the evolutionary history of genome streamlining as such a collection of evolutionarily related microbes from different habitats is rare in the microbial world.

Widespread Dominance of Kinetoplastids and Unexpected Presence of Diplonemids in Deep Freshwater Lakes

Kinetoplastid flagellates are generally abundant in the deep sea and recently they were even found to be dominant in the hypolimnion of a deep freshwater lake. Therefore, to understand the distribution of kinetoplastids in deep freshwater lakes, we have collected vertical samples from five lakes in Japan. The abundance of kinetoplastids was enumerated by Catalyzed Reporter Deposition-Fluorescence in situ Hybridization, and the diversity was determined by 18S amplicon sequencing using universal eukaryote and kinetoplastid-specific primers. Kinetoplastids were abundant in the deep waters of all the lakes, contributing up to 53.6% of total nanoeukaryotes. Despite this significant contribution, kinetoplastids remain undetected by amplicon sequencing using universal primers that are widely used in eukaryotic diversity studies. However, they were detected with specific primers, and the communities were characterized by both ubiquitous and lake-specific unique OTUs. Oligotyping of a ubiquitous and dominant OTU revealed the presence of lake-specific sequence types (oligotypes). Remarkably, we also detected diplonemids (a sister group of kinetoplastids and considered to be specific in the marine habitat) using kinetoplastid-specific primers, showing their presence in freshwaters. Underestimation of kinetoplastids and diplonemids using universal primers indicates that euglenozoan flagellates are overlooked in diversity studies worldwide. The present study highlighted the importance of kinetoplastids in the hypolimnion of deep lakes, thereby indicating their role in material cycling in deep waters.

Effects of Habitat Partitioning on the Distribution of Bacterioplankton in Deep Lakes

In deep lakes, many investigations highlighted the existence of exclusive groups of bacteria adapted to deep oxygenated and hypoxic and anoxic hypolimnia. Nevertheless, the extent of bacterial strain diversity has been much less scrutinized. This aspect is essential for an unbiased estimation of genetic variation, biodiversity, and population structure, which are essential for studying important research questions such as biogeographical patterns, temporal and spatial variability and the environmental factors affecting this variability. This study investigated the bacterioplankton community in the epilimnetic layers and in the oxygenated and hypoxic/anoxic hypolimnia of five large and deep lakes located at the southern border of the Alps using high throughput sequencing (HTS) analyses (16S rDNA) and identification of amplicon sequence variants (ASVs) resolving reads differing by as little as one nucleotide. The study sites, which included two oligomictic (Garda and Como) and three meromictic lakes (Iseo, Lugano, and Idro) with maximum depths spanning from 124 to 410 m, were chosen among large lakes to represent an oxic-hypoxic gradient. The analyses showed the existence of several unique ASVs in the three layers of the five lakes. In the case of cyanobacteria, this confirmed previous analyses made at the level of strains or based on oligotyping methods. As expected, the communities in the hypoxic/anoxic monimolimnia showed a strong differentiation from the oxygenated layer, with the exclusive presence in single lakes of several unique ASVs. In the meromictic lakes, results supported the hypothesis that the formation of isolated monimolimnia sustained the development of highly diversified bacterial communities through ecological selection, leading to the establishment of distinctive biodiversity zones. The genera identified in these layers are well-known to activate a wide range of redox reactions at low O₂ conditions. As inferred from 16S rDNA data, the highly diversified and coupled processes sustained by the monimolimnetic microbiota are essential ecosystem services that enhance mineralization of organic matter and formation of reduced compounds, and also abatement of undesirable greenhouse gasses.

Culturing the ubiquitous freshwater actinobacterial acI lineage by supplying a biochemical 'helper' catalase

The actinobacterial acI lineage is among the most successful and ubiquitous freshwater bacterioplankton found on all continents, often representing more than half of all microbial cells in the lacustrine environment and constituting multiple ecotypes. However, stably growing pure cultures of the acI lineage have not been established despite various cultivation efforts based on ecological and genomic studies on the lineage, which is in contrast to the ocean from which abundant microorganisms such as Prochlorococcus, Pelagibacter, and Nitrosopumilus have been isolated. Here, we report the first two pure cultures of the acI lineage successfully maintained by supplementing the growth media with catalase. Catalase was critical for stabilizing the growth of acI strains irrespective of the genomic presence of the catalase-peroxidase (katG) gene. The two strains, representing two novel species, displayed differential phenotypes and distinct preferences for reduced sulfurs and carbohydrates, some of which were difficult to predict based on genomic information. Our results suggest that culture of previously uncultured freshwater bacteria can be facilitated by a simple catalase-supplement method and indicate that genome-based metabolic prediction can be complemented by physiological analyses.

Niche-directed evolution modulates genome architecture in freshwater Planctomycetes

Freshwater environments teem with microbes that do not have counterparts in culture collections or genetic data available in genomic repositories. Currently, our apprehension of evolutionary ecology of freshwater bacteria is hampered by the difficulty to establish organism models for the most representative clades. To circumvent the bottlenecks inherent to the cultivation-based techniques, we applied ecogenomics approaches in order to unravel the evolutionary history and the processes that drive genome architecture in hallmark freshwater lineages from the phylum Planctomycetes. The evolutionary history inferences showed that sediment/soil Planctomycetes transitioned to aquatic environments, where they gave rise to new freshwater-specific clades. The most abundant lineage was found to have the most specialised lifestyle (increased regulatory genetic circuits, metabolism tuned for mineralization of proteinaceous sinking aggregates, psychrotrophic behaviour) within the analysed clades and to harbour the smallest freshwater Planctomycetes genomes, highlighting a genomic architecture shaped by niche-directed evolution (through loss of functions and pathways not needed in the newly acquired freshwater niche).

Metaepigenomic analysis reveals the unexplored diversity of DNA methylation in an environmental prokaryotic community

DNA methylation plays important roles in prokaryotes, and their genomic landscapes—prokaryotic epigenomes—have recently begun to be disclosed. However, our knowledge of prokaryotic methylation systems is focused on those of culturable microbes, which are rare in nature. Here, we used single-molecule real-time and circular consensus sequencing techniques to reveal the ‘metaepigenomes’ of a microbial community in the largest lake in Japan, Lake Biwa. We reconstructed 19 draft genomes from diverse bacterial and archaeal groups, most of which are yet to be cultured. The analysis of DNA chemical modifications in those genomes revealed 22 methylated motifs, nine of which were novel. We identified methyltransferase genes likely responsible for methylation of the novel motifs, and confirmed the catalytic specificities of four of them via transformation experiments using synthetic genes. Our study highlights metaepigenomics as a powerful approach for identification of the vast unexplored variety of prokaryotic DNA methylation systems in nature.

Our knowledge of DNA methylation systems in prokaryotes is mostly limited to those of culturable microbes. Here, Hiraoka et al. analyse DNA methylation patterns in metagenomic data from a microbial community, revealing new methylated motifs and experimentally validating the methyltransferases’ specificities.

The Broad Habitat Spectrum of the CL500-11 Lineage (Phylum Chloroflexi), a Dominant Bacterioplankton in Oxygenated Hypolimnia of Deep Freshwater Lakes

CL500-11 (phylum Chloroflexi) is one of the most ubiquitous and abundant bacterioplankton lineages in deep freshwater lakes inhabiting the oxygenated hypolimnion. While metagenomics predicted possible eco-physiological characteristics of this uncultured lineage, no consensus on their ecology has so far been reached, partly because their niche is not clearly understood due to a limited number of quantitative field observations. This study investigated the abundance and distribution of CL500-11 in seven deep perialpine lakes using catalyzed reporter deposition-fluorescence in situ hybridization (CARD-FISH). Samples were taken vertically (5–12 depths in each lake) and temporally (in two lakes) at the deepest point of the lakes located in Switzerland, Italy, and Austria with varying depth, trophic state, mixing regime, and water retention time. The results showed a dominance of CL500-11 in all the lakes; their proportion to total prokaryotes ranged from 4.3% (Mondsee) to 24.3% (Lake Garda) and their abundance ranged from 0.65 × 10⁵ (Mondsee) to 1.77 × 10⁵ (Lake Garda) cells mL^-1. By summarizing available information on CL500-11 occurrence to date, we demonstrated their broad habitat spectrum, ranging from ultra-oligotrophic to meso-eutrophic lakes, while low abundances or complete absence was observed in lakes with shallow depth, low pH, and/or short water retention time (<1 year). Together with available metagenomic and geochemical evidences from literatures, here we reviewed potential substrates supporting growth of CL500-11. Overall, the present study further endorsed ubiquity and quantitative significance of CL500-11 in deep freshwater systems and narrowed the focus on their physiological characteristics and ecological importance.

Hidden in plain sight-highly abundant and diverse planktonic freshwater Chloroflexi

BACKGROUND

Representatives of the phylum Chloroflexi, though reportedly highly abundant in the extensive deep water habitats of both marine (SAR202 up to 30% of total prokaryotes) and freshwater (CL500-11 up to 26% of total prokaryotes), remain uncultivated and uncharacterized. There are few metagenomic studies on marine Chloroflexi representatives, while the pelagic freshwater Chloroflexi community is largely unknown except for a single metagenome-assembled genome of CL500-11.

RESULTS

Here, we provide the first extensive examination of the community composition of this cosmopolitan phylum in a range of pelagic habitats (176 datasets) and highlight the impact of salinity and depth on their phylogenomic composition. Reconstructed genomes (53 in total) provide a perspective on the phylogeny, metabolism, and distribution of three novel classes and two family-level taxa within the phylum Chloroflexi. We unraveled a remarkable genomic diversity of pelagic freshwater Chloroflexi representatives that thrive not only in the hypolimnion as previously suspected, but also in the epilimnion. Our results suggest that the lake hypolimnion provides a globally stable habitat reflected in lower species diversity among hypolimnion-specific CL500-11 and TK10 clusters in distantly related lakes compared to a higher species diversity of the epilimnion-specific SL56 cluster. Cell volume analyses show that the CL500-11 are among the largest prokaryotic cells in the water column of deep lakes and with a biomass to abundance ratio of two they significantly contribute to the deep lake carbon flow. Metabolic insights indicate participation of JG30-KF-CM66 representatives in the global cobalamin production via cobinamide to cobalamin salvage pathway.

CONCLUSIONS

Extending phylogenomic comparisons to brackish and marine habitats suggests salinity as the major influencer of the community composition of the deep-dwelling Chloroflexi in marine (SAR202) and freshwater (CL500-11) habitats as both counterparts thrive in intermediate brackish salinity; however, freshwater habitats harbor the most phylogenetically diverse community of pelagic Chloroflexi representatives that reside both in epi- and hypolimnion.

Synchrony of Eukaryotic and Prokaryotic Planktonic Communities in Three Seasonally Sampled Austrian Lakes

Freshwater systems are characterized by an enormous diversity of eukaryotic protists and prokaryotic taxa. The community structures in different lakes are thereby influenced by factors such as habitat size, lake chemistry, biotic interactions, and seasonality. In our study, we used high throughput 454 sequencing to study the diversity and temporal changes of prokaryotic and eukaryotic planktonic communities in three Austrian lakes during the ice-free season. In the following year, one lake was sampled again with a reduced set of sampling dates to observe reoccurring patterns. Cluster analyses (based on SSU V9 (eukaryotic) and V4 (prokaryotic) OTU composition) grouped samples according to their origin followed by separation into seasonal clusters, indicating that each lake has a unique signature based on OTU composition. These results suggest a strong habitat-specificity of microbial communities and in particular of community patterns at the OTU level. A comparison of the prokaryotic and eukaryotic datasets via co-inertia analysis (CIA) showed a consistent clustering of prokaryotic and eukaryotic samples, probably reacting to the same environmental forces (e.g., pH, conductivity). In addition, the shifts in eukaryotic and bacterioplanktonic communities generally occurred at the same time and on the same scale. Regression analyses revealed a linear relationship between an increase in Bray-Curtis dissimilarities and elapsed time. Our study shows a pronounced coupling between bacteria and eukaryotes in seasonal samplings of the three analyzed lakes. However, our temporal resolution (biweekly sampling) and data on abiotic factors were insufficient to determine if this was caused by direct biotic interactions or by reacting to the same seasonally changing environmental forces.

Succession of Bacterial Communities in a Seasonally Stratified Lake with an Anoxic and Sulfidic Hypolimnion

Although bacteria play key roles in aquatic food webs and biogeochemical cycles, information on the seasonal succession of bacterial communities in lakes is still far from complete. Here, we report results of an integrative study on the successional trajectories of bacterial communities in a seasonally stratified lake with an anoxic hypolimnion. The bacterial community composition of epilimnion, metalimnion, and hypolimnion diverged during summer stratification and converged when the lake was mixed. In contrast, bacterial communities in the sediment remained relatively stable over the year. Phototrophic Cyanobacteria and heterotrophic Actinobacteria, Alphaproteobacteria and Planktomycetes were abundant in the aerobic epilimnion, Gammaproteobacteria (mainly Chromatiaceae) dominated in the metalimnion, and Chlorobi, Betaproteobacteria, Deltaproteobacteria, and Firmicutes were abundant in the anoxic sulfidic hypolimnion. Anoxic but nonsulfidic conditions expanded to the surface layer during fall turnover, when the epilimnion, metalimnion and upper hypolimnion mixed. During this period, phototrophic sulfur bacteria (Chromatiaceae and Chlorobi) disappeared, Polynucleobacter (Betaproteobacteria) and Methylobacter (Gammaproteobacteria) spread out from the former meta- and hypolimnion to the surface layer, and Epsilonproteobacteria dominated in the bottom water layer. Cyanobacteria and Planktomycetes regained dominance in early spring, after the oxygen concentration was restored by winter mixing. In total, these results show large spatio-temporal changes in bacterial community composition, especially during transitions from oxic to anoxic and from sulfidic to nonsulfidic conditions.

Ubiquity and quantitative significance of bacterioplankton lineages inhabiting the oxygenated hypolimnion of deep freshwater lakes

The oxygenated hypolimnion accounts for a volumetrically significant part of the global freshwater systems. Previous studies have proposed the presence of hypolimnion-specific bacterioplankton lineages that are distinct from those inhabiting the epilimnion. To date, however, no consensus exists regarding their ubiquity and abundance, which is necessary to evaluate their ecological importance. The present study investigated the bacterioplankton community in the oxygenated hypolimnia of 10 deep freshwater lakes. Despite the broad geochemical characteristics of the lakes, 16S rRNA gene sequencing demonstrated that the communities in the oxygenated hypolimnia were distinct from those in the epilimnia and identified several predominant lineages inhabiting multiple lakes. Catalyzed reporter deposition fluorescence in situ hybridization revealed that abundant hypolimnion-specific lineages, CL500-11 (Chloroflexi), CL500-3, CL500-37, CL500-15 (Planctomycetes) and Marine Group I (Thaumarchaeota), together accounted for 1.5-32.9% of all bacterioplankton in the hypolimnion of the lakes. Furthermore, an analysis of single-nucleotide variation in the partial 16S rRNA gene sequence (oligotyping) suggested the presence of different sub-populations between lakes and water layers among the lineages occurring in the entire water layer (for example, acI-B1 and acI-A7). Collectively, these results provide the first comprehensive overview of the bacterioplankton community in the oxygenated hypolimnion of deep freshwater lakes.

Spatiotemporal analysis of microbial community dynamics during seasonal stratification events in a freshwater lake (Grand Lake, OK, USA)

Many freshwater lakes undergo seasonal stratification, where the formation of phototrophic blooms in the epilimnion and subsequent sedimentation induces hypoxia/anoxia in the thermocline and hypolimnion. This autochthonously produced biomass represents a major seasonal organic input that impacts the entire ecosystem. While the limnological aspects of this process are fairly well documented, relatively little is known regarding the microbial community response to such events, especially in the deeper anoxic layers of the water column. Here, we conducted a spatiotemporal survey of the particle-associated and free-living microbial communities in a warm monomictic freshwater reservoir (Grand Lake O’ the Cherokees) in northeastern Oklahoma, USA. Pre-stratification samples (March) harbored a homogeneous community throughout the oxygenated water column dominated by typical oligotrophic aquatic lineages (acl clade within Actinobacteria, and Flavobacterium within the Bacteroidetes). The onset of phototrophic blooming in June induced the progression of this baseline community into two distinct trajectories. Within the oxic epilimnion, samples were characterized by the propagation of phototrophic (Prochlorococcus), and heterotrophic (Planctomycetes, Verrucomicrobia, and Beta-Proteobacteria) lineages. Within the oxygen-deficient thermocline and hypolimnion, the sedimentation of surface biomass induced the development of a highly diverse community, with the enrichment of Chloroflexi, “Latescibacteria”, Armatimonadetes, and Delta-Proteobacteria in the particle-associated fraction, and Gemmatimonadetes and “Omnitrophica” in the free-living fraction. Our work documents the development of multiple spatially and temporally distinct niches during lake stratification, and supports the enrichment of multiple yet-uncultured and poorly characterized lineages in the lake’s deeper oxygen-deficient layers, an ecologically relevant microbial niche that is often overlooked in lakes diversity surveys.

Competition and niche separation of pelagic bacteria in freshwater habitats

Freshwater bacterioplankton assemblages are composed of sympatric populations that can be delineated, for example, by ribosomal RNA gene relatedness and that differ in key ecophysiological properties. They may be free-living or attached, specialized for particular concentrations or subsets of substrates, or invest a variable amount of their resources in defence traits against protistan predators and viruses. Some may be motile and tactic whereas others are not, with far-reaching implications for their respective life styles and niche partitioning. The co-occurrence of competitors with overlapping growth requirements has profound consequences for the stability of community functions; it can to some extent be explained by habitat factors such as the microscale complexity and spatiotemporal variability of the lacustrine environments. On the other hand, the composition and diversity of freshwater microbial assemblages also reflects non-equilibrium states, dispersal and the stochasticity of community assembly processes. This review synoptically discusses the competition and niche separation of heterotrophic bacterial populations (defined at various levels of phylogenetic resolution) in the pelagic zone of inland surface waters from a variety of angles, focusing on habitat heterogeneity and the resulting biogeographic distribution patterns, the ecophysiological adaptations to the substrate field and the interactions of prokaryotes with predators and viruses.