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

Integrating depth-dependent protist dynamics and microbial interactions in spring succession of a freshwater reservoir

BACKGROUND

Protists are essential contributors to eukaryotic diversity and exert profound influence on carbon fluxes and energy transfer in freshwaters. Despite their significance, there is a notable gap in research on protistan dynamics, particularly in the deeper strata of temperate lakes. This study aimed to address this gap by integrating protists into the well-described spring dynamics of Římov reservoir, Czech Republic. Over a 2-month period covering transition from mixing to established stratification, we collected water samples from three reservoir depths (0.5, 10 and 30 m) with a frequency of up to three times per week. Microbial eukaryotic and prokaryotic communities were analysed using SSU rRNA gene amplicon sequencing and dominant protistan groups were enumerated by Catalysed Reporter Deposition-Fluorescence in situ Hybridization (CARD-FISH). Additionally, we collected samples for water chemistry, phyto- and zooplankton composition analyses.

RESULTS

Following the rapid changes in environmental and biotic parameters during spring, protistan and bacterial communities displayed swift transitions from a homogeneous community to distinct strata-specific communities. A prevalence of auto- and mixotrophic protists dominated by cryptophytes was associated with spring algal bloom-specialized bacteria in the epilimnion. In contrast, the meta- and hypolimnion showcased a development of a protist community dominated by putative parasitic Perkinsozoa, detritus or particle-associated ciliates, cercozoans, telonemids and excavate protists (Kinetoplastida), co-occurring with bacteria associated with lake snow.

CONCLUSIONS

Our high-resolution sampling matching the typical doubling time of microbes along with the combined microscopic and molecular approach and inclusion of all main components of the microbial food web allowed us to unveil depth-specific populations' successions and interactions in a deep lentic ecosystem.

Particle-attached bacteria act as gatekeepers in the decomposition of complex phytoplankton polysaccharides

BACKGROUND

Marine microalgae (phytoplankton) mediate almost half of the worldwide photosynthetic carbon dioxide fixation and therefore play a pivotal role in global carbon cycling, most prominently during massive phytoplankton blooms. Phytoplankton biomass consists of considerable proportions of polysaccharides, substantial parts of which are rapidly remineralized by heterotrophic bacteria. We analyzed the diversity, activity, and functional potential of such polysaccharide-degrading bacteria in different size fractions during a diverse spring phytoplankton bloom at Helgoland Roads (southern North Sea) at high temporal resolution using microscopic, physicochemical, biodiversity, metagenome, and metaproteome analyses.

RESULTS

Prominent active 0.2-3 µm free-living clades comprised Aurantivirga, "Formosa", Cd. Prosiliicoccus, NS4, NS5, Amylibacter, Planktomarina, SAR11 Ia, SAR92, and SAR86, whereas BD1-7, Stappiaceae, Nitrincolaceae, Methylophagaceae, Sulfitobacter, NS9, Polaribacter, Lentimonas, CL500-3, Algibacter, and Glaciecola dominated 3-10 µm and > 10 µm particles. Particle-attached bacteria were more diverse and exhibited more dynamic adaptive shifts over time in terms of taxonomic composition and repertoires of encoded polysaccharide-targeting enzymes. In total, 305 species-level metagenome-assembled genomes were obtained, including 152 particle-attached bacteria, 100 of which were novel for the sampling site with 76 representing new species. Compared to free-living bacteria, they featured on average larger metagenome-assembled genomes with higher proportions of polysaccharide utilization loci. The latter were predicted to target a broader spectrum of polysaccharide substrates, ranging from readily soluble, simple structured storage polysaccharides (e.g., laminarin, α-glucans) to less soluble, complex structural, or secreted polysaccharides (e.g., xylans, cellulose, pectins). In particular, the potential to target poorly soluble or complex polysaccharides was more widespread among abundant and active particle-attached bacteria.

CONCLUSIONS

Particle-attached bacteria represented only 1% of all bloom-associated bacteria, yet our data suggest that many abundant active clades played a pivotal gatekeeping role in the solubilization and subsequent degradation of numerous important classes of algal glycans. The high diversity of polysaccharide niches among the most active particle-attached clades therefore is a determining factor for the proportion of algal polysaccharides that can be rapidly remineralized during generally short-lived phytoplankton bloom events. Video Abstract.

