Species-resolved, single-cell respiration rates reveal dominance of sulfate reduction in a deep continental subsurface ecosystem

Rates of microbial processes are fundamental to understanding the significance of microbial impacts on environmental chemical cycling. However, it is often difficult to quantify rates or to link processes to specific taxa or individual cells, especially in environments where there are few cultured representatives with known physiology. Here, we describe the use of the redox-enzyme-sensitive molecular probe RedoxSensor™ Green to measure rates of anaerobic electron transfer physiology (i.e., sulfate reduction and methanogenesis) in individual cells and link those measurements to genomic sequencing of the same single cells. We used this method to investigate microbial activity in hot, anoxic, low-biomass (~103 cells mL-1) groundwater of the Death Valley Regional Flow System, California. Combining this method with electron donor amendment experiments and metatranscriptomics confirmed that the abundant spore formers including Candidatus Desulforudis audaxviator were actively reducing sulfate in this environment, most likely with acetate and hydrogen as electron donors. Using this approach, we measured environmental sulfate reduction rates at 0.14 to 26.9 fmol cell-1 h-1. Scaled to volume, this equates to a bulk environmental rate of ~103 pmol sulfate L-1 d-1, similar to potential rates determined with radiotracer methods. Despite methane in the system, there was no evidence for active microbial methanogenesis at the time of sampling. Overall, this method is a powerful tool for estimating species-resolved, single-cell rates of anaerobic metabolism in low-biomass environments while simultaneously linking genomes to phenomes at the single-cell level. We reveal active elemental cycling conducted by several species, with a large portion attributable to Ca. Desulforudis audaxviator.

Interplay between autotrophic and heterotrophic prokaryotic metabolism in the bathypelagic realm revealed by metatranscriptomic analyses

BACKGROUND

Heterotrophic microbes inhabiting the dark ocean largely depend on the settling of organic matter from the sunlit ocean. However, this sinking of organic materials is insufficient to cover their demand for energy and alternative sources such as chemoautotrophy have been proposed. Reduced sulfur compounds, such as thiosulfate, are a potential energy source for both auto- and heterotrophic marine prokaryotes.

METHODS

Seawater samples were collected from Labrador Sea Water (LSW, ~ 2000 m depth) in the North Atlantic and incubated in the dark at in situ temperature unamended, amended with 1 µM thiosulfate, or with 1 µM thiosulfate plus 10 µM glucose and 10 µM acetate (thiosulfate plus dissolved organic matter, DOM). Inorganic carbon fixation was measured in the different treatments and samples for metatranscriptomic analyses were collected after 1 h and 72 h of incubation.

RESULTS

Amendment of LSW with thiosulfate and thiosulfate plus DOM enhanced prokaryotic inorganic carbon fixation. The energy generated via chemoautotrophy and heterotrophy in the amended prokaryotic communities was used for the biosynthesis of glycogen and phospholipids as storage molecules. The addition of thiosulfate stimulated unclassified bacteria, sulfur-oxidizing Deltaproteobacteria (SAR324 cluster bacteria), Epsilonproteobacteria (Sulfurimonas sp.), and Gammaproteobacteria (SUP05 cluster bacteria), whereas, the amendment with thiosulfate plus DOM stimulated typically copiotrophic Gammaproteobacteria (closely related to Vibrio sp. and Pseudoalteromonas sp.).

CONCLUSIONS

The gene expression pattern of thiosulfate utilizing microbes specifically of genes involved in energy production via sulfur oxidation and coupled to CO2 fixation pathways coincided with the change in the transcriptional profile of the heterotrophic prokaryotic community (genes involved in promoting energy storage), suggesting a fine-tuned metabolic interplay between chemoautotrophic and heterotrophic microbes in the dark ocean. Video Abstract.

Life strategies for Aminicenantia in subseafloor oceanic crust

After decades studying the microbial "deep biosphere" in subseafloor oceanic crust, the growth and life strategies in this anoxic, low energy habitat remain poorly described. Using both single cell genomics and metagenomics, we reveal the life strategies of two distinct lineages of uncultivated Aminicenantia bacteria from the basaltic subseafloor oceanic crust of the eastern flank of the Juan de Fuca Ridge. Both lineages appear adapted to scavenge organic carbon, as each have genetic potential to catabolize amino acids and fatty acids, aligning with previous Aminicenantia reports. Given the organic carbon limitation in this habitat, seawater recharge and necromass may be important carbon sources for heterotrophic microorganisms inhabiting the ocean crust. Both lineages generate ATP via several mechanisms including substrate-level phosphorylation, anaerobic respiration, and electron bifurcation driving an Rnf ion translocation membrane complex. Genomic comparisons suggest these Aminicenantia transfer electrons extracellularly, perhaps to iron or sulfur oxides consistent with mineralogy of this site. One lineage, called JdFR-78, has small genomes that are basal to the Aminicenantia class and potentially use "primordial" siroheme biosynthetic intermediates for heme synthesis, suggesting this lineage retain characteristics of early evolved life. Lineage JdFR-78 contains CRISPR-Cas defenses to evade viruses, while other lineages contain prophage that may help prevent super-infection or no detectable viral defenses. Overall, genomic evidence points to Aminicenantia being well adapted to oceanic crust environments by taking advantage of simple organic molecules and extracellular electron transport.

Synthase-selected sorting approach identifies a beta-lactone synthase in a nudibranch symbiotic bacterium

BACKGROUND

Nudibranchs comprise a group of > 6000 marine soft-bodied mollusk species known to use secondary metabolites (natural products) for chemical defense. The full diversity of these metabolites and whether symbiotic microbes are responsible for their synthesis remains unexplored. Another issue in searching for undiscovered natural products is that computational analysis of genomes of uncultured microbes can result in detection of novel biosynthetic gene clusters; however, their in vivo functionality is not guaranteed which limits further exploration of their pharmaceutical or industrial potential. To overcome these challenges, we used a fluorescent pantetheine probe, which produces a fluorescent CoA-analog employed in biosynthesis of secondary metabolites, to label and capture bacterial symbionts actively producing these compounds in the mantle of the nudibranch Doriopsilla fulva.

RESULTS

We recovered the genome of Candidatus Doriopsillibacter californiensis from the Ca. Tethybacterales order, an uncultured lineage of sponge symbionts not found in nudibranchs previously. It forms part of the core skin microbiome of D. fulva and is nearly absent in its internal organs. We showed that crude extracts of D. fulva contained secondary metabolites that were consistent with the presence of a beta-lactone encoded in Ca. D. californiensis genome. Beta-lactones represent an underexplored group of secondary metabolites with pharmaceutical potential that have not been reported in nudibranchs previously.

CONCLUSIONS

Altogether, this study shows how probe-based, targeted sorting approaches can capture bacterial symbionts producing secondary metabolites in vivo. Video Abstract.

A compendium of bacterial and archaeal single-cell amplified genomes from oxygen deficient marine waters

Oxygen-deficient marine waters referred to as oxygen minimum zones (OMZs) or anoxic marine zones (AMZs) are common oceanographic features. They host both cosmopolitan and endemic microorganisms adapted to low oxygen conditions. Microbial metabolic interactions within OMZs and AMZs drive coupled biogeochemical cycles resulting in nitrogen loss and climate active trace gas production and consumption. Global warming is causing oxygen-deficient waters to expand and intensify. Therefore, studies focused on microbial communities inhabiting oxygen-deficient regions are necessary to both monitor and model the impacts of climate change on marine ecosystem functions and services. Here we present a compendium of 5,129 single-cell amplified genomes (SAGs) from marine environments encompassing representative OMZ and AMZ geochemical profiles. Of these, 3,570 SAGs have been sequenced to different levels of completion, providing a strain-resolved perspective on the genomic content and potential metabolic interactions within OMZ and AMZ microbiomes. Hierarchical clustering confirmed that samples from similar oxygen concentrations and geographic regions also had analogous taxonomic compositions, providing a coherent framework for comparative community analysis.

