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.

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.

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.

Guidelines for establishing a cytometry laboratory.

The purpose of this document is to provide guidance for establishing and maintaining growth and development of flow cytometry shared resource laboratories. While the best practices offered in this manuscript are not intended to be universal or exhaustive, they do outline key goals that should be prioritized to achieve operational excellence and meet the needs of the scientific community. Additionally, this document provides information on available technologies and software relevant to shared resource laboratories. This manuscript builds on the work of Barsky et al. 2016 published in Cytometry Part A and incorporates recent advancements in cytometric technology. A flow cytometer is a specialized piece of technology that require special care and consideration in its housing and operations. As with any scientific equipment, a thorough evaluation of the location, space requirements, auxiliary resources, and support is crucial for successful operation. This comprehensive resource has been written by past and present members of the International Society for Advancement of Cytometry (ISAC) Shared Resource Laboratory (SRL) Emerging Leaders Program https://isac-net.org/general/custom.asp?page=SRL-Emerging-Leaders with extensive expertise in managing flow cytometry SRLs from around the world in different settings including academia and industry. It is intended to assist in establishing a new flow cytometry SRL, re-purposing an existing space into such a facility, or adding a flow cytometer to an individual lab in academia or industry. This resource reviews the available cytometry technologies, the operational requirements, and best practices in SRL staffing and management.

Radiometric approach for the detection of picophytoplankton assemblages across oceanic fronts.

Cell abundances of Prochlorococcus, Synechococcus, and autotrophic picoeukaryotes were estimated in surface waters using principal component analysis (PCA) of hyperspectral and multispectral remote-sensing reflectance data. This involved the development of models that employed multilinear correlations between cell abundances across the Atlantic Ocean and a combination of PCA scores and sea surface temperatures. The models retrieve high Prochlorococcus abundances in the Equatorial Convergence Zone and show their numerical dominance in oceanic gyres, with decreases in Prochlorococcus abundances towards temperate waters where Synechococcus flourishes, and an emergence of picoeukaryotes in temperate waters. Fine-scale in-situ sampling across ocean fronts provided a large dynamic range of measurements for the training dataset, which resulted in the successful detection of fine-scale Synechococcus patches. Satellite implementation of the models showed good performance (R2 > 0.50) when validated against in-situ data from six Atlantic Meridional Transect cruises. The improved relative performance of the hyperspectral models highlights the importance of future high spectral resolution satellite instruments, such as the NASA PACE mission's Ocean Color Instrument, to extend our spatiotemporal knowledge about ecologically relevant phytoplankton assemblages.

Genome to phenome tools: In vivo and in vitro transfection of Crassostrea virginica hemocytes.

The sequencing of the Crassostrea virginica genome has brought back the interest for gene delivery and editing methodologies. Here, we report the expression in oyster hemocytes of two heterologous expression vectors under the CMV promoter delivered with dendrimers. Expression was monitored using confocal microscopy, flow cytometry, and immunofluorescence assay. C. virginica hemocytes were able to express the green fluorescence protein and Crassostrea gigas vascular endothelial growth factor under CMV viral promoter both in vivo and in vitro. These results provide the bases for interrogating the genome and adapting genome editing methodologies.

Imaging flow cytometry for phytoplankton analysis.

This review highlights the concepts and instrumentation of imaging flow cytometry technology and in particular its use for phytoplankton analysis. Imaging flow cytometry, a hybrid technology combining speed and statistical capabilities of flow cytometry with imaging features of microscopy, is rapidly advancing as a cell imaging platform that overcomes many of the limitations of current techniques and contributed significantly to the advancement of phytoplankton analysis in recent years. This review presents the various instrumentation relevant to the field and currently used for assessment of complex phytoplankton communities' composition and abundance, size structure determination, biovolume estimation, detection of harmful algal bloom species, evaluation of viability and metabolic activity and other applications. Also we present our data on viability and metabolic assessment of Aphanizomenon sp. cyanobacteria using Imagestream X Mark II imaging cytometer. Herein, we highlight the immense potential of imaging flow cytometry for microalgal research, but also discuss limitations and future developments.

Accumulation and enhanced cycling of polyphosphate by Sargasso Sea plankton in response to low phosphorus.

Phytoplankton alter their biochemical composition according to nutrient availability, such that their bulk elemental composition varies across oceanic provinces. However, the links between plankton biochemical composition and variation in biogeochemical cycling of nutrients remain largely unknown. In a survey of phytoplankton phosphorus stress in the western North Atlantic, we found that phytoplankton in the phosphorus-depleted subtropical Sargasso Sea were enriched in the biochemical polyphosphate (polyP) compared with nutrient-rich temperate waters, contradicting the canonical oceanographic view of polyP as a luxury phosphorus storage molecule. The enrichment in polyP coincided with enhanced alkaline phosphatase activity and substitution of sulfolipids for phospholipids, which are both indicators of phosphorus stress. Further, polyP appeared to be liberated preferentially over bulk phosphorus from sinking particles in the Sargasso Sea, thereby retaining phosphorus in shallow waters. Thus, polyP cycling may form a feedback loop that attenuates the export of phosphorus when it becomes scarce, contributes bioavailable P for primary production, and supports the export of carbon and nitrogen via sinking particles.

