Ecogenomics of the SAR11 clade

Microbial ecology of northern Gulf of Mexico estuarine waters

Estuarine and coastal ecosystems are of high economic and ecological importance, owing to their diverse communities and the disproportionate role they play in carbon cycling, particularly in carbon sequestration. Organisms inhabiting these environments must overcome strong natural fluctuations in salinity, nutrients, and turbidity, as well as numerous climate change-induced disturbances such as land loss, sea level rise, and, in some locations, increasingly severe tropical cyclones that threaten to disrupt future ecosystem health. The northern Gulf of Mexico (nGoM) along the Louisiana coast contains dozens of estuaries, including the Mississippi-Atchafalaya River outflow, which dramatically influence the region due to their vast upstream watershed. Nevertheless, the microbiology of these estuaries and surrounding coastal environments has received little attention. To improve our understanding of microbial ecology in the understudied coastal nGoM, we conducted a 16S rRNA gene amplicon survey at eight sites and multiple time points along the Louisiana coast and one inland swamp spanning freshwater to high brackish salinities, totaling 47 duplicated Sterivex (0.2-2.7 µm) and prefilter (>2.7 µm) samples. We cataloged over 13,000 Amplicon Sequence ariants (ASVs) from common freshwater and marine clades such as SAR11 (Alphaproteobacteria), Synechococcus (Cyanobacteria), and acI and Candidatus Actinomarina (Actinobacteria). We observed correlations with freshwater or marine habitats in many organisms and characterized a group of taxa with specialized distributions across brackish water sites, supporting the hypothesis of an endogenous brackish-water community. Additionally, we observed brackish-water associations for several aquatic clades typically considered marine or freshwater taxa, such as SAR11 subclade II, SAR324, and the acI Actinobacteria. The data presented here expand the geographic coverage of microbial ecology in estuarine communities, help delineate the native and transitory members of these environments, and provide critical aquatic microbiological baseline data for coastal and estuarine sites in the nGoM.IMPORTANCEEstuarine and coastal waters are diverse ecosystems influenced by tidal fluxes, interconnected wetlands, and river outflows, which are of high economic and ecological importance. Microorganisms play a pivotal role in estuaries as "first responders" and ecosystem architects, yet despite their ecological importance, they remain underrepresented in microbial studies compared to open ocean environments. This leads to substantial knowledge gaps that are important for understanding global biogeochemical cycling and making decisions about conservation and management strategies in these environments. Our study makes key contributions to the microbial ecology of estuarine and coastal habitats in the northern Gulf of Mexico. Our microbial community data support the concept of a globally distributed, core brackish microbiome and emphasize previously underrecognized brackish-water taxa. Given the projected worsening of land loss, oil spills, and natural disasters in this region, our results will serve as important baseline data for researchers investigating the microbial communities found across estuaries.

Gradients of bacteria in the oceanic water column reveal finely-resolved vertical distributions

Bacterial communities directly influence ecological processes in the ocean, and depth has a major influence due to the changeover in primary energy sources between the sunlit photic zone and dark ocean. Here, we examine the abundance and diversity of bacteria in Monterey Bay depth profiles collected from the surface to just above the sediments (e.g., 2000 m). Bacterial abundance in these Pacific Ocean samples decreased by >1 order of magnitude, from 1.22 ±0.69 ×106 cells ml-1 in the variable photic zone to 1.44 ± 0.25 ×105 and 6.71 ± 1.23 ×104 cells ml-1 in the mesopelagic and bathypelagic, respectively. V1-V2 16S rRNA gene profiling showed diversity increased sharply between the photic and mesopelagic zones. Weighted Gene Correlation Network Analysis clustered co-occurring bacterial amplicon sequence variants (ASVs) into seven subnetwork modules, of which five strongly correlated with depth-related factors. Within surface-associated modules there was a clear distinction between a 'copiotrophic' module, correlating with chlorophyll and dominated by e.g., Flavobacteriales and Rhodobacteraceae, and an 'oligotrophic' module dominated by diverse Oceanospirillales (such as uncultured JL-ETNP-Y6, SAR86) and Pelagibacterales. Phylogenetic reconstructions of Pelagibacterales and SAR324 using full-length 16S rRNA gene data revealed several additional subclades, expanding known microdiversity within these abundant lineages, including new Pelagibacterales subclades Ia.B, Id, and IIc, which comprised 4-10% of amplicons depending on the subclade and depth zone. SAR324 and Oceanospirillales dominated in the mesopelagic, with SAR324 clade II exhibiting its highest relative abundances (17±4%) in the lower mesopelagic (300-750 m). The two newly-identified SAR324 clades showed highest relative abundances in the photic zone (clade III), while clade IV was extremely low in relative abundance, but present across dark ocean depths. Hierarchical clustering placed microbial communities from 900 m samples with those from the bathypelagic, where Marinimicrobia was distinctively relatively abundant. The patterns resolved herein, through high resolution and statistical replication, establish baselines for marine bacterial abundance and taxonomic distributions across the Monterey Bay water column, against which future change can be assessed.

Decoding populations in the ocean microbiome

Understanding the characteristics and structure of populations is fundamental to comprehending ecosystem processes and evolutionary adaptations. While the study of animal and plant populations has spanned a few centuries, microbial populations have been under scientific scrutiny for a considerably shorter period. In the ocean, analyzing the genetic composition of microbial populations and their adaptations to multiple niches can yield important insights into ecosystem function and the microbiome's response to global change. However, microbial populations have remained elusive to the scientific community due to the challenges associated with isolating microorganisms in the laboratory. Today, advancements in large-scale metagenomics and metatranscriptomics facilitate the investigation of populations from many uncultured microbial species directly from their habitats. The knowledge acquired thus far reveals substantial genetic diversity among various microbial species, showcasing distinct patterns of population differentiation and adaptations, and highlighting the significant role of selection in structuring populations. In the coming years, population genomics is expected to significantly increase our understanding of the architecture and functioning of the ocean microbiome, providing insights into its vulnerability or resilience in the face of ongoing global change. Video Abstract.

Seasonality in land-ocean connectivity and local processes control sediment bacterial community structure and function in a High Arctic tidal flat

Climate change is altering patterns of precipitation, cryosphere thaw, and land-ocean influxes, affecting understudied Arctic estuarine tidal flats. These transitional zones between terrestrial and marine systems are hotspots for biogeochemical cycling, often driven by microbial processes. We investigated surface sediment bacterial community composition and function from May to September along a river-intertidal-subtidal-fjord gradient. We paired metabarcoding of in situ communities with in vitro carbon-source utilization assays. Bacterial communities differed in space and time, alongside varying environmental conditions driven by local seasonal processes and riverine inputs, with salinity emerging as the dominant structuring factor. Terrestrial and riverine taxa were found throughout the system, likely transported with runoff. In vitro assays revealed sediment bacteria utilized a broader range of organic matter substrates when incubated in fresh and brackish water compared to marine water. These results highlight the importance of salinity for ecosystem processes in these dynamic tidal flats, with the highest potential for utilization of terrestrially derived organic matter likely limited to tidal flat areas (and times) where sediments are permeated by freshwater. Our results demonstrate that intertidal flats must be included in future studies on impacts of increased riverine discharge and transport of terrestrial organic matter on coastal carbon cycling in a warming Arctic.

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.

