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

Aquirufa esocilacus sp. nov., Aquirufa originis sp. nov., Aquirufa avitistagni, and Aquirufa echingensis sp. nov. discovered in small freshwater habitats in Austria during a citizen science project

Within a citizen science project, 112 freshwater habitats in Austria were sampled to get bacterial cultures belonging to the genus Aquirufa using a strategy for targeted isolation. We focused on these bacteria because they are widespread and represent typical freshwater bacteria and, furthermore, the typic red pigmentation facilitates preselection. Among the 113 obtained Aquirufa strains were HETE-83DT, KTFRIE-69FT, OSTEICH-129VT and PLAD-142S6KT, originating from small ponds and a creek. Phylogenetic reconstructions with 16S rRNA gene sequences and genome-based analyses with amino acid sequences of 501 core genes showed that all four strains belonged to the A. antheringensis branch of the genus Aquirufa. Calculation of whole-genome average nucleotide identity values and digital DNA-DNA hybridization values revealed that they represent in each case a new species. The genome sizes of the four strains were between 2.5 and 2.8 Mbp and the G + C values were between 41.4 and 41.8%. Like all type strains of the genus Aquirufa, cells were rod-shaped, and liquid cultures and colonies on agar plates were red-pigmented, likely due to carotenoids. All strains except OSTEICH-129VT showed gliding motility on soft agar plates. All strains grew aerobically but only PLAD-142S6KT could grow weakly under anaerobic conditions. We propose here to establish the names Aquirufa esocilacus sp. nov. for strain HETE-83DT (= DSM 118087T = JCM 37094T), Aquirufa originis sp. nov. for KTFRIE-69FT (= DSM 117798T = JCM 37095T), Aquirufa avitistagni for OSTEICH-129VT (= DSM 118088T = JCM 37100T) and Aquirufa echingensis sp. nov. for PLAD-142S6KT (= DSM 117799T = JCM 37096T).

Extremely halophilic brine community manipulation shows higher robustness of microbiomes inhabiting human-driven solar saltern than naturally driven lake

Hypersaline ecosystems display taxonomically similar assemblages with low diversities and highly dense accompanying viromes. The ecological implications of viral infection on natural microbial populations remain poorly understood, especially at finer scales of diversity. Here, we sought to investigate the influence of changes in environmental physicochemical conditions and viral predation pressure by autochthonous and allochthonous viruses on host dynamics. For this purpose, we transplanted two microbiomes coming from distant hypersaline systems (solar salterns of Es Trenc in Spain and the thalassohaline lake of Aran-Bidgol lake in Iran), by exchanging the cellular fractions with the sterile-filtered accompanying brines with and without the free extracellular virus fraction. The midterm exposure (1 month) of the microbiomes to the new conditions showed that at the supraspecific taxonomic range, the assemblies from the solar saltern brine more strongly resisted the environmental changes and viral predation than that of the lake. The metagenome-assembled genomes (MAGs) analysis revealed an intraspecific transition at the ecotype level, mainly driven by changes in viral predation pressure, by both autochthonous and allochthonous viruses.

IMPORTANCE

Viruses greatly influence succession and diversification of their hosts, yet the effects of viral infection on the ecological dynamics of natural microbial populations remain poorly understood, especially at finer scales of diversity. By manipulating the viral predation pressure by autochthonous and allochthonous viruses, we uncovered potential phage-host interaction, and their important role in structuring the prokaryote community at an ecotype level.

Genome-resolved metagenomics reveals the effect of nutrient availability on bacterial genomic properties across 44 European freshwater lakes

Understanding intricate microbial interactions in the environment is crucial. This is especially true for the relationships between nutrients and bacteria, as phosphorus, nitrogen and organic carbon availability are known to influence bacterial population dynamics. It has been suggested that low nutrient conditions prompt the evolutionary process of genome streamlining. This process helps conserve scarce nutrients and allows for proliferation. Genome streamlining is associated with genomic properties such as %GC content, genes encoding sigma factors, percent coding regions, gene redundancy, and functional shifts in processes like cell motility and ATP binding cassette transporters, among others. The current study aims to unveil the impact of nutrition on the genome size, %GC content, and functional properties of pelagic freshwater bacteria. We do this at finer taxonomic resolutions for many metagenomically characterized communities. Our study confirms the interplay of trophic level and genomic properties. It also highlights that different nutrient types, particularly phosphorus and nitrogen, impact these properties differently. We observed a covariation of functional traits with genome size. Larger genomes exhibit enriched pathways for motility, environmental interaction, and regulatory genes. ABC transporter genes reflect the availability of nutrients in the environment, with small genomes presumably relying more on metabolites from other organisms. We also discuss the distinct strategies different phyla adopt to adapt to oligotrophic environments. The findings contribute to our understanding of genomic adaptations within complex microbial communities.

Water quality and microecosystem of water tanks in karst mountainous area, Southwest China

Karst mountainous areas in Southwest China, the world's largest bare karst area, are faced with growing water shortages. Rainwater harvesting plays an important role in alleviating water shortage. However, there remains a substantial gap in the research regarding the water quality of tanks. Water samples were seasonally collected from ten tanks to investigate the physicochemical properties, microbial communities, and their key influencing factors. The result showed that pH, turbidity, chroma, DOC, and CODMn exceeded drinking water guidelines. The alkaline pH value and the deterioration of sensory properties was the main feature of tank water, from which the over-standard rate of the uncleaned water tanks was higher. Moreover, principal component analyses suggested that tank water quality was influenced by human activities, catchment areas, and material cycling processes within the tanks, of which in-tank microbial activities were the most important driving factors in water quality variation. Proteobacteria, Actinobacteria, Bacteroidetes, Cyanobacteria, Firmicutes, and Verrucomicrobia were the predominant bacterial phyla in water tanks. Acinetobacter, Cyanobium-PCC-6307, CL500-29-marine-group, Candidatus-Aquiluna, and Exiguobacterium were the most abundant genera. The bacterial communities were significantly affected by the management practices. Higher relative abundance of Cyanobacteria and lower relative abundance of Proteobacteria was detected in the uncleaned tanks, which was a sign of tank water quality deterioration. The microbial community structure was closely related to the environmental factors. There was evidence that the water quality was affected by the existence of a microecosystem dominated by photosynthetic microorganisms in the water tanks. In addition, Acinetobacter, Enterobacter, Pseudomonas, and Legionella identified as the potential opportunistic pathogenic genera were frequently detected but the relative abundances except Acinetobacter were low in the tanks. Overall, our findings indicated that management style influences water quality and bacterial communities of tank water.

