Contrasting trends in distribution of four major planktonic betaproteobacterial groups along a pH gradient of epilimnia of 72 freshwater habitats

Fine-scale spatiotemporal variations in bacterial community diversity in agricultural pond water

Microbial communities in surface waters are affected by environmental conditions and can influence changes in water quality. To explore the hypothesis that the microbiome in agricultural waters associates with spatiotemporal variations in overall water quality and, in turn, has implications for resource monitoring and management, we characterized the relationships between the microbiota and physicochemical properties in a model irrigation pond as a factor of sampling time (i.e., 9:00, 12:00, 15:00) and location within the pond (i.e., bank vs. interior sites and cross-sectional depths at 0, 1, and 2 m). The microbial communities, which were defined by 16S rRNA gene sequencing analysis, significantly varied based on all sampling factors (PERMANOVA P < 0.05 for each). While the relative abundances of dominant phyla (e.g., Proteobacteria and Bacteroidetes) were relatively stable throughout the pond, subtle yet significant increases in α-diversity were observed as the day progressed (ANOVA P < 0.001). Key water quality properties that also increased between the morning and afternoon (i.e., pH, dissolved oxygen, and temperature) positively associated with relative abundances of Cyanobacteria, though were inversely proportional to Verrucomicrobia. These properties, among additional parameters such as bioavailable nutrients (e.g., NH3, NO3, PO4), chlorophyll, phycocyanin, conductivity, and colored dissolved organic matter, exhibited significant relationships with relative abundances of various bacterial genera as well. Further investigation of the microbiota in underlying sediments revealed significant differences between the bank and interior sites of the pond (P < 0.05 for α- and β-diversity). Overall, our findings emphasize the importance of accounting for time of day and water sampling location and depth when surveying the microbiomes of irrigation ponds and other small freshwater sources.

Tributary Inflows to a Regulated River Influence Bacterial Communities and Increase Bacterial Carbon Assimilation

Inflows from unregulated tributaries change the physical, chemical, and biotic conditions in receiving regulated rivers, impacting microbial community structure and metabolic function. Understanding how tributary inflows affect bacterial carbon production (BCP) is integral to understanding energy transfer in riverine ecosystems. To investigate the role of tributary inflows on bacterial community composition and BCP, a ~90th percentile natural flow event was sampled over 5 days along the Lachlan River and its tributaries within the Murray-Darling Basin of eastern Australia. Increased tributary inflows after rainfall corresponded with a significantly different and more diverse bacterial community in the regulated mainstem. The major contributor to this difference was an increase in relative abundance of bacterial groups with a potential metabolic preference for humic substances (Burkholderiaceae Polynucleobacter, Alcaligenaceae GKS98 freshwater group, Saccharimonadia) and a significant decrease in Spirosomaceae Pseudarcicella, known to metabolise algal exudates. Increases in orthophosphate and river discharge explained 31% of community change, suggesting a combination of resource delivery and microbial community coalescence as major drivers. BCP initially decreased significantly with tributary inflows, but the total load of carbon assimilated by bacteria increased by up to 20 times with flow due to increased water volume. The significant drivers of BCP were dissolved organic carbon, water temperature, and conductivity. Notably, BCP was not correlated with bacterial diversity or community composition. Tributary inflows were shown to alter mainstem bacterial community structure and metabolic function to take advantage of fresh terrestrial dissolved organic material, resulting in substantial changes to riverine carbon assimilation over small times scales.

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.

Fourteen new Polynucleobacter species: P. brandtiae sp. nov., P. kasalickyi sp. nov., P. antarcticus sp. nov., P. arcticus sp. nov., P. tropicus sp. nov., P. bastaniensis sp. nov., P. corsicus sp. nov., P. finlandensis sp. nov., P. ibericus sp. nov., P. hallstattensis sp. nov., P. alcilacus sp. nov., P. nymphae sp. nov., P. paludilacus sp. nov. and P. parvulilacunae sp. nov

Fourteen strains, all isolated from the surface of freshwater habitats, were genomically, phylogenetically and phenotypically characterized. The strains were obtained from geographically and climatically broadly scattered sites. This included two lakes in Antarctica, one arctic pond located on the Svalbard archipelago (Norway), a tropical habitat located in Uganda, some lakes in Southern Europe (Spain and France), lakes, ponds and a puddle in Central Europe (Austria, Czech Republic and Germany), and lakes in Northern Europe (Finland). Most of the investigated strains were characterized by rather small cell sizes and rather slow growth on media such as nutrient broth–soyotone–yeast extract (NSY) medium. Phylogenomic analyses indicated that all fourteen strains are affiliated with the genus Polynucleobacter ( Burkholderiaceae , Pseudomonadota ). Thirteen of the strains were found to be affiliated with subcluster PnecC of the genus. All these strains were characterized by genome sizes in the range of 1.7–2.3 Mbp and G+C values of 44.9–46.5 mol%. Furthermore, all PnecC-affiliated strains shared 16S rRNA gene sequence similarities >99 %. Only one strain characterized by a larger genome size of 2.9 Mbp and a lower G+C value of 41.0 mol% was found to be affiliated with subcluster PnecA. Whole genome sequence comparisons revealed that all 14 strains shared among each other, as well as with the type strains of the previously described 17 Polynucleobacter species, whole-genome average nucleotide identities values <95 %. This suggested that the 14 investigated strains represent 14 different new species. We propose the establishment of 14 new Polynucleobacter species represented by the following type strains: UB-Domo-W1T (=DSM 103491T=CIP 111598T=JCM 32562T), VK13T (=DSM 103488T=JCM 32564T), LimPoW16T (=DSM 24085T=CIP 111098T), UK-Long2-W17T (=DSM 103489T=CIP 111328T=JCM 32563T), UK-Pondora-W15T (=DSM 103423T=JCM 32939T), MWH-Mekk-B1T (=DSM 106779T=JCM 32556T), MWH-Mekk-C3T (=DSM 103415T=JCM 32557T), Ross1-W9T (=DSM 103416T=JCM 32561T), MWH-Hall10T (=DSM 107042T=JCM 32938T), AP-Basta-1000A-D1T (=DSM 107039T=JCM 32933T), AP-Melu-1000-A1T (=DSM 107036T=JCM 32935T), es-MAR-2T (=DSM 103424T=JCM 32554T), AP-Mumm-500A-B3T (=DSM 107037T=JCM 32936T), MWH-UH21BT (=DSM 23884T=LMG 29707T).