Disentangling microbial networks across pelagic zones in the tropical and subtropical global ocean

Microbial interactions are vital in maintaining ocean ecosystem function, yet their dynamic nature and complexity remain largely unexplored. Here, we use association networks to investigate possible ecological interactions in the marine microbiome among archaea, bacteria, and picoeukaryotes throughout different depths and geographical regions of the tropical and subtropical global ocean. Our findings reveal that potential microbial interactions change with depth and geographical scale, exhibiting highly heterogeneous distributions. A few potential interactions were global, meaning they occurred across regions at the same depth, while 11-36% were regional within specific depths. The bathypelagic zone had the lowest proportion of global associations, and regional associations increased with depth. Moreover, we observed that most surface water associations do not persist in deeper ocean layers despite microbial vertical dispersal. Our work contributes to a deeper understanding of the tropical and subtropical global ocean interactome, which is essential for addressing the challenges posed by global change.

Selective cell lysis pressure on rare and abundant prokaryotic taxa across a shelf-to-slope continuum in the Northern South China Sea

IMPORTANCE

Virus-induced host lysis contributes up to 40% of total prokaryotic mortality and plays crucial roles in shaping microbial composition and diversity in the ocean. Nonetheless, what taxon-specific cell lysis is caused by viruses remains to be studied. The present study, therefore, examined the taxon-specific cell lysis and estimated its contribution to the variations in the rare and abundant microbial taxa. The results demonstrate that taxon-specific mortality differed in surface and bottom of the coastal environment. In addition, active rare taxa are more susceptible to heightened lytic pressure and suggested the importance of viral lysis in regulating the microbial community composition. These results improve our understanding of bottom-up (abiotic environmental variables) and top-down (viral lysis) controls contributing to microbial community assembly in the ocean.

Phosphate-related genomic islands as drivers of environmental adaptation in the streamlined marine alphaproteobacterial HIMB59

IMPORTANCE

These results shed light on the evolutionary strategies of microbes with streamlined genomes to adapt and survive in the oligotrophic conditions that dominate the surface waters of the global ocean. At the individual level, these microbes have been subjected to evolutionary constraints that have led to a more efficient use of nutrients, removing non-essential genes named as "streamlining theory." However, at the population level, they conserve a highly diverse gene pool in flexible genomic islands resulting in polyclonal populations on the same genomic background as an evolutionary response to environmental pressures. Localization of these islands at equivalent positions in the genome facilitates horizontal transfer between clonal lineages. This high level of environmental genomic heterogeneity could explain their cosmopolitan distribution. In the case of the order HIMB59 within the class Alphaproteobacteria, two factors exert evolutionary pressure and determine this intraspecific diversity: phages and the concentration of P in the environment.

Metagenome-assembled genomes reveal greatly expanded taxonomic and functional diversification of the abundant marine Roseobacter RCA cluster

BACKGROUND

The RCA (Roseobacter clade affiliated) cluster belongs to the family Roseobacteracea and represents a major Roseobacter lineage in temperate to polar oceans. Despite its prevalence and abundance, only a few genomes and one described species, Planktomarina temperata, exist. To gain more insights into our limited understanding of this cluster and its taxonomic and functional diversity and biogeography, we screened metagenomic datasets from the global oceans and reconstructed metagenome-assembled genomes (MAG) affiliated to this cluster.

RESULTS

The total of 82 MAGs, plus five genomes of isolates, reveal an unexpected diversity and novel insights into the genomic features, the functional diversity, and greatly refined biogeographic patterns of the RCA cluster. This cluster is subdivided into three genera: Planktomarina, Pseudoplanktomarina, and the most deeply branching Candidatus Paraplanktomarina. Six of the eight Planktomarina species have larger genome sizes (2.44-3.12 Mbp) and higher G + C contents (46.36-53.70%) than the four Pseudoplanktomarina species (2.26-2.72 Mbp, 42.22-43.72 G + C%). Cand. Paraplanktomarina is represented only by one species with a genome size of 2.40 Mbp and a G + C content of 45.85%. Three novel species of the genera Planktomarina and Pseudoplanktomarina are validly described according to the SeqCode nomenclature for prokaryotic genomes. Aerobic anoxygenic photosynthesis (AAP) is encoded in three Planktomarina species. Unexpectedly, proteorhodopsin (PR) is encoded in the other Planktomarina and all Pseudoplanktomarina species, suggesting that this light-driven proton pump is the most important mode of acquiring complementary energy of the RCA cluster. The Pseudoplanktomarina species exhibit differences in functional traits compared to Planktomarina species and adaptations to more resource-limited conditions. An assessment of the global biogeography of the different species greatly expands the range of occurrence and shows that the different species exhibit distinct biogeographic patterns. They partially reflect the genomic features of the species.

CONCLUSIONS

Our detailed MAG-based analyses shed new light on the diversification, environmental adaptation, and global biogeography of a major lineage of pelagic bacteria. The taxonomic delineation and validation by the SeqCode nomenclature of prominent genera and species of the RCA cluster may be a promising way for a refined taxonomic identification of major prokaryotic lineages and sublineages in marine and other prokaryotic communities assessed by metagenomics approaches. Video Abstract.

Biogeographic Analysis Suggests Two Types of Planktonic Prokaryote Communities in the Barents Sea

The Barents Sea is one of the most rapidly changing Arctic regions, with an unprecedented sea ice decline and increase in water temperature and salinity. We have studied the diversity of prokaryotic communities using 16S metabarcoding in the western and northeastern parts of the Barents Sea along the Kola Section and the section from Novaya Zemlya to Franz Joseph Land. The hypothesis-independent clustering method revealed the existence of two distinct types of communities. The most common prokaryotic taxa were shared between two types of communities, but their relative abundance was different. It was found that the geographic location of the sampling sites explained more than 30% of the difference between communities, while no statistically significant correlation between environmental parameters and community composition was found. The representatives of the Psychrobacter, Sulfitobacter and Polaribacter genera were dominant in samples from both types of communities. The first type of community was also dominated by members of Halomonas, Pseudoalteromonas, Planococcaceae and an unclassified representative of the Alteromonadaceae family. The second type of community also had a significant proportion of Nitrincolaceae, SAR92, SAR11 Clade I, NS9, Cryomorphaceae and SUP05 representatives. The origin of these communities can be explained by the influence of environmental factors or by the different origins of water masses. This research highlights the importance of studying biogeographic patterns in the Barents Sea in comparison with those in the North Atlantic and Arctic Ocean prokaryote communities.

Partitioning of the denitrification pathway and other nitrite metabolisms within global oxygen deficient zones

