Single-cell genomics reveal low recombination frequencies in freshwater bacteria of the SAR11 clade

Freshwater genome-reduced bacteria exhibit pervasive episodes of adaptive stasis

The emergence of bacterial species is rooted in their inherent potential for continuous evolution and adaptation to an ever-changing ecological landscape. The adaptive capacity of most species frequently resides within the repertoire of genes encoding the secreted proteome (SP), as it serves as a primary interface used to regulate survival/reproduction strategies. Here, by applying evolutionary genomics approaches to metagenomics data, we show that abundant freshwater bacteria exhibit biphasic adaptation states linked to the eco-evolutionary processes governing their genome sizes. While species with average to large genomes adhere to the dominant paradigm of evolution through niche adaptation by reducing the evolutionary pressure on their SPs (via the augmentation of functionally redundant genes that buffer mutational fitness loss) and increasing the phylogenetic distance of recombination events, most of the genome-reduced species exhibit a nonconforming state. In contrast, their SPs reflect a combination of low functional redundancy and high selection pressure, resulting in significantly higher levels of conservation and invariance. Our findings indicate that although niche adaptation is the principal mechanism driving speciation, freshwater genome-reduced bacteria often experience extended periods of adaptive stasis. Understanding the adaptive state of microbial species will lead to a better comprehension of their spatiotemporal dynamics, biogeography, and resilience to global change.

Tools for microbial single-cell genomics for obtaining uncultured microbial genomes

The advent of next-generation sequencing technologies has facilitated the acquisition of large amounts of DNA sequence data at a relatively low cost, leading to numerous breakthroughs in decoding microbial genomes. Among the various genome sequencing activities, metagenomic analysis, which entails the direct analysis of uncultured microbial DNA, has had a profound impact on microbiome research and has emerged as an indispensable technology in this field. Despite its valuable contributions, metagenomic analysis is a "bulk analysis" technique that analyzes samples containing a wide diversity of microbes, such as bacteria, yielding information that is averaged across the entire microbial population. In order to gain a deeper understanding of the heterogeneous nature of the microbial world, there is a growing need for single-cell analysis, similar to its use in human cell biology. With this paradigm shift in mind, comprehensive single-cell genomics technology has become a much-anticipated innovation that is now poised to revolutionize microbiome research. It has the potential to enable the discovery of differences at the strain level and to facilitate a more comprehensive examination of microbial ecosystems. In this review, we summarize the current state-of-the-art in microbial single-cell genomics, highlighting the potential impact of this technology on our understanding of the microbial world. The successful implementation of this technology is expected to have a profound impact in the field, leading to new discoveries and insights into the diversity and evolution of microbes.

Geographic population structure and distinct intra-population dynamics of globally abundant freshwater bacteria

Implications of geographic separation and temporal dynamics on the evolution of free-living bacterial species are widely unclear. However, the vast amount of metagenome sequencing data generated during the last decades from various habitats around the world provides an unprecedented opportunity for such investigations. Here, we exploited publicly available and new freshwater metagenomes in combination with the genomes of abundant freshwater bacteria to reveal geographic and temporal population structure. We focused on species that were detected across broad geographic ranges at high enough sequence coverage for meaningful population genomic analyses, associated with the predominant freshwater taxa acI, LD12, Polynucleobacter, and Candidatus Methylopumilus. Despite the broad geographic ranges, each species appeared as a sequence-discrete cluster, in contrast to abundant marine taxa, for which continuous diversity structures were reported on a global scale. Population differentiation increased significantly with spatial distance in all species, but notable dispersal barriers (e.g. oceanic) were not apparent. Yet, the different species showed contrasting rates of geographic divergence and strikingly different intra-population dynamics in time series within individual habitats. The change in an LD12 population over 7 years was minor (FST = 0.04) compared to differentiation between lakes, whereas a Polynucleobacter population displayed strong changes within merely 2 months (FST up to 0.54), similar in scale to differentiation between populations separated by thousands of kilometers. The slowly and steadily evolving LD12 population showed high strain diversity, whereas the dynamic Polynucleobacter population exhibited alternating clonal expansions of mostly two strains only. Based on the contrasting population structures, we propose distinct models of speciation.

Genome reduction occurred in early Prochlorococcus with an unusually low effective population size

In the oligotrophic sunlit ocean, the most abundant free-living planktonic bacterial lineages evolve convergently through genome reduction. The cyanobacterium Prochlorococcus responsible for 10% global oxygen production is a prominent example. The dominant theory known as "genome streamlining" posits that they have extremely large effective population sizes (Ne) such that selection for metabolic efficiency acts to drive genome reduction. Because genome reduction largely took place anciently, this theory builds on the assumption that their ancestors' Ne was similarly large. Constraining Ne for ancient ancestors is challenging because experimental measurements of extinct organisms are impossible and alternatively reconstructing ancestral Ne with phylogenetic models gives large uncertainties. Here, we develop a new strategy that leverages agent-based modeling to simulate the changes in the genome-wide ratio of radical to conservative nonsynonymous nucleotide substitution rate (dR/dC) in a possible range of Ne in ancestral populations. This proxy shows expected increases with decreases of Ne only when Ne falls to about 10 k - 100 k or lower, magnitudes characteristic of Ne of obligate endosymbiont species where drift drives genome reduction. Our simulations therefore strongly support a scenario where the primary force of Prochlorococcus genome reduction is drift rather than selection.

Uncultured prokaryotic genomes in the spotlight: An examination of publicly available data from metagenomics and single-cell genomics

Owing to the ineffectiveness of traditional culture techniques for the vast majority of microbial species, culture-independent analyses utilizing next-generation sequencing and bioinformatics have become essential for gaining insight into microbial ecology and function. This mini-review focuses on two essential methods for obtaining genetic information from uncultured prokaryotes, metagenomics and single-cell genomics. We analyzed the registration status of uncultured prokaryotic genome data from major public databases and assessed the advantages and limitations of both the methods. Metagenomics generates a significant quantity of sequence data and multiple prokaryotic genomes using straightforward experimental procedures. However, in ecosystems with high microbial diversity, such as soil, most genes are presented as brief, disconnected contigs, and lack association of highly conserved genes and mobile genetic elements with individual species genomes. Although technically more challenging, single-cell genomics offers valuable insights into complex ecosystems by providing strain-resolved genomes, addressing issues in metagenomics. Recent technological advancements, such as long-read sequencing, machine learning algorithms, and in silico protein structure prediction, in combination with vast genomic data, have the potential to overcome the current technical challenges and facilitate a deeper understanding of uncultured microbial ecosystems and microbial dark matter genes and proteins. In light of this, it is imperative that continued innovation in both methods and technologies take place to create high-quality reference genome databases that will support future microbial research and industrial applications.