The Light-to-Nutrient Ratio in Alpine Lakes: Different Scenarios of Bacterial Nutrient Limitation and Community Structure in Lakes Above and Below the Treeline

The light-to-nutrient hypothesis proposes that under high light-to-nutrient conditions, bacteria tend to be limited by phosphorus (P), while under relatively low light-to-nutrient conditions, bacteria are likely driven towards carbon (C) limitation. Exploring whether this light-to-nutrient hypothesis is fitting for alpine lakes has profound implications for predicting the impacts of climatic and environmental changes on the structures and processes of aquatic ecosystems in climate-sensitive regions. We investigated the environmental conditions and bacterioplankton community compositions of 15 high-elevation lakes (7 above and 8 below treeline). High light-to-nutrient conditions (denoted by the reciprocal value of the attenuation coefficient (1/K) to total phosphorus (TP)), high chlorophyll a (Chl a) concentrations, low TP concentrations and low ratios of the dissolved organic carbon concentration to the dissolved total nitrogen concentration (DOC:DTN) were detected in above-treeline lakes. Significant positive correlations between the bacterioplankton community compositions with 1/K:TP ratios and Chl a concentrations indicated that not only high light energy but also nutrient competition between phytoplankton and bacteria could induce P limitation for bacteria. In contrast, low light-to-nutrient conditions and high allochthonous DOC input in below-treeline lakes lessen P limitation and C limitation. The most abundant genus, Polynucleobacter, was significantly enriched, and more diverse oligotypes of Polynucleobacter operational taxonomic units were identified in the below-treeline lakes, indicating the divergence of niche adaptations among Polynucleobacter oligotypes. The discrepancies in the light-to-P ratio and the components of organic matter between the above-treeline and below-treeline lakes have important implications for the nutrient limitation of bacterioplankton and their community compositions.

Cultivation of Dominant Freshwater Bacterioplankton Lineages Using a High-Throughput Dilution-to-Extinction Culturing Approach Over a 1-Year Period

Although many culture-independent molecular analyses have elucidated a great diversity of freshwater bacterioplankton, the ecophysiological characteristics of several abundant freshwater bacterial groups are largely unknown due to the scarcity of cultured representatives. Therefore, a high-throughput dilution-to-extinction culturing (HTC) approach was implemented herein to enable the culture of these bacterioplankton lineages using water samples collected at various seasons and depths from Lake Soyang, an oligotrophic reservoir located in South Korea. Some predominant freshwater bacteria have been isolated from Lake Soyang via HTC (e.g., the acI lineage); however, large-scale HTC studies encompassing different seasons and water depths have not been documented yet. In this HTC approach, bacterial growth was detected in 14% of 5,376 inoculated wells. Further, phylogenetic analyses of 16S rRNA genes from a total of 605 putatively axenic bacterial cultures indicated that the HTC isolates were largely composed of Actinobacteria, Bacteroidetes, Alphaproteobacteria, Betaproteobacteria, Gammaproteobacteria, and Verrucomicrobia. Importantly, the isolates were distributed across diverse taxa including phylogenetic lineages that are widely known cosmopolitan and representative freshwater bacterial groups such as the acI, acIV, LD28, FukuN57, MNG9, and TRA3-20 lineages. However, some abundant bacterial groups including the LD12 lineage, Chloroflexi, and Acidobacteria could not be domesticated. Among the 71 taxonomic groups in the HTC isolates, representative strains of 47 groups could either form colonies on agar plates or be revived from frozen glycerol stocks. Additionally, season and water depth significantly affected bacterial community structure, as demonstrated by 16S rRNA gene amplicon sequencing analyses. Therefore, our study successfully implemented a dilution-to-extinction cultivation strategy to cultivate previously uncultured or underrepresented freshwater bacterial groups, thus expanding the basis for future multi-omic studies.

Opening a next-generation black box: Ecological trends for hundreds of species-like taxa uncovered within a single bacterial >99% 16S rRNA operational taxonomic unit