Flotillin-associated rhodopsin (FArhodopsin), a widespread paralog of proteorhodopsin in aquatic bacteria with streamlined genomes

Microbial rhodopsins are found more than once in a single genome (paralogs) often have different functions. We screened a large dataset of open ocean single-amplified genomes (SAGs) for co-occurrences of multiple rhodopsin genes. Many such cases were found among Pelagibacterales (SAR11), HIMB59, and the Gammaproteobacteria Pseudothioglobus SAGs. These genomes always had a bona fide proteorhodopsin and a separate cluster of genes containing a second rhodopsin associated with a predicted flotillin coding gene and have thus been named flotillin-associated rhodopsins (FArhodopsins). Although they are members of the proteorhodopsin protein family, they form a separate clade within that family and are quite divergent from known proton-pumping proteorhodopsins. They contain either DTT, DTL, or DNI motifs in their key functional amino acids. FArhodopsins are mainly associated with the lower layers of the epipelagic zone. All marine FArhodopsins had the retinal binding lysine, but we found relatives in freshwater metagenomes lacking this key amino acid. AlphaFold predictions of marine FArhodopsins indicate that their retinal pocket might be very reduced or absent, hinting that they are retinal-less. Freshwater FArhodopsins were more diverse than marine ones, but we could not determine if there were other rhodopsins in the genome due to the lack of SAGs or isolates. Although the function of FArhodopsins could not be established, their conserved genomic context indicated involvement in the formation of membrane microdomains. The conservation of FArhodopsins in diverse and globally abundant microorganisms suggests that they may be important in the adaptation to the twilight zone of aquatic environments.IMPORTANCERhodopsins have been shown to play a key role in the ecology of aquatic microbes. Here, we describe a group of widespread rhodopsins in aquatic microbes associated with dim light conditions. Their characteristic genomic context found in both marine and freshwater environments indicates a novel potential involvement in membrane microstructure that could be important for the function of the coexisting proteorhodopsin proton pumps. The absence or reduction of the retinal binding pocket points to a drastically different physiological role.

Novel integrative elements and genomic plasticity in ocean ecosystems

Horizontal gene transfer accelerates microbial evolution. The marine picocyanobacterium Prochlorococcus exhibits high genomic plasticity, yet the underlying mechanisms are elusive. Here, we report a novel family of DNA transposons-"tycheposons"-some of which are viral satellites while others carry cargo, such as nutrient-acquisition genes, which shape the genetic variability in this globally abundant genus. Tycheposons share distinctive mobile-lifecycle-linked hallmark genes, including a deep-branching site-specific tyrosine recombinase. Their excision and integration at tRNA genes appear to drive the remodeling of genomic islands-key reservoirs for flexible genes in bacteria. In a selection experiment, tycheposons harboring a nitrate assimilation cassette were dynamically gained and lost, thereby promoting chromosomal rearrangements and host adaptation. Vesicles and phage particles harvested from seawater are enriched in tycheposons, providing a means for their dispersal in the wild. Similar elements are found in microbes co-occurring with Prochlorococcus, suggesting a common mechanism for microbial diversification in the vast oligotrophic oceans.

Decoupling of respiration rates and abundance in marine prokaryoplankton

The ocean-atmosphere exchange of CO2 largely depends on the balance between marine microbial photosynthesis and respiration. Despite vast taxonomic and metabolic diversity among marine planktonic bacteria and archaea (prokaryoplankton)1-3, their respiration usually is measured in bulk and treated as a 'black box' in global biogeochemical models4; this limits the mechanistic understanding of the global carbon cycle. Here, using a technology for integrated phenotype analyses and genomic sequencing of individual microbial cells, we show that cell-specific respiration rates differ by more than 1,000× among prokaryoplankton genera. The majority of respiration was found to be performed by minority members of prokaryoplankton (including the Roseobacter cluster), whereas cells of the most prevalent lineages (including Pelagibacter and SAR86) had extremely low respiration rates. The decoupling of respiration rates from abundance among lineages, elevated counts of proteorhodopsin transcripts in Pelagibacter and SAR86 cells and elevated respiration of SAR86 at night indicate that proteorhodopsin-based phototrophy3,5-7 probably constitutes an important source of energy to prokaryoplankton and may increase growth efficiency. These findings suggest that the dependence of prokaryoplankton on respiration and remineralization of phytoplankton-derived organic carbon into CO2 for its energy demands and growth may be lower than commonly assumed and variable among lineages.

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

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

Single cell genome sequencing of laboratory mouse microbiota improves taxonomic and functional resolution of this model microbial community

Laboratory mice are widely studied as models of mammalian biology, including the microbiota. However, much of the taxonomic and functional diversity of the mouse gut microbiome is missed in current metagenomic studies, because genome databases have not achieved a balanced representation of the diverse members of this ecosystem. Towards solving this problem, we used flow cytometry and low-coverage sequencing to capture the genomes of 764 single cells from the stool of three laboratory mice. From these, we generated 298 high-coverage microbial genome assemblies, which we annotated for open reading frames and phylogenetic placement. These genomes increase the gene catalog and phylogenetic breadth of the mouse microbiota, adding 135 novel species with the greatest increase in diversity to the Muribaculaceae and Bacteroidaceae families. This new diversity also improves the read mapping rate, taxonomic classifier performance, and gene detection rate of mouse stool metagenomes. The novel microbial functions revealed through our single-cell genomes highlight previously invisible pathways that may be important for life in the murine gastrointestinal tract.

Oceanic Crustal Fluid Single Cell Genomics Complements Metagenomic and Metatranscriptomic Surveys With Orders of Magnitude Less Sample Volume

Fluids circulating through oceanic crust play important roles in global biogeochemical cycling mediated by their microbial inhabitants, but studying these sites is challenged by sampling logistics and low biomass. Borehole observatories installed at the North Pond study site on the western flank of the Mid-Atlantic Ridge have enabled investigation of the microbial biosphere in cold, oxygenated basaltic oceanic crust. Here we test a methodology that applies redox-sensitive fluorescent molecules for flow cytometric sorting of cells for single cell genomic sequencing from small volumes of low biomass (approximately 103 cells ml-1) crustal fluid. We compare the resulting genomic data to a recently published paired metagenomic and metatranscriptomic analysis from the same site. Even with low coverage genome sequencing, sorting cells from less than one milliliter of crustal fluid results in similar interpretation of dominant taxa and functional profiles as compared to 'omics analysis that typically filter orders of magnitude more fluid volume. The diverse community dominated by Gammaproteobacteria, Bacteroidetes, Desulfobacterota, Alphaproteobacteria, and Zetaproteobacteria, had evidence of autotrophy and heterotrophy, a variety of nitrogen and sulfur cycling metabolisms, and motility. Together, results indicate fluorescence activated cell sorting methodology is a powerful addition to the toolbox for the study of low biomass systems or at sites where only small sample volumes are available for analysis.

Phylogenetically and functionally diverse microorganisms reside under the Ross Ice Shelf