Oxygen deficient zones (ODZs) account for about 30% of total oceanic fixed nitrogen loss via processes including denitrification, a microbially mediated pathway proceeding stepwise from NO3- to N2. This process may be performed entirely by complete denitrifiers capable of all four enzymatic steps, but many organisms possess only partial denitrification pathways, either producing or consuming key intermediates such as the greenhouse gas N2O. Metagenomics and marker gene surveys have revealed a diversity of denitrification genes within ODZs, but whether these genes co-occur within complete or partial denitrifiers and the identities of denitrifying taxa remain open questions. We assemble genomes from metagenomes spanning the ETNP and Arabian Sea, and map these metagenome-assembled genomes (MAGs) to 56 metagenomes from all three major ODZs to reveal the predominance of partial denitrifiers, particularly single-step denitrifiers. We find niche differentiation among nitrogen-cycling organisms, with communities performing each nitrogen transformation distinct in taxonomic identity and motility traits. Our collection of 962 MAGs presents the largest collection of pelagic ODZ microorganisms and reveals a clearer picture of the nitrogen cycling community within this environment.

In situ cell division and mortality rates of SAR11, SAR86, Bacteroidetes, and Aurantivirga during phytoplankton blooms reveal differences in population controls

Net growth of microbial populations, that is, changes in abundances over time, can be studied using 16S rRNA fluorescence in situ hybridization (FISH). However, this approach does not differentiate between mortality and cell division rates. We used FISH-based image cytometry in combination with dilution culture experiments to study net growth, cell division, and mortality rates of four bacterial taxa over two distinct phytoplankton blooms: the oligotrophs SAR11 and SAR86, and the copiotrophic phylum Bacteroidetes, and its genus Aurantivirga. Cell volumes, ribosome content, and frequency of dividing cells (FDC) co-varied over time. Among the three, FDC was the most suitable predictor to calculate cell division rates for the selected taxa. The FDC-derived cell division rates for SAR86 of up to 0.8/day and Aurantivirga of up to 1.9/day differed, as expected for oligotrophs and copiotrophs. Surprisingly, SAR11 also reached high cell division rates of up to 1.9/day, even before the onset of phytoplankton blooms. For all four taxonomic groups, the abundance-derived net growth (-0.6 to 0.5/day) was about an order of magnitude lower than the cell division rates. Consequently, mortality rates were comparably high to cell division rates, indicating that about 90% of bacterial production is recycled without apparent time lag within 1 day. Our study shows that determining taxon-specific cell division rates complements omics-based tools and provides unprecedented clues on individual bacterial growth strategies including bottom-up and top-down controls. IMPORTANCE The growth of a microbial population is often calculated from their numerical abundance over time. However, this does not take cell division and mortality rates into account, which are important for deriving ecological processes like bottom-up and top-down control. In this study, we determined growth by numerical abundance and calibrated microscopy-based methods to determine the frequency of dividing cells and subsequently calculate taxon-specific cell division rates in situ. The cell division and mortality rates of two oligotrophic (SAR11 and SAR86) and two copiotrophic (Bacteroidetes and Aurantivirga) taxa during two spring phytoplankton blooms showed a tight coupling for all four taxa throughout the blooms without any temporal offset. Unexpectedly, SAR11 showed high cell division rates days before the bloom while cell abundances remained constant, which is indicative of strong top-down control. Microscopy remains the method of choice to understand ecological processes like top-down and bottom-up control on a cellular level.

Evolutionary pathways for deep-sea adaptation in marine planktonic Actinobacteriota

The deep ocean, one of the largest ecosystems on earth, is dominated by microorganisms that are keystones in the regulation of biogeochemical cycles. However, the evolutionary pathways underlying the specific adaptations required (e.g., high pressure and low temperature) by this unique niche remain understudied. Here, we analyzed the first representatives belonging to the order Acidimicrobiales, a group of marine planktonic Actinobacteriota, that specifically inhabits the aphotic zone of the oceanic water column (>200 m). Compared with their epipelagic counterparts, deep-sea representatives showed the same evolution in genome architecture with higher GC content, longer intergenic spaces as well as higher nitrogen (N-ARSC) and lower carbon (C-ARSC) content in encoded amino acid residue side chains consistent with the higher nitrogen concentration and lower carbon concentration in deep waters compared to the photic zone. Metagenomic recruitment showed distribution patterns that allowed the description of different ecogenomic units within the three deep water-associated genera defined by our phylogenomic analyses (UBA3125, S20-B6 and UBA9410). The entire genus UBA3125 was found exclusively associated with oxygen minimum zones linked to the acquisition of genes involved in denitrification. Genomospecies of genus S20-B6 recruited in samples from both mesopelagic (200-1,000 m) and bathypelagic (1000-4,000 m) zones, including polar regions. Diversity in the genus UBA9410 was higher, with genomospecies widely distributed in temperate zones, others in polar regions, and the only genomospecies associated with abyssal zones (>4,000 m). At the functional level, groups beyond the epipelagic zone have a more complex transcriptional regulation including in their genomes a unique WhiB paralog. In addition, they showed higher metabolic potential for organic carbon and carbohydrate degradation as well as the ability to accumulate glycogen as a source of carbon and energy. This could compensate for energy metabolism in the absence of rhodopsins, which is only present in genomes associated with the photic zone. The abundance in deep samples of cytochrome P450 monooxygenases associated with the genomes of this order suggests an important role in remineralization of recalcitrant compounds throughout the water column.

Ecophysiology and genomics of the brackish water adapted SAR11 subclade IIIa

The Order Pelagibacterales (SAR11) is the most abundant group of heterotrophic bacterioplankton in global oceans and comprises multiple subclades with unique spatiotemporal distributions. Subclade IIIa is the primary SAR11 group in brackish waters and shares a common ancestor with the dominant freshwater IIIb (LD12) subclade. Despite its dominance in brackish environments, subclade IIIa lacks systematic genomic or ecological studies. Here, we combine closed genomes from new IIIa isolates, new IIIa MAGS from San Francisco Bay (SFB), and 460 highly complete publicly available SAR11 genomes for the most comprehensive pangenomic study of subclade IIIa to date. Subclade IIIa represents a taxonomic family containing three genera (denoted as subgroups IIIa.1, IIIa.2, and IIIa.3) that had distinct ecological distributions related to salinity. The expansion of taxon selection within subclade IIIa also established previously noted metabolic differentiation in subclade IIIa compared to other SAR11 subclades such as glycine/serine prototrophy, mosaic glyoxylate shunt presence, and polyhydroxyalkanoate synthesis potential. Our analysis further shows metabolic flexibility among subgroups within IIIa. Additionally, we find that subclade IIIa.3 bridges the marine and freshwater clades based on its potential for compatible solute transport, iron utilization, and bicarbonate management potential. Pure culture experimentation validated differential salinity ranges in IIIa.1 and IIIa.3 and provided detailed IIIa cell size and volume data. This study is an important step forward for understanding the genomic, ecological, and physiological differentiation of subclade IIIa and the overall evolutionary history of SAR11.

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