Large-scale phylogenomics of aquatic bacteria reveal molecular mechanisms for adaptation to salinity

The crossing of environmental barriers poses major adaptive challenges. Rareness of freshwater-marine transitions separates the bacterial communities, but how these are related to brackish counterparts remains elusive, as do the molecular adaptations facilitating cross-biome transitions. We conducted large-scale phylogenomic analysis of freshwater, brackish, and marine quality-filtered metagenome-assembled genomes (11,248). Average nucleotide identity analyses showed that bacterial species rarely existed in multiple biomes. In contrast, distinct brackish basins cohosted numerous species, but their intraspecific population structures displayed clear signs of geographic separation. We further identified the most recent cross-biome transitions, which were rare, ancient, and most commonly directed toward the brackish biome. Transitions were accompanied by systematic changes in amino acid composition and isoelectric point distributions of inferred proteomes, which evolved over millions of years, as well as convergent gains or losses of specific gene functions. Therefore, adaptive challenges entailing proteome reorganization and specific changes in gene content constrains the cross-biome transitions, resulting in species-level separation between aquatic biomes.

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

Postglacial adaptations enabled colonization and quasi-clonal dispersal of ammonia-oxidizing archaea in modern European large lakes

Ammonia-oxidizing archaea (AOA) play a key role in the aquatic nitrogen cycle. Their genetic diversity is viewed as the outcome of evolutionary processes that shaped ancestral transition from terrestrial to marine habitats. However, current genome-wide insights into AOA evolution rarely consider brackish and freshwater representatives or provide their divergence timeline in lacustrine systems. An unbiased global assessment of lacustrine AOA diversity is critical for understanding their origins, dispersal mechanisms, and ecosystem roles. Here, we leveraged continental-scale metagenomics to document that AOA species diversity in freshwater systems is remarkably low compared to marine environments. We show that the uncultured freshwater AOA, "Candidatus Nitrosopumilus limneticus," is ubiquitous and genotypically static in various large European lakes where it evolved 13 million years ago. We find that extensive proteome remodeling was a key innovation for freshwater colonization of AOA. These findings reveal the genetic diversity and adaptive mechanisms of a keystone species that has survived clonally in lakes for millennia.

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.

Phylogenetic divergence and adaptation of Nitrososphaeria across lake depths and freshwater ecosystems

Thaumarchaeota (now the class Nitrososphaeria in the phylum Thermoproteota in GTDB taxonomy) are abundant across marine and soil habitats; however, their genomic diversity and evolutionary history in freshwater environments remain elusive. Here, we reconstructed 17 high-quality metagenome-assembled genomes of Nitrososphaeria from a deep lake and two great rivers, and compared all available genomes between freshwater and marine habitats regarding their phylogenetic positions, relative abundance, and genomic content. We found that freshwater Nitrososphaeria were dominated by the family Nitrosopumilaceae and could be grouped into three distinct clades closely related to the genera Nitrosopumilus, Nitrosoarchaeum, and Nitrosotenuis. The Nitrosopumilus-like clade was exclusively from deep lakes, while the Nitrosoarchaeum-like clade was dominated by species from deep lakes and rivers, and the Nitrosotenuis-like clade was mainly from rivers, deep lakes, and estuaries. Interestingly, there was vertical niche separation between two clades in deep lakes, showing that the Nitrosopumilus-like species dominated shallow layers, whereas the relative abundance of the Nitrosoarchaeum-like clade increased toward deep waters. Phylogenetic clustering patterns in the Nitrosopumilaceae supported at least one freshwater-to-marine and two marine-to-freshwater transitions, the former of which refined the potential terrestrial-to-marine evolutionary path as previously proposed. The occurrence of the two marine-to-freshwater transitions were accompanied by horizontal transfer of the genes involved in nutrition regulation, osmoregulation, and cell motility during their colonization to freshwater habitats. Specifically, the Nitrosopumilus-like clade showed losses of genes encoding flagella assembly and ion transport, whereas the Nitrosoarchaeum-like clade had losses of intact genes involved in urea uptake and utilization and gains of genes encoding osmolarity-mediated mechanosensitive channels. Collectively, our results reveal for the first time the high genomic diversity of the class Nitrososphaeria across freshwater ecosystems and provide novel insights into their adaptive mechanisms and evolutionary histories.

Dissecting the dominant hot spring microbial populations based on community-wide sampling at single-cell genomic resolution

With advances in DNA sequencing and miniaturized molecular biology workflows, rapid and affordable sequencing of single-cell genomes has become a reality. Compared to 16S rRNA gene surveys and shotgun metagenomics, large-scale application of single-cell genomics to whole microbial communities provides an integrated snapshot of community composition and function, directly links mobile elements to their hosts, and enables analysis of population heterogeneity of the dominant community members. To that end, we sequenced nearly 500 single-cell genomes from a low diversity hot spring sediment sample from Dewar Creek, British Columbia, and compared this approach to 16S rRNA gene amplicon and shotgun metagenomics applied to the same sample. We found that the broad taxonomic profiles were similar across the three sequencing approaches, though several lineages were missing from the 16S rRNA gene amplicon dataset, likely the result of primer mismatches. At the functional level, we detected a large array of mobile genetic elements present in the single-cell genomes but absent from the corresponding same species metagenome-assembled genomes. Moreover, we performed a single-cell population genomic analysis of the three most abundant community members, revealing differences in population structure based on mutation and recombination profiles. While the average pairwise nucleotide identities were similar across the dominant species-level lineages, we observed differences in the extent of recombination between these dominant populations. Most intriguingly, the creek's Hydrogenobacter sp. population appeared to be so recombinogenic that it more closely resembled a sexual species than a clonally evolving microbe. Together, this work demonstrates that a randomized single-cell approach can be useful for the exploration of previously uncultivated microbes from community composition to population structure.

