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.

Procedures for Flow Cytometry-Based Sorting of Unfixed Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) Infected Cells and Other Infectious Agents

In response to the recent COVID‐19 pandemic, many laboratories are involved in research supporting SARS‐CoV‐2 vaccine development and clinical trials. Flow cytometry laboratories will be responsible for a large part of this effort by sorting unfixed antigen‐specific lymphocytes. Therefore, it is critical and timely that we have an understanding of risk assessment and established procedures of infectious cell sorting. Here we present procedures covering the biosafety aspects of sorting unfixed SARS‐CoV‐2‐infected cells and other infectious agents of similar risk level. These procedures follow the ISAC Biosafety Committee guidelines and were recently approved by the National Institutes of Health Institutional Biosafety Committee for sorting SARS‐CoV‐2‐infected cells. © 2020 International Society for Advancement of Cytometry

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.

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.

Novel Impactor and Microsphere-Based Assay Used to Measure Containment of Aerosols Generated in a Flow Cytometer Cell Sorter

Today's state-of-the-art cell sorting flow cytometers are equipped with aerosol containment systems designed to evacuate aerosols from the sort chamber during a sort. This biosafety device is especially important when the sort operator is sorting infectious or potentially infections samples. Hence, it is critical to evaluate the performance for this system in normal operation and in "failure" mode to determine the efficacy of containment. In the past decade, the most popular published method for evaluating containment has been the Glo-Germ bead procedure. These highly fluorescent and multisize particles can easily be detected on a microscope slide and enumerated using a fluorescent microscope. Collecting particles on this slide is accomplished using an Aerotech impactor. This sampler collects potentially escaping aerosols from the sort chamber before enumerating any particles. Although the Glo-Germ procedure has been adopted by many labs, there are several drawbacks with the procedure that have limited its adoption by cell sorter laboratories: The Aerotech impactor is a reusable device that requires rigorous cleaning between measurements. The surface area of the collection slide is large and difficult to scan on a fluorescence microscope. These beads produce a wide variation in sizes resulting in inconsistency in flow rates. Here, we describe a novel and replacement method utilizing a Cyclex-d impactor and Dragon Green beads. This method was compared for sensitivity of detection of escaped aerosols with a published method for aerosol detection which utilizes a UV-APS aerodynamic particle sizer and a UV-excitable dye. One of the advantages of the Cyclex-d system is the narrow-defined field of collection as compared to the standard Glo-Germ bead procedure, this means a smaller sampling area is used in the Cyclex-d impactor as compared to the AeroTech impactor. In addition, the sensitivity of detection was found to be better using the Cyclex-d collection device as compared to the standard Glo-Germ bead procedure. © 2018 International Society for Advancement of Cytometry.

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.

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.

The Role of Ocean Currents in the Temperature Selection of Plankton: Insights from an Individual-Based Model

Biogeography studies that correlate the observed distribution of organisms to environmental variables are typically based on local conditions. However, in cases with substantial translocation, like planktonic organisms carried by ocean currents, selection may happen upstream and local environmental factors may not be representative of those that shaped the local population. Here we use an individual-based model of microbes in the global surface ocean to explore this effect for temperature. We simulate up to 25 million individual cells belonging to up to 50 species with different temperature optima. Microbes are moved around the globe based on a hydrodynamic model, and grow and die based on local temperature. We quantify the role of currents using the “advective temperature differential” metric, which is the optimum temperature of the most abundant species from the model with advection minus that from the model without advection. This differential depends on the location and can be up to 4°C. Poleward-flowing currents, like the Gulf Stream, generally experience cooling and the differential is positive. We apply our results to three global datasets. For observations of optimum growth temperature of phytoplankton, accounting for the effect of currents leads to a slightly better agreement with observations, but there is large variability and the improvement is not statistically significant. For observed Prochlorococcus ecotype ratios and metagenome nucleotide divergence, accounting for advection improves the correlation significantly, especially in areas with relatively strong poleward or equatorward currents.

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.

Genomic and metabolic diversity of Marine Group I Thaumarchaeota in the mesopelagic of two subtropical gyres

Marine Group I (MGI) Thaumarchaeota are one of the most abundant and cosmopolitan chemoautotrophs within the global dark ocean. To date, no representatives of this archaeal group retrieved from the dark ocean have been successfully cultured. We used single cell genomics to investigate the genomic and metabolic diversity of thaumarchaea within the mesopelagic of the subtropical North Pacific and South Atlantic Ocean. Phylogenetic and metagenomic recruitment analysis revealed that MGI single amplified genomes (SAGs) are genetically and biogeographically distinct from existing thaumarchaea cultures obtained from surface waters. Confirming prior studies, we found genes encoding proteins for aerobic ammonia oxidation and the hydrolysis of urea, which may be used for energy production, as well as genes involved in 3-hydroxypropionate/4-hydroxybutyrate and oxidative tricarboxylic acid pathways. A large proportion of protein sequences identified in MGI SAGs were absent in the marine cultures Cenarchaeum symbiosum and Nitrosopumilus maritimus, thus expanding the predicted protein space for this archaeal group. Identifiable genes located on genomic islands with low metagenome recruitment capacity were enriched in cellular defense functions, likely in response to viral infections or grazing. We show that MGI Thaumarchaeota in the dark ocean may have more flexibility in potential energy sources and adaptations to biotic interactions than the existing, surface-ocean cultures.