Throughout coastal Antarctica, ice shelves separate oceanic waters from sunlight by hundreds of meters of ice. Historical studies have detected activity of nitrifying microorganisms in oceanic cavities below permanent ice shelves. However, little is known about the microbial composition and pathways that mediate these activities. In this study, we profiled the microbial communities beneath the Ross Ice Shelf using a multi-omics approach. Overall, beneath-shelf microorganisms are of comparable abundance and diversity, though distinct composition, relative to those in the open meso- and bathypelagic ocean. Production of new organic carbon is likely driven by aerobic lithoautotrophic archaea and bacteria that can use ammonium, nitrite, and sulfur compounds as electron donors. Also enriched were aerobic organoheterotrophic bacteria capable of degrading complex organic carbon substrates, likely derived from in situ fixed carbon and potentially refractory organic matter laterally advected by the below-shelf waters. Altogether, these findings uncover a taxonomically distinct microbial community potentially adapted to a highly oligotrophic marine environment and suggest that ocean cavity waters are primarily chemosynthetically-driven systems.

Clean room microbiome complexity impacts planetary protection bioburden

BACKGROUND

The Spacecraft Assembly Facility (SAF) at the NASA's Jet Propulsion Laboratory is the primary cleanroom facility used in the construction of some of the planetary protection (PP)-sensitive missions developed by NASA, including the Mars 2020 Perseverance Rover that launched in July 2020. SAF floor samples (n=98) were collected, over a 6-month period in 2016 prior to the construction of the Mars rover subsystems, to better understand the temporal and spatial distribution of bacterial populations (total, viable, cultivable, and spore) in this unique cleanroom.

RESULTS

Cleanroom samples were examined for total (living and dead) and viable (living only) microbial populations using molecular approaches and cultured isolates employing the traditional NASA standard spore assay (NSA), which predominantly isolated spores. The 130 NSA isolates were represented by 16 bacterial genera, of which 97% were identified as spore-formers via Sanger sequencing. The most spatially abundant isolate was Bacillus subtilis, and the most temporally abundant spore-former was Virgibacillus panthothenticus. The 16S rRNA gene-targeted amplicon sequencing detected 51 additional genera not found in the NSA method. The amplicon sequencing of the samples treated with propidium monoazide (PMA), which would differentiate between viable and dead organisms, revealed a total of 54 genera: 46 viable non-spore forming genera and 8 viable spore forming genera in these samples. The microbial diversity generated by the amplicon sequencing corresponded to ~86% non-spore-formers and ~14% spore-formers. The most common spatially distributed genera were Sphinigobium, Geobacillus, and Bacillus whereas temporally distributed common genera were Acinetobacter, Geobacilllus, and Bacillus. Single-cell genomics detected 6 genera in the sample analyzed, with the most prominent being Acinetobacter.

CONCLUSION

This study clearly established that detecting spores via NSA does not provide a complete assessment for the cleanliness of spacecraft-associated environments since it failed to detect several PP-relevant genera that were only recovered via molecular methods. This highlights the importance of a methodological paradigm shift to appropriately monitor bioburden in cleanrooms for not only the aeronautical industry but also for pharmaceutical, medical industries, etc., and the need to employ molecular sequencing to complement traditional culture-based assays. Video abstract.

Genomic Variation Influences Methanothermococcus Fitness in Marine Hydrothermal Systems

Hydrogenotrophic methanogens are ubiquitous chemoautotrophic archaea inhabiting globally distributed deep-sea hydrothermal vent ecosystems and associated subseafloor niches within the rocky subseafloor, yet little is known about how they adapt and diversify in these habitats. To determine genomic variation and selection pressure within methanogenic populations at vents, we examined five Methanothermococcus single cell amplified genomes (SAGs) in conjunction with 15 metagenomes and 10 metatranscriptomes from venting fluids at two geochemically distinct hydrothermal vent fields on the Mid-Cayman Rise in the Caribbean Sea. We observed that some Methanothermococcus lineages and their transcripts were more abundant than others in individual vent sites, indicating differential fitness among lineages. The relative abundances of lineages represented by SAGs in each of the samples matched phylogenetic relationships based on single-copy universal genes, and genes related to nitrogen fixation and the CRISPR/Cas immune system were among those differentiating the clades. Lineages possessing these genes were less abundant than those missing that genomic region. Overall, patterns in nucleotide variation indicated that the population dynamics of Methanothermococcus were not governed by clonal expansions or selective sweeps, at least in the habitats and sampling times included in this study. Together, our results show that although specific lineages of Methanothermococcus co-exist in these habitats, some outcompete others, and possession of accessory metabolic functions does not necessarily provide a fitness advantage in these habitats in all conditions. This work highlights the power of combining single-cell, metagenomic, and metatranscriptomic datasets to determine how evolution shapes microbial abundance and diversity in hydrothermal vent ecosystems.

Synthase-Selective Exploration of a Tunicate Microbiome by Activity-Guided Single-Cell Genomics

While thousands of environmental metagenomes have been mined for the presence of novel biosynthetic gene clusters, such computational predictions do not provide evidence of their in vivo biosynthetic functionality. Using fluorescent in situ enzyme assay targeting carrier proteins common to polyketide (PKS) and nonribosomal peptide synthetases (NRPS), we applied fluorescence-activated cell sorting to tunicate microbiome to enrich for microbes with active secondary metabolic capabilities. Single-cell genomics uncovered the genetic basis for a wide biosynthetic diversity in the enzyme-active cells and revealed a member of marine Oceanospirillales harboring a novel NRPS gene cluster with high similarity to phylogenetically distant marine and terrestrial bacteria. Interestingly, this synthase belongs to a larger class of siderophore biosynthetic gene clusters commonly associated with pestilence and disease. This demonstrates activity-guided single-cell genomics as a tool to guide novel biosynthetic discovery.

Genomic Analysis of the Yet-Uncultured Binatota Reveals Broad Methylotrophic, Alkane-Degradation, and Pigment Production Capacities

The recent leveraging of genome-resolved metagenomics has generated an enormous number of genomes from novel uncultured microbial lineages yet left many clades undescribed. Here, we present a global analysis of genomes belonging to Binatota (UBP10), a globally distributed, yet-uncharacterized bacterial phylum. All orders in Binatota encoded the capacity for aerobic methylotrophy using methanol, methylamine, sulfomethanes, and chloromethanes as the substrates. Methylotrophy in Binatota was characterized by order-specific substrate degradation preferences, as well as extensive metabolic versatility, i.e., the utilization of diverse sets of genes, pathways, and combinations to achieve a specific metabolic goal. The genomes also encoded multiple alkane hydroxylases and monooxygenases, potentially enabling growth on a wide range of alkanes and fatty acids. Pigmentation is inferred from a complete pathway for carotenoids (lycopene, β- and γ-carotenes, xanthins, chlorobactenes, and spheroidenes) production. Further, the majority of genes involved in bacteriochlorophyll a, c, and d biosynthesis were identified, although absence of key genes and failure to identify a photosynthetic reaction center preclude proposing phototrophic capacities. Analysis of 16S rRNA databases showed the preferences of Binatota to terrestrial and freshwater ecosystems, hydrocarbon-rich habitats, and sponges, supporting their potential role in mitigating methanol and methane emissions, breakdown of alkanes, and their association with sponges. Our results expand the lists of methylotrophic, aerobic alkane-degrading, and pigment-producing lineages. We also highlight the consistent encountering of incomplete biosynthetic pathways in microbial genomes, a phenomenon necessitating careful assessment when assigning putative functions based on a set-threshold of pathway completion.IMPORTANCE A wide range of microbial lineages remain uncultured, yet little is known regarding their metabolic capacities, physiological preferences, and ecological roles in various ecosystems. We conducted a thorough comparative genomic analysis of 108 genomes belonging to the Binatota (UBP10), a globally distributed, yet-uncharacterized bacterial phylum. We present evidence that members of the order Binatota specialize in methylotrophy and identify an extensive repertoire of genes and pathways mediating the oxidation of multiple one-carbon (C1) compounds in Binatota genomes. The occurrence of multiple alkane hydroxylases and monooxygenases in these genomes was also identified, potentially enabling growth on a wide range of alkanes and fatty acids. Pigmentation is inferred from a complete pathway for carotenoids production. We also report on the presence of incomplete chlorophyll biosynthetic pathways in all genomes and propose several evolutionary-grounded scenarios that could explain such a pattern. Assessment of the ecological distribution patterns of the Binatota indicates preference of its members to terrestrial and freshwater ecosystems characterized by high methane and methanol emissions, as well as multiple hydrocarbon-rich habitats and marine sponges.

