Chloroflexi CL500-11 Populations That Predominate Deep-Lake Hypolimnion Bacterioplankton Rely on Nitrogen-Rich Dissolved Organic Matter Metabolism and C1 Compound Oxidation

Niche-Aware Metagenomic Screening for Enzyme Methioninase Illuminates Its Contribution to Metabolic Syntrophy

The single-step methioninase-mediated degradation of methionine (as a sulfur containing amino acid) is a reaction at the interface of carbon, nitrogen, sulfur, and methane metabolism in microbes. This enzyme also has therapeutic application due to its role in starving auxotrophic cancer cells. Applying our refined in silico screening pipeline on 33,469 publicly available genome assemblies and 1878 metagenome assembled genomes/single-cell amplified genomes from brackish waters of the Caspian Sea and the Fennoscandian Shield deep groundwater resulted in recovering 1845 methioninases. The majority of recovered methioninases belong to representatives of phyla Proteobacteria (50%), Firmicutes (29%), and Firmicutes_A (13%). Prevalence of methioninase among anaerobic microbes and in the anoxic deep groundwater together with the relevance of its products for energy conservation in anaerobic metabolism highlights such environments as desirable targets for screening novel methioninases and resolving its contribution to microbial metabolism and interactions. Among archaea, majority of detected methioninases are from representatives of Methanosarcina that are able to use methanethiol, the sulfur containing product from methionine degradation, as a precursor for methanogenesis. Branching just outside these archaeal methioninases in the phylogenetic tree, we recovered three methioninases belonging to representatives of Patescibacteria reconstructed from deep groundwater metagenomes. We hypothesize that methioninase in Patescibacteria could contribute to their syntrophic interactions where their methanogenic partners/hosts benefit from the produced 2-oxobutyrate and methanethiol. Our results underscore the significance of accounting for specific ecological niche in screening for enzyme variates with desired characteristics. Finally, complementing of our findings with experimental validation of methioninase activity confirms the potential of our in silico screening in clarifying the peculiar ecological role of methioninase in anoxic environments.

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.

Taxonomic Re-Classification and Expansion of the Phylum Chloroflexota Based on over 5000 Genomes and Metagenome-Assembled Genomes

The phylum Chloroflexota (formerly Chloroflexi) encompasses metabolically diverse bacteria that often have high prevalence in terrestrial and aquatic habitats, some even with biotechnological application. However, there is substantial disagreement in public databases which lineage should be considered a member of the phylum and at what taxonomic level. Here, we addressed these issues through extensive phylogenomic analyses. The analyses were based on a collection of >5000 Chloroflexota genomes and metagenome-assembled genomes (MAGs) from public databases, novel environmental sites, as well as newly generated MAGs from publicly available sequence reads via an improved binning approach incorporating covariance information. Based on calculated relative evolutionary divergence, we propose that Candidatus Dormibacterota should be listed as a class (i.e., Ca. Dormibacteria) within Chloroflexota together with the classes Anaerolineae, Chloroflexia, Dehalococcoidia, Ktedonobacteria, Ca. Limnocylindria, Thermomicrobia, and two other classes containing only uncultured members. All other Chloroflexota lineages previously listed at the class rank appear to be rather orders or families in the Anaerolineae and Dehalococcoidia, which contain the vast majority of genomes and exhibited the strongest phylogenetic radiation within the phylum. Furthermore, the study suggests that a common ecophysiological capability of members of the phylum is to successfully cope with low energy fluxes.

Thermophilic Dehalococcoidia with unusual traits shed light on an unexpected past

Although the phylum Chloroflexota is ubiquitous, its biology and evolution are poorly understood due to limited cultivability. Here, we isolated two motile, thermophilic bacteria from hot spring sediments belonging to the genus Tepidiforma and class Dehalococcoidia within the phylum Chloroflexota. A combination of cryo-electron tomography, exometabolomics, and cultivation experiments using stable isotopes of carbon revealed three unusual traits: flagellar motility, a peptidoglycan-containing cell envelope, and heterotrophic activity on aromatics and plant-associated compounds. Outside of this genus, flagellar motility has not been observed in Chloroflexota, and peptidoglycan-containing cell envelopes have not been described in Dehalococcoidia. Although these traits are unusual among cultivated Chloroflexota and Dehalococcoidia, ancestral character state reconstructions showed flagellar motility and peptidoglycan-containing cell envelopes were ancestral within the Dehalococcoidia, and subsequently lost prior to a major adaptive radiation of Dehalococcoidia into marine environments. However, despite the predominantly vertical evolutionary histories of flagellar motility and peptidoglycan biosynthesis, the evolution of enzymes for degradation of aromatics and plant-associated compounds was predominantly horizontal and complex. Together, the presence of these unusual traits in Dehalococcoidia and their evolutionary histories raise new questions about the timing and selective forces driving their successful niche expansion into global oceans.

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'.

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.