Structure-informed microbial population genetics elucidate selective pressures that shape protein evolution

Comprehensive sampling of natural genetic diversity with metagenomics enables highly resolved insights into the interplay between ecology and evolution. However, resolving adaptive, neutral, or purifying processes of evolution from intrapopulation genomic variation remains a challenge, partly due to the sole reliance on gene sequences to interpret variants. Here, we describe an approach to analyze genetic variation in the context of predicted protein structures and apply it to a marine microbial population within the SAR11 subclade 1a.3.V, which dominates low-latitude surface oceans. Our analyses reveal a tight association between genetic variation and protein structure. In a central gene in nitrogen metabolism, we observe decreased occurrence of nonsynonymous variants from ligand-binding sites as a function of nitrate concentrations, revealing genetic targets of distinct evolutionary pressures maintained by nutrient availability. Our work yields insights into the governing principles of evolution and enables structure-aware investigations of microbial population genetics.

Single-amplified genomes reveal most streamlined free-living marine bacteria

Evolutionary adaptations of prokaryotes to the environment sometimes result in genome reduction. Our knowledge of this phenomenon among free-living bacteria remains scarce. We address the dynamics and limits of genome reduction by examining one of the most abundant bacteria in the ocean, the SAR86 clade. Despite its abundance, comparative genomics has been limited by the absence of pure cultures and the poor representation in metagenome-assembled genomes. We co-assembled multiple previously available single-amplified genomes to obtain the first complete genomes from members of the four families. All families showed a convergent evolutionary trajectory with characteristic features of streamlined genomes, most pronounced in the TMED112 family. This family has a genome size of ca. 1 Mb and only 1 bp as median intergenic distance, exceeding values found in other abundant microbes such as SAR11, OM43 and Prochlorococcus. This genomic simplification led to a reduction in the biosynthesis of essential molecules, DNA repair-related genes, and the ability to sense and respond to environmental factors, which could suggest an evolutionary dependence on other co-occurring microbes for survival (Black Queen hypothesis). Therefore, these reconstructed genomes within the SAR86 clade provide new insights into the limits of genome reduction in free-living marine bacteria.

SAR11 clade microdiversity and activity during the early spring blooms off Kerguelen Island, Southern Ocean

The ecology of the SAR11 clade, the most abundant bacterial group in the ocean, has been intensively studied in temperate and tropical regions, but its distribution remains largely unexplored in the Southern Ocean. Through amplicon sequencing of the 16S rRNA gene, we assessed the contribution of the SAR11 clade to bacterial community composition in the naturally iron fertilized region off Kerguelen Island. We investigated the upper 300 m at seven sites located in early spring phytoplankton blooms and at one high-nutrient low-chlorophyll site. Despite pronounced vertical patterns of the bacterioplankton assemblages, the SAR11 clade had high relative abundances at all depths and sites, averaging 40% (±15%) of the total community relative abundance. Micro-autoradiography combined with CARD-FISH further revealed that the clade had an overall stable contribution (45%-60% in surface waters) to bacterial biomass production (determined by ³ H-leucine incorporation) during different early bloom stages. The spatio-temporal partitioning of some of the SAR11 subclades suggests a niche specificity and periodic selection of different subclades in response to the fluctuating extreme conditions of the Southern Ocean. These observations improve our understanding of the ecology of the SAR11 clade and its implications in biogeochemical cycles in the rapidly changing Southern Ocean.

Influence of short and long term processes on SAR11 communities in open ocean and coastal systems

SAR11 bacteria dominate the surface ocean and are major players in converting fixed carbon back to atmospheric carbon dioxide. The SAR11 clade is comprised of niche-specialized ecotypes that display distinctive spatiotemporal transitions. We analyzed SAR11 ecotype seasonality in two long-term 16S rRNA amplicon time series representing different North Atlantic regimes: the Sargasso Sea (subtropical ocean-gyre; BATS) and the temperate coastal Western English Channel (WEC). Using phylogenetically resolved amplicon sequence variants (ASVs), we evaluated seasonal environmental constraints on SAR11 ecotype periodicity. Despite large differences in temperature and nutrient availability between the two sites, at both SAR11 succession was defined by summer and winter clusters of ASVs. The summer cluster was dominated by ecotype Ia.3 in both sites. Winter clusters were dominated by ecotypes Ib and IIa.A at BATS and Ia.1 and IIa.B at WEC. A 2-year weekly analysis within the WEC time series showed that the response of SAR11 communities to short-term environmental fluctuations was variable. In 2016, community shifts were abrupt and synchronized to environmental shifts. However, in 2015, changes were gradual and decoupled from environmental fluctuations, likely due to increased mixing from strong winds. We demonstrate that interannual weather variability disturb the pace of SAR11 seasonal progression.

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.

Spatial Variability of Bacterial Community Compositions in the Mariana Trench

Hadal microorganisms play an important role in the biogeochemical processes in marine ecosystems and act as a valuable resource for industrial applications. This paper presents the bacterial community analysis of samples taken from the Challenger Deep within the Mariana Trench, which is the deepest site in the ocean. High-throughput 16S rRNA gene amplicon sequencing was used to reveal that the vertically sampled bacterial populations at eight stations varied at the surface to 10 km depth. The surface water samples harbored a distinct bacterial assemblage, while the mesopelagic and bathyal samples manifested different bacterial community composition, which was not consistent with previous studies. Gammaproteobacteria was the most abundant bacteria in the bathyal and hadal water. The hadal bacterial community consisted mostly of Alteromonadales and Oceanospirillales. The former was widely spread in the water column, which might suggest habitat partitioning at the genus and OTU levels, while the latter might represent hadal-enriched hydrocarbon degraders. The present work complements the current knowledge and understanding of the bathyal and hadal bacterial communities of the Mariana Trench.

Composition and Biogeography of Planktonic Pro- and Eukaryotic Communities in the Atlantic Ocean: Primer Choice Matters

Basin-scale biogeographic observations of marine pelagic pro- and eukaryotic communities are necessary to understand forces driving community composition and for providing a baseline to monitor global change. Deep sequencing of rRNA genes provides community composition at high resolution; yet, it is unclear how the choice of primers affects biogeographic patterns. Here, we re-amplified 16S rRNA genes from DNA sampled during R/V Polarstern Cruise ANT28-5 over a latitudinal transect across the Atlantic Ocean from 52°S to 47°N using universal V4-V5 primers and compared the results with those obtained previously with V5-V6 bacteria-specific primers. For validation of our results, we inferred community composition based on 16S rRNA genes of metagenomes from the same stations and single amplified genomes (SAGs) from the Global Ocean Reference Genome (GORG) database. We found that the universal V4-V5 primers retrieved SAR11 clades with similar relative proportions as those found in the GORG database while the V5-V6 primers recovered strongly diverging clade abundances. We confirmed an inverse bell-shaped distance-decay relationship and a latitudinal diversity gradient that did not decline linearly with absolute latitude in the Atlantic Ocean. Patterns were modified by sampling depth, sequencing depth, choice of primers, and abundance filtering. Especially richness patterns were not robust to methodological change. This study offers a detailed picture of the Atlantic Ocean microbiome using a universal set of PCR primers that allow for the conjunction of biogeographical patterns among organisms from different domains of life.