Description of a new freshwater bacterium Aquirufa regiilacus sp. nov., classification of the genera Aquirufa, Arundinibacter, Sandaracinomonas, and Tellurirhabdus to the family Spirosomataceae, classification of the genus Chryseotalea to the family Fulvivirgaceae and Litoribacter to the family Cyclobacteriaceae, as well as classification of Litoribacter alkaliphilus as a later heterotypic synonym of Litoribacter ruber

Strains LEOWEIH-7CT and LEPPI-3A were isolated from the Leopoldskroner Weiher, a lake located in the city of Salzburg, Austria. 16S rRNA gene similarities and phylogenetic reconstructions with 16S rRNA gene sequences as well as based on genome sequences revealed that the new strains belong to the A. antheringensis branch of the genus Aquirufa. Calculated whole-genome average nucleotide identity (gANI) and digital DNA-DNA hybridization (dDDH) values with the closely related type strains showed that the two strains represent a single new species. The strains grew aerobically and chemoorganotrophically, and the cells were rod shaped, on average 0.8 µm long and 0.3 µm wide, red pigmented and motile by gliding. The genome size of both strains was 2.6 Mbp and the G+C value was 41.9%. The genomes comprised genes predicted for the complete light-harvesting rhodopsin system and various carotenoids. We proposed to establish the name Aquirufa regiilacus sp. nov. for strain LEOWEIH-7CT (=DSM 116390T = JCM 36347T) as the type strain. Strain LEPPI-3A (=DSM 116391 = JCM 36348) also belongs to this new species. The calculated genome-based phylogenetic tree revealed that Aquirufa and some other genera currently allocated in the family Cytophagaceae need a reclassification. Aquirufa, Arundinibacter, Sandaracinomonas, and Tellurirhabdus should be designated to the family Spirosomataceae, the genus Chryseotalea to the family Fulvivirgaceae, and the genus Litoribacter to the family Cyclobacteriaceae. Furthermore, based on calculated gANI and dDDH values, Litoribacter alkaliphilus should be reclassified as a later heterotypic synonym of Litoribacter ruber.

Flexible genomic island conservation across freshwater and marine Methylophilaceae

The evolutionary trajectory of Methylophilaceae includes habitat transitions from freshwater sediments to freshwater and marine pelagial that resulted in genome reduction (genome-streamlining) of the pelagic taxa. However, the extent of genetic similarities in the genomic structure and microdiversity of the two genome-streamlined pelagic lineages (freshwater "Ca. Methylopumilus" and the marine OM43 lineage) has so far never been compared. Here, we analyzed complete genomes of 91 "Ca. Methylopumilus" strains isolated from 14 lakes in Central Europe and 12 coastal marine OM43 strains. The two lineages showed a remarkable niche differentiation with clear species-specific differences in habitat preference and seasonal distribution. On the other hand, we observed a synteny preservation in their genomes by having similar locations and types of flexible genomic islands (fGIs). Three main fGIs were identified: a replacement fGI acting as phage defense, an additive fGI harboring metabolic and resistance-related functions, and a tycheposon containing nitrogen-, thiamine-, and heme-related functions. The fGIs differed in relative abundances in metagenomic datasets suggesting different levels of variability ranging from strain-specific to population-level adaptations. Moreover, variations in one gene seemed to be responsible for different growth at low substrate concentrations and a potential biogeographic separation within one species. Our study provides a first insight into genomic microdiversity of closely related taxa within the family Methylophilaceae and revealed remarkably similar dynamics involving mobile genetic elements and recombination between freshwater and marine family members.

Ecological success of extreme halophiles subjected to recurrent osmotic disturbances is primarily driven by congeneric species replacement

To understand how extreme halophiles respond to recurrent disturbances, we challenged the communities thriving in salt-saturated (~36% salts) ~230 L brine mesocosms to repeated dilutions down to 13% (D13 mesocosm) or 20% (D20 mesocosm) salts each time mesocosms reached salt saturation due to evaporation (for 10 and 17 cycles, respectively) over 813 days. Depending on the magnitude of dilution, the most prevalent species, Haloquadratum walsbyi and Salinibacter ruber, either increased in dominance by replacing less competitive populations (for D20, moderate stress conditions), or severely decreased in abundance and were eventually replaced by other congeneric species better adapted to the higher osmotic stress (for D13, strong stress conditions). Congeneric species replacement was commonly observed within additional abundant genera in response to changes in environmental or biological conditions (e.g. phage predation) within the same system and under a controlled perturbation of a relevant environmental parameter. Therefore, a genus is an ecologically important level of diversity organization, not just a taxonomic rank, that persists in the environment based on congeneric species replacement due to relatively high functional overlap (gene sharing), with important consequences for the success of the lineage, and similar to the success of a species via strain-replacement. Further, our results showed that successful species were typically accompanied by the emergence of their own viral cohorts, whose intra-cohort diversity appeared to strongly covary with, and likely drive, the intra-host diversity. Collectively, our results show that brine communities are ecologically resilient and continuously adapting to changing environments by transitioning to alternative stable states.

Exploring environmental intra-species diversity through non-redundant pangenome assemblies