Diversity, distribution, and expression of opsin genes in freshwater lakes

Microbial rhodopsins are widely distributed in aquatic environments and may significantly contribute to phototrophy and energy budgets in global oceans. However, the study of freshwater rhodopsins has been largely limited. Here, we explored the diversity, ecological distribution, and expression of opsin genes that encode the apoproteins of type I rhodopsins in humic and clearwater lakes with contrasting physicochemical and optical characteristics. Using metagenomes and metagenome-assembled genomes, we recovered opsin genes from a wide range of taxa, mostly predicted to encode green light-absorbing proton pumps. Viral opsin and novel bacterial opsin clades were recovered. Opsin genes occurred more frequently in taxa from clearwater than from humic water, and opsins in some taxa have nontypical ion-pumping motifs that might be associated with physicochemical conditions of these two freshwater types. Analyses of the surface layer of 33 freshwater systems revealed an inverse correlation between opsin gene abundance and lake dissolved organic carbon (DOC). In humic water with high terrestrial DOC and light-absorbing humic substances, opsin gene abundance was low and dramatically declined within the first few meters, whereas the abundance remained relatively high along the bulk water column in clearwater lakes with low DOC, suggesting opsin gene distribution is influenced by lake optical properties and DOC. Gene expression analysis confirmed the significance of rhodopsin-based phototrophy in clearwater lakes and revealed different diel expressional patterns among major phyla. Overall, our analyses revealed freshwater opsin diversity, distribution and expression patterns, and suggested the significance of rhodopsin-based phototrophy in freshwater energy budgets, especially in clearwater lakes.

Novel nitrite reductase domain structure suggests a chimeric denitrification repertoire in the phylum Chloroflexi

Denitrification plays a central role in the global nitrogen cycle, reducing and removing nitrogen from marine and terrestrial ecosystems. The flux of nitrogen species through this pathway has a widespread impact, affecting ecological carrying capacity, agriculture, and climate. Nitrite reductase (Nir) and nitric oxide reductase (NOR) are the two central enzymes in this pathway. Here we present a previously unreported Nir domain architecture in members of phylum Chloroflexi. Phylogenetic analyses of protein domains within Nir indicate that an ancestral horizontal transfer and fusion event produced this chimeric domain architecture. We also identify an expanded genomic diversity of a rarely reported NOR subtype, eNOR. Together, these results suggest a greater diversity of denitrification enzyme arrangements exist than have been previously reported.

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.

Bacterium Lacking a Known Gene for Retinal Biosynthesis Constructs Functional Rhodopsins

Microbial rhodopsins, comprising a protein moiety (rhodopsin apoprotein) bound to the light-absorbing chromophore retinal, function as ion pumps, ion channels, or light sensors. However, recent genomic and metagenomic surveys showed that some rhodopsin-possessing prokaryotes lack the known genes for retinal biosynthesis. Since rhodopsin apoproteins cannot absorb light energy, rhodopsins produced by prokaryotic strains lacking genes for retinal biosynthesis are hypothesized to be non-functional in cells. In the present study, we investigated whether Aurantimicrobium minutum KNC^T, which is widely distributed in terrestrial environments and lacks any previously identified retinal biosynthesis genes, possesses functional rhodopsin. We initially measured ion transport activity in cultured cells. A light-induced pH change in a cell suspension of rhodopsin-possessing bacteria was detected in the absence of exogenous retinal. Furthermore, spectroscopic analyses of the cell lysate and HPLC-MS/MS analyses revealed that this strain contained an endogenous retinal. These results confirmed that A. minutum KNC^T possesses functional rhodopsin and, hence, produces retinal via an unknown biosynthetic pathway. These results suggest that rhodopsin-possessing prokaryotes lacking known retinal biosynthesis genes also have functional rhodopsins.

Genomics, Exometabolomics, and Metabolic Probing Reveal Conserved Proteolytic Metabolism of Thermoflexus hugenholtzii and Three Candidate Species From China and Japan