Current knowledge on environmental distribution and taxon richness of free-living bacteria is mainly based on cultivation-independent investigations employing 16S rRNA gene sequencing methods. Yet, 16S rRNA genes are evolutionarily rather conserved, resulting in limited taxonomic and ecological resolutions provided by this marker. The faster evolving protein-encoding gene priB was used to reveal ecological patterns hidden within a single operational taxonomic unit (OTU) defined by >99% 16S rRNA sequence similarity. The studied subcluster PnecC of the genus Polynucleobacter represents a ubiquitous group of abundant freshwater bacteria with cosmopolitan distribution, which is very frequently detected by diversity surveys of freshwater systems. Based on genome taxonomy and a large set of genome sequences, a sequence similarity threshold for delineation of species-like taxa could be established. In total, 600 species-like taxa were detected in 99 freshwater habitats scattered across three regions representing a latitudinal range of 3,400 km (42°N to 71°N) and a pH gradient of 4.2 to 8.6. In addition to the unexpectedly high richness, the increased taxonomic resolution revealed structuring of Polynucleobacter communities by a couple of macroecological trends, which was previously only demonstrated for phylogenetically much broader groups of bacteria. An unexpected pattern was the almost complete compositional separation of Polynucleobacter communities of Ca²⁺ -rich and Ca²⁺ -poor habitats. This compositional pattern strongly resembled the vicariance of plant species on silicate and limestone soils. The new cultivation-independent approach presented opened a window to an incredible, previously unseen diversity, and enables investigations aiming on deeper understanding of how environmental conditions shape bacterial communities and drive evolution of free-living bacteria.

Continental-Scale Gene Flow Prevents Allopatric Divergence of Pelagic Freshwater Bacteria

Allopatric divergence is one of the principal mechanisms for speciation of macro-organisms. Microbes by comparison are assumed to disperse more freely and to be less limited by dispersal barriers. However, thermophilic prokaryotes restricted to geothermal springs have shown clear signals of geographic isolation, but robust studies on this topic for microbes with less strict habitat requirements are scarce. Furthermore, it has only recently been recognized that homologous recombination among conspecific individuals provides species coherence in a wide range of prokaryotes. Recombination barriers thus may define prokaryotic species boundaries, yet, the extent to which geographic distance between populations gives rise to such barriers is an open question. Here, we investigated gene flow and population structure in a widespread species of pelagic freshwater bacteria, Polynucleobacter paneuropaeus. Through comparative genomics of 113 conspecific strains isolated from freshwater lakes and ponds located across a North–South range of more than 3,000 km, we were able to reconstruct past gene flow events. The species turned out to be highly recombinogenic as indicated by significant signs of gene transfer and extensive genome mosaicism. Although genomic differences increased with spatial distance on a regional scale (<170 km), such correlations were mostly absent on larger scales up to 3,400 km. We conclude that allopatric divergence in European P. paneuropaeus is minor, and that effective gene flow across the sampled geographic range in combination with a high recombination efficacy maintains species coherence.

Ecological Differentiation in Two Major Freshwater Bacterial Taxa Along Environmental Gradients

Polynucleobacter (Burkholderiaceae, Betaproteobacteria) and Limnohabitans (Comamonadaceae, Betaproteobacteria) are abundant freshwater bacteria comprising large genetic and taxonomic diversities, with species adapted to physico-chemically distinct types of freshwater systems. The relative importance of environmental drivers, i.e., physico-chemistry, presence of microeukaryotes and geographic position for the diversity and prevalence has not been investigated for both taxa before. Here, we present the first pan-European study on this topic, comprising 255 freshwater lakes. We investigated Limnohabitans and Polynucleobacter using an amplicon sequencing approach of partial 16S rRNA genes along environmental gradients. We show that physico-chemical factors had the greatest impact on both genera. Analyses on environmental gradients revealed an exceptionally broad ecological spectrum of operational taxonomic units (OTUs). Despite the coarse resolution of the genetic marker, we found OTUs with contrasting environmental preferences within Polynucleobacter and Limnohabitans subclusters. Such an ecological differentiation has been characterized for PnecC and LimC before but was so far unknown for less well studied subclusters such as PnecA and PnecB. Richness and abundance of OTUs are geographically clustered, suggesting that geographic diversity patterns are attributable to region-specific physico-chemical characteristics (e.g., pH and temperature) rather than latitudinal gradients or lake sizes.

Disparate effects of antibiotic-induced microbiome change and enhanced fitness in Daphnia magna

It is a common view that an organism’s microbiota has a profound influence on host fitness; however, supporting evidence is lacking in many organisms. We manipulated the gut microbiome of Daphnia magna by chronic exposure to different concentrations of the antibiotic Ciprofloxacin (0.01–1 mg L^-1), and evaluated whether this affected the animals fitness and antioxidant capacity. In line with our expectations, antibiotic exposure altered the microbiome in a concentration-dependent manner. However, contrary to these expectations, the reduced diversity of gut bacteria was not associated with any fitness detriment. Moreover, the growth-related parameters correlated negatively with microbial diversity; and, in the daphnids exposed to the lowest Ciprofloxacin concentrations, the antioxidant capacity, growth, and fecundity were even higher than in control animals. These findings suggest that Ciprofloxacin exerts direct stimulatory effects on growth and reproduction in the host, while microbiome- mediated effects are of lesser importance. Thus, although microbiome profiling of Daphnia may be a sensitive tool to identify early effects of antibiotic exposure, disentangling direct and microbiome-mediated effects on the host fitness is not straightforward.

Co-occurrence Networks Among Bacteria and Microbial Eukaryotes of Lake Baikal During a Spring Phytoplankton Bloom