At the genome level, microorganisms are highly adaptable both in terms of allele and gene composition. Such heritable traits emerge in response to different environmental niches and can have a profound influence on microbial community dynamics. As a consequence, any individual genome or population will contain merely a fraction of the total genetic diversity of any operationally defined "species", whose ecological potential can thus be only fully understood by studying all of their genomes and the genes therein. This concept, known as the pangenome, is valuable for studying microbial ecology and evolution, as it partitions genomes into core (present in all the genomes from a species, and responsible for housekeeping and species-level niche adaptation among others) and accessory regions (present only in some, and responsible for intra-species differentiation). Here we present SuperPang, an algorithm producing pangenome assemblies from a set of input genomes of varying quality, including metagenome-assembled genomes (MAGs). SuperPang runs in linear time and its results are complete, non-redundant, preserve gene ordering and contain both coding and non-coding regions. Our approach provides a modular view of the pangenome, identifying operons and genomic islands, and allowing to track their prevalence in different populations. We illustrate this by analysing intra-species diversity in Polynucleobacter, a bacterial genus ubiquitous in freshwater ecosystems, characterized by their streamlined genomes and their ecological versatility. We show how SuperPang facilitates the simultaneous analysis of allelic and gene content variation under different environmental pressures, allowing us to study the drivers of microbial diversification at unprecedented resolution.

Bacterioplankton Communities in Dissolved Organic Carbon-Rich Amazonian Black Water

The Amazon River basin sustains dramatic hydrochemical gradients defined by three water types: white, clear, and black waters. In black water, important loads of allochthonous humic dissolved organic matter (DOM) result from the bacterioplankton degradation of plant lignin. However, the bacterial taxa involved in this process remain unknown, since Amazonian bacterioplankton has been poorly studied. Its characterization could lead to a better understanding of the carbon cycle in one of the Earth's most productive hydrological systems. Our study characterized the taxonomic structure and functions of Amazonian bacterioplankton to better understand the interplay between this community and humic DOM. We conducted a field sampling campaign comprising 15 sites distributed across the three main Amazonian water types (representing a gradient of humic DOM), and a 16S rRNA metabarcoding analysis based on bacterioplankton DNA and RNA extracts. Bacterioplankton functions were inferred using 16S rRNA data in combination with a tailored functional database from 90 Amazonian basin shotgun metagenomes from the literature. We discovered that the relative abundances of fluorescent DOM fractions (humic-, fulvic-, and protein-like) were major drivers of bacterioplankton structure. We identified 36 genera for which the relative abundance was significantly correlated with humic DOM. The strongest correlations were found in the Polynucleobacter, Methylobacterium, and Acinetobacter genera, three low abundant but omnipresent taxa that possessed several genes involved in the main steps of the β-aryl ether enzymatic degradation pathway of diaryl humic DOM residues. Overall, this study identified key taxa with DOM degradation genomic potential, the involvement of which in allochthonous Amazonian carbon transformation and sequestration merits further investigation. IMPORTANCE The Amazon basin discharge carries an important load of terrestrially derived dissolved organic matter (DOM) to the ocean. The bacterioplankton from this basin potentially plays important roles in transforming this allochthonous carbon, which has consequences on marine primary productivity and global carbon sequestration. However, the structure and function of Amazonian bacterioplanktonic communities remain poorly studied, and their interactions with DOM are unresolved. In this study, we (i) sampled bacterioplankton in all the main Amazon tributaries, (ii) combined information from the taxonomic structure and functional repertory of Amazonian bacterioplankton communities to understand their dynamics, (iii) identified the main physicochemical parameters shaping bacterioplanktonic communities among a set of >30 measured environmental parameters, and (iv) characterized how bacterioplankton structure varies according to the relative abundance of humic compounds, a by-product from the bacterial degradation process of allochthonous DOM.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Adaptive genetic traits in pelagic freshwater microbes

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

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

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

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

Aquirufa lenticrescens sp. nov. and Aquirufa aurantiipilula sp. nov.: two new species of a lineage of widespread freshwater bacteria

Two bacterial strains, 9H-EGSE^T and 15D-MOB^T, were isolated from small freshwater habitats located near Salzburg, Austria. They showed the highest 16S rRNA sequence similarities of 100% and 99.9%, respectively, with type strains of species of the genus Aquirufa (Bacteroidota). Genome-based phylogenetic reconstructions with 119 amino acid sequences assigned the new taxa to the two distinct branches of the genus Aquirufa. Whole-genome average nucleotide identities were calculated with all possible pairs belonging to the genus. Values between 75.4% and 88.6% revealed that the two new strains represent each a new species. Like all, so far described members of the genus, they grew aerobically and chemoorganotrophically, were rod-shaped, red-pigmented, and motile by gliding, and showed genome sizes of about 3 Mbp and G + C values of about 40%. They could be distinguished by some phenotypic and chemotaxonomic features from their nearest related species. Until now, strain 9H-EGSE^T is the only one among the Aquirufa strains which contained traces of MK8 as respiratory quinone, and strain 15D-MOB^T is the only one that formed tiny orange globules in liquid medium. The genome of strain 9H-EGSE^T comprised genes for the complete light-harvesting rhodopsin / retinal system, in the case of 15D-MOB^T genes predicted for a nitrous oxide reductase were present. For the two new species of the genus Aquirufa, we propose to establish the names Aquirufa lenticrescens for strain 9H-EGSE^T (= JCM 34077 ^T = CIP 111926 ^T) and Aquirufa aurantiipilula for strain 15D-MOB^T (= JCM 34078 ^T = CIP 111925 ^T).

Water Column Microbial Communities Vary along Salinity Gradients in the Florida Coastal Everglades Wetlands

Planktonic microbial communities mediate many vital biogeochemical processes in wetland ecosystems, yet compared to other aquatic ecosystems, like oceans, lakes, rivers or estuaries, they remain relatively underexplored. Our study site, the Florida Everglades (USA)—a vast iconic wetland consisting of a slow-moving system of shallow rivers connecting freshwater marshes with coastal mangrove forests and seagrass meadows—is a highly threatened model ecosystem for studying salinity and nutrient gradients, as well as the effects of sea level rise and saltwater intrusion. This study provides the first high-resolution phylogenetic profiles of planktonic bacterial and eukaryotic microbial communities (using 16S and 18S rRNA gene amplicons) together with nutrient concentrations and environmental parameters at 14 sites along two transects covering two distinctly different drainages: the peat-based Shark River Slough (SRS) and marl-based Taylor Slough/Panhandle (TS/Ph). Both bacterial as well as eukaryotic community structures varied significantly along the salinity gradient. Although freshwater communities were relatively similar in both transects, bacterioplankton community composition at the ecotone (where freshwater and marine water mix) differed significantly. The most abundant taxa in the freshwater marshes include heterotrophic Polynucleobacter sp. and potentially phagotrophic cryptomonads of the genus Chilomonas, both of which could be key players in the transfer of detritus-based biomass to higher trophic levels.

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.