Genomic features underlying the evolutionary transitions of Apibacter to honey bee gut symbionts

The gut bacteria of honey bee recognized as a mutualistic partner with the insect host might have originated from a free-living or parasitic lifestyle. However, little is known about the genomic features underlying this lifestyle transition. Here we compared the genomes of bee gut bacteria Apibacter with their close relatives living in different lifestyles. We found that despite general reduction in the Apibacter genome, genes involved in amino acid synthesis and monosaccharide detoxification were retained, which is putatively beneficial to the host. Interestingly, the microaerobic Apibacter species specifically acquired genes encoding for the nitrate respiration (NAR). These together with nitrate transporter and enzymatic cofactor synthesis genes were found clustered in the genomes. The NAR system is also conserved in the cohabitating bee gut microbe Snodgrassella, although with a different structure. This convergence suggests a key role of respiratory nitrate reduction for microaerophilic microbiomes to colonize bee gut epithelium. Genes involved in lipid, histidine degradation were found partially or completely lost in Apibacter. Particularly, genes encoding for the conversion to the toxic intermediates in phenylacetate degradation, as well as other potential virulence factors, are specifically lost in Apibacter group. Antibiotic resistance genes are only sporadically distributed among Apibacter species, but are prevalent in their relatives, which may be related to the remotely living feature and less exposure to antibiotics of their bee hosts. Collectively, this study advanced our knowledge of genomic features specialized to bee gut symbionts.

Heterozygous, Polyploid, Giant Bacterium, Achromatium, Possesses an Identical Functional Inventory Worldwide across Drastically Different Ecosystems

Achromatium is large, hyperpolyploid and the only known heterozygous bacterium. Single cells contain approximately 300 different chromosomes with allelic diversity far exceeding that typically harbored by single bacteria genera. Surveying all publicly available sediment sequence archives, we show that Achromatium is common worldwide, spanning temperature, salinity, pH, and depth ranges normally resulting in bacterial speciation. Although saline and freshwater Achromatium spp. appear phylogenetically separated, the genus Achromatium contains a globally identical, complete functional inventory regardless of habitat. Achromatium spp. cells from differing ecosystems (e.g., from freshwater to saline) are, unexpectedly, equally functionally equipped but differ in gene expression patterns by transcribing only relevant genes. We suggest that environmental adaptation occurs by increasing the copy number of relevant genes across the cell's hundreds of chromosomes, without losing irrelevant ones, thus maintaining the ability to survive in any ecosystem type. The functional versatility of Achromatium and its genomic features reveal alternative genetic and evolutionary mechanisms, expanding our understanding of the role and evolution of polyploidy in bacteria while challenging the bacterial species concept and drivers of bacterial speciation.

Microbial Diversity and Phage-Host Interactions in the Georgian Coastal Area of the Black Sea Revealed by Whole Genome Metagenomic Sequencing

Viruses have the greatest abundance and highest genetic diversity in marine ecosystems. The interactions between viruses and their hosts is one of the hot spots of marine ecology. Besides their important role in various ecosystems, viruses, especially bacteriophages and their gene pool, are of enormous interest for the development of new gene products with high innovation value. Various studies have been conducted in diverse ecosystems to understand microbial diversity and phage-host interactions; however, the Black Sea, especially the Eastern coastal area, remains among the least studied ecosystems in this regard. This study was aimed at to fill this gap by analyzing microbial diversity and bacteriophage-host interactions in the waters of Eastern Black Sea using a metagenomic approach. To this end, prokaryotic and viral metagenomic DNA from two sampling sites, Poti and Gonio, were sequenced on the Illumina Miseq platform and taxonomic and functional profiles of the metagenomes were obtained using various bioinformatics tools. Our metagenomics analyses allowed us to identify the microbial communities, with Proteobacteria, Cyanobacteria, Actinibacteria, and Firmicutes found to be the most dominant bacterial phyla and Synechococcus and Candidatus Pelagibacter phages found to be the most dominant viral groups in the Black Sea. As minor groups, putative phages specific to human pathogens were identified in the metagenomes. We also characterized interactions between the phages and prokaryotic communities by determining clustered regularly interspaced short palindromic repeats (CRISPR), prophage-like sequences, and integrase/excisionase sequences in the metagenomes, along with identification of putative horizontally transferred genes in the viral contigs. In addition, in the viral contig sequences related to peptidoglycan lytic activity were identified as well. This is the first study on phage and prokaryote diversity and their interactions in the Eastern coastal area of the Black Sea using a metagenomic approach.

The Evolutionary Success of the Marine Bacterium SAR11 Analyzed through a Metagenomic Perspective

As the most abundant bacteria in oceans, the Pelagibacterales order (here SAR11) plays an important role in the global carbon cycle, but the study of the evolutionary forces driving its evolution has lagged considerably due to the inherent difficulty of obtaining pure cultures. Multiple evolutionary models have been proposed to explain the diversification of distinct lineages within a population; however, the identification of many of these patterns in natural populations remains mostly enigmatic. We have used a metagenomic approach to explore microdiversity patterns in their natural habitats. Comparison with a collection of bacterial and archaeal groups from the same environments shows that SAR11 populations have a different evolutionary regime, where multiple genotypes coexist within the same population and remain stable over time. Widespread homologous recombination could be one of the main driving factors of this homogenization.