Thermoflexus hugenholtzii JAD2^T, the only cultured representative of the Chloroflexota order Thermoflexales, is abundant in Great Boiling Spring (GBS), NV, United States, and close relatives inhabit geothermal systems globally. However, no defined medium exists for T. hugenholtzii JAD2^T and no single carbon source is known to support its growth, leaving key knowledge gaps in its metabolism and nutritional needs. Here, we report comparative genomic analysis of the draft genome of T. hugenholtzii JAD2^T and eight closely related metagenome-assembled genomes (MAGs) from geothermal sites in China, Japan, and the United States, representing “Candidatus Thermoflexus japonica,” “Candidatus Thermoflexus tengchongensis,” and “Candidatus Thermoflexus sinensis.” Genomics was integrated with targeted exometabolomics and ¹³C metabolic probing of T. hugenholtzii. The Thermoflexus genomes each code for complete central carbon metabolic pathways and an unusually high abundance and diversity of peptidases, particularly Metallo- and Serine peptidase families, along with ABC transporters for peptides and some amino acids. The T. hugenholtzii JAD2^T exometabolome provided evidence of extracellular proteolytic activity based on the accumulation of free amino acids. However, several neutral and polar amino acids appear not to be utilized, based on their accumulation in the medium and the lack of annotated transporters. Adenine and adenosine were scavenged, and thymine and nicotinic acid were released, suggesting interdependency with other organisms in situ. Metabolic probing of T. hugenholtzii JAD2^T using ¹³C-labeled compounds provided evidence of oxidation of glucose, pyruvate, cysteine, and citrate, and functioning glycolytic, tricarboxylic acid (TCA), and oxidative pentose-phosphate pathways (PPPs). However, differential use of position-specific ¹³C-labeled compounds showed that glycolysis and the TCA cycle were uncoupled. Thus, despite the high abundance of Thermoflexus in sediments of some geothermal systems, they appear to be highly focused on chemoorganotrophy, particularly protein degradation, and may interact extensively with other microorganisms in situ.

From ecophysiology to cultivation methodology: filling the knowledge gap between uncultured and cultured microbes

Earth is dominated by a myriad of microbial communities, but the majority fails to grow under in situ laboratory conditions. The basic cause of unculturability is that bacteria dominantly occur as biofilms in natural environments. Earlier improvements in the culture techniques are mostly done by optimizing media components. However, with technological advancement particularly in the field of genome sequencing and cell imagining techniques, new tools have become available to understand the ecophysiology of microbial communities. Hence, it becomes easier to mimic environmental conditions in the culture plate. Other methods include co-culturing, emendation of growth factors, and cultivation after physical cell sorting. Most recently, techniques have been proposed for bacterial cultivation by employing genomic data to understand either microbial interactions (network-directed targeted bacterial isolation) or ecosystem engineering (reverse genomics). Hopefully, these techniques may be applied to almost all environmental samples, and help fill the gaps between the cultured and uncultured microbial communities.

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.

Bacterial Diversity in a Dynamic and Extreme Sub-Arctic Watercourse (Pasvik River, Norwegian Arctic)

Microbial communities promptly respond to the environmental perturbations, especially in the Arctic and sub-Arctic systems that are highly impacted by climate change, and fluctuations in the diversity level of microbial assemblages could give insights on their expected response. 16S rRNA gene amplicon sequencing was applied to describe the bacterial community composition in water and sediment through the sub-Arctic Pasvik River. Our results showed that river water and sediment harbored distinct communities in terms of diversity and composition at genus level. The distribution of the bacterial communities was mainly affected by both salinity and temperature in sediment samples, and by oxygen in water samples. Glacial meltwaters and runoff waters from melting ice probably influenced the composition of the bacterial community at upper and middle river sites. Interestingly, marine-derived bacteria consistently accounted for a small proportion of the total sequences and were also more prominent in the inner part of the river. Results evidenced that particular conditions occurring at sampling sites (such as algal blooms, heavy metal contamination and anaerobiosis) may select species at local scale from a shared bacterial pool, thus favoring certain bacterial taxa. Conversely, the few phylotypes specifically detected in some sites are probably due to localized external inputs introducing allochthonous microbial groups.

Freshwater microbial community diversity in a rapidly changing High Arctic watershed

