Genomic insights to SAR86, an abundant and uncultivated marine bacterial lineage

Bacteria in the 16S rRNA clade SAR86 are among the most abundant uncultivated constituents of microbial assemblages in the surface ocean for which little genomic information is currently available. Bioinformatic techniques were used to assemble two nearly complete genomes from marine metagenomes and single-cell sequencing provided two more partial genomes. Recruitment of metagenomic data shows that these SAR86 genomes substantially increase our knowledge of non-photosynthetic bacteria in the surface ocean. Phylogenomic analyses establish SAR86 as a basal and divergent lineage of γ-proteobacteria, and the individual genomes display a temperature-dependent distribution. Modestly sized at 1.25–1.7 Mbp, the SAR86 genomes lack several pathways for amino-acid and vitamin synthesis as well as sulfate reduction, trends commonly observed in other abundant marine microbes. SAR86 appears to be an aerobic chemoheterotroph with the potential for proteorhodopsin-based ATP generation, though the apparent lack of a retinal biosynthesis pathway may require it to scavenge exogenously-derived pigments to utilize proteorhodopsin. The genomes contain an expanded capacity for the degradation of lipids and carbohydrates acquired using a wealth of tonB-dependent outer membrane receptors. Like the abundant planktonic marine bacterial clade SAR11, SAR86 exhibits metabolic streamlining, but also a distinct carbon compound specialization, possibly avoiding competition.     

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.     

The unique metabolism of SAR11 aquatic bacteria

The deeply branching clade of abundant, globally distributed aquatic α-Proteobacteria known as “SAR11”, are adapted to nutrient-poor environments such as the surface waters of the open ocean. Unknown prior to 1990, uncultured until 2002, members of the SAR11 clade can now be cultured in artificial, defined media to densities three orders of magnitude higher than in unamended natural media. Cultivation in natural and defined media has confirmed genomic and metagenomic predictions such as an inability to reduce sulfate to sulfide, a requirement for pyruvate, an ability to oxidize a wide variety of methylated and one-carbon compounds for energy, and an unusual form of conditional glycine auxotrophy. Here we describe the metabolism of the SAR11 type strain Candidatus “Pelagibacter ubique” str. HTCC1062, as revealed by genome-assisted studies of laboratory cultures. We also describe the discovery of SAR11 and field studies that have been done on natural populations.     

Temporal variability and coherence of euphotic zone bacterial communities over a decade in the Southern California Bight

Time-series are critical to understanding long-term natural variability in the oceans. Bacterial communities in the euphotic zone were investigated for over a decade at the San Pedro Ocean Time-series station (SPOT) off southern California. Community composition was assessed by Automated Ribosomal Intergenic Spacer Analysis (ARISA) and coupled with measurements of oceanographic parameters for the surface ocean (0–5 m) and deep chlorophyll maximum (DCM, average depth ∼30 m). SAR11 and cyanobacterial ecotypes comprised typically more than one-third of the measured community; diversity within both was temporally variable, although a few operational taxonomic units (OTUs) were consistently more abundant. Persistent OTUs, mostly Alphaproteobacteria (SAR11 clade), Actinobacteria and Flavobacteria, tended to be abundant, in contrast to many rarer yet intermittent and ephemeral OTUs. Association networks revealed potential niches for key OTUs from SAR11, cyanobacteria, SAR86 and other common clades on the basis of robust correlations. Resilience was evident by the average communities drifting only slightly as years passed. Average Bray-Curtis similarity between any pair of dates was ∼40%, with a slight decrease over the decade and obvious near-surface seasonality; communities 8–10 years apart were slightly more different than those 1–4 years apart with the highest rate of change at 0–5 m between communities <4 years apart. The surface exhibited more pronounced seasonality than the DCM. Inter-depth Bray-Curtis similarities repeatedly decreased as the water column stratified each summer. Environmental factors were better predictors of shifts in community composition than months or elapsed time alone; yet, the best predictor was community composition at the other depth (that is, 0–5 m versus DCM).     

Evolutionary analysis of a streamlined lineage of surface ocean Roseobacters