The SAR11 clade of Alphaproteobacteria is the most abundant group of planktonic cells in the near-surface epipelagic waters of the ocean, but the mechanisms underlying its exceptional success have not been fully elucidated. Here, we applied a metagenomic approach to explore microdiversity patterns by measuring the accumulation of synonymous and nonsynonymous mutations as well as homologous recombination in populations of SAR11 from different aquatic habitats (marine epipelagic, bathypelagic, and surface freshwater). The patterns of mutation accumulation and recombination were compared to those of other groups of representative marine microbes with multiple ecological strategies that share the same marine habitat, namely, Cyanobacteria (Prochlorococcus and Synechococcus), Archaea (“Candidatus Nitrosopelagicus” and Marine Group II Thalassoarchaea), and some heterotrophic marine bacteria (Alteromonas and Erythrobacter). SAR11 populations showed widespread recombination among distantly related members, preventing divergence leading to a genetically stable population. Moreover, their high intrapopulation sequence diversity with an enrichment in synonymous replacements supports the idea of a very ancient divergence and the coexistence of multiple different clones. However, other microbes analyzed seem to follow different evolutionary dynamics where processes of diversification driven by geographic and ecological instability produce a higher number of nonsynonymous replacements and lower intrapopulation sequence diversity. Together, these data shed light on some of the evolutionary and ecological processes that lead to the large genomic diversity in SAR11. Furthermore, this approach can be applied to other similar microbes that are difficult to culture in the laboratory, but abundant in nature, to investigate the underlying dynamics of their genomic evolution.

IMPORTANCE As the most abundant bacteria in oceans, the Pelagibacterales order (here SAR11) plays an important role in the global carbon cycle, but the study of the evolutionary forces driving its evolution has lagged considerably due to the inherent difficulty of obtaining pure cultures. Multiple evolutionary models have been proposed to explain the diversification of distinct lineages within a population; however, the identification of many of these patterns in natural populations remains mostly enigmatic. We have used a metagenomic approach to explore microdiversity patterns in their natural habitats. Comparison with a collection of bacterial and archaeal groups from the same environments shows that SAR11 populations have a different evolutionary regime, where multiple genotypes coexist within the same population and remain stable over time. Widespread homologous recombination could be one of the main driving factors of this homogenization.

Streamlined and Abundant Bacterioplankton Thrive in Functional Cohorts

This study examines evolutionary and ecological relationships of three of the most ubiquitous and abundant freshwater bacterial genera: “Ca. Planktophila” (acI-A), “Ca. Nanopelagicus” (acI-B), and “Ca. Fonsibacter” (LD12). Due to high abundance, these genera might have a significant influence on nutrient cycling in freshwaters worldwide, and this study adds a layer of understanding to how seemingly competing clades of bacteria can coexist by having different cooperation strategies. Our synthesis ties together network and ecological theory with empirical evidence and lays out a framework for how the functioning of populations within complex microbial communities can be studied.

While fastidious microbes can be abundant and ubiquitous in their natural communities, many fail to grow axenically in laboratories due to auxotrophies or other dependencies. To overcome auxotrophies, these microbes rely on their surrounding cohort. A cohort may consist of kin (ecotypes) or more distantly related organisms (community) with the cooperation being reciprocal or nonreciprocal and expensive (Black Queen hypothesis) or costless (by-product). These metabolic partnerships (whether at single species population or community level) enable dominance by and coexistence of these lineages in nature. Here we examine the relevance of these cooperation models to explain the abundance and ubiquity of the dominant fastidious bacterioplankton of a dimictic mesotrophic freshwater lake. Using both culture-dependent (dilution mixed cultures) and culture-independent (small subunit [SSU] rRNA gene time series and environmental metagenomics) methods, we independently identified the primary cohorts of actinobacterial genera “Candidatus Planktophila” (acI-A) and “Candidatus Nanopelagicus” (acI-B) and the proteobacterial genus “Candidatus Fonsibacter” (LD12). While “Ca. Planktophila” and “Ca. Fonsibacter” had no correlation in their natural habitat, they have the potential to be complementary in laboratory settings. We also investigated the bifunctional catalase-peroxidase enzyme KatG (a common good which “Ca. Planktophila” is dependent upon) and its most likely providers in the lake. Further, we found that while ecotype and community cooperation combined may explain “Ca. Planktophila” population abundance, the success of “Ca. Nanopelagicus” and “Ca. Fonsibacter” is better explained as a community by-product. Ecotype differentiation of “Ca. Fonsibacter” as a means of escaping predation was supported but not for overcoming auxotrophies.

IMPORTANCE This study examines evolutionary and ecological relationships of three of the most ubiquitous and abundant freshwater bacterial genera: “Ca. Planktophila” (acI-A), “Ca. Nanopelagicus” (acI-B), and “Ca. Fonsibacter” (LD12). Due to high abundance, these genera might have a significant influence on nutrient cycling in freshwaters worldwide, and this study adds a layer of understanding to how seemingly competing clades of bacteria can coexist by having different cooperation strategies. Our synthesis ties together network and ecological theory with empirical evidence and lays out a framework for how the functioning of populations within complex microbial communities can be studied.

Metagenomics Uncovers a Core SAR11 Population in Brackish Surface Waters of the Baltic Sea

The Baltic Sea represents one of the largest brackish ecosystems where various environmental factors control dynamic seasonal shifts in the structure, diversity, and function of the planktonic microbial communities. In this study, despite seasonal fluctuations, several bacterial populations (<2% of the total OTUs) that are highly dominant (25% of relative abundance) and highly frequently occurring (>85% of occurrence) over four seasons were identified. Mathematical models using occurrence frequency and relative abundance data were able to describe community assembly persisting over time. Further, this work uncovered one of the core bacterial populations phylogenetically affiliated to SAR11 subclade IIIa. The analysis of the hypervariable region of 16S rRNA gene and single copy housekeeping genes recovered from metagenomic datasets suggested that the population was unexpectedly evolutionarily closely related to those inhabiting a mesosaline lacustrine ecosystem rather than other marine/coastal members. Our metagenomic results further revealed that the newly-identified population was the major driver facilitating the seasonal shifts in the overall community structure over the brackish waters of the Baltic Sea. The core community uncovered in this study supports the presence of a brackish water microbiome distinguishable from other marine and freshwater counterparts and will be a useful sentinel for monitoring local/global environmental changes posed on brackish surface waters.

Ecogenomics of the SAR11 clade