Long-Read Metagenomics Improves the Recovery of Viral Diversity from Complex Natural Marine Samples

The recovery of DNA from viromes is a major obstacle in the use of long-read sequencing to study their genomes. For this reason, the use of cellular metagenomes (>0.2-μm size range) emerges as an interesting complementary tool, since they contain large amounts of naturally amplified viral genomes from prelytic replication. We have applied second-generation (Illumina NextSeq; short reads) and third-generation (PacBio Sequel II; long reads) sequencing to compare the diversity and features of the viral community in a marine sample obtained from offshore waters of the western Mediterranean. We found that a major wedge of the expected marine viral diversity was directly recovered by the raw PacBio circular consensus sequencing (CCS) reads. More than 30,000 sequences were detected only in this data set, with no homologues in the long- and short-read assembly, and ca. 26,000 had no homologues in the large data set of the Global Ocean Virome 2 (GOV2), highlighting the information gap created by the assembly bias. At the level of complete viral genomes, the performance was similar in both approaches. However, the hybrid long- and short-read assembly provided the longest average length of the sequences and improved the host assignment. Although no novel major clades of viruses were found, there was an increase in the intraclade genomic diversity recovered by long reads that produced an enriched assessment of the real diversity and allowed the discovery of novel genes with biotechnological potential (e.g., endolysin genes).

IMPORTANCE We explored the vast genetic diversity of environmental viruses by using a combination of cellular metagenome (as opposed to virome) sequencing using high-fidelity long-read sequences (in this case, PacBio CCS). This approach resulted in the recovery of a representative sample of the viral population, and it performed better (more phage contigs, larger average contig size) than Illumina sequencing applied to the same sample. By this approach, the many biases of assembly are avoided, as the CCS reads recovers (typically around 5 kb) complete genes and even operons, resulting in a better discovery of the viral gene diversity based on viral marker proteins. Thus, biotechnologically promising genes, such as endolysin genes, can be very efficiently searched with this approach. In addition, hybrid assembly produces more complete and longer contigs, which is particularly important for studying little-known viral groups such as the nucleocytoplasmic large DNA viruses (NCLDV).

Compositional Changes and Co-Occurrence Patterns of Planktonic Bacteria and Microeukaryotes in a Subtropical Estuarine Ecosystem, the Pearl River Delta

Planktonic microorganisms in aquatic ecosystems form complex assemblages of highly interactive taxa and play key roles in biogeochemical cycles. However, the microbial interactions within bacterial and microeukaryotic communities, and the mechanisms underpinning the responses of abundant and rare microbial taxa to environmental disturbances in the river estuary remain unknown. Here, 16S and 18S rRNA gene sequencing were used to investigate the compositional changes and the co-occurrence patterns of bacterial and microeukaryotic communities. The results showed that the rare taxa in the bacterial communities were more prevalent than those in the microeukaryotic communities and may influence the resilience and resistance of microorganisms to environmental variations in estuarine ecosystems. The environmental variations had strong effects on the microeukaryotic communities and their assembly mechanisms but not on the bacterial communities in our studied area. However, based on co-occurrence network analyses, the bacterial communities had stronger links and more complex interactions than microeukaryotic communities, suggesting that bacterial networks may help improve the buffering capacities of the estuarine ecosystem against environmental change. The keystone taxa of bacteria mainly belonged to rare subcommunities, which further illustrates that rare taxa may play fundamental roles in network persistence. Overall, these results provide insights into the microbial responses of aquatic ecosystems to environmental heterogeneity.

Diverse DNA modification in marine prokaryotic and viral communities

DNA chemical modifications, including methylation, are widespread and play important roles in prokaryotes and viruses. However, current knowledge of these modification systems is severely biased towards a limited number of culturable prokaryotes, despite the fact that a vast majority of microorganisms have not yet been cultured. Here, using single-molecule real-time sequencing, we conducted culture-independent 'metaepigenomic' analyses (an integrated analysis of metagenomics and epigenomics) of marine microbial communities. A total of 233 and 163 metagenomic-assembled genomes (MAGs) were constructed from diverse prokaryotes and viruses, respectively, and 220 modified motifs and 276 DNA methyltransferases (MTases) were identified. Most of the MTase genes were not genetically linked with the endonuclease genes predicted to be involved in defense mechanisms against extracellular DNA. The MTase-motif correspondence found in the MAGs revealed 10 novel pairs, 5 of which showed novel specificities and experimentally confirmed the catalytic specificities of the MTases. We revealed novel alternative specificities in MTases that are highly conserved in Alphaproteobacteria, which may enhance our understanding of the co-evolutionary history of the methylation systems and the genomes. Our findings highlight diverse unexplored DNA modifications that potentially affect the ecology and evolution of prokaryotes and viruses in nature.

Searching Metagenomes for New Rhodopsins

Most microbial groups have not been cultivated yet, and the only way to approach the enormous diversity of rhodopsins that they contain in a sensible timeframe is through the analysis of their genomes. High-throughput sequencing technologies have allowed the release of community genomics (metagenomics) of many habitats in the photic zones of the ocean and lakes. Already the harvest is impressive and included from the first bacterial rhodopsin (proteorhodopsin) to the recent discovery of heliorhodopsin by functional metagenomics. However, the search continues using bioinformatic or biochemical routes.

Spatial and temporal dynamics of SAR11 marine bacteria across a nearshore to offshore transect in the tropical Pacific Ocean

Surveys of microbial communities across transitions coupled with contextual measures of the environment provide a useful approach to dissect the factors determining distributions of microorganisms across ecological niches. Here, monthly time-series samples of surface seawater along a transect spanning the nearshore coastal environment within Kāneʻohe Bay on the island of Oʻahu, Hawaiʻi, and the adjacent offshore environment were collected to investigate the diversity and abundance of SAR11 marine bacteria (order Pelagibacterales) over a 2-year time period. Using 16S ribosomal RNA gene amplicon sequencing, the spatiotemporal distributions of major SAR11 subclades and exact amplicon sequence variants (ASVs) were evaluated. Seven of eight SAR11 subclades detected in this study showed distinct subclade distributions across the coastal to offshore environments. The SAR11 community was dominated by seven (of 106 total) SAR11 ASVs that made up an average of 77% of total SAR11. These seven ASVs spanned five different SAR11 subclades (Ia, Ib, IIa, IV, and Va), and were recovered from all samples collected from either the coastal environment, the offshore, or both. SAR11 ASVs were more often restricted spatially to coastal or offshore environments (64 of 106 ASVs) than they were shared among coastal, transition, and offshore environments (39 of 106 ASVs). Overall, offshore SAR11 communities contained a higher diversity of SAR11 ASVs than their nearshore counterparts, with the highest diversity within the little-studied subclade IIa. This study reveals ecological differentiation of SAR11 marine bacteria across a short physiochemical gradient, further increasing our understanding of how SAR11 genetic diversity partitions into distinct ecological units.