Aquiluna borgnonia gen. nov., sp. nov., a member of a Microbacteriaceae lineage of freshwater bacteria with small genome sizes

The actinobacterial strain 15G-AUS-rot^T was isolated from an artificial pond located near Salzburg, Austria. The strain showed 16S rRNA gene sequence similarities of 98.7 % to Candidatus Aquiluna rubra and of 96.6 and 96.7 % to the two validly described species of the genus Rhodoluna. Phylogenetic reconstructions based on 16S rRNA gene sequences and genome-based on amino acid sequences of 118 single copy genes referred strain 15G-AUS-rot^T to the family Microbacteriaceae and therein to the so-called subcluster Luna-1. The genome-based phylogenetic tree showed that the new strain represents a putative new genus. Cultures of strain 15G-AUS-rot^T were light red pigmented and comprised very small, rod-shaped cells. They metabolized a broad variety of substrates. Major fatty acids (>10 %) of cells were iso-C_16 : 0, antiso-C_15 : 0 and iso-C_14 : 0. The major respiratory quinone was MK-11 and a minor component was MK-10. The peptidoglycan structure belonged to an unusual B type. The closed genome sequence of the strain was very small (1.4 Mbp) and had a DNA G+C content of 54.8 mol%. An interesting feature was the presence of genes putatively encoding the complete light-driven proton pumping actinorhodopsin/retinal system, which were located at three different positions of the genome. Based on the characteristics of the strain, a new genus and a new species termed Aquiluna borgnonia is proposed for strain 15G-AUS-rot^T (=DSM 107803^T=JCM 32974^T).

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.

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.

Martelella alba sp. nov., isolated from mangrove rhizosphere soil within the Beibu Gulf

Strain BGMRC 2036^T was isolated from rhizosphere soil of Bruguiear gymnorrhiza collected from the Beibu Gulf of China. Optimum growth occurred at 28 °C, pH 7.0, and under the conditions of 3-5% (w/v) NaCl. The phylogenetic comparisons of 16S rRNA gene sequences displayed that strain BGMRC 2036^T was closely related to Martelella limonii NBRC109441^T (96.6% sequence similarity), M. mediterranea CGMCC 1.12224^T (96.5%), M. lutilitoris GH2-6^T (96.5%), M. radicis BM5-7^T (96.2%), and M. mangrove BM9-1^T (95.9%), M. suaedae NBRC109440^T (95.8%). The phylogenomic tree based on the up-to-date bacterial core gene set indicated that the strain BGMRC 2036^T form a clade formed with members of the genera Martelella. The major polar lipids include phosphatidylmethylethanolamine, phosphatidylglycerol, phosphatidylcholine, phosphotidylinositol, two unidentified phospholipids, and three unidentified ninhydrin positive phospholipids. The major respiratory quinone is Q-10, which is similar to those of genera Martelella. The main cellular fatty acids are C_18:1 ω7c, C_16:0, and C_12:0 aldehyde. Genome sequencing revealed a genome size of 4.99 Mbp and a G + C content of 62.3 mol%. Pairwise comparison of the genomes of the new strain BGMRC 2036^T and the three reference strains M. endophytica YC 6887^T, M. mediterranea CGMCC 1.12224^T, and M. mangrovi USBA-857 indicated that gANI value was lower than 81% and a digital DNA-DNA hybridization value was lower than 27%. The strain BGMRC 2036^T possessed genes putatively encoding riboflavin synthesis and flavodoxin A polyphasic taxonomic study suggested that strain BGMRC 2036^T represented a novel species belonging to the genus Martelella, and it was named Martelella alba sp. nov. The type strain is BGMRC 2036^T (=KCTC 52121^T =NBRC 111908^T).

Combined Methylome, Transcriptome and Proteome Analyses Document Rapid Acclimatization of a Bacterium to Environmental Changes

Polynucleobacter asymbioticus strain QLW-P1DMWA-1^T represents a group of highly successful heterotrophic ultramicrobacteria that is frequently very abundant (up to 70% of total bacterioplankton) in freshwater habitats across all seven continents. This strain was originally isolated from a shallow Alpine pond characterized by rapid changes in water temperature and elevated UV radiation due to its location at an altitude of 1300 m. To elucidate the strain’s adjustment to fluctuating environmental conditions, we recorded changes occurring in its transcriptomic and proteomic profiles under contrasting experimental conditions by simulating thermal conditions in winter and summer as well as high UV irradiation. To analyze the potential connection between gene expression and regulation via methyl group modification of the genome, we also analyzed its methylome. The methylation pattern differed between the three treatments, pointing to its potential role in differential gene expression. An adaptive process due to evolutionary pressure in the genus was deduced by calculating the ratios of non-synonymous to synonymous substitution rates for 20 Polynucleobacter spp. genomes obtained from geographically diverse isolates. The results indicate purifying selection.

Mycoavidus sp. Strain B2-EB: Comparative Genomics Reveals Minimal Genomic Features Required by a Cultivable Burkholderiaceae-Related Endofungal Bacterium