The vast majority of surface ocean bacteria are uncultivated. Compared with their cultured relatives, they frequently exhibit a streamlined genome, reduced G+C content and distinct gene repertoire. These genomic traits are relevant to environmental adaptation, and have generally been thought to become fixed in marine bacterial populations through selection. Using single-cell genomics, we sequenced four uncultivated cells affiliated with the ecologically relevant Roseobacter clade and used a composition-heterogeneous Bayesian phylogenomic model to resolve these single-cell genomes into a new clade. This lineage has no representatives in culture, yet accounts for ∼35% of Roseobacters in some surface ocean waters. Analyses of multiple genomic traits, including genome size, G+C content and percentage of noncoding DNA, suggest that these single cells are representative of oceanic Roseobacters but divergent from isolates. Population genetic analyses showed that substitution of physicochemically dissimilar amino acids and replacement of G+C-rich to G+C-poor codons are accelerated in the uncultivated clade, processes that are explained equally well by genetic drift as by the more frequently invoked explanation of natural selection. The relative importance of drift vs selection in this clade, and perhaps in other marine bacterial clades with streamlined G+C-poor genomes, remains unresolved until more evidence is accumulated.     

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.     

Global biogeography of SAR11 marine bacteria

The ubiquitous SAR11 bacterial clade is the most abundant type of organism in the world's oceans, but the reasons for its success are not fully elucidated. We analysed 128 surface marine metagenomes, including 37 new Antarctic metagenomes. The large size of the data set enabled internal transcribed spacer (ITS) regions to be obtained from the Southern polar region, enabling the first global characterization of the distribution of SAR11, from waters spanning temperatures ?2 to 30°C. Our data show a stable co‐occurrence of phylotypes within both ‘tropical’ (>20°C) and ‘polar’ (<10°C) biomes, highlighting ecological niche differentiation between major SAR11 subgroups. All phylotypes display transitions in abundance that are strongly correlated with temperature and latitude. By assembling SAR11 genomes from Antarctic metagenome data, we identified specific genes, biases in gene functions and signatures of positive selection in the genomes of the polar SAR11—genomic signatures of adaptive radiation. Our data demonstrate the importance of adaptive radiation in the organism's ability to proliferate throughout the world's oceans, and describe genomic traits characteristic of different phylotypes in specific marine biomes.     

Phaeobacter gallaeciensis genomes from globally opposite locations reveal high similarity of adaptation to surface life

Phaeobacter gallaeciensis, a member of the abundant marine Roseobacter clade, is known to be an effective colonizer of biotic and abiotic marine surfaces. Production of the antibiotic tropodithietic acid (TDA) makes P. gallaeciensis a strong antagonist of many bacteria, including fish and mollusc pathogens. In addition to TDA, several other secondary metabolites are produced, allowing the mutualistic bacterium to also act as an opportunistic pathogen. Here we provide the manually annotated genome sequences of the P. gallaeciensis strains DSM 17395 and 2.10, isolated at the Atlantic coast of north western Spain and near Sydney, Australia, respectively. Despite their isolation sites from the two different hemispheres, the genome comparison demonstrated a surprisingly high level of synteny (only 3% nucleotide dissimilarity and 88% and 93% shared genes). Minor differences in the genomes result from horizontal gene transfer and phage infection. Comparison of the P. gallaeciensis genomes with those of other roseobacters revealed unique genomic traits, including the production of iron-scavenging siderophores. Experiments supported the predicted capacity of both strains to grow on various algal osmolytes. Transposon mutagenesis was used to expand the current knowledge on the TDA biosynthesis pathway in strain DSM 17395. This first comparative genomic analysis of finished genomes of two closely related strains belonging to one species of the Roseobacter clade revealed features that provide competitive advantages and facilitate surface attachment and interaction with eukaryotic hosts.     

Distribution of Roseobacter RCA and SAR11 lineages in the North Sea and characteristics of an abundant RCA isolate

The Roseobacter group and SAR11 clade constitute high proportions of the marine bacterioplankton, but only scarce information exists on the abundance of distinct populations of either lineage. Therefore, we quantified the abundance of the largest cluster of the Roseobacter group, the RCA (Roseobacter clade affiliated) cluster together with the SAR11 clade by quantitative PCR in the southern and eastern North Sea. The RCA cluster constituted up to 15 and 21% of total bacterial 16S ribosomal RNA (rRNA) genes in September 2005 and May 2006, respectively. At a few stations, the RCA cluster exceeded the SAR11 clade, whereas at most stations, SAR11 constituted higher fractions with maxima of 37%. In most samples, only one RCA ribotype was detected. RCA abundance was positively correlated with phaeopigments, chlorophyll, dissolved and particulate organic carbon (POC), turnover rates of dissolved free amino acids (DFAAs), temperature, and negatively correlated with salinity. The SAR11 clade was only correlated with POC (negatively, May) and with DFAA turnover rates (positively, September). An abundant RCA strain, ‘Candidatus Planktomarina temperata'', was isolated from the southern North Sea. This strain has an identical 16S rRNA gene sequence to the dominant RCA ribotype. Detection of the pufM gene, coding for a subunit of the reaction center of bacteriochlorophyll a, indicates the potential of the isolate for aerobic anoxygenic photosynthesis. Our study shows that a distinct population of the RCA cluster constitutes an abundant bacterioplankton group in a neritic sea of the temperate zone and indicates that this population has an important role during decaying phytoplankton blooms.     