Current models predict increases in High Arctic temperatures and precipitation that will have profound impacts on the Arctic hydrological cycle, including enhanced glacial melt and thawing of active layer soils. However, it remains uncertain how these changes will impact the structure of downstream resident freshwater microbial communities and ensuing microbially driven freshwater ecosystem services. Using the Lake Hazen watershed (Nunavut, Canada; 82°N, 71°W) as a sentinel system, we related microbial community composition (16S rRNA gene sequencing) to physicochemical parameters (e.g. dissolved oxygen and nutrients) over an annual hydrological cycle in three freshwater compartments within the watershed: (i) glacial rivers; (ii) active layer thaw-fed streams and waterbodies and (iii) Lake Hazen, into which (i) and (ii) drain. Microbial communities throughout these freshwater compartments were strongly interconnected, hydrologically, and often correlated with the presence of melt-sourced chemicals (e.g. dissolved inorganic carbon) as the melt season progressed. Within Lake Hazen itself, water column microbial communities were generally stable over spring and summer, despite fluctuating lake physicochemistry, indicating that these communities and the potential ecosystem services they provide therein may be resilient to environmental change. This work helps to establish a baseline understanding of how microbial communities and the ecosystem services they provide in Arctic watersheds might respond to future climate change.

Temperature and Nutrient Levels Correspond with Lineage-Specific Microdiversification in the Ubiquitous and Abundant Freshwater Genus Limnohabitans

Most freshwater bacterial communities are characterized by a few dominant taxa that are often ubiquitous across freshwater biomes worldwide. Our understanding of the genomic diversity within these taxonomic groups is limited to a subset of taxa. Here, we investigated the genomic diversity that enables Limnohabitans, a freshwater genus key in funneling carbon from primary producers to higher trophic levels, to achieve abundance and ubiquity. We reconstructed eight putative Limnohabitans metagenome-assembled genomes (MAGs) from stations located along broad environmental gradients existing in Lake Michigan, part of Earth's largest surface freshwater system. De novo strain inference analysis resolved a total of 23 strains from these MAGs, which strongly partitioned into two habitat-specific clusters with cooccurring strains from different lineages. The largest number of strains belonged to the abundant LimB lineage, for which robust in situ strain delineation had not previously been achieved. Our data show that temperature and nutrient levels may be important environmental parameters associated with microdiversification within the Limnohabitans genus. In addition, strains predominant in low- and high-phosphorus conditions had larger genomic divergence than strains abundant under different temperatures. Comparative genomics and gene expression analysis yielded evidence for the ability of LimB populations to exhibit cellular motility and chemotaxis, a phenotype not yet associated with available Limnohabitans isolates. Our findings broaden historical marker gene-based surveys of Limnohabitans microdiversification and provide in situ evidence of genome diversity and its functional implications across freshwater gradients.IMPORTANCE Limnohabitans is an important bacterial taxonomic group for cycling carbon in freshwater ecosystems worldwide. Here, we examined the genomic diversity of different Limnohabitans lineages. We focused on the LimB lineage of this genus, which is globally distributed and often abundant, and its abundance has shown to be largely invariant to environmental change. Our data show that the LimB lineage is actually comprised of multiple cooccurring populations for which the composition and genomic characteristics are associated with variations in temperature and nutrient levels. The gene expression profiles of this lineage suggest the importance of chemotaxis and motility, traits that had not yet been associated with the Limnohabitans genus, in adapting to environmental conditions.

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.

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.

Wide Distribution of Phage That Infect Freshwater SAR11 Bacteria

Fonsibacter (LD12 subclade) is among the most abundant bacterioplankton in freshwater ecosystems. These bacteria belong to the order Pelagibacterales (SAR11) and are related to Pelagibacter (marine SAR11), which dominates many marine habitats. Although a few Pelagibacter phage (Pelagiphage) have been described, no phage that infect Fonsibacter have been reported. In this study, we describe two groups of Podoviridae phage that infect Fonsibacter A complete Fonsibacter genome containing a prophage was reconstructed from metagenomic data. A circularized and complete genome related to the prophage, referred to as uv-Fonsiphage-EPL (lysogenic strategy), shows high similarity to marine Pelagiphage HTVC025P. Additionally, we reconstructed three complete genomes and one draft genome of phage related to marine Pelagiphage HTVC010P and predicted a lytic strategy. The similarity in codon usage and cooccurrence patterns of HTVC010P-related phage and Fonsibacter suggested that these phage infect Fonsibacter Similar phage were detected in Lake Mendota, Wisconsin, where Fonsibacter is also present. A search of related phage revealed the worldwide distribution of some genotypes in freshwater ecosystems, suggesting their substantial role in shaping indigenous microbial assemblages and influence on biogeochemical cycling. However, the uv-Fonsiphage-EPL and one group of HTVC010P-related phage have a more limited distribution in freshwater ecosystems. Overall, the findings provide insights into the genomic features of phage that infect Fonsibacter and expand understanding of the ecology and evolution of these important bacteria.IMPORTANCE Fonsibacter represents a significant microbial group of freshwater ecosystems. Although the genomic and metabolic features of these bacteria have been well studied, no phage infecting them has been reported. In this study, we reconstructed complete genomes of Fonsibacter and infecting phage and revealed their close relatedness to the phage infecting marine SAR11 members. Also, we illustrated that phage that infect Fonsibacter are widely distributed in freshwater habitats. In summary, the results contribute new insights into the ecology and evolution of Fonsibacter and phage.