The pelagic zone of Lake Baikal is an ecological niche where phytoplankton bloom causes increasing microbial abundance in spring which plays a key role in carbon turnover in the freshwater lake. Co-occurrence patterns revealed among different microbes can be applied to predict interactions between the microbes and environmental conditions in the ecosystem. We used 454 pyrosequencing of 16S rRNA and 18S rRNA genes to study bacterial and microbial eukaryotic communities and their co-occurrence patterns at the pelagic zone of Lake Baikal during a spring phytoplankton bloom. We found that microbes within one domain mostly correlated positively with each other and are highly interconnected. The highly connected taxa in co-occurrence networks were operational taxonomic units (OTUs) of Actinobacteria, Bacteroidetes, Alphaproteobacteria, and autotrophic and unclassified Eukaryota which might be analogous to microbial keystone taxa. Constrained correspondence analysis revealed the relationships of bacterial and microbial eukaryotic communities with geographical location.

Polynucleobacter paneuropaeus sp. nov., characterized by six strains isolated from freshwater lakes located along a 3000 km north-south cross-section across Europe

Six Polynucleobacter (Burkholderiaceae, Betaproteobacteria) strains isolated from different freshwater lakes located across Europe were taxonomically investigated. Phylogenetic analyses based on 16S rRNA gene sequences assigns all six strains to the cryptic species complex PnecC within the genus Polynucleobacter. Analyses of partial glutamine synthetase (glnA) genes suggests that all six strains belong to the species-like taxon designated F15 in previous papers. Comparative genome analyses reveal that the six strains form a genomically coherent group characterized by whole-genome average nucleotide identity (gANI) values of >98 % but separated by gANI values of <88 % from the type strains and representatives of the 16 previously described Polynucleobacter species. In phylogenetic analyses based on nucleotide sequences of 319 orthologous genes, the six strains represent a monophyletic cluster that is clearly separated from all other described species. Genome sizes of the six strains range from 1.61 to 1.83 Mbp, which is smaller than genome sizes of the majority of type strains representing previously described Polynucleobacter species. By contrast, the G+C content of the DNA of the strains is well in the range of 44.8-46.6 mol% previously found for other type strains of species affiliated with the subgroup PnecC. Variation among the six strains representing the new species is evident in a number of traits. These include gene content differences, for instance regarding a gene cluster encoding anoxygenic photosynthesis, as well as phenotypic traits. We propose to name the new species represented by the six strains Polynucleobacter paneuropaeus sp. nov. and designate strain MG-25-Pas1-D2^T (=DSM 103454^T =CIP 111323^T) as the type strain.

Cryptophyta as major bacterivores in freshwater summer plankton

Small bacterivorous eukaryotes play a cardinal role in aquatic food webs and their taxonomic classification is currently a hot topic in aquatic microbial ecology. Despite increasing interest in their diversity, core questions regarding predator–prey specificity remain largely unanswered, e.g., which heterotrophic nanoflagellates (HNFs) are the main bacterivores in freshwaters and which prokaryotes support the growth of small HNFs. To answer these questions, we fed natural communities of HNFs from Římov reservoir (Czech Republic) with five different bacterial strains of the ubiquitous betaproteobacterial genera Polynucleobacter and Limnohabitans. We combined amplicon sequencing and catalyzed reporter deposition fluorescence in situ hybridization (CARD-FISH) targeting eukaryotic 18 S rRNA genes to track specific responses of the natural HNF community to prey amendments. While amplicon sequencing provided valuable qualitative data and a basis for designing specific probes, the number of reads was insufficient to accurately quantify certain eukaryotic groups. We also applied a double-hybridization technique that allows simultaneous phylogenetic identification of both predator and prey. Our results show that community composition of HNFs is strongly dependent upon prey type. Surprisingly, Cryptophyta were the most abundant bacterivores, although this phylum has been so far assumed to be mainly autotrophic. Moreover, the growth of a small lineage of Cryptophyta (CRY1 clade) was strongly stimulated by one Limnohabitans strain in our experiment. Thus, our study is the first report that colorless Cryptophyta are major bacterivores in summer plankton samples and can play a key role in the carbon transfer from prokaryotes to higher trophic levels.

Differences in planktonic microbial communities associated with three types of macrophyte stands in a shallow lake

Little is known about how various substances from living and decomposing aquatic macrophytes affect the horizontal patterns of planktonic bacterial communities. Study sites were located within Lake Kolon, which is a freshwater marsh and can be characterised by open-water sites and small ponds with different macrovegetation (Phragmites australis, Nymphea alba and Utricularia vulgaris). Our aim was to reveal the impact of these macrophytes on the composition of the planktonic microbial communities using comparative analysis of environmental parameters, microscopy and pyrosequencing data. Bacterial 16S rRNA gene sequences were dominated by members of phyla Proteobacteria (36%-72%), Bacteroidetes (12%-33%) and Actinobacteria (5%-26%), but in the anoxic sample the ratio of Chlorobi (54%) was also remarkable. In the phytoplankton community, Cryptomonas sp., Dinobryon divergens, Euglena acus and chrysoflagellates had the highest proportion. Despite the similarities in most of the measured environmental parameters, the inner ponds had different bacterial and algal communities, suggesting that the presence and quality of macrophytes directly and indirectly controlled the composition of microbial plankton.

Metagenomic analysis of soil and freshwater from zoo agricultural area with organic fertilization