Members of the SAR11 clade, despite their high abundance, are often poorly represented by metagenome‐assembled genomes. This fact has hampered our knowledge about their ecology and genetic diversity. Here we examined 175 SAR11 genomes, including 47 new single‐amplified genomes. The presence of the first genomes associated with subclade IV suggests that, in the same way as subclade V, they might be outside the proposed Pelagibacterales order. An expanded phylogenomic classification together with patterns of metagenomic recruitment at a global scale have allowed us to define new ecogenomic units of classification (genomospecies), appearing at different, and sometimes restricted, metagenomic data sets. We detected greater microdiversity across the water column at a single location than in samples collected from similar depth across the global ocean, suggesting little influence of biogeography. In addition, pangenome analysis revealed that the flexible genome was essential to shape genomospecies distribution. In one genomospecies preferentially found within the Mediterranean, a set of genes involved in phosphonate utilization was detected. While another, with a more cosmopolitan distribution, was unique in having an aerobic purine degradation pathway. Together, these results provide a glimpse of the enormous genomic diversity within this clade at a finer resolution than the currently defined clades.

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.

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

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

Wide Distribution of Phage That Infect Freshwater SAR11 Bacteria

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

Wide Distribution of Phage That Infect Freshwater SAR11 Bacteria

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

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

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

Single-amino acid variants reveal evolutionary processes that shape the biogeography of a global SAR11 subclade

Members of the SAR11 order Pelagibacterales dominate the surface oceans. Their extensive diversity challenges emerging operational boundaries defined for microbial 'species' and complicates efforts of population genetics to study their evolution. Here, we employed single-amino acid variants (SAAVs) to investigate ecological and evolutionary forces that maintain the genomic heterogeneity within ubiquitous SAR11 populations we accessed through metagenomic read recruitment using a single isolate genome. Integrating amino acid and protein biochemistry with metagenomics revealed that systematic purifying selection against deleterious variants governs non-synonymous variation among very closely related populations of SAR11. SAAVs partitioned metagenomes into two main groups matching large-scale oceanic current temperatures, and six finer proteotypes that connect distant oceanic regions. These findings suggest that environmentally-mediated selection plays a critical role in the journey of cosmopolitan surface ocean microbial populations, and the idea 'everything is everywhere but the environment selects' has credence even at the finest resolutions.

Marine-freshwater prokaryotic transitions require extensive changes in the predicted proteome

BACKGROUND

The adaptation of a marine prokaryote to live in freshwater environments or vice versa is generally believed to be an unusual and evolutionary demanding process. However, the reasons are not obvious given the similarity of both kinds of habitats.

RESULTS

We have found major differences at the level of the predicted metaproteomes of marine and freshwater habitats with more acidic values of the isoelectric points (pI) in marine microbes. Furthermore, by comparing genomes of marine-freshwater phylogenetic relatives, we have found higher pI values (basic shift) in the freshwater ones. This difference was sharper in secreted > cytoplasmic > membrane proteins. The changes are concentrated on the surface of soluble proteins. It is also detectable at the level of total amino acid composition and involves similarly core and flexible genome- encoded proteins.

CONCLUSIONS

The marked changes at the level of protein amino acid composition and pI provide a tool to predict the preferred habitat of a culture or a metagenome-assembled genome (MAG). The exact physiological explanation for such variations in the pIs and electrostatic surface potentials is not known yet. However, these changes might reflect differences in membrane bioenergetics derived from the absence of significant Na⁺ concentrations in most freshwater habitats. In any case, the changes in amino acid composition in most proteins imply that a long evolutionary time is required to adapt from one type of habitat to the other.

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.

Repeated evolutionary transitions of flavobacteria from marine to non-marine habitats

The taxonomy of marine and non-marine organisms rarely overlap, but the mechanisms underlying this distinction are often unknown. Here, we predicted three major ocean-to-land transitions in the evolutionary history of Flavobacteriaceae, a family known for polysaccharide and peptide degradation. These unidirectional transitions were associated with repeated losses of marine signature genes and repeated gains of non-marine adaptive genes. This included various Na⁺ -dependent transporters, osmolyte transporters and glycoside hydrolases (GH) for sulfated polysaccharide utilization in marine descendants, and in non-marine descendants genes for utilizing the land plant material pectin and genes facilitating terrestrial host interactions. The K⁺ scavenging ATPase was repeatedly gained whereas the corresponding low-affinity transporter repeatedly lost upon transitions, reflecting K⁺ ions are less available to non-marine bacteria. Strikingly, the central metabolism Na⁺ -translocating NADH: quinone dehydrogenase gene was repeatedly gained in marine descendants, whereas the H⁺ -translocating counterpart was repeatedly gained in non-marine lineages. Furthermore, GH genes were depleted in isolates colonizing animal hosts but abundant in bacteria inhabiting other non-marine niches; thus relative abundances of GH versus peptidase genes among Flavobacteriaceae lineages were inconsistent with the marine versus non-marine dichotomy. We suggest that phylogenomic analyses can cast novel light on mechanisms explaining the distribution and ecology of key microbiome components.

Reevaluating the Salty Divide: Phylogenetic Specificity of Transitions between Marine and Freshwater Systems

Marine and freshwater microbial communities are phylogenetically distinct, and transitions between habitat types are thought to be infrequent. We compared the phylogenetic diversity of marine and freshwater microorganisms and identified specific lineages exhibiting notably high or low similarity between marine and freshwater ecosystems using a meta-analysis of 16S rRNA gene tag-sequencing data sets. As expected, marine and freshwater microbial communities differed in the relative abundance of major phyla and contained habitat-specific lineages. At the same time, and contrary to expectations, many shared taxa were observed in both habitats. Based on several metrics, we found that Gammaproteobacteria, Alphaproteobacteria, Bacteroidetes, and Betaproteobacteria contained the highest number of closely related marine and freshwater sequences, suggesting comparatively recent habitat transitions in these groups. Using the abundant alphaproteobacterial group SAR11 as an example, we found evidence that new lineages, beyond the recognized LD12 clade, are detected in freshwater at low but reproducible abundances; this evidence extends beyond the 16S rRNA locus to core genes throughout the genome. Our results suggest that shared taxa are numerous, but tend to occur sporadically and at low relative abundance in one habitat type, leading to an underestimation of transition frequency between marine and freshwater habitats. Rare taxa with abundances near or below detection, including lineages that appear to have crossed the salty divide relatively recently, may possess adaptations enabling them to exploit opportunities for niche expansion when environments are disturbed or conditions change. IMPORTANCE The distribution of microbial diversity across environments yields insight into processes that create and maintain this diversity as well as potential to infer how communities will respond to future environmental changes. We integrated data sets from dozens of freshwater lake and marine samples to compare diversity across open water habitats differing in salinity. Our novel combination of sequence-based approaches revealed lineages that likely experienced a recent transition across habitat types. These taxa are promising targets for studying physiological constraints on salinity tolerance. Our findings contribute to understanding the ecological and evolutionary controls on microbial distributions, and open up new questions regarding the plasticity and adaptability of particular lineages.

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

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