Obligate bacterial endosymbionts are critical to the existence of many eukaryotes. Such endobacteria are usually characterized by reduced genomes and metabolic dependence on the host, which may cause difficulty in isolating them in pure cultures. Family Burkholderiaceae-related endofungal bacteria affiliated with the Mycoavidus-Glomeribacter clade can be associated with the fungal subphyla Mortierellomycotina and Glomeromycotina. In this study, a cultivable endosymbiotic bacterium, Mycoavidus sp. strain B2-EB, present in the fungal host Mortierella parvispora was obtained successfully. The B2-EB genome (1.88 Mb) represents the smallest genome among the endofungal bacterium Mycoavidus cysteinexigens (2.64-2.80 Mb) of Mortierella elongata and the uncultured endosymbiont "Candidatus Glomeribacter gigasporarum" (1.37 to 2.36 Mb) of arbuscular mycorrhizal fungi. Despite a reduction in genome size, strain B2-EB displays a high genome completeness, suggesting a nondegenerative reduction in the B2-EB genome. Compared with a large proportion of transposable elements (TEs) in other known Mycoavidus genomes (7.2 to 11.5% of the total genome length), TEs accounted for only 2.4% of the B2-EB genome. This pattern, together with a high proportion of single-copy genes in the B2-EB genome, suggests that the B2-EB genome reached a state of relative evolutionary stability. These results represent the most streamlined structure among the cultivable endofungal bacteria and suggest the minimal genome features required by both an endofungal lifestyle and artificial culture. This study allows us to understand the genome evolution of Burkholderiaceae-related endosymbionts and to elucidate microbiological interactions.IMPORTANCE This study attempted the isolation of a novel endobacterium, Mycoavidus sp. B2-EB (JCM 33615), harbored in the fungal host Mortierella parvispora E1425 (JCM 39028). We report the complete genome sequence of this strain, which possesses a reduced genome size with relatively high genome completeness and a streamlined genome structure. The information indicates the minimal genomic features required by both the endofungal lifestyle and artificial cultivation, which furthers our understanding of genome reduction in fungal endosymbionts and extends the culture resources for biotechnological development on engineering synthetic microbiomes.

Cryptic speciation of a pelagic Roseobacter population varying at a few thousand nucleotide sites

A drop of seawater contains numerous microspatial niches at the scale relevant to microbial activities. Examples are abiotic niches such as detrital particles that show different sizes and organic contents, and biotic niches resulting from bacteria-phage and bacteria-phytoplankton interactions. A common practice to investigate the impact of microenvironments on bacterial evolution is to separate the microenvironments physically and compare the bacterial inhabitants from each. It remains poorly understood, however, which microenvironment primarily drives bacterioplankton evolution in the pelagic ocean. By applying a dilution cultivation approach to an undisturbed coastal water sample, we isolate a bacterial population affiliated with the globally dominant Roseobacter group. Although varying at just a few thousand nucleotide sites across the whole genomes, members of this clonal population are diverging into two genetically separated subspecies. Genes underlying speciation are not unique to subspecies but instead clustered at the shared regions that represent ~6% of the genomic DNA. They are primarily involved in vitamin synthesis, motility, oxidative defense, carbohydrate, and amino acid utilization, consistent with the known strategies that roseobacters take to interact with phytoplankton and particles. Physiological assays corroborate that one subspecies outcompetes the other in these traits. Our results indicate that the microenvironments in the pelagic ocean represented by phytoplankton and organic particles are likely important niches that drive the cryptic speciation of the Roseobacter population, though microhabitats contributed by other less abundant pelagic hosts cannot be ruled out.

Aquirufa ecclesiirivi sp. nov. and Aquirufa beregesia sp. nov., isolated from a small creek and classification of Allopseudarcicella aquatilis as a later heterotypic synonym of Aquirufa nivalisilvae

Two bacterial strains, 50A-KIRBA^T and 50C-KIRBA^T, were isolated from the same freshwater creek located near Salzburg, Austria. They showed 16S rRNA gene sequence similarities to Aquirufa nivalisilvae of 100 and 99.9 %, respectively. A genome-based phylogenetic reconstruction with amino acid sequences of 119 single-copy genes suggested that the new strains represent two new species of the genus Aquirufa. Pairwise calculated whole-genome average nucleotide identity (gANI) values ranging from 85.4 to 87.5 % confirmed this conclusion. Phenotypic, chemotaxonomic and genomic traits were investigated. Like strains of other Aquirufa species, 50A-KIRBA^T and 50C-KIRBA^T grew aerobically and chemoorganotrophically, were rod-shaped, red-pigmented and motile, most likely by gliding. They could be distinguished by slight differences in the chemotaxonomic features. We propose to establish for strain 50A-KIRBA^T (=CIP 111735^T=LMG 31080^T) as type strain the name Aquirufa ecclesiirivi and for strain 50C-KIRBA^T (=CIP 111736^T=LMG 31501^T) as type strain the name Aquirufa beregesia. Furthermore, the relationship between the type strains of Aquirufa nivalisilvae (59G-WUEMPEL^T) and Allopseudarcicella aquatilis (HME7025^T) was investigated. Results of polyphasic analyses, especially a gANI value of 97.6 %, as well as the genome-based phylogenetic reconstruction, suggested that Allopseudarcicella aquatilis is a heterotypic synonym of Aquirufa nivalisilvae. According to rule 24b of the International Code of Nomenclature of Prokaryotes we propose to classify strain HME7025 as Aquirufa nivalisilvae and provide an emended description for the latter.

Fluviispira multicolorata gen. nov., sp. nov. and Silvanigrella paludirubra sp. nov., isolated from freshwater habitats

Strain 33A1-SZDP^T was isolated from a small creek located in Puch, Austria. Strain SP-Ram-0.45-NSY-1^T was obtained from a small pond located in Schönramer Moor, Germany. 16S rRNA gene sequence similarities between the type strain of Silvanigrella aquatica, currently the only member of the family Silvanigrellaceae, and strains 33A1-SZDP^T and SP-Ram-0.45-NSY-1^T of 94.1 and 99.1 %, respectively, suggested affiliation of the two strains with this family. Phylogenetic reconstructions with 16S rRNA gene sequences and phylogenomic analyses with amino acid sequences obtained from 103 single-copy genes suggested that the strains represent a new genus and a new species in the case of strain 33A1-SZDP^T (=JCM 32978^T=DSM 107810^T), and a new species within the genus Silvanigrella in the case of strain SP-Ram-0.45-NSY-1^T (=JCM 32975^T=DSM 107809^T). Cells of strain 33A1-SZDP^T were motile, pleomorphic, purple-pigmented on agar plates, putatively due to violacein, and showed variable pigmentation in liquid media. They grew chemoorganotrophically and aerobically and tolerated salt concentrations up to 1.2 % NaCl (v/w). The genome size of strain 33A1-SZDP^T was 3.4 Mbp and the G+C content was 32.2 mol%. For this new genus and new species, we propose the name Fluviispira multicolorata gen. nov., sp. nov. Cells of strain SP-Ram-0.45-NSY-1^T were motile, pleomorphic, red-pigmented and grew chemoorganotrophically and aerobically. They tolerated salt concentrations up to 1.1 % NaCl (v/w). The genome size of strain SP-Ram-0.45-NSY-1^T was 3.9 Mbp and the G+C content 29.3 mol%. For the new species within the genus Silvanigrella we propose the name Silvanigrella paludirubra sp. nov.

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.