Phylogenomics of SAR116 Clade Reveals Two Subclades with Different Evolutionary Trajectories and an Important Role in the Ocean Sulfur Cycle

The SAR116 clade within the class Alphaproteobacteria represents one of the most abundant groups of heterotrophic bacteria inhabiting the surface of the ocean. The small number of cultured representatives of SAR116 (only two to date) is a major bottleneck that has prevented an in-depth study at the genomic level to understand the relationship between genome diversity and its role in the marine environment. In this study, we use all publicly available genomes to provide a genomic overview of the phylogeny, metabolism, and biogeography within the SAR116 clade. This increased genomic diversity has led to the discovery of two subclades that, despite coexisting in the same environment, display different properties in their genomic makeup. One represents a novel subclade for which no pure cultures have been isolated and is composed mainly of single-amplified genomes (SAGs). Genomes within this subclade showed convergent evolutionary trajectories with more streamlined features, such as low GC content (ca. 30%), short intergenic spacers (<22 bp), and strong purifying selection (low ratio of nonsynonymous to synonymous polymorphisms [dN/dS]). Besides, they were more abundant in metagenomic databases recruiting at the deep chlorophyll maximum. Less abundant and restricted to the upper photic layers of the global ocean, the other subclade of SAR116, enriched in metagenome-assembled genomes (MAGs), included the only two pure cultures. Genomic analysis suggested that both clades have a significant role in the sulfur cycle with differences in the way both clades can metabolize dimethylsulfoniopropionate (DMSP). IMPORTANCE The SAR116 clade of Alphaproteobacteria is a ubiquitous group of heterotrophic bacteria inhabiting the surface of the ocean, but the information about their ecology and population genomic diversity is scarce due to the difficulty of getting pure culture isolates. The combination of single-cell genomics and metagenomics has become an alternative approach to study these kinds of microbes. Our results expand the understanding of the genomic diversity, distribution, and lifestyles within this clade and provide evidence of different evolutionary trajectories in the genomic makeup of the two subclades that could serve to illustrate how evolutionary pressure can drive different adaptations to the same environment. Therefore, the SAR116 clade represents an ideal model organism for the study of the evolutionary streamlining of genomes in microbes that have relatively close relatedness to each other.

Prokaryotic responses to a warm temperature anomaly in northeast subarctic Pacific waters

Recent studies on marine heat waves describe water temperature anomalies causing changes in food web structure, bloom dynamics, biodiversity loss, and increased plant and animal mortality. However, little information is available on how water temperature anomalies impact prokaryotes (bacteria and archaea) inhabiting ocean waters. This is a nontrivial omission given their integral roles in driving major biogeochemical fluxes that influence ocean productivity and the climate system. Here we present a time-resolved study on the impact of a large-scale warm water surface anomaly in the northeast subarctic Pacific Ocean, colloquially known as the Blob, on prokaryotic community compositions. Multivariate statistical analyses identified significant depth- and season-dependent trends that were accentuated during the Blob. Moreover, network and indicator analyses identified shifts in specific prokaryotic assemblages from typically particle-associated before the Blob to taxa considered free-living and chemoautotrophic during the Blob, with potential implications for primary production and organic carbon conversion and export.

Enhanced Recovery of Microbial Genes and Genomes From a Marine Water Column Using Long-Read Metagenomics

Third-generation sequencing has penetrated little in metagenomics due to the high error rate and dependence for assembly on short-read designed bioinformatics. However, second-generation sequencing metagenomics (mostly Illumina) suffers from limitations, particularly in the assembly of microbes with high microdiversity and retrieval of the flexible (adaptive) fraction of prokaryotic genomes. Here, we have used a third-generation technique to study the metagenome of a well-known marine sample from the mixed epipelagic water column of the winter Mediterranean. We have compared PacBio Sequel II with the classical approach using Illumina Nextseq short reads followed by assembly to study the metagenome. Long reads allow for efficient direct retrieval of complete genes avoiding the bias of the assembly step. Besides, the application of long reads on metagenomic assembly allows for the reconstruction of much more complete metagenome-assembled genomes (MAGs), particularly from microbes with high microdiversity such as Pelagibacterales. The flexible genome of reconstructed MAGs was much more complete containing many adaptive genes (some with biotechnological potential). PacBio Sequel II CCS appears particularly suitable for cellular metagenomics due to its low error rate. For most applications of metagenomics, from community structure analysis to ecosystem functioning, long reads should be applied whenever possible. Specifically, for in silico screening of biotechnologically useful genes, or population genomics, long-read metagenomics appears presently as a very fruitful approach and can be analyzed from raw reads before a computationally demanding (and potentially artifactual) assembly step.

Metapangenomics reveals depth-dependent shifts in metabolic potential for the ubiquitous marine bacterial SAR324 lineage

BACKGROUND

Oceanic microbiomes play a pivotal role in the global carbon cycle and are central to the transformation and recycling of carbon and energy in the ocean's interior. SAR324 is a ubiquitous but poorly understood uncultivated clade of Deltaproteobacteria that inhabits the entire water column, from ocean surface waters to its deep interior. Although some progress has been made in elucidating potential metabolic traits of SAR324 in the dark ocean, very little is known about the ecology and the metabolic capabilities of this group in the euphotic and twilight zones. To investigate the comparative genomics, ecology, and physiological potential of the SAR324 clade, we examined the distribution and variability of key genomic features and metabolic pathways in this group from surface waters to the abyss in the North Pacific Subtropical Gyre, one of the largest biomes on Earth.

RESULTS

We leveraged a pangenomic ecological approach, combining spatio-temporally resolved single-amplified genome, metagenomic, and metatranscriptomic datasets. The data revealed substantial genomic diversity throughout the SAR324 clade, with distinct depth and temporal distributions that clearly differentiated ecotypes. Phylogenomic subclade delineation, environmental distributions, genomic feature similarities, and metabolic capacities revealed strong congruence. The four SAR324 ecotypes delineated in this study revealed striking divergence from one another with respect to their habitat-specific metabolic potentials. The ecotypes living in the dark or twilight oceans shared genomic features and metabolic capabilities consistent with a sulfur-based chemolithoautotrophic lifestyle. In contrast, those inhabiting the sunlit ocean displayed higher plasticity energy-related metabolic pathways, supporting a presumptive photoheterotrophic lifestyle. In epipelagic SAR324 ecotypes, we observed the presence of two types of proton-pumping rhodopsins, as well as genomic, transcriptomic, and ecological evidence for active photoheterotrophy, based on xanthorhodopsin-like light-harvesting proteins.