Microbial communities drive biogeochemical cycles in agricultural areas by decomposing organic materials and converting essential nutrients. Organic amendments improve soil quality by increasing the load of essential nutrients and enhancing the productivity. Additionally, fresh water used for irrigation can affect soil quality of agricultural soils, mainly due to the presence of microbial contaminants and pathogens. In this study, we investigated how microbial communities in irrigation water might contribute to the microbial diversity and function of soil. Whole-metagenomic sequencing approaches were used to investigate the taxonomic and the functional profiles of microbial communities present in fresh water used for irrigation, and in soil from a vegetable crop, which received fertilization with organic compost made from animal carcasses. The taxonomic analysis revealed that the most abundant genera were Polynucleobacter (~8% relative abundance) and Bacillus (~10%) in fresh water and soil from the vegetable crop, respectively. Low abundance (0.38%) of cyanobacterial groups were identified. Based on functional gene prediction, denitrification appears to be an important process in the soil community analysed here. Conversely, genes for nitrogen fixation were abundant in freshwater, indicating that the N-fixation plays a crucial role in this particular ecosystem. Moreover, pathogenicity islands, antibiotic resistance and potential virulence related genes were identified in both samples, but no toxigenic genes were detected. This study provides a better understanding of the community structure of an area under strong agricultural activity with regular irrigation and fertilization with an organic compost made from animal carcasses. Additionally, the use of a metagenomic approach to investigate fresh water quality proved to be a relevant method to evaluate its use in an agricultural ecosystem.

The Limnohabitans Genus Harbors Generalistic and Opportunistic Subtypes: Evidence from Spatiotemporal Succession in a Canyon-Shaped Reservoir

We studied the diversity of Limnohabitans using reverse line blot hybridization with Limnohabitans lineage-specific probes in the freshwater canyon-shaped Římov reservoir (Czech Republic). To examine the succession of distinct lineages, we performed (i) a study of an intensive spring sampling program at the lacustrine part of the Římov reservoir (from ice melt through a phytoplankton peak to the clear-water phase), and (ii) a seasonal study (April to November) when the occurrence of distinct Limnohabitans lineages was related to the inherent longitudinal heterogeneity of the reservoir. Significant spatiotemporal changes in the compositions of distinct Limnohabitans lineages allowed for the identification of "generalists" that were always present throughout the whole season as well as "specialists" that appeared in the reservoir only for limited periods of time or irregularly. Our results indicate that some phytoplankton groups, such as cryptophytes or cyanobacteria, and zooplankton composition were the major factors modulating the distribution and dynamics of distinct Limnohabitans lineages. The highest Limnohabitans diversity was observed during the spring algal bloom, whereas the lowest was during the summer cyanobacterial bloom. The microdiversity also markedly increased upstream in the reservoir, being highest at the inflow, and thus likely reflecting strong influences of the watershed.IMPORTANCE The genus Limnohabitans is a typical freshwater bacterioplankton and is believed to play a significant role in inland freshwater habitats. This work is unique in detecting and tracing different closely related lineages of this bacterial genus in its natural conditions using the semiquantitative reverse line blot hybridization method and in discovering the factors influencing the microdiversity, subtype alternations, and seasonality.

Genomic divergence and cohesion in a species of pelagic freshwater bacteria

Background

In many prokaryotic genera a clustered phylogeny is observed, akin to the occurrence of species in sexually reproducing organisms. For some taxa, homologous recombination has been invoked as the underlying mechanism providing genomic cohesion among conspecific individuals. Whether this mechanism is applicable to prokaryotes in freshwaters with low habitat connectivity – i.e. elevated geographic barriers to gene flow – is unclear. To investigate further we studied genomic trends within the globally abundant PnecC cluster (genus Polynucleobacter, Betaproteobacteria) and analyzed homologous recombination within the affiliated species P. asymbioticus.

Results

Comparisons among 20 PnecC genomes revealed a clearly discontinuous distribution of nucleotide sequence similarities. Among the nine conspecific individuals (P. asymbioticus) all average nucleotide identity (ANI) values were greater than 97%, whereas all other comparisons exhibited ANI values lower than 85%. The reconstruction of recombination and mutation events for the P. asymbioticus core genomes yielded an r/m ratio of 7.4, which is clearly above estimated thresholds for recombination to act as a cohesive force. Hotspots of recombination were found to be located in the flanking regions of genomic islands. Even between geographically separated habitats a high flux of recombination was evident. While a biogeographic population structure was suggested from MLST data targeting rather conserved loci, such a structure was barely visible when whole genome data was considered. However, both MLST and whole genome data showed evidence of differentiation between two lineages of P. asymbioticus. The ratios of non-synonymous to synonymous substitution rates as well as growth rates in transplantation experiments suggested that this divergence was not selectively neutral.

Conclusions

The high extent of homologous recombination among P. asymbioticus bacteria can act as a cohesive force that effectively counteracts genetic divergence. At least on a regional scale, homologous recombination can act across geographically separated ecosystems and therefore plays an important role in the evolution and consistency of bacterial freshwater species. A species model akin to the biological species concept may be applicable for P. asymbioticus. Nonetheless, two genetically distinct lineages have emerged and further research may clarify if their divergence has been initiated by reinforced geographical barriers or has been evolving in sympatry.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4199-z) contains supplementary material, which is available to authorized users.

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.

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.

Distribution and ecological preferences of the freshwater lineage LimA (genus Limnohabitans) revealed by a new double hybridization approach

