Search

Metagenomics and single‐cell genomics have enabled the discovery of relevant uncultured microbes. Recently, single‐virus genomics (SVG), although still in an incipient stage, has opened new avenues in viral ecology by allowing the sequencing of one single virus at a time. The investigation of methodological alternatives and optimization of existing procedures for SVG is paramount to deliver high‐quality genomic data. We report a sequencing dataset of viral single‐amplified genomes (vSAGs) from cultured and uncultured viruses obtained by applying different conditions in each SVG step, from viral preservation and novel whole‐genome amplification (WGA) to sequencing platforms and genome assembly. Sequencing data showed that cryopreservation and mild fixation were compatible with WGA, although fresh samples delivered better genome quality data. The novel TruPrime WGA, based on primase‐polymerase features, and WGA‐X employing a thermostable phi29 polymerase, were proven to be with sufficient sensitivity in SVG. The Oxford Nanopore (ON) sequencing platform did not provide a significant improvement of vSAG assembly compared to Illumina alone. Finally, the SPAdes assembler performed the best. Overall, our results represent a valuable genomic dataset that will help to standardized and advance new tools in viral ecology. This work has been supported by Gordon and Betty Moore Foundation (grant 5334) and Spanish Ministry of Economy and Competitiveness (refs CGL2013‐40564‐R, RTI2018‐094248‐B‐I00 and SAF2013‐49267‐EXP). Work at CRG, BIST and UPF was in part funded by the Spanish Ministry of Economy and Competitiveness, ‘Centro de Excelencia Severo Ochoa 2013‐2017’ and the Spanish Ministry of Economy and Competitiveness, ‘Centro de Excelencia Maria de Maeztu 2016‐2019’.

Members of the SAR11 clade, despite their high abundance, are often poorly represented by metagenome‐assembled genomes. This fact has hampered our knowledge about their ecology and genetic diversity. Here we examined 175 SAR11 genomes, including 47 new single‐amplified genomes. The presence of the first genomes associated with subclade IV suggests that, in the same way as subclade V, they might be outside the proposed Pelagibacterales order. An expanded phylogenomic classification together with patterns of metagenomic recruitment at a global scale have allowed us to define new ecogenomic units of classification (genomospecies), appearing at different, and sometimes restricted, metagenomic data sets. We detected greater microdiversity across the water column at a single location than in samples collected from similar depth across the global ocean, suggesting little influence of biogeography. In addition, pangenome analysis revealed that the flexible genome was essential to shape genomospecies distribution. In one genomospecies preferentially found within the Mediterranean, a set of genes involved in phosphonate utilization was detected. While another, with a more cosmopolitan distribution, was unique in having an aerobic purine degradation pathway. Together, these results provide a glimpse of the enormous genomic diversity within this clade at a finer resolution than the currently defined clades. This work was supported by grant ‘VIREVO’ CGL2016-76273-P [AEI/FEDER, EU] (cofounded with FEDER funds) from the Spanish Ministerio de Economía, Industria y Competitividad to FRV, and grants CGL2013-40564-R and SAF2013-49267-EXP from the Spanish Ministerio de Economía, Industria y Competitividad, grant ACIF/2015/332 from Generalitat Valenciana and grant 5334 from the Betty Moore Foundation to MMG. FRV was also a beneficiary of the 5top100-program of the Ministry for Science and Education of Russia. JHM was supported by a Ph.D. fellowship from the Spanish Ministerio de Economía y Competitividad (BES-2014-067828). MLP was supported by a postdoctoral fellowship from the Spanish Ministerio de Economía, Industria y Competitividad (IJCI-2017-34002).

Absolute abundances of prokaryotes are typically determined by FISH. Due to the lack of a universal conserved gene among all viruses, metagenomic fragment recruitment is commonly used to estimate the relative viral abundance. However, the paucity of absolute virus abundance data hinders our ability to fully understand how viruses drive global microbial populations. The cosmopolitan marine Pelagibacter ubique is host for the highly widespread HTVC010P pelagiphage isolate and the extremely abundant uncultured virus vSAG 37-F6 recently discovered by single-virus genomics. Here we applied droplet digital PCR (ddPCR) to calculate the absolute abundance of these pelagiphage genotypes in the Mediterranean Sea and the Gulf of Maine. Abundances were between 360 and 8,510 virus mL-1 and 1,270–14,400 virus mL-1 for vSAG 37-F6 and HTVC010P, respectively. Illumina PCR-amplicon sequencing corroborated the absence of ddPCR non-specific amplifications for vSAG 37-F6, but showed an overestimation of 6% for HTVC010P from off-targets, genetically unrelated viruses. Absolute abundances of both pelagiphages, two of the most abundance marine viruses, suggest a large viral pelagiphage diversity in marine environments, and show the efficiency and power of ddPCR to disentangle the structure of marine viral communities. Results also highlight the need for a standardized workflow to obtain accurate quantification that allows cross data comparison. This work has been supported by Spanish Ministry of Economy and Competitiveness (Ref. RTI2018-094248-B-I00), Generalitat Valenciana (Ref. ACOM/2015/133 and ACIF/2015/332), and Gordon and Betty Moore Foundation (Grant 5334).