Comparing effective population sizes of dominant marine alphaproteobacteria lineages

A fundamental question in marine microbial ecology is how microbes adapt to ocean environments. Although numerically dominant populations are typically considered more successful, higher census population sizes (Nc) do not equate directly to a greater capability for adaptation. Instead, effective population size (Ne) determines the fate of deleterious and favourable mutations, and thus is a key parameter for determining the adaptive potential of a population. In the case of the SAR11 and Roseobacter lineages, two abundant heterotrophic bacteria in ocean surface waters with contrasting life history strategies, culture-independent population surveys suggest that SAR11s have greater Nc than Roseobacters. To determine relative Ne, we compared the ratio of nonsynonymous to synonymous substitution rates (ω) of recently diverged lineages of these taxa. Values of ω associated with several of the Roseobacter subclades were lower than for SAR11 subclades, suggesting greater Ne in these cases. Most Roseobacter lineages also had smaller ω values compared with an atypical basal Roseobacter lineage with a large Nc. This finding provides insight into variability in Ne across two important marine bacterial lineages, and provides an evolutionary context for considering how heterotrophic marine bacteria may differ in their ability to adapt to changing ocean habitats.

Single-cell enabled comparative genomics of a deep ocean SAR11 bathytype

Bacterioplankton of the SAR11 clade are the most abundant microorganisms in marine systems, usually representing 25% or more of the total bacterial cells in seawater worldwide. SAR11 is divided into subclades with distinct spatiotemporal distributions (ecotypes), some of which appear to be specific to deep water. Here we examine the genomic basis for deep ocean distribution of one SAR11 bathytype (depth-specific ecotype), subclade Ic. Four single-cell Ic genomes, with estimated completeness of 55%-86%, were isolated from 770 m at station ALOHA and compared with eight SAR11 surface genomes and metagenomic datasets. Subclade Ic genomes dominated metagenomic fragment recruitment below the euphotic zone. They had similar COG distributions, high local synteny and shared a large number (69%) of orthologous clusters with SAR11 surface genomes, yet were distinct at the 16S rRNA gene and amino-acid level, and formed a separate, monophyletic group in phylogenetic trees. Subclade Ic genomes were enriched in genes associated with membrane/cell wall/envelope biosynthesis and showed evidence of unique phage defenses. The majority of subclade Ic-specfic genes were hypothetical, and some were highly abundant in deep ocean metagenomic data, potentially masking mechanisms for niche differentiation. However, the evidence suggests these organisms have a similar metabolism to their surface counterparts, and that subclade Ic adaptations to the deep ocean do not involve large variations in gene content, but rather more subtle differences previously observed deep ocean genomic data, like preferential amino-acid substitutions, larger coding regions among SAR11 clade orthologs, larger intergenic regions and larger estimated average genome size.

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.

Marine viruses, a genetic reservoir revealed by targeted viromics

Metagenomics has opened new windows on investigating viral diversity and functions. Viromic studies typically require large sample volumes and filtration through 0.2 μm pore-size filters, consequently excluding or under-sampling tailed and very large viruses. We have optimized a targeted viromic approach that employs fluorescence-activated sorting and whole genome amplification to produce dsDNA-enriched libraries from discrete viral populations from a 1-ml water sample. Using this approach on an environmental sample from the Patagonian Shelf, we produced three distinct libraries. One of the virus libraries was dominated (79.65% of sequences with known viral homology) by giant viruses from the Mimiviridae and Phycodnaviridae families, while the two other viromes were dominated by smaller phycodnaviruses, cyanophages and other bacteriophages. The estimated genotypic richness and diversity in our sorted viromes, with 52-163 estimated genotypes, was much lower than in previous virome reports. Fragment recruitment of metagenome reads to selected reference viral genomes yields high genome coverage, suggesting little amplification and sequencing bias against some genomic regions. These results underscore the value of our approach as an effective way to target and investigate specific virus groups. In particular, it will help reveal the diversity and abundance of giant viruses in marine ecosystems.

Single-cell genomics shedding light on marine Thaumarchaeota diversification

Previous studies based on analysis of amoA, 16S ribosomal RNA or accA gene sequences have established that marine Thaumarchaeota fall into two phylogenetically distinct groups corresponding to shallow- and deep-water clades, but it is not clear how water depth interacts with other environmental factors, including light, temperature and location, to affect this pattern of diversification. Earlier studies focused on single-gene distributions were not able to link phylogenetic structure to other aspects of functional adaptation. Here, we analyzed the genome content of 46 uncultivated single Thaumarchaeota cells sampled from epi- and mesopelagic waters of subtropical, temperate and polar oceans. Phylogenomic analysis showed that populations diverged by depth, as expected, and that mesopelagic populations from different locations were well mixed. Functional analysis showed that some traits, including putative DNA photolyase and catalase genes that may be related to adaptive mechanisms to reduce light-induced damage, were found exclusively in members of the epipelagic clade. Our analysis of partial genomes has thus confirmed the depth differentiation of Thaumarchaeota populations observed previously, consistent with the distribution of putative mechanisms to reduce light-induced damage in shallow- and deep-water populations.