Diversification of methanogens into hyperalkaline serpentinizing environments through adaptations to minimize oxidant limitation

Metagenome assembled genomes (MAGs) and single amplified genomes (SAGs) affiliated with two distinct Methanobacterium lineages were recovered from subsurface fracture waters of the Samail Ophiolite, Sultanate of Oman. Lineage Type I was abundant in waters with circumneutral pH, whereas lineage Type II was abundant in hydrogen rich, hyperalkaline waters. Type I encoded proteins to couple hydrogen oxidation to CO₂ reduction, typical of hydrogenotrophic methanogens. Surprisingly, Type II, which branched from the Type I lineage, lacked homologs of two key oxidative [NiFe]-hydrogenases. These functions were presumably replaced by formate dehydrogenases that oxidize formate to yield reductant and cytoplasmic CO₂ via a pathway that was unique among characterized Methanobacteria, allowing cells to overcome CO₂/oxidant limitation in high pH waters. This prediction was supported by microcosm-based radiotracer experiments that showed significant biological methane generation from formate, but not bicarbonate, in waters where the Type II lineage was detected in highest relative abundance. Phylogenetic analyses and variability in gene content suggested that recent and ongoing diversification of the Type II lineage was enabled by gene transfer, loss, and transposition. These data indicate that selection imposed by CO₂/oxidant availability drove recent methanogen diversification into hyperalkaline waters that are heavily impacted by serpentinization.

Single Cell Genomics Reveals Viruses Consumed by Marine Protists

The predominant model of the role of viruses in the marine trophic web is that of the “viral shunt,” where viral infection funnels a substantial fraction of the microbial primary and secondary production back to the pool of dissolved organic matter. Here, we analyzed the composition of non-eukaryotic DNA associated with individual cells of small, planktonic protists in the Gulf of Maine (GoM) and the Mediterranean Sea. We found viral DNA associated with a substantial fraction cells from the GoM (51%) and the Mediterranean Sea (35%). While Mediterranean SAGs contained a larger proportion of cells containing bacterial sequences (49%), a smaller fraction of cells contained bacterial sequences in the GoM (19%). In GoM cells, nearly identical bacteriophage and ssDNA virus sequences where found across diverse lineages of protists, suggesting many of these viruses are non-infective. The fraction of cells containing viral DNA varied among protistan lineages and reached 100% in Picozoa and Choanozoa. These two groups also contained significantly higher numbers of viral sequences than other identified taxa. We consider mechanisms that may explain the presence of viral DNA in protistan cells and conclude that protistan predation on free viral particles contributed to the observed patterns. These findings confirm prior experiments with protistan isolates and indicate that the viral shunt is complemented by a viral link in the marine microbial food web. This link may constitute a sink of viral particles in the ocean and has implications for the flow of carbon through the microbial food web.

Ancestral Absence of Electron Transport Chains in Patescibacteria and DPANN

Recent discoveries suggest that the candidate superphyla Patescibacteria and DPANN constitute a large fraction of the phylogenetic diversity of Bacteria and Archaea. Their small genomes and limited coding potential have been hypothesized to be ancestral adaptations to obligate symbiotic lifestyles. To test this hypothesis, we performed cell-cell association, genomic, and phylogenetic analyses on 4,829 individual cells of Bacteria and Archaea from 46 globally distributed surface and subsurface field samples. This confirmed the ubiquity and abundance of Patescibacteria and DPANN in subsurface environments, the small size of their genomes and cells, and the divergence of their gene content from other Bacteria and Archaea. Our analyses suggest that most Patescibacteria and DPANN in the studied subsurface environments do not form specific physical associations with other microorganisms. These data also suggest that their unusual genomic features and prevalent auxotrophies may be a result of ancestral, minimal cellular energy transduction mechanisms that lack respiration, thus relying solely on fermentation for energy conservation.

Insights into the dynamics between viruses and their hosts in a hot spring microbial mat

Our current knowledge of host-virus interactions in biofilms is limited to computational predictions based on laboratory experiments with a small number of cultured bacteria. However, natural biofilms are diverse and chiefly composed of uncultured bacteria and archaea with no viral infection patterns and lifestyle predictions described to date. Herein, we predict the first DNA sequence-based host-virus interactions in a natural biofilm. Using single-cell genomics and metagenomics applied to a hot spring mat of the Cone Pool in Mono County, California, we provide insights into virus-host range, lifestyle and distribution across different mat layers. Thirty-four out of 130 single cells contained at least one viral contig (26%), which, together with the metagenome-assembled genomes, resulted in detection of 59 viruses linked to 34 host species. Analysis of single-cell amplification kinetics revealed a lack of active viral replication on the single-cell level. These findings were further supported by mapping metagenomic reads from different mat layers to the obtained host-virus pairs, which indicated a low copy number of viral genomes compared to their hosts. Lastly, the metagenomic data revealed high layer specificity of viruses, suggesting limited diffusion to other mat layers. Taken together, these observations indicate that in low mobility environments with high microbial abundance, lysogeny is the predominant viral lifestyle, in line with the previously proposed "Piggyback-the-Winner" theory.

Hiding in Plain Sight: The Globally Distributed Bacterial Candidate Phylum PAUC34f

Bacterial candidate phylum PAUC34f was originally discovered in marine sponges and is widely considered to be composed of sponge symbionts. Here, we report 21 single amplified genomes (SAGs) of PAUC34f from a variety of environments, including the dark ocean, lake sediments, and a terrestrial aquifer. The diverse origins of the SAGs and the results of metagenome fragment recruitment suggest that some PAUC34f lineages represent relatively abundant, free-living cells in environments other than sponge microbiomes, including the deep ocean. Both phylogenetic and biogeographic patterns, as well as genome content analyses suggest that PAUC34f associations with hosts evolved independently multiple times, while free-living lineages of PAUC34f are distinct and relatively abundant in a wide range of environments.

Charting the Complexity of the Marine Microbiome through Single-Cell Genomics

Marine bacteria and archaea play key roles in global biogeochemistry. To improve our understanding of this complex microbiome, we employed single-cell genomics and a randomized, hypothesis-agnostic cell selection strategy to recover 12,715 partial genomes from the tropical and subtropical euphotic ocean. A substantial fraction of known prokaryoplankton coding potential was recovered from a single, 0.4 mL ocean sample, which indicates that genomic information disperses effectively across the globe. Yet, we found each genome to be unique, implying limited clonality within prokaryoplankton populations. Light harvesting and secondary metabolite biosynthetic pathways were numerous across lineages, highlighting the value of single-cell genomics to advance the identification of ecological roles and biotechnology potential of uncultured microbial groups. This genome collection enabled functional annotation and genus-level taxonomic assignments for >80% of individual metagenome reads from the tropical and subtropical surface ocean, thus offering a model to improve reference genome databases for complex microbiomes.

Correction: Single-cell genomics shedding light on marine Thaumarchaeota diversification

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Single Cell Genomics-Based Analysis of Gene Content and Expression of Prophages in a Diffuse-Flow Deep-Sea Hydrothermal System