Wide Distribution of Phage That Infect Freshwater SAR11 Bacteria

Fonsibacter represents a significant microbial group of freshwater ecosystems. Although the genomic and metabolic features of these bacteria have been well studied, no phage infecting them has been reported. In this study, we reconstructed complete genomes of Fonsibacter and infecting phage and revealed their close relatedness to the phage infecting marine SAR11 members. Also, we illustrated that phage that infect Fonsibacter are widely distributed in freshwater habitats. In summary, the results contribute new insights into the ecology and evolution of Fonsibacter and phage.

Fonsibacter (LD12 subclade) is among the most abundant bacterioplankton in freshwater ecosystems. These bacteria belong to the order Pelagibacterales (SAR11) and are related to Pelagibacter (marine SAR11), which dominates many marine habitats. Although a few Pelagibacter phage (Pelagiphage) have been described, no phage that infect Fonsibacter have been reported. In this study, we describe two groups of Podoviridae phage that infect Fonsibacter. A complete Fonsibacter genome containing a prophage was reconstructed from metagenomic data. A circularized and complete genome related to the prophage, referred to as uv-Fonsiphage-EPL (lysogenic strategy), shows high similarity to marine Pelagiphage HTVC025P. Additionally, we reconstructed three complete genomes and one draft genome of phage related to marine Pelagiphage HTVC010P and predicted a lytic strategy. The similarity in codon usage and cooccurrence patterns of HTVC010P-related phage and Fonsibacter suggested that these phage infect Fonsibacter. Similar phage were detected in Lake Mendota, Wisconsin, where Fonsibacter is also present. A search of related phage revealed the worldwide distribution of some genotypes in freshwater ecosystems, suggesting their substantial role in shaping indigenous microbial assemblages and influence on biogeochemical cycling. However, the uv-Fonsiphage-EPL and one group of HTVC010P-related phage have a more limited distribution in freshwater ecosystems. Overall, the findings provide insights into the genomic features of phage that infect Fonsibacter and expand understanding of the ecology and evolution of these important bacteria.

IMPORTANCEFonsibacter represents a significant microbial group of freshwater ecosystems. Although the genomic and metabolic features of these bacteria have been well studied, no phage infecting them has been reported. In this study, we reconstructed complete genomes of Fonsibacter and infecting phage and revealed their close relatedness to the phage infecting marine SAR11 members. Also, we illustrated that phage that infect Fonsibacter are widely distributed in freshwater habitats. In summary, the results contribute new insights into the ecology and evolution of Fonsibacter and phage.

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.

Gene Expansion and Positive Selection as Bacterial Adaptations to Oligotrophic Conditions

By combining a genome-centric metagenomic approach with a culture-based approach, we investigated the genomic adaptations of prevalent populations in an engineered oligotrophic freshwater system. We found evidence for widespread positive selection on genes involved in phosphorus and carbon scavenging pathways and for gene expansions in motility and environmental sensing to be important genomic adaptations of the abundant taxon in this system. In addition, microscopic and flow cytometric analysis of the first freshwater representative of this population (Ramlibacter aquaticus LMG 30558^T) demonstrated phenotypic plasticity, possibly due to the metabolic versatility granted by its larger genome, to be a strategy to cope with nutrient limitation. Our study clearly demonstrates the need for the use of a broad set of genomic tools combined with culture-based physiological characterization assays to investigate and validate genomic adaptations.