PAR-CLIP for Discovering Target Sites of RNA-Binding Proteins

RNA-binding proteins (RBPs) establish posttranscriptional gene regulation (PTGR) by coordinating the maturation, editing, transport, stability, and translation of cellular RNAs. A variety of experimental approaches have been developed to characterize the RNAs associated with RBPs in vitro as well as in vivo. Our laboratory developed Photoactivatable-Ribonucleoside-Enhanced Cross-Linking and Immunoprecipitation (PAR-CLIP), which in combination with next-generation sequencing enables the identification of RNA targets of RBPs at a nucleotide-level resolution. Here we present an updated and condensed step-by-step PAR-CLIP protocol followed by the description of our RNA-seq data analysis pipeline.

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.

Contrasting patterns of genome-level diversity across distinct co-occurring bacterial populations

To understand the forces driving differentiation and diversification in wild bacterial populations, we must be able to delineate and track ecologically relevant units through space and time. Mapping metagenomic sequences to reference genomes derived from the same environment can reveal genetic heterogeneity within populations, and in some cases, be used to identify boundaries between genetically similar, but ecologically distinct, populations. Here we examine population-level heterogeneity within abundant and ubiquitous freshwater bacterial groups such as the acI Actinobacteria and LD12 Alphaproteobacteria (the freshwater sister clade to the marine SAR11) using 33 single-cell genomes and a 5-year metagenomic time series. The single-cell genomes grouped into 15 monophyletic clusters (termed “tribes”) that share at least 97.9% 16S rRNA identity. Distinct populations were identified within most tribes based on the patterns of metagenomic read recruitments to single-cell genomes representing these tribes. Genetically distinct populations within tribes of the acI Actinobacterial lineage living in the same lake had different seasonal abundance patterns, suggesting these populations were also ecologically distinct. In contrast, sympatric LD12 populations were less genetically differentiated. This suggests that within one lake, some freshwater lineages harbor genetically discrete (but still closely related) and ecologically distinct populations, while other lineages are composed of less differentiated populations with overlapping niches. Our results point at an interplay of evolutionary and ecological forces acting on these communities that can be observed in real time.

Analysis of single-cell genome sequences of bacteria and archaea

Single-cell genome sequencing of individual archaeal and bacterial cells is a vital approach to decipher the genetic makeup of uncultured microorganisms. With this review, we describe single-cell genome analysis with a focus on the unique properties of single-cell sequence data and with emphasis on quality assessment and assurance.

Reconstruction of Diverse Verrucomicrobial Genomes from Metagenome Datasets of Freshwater Reservoirs

The phylum Verrucomicrobia contains freshwater representatives which remain poorly studied at the genomic, taxonomic, and ecological levels. In this work we present eighteen new reconstructed verrucomicrobial genomes from two freshwater reservoirs located close to each other (Tous and Amadorio, Spain). These metagenome-assembled genomes (MAGs) display a remarkable taxonomic diversity inside the phylum and comprise wide ranges of estimated genome sizes (from 1.8 to 6 Mb). Among all Verrucomicrobia studied we found some of the smallest genomes of the Spartobacteria and Opitutae classes described so far. Some of the Opitutae family MAGs were small, cosmopolitan, with a general heterotrophic metabolism with preference for carbohydrates, and capable of xylan, chitin, or cellulose degradation. Besides, we assembled large copiotroph genomes, which contain a higher number of transporters, polysaccharide degrading pathways and in general more strategies for the uptake of nutrients and carbohydrate-based metabolic pathways in comparison with the representatives with the smaller genomes. The diverse genomes revealed interesting features like green-light absorbing rhodopsins and a complete set of genes involved in nitrogen fixation. The large diversity in genome sizes and physiological properties emphasize the diversity of this clade in freshwaters enlarging even further the already broad eco-physiological range of these microbes.

Genomic variation in microbial populations inhabiting the marine subseafloor at deep-sea hydrothermal vents

Little is known about evolutionary drivers of microbial populations in the warm subseafloor of deep-sea hydrothermal vents. Here we reconstruct 73 metagenome-assembled genomes (MAGs) from two geochemically distinct vent fields in the Mid-Cayman Rise to investigate patterns of genomic variation within subseafloor populations. Low-abundance populations with high intra-population diversity coexist alongside high-abundance populations with low genomic diversity, with taxonomic differences in patterns of genomic variation between the mafic Piccard and ultramafic Von Damm vent fields. Populations from Piccard are significantly enriched in nonsynonymous mutations, suggesting stronger purifying selection in Von Damm relative to Piccard. Comparison of nine Sulfurovum MAGs reveals two high-coverage, low-diversity MAGs from Piccard enriched in unique genes related to the cellular membrane, suggesting these populations were subject to distinct evolutionary pressures that may correlate with genes related to nutrient uptake, biofilm formation, or viral invasion. These results are consistent with distinct evolutionary histories between geochemically different vent fields, with implications for understanding evolutionary processes in subseafloor microbial populations.

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.

Excess of non-conservative amino acid changes in marine bacterioplankton lineages with reduced genomes

Surface ocean waters are dominated by planktonic bacterial lineages with highly reduced genomes. The best examples are the cyanobacterial genus Prochlorococcus, the alphaproteobacterial clade SAR11 and the gammaproteobacterial clade SAR86, which together represent over 50% of the cells in surface oceans. Several studies have identified signatures of selection on these lineages in today's ocean and have postulated selection as the primary force throughout their evolutionary history. However, massive loss of genomic DNA in these lineages often occurred in the distant past, and the selective pressures underlying these ancient events have not been assessed. Here, we probe ancient selective pressures by computing %GC-corrected rates of conservative and radical nonsynonymous nucleotide substitutions. Surprisingly, we found an excess of radical changes in several of these lineages in comparison to their relatives with larger genomes. Furthermore, analyses of allelic genome sequences of several populations within these lineages consistently supported that radical replacements are more likely to be deleterious than conservative changes. Our results suggest coincidence of massive genomic DNA losses and increased power of genetic drift, but we also suggest that additional evidence independent of the nucleotide substitution analyses is needed to support a primary role of genetic drift driving ancient genome reduction of marine bacterioplankton lineages.

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

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