Small-scale spatial beta diversity of bacteria in the mixed upper layer of a lake

Bacterial community composition among individual, experimentally generated 'lake snow' particles may be highly variable. Since such aggregates are seasonally abundant in the mixed upper layer of lakes, we hypothesized that particle-attached (PA) bacteria disproportionally contribute to the small-scale spatial beta diversity of pelagic communities. Community composition was analysed in sets of small (10 mL) samples collected from a pre-alpine lake in May, July and October 2018. Bacteria were classified as free-living (FL) or PA depending on their presence in large, 5-μm pre-filtered reference samples. FL exhibited clear seasonal differences in community composition and assembly. They were spatially uniform in May and July, and only a few FL taxa exhibited significant spatial variability. Spatial heterogeneity of FL in October was caused by high alpha and beta diversity of rare taxa, many with a presumably 'tychoplanktic' (alternating attached and free-living) lifestyle. The spatial beta diversity of PA was always high, and only about 10% of their seasonal richness was present in any single sample. Thus, most compositional variability of pelagic bacteria at spatial scales of cm to m either directly or indirectly originated from PA. On a functional level, this genotypic heterogeneity might affect the spatial distribution of rare metabolic traits.

In the right place, at the right time: the integration of bacteria into the Plankton Ecology Group model

BACKGROUND

Planktonic microbial communities have critical impacts on the pelagic food web and water quality status in freshwater ecosystems, yet no general model of bacterial community assembly linked to higher trophic levels and hydrodynamics has been assessed. In this study, we utilized a 2-year survey of planktonic communities from bacteria to zooplankton in three freshwater reservoirs to investigate their spatiotemporal dynamics.

RESULTS

We observed site-specific occurrence and microdiversification of bacteria in lacustrine and riverine environments, as well as in deep hypolimnia. Moreover, we determined recurrent bacterial seasonal patterns driven by both biotic and abiotic conditions, which could be integrated into the well-known Plankton Ecology Group (PEG) model describing primarily the seasonalities of larger plankton groups. Importantly, bacteria with different ecological potentials showed finely coordinated successions affiliated with four seasonal phases, including the spring bloom dominated by fast-growing opportunists, the clear-water phase associated with oligotrophic ultramicrobacteria, the summer phase characterized by phytoplankton bloom-associated bacteria, and the fall/winter phase driven by decay-specialists.

CONCLUSIONS

Our findings elucidate the major principles driving the spatiotemporal microbial community distribution in freshwater ecosystems. We suggest an extension to the original PEG model by integrating new findings on recurrent bacterial seasonal trends. Video Abstract.

Adaptive genetic traits in pelagic freshwater microbes

Pelagic microbes have adopted distinct strategies to inhabit the pelagial of lakes and oceans and can be broadly categorized in two groups: free-living, specialized oligotrophs and patch-associated generalists or copiotrophs. In this review, we aim to identify genomic traits that enable pelagic freshwater microbes to thrive in their habitat. To do so, we discuss the main genetic differences of pelagic marine and freshwater microbes that are both dominated by specialized oligotrophs and the difference to freshwater sediment microbes, where copiotrophs are more prevalent. We phylogenomically analysed a collection of >7700 metagenome-assembled genomes, classified habitat preferences on different taxonomic levels, and compared the metabolic traits of pelagic freshwater, marine, and freshwater sediment microbes. Metabolic differences are mainly associated with transport functions, environmental information processing, components of the electron transport chain, osmoregulation and the isoelectric point of proteins. Several lineages with known habitat transitions (Nitrososphaeria, SAR11, Methylophilaceae, Synechococcales, Flavobacteriaceae, Planctomycetota) and the underlying mechanisms in this process are discussed in this review. Additionally, the distribution, ecology and genomic make-up of the most abundant freshwater prokaryotes are described in details in separate chapters for Actinobacteriota, Bacteroidota, Burkholderiales, Verrucomicrobiota, Chloroflexota, and 'Ca. Patescibacteria'.

A Fast and Easy Method to Co-extract DNA and RNA from an Environmental Microbial Sample