We examined the genomic adaptations of prevalent bacterial taxa in a highly nutrient- and ion-depleted freshwater environment located in the secondary cooling water system of a nuclear research reactor. Using genome-centric metagenomics, we found that none of the prevalent bacterial taxa were related to typical freshwater bacterial lineages. We also did not identify strong signatures of genome streamlining, which has been shown to be one of the ecoevolutionary forces shaping the genome characteristics of bacterial taxa in nutrient-depleted environments. Instead, focusing on the dominant taxon, a novel Ramlibacter sp. which we propose to name Ramlibacter aquaticus, we detected extensive positive selection on genes involved in phosphorus and carbon scavenging pathways. These genes were involved in the high-affinity phosphate uptake and storage into polyphosphate granules, metabolism of nitrogen-rich organic matter, and carbon/energy storage into polyhydroxyalkanoate. In parallel, comparative genomics revealed a high number of paralogs and an accessory genome significantly enriched in environmental sensing pathways (i.e., chemotaxis and motility), suggesting extensive gene expansions in R. aquaticus. The type strain of R. aquaticus (LMG 30558^T) displayed optimal growth kinetics and productivity at low nutrient concentrations, as well as substantial cell size plasticity. Our findings with R. aquaticus LMG 30558^T demonstrate that positive selection and gene expansions may represent successful adaptive strategies to oligotrophic environments that preserve high growth rates and cellular productivity.

IMPORTANCE By combining a genome-centric metagenomic approach with a culture-based approach, we investigated the genomic adaptations of prevalent populations in an engineered oligotrophic freshwater system. We found evidence for widespread positive selection on genes involved in phosphorus and carbon scavenging pathways and for gene expansions in motility and environmental sensing to be important genomic adaptations of the abundant taxon in this system. In addition, microscopic and flow cytometric analysis of the first freshwater representative of this population (Ramlibacter aquaticus LMG 30558^T) demonstrated phenotypic plasticity, possibly due to the metabolic versatility granted by its larger genome, to be a strategy to cope with nutrient limitation. Our study clearly demonstrates the need for the use of a broad set of genomic tools combined with culture-based physiological characterization assays to investigate and validate genomic adaptations.

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.

Effects of Phytoremediation Treatment on Bacterial Community Structure and Diversity in Different Petroleum-Contaminated Soils

Increased exploitation and use of petroleum resources is leading to increased risk of petroleum contamination of soil and groundwater. Although phytoremediation is a widely-used and cost-effective method for rehabilitating soils polluted by petroleum, bacterial community structure and diversity in soils undergoing phytoremediation is poorly understood. We investigate bacterial community response to phytoremediation in two distinct petroleum-contaminated soils (add prepared petroleum-contaminated soils) from northwest China, Weihe Terrace soil and silty loam from loess tableland. High-throughput sequencing technology was used to compare the bacterial communities in 24 different samples, yielding 18,670 operational taxonomic units (OTUs). The dominant bacterial groups, Proteobacteria (31.92%), Actinobacteria (16.67%), Acidobacteria (13.29%) and Bacteroidetes (6.58%), increased with increasing petroleum concentration from 3000 mg/kg–10,000 mg/kg, while Crenarchaeota (13.58%) and Chloroflexi (4.7%) decreased. At the order level, RB41, Actinomycetales, Cytophagales, envOPS12, Rhodospirillales, MND1 and Xanthomonadales, except Nitrososphaerales, were dominant in Weihe Terrace soil. Bacterial community structure and diversity in the two soils were significantly different at similar petroleum concentrations. In addition, the dominant genera were affected by available nitrogen, which is strongly associated with the plants used for remediation. Overall, the bacterial community structure and diversity were markedly different in the two soils, depending on the species of plants used and the petroleum concentration.

Cultivation and genomics of the first freshwater SAR11 (LD12) isolate

Evolutionary transitions between fresh and salt water happen infrequently among bacterioplankton. Within the ubiquitous and highly abundant heterotrophic Alphaproteobacteria order Pelagibacterales (SAR11), most members live in marine habitats, but the LD12 subclade has evolved as a unique freshwater lineage. LD12 cells occur as some of the most dominant freshwater bacterioplankton, yet this group has remained elusive to cultivation, hampering a more thorough understanding of its biology. Here, we report the first successful isolation of an LD12 representative, strain LSUCC0530, using high-throughput dilution-to-extinction cultivation methods, and its complete genome sequence. Growth experiments corroborate ecological data suggesting active populations of LD12 in brackish water up to salinities of ~5. LSUCC0530 has the smallest closed genome thus far reported for a SAR11 strain (1.16 Mbp). The genome affirms many previous metabolic predictions from cultivation-independent analyses, like a complete Embden–Meyerhof–Parnas glycolysis pathway, but also provides novel insights, such as the first isocitrate dehydrogenase in LD12, a likely homologous recombination of malate synthase from outside of the SAR11 clade, and analogous substitutions of ion transporters with others that occur throughout the rest of the SAR11 clade. Growth data support metagenomic recruitment results suggesting temperature-based ecotype diversification within LD12. Key gene losses for osmolyte uptake provide a succinct hypothesis for the evolutionary transition of LD12 from salt to freshwater. For strain LSUCC0530, we propose the provisional nomenclature Candidatus fonsibacter ubiquis.