A new class of marine Euryarchaeota group II from the mediterranean deep chlorophyll maximum

We have analyzed metagenomic fosmid clones from the deep chlorophyll maximum (DCM), which, by genomic parameters, correspond to the 16S ribosomal RNA (rRNA)-defined marine Euryarchaeota group IIB (MGIIB). The fosmid collections associated with this group add up to 4 Mb and correspond to at least two species within this group. From the proposed essential genes contained in the collections, we infer that large sections of the conserved regions of the genomes of these microbes have been recovered. The genomes indicate a photoheterotrophic lifestyle, similar to that of the available genome of MGIIA (assembled from an estuarine metagenome in Puget Sound, Washington Pacific coast), with a proton-pumping rhodopsin of the same kind. Several genomic features support an aerobic metabolism with diversified substrate degradation capabilities that include xenobiotics and agar. On the other hand, these MGIIB representatives are non-motile and possess similar genome size to the MGIIA-assembled genome, but with a lower GC content. The large phylogenomic gap with other known archaea indicates that this is a new class of marine Euryarchaeota for which we suggest the name Thalassoarchaea. The analysis of recruitment from available metagenomes indicates that the representatives of group IIB described here are largely found at the DCM (ca. 50 m deep), in which they are abundant (up to 0.5% of the reads), and at the surface mostly during the winter mixing, which explains formerly described 16S rRNA distribution patterns. Their uneven representation in environmental samples that are close in space and time might indicate sporadic blooms.     

Discovery of a SAR11 growth requirement for thiamin's pyrimidine precursor and its distribution in the Sargasso Sea

Vitamin traffic, the production of organic growth factors by some microbial community members and their use by other taxa, is being scrutinized as a potential explanation for the variation and highly connected behavior observed in ocean plankton by community network analysis. Thiamin (vitamin B₁), a cofactor in many essential biochemical reactions that modify carbon–carbon bonds of organic compounds, is distributed in complex patterns at subpicomolar concentrations in the marine surface layer (0–300 m). Sequenced genomes from organisms belonging to the abundant and ubiquitous SAR11 clade of marine chemoheterotrophic bacteria contain genes coding for a complete thiamin biosynthetic pathway, except for thiC, encoding the 4-amino-5-hydroxymethyl-2-methylpyrimidine (HMP) synthase, which is required for de novo synthesis of thiamin''s pyrimidine moiety. Here we demonstrate that the SAR11 isolate ‘Candidatus Pelagibacter ubique'', strain HTCC1062, is auxotrophic for the thiamin precursor HMP, and cannot use exogenous thiamin for growth. In culture, strain HTCC1062 required 0.7 zeptomoles per cell (ca. 400 HMP molecules per cell). Measurements of dissolved HMP in the Sargasso Sea surface layer showed that HMP ranged from undetectable (detection limit: 2.4 pM) to 35.7 pM, with maximum concentrations coincident with the deep chlorophyll maximum. In culture, some marine cyanobacteria, microalgae and bacteria exuded HMP, and in the Western Sargasso Sea, HMP profiles changed between the morning and evening, suggesting a dynamic biological flux from producers to consumers.     

Abundance and Distribution of Dimethylsulfoniopropionate Degradation Genes and the Corresponding Bacterial Community Structure at Dimethyl Sulfide Hot Spots in the Tropical and Subtropical Pacific Ocean

Dimethylsulfoniopropionate (DMSP) is mainly produced by marine phytoplankton but is released into the microbial food web and degraded by marine bacteria to dimethyl sulfide (DMS) and other products. To reveal the abundance and distribution of bacterial DMSP degradation genes and the corresponding bacterial communities in relation to DMS and DMSP concentrations in seawater, we collected surface seawater samples from DMS hot spot sites during a cruise across the Pacific Ocean. We analyzed the genes encoding DMSP lyase (dddP) and DMSP demethylase (dmdA), which are responsible for the transformation of DMSP to DMS and DMSP assimilation, respectively. The averaged abundance (±standard deviation) of these DMSP degradation genes relative to that of the 16S rRNA genes was 33% ± 12%. The abundances of these genes showed large spatial variations. dddP genes showed more variation in abundances than dmdA genes. Multidimensional analysis based on the abundances of DMSP degradation genes and environmental factors revealed that the distribution pattern of these genes was influenced by chlorophyll a concentrations and temperatures. dddP genes, dmdA subclade C/2 genes, and dmdA subclade D genes exhibited significant correlations with the marine Roseobacter clade, SAR11 subgroup Ib, and SAR11 subgroup Ia, respectively. SAR11 subgroups Ia and Ib, which possessed dmdA genes, were suggested to be the main potential DMSP consumers. The Roseobacter clade members possessing dddP genes in oligotrophic subtropical regions were possible DMS producers. These results suggest that DMSP degradation genes are abundant and widely distributed in the surface seawater and that the marine bacteria possessing these genes influence the degradation of DMSP and regulate the emissions of DMS in subtropical gyres of the Pacific Ocean.     