Wastewater treatment plants effluents are considered as hotspots for the dispersion of antibiotic resistance genes (ARGs) into natural ecosystems. The bacterial resistome (ARG collection in a metagenome) analyses have provided clues on antibacterial resistance dynamics. However, viruses and vesicles are frequently ignored. Here, we addressed the bacterial, viral and vesicle resistomes from a representative wastewater effluent in natural conditions and amended with polymyxin, which is used as a last resort antibiotic. Metagenomics showed that the natural prokaryotic resistome was vast (9000 ARG hits/Gb metagenome) and diverse, while viral resistome was two orders of magnitude lower (50 ARG hits/Gb metagenome) suggesting that viruses rarely encoded ARGs. After polymyxin amendment, data showed no ARG enrichment – including to polymyxin – in the microbiome. Remarkably, microbiomes responded to polymyxin with a vast release of putative vesicles (threefold increase compared with the control), which might be used as 'traps' to decrease the antibiotic concentration. Intriguingly, although polymyxin resistance genes (PRGs) were rare in the microbiome (0.018% of total ARG found), in the viral and vesicle fractions, PRGs were more abundant (0.5%–0.8% of total ARG found). Our data suggest that vesicles could have a more active role in the context of transmission of antibiotic resistances. This work has been supported by Ministerio de Ciencia, Innovación y Universidades (ref. RTI2018-094248-B-I00), Generalitat Valenciana (refs. ACOM/2015/133 and ACIF/2015/332), and Gordon and Betty Moore Foundation (grant 5334).

The identification of relevant virus–host pairs that globally account for a large pool of carbon and nutrients in the ocean is paramount to build accurate ecological models. A previous work using single-virus genomics led to the discovery of the uncultured single-virus vSAG 37-F6, originally sorted from the Mediterranean Sea (Blanes Bay Microbial Observatory), that represents one of the most abundant dsDNA viral population in the marine surface virosphere. Here, from same sampling site, we report that a Pelagibacter single-cell contained a viral member of vSAG 37-F6 population, by means of PCR screening of sorted, genome-amplified single cells with vSAG 37-F6-specific primers and whole-genome sequencing. Furthermore, viruses from this population were also found in three other Pelagibacter single cells from the South Pacific and Atlantic oceans. These new uncultured pelagiphages were genetically different from the previously characterized pelagiphage isolates. Data showed that the uncultured vSAG 37-F6 population represents the Pelagibacter phages that inhabit the sunlit ocean better, and contains a vast unrecognized microdiversity. This work has been supported by Spanish Ministry of Economy and Competitiveness (refs. CGL2013-40564-R and SAF2013-49267-EXP), Generalitat Valenciana (refs. ACOM/2015/133 and ACIF/2015/332), and Gordon and Betty Moore Foundation (grant 5334).

Single-cell genomics has unveiled the metabolic potential of dominant microbes inhabiting different environments, including the human body. The lack of genomic information for predominant microbes of the human body, such as bacteriophages, hinders our ability to answer fundamental questions about our viral communities. Here, we applied single-virus genomics (SVGs) to natural human salivary samples in combination with viral metagenomics to gain some insights into the viral community structure of the oral cavity. Saliva samples were processed for viral metagenomics (n = 15) and SVGs (n = 3). A total of 1328 uncultured single viruses were sorted by fluorescence-activated virus sorting followed by whole genome amplification. Sequencing of 24 viral single amplified genomes (vSAGs) showed that half of the vSAGs contained viral hallmark genes. Among those bona fide viruses, the uncultured single virus 92-C13 putatively infecting oral Streptococcus-like species was within the top ≈10 most abundant viruses in the oral virome. Viral gene network and viral metagenomics analyses of 439 oral viruses from cultures, metagenomics, and SVGs revealed that salivary viruses were tentatively structured into ≈200 major viral clusters, corresponding to approximately genus-level groupings. Data showed that none of the publicly available viral isolates, excepting an Actinomyces phage, were significantly abundant in the oral viromes. In addition, none of the obtained viral contigs and vSAGs from this study were present in all viromes. Overall, the data demonstrates that most viral isolates are not naturally abundant in saliva, and furthermore, the predominant viruses in the oral cavity are yet uncharacterized. Results suggest a variable, complex, and interpersonal viral profile. Finally, we demonstrated the power of SVGs in combination with viral metagenomics to unveil the genetic information of the uncultured viruses of the human virome. This work has been supported by the Spanish Ministry of Economy and Competitiveness (refs. CGL2013-40564-R and SAF2013-49267-EXP), the Generalitat Valenciana (refs. ACOM/2015/133 and ACIF/2015/332), and the Gordon and Betty Moore Foundation (grant 5334). We thank the English editor Karen Neller.

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.

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. This research was supported by the NSF grants DEB-841933 and OCE-821374 to RS, OCE-0848703 to CA and by a Maine Technology Institute research infrastructure grant to the Bigelow Laboratory. The Los Alamos National Laboratory researchers were supported in part by the U.S. Department of Energy Joint Genome Institute through the Office of Science of the U.S. Department of Energy under Contract Number DE-AC02-05CH11231 and grants from the U.S. Defense Threat Reduction Agency under contract numbers B104153I and B084531I.

Heterotrophic protists are a highly diverse and biogeochemically significant component of marine ecosystems, yet little is known about their species-specific prey preferences and symbiotic interactions in situ. Here we demonstrate how these previously unresolved questions can be addressed by sequencing the eukaryote and bacterial SSU rRNA genes from individual, uncultured protist cells collected from their natural marine environment and sorted by flow cytometry. We detected Pelagibacter ubique in association with a MAST-4 protist, an actinobacterium in association with a chrysophyte and three bacteroidetes in association with diverse protist groups. The presence of identical phylotypes among the putative prey and the free bacterioplankton in the same sample provides evidence for predator–prey interactions. Our results also suggest a discovery of novel symbionts, distantly related to Rickettsiales and the candidate divisions ZB3 and TG2, associated with Cercozoa and Chrysophyta cells. This study demonstrates the power of single cell sequencing to untangle ecological interactions between uncultured protists and prokaryotes.

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.