Phage–host interactions likely play a major role in the composition and functioning of many microbiomes, yet remain poorly understood. Here, we employed single cell genomics to investigate phage–host interactions in a diffuse-flow, low-temperature hydrothermal vent that may be reflective of a broadly distributed biosphere in the subseafloor. We identified putative prophages in 13 of 126 sequenced single amplified genomes (SAGs), with no evidence for lytic infections, which is in stark contrast to findings in the surface ocean. Most were distantly related to known prophages, while their hosts included bacterial phyla Campylobacterota, Bacteroidetes, Chlorobi, Proteobacteria, Lentisphaerae, Spirochaetes, and Thermotogae. Our results suggest the predominance of lysogeny over lytic interaction in diffuse-flow, deep-sea hydrothermal vents, despite the high activity of the dominant Campylobacteria that would favor lytic infections. We show that some of the identified lysogens have co-evolved with their host over geological time scales and that their genes are transcribed in the environment. Functional annotations of lysogeny-related genes suggest involvement in horizontal gene transfer enabling host’s protection against toxic metals and antibacterial compounds.

Phylogenomics suggests oxygen availability as a driving force in Thaumarchaeota evolution

Ammonia-oxidizing archaea (AOA) of the phylum Thaumarchaeota are widespread in marine and terrestrial habitats, playing a major role in the global nitrogen cycle. However, their evolutionary history remains unexplored, which limits our understanding of their adaptation mechanisms. Here, our comprehensive phylogenomic tree of Thaumarchaeota supports three sequential events: origin of AOA from terrestrial non-AOA ancestors, colonization of the shallow ocean, and expansion to the deep ocean. Careful molecular dating suggests that these events coincided with the Great Oxygenation Event around 2300 million years ago (Mya), and oxygenation of the shallow and deep ocean around 800 and 635-560 Mya, respectively. The first transition was likely enabled by the gain of an aerobic pathway for energy production by ammonia oxidation and biosynthetic pathways for cobalamin and biotin that act as cofactors in aerobic metabolism. The first transition was also accompanied by the loss of dissimilatory nitrate and sulfate reduction, loss of oxygen-sensitive pyruvate oxidoreductase, which reduces pyruvate to acetyl-CoA, and loss of the Wood-Ljungdahl pathway for anaerobic carbon fixation. The second transition involved gain of a K⁺ transporter and of the biosynthetic pathway for ectoine, which may function as an osmoprotectant. The third transition was accompanied by the loss of the uvr system for repairing ultraviolet light-induced DNA lesions. We conclude that oxygen availability drove the terrestrial origin of AOA and their expansion to the photic and dark oceans, and that the stressors encountered during these events were partially overcome by gene acquisitions from Euryarchaeota and Bacteria, among other sources.

Single cell genomics yields a wide diversity of small planktonic protists across major ocean ecosystems

Marine planktonic protists are critical components of ocean ecosystems and are highly diverse. Molecular sequencing methods are being used to describe this diversity and reveal new associations and metabolisms that are important to how these ecosystems function. We describe here the use of the single cell genomics approach to sample and interrogate the diversity of the smaller (pico- and nano-sized) protists from a range of oceanic samples. We created over 900 single amplified genomes (SAGs) from 8 Tara Ocean samples across the Indian Ocean and the Mediterranean Sea. We show that flow cytometric sorting of single cells effectively distinguishes plastidic and aplastidic cell types that agree with our understanding of protist phylogeny. Yields of genomic DNA with PCR-identifiable 18S rRNA gene sequence from single cells was low (15% of aplastidic cell sorts, and 7% of plastidic sorts) and tests with alternate primers and comparisons to metabarcoding did not reveal phylogenetic bias in the major protist groups. There was little evidence of significant bias against or in favor of any phylogenetic group expected or known to be present. The four open ocean stations in the Indian Ocean had similar communities, despite ranging from 14°N to 20°S latitude, and they differed from the Mediterranean station. Single cell genomics of protists suggests that the taxonomic diversity of the dominant taxa found in only several hundreds of microliters of surface seawater is similar to that found in molecular surveys where liters of sample are filtered.

Emergence of trait variability through the lens of nitrogen assimilation in Prochlorococcus

Intraspecific trait variability has important consequences for the function and stability of marine ecosystems. Here we examine variation in the ability to use nitrate across hundreds of Prochlorococcus genomes to better understand the modes of evolution influencing intraspecific allocation of ecologically important functions. Nitrate assimilation genes are absent in basal lineages but occur at an intermediate frequency that is randomly distributed within recently emerged clades. The distribution of nitrate assimilation genes within clades appears largely governed by vertical inheritance, gene loss, and homologous recombination. By mapping this process onto a model of Prochlorococcus’ macroevolution, we propose that niche-constructing adaptive radiations and subsequent niche partitioning set the stage for loss of nitrate assimilation genes from basal lineages as they specialized to lower light levels. Retention of these genes in recently emerged lineages has likely been facilitated by selection as they sequentially partitioned into niches where nitrate assimilation conferred a fitness benefit.

Genomics insights into ecotype formation of ammonia-oxidizing archaea in the deep ocean

Various lineages of ammonia-oxidizing archaea (AOA) are present in deep waters, but the mechanisms that determine ecotype formation are obscure. We studied 18 high-quality genomes of the marine group I AOA lineages (alpha, gamma and delta) from the Mariana and Ogasawara trenches. The genomes of alpha AOA resembled each other, while those of gamma and delta lineages were more divergent and had even undergone insertion of some phage genes. The instability of the gamma and delta AOA genomes could be partially due to the loss of DNA polymerase B (polB) and methyladenine DNA glycosylase (tag) genes responsible for the repair of point mutations. The alpha AOA genomes harbour genes encoding a thrombospondin-like outer membrane structure that probably serves as a barrier to gene flow. Moreover, the gamma and alpha AOA lineages rely on vitamin B₁₂ -independent MetE and B₁₂ -dependent MetH, respectively, for methionine synthesis. The delta AOA genome contains genes involved in uptake of sugar and peptide perhaps for heterotrophic lifestyle. Our study provides insights into co-occurrence of cladogenesis and anagenesis in the formation of AOA ecotypes that perform differently in nitrogen and carbon cycling in dark oceans.

Single-cell genomics of co-sorted Nanoarchaeota suggests novel putative host associations and diversification of proteins involved in symbiosis

BACKGROUND

Nanoarchaeota are obligate symbionts of other Archaea first discovered 16 years ago, yet little is known about this largely uncultivated taxon. While Nanoarchaeota diversity has been detected in a variety of habitats using 16S rRNA gene surveys, genome sequences have been available for only three Nanoarchaeota and their hosts. The host range and adaptation of Nanoarchaeota to a wide range of environmental conditions has thus largely remained elusive. Single-cell genomics is an ideal approach to address these questions as Nanoarchaeota can be isolated while still attached to putative hosts, enabling the exploration of cell-cell interactions and fine-scale genomic diversity.

RESULTS

From 22 single amplified genomes (SAGs) from three hot springs in Yellowstone National Park, we derived a genome-based phylogeny of the phylum Nanoarchaeota, linking it to global 16S rRNA gene diversity. By exploiting sequencing of co-sorted tightly attached cells, we associated Nanoarchaeota with 6 novel putative hosts, 2 of which were found in multiple SAGs, and showed that the same host species may associate with multiple species of Nanoarchaeota. Comparison of single nucleotide polymorphisms (SNPs) within a population of Nanoarchaeota SAGs indicated that Nanoarchaeota attached to a single host cell in situ are likely clonal. In addition to an overall pattern of purifying selection, we found significantly higher densities of non-synonymous SNPs in hypothetical cell surface proteins, as compared to other functional categories. Genes implicated in interactions in other obligate microbe-microbe symbioses, including those encoding a cytochrome bd-I ubiquinol oxidase and a FlaJ/TadC homologue possibly involved in type IV pili production, also had relatively high densities of non-synonymous SNPs.