Competition and niche separation of pelagic bacteria in freshwater habitats

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

Evolution and adaptation of SAR11 and Cyanobium in a saline Tibetan lake

Lake Qinghai is a unique lacustrine ecosystem located on the Tibetan Plateau and exhibits oligotrophic, alkaline, and saline conditions. Previous studies have focused on the community phylogenetic diversity of bacterioplankton in the ecosystem. This study aimed to address the ecotype diversity of bacterioplankton populations in the unique microbial habitat, using metagenomic sequencing and analysis. Phylogenetic analysis revealed two major bacterial populations: SAR11 IIIa (14% of the total) and Cyanobium (14%). Although the two populations shared high 16S rRNA gene sequence identity (> 98% identity) with their closest marine counterparts, they displayed substantial genomic divergence (≤ 80% average amino acid sequence identity). Comparative genomic analysis identified conservation of carbon and energy storage metabolism (biosynthesis of polyphosphate and polyhydroxyalkanoate) gene operons in the SAR11 IIIa and a cyanate (potential nitrogen source in alkaline conditions) transporter gene operon in the Cyanobium. We further identified genetic signature of positive selection acting on an exodeoxyribonuclease gene of the SAR11 IIIa population, which is potentially associated with DNA repair responsive to strong UV radiation on the high altitude mountain. Taken together, our results revealed the ecosystem-specific gene content of the bacterioplankton populations and provided new insights into their adaptations unique to the Tibetan lake.

Spatial and temporal patterns in the Pelagibacteraceae across an estuarine gradient

Marine bacterial communities show strong spatial and seasonal patterns, often characterized by changes at high taxonomic levels. The Pelagibacteraceae are common members of bacterial communities, with well-documented biogeography at the subclade level. To identify patterns within the subclades, the abundance and diversity of Pelagibacteraceae were analyzed over a two-year period at four stations across an estuarine gradient. Pelagibacteraceae was the most abundant bacterial family, averaging 27% of the community, but varying from 1% to 57% in any one sample. Highest abundances were detected in autumn and winter. Pelagibacteraceae richness was lowest at the most inshore site, and highest in autumn and winter at all sites. Shannon diversity decreased in winter, when a few OTUs dominated the community. Dissolved oxygen, dissolved silicate and prokaryote abundance explained most of the variability in the Pelagibacteraceae communities, with salinity differentiating low salinity communities. The 10 most abundant OTUs included OTUs that varied across sites, with little seasonality as well as those with small site effects, but strong seasonal patterns indicating differences in the niches of individual OTUs. While salinity was important in structuring low salinity communities, higher salinity communities appear to be responding to additional environmental parameters including oxygen, nutrients and other organisms.

Ancestral genome reconstruction identifies the evolutionary basis for trait acquisition in polyphosphate accumulating bacteria

The evolution of complex traits is hypothesized to occur incrementally. Identifying the transitions that lead to extant complex traits may provide a better understanding of the genetic nature of the observed phenotype. A keystone functional group in wastewater treatment processes are polyphosphate accumulating organisms (PAOs), however the evolution of the PAO phenotype has yet to be explicitly investigated and the specific metabolic traits that discriminate non-PAO from PAO are currently unknown. Here we perform the first comprehensive investigation on the evolution of the PAO phenotype using the model uncultured organism Candidatus Accumulibacter phosphatis (Accumulibacter) through ancestral genome reconstruction, identification of horizontal gene transfer, and a kinetic/stoichiometric characterization of Accumulibacter Clade IIA. The analysis of Accumulibacter's last common ancestor identified 135 laterally derived genes, including genes involved in glycogen, polyhydroxyalkanoate, pyruvate and NADH/NADPH metabolisms, as well as inorganic ion transport and regulatory mechanisms. In contrast, pathways such as the TCA cycle and polyphosphate metabolism displayed minimal horizontal gene transfer. We show that the transition from non-PAO to PAO coincided with horizontal gene transfer within Accumulibacter's core metabolism; likely alleviating key kinetic and stoichiometric bottlenecks, such as anaerobically linking glycogen degradation to polyhydroxyalkanoate synthesis. These results demonstrate the utility of investigating the derived genome of a lineage to identify key transitions leading to an extant complex phenotype.

Tuning fresh: radiation through rewiring of central metabolism in streamlined bacteria

Most free-living planktonic cells are streamlined and in spite of their limitations in functional flexibility, their vast populations have radiated into a wide range of aquatic habitats. Here we compared the metabolic potential of subgroups in the Alphaproteobacteria lineage SAR11 adapted to marine and freshwater habitats. Our results suggest that the successful leap from marine to freshwaters in SAR11 was accompanied by a loss of several carbon degradation pathways and a rewiring of the central metabolism. Examples for these are C1 and methylated compounds degradation pathways, the Entner–Doudouroff pathway, the glyoxylate shunt and anapleuretic carbon fixation being absent from the freshwater genomes. Evolutionary reconstructions further suggest that the metabolic modules making up these important freshwater metabolic traits were already present in the gene pool of ancestral marine SAR11 populations. The loss of the glyoxylate shunt had already occurred in the common ancestor of the freshwater subgroup and its closest marine relatives, suggesting that the adaptation to freshwater was a gradual process. Furthermore, our results indicate rapid evolution of TRAP transporters in the freshwater clade involved in the uptake of low molecular weight carboxylic acids. We propose that such gradual tuning of metabolic pathways and transporters toward locally available organic substrates is linked to the formation of subgroups within the SAR11 clade and that this process was critical for the freshwater clade to find and fix an adaptive phenotype.

Genome Reconstruction from Metagenomic Data Sets Reveals Novel Microbes in the Brackish Waters of the Caspian Sea