ezTree: an automated pipeline for identifying phylogenetic marker genes and inferring evolutionary relationships among uncultivated prokaryotic draft genomes

BACKGROUND

Inferring phylogenetic trees for newly recovered genomes from metagenomic samples is very useful in determining the identities of uncultivated microorganisms. Even though 16S ribosomal RNA small subunit genes have been established as "gold standard" markers for inferring phylogenetic trees, they usually cannot be assembled very well in metagenomes due to shared regions among 16S genes. Using single-copy marker genes to build genome trees has become increasingly popular for uncultivated species. Predefined marker gene sets were discovered and have been applied in various genomic studies; however these gene sets might not be adequate for novel, uncultivated, draft, or incomplete genomes. The automatic identification of marker gene sets among a set of genomes with different assembly qualities has thus become a very important task for inferring reliable phylogenetic relationships for microbial populations.

RESULTS

A computational pipeline, ezTree, was developed to automatically identify single-copy marker genes for a group of genomes and build phylogenetic trees from the marker genes. Testing ezTree on a group of proteobacteria species revealed that ezTree was highly effective in pinpointing marker genes and constructing reliable trees for different groups of bacterial genomes. Applying ezTree to genomes that were recently recovered from metagenomes also showed that ezTree can help elucidate taxonomic relationships among newly recovered genomes and existing ones.

CONCLUSIONS

The development of ezTree can help scientists build reliable phylogenetic trees for uncultivated species retrieved from environmental samples. The uncovered single-copy marker genes may also provide crucial hints for understanding shared features of a group of microbes. The ezTree pipeline is freely available at https://github.com/yuwwu/ezTree under a GNU GPLv3 license.

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.

Lake Bacterial Assemblage Composition Is Sensitive to Biological Disturbance Caused by an Invasive Filter Feeder

Freshwater bacteria play fundamental roles in global elemental cycling and are an intrinsic part of local food webs. Human activities are altering freshwater environments, and much has been learned regarding the sensitivity of bacterial assemblages to a variety of these disturbances. Yet, relatively few studies have focused on how species invasion, which is one of the most important aspects of anthropogenic global change, affects freshwater bacterial assemblages. This study focuses on the impact of invasive dreissenid mussels (IDMs), a globally distributed group of invasive species with large impacts on freshwater phyto- and zooplankton assemblages. We show that IDMs have direct effects on lake bacterioplankton abundance, taxonomic composition, and inferred bacterial functional group representation.

One approach to improve forecasts of how global change will affect ecosystem processes is to better understand how anthropogenic disturbances alter bacterial assemblages that drive biogeochemical cycles. Species invasions are important contributors to global change, but their impacts on bacterial community ecology are rarely investigated. Here, we studied direct impacts of invasive dreissenid mussels (IDMs), one of many invasive filter feeders, on freshwater lake bacterioplankton. We demonstrated that direct effects of IDMs reduced bacterial abundance and altered assemblage composition by preferentially removing larger and particle-associated bacteria. While this increased the relative abundances of many free-living bacterial taxa, some were susceptible to filter feeding, in line with efficient removal of phytoplankton cells of <2 μm. This selective removal of particle-associated and larger bacteria by IDMs altered inferred bacterial functional group representation, defined by carbon and energy source utilization. Specifically, we inferred an increased relative abundance of chemoorganoheterotrophs predicted to be capable of rhodopsin-dependent energy generation. In contrast to the few previous studies that have focused on the longer-term combined direct and indirect effects of IDMs on bacterioplankton, our study showed that IDMs act directly as a biological disturbance to which freshwater bacterial assemblages are sensitive. The negative impacts on particle-associated bacteria, which have been shown to be more active than free-living bacteria, and the inferred shifts in functional group representation raise the possibility that IDMs may directly alter bacterially mediated ecosystem functions.