CONCLUSIONS

This population genetics study of Nanoarchaeota greatly expands the known potential host range of the phylum and hints at what genes may be involved in adaptation to diverse environments or different hosts. We provide the first evidence that Nanoarchaeota cells attached to the same host cell are clonal and propose a hypothesis for how clonality may occur despite diverse symbiont populations.

Single cell genomes of Prochlorococcus, Synechococcus, and sympatric microbes from diverse marine environments

Prochlorococcus and Synechococcus are the dominant primary producers in marine ecosystems and perform a significant fraction of ocean carbon fixation. These cyanobacteria interact with a diverse microbial community that coexists with them. Comparative genomics of cultivated isolates has helped address questions regarding patterns of evolution and diversity among microbes, but the fraction that can be cultivated is miniscule compared to the diversity in the wild. To further probe the diversity of these groups and extend the utility of reference sequence databases, we report a data set of single cell genomes for 489 Prochlorococcus, 50 Synechococcus, 9 extracellular virus particles, and 190 additional microorganisms from a diverse range of bacterial, archaeal, and viral groups. Many of these uncultivated single cell genomes are derived from samples obtained on GEOTRACES cruises and at well-studied oceanographic stations, each with extensive suites of physical, chemical, and biological measurements. The genomic data reported here greatly increases the number of available Prochlorococcus genomes and will facilitate studies on evolutionary biology, microbial ecology, and biological oceanography.

Major role of nitrite-oxidizing bacteria in dark ocean carbon fixation

Carbon fixation by chemoautotrophic microorganisms in the dark ocean has a major impact on global carbon cycling and ecological relationships in the ocean's interior, but the relevant taxa and energy sources remain enigmatic. We show evidence that nitrite-oxidizing bacteria affiliated with the Nitrospinae phylum are important in dark ocean chemoautotrophy. Single-cell genomics and community metagenomics revealed that Nitrospinae are the most abundant and globally distributed nitrite-oxidizing bacteria in the ocean. Metaproteomics and metatranscriptomics analyses suggest that nitrite oxidation is the main pathway of energy production in Nitrospinae. Microautoradiography, linked with catalyzed reporter deposition fluorescence in situ hybridization, indicated that Nitrospinae fix 15 to 45% of inorganic carbon in the mesopelagic western North Atlantic. Nitrite oxidation may have a greater impact on the carbon cycle than previously assumed.

Rokubacteria: Genomic Giants among the Uncultured Bacterial Phyla

Recent advances in single-cell genomic and metagenomic techniques have facilitated the discovery of numerous previously unknown, deep branches of the tree of life that lack cultured representatives. Many of these candidate phyla are composed of microorganisms with minimalistic, streamlined genomes lacking some core metabolic pathways, which may contribute to their resistance to growth in pure culture. Here we analyzed single-cell genomes and metagenome bins to show that the "Candidate phylum Rokubacteria," formerly known as SPAM, represents an interesting exception, by having large genomes (6-8 Mbps), high GC content (66-71%), and the potential for a versatile, mixotrophic metabolism. We also observed an unusually high genomic heterogeneity among individual Rokubacteria cells in the studied samples. These features may have contributed to the limited recovery of sequences of this candidate phylum in prior cultivation and metagenomic studies. Our analyses suggest that Rokubacteria are distributed globally in diverse terrestrial ecosystems, including soils, the rhizosphere, volcanic mud, oil wells, aquifers, and the deep subsurface, with no reports from marine environments to date.

Characterization of Two Cryptic Plasmids Isolated in Haiti from Clinical Vibrio cholerae Non-O1/Non-O139

We report the complete sequence of two novel plasmids, pSDH-1 and pSDH-2, isolated from clinical Vibrio cholerae non-O1/non-O139 during the early phase of the 2010 Haitian cholera epidemic. Plasmids were revealed by employing single-cell genomics and their genome content suggests self-mobilization and, for pSDH-2, a toxin-antitoxin (TA) system for plasmid stabilization was identified. The putative origin of replication of pSDH-2 was mapped suggesting it replicates following the ColE1 model of plasmid replication. pSDH-1 and pSDH-2 were widespread among environmental V. cholerae non-O1/non-O139 with variable prevalence in four Haitian Departments. pSDH-2 was the most common element, either alone or with pSDH-1. The two plasmids detection adds to the composite scenario of mobile genetic elements (MGEs) observed in V. cholerae in Haiti. The role these small cryptic plasmids circulating in Vibrio spp. play in bacterial fitness or pathogenicity merits further investigation.

Genomic Comparison of Two Family-Level Groups of the Uncultivated NAG1 Archaeal Lineage from Chemically and Geographically Disparate Hot Springs

Recent progress based on single-cell genomics and metagenomic investigations of archaea in a variety of extreme environments has led to significant advances in our understanding of the diversity, evolution, and metabolic potential of archaea, yet the vast majority of archaeal diversity remains undersampled. In this work, we coordinated single-cell genomics with metagenomics in order to construct a near-complete genome from a deeply branching uncultivated archaeal lineage sampled from Great Boiling Spring (GBS) in the U.S. Great Basin, Nevada. This taxon is distantly related (distinct families) to an archaeal genome, designated "Novel Archaeal Group 1" (NAG1), which was extracted from a metagenome recovered from an acidic iron spring in Yellowstone National Park (YNP). We compared the metabolic predictions of the NAG1 lineage to better understand how these archaea could inhabit such chemically distinct environments. Similar to the NAG1 population previously studied in YNP, the NAG1 population from GBS is predicted to utilize proteins as a primary carbon source, ferment simple carbon sources, and use oxygen as a terminal electron acceptor under oxic conditions. However, GBS NAG1 populations contained distinct genes involved in central carbon metabolism and electron transfer, including nitrite reductase, which could confer the ability to reduce nitrite under anaerobic conditions. Despite inhabiting chemically distinct environments with large variations in pH, GBS NAG1 populations shared many core genomic and metabolic features with the archaeon identified from YNP, yet were able to carve out a distinct niche at GBS.

Minimum information about a single amplified genome (MISAG) and a metagenome-assembled genome (MIMAG) of bacteria and archaea

Standards for sequencing the microbial 'uncultivated majority', namely bacterial and archaeal single-cell genome sequences, and genome sequences from metagenomic datasets, are proposed.

We present two standards developed by the Genomic Standards Consortium (GSC) for reporting bacterial and archaeal genome sequences. Both are extensions of the Minimum Information about Any (x) Sequence (MIxS). The standards are the Minimum Information about a Single Amplified Genome (MISAG) and the Minimum Information about a Metagenome-Assembled Genome (MIMAG), including, but not limited to, assembly quality, and estimates of genome completeness and contamination. These standards can be used in combination with other GSC checklists, including the Minimum Information about a Genome Sequence (MIGS), Minimum Information about a Metagenomic Sequence (MIMS), and Minimum Information about a Marker Gene Sequence (MIMARKS). Community-wide adoption of MISAG and MIMAG will facilitate more robust comparative genomic analyses of bacterial and archaeal diversity.

Viral to metazoan marine plankton nucleotide sequences from the Tara Oceans expedition