Characterizing the phytoplankton soup: pump and plumbing effects on the particle assemblage in underway optical seawater systems.

Many optical and biogeochemical data sets, crucial for algorithm development and satellite data validation, are collected using underway seawater systems over the course of research cruises. Phytoplankton and particle size distribution (PSD) in the ocean is a key measurement, required in oceanographic research and ocean optics. Using a data set collected in the North Atlantic, spanning different oceanic water types, we outline the differences observed in concurrent samples collected from two different flow-through systems: a permanently plumbed science seawater supply with an impeller pump, and an independent system with shorter, clean tubing runs and a diaphragm pump. We observed an average of 40% decrease in phytoplankton counts, and significant changes to the PSD in 10-45 µm range, when comparing impeller and diaphragm pump systems. Change in PSD seems to be more dependent on the type of the phytoplankton, than the size, with photosynthetic ciliates displaying the largest decreases in cell counts (78%). Comparison of chlorophyll concentrations across the two systems demonstrated lower sensitivity to sampling system type. Observed changes in several measured biogeochemical parameters (associated with phytoplankton size distribution) using the two sampling systems, should be used as a guide towards building best practices when it comes to the deployment of flow-through systems in the field for examining optics and biogeochemistry. Using optical models, we evaluated potential impact of the observed change in measured phytoplankton size spectra onto scattering measurements, resulting in significant differences between modeled optical properties across systems (~40%). Researchers should be aware of the methods used with previously collected data sets, and take into consideration the potentially significant and highly variable ecosystem-dependent biases in designing field studies in the future.

Prevalent genome streamlining and latitudinal divergence of planktonic bacteria in the surface ocean.

Planktonic bacteria dominate surface ocean biomass and influence global biogeochemical processes, but remain poorly characterized owing to difficulties in cultivation. Using large-scale single cell genomics, we obtained insight into the genome content and biogeography of many bacterial lineages inhabiting the surface ocean. We found that, compared with existing cultures, natural bacterioplankton have smaller genomes, fewer gene duplications, and are depleted in guanine and cytosine, noncoding nucleotides, and genes encoding transcription, signal transduction, and noncytoplasmic proteins. These findings provide strong evidence that genome streamlining and oligotrophy are prevalent features among diverse, free-living bacterioplankton, whereas existing laboratory cultures consist primarily of copiotrophs. The apparent ubiquity of metabolic specialization and mixotrophy, as predicted from single cell genomes, also may contribute to the difficulty in bacterioplankton cultivation. Using metagenome fragment recruitment against single cell genomes, we show that the global distribution of surface ocean bacterioplankton correlates with temperature and latitude and is not limited by dispersal at the time scales required for nucleotide substitution to exceed the current operational definition of bacterial species. Single cell genomes with highly similar small subunit rRNA gene sequences exhibited significant genomic and biogeographic variability, highlighting challenges in the interpretation of individual gene surveys and metagenome assemblies in environmental microbiology. Our study demonstrates the utility of single cell genomics for gaining an improved understanding of the composition and dynamics of natural microbial assemblages.

Osmotrophy of dissolved organic compounds by coccolithophore populations: Fixation into particulate organic and inorganic carbon.

Coccolithophores are typically thought of as photoautotrophs, yet a few genera inhabit sub-euphotic environments with insufficient light for photosynthesis, suggesting that other carbon acquisition strategies are likely. Field experiments were performed in the northwest Atlantic (a region with potentially abundant coccolithophores). Phytoplankton populations were incubated with 14C-labeled dissolved organic carbon (DOC) compounds, acetate, mannitol, and glycerol. Coccolithophores were sorted from these populations 24 hours later using flow cytometry, and DOC uptake was measured. DOC uptake rates were as high as 10-15 moles cell-1 day-1, slow relative to photosynthesis rates (10-12 moles cell-1 day-1). Growth rates on the organic compounds were low, suggesting that osmotrophy plays more of a survival strategy in low-light situations. Assimilated DOC was found in both particulate organic carbon and calcite coccoliths (particulate inorganic carbon), suggesting that osmotrophic uptake of DOC into coccolithophore calcite is a small but notable part of the biological carbon pump and alkalinity pump paradigms.

Mixoplankton and mixotrophy: future research priorities.