We present here the findings from a study of the microbiome of the southern basin of the Caspian Sea, the largest water body on Earth disconnected from any ocean and a brackish inland sea. By high-throughput metagenomics, we were able to reconstruct the genomes of representative microbes. The gross community structure (at the phylum level) was different from the structure of typical marine and freshwater communities in temperate open oceans, with the Caspian Sea having freshwater-like amounts of Actinobacteria and Alphaproteobacteria, while Gammaproteobacteria and Betaproteobacteria were present at intermediate levels. We assembled the genomes of several groups and provide detailed descriptions of partial genomes from Actinobacteria, Thaumarchaea, and Alphaproteobacteria. Most belonged to hitherto unknown groups, although they were related to either marine or freshwater groups. The phylogenetic placement of the Caspian genomes indicates that the organisms have multiple and separate phylogenetic origins and that they are related to organisms with both freshwater and marine lineages. Comparative recruitment from global aquatic metagenomes indicated that most Caspian microbes are endemic. However, some Caspian genomes were recruited significantly from either marine water (a member of the Alphaproteobacteria) or freshwater (a member of the Actinobacteria). Reciprocally, some genomes of other origins, such as the marine thaumarchaeon " Candidatus Nitrosopelagicus" or the actinobacterium "Candidatus Actinomarina," were recruited from the Caspian Sea, indicating some degree of overlap with the microbiota of other water bodies. Some of these microbes seem to have a remarkably widespread geographic and environmental distribution.

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

The Chloroflexi CL500-11 clade contributes a large proportion of the bacterial biomass in the oxygenated hypolimnia of deep lakes worldwide, including the world's largest freshwater system, the Laurentian Great Lakes. Traits that allow CL500-11 to thrive and its biogeochemical role in these environments are currently unknown. Here, we found that a CL500-11 population was present mostly in offshore waters along a transect in ultraoligotrophic Lake Michigan (a Laurentian Great Lake). It occurred throughout the water column in spring and only in the hypolimnion during summer stratification, contributing up to 18.1% of all cells. Genome reconstruction from metagenomic data suggested an aerobic, motile, heterotrophic lifestyle, with additional energy being gained through carboxidovory and methylovory. Comparisons to other available streamlined freshwater genomes revealed that the CL500-11 genome contained a disproportionate number of cell wall/capsule biosynthesis genes and the most diverse spectrum of genes involved in the uptake of dissolved organic matter (DOM) substrates, particularly peptides. In situ expression patterns indicated the importance of DOM uptake and protein/peptide turnover, as well as type I and type II carbon monoxide dehydrogenase and flagellar motility. Its location in the water column influenced its gene expression patterns the most. We observed increased bacteriorhodopsin gene expression and a response to oxidative stress in surface waters compared to its response in deep waters. While CL500-11 carries multiple adaptations to an oligotrophic lifestyle, its investment in motility, its large cell size, and its distribution in both oligotrophic and mesotrophic lakes indicate its ability to thrive under conditions where resources are more plentiful. Our data indicate that CL500-11 plays an important role in nitrogen-rich DOM mineralization in the extensive deep-lake hypolimnion habitat.

Specificity and nonspecificity in RNA-protein interactions

To fully understand the regulation of gene expression, it is critical to quantitatively define whether and how RNA-binding proteins (RBPs) discriminate between alternative binding sites in RNAs. Here, we describe new methods that measure protein binding to large numbers of RNA variants, and ways to analyse and interpret data obtained by these approaches, including affinity distributions and free energy landscapes. We discuss how the new methodologies and the associated concepts enable the development of inclusive, quantitative models for RNA-protein interactions that transcend the traditional binary classification of RBPs as either specific or nonspecific.

How do divergent ecological strategies emerge among marine bacterioplankton lineages?

Heterotrophic bacteria in pelagic marine environments are frequently categorized into two canonical ecological groups: patch-associated and free-living. This framework provides a conceptual basis for understanding bacterial utilization of oceanic organic matter. Some patch-associated bacteria are ecologically linked with eukaryotic phytoplankton, and this observation fits with predicted coincidence of their genome expansion with marine phytoplankton diversification. By contrast, free-living bacteria in today's oceans typically live singly with streamlined metabolic and regulatory functions that allow them to grow in nutrient-poor seawater. Recent analyses of marine Alphaproteobacteria suggest that some free-living bacterioplankton lineages evolved from patch-associated ancestors up to several hundred million years ago. While evolutionary analyses agree with the hypothesis that natural selection has maintained these distinct ecological strategies and genomic traits in present-day populations, they do not rule out a major role for genetic drift in driving ancient ecological switches. These two evolutionary forces may have acted on ocean bacteria at different geological time scales and under different geochemical constraints, with possible implications for future adaptations to a changing ocean. New evolutionary models and genomic data are leading to a more comprehensive understanding of marine bacterioplankton evolutionary history.

Single-cell genomics of a rare environmental alphaproteobacterium provides unique insights into Rickettsiaceae evolution

The bacterial family Rickettsiaceae includes a group of well-known etiological agents of many human and vertebrate diseases, including epidemic typhus-causing pathogen Rickettsia prowazekii. Owing to their medical relevance, rickettsiae have attracted a great deal of attention and their host-pathogen interactions have been thoroughly investigated. All known members display obligate intracellular lifestyles, and the best-studied genera, Rickettsia and Orientia, include species that are hosted by terrestrial arthropods. Their obligate intracellular lifestyle and host adaptation is reflected in the small size of their genomes, a general feature shared with all other families of the Rickettsiales. Yet, despite that the Rickettsiaceae and other Rickettsiales families have been extensively studied for decades, many details of the origin and evolution of their obligate host-association remain elusive. Here we report the discovery and single-cell sequencing of 'Candidatus Arcanobacter lacustris', a rare environmental alphaproteobacterium that was sampled from Damariscotta Lake that represents a deeply rooting sister lineage of the Rickettsiaceae. Intriguingly, phylogenomic and comparative analysis of the partial 'Candidatus Arcanobacter lacustris' genome revealed the presence chemotaxis genes and vertically inherited flagellar genes, a novelty in sequenced Rickettsiaceae, as well as several host-associated features. This finding suggests that the ancestor of the Rickettsiaceae might have had a facultative intracellular lifestyle. Our study underlines the efficacy of single-cell genomics for studying microbial diversity and evolution in general, and for rare microbial cells in particular.