Capturing single cell genomes of active polysaccharide degraders: an unexpected contribution of Verrucomicrobia

Microbial hydrolysis of polysaccharides is critical to ecosystem functioning and is of great interest in diverse biotechnological applications, such as biofuel production and bioremediation. Here we demonstrate the use of a new, efficient approach to recover genomes of active polysaccharide degraders from natural, complex microbial assemblages, using a combination of fluorescently labeled substrates, fluorescence-activated cell sorting, and single cell genomics. We employed this approach to analyze freshwater and coastal bacterioplankton for degraders of laminarin and xylan, two of the most abundant storage and structural polysaccharides in nature. Our results suggest that a few phylotypes of Verrucomicrobia make a considerable contribution to polysaccharide degradation, although they constituted only a minor fraction of the total microbial community. Genomic sequencing of five cells, representing the most predominant, polysaccharide-active Verrucomicrobia phylotype, revealed significant enrichment in genes encoding a wide spectrum of glycoside hydrolases, sulfatases, peptidases, carbohydrate lyases and esterases, confirming that these organisms were well equipped for the hydrolysis of diverse polysaccharides. Remarkably, this enrichment was on average higher than in the sequenced representatives of Bacteroidetes, which are frequently regarded as highly efficient biopolymer degraders. These findings shed light on the ecological roles of uncultured Verrucomicrobia and suggest specific taxa as promising bioprospecting targets. The employed method offers a powerful tool to rapidly identify and recover discrete genomes of active players in polysaccharide degradation, without the need for cultivation.

Potential for chemolithoautotrophy among ubiquitous bacteria lineages in the dark ocean

Recent studies suggest that unidentified prokaryotes fix inorganic carbon at globally significant rates in the immense dark ocean. Using single-cell sorting and whole-genome amplification of prokaryotes from two subtropical gyres, we obtained genomic DNA from 738 cells representing most cosmopolitan lineages. Multiple cells of Deltaproteobacteria cluster SAR324, Gammaproteobacteria clusters ARCTIC96BD-19 and Agg47, and some Oceanospirillales from the lower mesopelagic contained ribulose-1,5-bisphosphate carboxylase-oxygenase and sulfur oxidation genes. These results corroborated community DNA and RNA profiling from diverse geographic regions. The SAR324 genomes also suggested C(1) metabolism and a particle-associated life-style. Microautoradiography and fluorescence in situ hybridization confirmed bicarbonate uptake and particle association of SAR324 cells. Our study suggests potential chemolithoautotrophy in several uncultured Proteobacteria lineages that are ubiquitous in the dark oxygenated ocean and provides new perspective on carbon cycling in the ocean's largest habitat.

High-throughput single-cell sequencing identifies photoheterotrophs and chemoautotrophs in freshwater bacterioplankton

Recent discoveries suggest that photoheterotrophs (rhodopsin-containing bacteria (RBs) and aerobic anoxygenic phototrophs (AAPs)) and chemoautotrophs may be significant for marine and freshwater ecosystem productivity. However, their abundance and taxonomic identities remain largely unknown. We used a combination of single-cell and metagenomic DNA sequencing to study the predominant photoheterotrophs and chemoautotrophs inhabiting the euphotic zone of temperate, physicochemically diverse freshwater lakes. Multi-locus sequencing of 712 single amplified genomes, generated by fluorescence-activated cell sorting and whole genome multiple displacement amplification, showed that most of the cosmopolitan freshwater clusters contain photoheterotrophs. These comprised at least 10-23% of bacterioplankton, and RBs were the dominant fraction. Our data demonstrate that Actinobacteria, including clusters acI, Luna and acSTL, are the predominant freshwater RBs. We significantly broaden the known taxonomic range of freshwater RBs, to include Alpha-, Beta-, Gamma- and Deltaproteobacteria, Verrucomicrobia and Sphingobacteria. By sequencing single cells, we found evidence for inter-phyla horizontal gene transfer and recombination of rhodopsin genes and identified specific taxonomic groups involved in these evolutionary processes. Our data suggest that members of the ubiquitous betaproteobacteria Polynucleobacter spp. are the dominant AAPs in temperate freshwater lakes. Furthermore, the RuBisCO (ribulose 1,5-bisphosphate carboxylase/oxygenase) gene was found in several single cells of Betaproteobacteria, Bacteroidetes and Gammaproteobacteria, suggesting that chemoautotrophs may be more prevalent among aerobic bacterioplankton than previously thought. This study demonstrates the power of single-cell DNA sequencing addressing previously unresolved questions about the metabolic potential and evolutionary histories of uncultured microorganisms, which dominate most natural environments.