The ecological relevance and factors shaping dynamics of Limnohabitans sp. have been largely studied by fluorescence in situ hybridization with a 16S rRNA probe targeting the R-BT group (lineages LimBCDE), but not lineage LimA. Consequently, ecology and distribution of LimA remained unknown. We developed a double hybridization strategy using a novel 23S rRNA probe specifically targeting LimA and LimE that in combination with the existing R-BT probe can discriminate LimA populations. This technique was applied for more than 1000 samples from 46 freshwater sites including long-term data sets from oligo-mesotrophic Lake Zurich, CH and meso-eutrophic Římov reservoir, CZ. LimA was ubiquitously distributed and highly abundant. Observed ecological preferences of LimA in Lake Zurich were in general similar to already reported for Limnohabitans with highest numbers in surface waters during growing seasons. Three times higher densities of LimA were detected in Římov reservoir, where they were significantly more abundant at the riverine zone especially after flood events that introduced fresh terrestrial DOM (dissolved organic matter). Moreover, statistical analyses of biological and physicochemical parameters obtained from small dynamic water bodies confirmed a correspondence between LimA and allochthonous DOM, in opposite to R-BT that was more related to algal primary production.

Microdiversification of a Pelagic Polynucleobacter Species Is Mainly Driven by Acquisition of Genomic Islands from a Partially Interspecific Gene Pool

Microdiversification of a planktonic freshwater bacterium was studied by comparing 37 Polynucleobacter asymbioticus strains obtained from three geographically separated sites in the Austrian Alps. Genome comparison of nine strains revealed a core genome of 1.8 Mb, representing 81% of the average genome size. Seventy-five percent of the remaining flexible genome is clustered in genomic islands (GIs). Twenty-four genomic positions could be identified where GIs are potentially located. These positions are occupied strain specifically from a set of 28 GI variants, classified according to similarities in their gene content. One variant, present in 62% of the isolates, encodes a pathway for the degradation of aromatic compounds, and another, found in 78% of the strains, contains an operon for nitrate assimilation. Both variants were shown in ecophysiological tests to be functional, thus providing the potential for microniche partitioning. In addition, detected interspecific horizontal exchange of GIs indicates a large gene pool accessible to Polynucleobacter species. In contrast to core genes, GIs are spread more successfully across spatially separated freshwater habitats. The mobility and functional diversity of GIs allow for rapid evolution, which may be a key aspect for the ubiquitous occurrence of Polynucleobacter bacteria.

IMPORTANCE

Assessing the ecological relevance of bacterial diversity is a key challenge for current microbial ecology. The polyphasic approach which was applied in this study, including targeted isolation of strains, genome analysis, and ecophysiological tests, is crucial for the linkage of genetic and ecological knowledge. Particularly great importance is attached to the high number of closely related strains which were investigated, represented by genome-wide average nucleotide identities (ANI) larger than 97%. The extent of functional diversification found on this narrow phylogenetic scale is compelling. Moreover, the transfer of metabolically relevant genomic islands between more distant members of the Polynucleobacter community provides important insights toward a better understanding of the evolution of these globally abundant freshwater bacteria.

Acidification and warming affect prominent bacteria in two seasonal phytoplankton bloom mesocosms

In contrast to clear stimulatory effects of rising temperature, recent studies of the effects of CO₂ on planktonic bacteria have reported conflicting results. To better understand the potential impact of predicted climate scenarios on the development and performance of bacterial communities, we performed bifactorial mesocosm experiments (pCO₂ and temperature) with Baltic Sea water, during a diatom dominated bloom in autumn and a mixed phytoplankton bloom in summer. The development of bacterial community composition (BCC) followed well-known algal bloom dynamics. A principal coordinate analysis (PCoA) of bacterial OTUs (operational taxonomic units) revealed that phytoplankton succession and temperature were the major variables structuring the bacterial community whereas the impact of pCO₂ was weak. Prokaryotic abundance and carbon production, and organic matter concentration and composition were partly affected by temperature but not by increased pCO₂ . However, pCO₂ did have significant and potentially direct effects on the relative abundance of several dominant OTUs; in some cases, these effects were accompanied by an antagonistic impact of temperature. Our results suggest the necessity of high-resolution BCC analyses and statistical analyses at the OTU level to detect the strong impact of CO₂ on specific bacterial groups, which in turn might also influence specific organic matter degradation processes.

Microbial Community Profile and Water Quality in a Protected Area of the Caatinga Biome

The Caatinga is a semi-arid biome in northeast Brazil. The Paraguaçú River is located in the Caatinga biome, and part of its course is protected by the National Park of Chapada Diamantina (PNCD). In this study we evaluated the effect of PNCD protection on the water quality and microbial community diversity of this river by analyzing water samples obtained from points located inside and outside the PNCD in both wet and dry seasons. Results of water quality analysis showed higher levels of silicate, ammonia, particulate organic carbon, and nitrite in samples from the unprotected area compared with those from protected areas. Pyrosequencing of the 16S rRNA genes revealed that Burkholderiales was abundant in samples from all three sites during both seasons and was represented primarily by the genus Polynucleobacter and members of the Comamonadaceae family (e.g., genus Limnohabitans). During the dry season, the unprotected area showed a higher abundance of Flavobacterium sp. and Arthrobacter sp., which are frequently associated with the presence and/or degradation of arsenic and pesticide compounds. In addition, genes that appear to be related to agricultural impacts on the environment, as well as those involved in arsenic and cadmium resistance, copper homeostasis, and propanediol utilization, were detected in the unprotected areas by metagenomic sequencing. Although PNCD protection improves water quality, agricultural activities around the park may affect water quality within the park and may account for the presence of bacteria capable of pesticide degradation and assimilation, evidencing possible anthropogenic impacts on the Caatinga.

Distributions of Virus-Like Particles and Prokaryotes within Microenvironments