A unique collection of oceanic samples was gathered by the Tara Oceans expeditions (2009–2013), targeting plankton organisms ranging from viruses to metazoans, and providing rich environmental context measurements. Thanks to recent advances in the field of genomics, extensive sequencing has been performed for a deep genomic analysis of this huge collection of samples. A strategy based on different approaches, such as metabarcoding, metagenomics, single-cell genomics and metatranscriptomics, has been chosen for analysis of size-fractionated plankton communities. Here, we provide detailed procedures applied for genomic data generation, from nucleic acids extraction to sequence production, and we describe registries of genomics datasets available at the European Nucleotide Archive (ENA, www.ebi.ac.uk/ena). The association of these metadata to the experimental procedures applied for their generation will help the scientific community to access these data and facilitate their analysis. This paper complements other efforts to provide a full description of experiments and open science resources generated from the Tara Oceans project, further extending their value for the study of the world’s planktonic ecosystems.

Exploring Microdiversity in Novel Kordia sp. (Bacteroidetes) with Proteorhodopsin from the Tropical Indian Ocean via Single Amplified Genomes

Marine Bacteroidetes constitute a very abundant bacterioplankton group in the oceans that plays a key role in recycling particulate organic matter and includes several photoheterotrophic members containing proteorhodopsin. Relatively few marine Bacteroidetes species have been described and, moreover, they correspond to cultured isolates, which in most cases do not represent the actual abundant or ecologically relevant microorganisms in the natural environment. In this study, we explored the microdiversity of 98 Single Amplified Genomes (SAGs) retrieved from the surface waters of the underexplored North Indian Ocean, whose most closely related isolate is Kordia algicida OT-1. Using Multi Locus Sequencing Analysis (MLSA) we found no microdiversity in the tested conserved phylogenetic markers (16S rRNA and 23S rRNA genes), the fast-evolving Internal Transcribed Spacer and the functional markers proteorhodopsin and the beta-subunit of RNA polymerase. Furthermore, we carried out a Fragment Recruitment Analysis (FRA) with marine metagenomes to learn about the distribution and dynamics of this microorganism in different locations, depths and size fractions. This analysis indicated that this taxon belongs to the rare biosphere, showing its highest abundance after upwelling-induced phytoplankton blooms and sinking to the deep ocean with large organic matter particles. This uncultured Kordia lineage likely represents a novel Kordia species (Kordia sp. CFSAG39SUR) that contains the proteorhodopsin gene and has a widespread spatial and vertical distribution. The combination of SAGs and MLSA makes a valuable approach to infer putative ecological roles of uncultured abundant microorganisms.

Unexpected genomic features in widespread intracellular bacteria: evidence for motility of marine chlamydiae

Chlamydiae are obligate intracellular bacteria comprising important human pathogens and symbionts of protists. Molecular evidence indicates a tremendous diversity of chlamydiae particularly in marine environments, yet our current knowledge is based mainly on terrestrial representatives. Here we provide first insights into the biology of marine chlamydiae representing three divergent clades. Our analysis of single-cell amplified genomes revealed hallmarks of the chlamydial lifestyle, supporting the ancient origin of their characteristic developmental cycle and major virulence mechanisms. Surprisingly, these chlamydial genomes encode a complete flagellar apparatus, a previously unreported feature. We show that flagella are an ancient trait that was subject to differential gene loss among extant chlamydiae. Together with a chemotaxis system, these marine chlamydiae are likely motile, with flagella potentially playing a role during host cell infection. This study broadens our view on chlamydial biology and indicates a largely underestimated potential to adapt to different hosts and environments.

SAR202 Genomes from the Dark Ocean Predict Pathways for the Oxidation of Recalcitrant Dissolved Organic Matter

Deep-ocean regions beyond the reach of sunlight contain an estimated 615 Pg of dissolved organic matter (DOM), much of which persists for thousands of years. It is thought that bacteria oxidize DOM until it is too dilute or refractory to support microbial activity. We analyzed five single-amplified genomes (SAGs) from the abundant SAR202 clade of dark-ocean bacterioplankton and found they encode multiple families of paralogous enzymes involved in carbon catabolism, including several families of oxidative enzymes that we hypothesize participate in the degradation of cyclic alkanes. The five partial genomes encoded 152 flavin mononucleotide/F420-dependent monooxygenases (FMNOs), many of which are predicted to be type II Baeyer-Villiger monooxygenases (BVMOs) that catalyze oxygen insertion into semilabile alicyclic alkanes. The large number of oxidative enzymes, as well as other families of enzymes that appear to play complementary roles in catabolic pathways, suggests that SAR202 might catalyze final steps in the biological oxidation of relatively recalcitrant organic compounds to refractory compounds that persist.IMPORTANCE Carbon in the ocean is massively sequestered in a complex mixture of biologically refractory molecules that accumulate as the chemical end member of biological oxidation and diagenetic change. However, few details are known about the biochemical machinery of carbon sequestration in the deep ocean. Reconstruction of the metabolism of a deep-ocean microbial clade, SAR202, led to postulation of new biochemical pathways that may be the penultimate stages of DOM oxidation to refractory forms that persist. These pathways are tied to a proliferation of oxidative enzymes. This research illuminates dark-ocean biochemistry that is broadly consequential for reconstructing the global carbon cycle.

Genomic comparisons of a bacterial lineage that inhabits both marine and terrestrial deep subsurface systems

It is generally accepted that diverse, poorly characterized microorganisms reside deep within Earth’s crust. One such lineage of deep subsurface-dwelling bacteria is an uncultivated member of the Firmicutes phylum that can dominate molecular surveys from both marine and continental rock fracture fluids, sometimes forming the sole member of a single-species microbiome. Here, we reconstructed a genome from basalt-hosted fluids of the deep subseafloor along the eastern Juan de Fuca Ridge flank and used a phylogenomic analysis to show that, despite vast differences in geographic origin and habitat, it forms a monophyletic clade with the terrestrial deep subsurface genome of “Candidatus Desulforudis audaxviator” MP104C. While a limited number of differences were observed between the marine genome of “Candidatus Desulfopertinax cowenii” modA32 and its terrestrial relative that may be of potential adaptive importance, here it is revealed that the two are remarkably similar thermophiles possessing the genetic capacity for motility, sporulation, hydrogenotrophy, chemoorganotrophy, dissimilatory sulfate reduction, and the ability to fix inorganic carbon via the Wood-Ljungdahl pathway for chemoautotrophic growth. Our results provide insights into the genetic repertoire within marine and terrestrial members of a bacterial lineage that is widespread in the global deep subsurface biosphere, and provides a natural means to investigate adaptations specific to these two environments.

De novo synthesis and functional analysis of the phosphatase-encoding gene acI-B of uncultured Actinobacteria from Lake Stechlin (NE Germany)