CONCLUSIONS

Combining pangenomic and both metagenomic and metatranscriptomic profiling revealed a striking divergence in the vertical distribution, genomic composition, metabolic potential, and predicted lifestyle strategies of geographically co-located members of the SAR324 bacterial clade. The results highlight the utility of metapangenomic approaches employed across environmental gradients, to decipher the properties and variation in function and ecological traits of specific phylogenetic clades within complex microbiomes. Video abstract.

Adaptability of a caproate-producing bacterium contributes to its dominance in an anaerobic fermentation system

The transformation of diverse feedstocks into medium-chain fatty acids (MCFAs) by mixed cultures is a promising biorefinery route because of the high value of MCFAs. A particular concern is how to maintain the microbial consortia in mixed cultures to achieve stable MCFA production. Chinese strong aroma-type liquor (Baijiu) fermentation system continually produces caproic acid for decades through a spontaneous inoculation of anaerobes from pit mud into fermented grains. Therefore, illuminating the dominant caproate-producing bacterium (CPB) in pit mud and how the CPB sustains in the spontaneous fermentation system will benefit to reveal the microbiological mechanisms of the stable caproate production. Here, we examined pit mud samples across four Chinese strong aroma-type Baijiu producing areas and found that a caproate-producing Caproicibacterium sp. was widely distributed in these distilleries with relative abundance ranging from 1.4% to 35.5% and an average abundance of 11.4%. Through controlling carbon source availability, we achieved different simplified caproate-producing consortia and found that the growth advantage of Caproicibacterium sp. was highly dependent on glucose. Then two strains, named Caproicibacterium sp. LBM19010 and Caproicibacterium sp. JNU-WLY1368, were isolated from pit mud of two regions. The metabolic versatility of this bacterium utilizing starch, maltose, glucose and lactate reflected its adaptability to the fermentation environment where these carbon sources coexist. The simultaneous utilization of glucose and lactate contributed to the balance between cell growth and pH homeostasis. This study reveals that multiple adaptation strategies employed by the predominant CPB promotes its stability and dominance in a saccharide- and lactate-rich anaerobic habitat. IMPORTANCE Chinese strong aroma-type liquor (Baijiu) fermentation environment is a typical medium-chain fatty acid producing system with complex nutrients. Although several studies have revealed the correlation between microbial community composition and abiotic factors, the adaptation mechanisms of dominant species to abiotic environment are still unknown in this special anaerobic habitat. This study identified the predominant CPB in Chinese strong aroma-type Baijiu fermentation system. Metabolic versatility and flexibility of the dominant CPB with a small-size genome indicated that this bacterium can effectively exploit available carbon and nitrogen sources, which could be a key factor to promote its ecological success in a multi-species environment. The understanding of growth and metabolic features of CPB responsible for its dominance in microbial community will not only contribute to the improvement of Chinese strong aroma-type Baijiu production but also expand its potential industrial applications in caproate production.

The microbiome of the Black Sea water column analyzed by shotgun and genome centric metagenomics

BACKGROUND

The Black Sea is the largest brackish water body in the world, although it is connected to the Mediterranean Sea and presents an upper water layer similar to some regions of the former, albeit with lower salinity and temperature. Despite its well-known hydrology and physicochemical features, this enormous water mass remains poorly studied at the microbial genomics level.

RESULTS

We have sampled its different water masses and analyzed the microbiome by shotgun and genome-resolved metagenomics, generating a large number of metagenome-assembled genomes (MAGs) from them. We found various similarities with previously described Black Sea metagenomic datasets, that show remarkable stability in its microbiome. Our datasets are also comparable to other marine anoxic water columns like the Cariaco Basin. The oxic zone resembles to standard marine (e.g. Mediterranean) photic zones, with Cyanobacteria (Synechococcus but a conspicuously absent Prochlorococcus), and photoheterotrophs domination (largely again with marine relatives). The chemocline presents very different characteristics from the oxic surface with many examples of chemolithotrophic metabolism (Thioglobus) and facultatively anaerobic microbes. The euxinic anaerobic zone presents, as expected, features in common with the bottom of meromictic lakes with a massive dominance of sulfate reduction as energy-generating metabolism, a few (but detectable) methanogenesis marker genes, and a large number of "dark matter" streamlined genomes of largely unpredictable ecology.

CONCLUSIONS

The Black Sea oxic zone presents many similarities to the global ocean while the redoxcline and euxinic water masses have similarities to other similar aquatic environments of marine (Cariaco Basin or other Black Sea regions) or freshwater (meromictic monimolimnion strata) origin. The MAG collection represents very well the different types of metabolisms expected in this kind of environment. We are adding critical information about this unique and important ecosystem and its microbiome.

Description of Candidatus Mesopelagibacter carboxydoxydans and Candidatus Anoxipelagibacter denitrificans: Nitrate-reducing SAR11 genera that dominate mesopelagic and anoxic marine zones

The diverse and ubiquitous members of the SAR11 lineage (Alphaproteobacteria) represent up to 30-40% of the surface and mesopelagic oceanic microbial communities. However, the molecular and ecological mechanisms that differentiate closely related, yet distinct, SAR11 members that often co-occur under similar environmental conditions remain speculative. Recently, two mesopelagic and oxygen minimum zone (OMZ)-associated subclades of SAR11 (Ic and IIa.A) were described using single-cell amplified genomes (SAGs) linked to nitrate reduction in OMZs. In this current study, the collection of genomes belonging to these two subclades was expanded with thirteen new metagenome-assembled genomes (MAGs), thus providing a more detailed phylogenetic and functional characterization of these subclades. Gene content-based predictions of metabolic functions revealed similarities in central carbon metabolism between subclades Ic and IIa.A and surface SAR11 clades, with small variations in central pathways. These variations included more versatile sulfur assimilation pathways, as well as a previously predicted capacity for nitrate reduction that conferred unique versatility on mesopelagic-adapted clades compared to their surface counterparts. Finally, consistent with previously reported abundances of carbon monoxide (CO) in surface and mesopelagic waters, subclades Ia (surface) and Ic (mesopelagic) have the genetic potential to oxidize carbon monoxide (CO), presumably taking advantage of this abundant compound as an electron donor. Based on genomic analyses, environmental distribution and metabolic reconstruction, we propose two new SAR11 genera, Ca. Mesopelagibacter carboxydoxydans (subclade Ic) and Ca. Anoxipelagibacter denitrificans (subclade IIa.A), which represent members of the mesopelagic and OMZ-adapted SAR11 clades.

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.

Isolation of SAR11 Marine Bacteria from Cryopreserved Seawater