One carbon metabolism in SAR11 pelagic marine bacteria

The SAR11 Alphaproteobacteria are the most abundant heterotrophs in the oceans and are believed to play a major role in mineralizing marine dissolved organic carbon. Their genomes are among the smallest known for free-living heterotrophic cells, raising questions about how they successfully utilize complex organic matter with a limited metabolic repertoire. Here we show that conserved genes in SAR11 subgroup Ia (Candidatus Pelagibacter ubique) genomes encode pathways for the oxidation of a variety of one-carbon compounds and methyl functional groups from methylated compounds. These pathways were predicted to produce energy by tetrahydrofolate (THF)-mediated oxidation, but not to support the net assimilation of biomass from C1 compounds. Measurements of cellular ATP content and the oxidation of (14)C-labeled compounds to (14)CO(2) indicated that methanol, formaldehyde, methylamine, and methyl groups from glycine betaine (GBT), trimethylamine (TMA), trimethylamine N-oxide (TMAO), and dimethylsulfoniopropionate (DMSP) were oxidized by axenic cultures of the SAR11 strain Ca. P. ubique HTCC1062. Analyses of metagenomic data showed that genes for C1 metabolism occur at a high frequency in natural SAR11 populations. In short term incubations, natural communities of Sargasso Sea microbial plankton expressed a potential for the oxidation of (14)C-labeled formate, formaldehyde, methanol and TMAO that was similar to cultured SAR11 cells and, like cultured SAR11 cells, incorporated a much larger percentage of pyruvate and glucose (27-35%) than of C1 compounds (2-6%) into biomass. Collectively, these genomic, cellular and environmental data show a surprising capacity for demethylation and C1 oxidation in SAR11 cultures and in natural microbial communities dominated by SAR11, and support the conclusion that C1 oxidation might be a significant conduit by which dissolved organic carbon is recycled to CO(2) in the upper ocean.     

Diversity and abundance of “Pelagibacterales” (SAR11) in the Baltic Sea salinity gradient

The candidate order “Pelagibacterales” (SAR11) is one of the most abundant bacterial orders in ocean surface waters and, periodically, in freshwater lakes. The presence of several stable phylogenetic lineages comprising “Pelagibacterales” correlates with the physico-chemical parameters in aquatic environments. A previous amplicon sequencing study covering the bacterial community in the salinity gradient of the Baltic Sea suggested that pelagibacteral subclade SAR11-I was replaced by SAR11-IIIa in the mesohaline region of the Baltic Sea. In this current study, we investigated the cellular abundances of “Pelagibacterales” subclades along the Baltic Sea salinity gradient using catalyzed reporter deposition fluorescence in situ hybridization (CARD-FISH). The results obtained with a newly designed probe, which exclusively detected SAR11-IIIa, were compared to CARD-FISH abundances of the marine SAR11-I/II subclade and the freshwater lineage SAR11-IIIb (LD12). The results showed that SAR11-IIIa was abundant in oligohaline–mesohaline conditions (salinities 2.7–13.3), with maximal abundances at a salinity of 7 (up to 35% of total Bacteria, quantified with a universal bacterial probe EUB). As expected, SAR11-I/II was abundant (27% of EUB) in the marine parts of the Baltic Sea, whereas counts of the freshwater lineage SAR11-IIIb were below the detection limit at all stations. The shift from SAR11-IIIa to SAR11-I/II was confirmed in the vertical salinity gradient in the deeper basins of the Baltic Sea. These findings were consistent with an overlapping but defined distribution of SAR11-I/II and SAR11-IIIa in the salinity gradient of the Baltic Sea and suggested the adaptation of SAR11-IIIa for growth and survival in mesohaline conditions.     

Comparable light stimulation of organic nutrient uptake by SAR11 and Prochlorococcus in the North Atlantic subtropical gyre