We herein propose a fast and easy DNA and RNA co-extraction method for environmental microbial samples. It combines bead beating and phenol-chloroform phase separation followed by the separation and purification of DNA and RNA using the Qiagen AllPrep DNA/RNA mini kit. With a handling time of ~3 h, our method simultaneously extracted high-quality DNA (peak size >10-15‍ ‍kb) and RNA (RNA integrity number >6) from lake bacterioplankton filtered samples. The method is also applicable to low-biomass samples (expected DNA or RNA yield <50‍ ‍ng) and eukaryotic microbial samples, providing an easy option for more versatile eco-genomic applications.

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.

Genome Streamlining, Proteorhodopsin, and Organic Nitrogen Metabolism in Freshwater Nitrifiers

Microbial nitrification is a critical process governing nitrogen availability in aquatic systems. Freshwater nitrifiers have received little attention, leaving many unanswered questions about their taxonomic distribution, functional potential, and ecological interactions. Here, we reconstructed genomes to infer the metabolism and ecology of free-living picoplanktonic nitrifiers across the Laurentian Great Lakes, a connected series of five of Earth’s largest lakes. Surprisingly, ammonia-oxidizing bacteria (AOB) related to Nitrosospira dominated over ammonia-oxidizing archaea (AOA) at nearly all stations, with distinct ecotypes prevailing in the transparent, oligotrophic upper lakes compared to Lakes Erie and Ontario. Unexpectedly, one ecotype of Nitrosospira encodes proteorhodopsin, which could enhance survival under conditions where ammonia oxidation is inhibited or substrate limited. Nitrite-oxidizing bacteria (NOB) “Candidatus Nitrotoga” and Nitrospira fluctuated in dominance, with the latter prevailing in deeper, less-productive basins. Genome reconstructions reveal highly reduced genomes and features consistent with genome streamlining, along with diverse adaptations to sunlight and oxidative stress and widespread capacity for organic nitrogen use. Our findings expand the known functional diversity of nitrifiers and establish their ecological genomics in large lake ecosystems. By elucidating links between microbial biodiversity and biogeochemical cycling, our work also informs ecosystem models of the Laurentian Great Lakes, a critical freshwater resource experiencing rapid environmental change.

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 (&gt;25 m) and mostly oxygenated lakes showed highly similar community assemblages that were distinct from those of two shallower lakes (&lt;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.

Remarkable characteristics and distinct community of biofilms on the photoaged polyethylene films in riverine microcosms

Recalcitrant plastics in the environment are gradually fragmented into weathered debris distinguished from their original state by the integrative action of influencing factors, such as UV light, heating and physical abrasion. As new artificial carbon-source substrates in aquatic ecosystems, plastic products can be colonized by biofilms and even utilized by microorganisms. To investigate the influences of weathering of plastics on the colonized biofilms, freshwater samples from the Yangtze River (Nanjing, China) were collected for biofilm incubation. Based on the characterization of plastics and biofilms, the effects of plastic surface properties on biofilm characteristics were revealed by the analysis of partial least squares regression (PLSR). Roughness was the principal influencing factor, while rigidity had the opposite effect to it. 16S rRNA gene high-throughput sequencing results indicated the high relative abundance of Cyanobacteria and rising proportion of harmful components (e.g., Flavobacterium) on photoaged polyethylene plastics. The microbial functional profiles (KEGG) predicted by Tax4Fun showed that the functions (e.g., membrane transport, energy metabolism, etc.) of biofilm on photoaged plastics were dissimilar with those on original ones. These findings suggested that the distinct microbial community and the adverse functional changes in biofilms on photoaged plastics potentially enhanced their environmental risks. On the other hand, 28-day cultured biofilms on original low-density polyethylene (LDPE) films were dominated by Exiguobacterium. The previously ignored potentials of this microorganism in rapidly accommodating to a hydrophobic substrate and its plastic degrading ability were both worthy of attention. Therefore, it is necessary to consider the weathering process of plastics in exploring the "plastisphere", and to give further insights into the double-edged nature of the "plastisphere".

Drivers and implications of distance decay differ for ectomycorrhizal and foliar endophytic fungi across an anciently fragmented landscape

Fungal communities associated with plants often decrease in similarity as the distance between sampling sites increases (i.e., they demonstrate distance decay). In the southwestern USA, forests occur in highlands separated from one another by warmer, drier biomes with plant and fungal communities that differ from those at higher elevations. These disjunct forests are broadly similar in climate to one another, offering an opportunity to examine drivers of distance decay in plant-associated fungi across multiple ecologically similar yet geographically disparate landscapes. We examined ectomycorrhizal and foliar endophytic fungi associated with a dominant forest tree (Pinus ponderosa) in forests across ca. 550 km of geographic distance from northwestern to southeastern Arizona (USA). Both guilds of fungi showed distance decay, but drivers differed for each: ectomycorrhizal fungi are constrained primarily by dispersal limitation, whereas foliar endophytes are constrained by specific environmental conditions. Most ectomycorrhizal fungi were found in only a single forested area, as were many endophytic fungi. Such regional-scale perspectives are needed for baseline estimates of fungal diversity associated with forest trees at a landscape scale, with attention to the sensitivity of different guilds of fungal symbionts to decreasing areas of suitable habitat, increasing disturbance, and related impacts of climate change.

Seasonal Variation in Viral Infection Rates and Cell Sizes of Infected Prokaryotes in a Large and Deep Freshwater Lake (Lake Biwa, Japan)

As viruses regulate prokaryotic abundance and the carbon cycle by infecting and lysing their prokaryotic hosts, the volume of infected prokaryotes is an important parameter for understanding the impact of viruses on aquatic environments. However, literature regarding the seasonal and spatial variations in the cell volume of infected prokaryotes is limited, despite the volume of the prokaryotic community varying dynamically with season and water column depth. Here, we conducted a field survey for two annual cycles in a large and deep freshwater lake (Lake Biwa, Japan), where large prokaryotes inhabit the deeper layer during the stratified period. We used transmission electron microscopy to reveal the seasonal and spatial variation in the frequency of viral infection and cell volume of infected prokaryotes. We found that the viral infection rate in the surface layer increased when estimated contact rates increased during the middle of the stratified period, whereas the infection rate in the deeper layer increased despite low estimated contact rates during the end of the stratified period. In addition, in the deeper layer, the fraction of large prokaryotes in the total and infected prokaryotic communities increased progressively while the number of intracellular viral particles increased. We suggest different ways in which the viral abundance is maintained in the two water layers. In the surface layer, it is speculated that viral abundance is supported by the high viral infection rate because of the high activity of prokaryotes, whereas in the deeper layer, it might be supported by the larger number of intracellular viral particles released from large prokaryotes. Moreover, large prokaryotes could contribute as important sources of organic substrates via viral lysis in the deeper layer, where labile dissolved organic matter is depleted.

Metagenomics combined with DNA-based stable isotope probing provide comprehensive insights of active triclosan-degrading bacteria in wastewater treatment

The biotransformation of triclosan (TCS) during wastewater treatment occurred frequently, while little researches are known the identity of microorganisms involved in the biodegradation process. In this work, DNA-based stable isotope probing (DNA-SIP) was occupied to investigate the TCS assimilation microbes originated from a full-scale cyclic activated sludge system in Beijing. Results of TCS removal pathway showed that the TCS removal in nitrification process was mainly contributed by the metabolism of heterotrophic bacteria, accounting for about 18.54%. DNA-SIP assay indicated that Sphingobium dominated the degradation of TCS. Oligotyping analysis further indicated that oligotype GCTAAT and ATGTTA of Sphingobium played important roles in degrading TCS. Furthermore, the Kyoto Encyclopedia of Genes and Genomes functional abundance statistics based on PICRUSt2 showed that glutathione transferase was the most prevalent enzyme involved in TCS metabolism, and TCS might be removed through microbial carbon metabolism. Metagenomics made clear that Sphingobium might play irrelevant role on the propagation of antibiotics resistance genes (ARGs), even though, it could degrade TCS. Thauera and Dechloromonas were identified as the key hosts of most ARGs. This study revealed the potential metabolic pathway and microbial ecology of TCS biodegradation in nitrification process of wastewater treatment system.

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.

Metagenomics Approaches for the Detection and Surveillance of Emerging and Recurrent Plant Pathogens

Globalization has a dramatic effect on the trade and movement of seeds, fruits and vegetables, with a corresponding increase in economic losses caused by the introduction of transboundary plant pathogens. Current diagnostic techniques provide a useful and precise tool to enact surveillance protocols regarding specific organisms, but this approach is strictly targeted, while metabarcoding and shotgun metagenomics could be used to simultaneously detect all known pathogens and potentially new ones. This review aims to present the current status of high-throughput sequencing (HTS) diagnostics of fungal and bacterial plant pathogens, discuss the challenges that need to be addressed, and provide direction for the development of methods for the detection of a restricted number of related taxa (specific surveillance) or all of the microorganisms present in a sample (general surveillance). HTS techniques, particularly metabarcoding, could be useful for the surveillance of soilborne, seedborne and airborne pathogens, as well as for identifying new pathogens and determining the origin of outbreaks. Metabarcoding and shotgun metagenomics still suffer from low precision, but this issue can be limited by carefully choosing primers and bioinformatic algorithms. Advances in bioinformatics will greatly accelerate the use of metagenomics to address critical aspects related to the detection and surveillance of plant pathogens in plant material and foodstuffs.

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.

Bacterial and Eukaryotic Small-Subunit Amplicon Data Do Not Provide a Quantitative Picture of Microbial Communities, but They Are Reliable in the Context of Ecological Interpretations

High-throughput sequencing (HTS) of amplified fragments of rRNA genes provides unprecedented insight into the diversity of prokaryotic and eukaryotic microorganisms. Unfortunately, HTS data are prone to quantitative biases, which may lead to an erroneous picture of microbial community composition and thwart efforts to advance its understanding. These concerns motivated us to investigate how accurately HTS data characterize the variability of microbial communities, the relative abundances of specific phylotypes, and their relationships with environmental factors in comparison to an established microscopy-based method. We compared results obtained by HTS and catalyzed reporter deposition-fluorescence in situ hybridization (CARD-FISH) from three independent aquatic time series for both prokaryotic and eukaryotic microorganisms (almost 900 data points, the largest obtained with both methods so far). HTS and CARD-FISH data disagree with regard to relative abundances of bacterial and eukaryotic phylotypes but identify similar environmental drivers shaping bacterial and eukaryotic communities.

High-throughput sequencing (HTS) of gene amplicons is a preferred method of assessing microbial community composition, because it rapidly provides information from a large number of samples at high taxonomic resolution and low costs. However, mock community studies show that HTS data poorly reflect the actual relative abundances of individual phylotypes, casting doubt on the reliability of subsequent statistical analysis and data interpretation. We investigated how accurately HTS data reflect the variability of bacterial and eukaryotic community composition and their relationship with environmental factors in natural samples. For this, we compared results of HTS from three independent aquatic time series (n = 883) with those from an established, quantitative microscopic method (catalyzed reporter deposition-fluorescence in situ hybridization [CARD-FISH]). Relative abundances obtained by CARD-FISH and HTS disagreed for most bacterial and eukaryotic phylotypes. Nevertheless, the two methods identified the same environmental drivers to shape bacterial and eukaryotic communities. Our results show that amplicon data do provide reliable information for their ecological interpretations. Yet, when studying specific phylogenetic groups, it is advisable to combine HTS with quantification using microscopy and/or the addition of internal standards.

IMPORTANCE High-throughput sequencing (HTS) of amplified fragments of rRNA genes provides unprecedented insight into the diversity of prokaryotic and eukaryotic microorganisms. Unfortunately, HTS data are prone to quantitative biases, which may lead to an erroneous picture of microbial community composition and thwart efforts to advance its understanding. These concerns motivated us to investigate how accurately HTS data characterize the variability of microbial communities, the relative abundances of specific phylotypes, and their relationships with environmental factors in comparison to an established microscopy-based method. We compared results obtained by HTS and catalyzed reporter deposition-fluorescence in situ hybridization (CARD-FISH) from three independent aquatic time series for both prokaryotic and eukaryotic microorganisms (almost 900 data points, the largest obtained with both methods so far). HTS and CARD-FISH data disagree with regard to relative abundances of bacterial and eukaryotic phylotypes but identify similar environmental drivers shaping bacterial and eukaryotic communities.

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.

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.

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.

Influence of salinity cycles in bioreactor performance and microbial community structure of membrane-based tidal-like variable salinity wastewater treatment systems

A membrane bioreactor and two hybrid moving bed bioreactor-membrane bioreactors were operated for the treatment of variable salinity wastewater, changing in cycles of 6-h wastewater base salinity and 6-h maximum salinity (4.5 and 8.5 mS cm^-1 electric conductivity, which relate to 2.4 and 4.8 g L^-1 NaCl, respectively), under different hydraulic retention times (6, 9.5, and 12 h) and total solids concentrations (2500 and 3500 mg L^-1). The evaluation of the performance of the systems showed that COD removal performance was unaffected by salinity conditions, while BOD₅ and TN removals were significantly higher in the low-salinity scenario. The microbial community structure showed differences with respect to salinity conditions for Eukarya, suggesting their higher sensitivity for salinity with respect to Prokarya, which were similar at both salinity scenarios. Nevertheless, the intra-OTU distribution of consistently represented OTUs of Eukarya and Prokarya was affected by the different salinity maximums. Multivariate redundancy analyses showed that several genera such as Amphiplicatus (0.01-5.90%), Parvibaculum (0.27-1.19%), Thiothrix (0.30-1.19%), Rhodanobacter (2.81-5.85%), Blastocatella (0.21-2.01%), and Nitrobacter (0.80-0.99%) were positively correlated with BOD₅ and TN removal, and the ecological roles of these were proposed. All these genera were substantially more represented under low-salinity conditions (10-500% higher relative abundance), demonstrating that they might be of importance for the treatment of variable salinity wastewater. Evaluation of Eukarya OTUs showed that many of them lack a consistent taxonomic classification, which highlights the lack of knowledge of the diversity and ecological role of Eukaryotes in saline wastewater treatment processes. The results obtained will be of interest for future design and operation of salinity wastewater treatment systems particularly because little is known on the effect of variable salinity conditions in wastewater treatment.

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.

Diversity and Cyclical Seasonal Transitions in the Bacterial Community in a Large and Deep Perialpine Lake

High-throughput sequencing (HTS) was used to analyze the seasonal variations in the bacterioplankton community composition (BCC) in the euphotic layer of a large and deep lake south of the Alps (Lake Garda). The BCC was analyzed throughout two annual cycles by monthly samplings using the amplification and sequencing of the V3-V4 hypervariable region of the 16S rRNA gene by the MiSeq Illumina platform. The dominant and most diverse bacterioplankton phyla were among the more frequently reported in freshwater ecosystems, including the Proteobacteria, Cyanobacteria, Bacteroidetes, Verrucomicrobia, Actinobacteria, and Planctomycetes. As a distinctive feature, the development of the BCC showed a cyclical temporal pattern in the two analyzed years and throughout the euphotic layer. The recurring temporal development was controlled by the strong seasonality in water temperature and thermal stratification, and by cyclical temporal changes in nutrients and, possibly, by the remarkable annual cyclical development of cyanobacteria and eukaryotic phytoplankton hosting bacterioplankton that characterizes Lake Garda. Further downstream analyses of operational taxonomic units associated to cyanobacteria allowed confirming the presence of the most abundant taxa previously identified by microscopy and/or phylogenetic analyses, as well as the presence of other small Synechococcales/Chroococcales and rare Nostocales never identified so far in the deep lakes south of the Alps. The implications of the high diversity and strong seasonality are relevant, opening perspectives for the definition of common and discriminating patterns characterizing the temporal and spatial distribution in the BCC, and for the application of the new sequencing technologies in the monitoring of water quality in large and deep lakes.

Watershed-Induced Limnological and Microbial Status in Two Oligotrophic Andean Lakes Exposed to the Same Climatic Scenario

Laguna Negra and Lo Encañado are two oligotrophic Andean lakes forming part of the system fed by meltwater from distinct glacial tongues of the Echaurren glacier in central Chile, which is in a recession period. The recent increase in temperature and decline in precipitation have led to an increase of glacial meltwater and sediments entering these lakes. Although the lacustrine systems are also hydrogeologically connected, the limnology of the lakes is strongly controlled by the surface processes related to the respective sub-watersheds and hydrology. Watershed characteristics (area and length, slope, lithology, resistance to erosion, among others) affect the chemical and physical characteristics of both lakes (e.g., nutrient concentration and turbidity). We studied physical and chemical variables and performed 16S rRNA amplicon sequencing to determine the specific microbial signature of the lakes. The transparency, temperature, turbidity and concentrations of chlorophyll-a, dissolved organic matter, nutrients and the total number of cells, revealed the different status of both lakes at the time of sampling. The predominant bacterial groups in both lakes were Proteobacteria, Verrucomicrobia, and Bacteroidetes. Interestingly, the contribution of phototrophs was significantly higher in LN compared to LE (13 and 4% respectively) and the major fraction corresponded to Anoxygenic Phototrophs (AP) represented by Chloroflexi, Alpha, and Betaproteobacteria. Multivariate analyses showed that the nutrient levels and the light availability of both lakes, which finally depend on the hydrological characteristics of the respective watersheds, explain the differential community composition/function. The abundance of a diverse photoheterotrophic bacterioplankton community suggests that the ability to utilize solar energy along with organic and inorganic substrates is a key function in these oligotrophic mountain lakes.

Microbial ecology of full-scale wastewater treatment systems in the Polar Arctic Circle: Archaea, Bacteria and Fungi

Seven full-scale biological wastewater treatment systems located in the Polar Arctic Circle region in Finland were investigated to determine their Archaea, Bacteria and Fungi community structure, and their relationship with the operational conditions of the bioreactors by the means of quantitative PCR, massive parallel sequencing and multivariate redundancy analysis. The results showed dominance of Archaea and Bacteria members in the bioreactors. The activated sludge systems showed strong selection of Bacteria but not for Archaea and Fungi, as suggested by diversity analyses. Core OTUs in influent and bioreactors were classified as Methanobrevibacter, Methanosarcina, Terrestrial Group Thaumarchaeota and unclassified Euryarchaeota member for Archaea; Trichococcus, Leptotrichiaceae and Comamonadaceae family, and Methylorosula for Bacteria and Trichosporonaceae family for Fungi. All influents shared core OTUs in all domains, but in bioreactors this did not occur for Bacteria. Oligotype structure of core OTUs showed several ubiquitous Fungi oligotypes as dominant in sewage and bioreactors. Multivariate redundancy analyses showed that the majority of core OTUs were related to organic matter and nutrients removal. Also, there was evidence of competition among Archaea and Fungi core OTUs, while all Bacteria OTUs were positively correlated among them. The results obtained highlighted interesting features of extremely cold temperature bioreactors.

Microbial community diversity patterns are related to physical and chemical differences among temperate lakes near Beaver Island, MI

Lakes are dynamic and complex ecosystems that can be influenced by physical, chemical, and biological processes. Additionally, individual lakes are often chemically and physically distinct, even within the same geographic region. Here we show that differences in physicochemical conditions among freshwater lakes located on (and around) the same island, as well as within the water column of each lake, are significantly related to aquatic microbial community diversity. Water samples were collected over time from the surface and bottom-water within four freshwater lakes located around Beaver Island, MI within the Laurentian Great Lakes region. Three of the sampled lakes experienced seasonal lake mixing events, impacting either O₂, pH, temperature, or a combination of the three. Microbial community alpha and beta diversity were assessed and individual microbial taxa were identified via high-throughput sequencing of the 16S rRNA gene. Results demonstrated that physical and chemical variability (temperature, dissolved oxygen, and pH) were significantly related to divergence in the beta diversity of surface and bottom-water microbial communities. Despite its correlation to microbial community structure in unconstrained analyses, constrained analyses demonstrated that dissolved organic carbon (DOC) concentration was not strongly related to microbial community structure among or within lakes. Additionally, several taxa were correlated (either positively or negatively) to environmental variables, which could be related to aerobic and anaerobic metabolisms. This study highlights the measurable relationships between environmental conditions and microbial communities within freshwater temperate lakes around the same island.