While marine microorganisms are frequently studied in their natural environment, isolated strains are invaluable resources that can be used in controlled experiments to expand upon direct observations from natural systems. Here, we sought a means to enhance culture collections of SAR11 marine bacteria by testing the use of seawater cryopreserved with glycerol as an inoculum. Using a raw seawater sample collected from the tropical Pacific Ocean, a subsample was diluted in seawater growth medium to create 576 2-ml dilution cultures containing 5 cells each and incubated for a high-throughput culturing (HTC) experiment, while another portion was cryopreserved in 10% glycerol. After 10 months, a cryopreserved aliquot was thawed and used to create a second cultivation experiment of 480 2-ml cultures containing 5 cells each and 470 cultures containing 105 cells each. The raw seawater cultivation experiment resulted in the successful isolation of 54 monocultures and 29 mixed cultures, while cryopreserved seawater resulted in 59 monocultures and 29 mixed cultures. Combined, the cultures included 51 SAR11 isolates spanning 11 unique 16S rRNA gene amplicon sequence variants (ASVs) from the raw seawater inoculum and 74 SAR11 isolates spanning 13 unique ASVs from cryopreserved seawater. A vast majority (92%) of SAR11 isolates from the two HTC experiments were members of SAR11 subclade Ia, though subclades IIIa and Va were also recovered from cryopreserved seawater and subclade Ib was recovered from both. The four most abundant SAR11 subclade Ia ASVs found in the initial seawater environmental sample were isolated by both approaches.IMPORTANCE High-throughput dilution culture has proved to be a successful approach to bring some difficult-to-isolate planktonic microorganisms into culture, including the highly abundant SAR11 lineage of marine bacteria. While the long-term preservation of bacterial isolates by freezing them in the presence of cryoprotectants, such as glycerol, has been shown to be an effective method of storing viable cells over long time periods (i.e., years), to our knowledge it had not previously been tested for its efficacy in preserving raw seawater for later use as an inoculum for high-throughput cultivation experiments. We found that SAR11 and other abundant marine bacteria could be isolated from seawater that was previously cryopreserved for nearly 10 months at a rate of culturability similar to that of the same seawater used fresh, immediately after collection. Our findings (i) expand the potential of high-throughput cultivation experiments to include testing when immediate isolation experiments are impractical, (ii) allow for targeted isolation experiments from specific samples based on analyses such as microbial community structure, and (iii) enable cultivation experiments across a wide range of other conditions that would benefit from having source inocula available over extended periods of time.

Genomes of the "Candidatus Actinomarinales" Order: Highly Streamlined Marine Epipelagic Actinobacteria

Microbiology is in a new age in which sequence databases are primary sources of information about many microbes. However, in-depth analysis of environmental genomes thus retrieved is essential to substantiate the new knowledge.

“Candidatus Actinomarinales” was defined as a subclass of exclusively marine Actinobacteria with small cells and genomes. We have collected all the available genomes in databases to assess the diversity included in this group and analyzed it by comparative genomics. We have found the equivalent of five genera and 18 genomospecies. They have genome reduction parameters equal to those of freshwater actinobacterial “Candidatus Nanopelagicales” or marine alphaproteobacterial Pelagibacterales. Genome recruitment shows that they are found only in the photic zone and mainly in surface waters, with only one genus that is found preferentially at or below the deep chlorophyll maximum. “Ca. Actinomarinales” show a highly conserved core genome (80% of the gene families conserved for the whole order) with a saturation of genomic diversity of the flexible genome at the genomospecies level. We found only a flexible genomic island preserved throughout the order; it is related to the sugar decoration of the envelope and uses several tRNAs as hot spots to increase its genomic diversity. Populations had a discrete level of sequence diversity similar to other marine microbes but drastically different from the much higher levels found for Pelagibacterales. Genomic analysis suggests that they are all aerobic photoheterotrophs with one type 1 rhodopsin and a heliorhodopsin. Like other actinobacteria, they possess the F420 coenzyme biosynthesis pathway, and its lower reduction potential could provide access to an increased range of redox chemical transformations. Last, sequence analysis revealed the first “Ca. Actinomarinales” phages, including a prophage, with metaviromic islands related to sialic acid cleavage.

IMPORTANCE Microbiology is in a new age in which sequence databases are primary sources of information about many microbes. However, in-depth analysis of environmental genomes thus retrieved is essential to substantiate the new knowledge. Here, we study 182 genomes belonging to the only known exclusively marine pelagic group of the phylum Actinobacteria. The aquatic branch of this phylum is largely known from environmental sequencing studies (single-amplified genomes [SAGs] and metagenome-assembled genomes [MAGs]), and we have collected and analyzed the available information present in databases about the “Ca. Actinomarinales.” They are among the most streamlined microbes to live in the epipelagic zone of the ocean, and their study is critical to obtain a proper view of the diversity of Actinobacteria and their role in aquatic ecosystems.

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

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

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.

Metagenome Mining Reveals Hidden Genomic Diversity of Pelagimyophages in Aquatic Environments

SAR11 clade members are among the most abundant bacteria on Earth. Their study is complicated by their great diversity and difficulties in being grown and manipulated in the laboratory. On the other hand, and due to their extraordinary abundance, metagenomic data sets provide enormous richness of information about these microbes. Given the major role played by phages in the lifestyle and evolution of prokaryotic cells, the contribution of several new bacteriophage genomes preying on this clade opens windows into the infection strategies and life cycle of its viruses. Such strategies could provide models of attack of large-genome phages preying on streamlined aquatic microbes.

The SAR11 clade is one of the most abundant bacterioplankton groups in surface waters of most of the oceans and lakes. However, only 15 SAR11 phages have been isolated thus far, and only one of them belongs to the Myoviridae family (pelagimyophages). Here, we have analyzed 26 sequences of myophages that putatively infect the SAR11 clade. They have been retrieved by mining ca. 45 Gbp aquatic assembled cellular metagenomes and viromes. Most of the myophages were obtained from the cellular fraction (0.2 μm), indicating a bias against this type of virus in viromes. We have found the first myophages that putatively infect Candidatus Fonsibacter (freshwater SAR11) and another group putatively infecting bathypelagic SAR11 phylogroup Ic. The genomes have similar sizes and maintain overall synteny in spite of low average nucleotide identity values, revealing high similarity to marine cyanomyophages. Pelagimyophages recruited metagenomic reads widely from several locations but always much more from cellular metagenomes than from viromes, opposite to what happens with pelagipodophages. Comparing the genomes resulted in the identification of a hypervariable island that is related to host recognition. Interestingly, some genes in these islands could be related to host cell wall synthesis and coinfection avoidance. A cluster of curli-related proteins was widespread among the genomes, although its function is unclear.

IMPORTANCE SAR11 clade members are among the most abundant bacteria on Earth. Their study is complicated by their great diversity and difficulties in being grown and manipulated in the laboratory. On the other hand, and due to their extraordinary abundance, metagenomic data sets provide enormous richness of information about these microbes. Given the major role played by phages in the lifestyle and evolution of prokaryotic cells, the contribution of several new bacteriophage genomes preying on this clade opens windows into the infection strategies and life cycle of its viruses. Such strategies could provide models of attack of large-genome phages preying on streamlined aquatic microbes.