Phago-mixotrophy, the combination of photoautotrophy and phagotrophy in mixoplankton, organisms that can combine both trophic strategies, have gained increasing attention over the past decade. It is now recognized that a substantial number of protistan plankton species engage in phago-mixotrophy to obtain nutrients for growth and reproduction under a range of environmental conditions. Unfortunately, our current understanding of mixoplankton in aquatic systems significantly lags behind our understanding of zooplankton and phytoplankton, limiting our ability to fully comprehend the role of mixoplankton (and phago-mixotrophy) in the plankton food web and biogeochemical cycling. Here, we put forward five research directions that we believe will lead to major advancement in the field: (i) evolution: understanding mixotrophy in the context of the evolutionary transition from phagotrophy to photoautotrophy; (ii) traits and trade-offs: identifying the key traits and trade-offs constraining mixotrophic metabolisms; (iii) biogeography: large-scale patterns of mixoplankton distribution; (iv) biogeochemistry and trophic transfer: understanding mixoplankton as conduits of nutrients and energy; and (v) in situ methods: improving the identification of in situ mixoplankton and their phago-mixotrophic activity.

Cryptic and abundant marine viruses at the evolutionary origins of Earth's RNA virome.

Whereas DNA viruses are known to be abundant, diverse, and commonly key ecosystem players, RNA viruses are insufficiently studied outside disease settings. In this study, we analyzed ≈28 terabases of Global Ocean RNA sequences to expand Earth's RNA virus catalogs and their taxonomy, investigate their evolutionary origins, and assess their marine biogeography from pole to pole. Using new approaches to optimize discovery and classification, we identified RNA viruses that necessitate substantive revisions of taxonomy (doubling phyla and adding >50% new classes) and evolutionary understanding. "Species"-rank abundance determination revealed that viruses of the new phyla "Taraviricota," a missing link in early RNA virus evolution, and "Arctiviricota" are widespread and dominant in the oceans. These efforts provide foundational knowledge critical to integrating RNA viruses into ecological and epidemiological models.

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.

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.

Single cell ecogenomics reveals mating types of individual cells and ssDNA viral infections in the smallest photosynthetic eukaryotes.

Planktonic photosynthetic organisms of the class Mamiellophyceae include the smallest eukaryotes (less than 2 µm), are globally distributed and form the basis of coastal marine ecosystems. Eight complete fully annotated 13-22 Mb genomes from three genera, Ostreococcus, Bathycoccus and Micromonas, are available from previously isolated clonal cultured strains and provide an ideal resource to explore the scope and challenges of analysing single cell amplified genomes (SAGs) isolated from a natural environment. We assembled data from 12 SAGs sampled during the Tara Oceans expedition to gain biological insights about their in situ ecology, which might be lost by isolation and strain culture. Although the assembled nuclear genomes were incomplete, they were large enough to infer the mating types of four Ostreococcus SAGs. The systematic occurrence of sequences from the mitochondria and chloroplast, representing less than 3% of the total cell's DNA, intimates that SAGs provide suitable substrates for detection of non-target sequences, such as those of virions. Analysis of the non-Mamiellophyceae assemblies, following filtering out cross-contaminations during the sequencing process, revealed two novel 1.6 and 1.8 kb circular DNA viruses, and the presence of specific Bacterial and Oomycete sequences suggests that these organisms might co-occur with the Mamiellales. This article is part of a discussion meeting issue 'Single cell ecology'.

Temporal dynamics of eukaryotic microbial diversity at a coastal Pacific site.

High-throughput sequencing of ocean biomes has revealed vast eukaryotic microbial diversity, a significant proportion of which remains uncharacterized. Here we use a temporal approach to understanding eukaryotic diversity at the Scripps Pier, La Jolla, California, USA, via high-throughput amplicon sequencing of the 18S rRNA gene, the abundances of both Synechococcus and Synechococcus grazers, and traditional oceanographic parameters. We also exploit our ability to track operational taxonomic units (OTUs) temporally to evaluate the ability of 18S sequence-based OTU assignments to meaningfully reflect ecological dynamics. The eukaryotic community is highly dynamic in terms of both species richness and composition, although proportional representation of higher-order taxa remains fairly consistent over time. Synechococcus abundance fluctuates throughout the year. OTUs unique to dates of Synechococcus blooms and crashes or enriched in Synechococcus addition incubation experiments suggest that the prasinophyte Tetraselmis sp. and Gymnodinium-like dinoflagellates are likely Synechococcus grazers under certain conditions, and may play an important role in their population fluctuations.

Single-cell genomics of multiple uncultured stramenopiles reveals underestimated functional diversity across oceans.

Single-celled eukaryotes (protists) are critical players in global biogeochemical cycling of nutrients and energy in the oceans. While their roles as primary producers and grazers are well appreciated, other aspects of their life histories remain obscure due to challenges in culturing and sequencing their natural diversity. Here, we exploit single-cell genomics and metagenomics data from the circumglobal Tara Oceans expedition to analyze the genome content and apparent oceanic distribution of seven prevalent lineages of uncultured heterotrophic stramenopiles. Based on the available data, each sequenced genome or genotype appears to have a specific oceanic distribution, principally correlated with water temperature and depth. The genome content provides hypotheses for specialization in terms of cell motility, food spectra, and trophic stages, including the potential impact on their lifestyles of horizontal gene transfer from prokaryotes. Our results support the idea that prominent heterotrophic marine protists perform diverse functions in ocean ecology.