Microbial interactions are important for ecosystem function, but occur at the microscale and so are difficult to observe. Previous studies in marine systems have shown significant shifts in microbial community abundance and composition over scales of micrometres to centimetres. This study investigates the microscale abundance distributions of virus-like particles (VLPs) and prokaryotes in the lower reaches of a river to determine the extent to which microscale microbial patchiness exists in freshwater systems. Here we report local hotspots surrounded by gradients that reach a maximum 80 and 107 fold change in abundance over 0.9 cm for prokaryotic and VLP subpopulations. Changes in prokaryotic and VLP hotspots were tightly coupled. There were no gradients at tens of centimetres across the boundary layers, which is consistent with strong mixing and turbulence-driven aggregation found in river systems. Quantification of the patchiness shows a marked asymmetry with patches 10 times greater than background common, but depletions being rare or absent in most samples. This consistent asymmetry suggests that coldspots depleted by grazing and lysis are rapidly mixed to background concentrations, while the prevalence of hotspots indicates persistence against disruption. The hotspot to coldspot relative abundance may be useful for understanding microbial river dynamics. The patchiness indicates that the mean- field approach of bulk phase sampling misses the microbially relevant community variation and may underestimate the concentrations of these important microbial groups.

Prey-Specific Growth Responses of Freshwater Flagellate Communities Induced by Morphologically Distinct Bacteria from the Genus Limnohabitans

Because their large growth potential is counterbalanced with grazing by heterotrophic nanoflagellates (HNF), bacteria of the genus Limnohabitans, which are common in many freshwater habitats, represent a valuable model for examining bacterial carbon flow to the grazer food chain. We conducted experiments with natural HNF communities taken from two distinct habitats, the meso-eutrophic Římov Reservoir and the oligo-mesotrophic Lake Cep (South Bohemia). HNF communities from each habitat at distinct seasonal phases, a late April algal bloom and a late May clear water phase, were each fed 3 Limnohabitans strains of differing cell sizes. Water samples were prefiltered (5 μm) to release natural HNF communities from zooplankton control and then amended with the Limnohabitans strains L. planktonicus II-D5 (medium sized, rod shaped), Limnohabitans sp. strain T6-5 (thin, long, curved rod), and Limnohabitans sp. strain 2KL-3 (large solenoid). Using temporal sampling and prey treatment, we determined HNF growth parameters such as doubling time, growth efficiency, and length of lag phase prior starting to exponential growth. All three Limnohabitans strains supported HNF growth but in significant prey-, site-, and season-dependent fashions. For instance, addition of the moderately large T6-5 strain yielded very rapid HNF growth with a short lag phase. In contrast, the curved morphology and larger cell size of strain 2KL-3 made this prey somewhat protected against grazing by smaller HNF, resulting in slower HNF growth and longer lag phases. These trends were particularly pronounced during the late May clear-water phase, which was dominated by smaller HNF cells. This may indicate a longer "adaptation time" for the flagellate communities toward the large prey size offered.

pH influences the importance of niche-related and neutral processes in lacustrine bacterioplankton assembly

pH is an important factor that shapes the structure of bacterial communities. However, we have very limited information about the patterns and processes by which overall bacterioplankton communities assemble across wide pH gradients in natural freshwater lakes. Here, we used pyrosequencing to analyze the bacterioplankton communities in 25 discrete freshwater lakes in Denmark with pH levels ranging from 3.8 to 8.8. We found that pH was the key factor impacting lacustrine bacterioplankton community assembly. More acidic lakes imposed stronger environmental filtering, which decreased the richness and evenness of bacterioplankton operational taxonomic units (OTUs) and largely shifted community composition. Although environmental filtering was determined to be the most important determinant of bacterioplankton community assembly, the importance of neutral assembly processes must also be considered, notably in acidic lakes, where the species (OTU) diversity was low. We observed that the strong effect of environmental filtering in more acidic lakes was weakened by the enhanced relative importance of neutral community assembly, and bacterioplankton communities tended to be less phylogenetically clustered in more acidic lakes. In summary, we propose that pH is a major environmental determinant in freshwater lakes, regulating the relative importance and interplay between niche-related and neutral processes and shaping the patterns of freshwater lake bacterioplankton biodiversity.

Bacterial community composition in three freshwater reservoirs of different alkalinity and trophic status

In order to investigate the factors controlling the bacterial community composition (BCC) in reservoirs, we sampled three freshwater reservoirs with contrasted physical and chemical characteristics and trophic status. The BCC was analysed by 16S rRNA gene amplicon 454 pyrosequencing. In parallel, a complete dataset of environmental parameters and phytoplankton community composition was also collected. BCC in the analysed reservoirs resembled that of epilimnetic waters of natural freshwater lakes with presence of Actinobacteria, Alpha- and Betaproteobacteria, Cytophaga–Flavobacteria–Bacteroidetes (CFB) and Verrucomicrobia groups. Our results evidenced that the retrieved BCC in the analysed reservoirs was strongly influenced by pH, alkalinity and organic carbon content, whereas comparatively little change was observed among layers in stratified conditions.

Global phylogeography of pelagic Polynucleobacter bacteria: restricted geographic distribution of subgroups, isolation by distance and influence of climate