The National Center for Biotechnology Information [http://www.ncbi.nlm.nih.gov/guide/taxonomy/] database enlists more than 15,500 bacterial species. But this also includes a plethora of uncultured bacterial representations. Owing to their metabolism, they directly influence biogeochemical cycles, which underscores the the important status of bacteria on our planet. To study the function of a gene from an uncultured bacterium, we have undertaken a de novo gene synthesis approach. Actinobacteria of the acI-B subcluster are important but yet uncultured members of the bacterioplankton in temperate lakes of the northern hemisphere such as oligotrophic Lake Stechlin (NE Germany). This lake is relatively poor in phosphate (P) and harbors on average ~1.3 x 10 6 bacterial cells/ml, whereby Actinobacteria of the ac-I lineage can contribute to almost half of the entire bacterial community depending on seasonal variability. Single cell genome analysis of Actinobacterium SCGC AB141-P03, a member of the acI-B tribe in Lake Stechlin has revealed several phosphate-metabolizing genes. The genome of acI-B Actinobacteria indicates potential to degrade polyphosphate compound. To test for this genetic potential, we targeted the exoP-annotated gene potentially encoding polyphosphatase and synthesized it artificially to examine its biochemical role. Heterologous overexpression of the gene in Escherichia coli and protein purification revealed phosphatase activity. Comparative genome analysis suggested that homologs of this gene should be also present in other Actinobacteria of the acI lineages. This strategic retention of specialized genes in their genome provides a metabolic advantage over other members of the aquatic food web in a P-limited ecosystem. [Int Microbiol 2016; 19(1):39-47].

Single-Cell Genome and Group-Specific dsrAB Sequencing Implicate Marine Members of the Class Dehalococcoidia (Phylum Chloroflexi) in Sulfur Cycling

The marine subsurface sediment biosphere is widely inhabited by bacteria affiliated with the class Dehalococcoidia (DEH), phylum Chloroflexi, and yet little is known regarding their metabolisms. In this report, genomic content from a single DEH cell (DEH-C11) with a 16S rRNA gene that was affiliated with a diverse cluster of 16S rRNA gene sequences prevalent in marine sediments was obtained from sediments of Aarhus Bay, Denmark. The distinctive gene content of this cell suggests metabolic characteristics that differ from those of known DEH and Chloroflexi The presence of genes encoding dissimilatory sulfite reductase (Dsr) suggests that DEH could respire oxidized sulfur compounds, although Chloroflexi have never been implicated in this mode of sulfur cycling. Using long-range PCR assays targeting DEH dsr loci, dsrAB genes were amplified and sequenced from various marine sediments. Many of the amplified dsrAB sequences were affiliated with the DEH Dsr clade, which we propose equates to a family-level clade. This provides supporting evidence for the potential for sulfite reduction by diverse DEH species. DEH-C11 also harbored genes encoding reductases for arsenate, dimethyl sulfoxide, and halogenated organics. The reductive dehalogenase homolog (RdhA) forms a monophyletic clade along with RdhA sequences from various DEH-derived contigs retrieved from available metagenomes. Multiple facts indicate that this RdhA may not be a terminal reductase. The presence of other genes indicated that nutrients and energy may be derived from the oxidation of substituted homocyclic and heterocyclic aromatic compounds. Together, these results suggest that marine DEH play a previously unrecognized role in sulfur cycling and reveal the potential for expanded catabolic and respiratory functions among subsurface DEH.

IMPORTANCE

Sediments underlying our oceans are inhabited by microorganisms in cell numbers similar to those estimated to inhabit the oceans. Microorganisms in sediments consist of various diverse and uncharacterized groups that contribute substantially to global biogeochemical cycles. Since most subsurface microorganisms continue to evade cultivation, possibly due to very slow growth, we obtained and analyzed genomic information from a representative of one of the most widespread and abundant, yet uncharacterized bacterial groups of the marine subsurface. We describe several key features that may contribute to their widespread distribution, such as respiratory flexibility and the potential to use oxidized sulfur compounds, which are abundant in marine environments, as electron acceptors. Together, these data provide important information that can be used to assist in designing enrichment strategies or other postgenomic studies, while also improving our understanding of the diversity and distribution of dsrAB genes, which are widely used functional marker genes for sulfur-cycling microbes.

Diversification and niche adaptations of Nitrospina-like bacteria in the polyextreme interfaces of Red Sea brines

Nitrite-oxidizing bacteria (NOB) of the genus Nitrospina have exclusively been found in marine environments. In the brine–seawater interface layer of Atlantis II Deep (Red Sea), Nitrospina-like bacteria constitute up to one-third of the bacterial 16S ribosomal RNA (rRNA) gene sequences. This is much higher compared with that reported in other marine habitats (~10% of all bacteria), and was unexpected because no NOB culture has been observed to grow above 4.0% salinity, presumably due to the low net energy gained from their metabolism that is insufficient for both growth and osmoregulation. Using phylogenetics, single-cell genomics and metagenomic fragment recruitment approaches, we document here that these Nitrospina-like bacteria, designated as Candidatus Nitromaritima RS, are not only highly diverged from the type species Nitrospina gracilis (pairwise genome identity of 69%) but are also ubiquitous in the deeper, highly saline interface layers (up to 11.2% salinity) with temperatures of up to 52 °C. Comparative pan-genome analyses revealed that less than half of the predicted proteome of Ca. Nitromaritima RS is shared with N. gracilis. Interestingly, the capacity for nitrite oxidation is also conserved in both genomes. Although both lack acidic proteomes synonymous with extreme halophiles, the pangenome of Ca. Nitromaritima RS specifically encodes enzymes with osmoregulatory and thermoprotective roles (i.e., ectoine/hydroxyectoine biosynthesis) and of thermodynamic importance (i.e., nitrate and nitrite reductases). Ca. Nitromaritima RS also possesses many hallmark traits of microaerophiles and high-affinity NOB. The abundance of the uncultured Ca. Nitromaritima lineage in marine oxyclines suggests their unrecognized ecological significance in deoxygenated areas of the global ocean.

Single cell genomics indicates horizontal gene transfer and viral infections in a deep subsurface Firmicutes population

A major fraction of Earth's prokaryotic biomass dwells in the deep subsurface, where cellular abundances per volume of sample are lower, metabolism is slower, and generation times are longer than those in surface terrestrial and marine environments. How these conditions impact biotic interactions and evolutionary processes is largely unknown. Here we employed single cell genomics to analyze cell-to-cell genome content variability and signatures of horizontal gene transfer (HGT) and viral infections in five cells of Candidatus Desulforudis audaxviator, which were collected from a 3 km-deep fracture water in the 2.9 Ga-old Witwatersrand Basin of South Africa. Between 0 and 32% of genes recovered from single cells were not present in the original, metagenomic assembly of Desulforudis, which was obtained from a neighboring subsurface fracture. We found a transposable prophage, a retron, multiple clustered regularly interspaced short palindromic repeats (CRISPRs) and restriction-modification systems, and an unusually high frequency of transposases in the analyzed single cell genomes. This indicates that recombination, HGT and viral infections are prevalent evolutionary events in the studied population of microorganisms inhabiting a highly stable deep subsurface environment.

Single-cell genomics-based analysis of virus-host interactions in marine surface bacterioplankton

Viral infections dynamically alter the composition and metabolic potential of marine microbial communities and the evolutionary trajectories of host populations with resulting feedback on biogeochemical cycles. It is quite possible that all microbial populations in the ocean are impacted by viral infections. Our knowledge of virus–host relationships, however, has been limited to a minute fraction of cultivated host groups. Here, we utilized single-cell sequencing to obtain genomic blueprints of viruses inside or attached to individual bacterial and archaeal cells captured in their native environment, circumventing the need for host and virus cultivation. A combination of comparative genomics, metagenomic fragment recruitment, sequence anomalies and irregularities in sequence coverage depth and genome recovery were utilized to detect viruses and to decipher modes of virus–host interactions. Members of all three tailed phage families were identified in 20 out of 58 phylogenetically and geographically diverse single amplified genomes (SAGs) of marine bacteria and archaea. At least four phage–host interactions had the characteristics of late lytic infections, all of which were found in metabolically active cells. One virus had genetic potential for lysogeny. Our findings include first known viruses of Thaumarchaeota, Marinimicrobia, Verrucomicrobia and Gammaproteobacteria clusters SAR86 and SAR92. Viruses were also found in SAGs of Alphaproteobacteria and Bacteroidetes. A high fragment recruitment of viral metagenomic reads confirmed that most of the SAG-associated viruses are abundant in the ocean. Our study demonstrates that single-cell genomics, in conjunction with sequence-based computational tools, enable in situ, cultivation-independent insights into host–virus interactions in complex microbial communities.