Genomic, metabolic and phenotypic variability shapes ecological differentiation and intraspecies interactions of Alteromonas macleodii

Ecological differentiation between strains of bacterial species is shaped by genomic and metabolic variability. However, connecting genotypes to ecological niches remains a major challenge. Here, we linked bacterial geno- and phenotypes by contextualizing pangenomic, exometabolomic and physiological evidence in twelve strains of the marine bacterium Alteromonas macleodii, illuminating adaptive strategies of carbon metabolism, microbial interactions, cellular communication and iron acquisition. In A. macleodii strain MIT1002, secretion of amino acids and the unique capacity for phenol degradation may promote associations with Prochlorococcus cyanobacteria. Strain 83-1 and three novel Pacific isolates, featuring clonal genomes despite originating from distant locations, have profound abilities for algal polysaccharide utilization but without detrimental implications for Ecklonia macroalgae. Degradation of toluene and xylene, mediated via a plasmid syntenic to terrestrial Pseudomonas, was unique to strain EZ55. Benzoate degradation by strain EC673 related to a chromosomal gene cluster shared with the plasmid of A. mediterranea EC615, underlining that mobile genetic elements drive adaptations. Furthermore, we revealed strain-specific production of siderophores and homoserine lactones, with implications for nutrient acquisition and cellular communication. Phenotypic variability corresponded to different competitiveness in co-culture and geographic distribution, indicating linkages between intraspecific diversity, microbial interactions and biogeography. The finding of "ecological microdiversity" helps understanding the widespread occurrence of A. macleodii and contributes to the interpretation of bacterial niche specialization, population ecology and biogeochemical roles.

Rariglobus hedericola gen. nov., sp. nov., belonging to the Verrucomicrobia, isolated from a temperate freshwater habitat

The bacterial strain 53C-WASEF was isolated from a small freshwater ditch located in Eugendorf, Austria. Phylogenetic reconstructions with 16S rRNA gene sequences and genome based, with amino acid sequences obtained from 105 single copy genes, suggested that the strain represents a new genus and a new species within the family Opitutaceae, which belongs to the class Opitutae of the phylum Verrucomicrobia. Comparisons of the 16S rRNA gene sequence of strain 53C-WASEF with those of related type strains revealed a highest sequence similarity of 93.5 % to Nibricoccus aquaticus and of 92.9 % to Geminisphaera colitermitum. Interestingly, phylogentic trees indicated the latter as being the closest known relative of the new strain. Phenotypic, chemotaxonomic and genomic traits were investigated. Cells were observed to be small, spherical, motile and unpigmented, and grew chemoorganotrophically and aerobically. The respiratory quinone was MK-7, the predominant fatty acids were anteiso-C_15 : 0, C_16 : 1ω5c and C_16 : 0. The identified polar lipids were phosphatidylethanolamine, phosphatidylglycerol and diphosphatidylglycerol. Genome sequencing revealed genes putatively encoding for flagella synthesis and cellulose degradation. The genome size was 4.1 Mbp and the G+C content 60.6 mol%. For the new genus and the new species, we propose the name Rariglobus hedericola gen. nov., sp. nov. (=CIP 111665^T=DSM 109123^T).

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.

Gene gain and loss push prokaryotes beyond the homologous recombination barrier and accelerate genome sequence divergence

Bacterial and archaeal evolution involve extensive gene gain and loss. Thus, phylogenetic trees of prokaryotes can be constructed both by traditional sequence-based methods (gene trees) and by comparison of gene compositions (genome trees). Comparing the branch lengths in gene and genome trees with identical topologies for 34 clusters of closely related bacterial and archaeal genomes, we show here that terminal branches of gene trees are systematically compressed compared to those of genome trees. Thus, sequence evolution is delayed compared to genome evolution by gene gain and loss. The extent of this delay differs widely among bacteria and archaea. Mathematical modeling shows that the divergence delay can result from sequence homogenization by homologous recombination. The model explains how homologous recombination maintains the cohesiveness of the core genome of a species while allowing extensive gene gain and loss within the accessory genome. Once evolving genomes become isolated by barriers impeding homologous recombination, gene and genome evolution processes settle into parallel trajectories, and genomes diverge, resulting in speciation.

A significant proportion of the molecular evolution of bacteria and archaea occurs through gene gain and loss. Here Iranzo et al. develop a mathematical model that explains observed differential patterns of sequence evolution vs. gene content evolution as a consequence of homologous recombination.

Microdiversity ensures the maintenance of functional microbial communities under changing environmental conditions

Microdiversity can lead to different ecotypes within the same species. These are assumed to provide stability in time and space to those species. However, the role of microdiversity in the stability of whole microbial communities remains underexplored. Understanding the drivers of microbial community stability is necessary to predict community response to future disturbances. Here, we analyzed 16S rRNA gene amplicons from eight different temperate bog lakes at the 97% OTU and amplicon sequence variant (ASV) levels and found ecotypes within the same OTU with different distribution patterns in space and time. We observed that these ecotypes are adapted to different values of environmental factors such as water temperature and oxygen concentration. Our results showed that the existence of several ASVs within a OTU favored its persistence across changing environmental conditions. We propose that microdiversity aids the stability of microbial communities in the face of fluctuations in environmental factors.

Ecogenomic characterization of widespread, closely-related SAR11 clades of the freshwater genus "Candidatus Fonsibacter" and proposal of Ca. Fonsibacter lacus sp. nov