The free-living planktonic freshwater bacterium Polynucleobacter necessarius subspecies asymbioticus (> 99% 16S rRNA similarity) represents a taxon with a cosmopolitan distribution and apparently ubiquitous occurrence in lentic freshwater habitats. We tested for intra-taxon biogeographic patterns by combining cultivation-independent and cultivation methods. A culture collection of 204 strains isolated from globally distributed freshwater habitats (Arctic to Antarctica) was investigated for phylogeographic patterns based on sequences of two markers, the 16S–23S internal transcribed spacers and the glutamine synthetase gene (glnA). Genetic distance between isolates showed significant geographic distance-decay patterns for both markers, suggesting that an isolation-by-distance mechanism influences the global phylogeography. Furthermore, a couple of subgroups showed restricted geographic distributions. Strains of one subgroup were exclusively obtained from tropical sites on four continents (pantropical subgroup). Cultivation-independent methods were used to confirm the restricted geographic distributions of two subgroups. The pantropical taxon could be detected in 63% of investigated tropical habitats but not in any of 121 European freshwater samples. Physiological tests indicated that almost all strains of the pantropical subgroup failed to grow at temperatures of 4°C, while strains affiliated with other subgroups showed good growth at this temperature. This suggests that thermal adaptation is involved in phylogeographic structuring of the global Polynucleobacter population.

Phytoplankton succession affects the composition of Polynucleobacter subtypes in humic lakes

Phytoplankton influence the composition of bacterial communities, but the taxonomic specificity of algal-bacterial interactions is unclear due to the aggregation of ecologically distinct bacterial populations by community characterization methods. Here we examine whether phytoplankton seasonal succession affects the composition of subtypes within the cosmopolitan freshwater bacterial genus Polynucleobacter. Changes in the composition of Polynucleobacter subtypes were characterized in samples collected weekly from May to August in 2003 and 2008 from three humic lakes using terminal restriction fragment length polymorphism fingerprinting of the protein-encoding cytochrome c oxidase ccoN gene. Changes in phytoplankton population abundances explained, on average, 30% of temporal variation in the composition of Polynucleobacter subtypes and the interaction between phytoplankton and the environment explained an additional 18% of temporal variation. The effect of phytoplankton on specific Polynucleobacter subtypes was experimentally confirmed by changes in Polynucleobacter subtype composition following incubation with different phytoplankton assemblages or a no-phytoplankton control. Phytoplankton-associated subtypes and differentiation in substrate use among subtypes likely contribute to the effects of phytoplankton on Polynucleobacter subtype composition. Interactions between unique Polynucleobacter populations and phytoplankton highlight the ecological significance and specificity of species interactions in freshwater communities.

Differential freshwater flagellate community response to bacterial food quality with a focus on Limnohabitans bacteria

Different bacterial strains can have different value as food for heterotrophic nanoflagellates (HNF), thus modulating HNF growth and community composition. We examined the influence of prey food quality using four Limnohabitans strains, one Polynucleobacter strain and one freshwater actinobacterial strain on growth (growth rate, length of lag phase and growth efficiency) and community composition of a natural HNF community from a freshwater reservoir. Pyrosequencing of eukaryotic small subunit rRNA amplicons was used to assess time-course changes in HNF community composition. All four Limnohabitans strains and the Polynucleobacter strain yielded significant HNF community growth while the actinobacterial strain did not although it was detected in HNF food vacuoles. Notably, even within the Limnohabitans strains we found significant prey-related differences in HNF growth parameters, which could not be related only to size of the bacterial prey. Sequence data characterizing the HNF communities showed also that different bacterial prey items induced highly significant differences in community composition of flagellates. Generally, Stramenopiles dominated the communities and phylotypes closely related to Pedospumella (Chrysophyceae) were most abundant bacterivorous flagellates rapidly reacting to addition of the bacterial prey of high food quality.

Patterns of Limnohabitans microdiversity across a large set of freshwater habitats as revealed by Reverse Line Blot Hybridization

Among abundant freshwater Betaproteobacteria, only few groups are considered to be of central ecological importance. One of them is the well-studied genus Limnohabitans and mainly its R-BT subcluster, investigated previously mainly by fluorescence in situ hybridization methods. We designed, based on sequences from a large Limnohabitans culture collection, 18 RLBH (Reverse Line Blot Hybridization) probes specific for different groups within the genus Limnohabitans by targeting diagnostic sequences on their 16 S–23 S rRNA ITS regions. The developed probes covered in sum 92% of the available isolates. This set of probes was applied to environmental DNA originating from 161 different European standing freshwater habitats to reveal the microdiversity (intra-genus) patterns of the Limnohabitans genus along a pH gradient. Investigated habitats differed in various physicochemical parameters, and represented a very broad range of standing freshwater habitats. The Limnohabitans microdiversity, assessed as number of RLBH-defined groups detected, increased significantly along the gradient of rising pH of habitats. 14 out of 18 probes returned detection signals that allowed predictions on the distribution of distinct Limnohabitans groups. Most probe-defined Limnohabitans groups showed preferences for alkaline habitats, one for acidic, and some seemed to lack preferences. Complete niche-separation was indicated for some of the probe-targeted groups. Moreover, bimodal distributions observed for some groups of Limnohabitans, suggested further niche separation between genotypes within the same probe-defined group. Statistical analyses suggested that different environmental parameters such as pH, conductivity, oxygen and altitude influenced the distribution of distinct groups. The results of our study do not support the hypothesis that the wide ecological distribution of Limnohabitans bacteria in standing freshwater habitats results from generalist adaptations of these bacteria. Instead, our observations suggest that the genus Limnohabitans, as well as its R-BT subgroup, represent ecologically heterogeneous taxa, which underwent pronounced ecological diversification.