Subtropical oceanic gyres are the most extensive biomes on Earth where SAR11 and Prochlorococcus bacterioplankton numerically dominate the surface waters depleted in inorganic macronutrients as well as in dissolved organic matter. In such nutrient poor conditions bacterioplankton could become photoheterotrophic, that is, potentially enhance uptake of scarce organic molecules using the available solar radiation to energise appropriate transport systems. Here, we assessed the photoheterotrophy of the key microbial taxa in the North Atlantic oligotrophic gyre and adjacent regions using ³³P-ATP, ³H-ATP and ³⁵S-methionine tracers. Light-stimulated uptake of these substrates was assessed in two dominant bacterioplankton groups discriminated by flow cytometric sorting of tracer-labelled cells and identified using catalysed reporter deposition fluorescence in situ hybridisation. One group of cells, encompassing 48% of all bacterioplankton, were identified as members of the SAR11 clade, whereas the other group (24% of all bacterioplankton) was Prochlorococcus. When exposed to light, SAR11 cells took 31% more ATP and 32% more methionine, whereas the Prochlorococcus cells took 33% more ATP and 34% more methionine. Other bacterioplankton did not demonstrate light stimulation. Thus, the SAR11 and Prochlorococcus groups, with distinctly different light-harvesting mechanisms, used light equally to enhance, by approximately one-third, the uptake of different types of organic molecules. Our findings indicate the significance of light-driven uptake of essential organic nutrients by the dominant bacterioplankton groups in the surface waters of one of the less productive, vast regions of the world''s oceans—the oligotrophic North Atlantic subtropical gyre.     

Efficient dilution-to-extinction isolation of novel virus–host model systems for fastidious heterotrophic bacteria

Microbes and their associated viruses are key drivers of biogeochemical processes in marine and soil biomes. While viruses of phototrophic cyanobacteria are well-represented in model systems, challenges of isolating marine microbial heterotrophs and their viruses have hampered experimental approaches to quantify the importance of viruses in nutrient recycling. A resurgence in cultivation efforts has improved the availability of fastidious bacteria for hypothesis testing, but this has not been matched by similar efforts to cultivate their associated bacteriophages. Here, we describe a high-throughput method for isolating important virus–host systems for fastidious heterotrophic bacteria that couples advances in culturing of hosts with sequential enrichment and isolation of associated phages. Applied to six monthly samples from the Western English Channel, we first isolated one new member of the globally dominant bacterial SAR11 clade and three new members of the methylotrophic bacterial clade OM43. We used these as bait to isolate 117 new phages, including the first known siphophage-infecting SAR11, and the first isolated phage for OM43. Genomic analyses of 13 novel viruses revealed representatives of three new viral genera, and infection assays showed that the viruses infecting SAR11 have ecotype-specific host ranges. Similar to the abundant human-associated phage ɸCrAss001, infection dynamics within the majority of isolates suggested either prevalent lysogeny or chronic infection, despite a lack of associated genes, or host phenotypic bistability with lysis putatively maintained within a susceptible subpopulation. Broader representation of important virus–host systems in culture collections and genomic databases will improve both our understanding of virus–host interactions, and accuracy of computational approaches to evaluate ecological patterns from metagenomic data.Subject terms: Bacteriophages, Microbial ecology     

Structure,fluctuation and magnitude of a natural grassland soil metagenome

The soil ecosystem is critical for human health, affecting aspects of the environment from key agricultural and edaphic parameters to critical influence on climate change. Soil has more unknown biodiversity than any other ecosystem. We have applied diverse DNA extraction methods coupled with high throughput pyrosequencing to explore 4.88 × 10⁹ bp of metagenomic sequence data from the longest continually studied soil environment (Park Grass experiment at Rothamsted Research in the UK). Results emphasize important DNA extraction biases and unexpectedly low seasonal and vertical soil metagenomic functional class variations. Clustering-based subsystems and carbohydrate metabolism had the largest quantity of annotated reads assigned although <50% of reads were assigned at an E value cutoff of 10⁻⁵. In addition, with the more detailed subsystems, cAMP signaling in bacteria (3.24±0.27% of the annotated reads) and the Ton and Tol transport systems (1.69±0.11%) were relatively highly represented. The most highly represented genome from the database was that for a Bradyrhizobium species. The metagenomic variance created by integrating natural and methodological fluctuations represents a global picture of the Rothamsted soil metagenome that can be used for specific questions and future inter-environmental metagenomic comparisons. However, only 1% of annotated sequences correspond to already sequenced genomes at 96% similarity and E values of <10⁻⁵, thus, considerable genomic reconstructions efforts still have to be performed.