The ubiquitous alpha-proteobacteria of the order "Candidatus Pelagibacterales" (SAR11) are highly abundant in aquatic environments, and among them, members of the monophyletic lineage LD12 (also known as SAR11 clade IIIb) are specifically found in lacustrine ecosystems. Clade IIIb bacteria are some of the most prominent members of freshwater environments, but little is known about their biology due to the lack of genome representatives. Only recently, the first non-marine isolate was cultured and described as "Candidatus Fonsibacter ubiquis". Here, we expand the collection of freshwater IIIb representatives and describe a new IIIb species of the genus "Ca. Fonsibacter". Specifically, we assembled a collection of 67 freshwater metagenomic datasets from the interconnected lakes of the Chattahoochee River basin (GA, USA) and obtained nearly complete metagenome-assembled genomes (MAGs) representing 5 distinct IIIb subclades, roughly equivalent to species based on genomic standards, including the previously described "Ca. F. ubiquis". Genomic comparisons between members of the IIIb species revealed high similarity in gene content. However, when comparing their abundance profiles in the Chattahoochee basin and various aquatic environments, differences in temporal and spatial distributions among the distinct species were observed implying niche differentiation might be underlying the coexistence of the highly functionally similar representatives. The name Ca. Fonsibacter lacus sp. nov. is proposed for the most abundant and widespread species in the Chattahoochee River basin and various freshwater ecosystems.

Comparative analysis reveals the Genomic Islands in Pasteurella multocida population genetics: on Symbiosis and adaptability

BACKGROUND

Pasteurella multocida (P. multocida) is a widespread opportunistic pathogen that infects human and various animals. Genomic Islands (GIs) are one of the most important mobile components that quickly help bacteria acquire large fragments of foreign genes. However, the effects of GIs on P. multocida are unknown in the evolution of bacterial populations.

RESULTS

Ten avian-sourced P. multocida obtained through high-throughput sequencing together with 104 publicly available P. multocida genomes were used to analyse their population genetics, thus constructed a pan-genome containing 3948 protein-coding genes. Through the pan-genome, the open evolutionary pattern of P. multocida was revealed, and the functional components of 944 core genes, 2439 accessory genes and 565 unique genes were analysed. In addition, a total of 280 GIs were predicted in all strains. Combined with the pan-genome of P. multocida, the GIs accounted for 5.8% of the core genes in the pan-genome, mainly related to functional metabolic activities; the accessory genes accounted for 42.3%, mainly for the enrichment of adaptive genes; and the unique genes accounted for 35.4%, containing some defence mechanism-related genes.

CONCLUSIONS

The effects of GIs on the population genetics of P. multocida evolution and adaptation to the environment are reflected by the proportion and function of the pan-genome acquired from GIs, and the large quantities of GI data will aid in additional population genetics studies.

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.

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.

Polynucleobacter hirudinilacicola sp. nov. and Polynucleobacter campilacus sp. nov., both isolated from freshwater systems

Strains MWH-EgelM1-30-B4^T and MWH-Feld-100^T were isolated from the water columns of two freshwater systems. Both strains represent delicate bacteria not easy to work with in laboratory experiments. Phylogenetic analyses of the 16S rRNA genes suggested that both strains were affiliated with the genus Polynucleobacter. Both strains share 16S rRNA gene sequence similarities of >99 % with eight free-living Polynucleobacter type strains, all affiliated with the cryptic species complex PnecC. The full-length 16S rRNA gene sequences of the two strains differ only in two and three positions, respectively, from the sequence of the closest related Polynucleobacter type strain. Genome sequencing of both strains revealed relatively small genome sizes of 2.0 Mbp and G+C contents of 45 mol%. Phylogenetic analyses based on nucleotide sequences of 319 shared protein-encoding genes consistently placed the two strains in taxon PnecC but did not suggest an affiliation with one of the previously described species. Pairwise analyses of whole genome average nucleotide identities (gANI) with representatives of all previously described Polynucleobacter species resulted in both cases throughout in values <80 %. Pairwise comparison of the genomes of the two new strains resulted in gANI values of 83.3 %. All gANI analyses clearly suggested that strains MWH-EgelM1-30-B4^T and MWH-Feld-100^T represent two novel Polynucleobacter species. We propose for these novel species the names Polynucleobacter hirudinilacicola sp. nov. and Polynucleobacter campilacus sp. nov. and strains MWH-EgelM1-30-B4^T (=DSM 23911^T=LMG 30144^T) and MWH-Feld-100^T (=DSM 24007^T=LMG 29705^T) as the type strains, respectively.

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.

Genomes of Novel Microbial Lineages Assembled from the Sub-Ice Waters of Lake Baikal

We present a metagenomic study of Lake Baikal (East Siberia). Two samples obtained from the water column under the ice cover (5 and 20 m deep) in March 2016 have been deep sequenced and the reads assembled to generate metagenome-assembled genomes (MAGs) that are representative of the microbes living in this special environment. Compared with freshwater bodies studied around the world, Lake Baikal had an unusually high fraction of Verrucomicrobia Other groups, such as Actinobacteria and Proteobacteria, were in proportions similar to those found in other lakes. The genomes (and probably cells) tended to be small, presumably reflecting the extremely oligotrophic and cold prevalent conditions. Baikal microbes are novel lineages recruiting very little from other water bodies and are distantly related to other freshwater microbes. Despite their novelty, they showed the closest relationship to genomes discovered by similar approaches from other freshwater lakes and reservoirs. Some of them were particularly similar to MAGs from the Baltic Sea, which, although it is brackish, connected to the ocean, and much more eutrophic, has similar climatological conditions. Many of the microbes contained rhodopsin genes, indicating that, in spite of the decreased light penetration allowed by the thick ice/snow cover, photoheterotrophy could be widespread in the water column, either because enough light penetrates or because the microbes are already adapted to the summer ice-less conditions. We have found a freshwater SAR11 subtype I/II representative showing striking synteny with Pelagibacterubique strains, as well as a phage infecting the widespread freshwater bacterium PolynucleobacterIMPORTANCE Despite the increasing number of metagenomic studies on different freshwater bodies, there is still a missing component in oligotrophic cold lakes suffering from long seasonal frozen cycles. Here, we describe microbial genomes from metagenomic assemblies that appear in the upper water column of Lake Baikal, the largest and deepest freshwater body on Earth. This lake is frozen from January to May, which generates conditions that include an inverted temperature gradient (colder up), decrease in light penetration due to ice, and, especially, snow cover, and oligotrophic conditions more similar to the open-ocean and high-altitude lakes than to other freshwater or brackish systems. As could be expected, most reconstructed genomes are novel lineages distantly related to others in cold environments, like the Baltic Sea and other freshwater lakes. Among them, there was a broad set of streamlined microbes with small genomes/intergenic spacers, including a new nonmarine Pelagibacter-like (subtype I/II) genome.