IMPORTANCE Freshwater bacteria play fundamental roles in global elemental cycling and are an intrinsic part of local food webs. Human activities are altering freshwater environments, and much has been learned regarding the sensitivity of bacterial assemblages to a variety of these disturbances. Yet, relatively few studies have focused on how species invasion, which is one of the most important aspects of anthropogenic global change, affects freshwater bacterial assemblages. This study focuses on the impact of invasive dreissenid mussels (IDMs), a globally distributed group of invasive species with large impacts on freshwater phyto- and zooplankton assemblages. We show that IDMs have direct effects on lake bacterioplankton abundance, taxonomic composition, and inferred bacterial functional group representation.

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.

Freshwater bacteria release methane as a byproduct of phosphorus acquisition

Freshwater lakes emit large amounts of methane, some of which is produced in oxic surface waters. Two potential pathways for aerobic methane production exist: methanogenesis in oxygenated water, which has been observed in some lakes, or demethylation of small organic molecules. Although methane is produced via demethylation in oxic marine environments, this mechanism of methane release has not yet been demonstrated in freshwater systems. Genes related to the C-P lyase pathway, which cleaves C-P bonds in phosphonate compounds, were found in a metagenomic survey of the surface water of Lake Matano, which is chronically P-starved and methane-rich. We demonstrate that four bacterial isolates from Lake Matano obtain P from methylphosphonate and release methane, and that this activity is repressed by phosphate. We further demonstrate that expression of phnJ, which encodes the enzyme that releases methane, is higher in the presence of methylphosphonate and lower when both methylphosphonate and phosphate are added. This gene is also found in most of the metagenomic data sets from freshwater environments. These experiments link methylphosphonate degradation and methane production with gene expression and phosphate availability in freshwater organisms, and suggest that some of the excess methane in the Lake Matano surface water, and in other methane-rich lakes, may be produced by P-starved bacteria.

IMPORTANCE

Methane is an important greenhouse gas, and contributes substantially to global warming. Although freshwater environments are known to release methane into the atmosphere, estimates of the amount of methane emitted by freshwater lakes vary from 8 to 73 Tg per year. Methane emissions are difficult to predict in part because the source of the methane can vary: it is the end product of the energy-conserving pathway in methanogenic archaea, which predominantly live in anoxic sediments or waters, but have also been identified in some oxic freshwater environments. More recently, methane release from small organic molecules has been observed in oxic marine environments. Here we show that demethylation of methylphosphonate may also contribute to methane release from lakes, and that phosphate can repress this activity. Since lakes are typically phosphorus-limited, some methane release in these environments may be a byproduct of phosphorus metabolism, rather than carbon or energy metabolism. Methane emissions from lakes are currently predicted using primary production, eutrophication status, extent of anoxia, and the shape and size of the lake; to improve prediction of methane emissions, phosphorus availability and sources may also need to be included in these models.

Seasonal Succession Leads to Habitat-Dependent Differentiation in Ribosomal RNA:DNA Ratios among Freshwater Lake Bacteria

Relative abundance profiles of bacterial populations measured by sequencing DNA or RNA of marker genes can widely differ. These differences, made apparent when calculating ribosomal RNA:DNA ratios, have been interpreted as variable activities of bacterial populations. However, inconsistent correlations between ribosomal RNA:DNA ratios and metabolic activity or growth rates have led to a more conservative interpretation of this metric as the cellular protein synthesis potential (PSP). Little is known, particularly in freshwater systems, about how PSP varies for specific taxa across temporal and spatial environmental gradients and how conserved PSP is across bacterial phylogeny. Here, we generated 16S rRNA gene sequencing data using simultaneously extracted DNA and RNA from fractionated (free-living and particulate) water samples taken seasonally along a eutrophic freshwater estuary to oligotrophic pelagic transect in Lake Michigan. In contrast to previous reports, we observed frequent clustering of DNA and RNA data from the same sample. Analysis of the overlap in taxa detected at the RNA and DNA level indicated that microbial dormancy may be more common in the estuary, the particulate fraction, and during the stratified period. Across spatiotemporal gradients, PSP was often conserved at the phylum and class levels. PSPs for specific taxa were more similar across habitats in spring than in summer and fall. This was most notable for PSPs of the same taxa when located in the free-living or particulate fractions, but also when contrasting surface to deep, and estuary to Lake Michigan communities. Our results show that community composition assessed by RNA and DNA measurements are more similar than previously assumed in freshwater systems. However, the similarity between RNA and DNA measurements and taxa-specific PSPs that drive community-level similarities are conditional on spatiotemporal factors.