Reclassification of a Polynucleobacter cosmopolitanus strain isolated from tropical Lake Victoria as Polynucleobacter victoriensis sp. nov

The genus Polynucleobacter (family Burkholderiaceae) is phylogenetically subdivided into at least four subclusters. One of those, subcluster PnecC, was recognized as a cryptic species complex. Here we test by comparative genome analyses whether subcluster PnecD, currently solely represented by the species Polynucleobactercosmopolitanus, also represents such a cryptic species complex. The genome sequences of the two P. cosmopolitanus strains, MWH-MoIso2^T and MWH-VicM1, were determined. The latter strain was also characterized in the previous description of P. cosmopolitanus. These two strains originate from a temperate lake located in Austria and from the large tropical Lake Victoria located in East Africa, respectively. Strains MWH-MoIso2^T and MWH-VicM1 possess quite small genomes of 1.78 and 1.63 Mbp, respectively, and share similar G+C values of 44.1 and 43.1 mol%, respectively. Both strains encode only a single copy of the ribosomal operon, and their 16S rRNA genes differ only in four positions, equalling a sequence similarity of 99.74 %. Both genomes possess characteristics indicating evolutionary genome streamlining, such as high coding densities of 93.9 and 94.6 % of bases, respectively. Average nucleotide identity (ANI) comparisons of the genomes of the two strains resulted in a value of 78.4 %, suggesting that each of the strains represents a separate species. Our investigation suggests that PnecD represents an additional cryptic species complex within the genus Polynucleobacter that was not resolved by 16S rRNA gene sequence analyses. We propose reclassification of strain MWH-VicM1 as Polynucleobacter victoriensis sp. nov., with type strain MWH-VicM1^T(=DSM 21486^T=JCM 32005^T).

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.

Microdiversification in genome-streamlined ubiquitous freshwater Actinobacteria

Actinobacteria of the acI lineage are the most abundant microbes in freshwater systems, but there are so far no pure living cultures of these organisms, possibly because of metabolic dependencies on other microbes. This, in turn, has hampered an in-depth assessment of the genomic basis for their success in the environment. Here we present genomes from 16 axenic cultures of acI Actinobacteria. The isolates were not only of minute cell size, but also among the most streamlined free-living microbes, with extremely small genome sizes (1.2-1.4 Mbp) and low genomic GC content. Genome reduction in these bacteria might have led to auxotrophy for various vitamins, amino acids and reduced sulphur sources, thus creating dependencies to co-occurring organisms (the 'Black Queen' hypothesis). Genome analyses, moreover, revealed a surprising degree of inter- and intraspecific diversity in metabolic pathways, especially of carbohydrate transport and metabolism, and mainly encoded in genomic islands. The striking genotype microdiversification of acI Actinobacteria might explain their global success in highly dynamic freshwater environments with complex seasonal patterns of allochthonous and autochthonous carbon sources. We propose a new order within Actinobacteria ('Candidatus Nanopelagicales') with two new genera ('Candidatus Nanopelagicus' and 'Candidatus Planktophila') and nine new species.

Polynucleobacter aenigmaticus sp. nov. isolated from the permanently anoxic monimolimnion of a temperate meromictic lake

The bacterial strain MWH-K35W1^T was isolated from a permanently anoxic water layer of a meromictic lake located in the Austrian Salzkammergut area. The basically chemo-organoheterotrophic strain was isolated and maintained under aerobic conditions. Phylogenetic analyses of the 16S rRNA gene and the glutamine synthetase gene (glnA) of the strain suggested an affiliation to the genus Polynucleobacter and the cryptic species complex PnecC. Strain MWH-K35W1^T shares with the type strains of the six free-living species of the genus Polynucleobacter affiliated with this species complex 16S rRNA gene sequence similarities of 99.6-99.9 %, while the type material of the obligate endosymbiont Polynucleobacternecessarius, which is also affiliated with this species complex, shares a gene sequence similarity of 99.1 %. Genome sequencing resulted in a genome size of 2.14 Mbp and a DNA G+C content of 45.98 mol%. Major fatty acids were C16 : 1ω7c, C18 : 1ω7c and C16 : 0. This strain is the first strain of the genus Polynucleobacter found to encode a proteorhodopsin-like protein but, in contrast to some other strains affiliated to this genus, it does not encode a putative anoxygenic photosynthesis system. Multilocus sequence analysis based on partial sequences of eight housekeeping genes, as well as average nucleotide identity (ANI) analyses, did not suggest that strain MWH-K35W1^T belongs to a previously described species. We propose the name Polynucleobacter aenigmaticus for a novel species with strain MWH-K35W1^T (=DSM 24006^T=LMG 29706^T) as the type strain.

Parallel genome reduction in symbionts descended from closely related free-living bacteria

Endosymbiosis plays an important role in ecology and evolution, but fundamental aspects of the origin of intracellular symbionts remain unclear. The extreme age of many symbiotic relationships, lack of data on free-living ancestors and uniqueness of each event hinder investigations. Here, we describe multiple strains of the bacterium Polynucleobacter that evolved independently and under similar conditions from closely related, free-living ancestors to become obligate endosymbionts of closely related ciliate hosts. As these genomes reduced in parallel from similar starting states, they provide unique glimpses into the mechanisms underlying genome reduction in symbionts. We found that gene loss is contingently lineage-specific, with no evidence for ordered streamlining. However, some genes in otherwise disrupted pathways are retained, possibly reflecting cryptic genetic network complexity. We also measured substitution rates between many endosymbiotic and free-living pairs for hundreds of genes, which showed that genetic drift, and not mutation pressure, is the main non-selective factor driving molecular evolution in endosymbionts.