Expanding the marine virosphere using metagenomics

PhagePhisher: a pipeline for the discovery of covert viral sequences in complex genomic datasets

Obtaining meaningful viral information from large sequencing datasets presents unique challenges distinct from prokaryotic and eukaryotic sequencing efforts. The difficulties surrounding this issue can be ascribed in part to the genomic plasticity of viruses themselves as well as the scarcity of existing information in genomic databases. The open-source software PhagePhisher (http://www.putonti-lab.com/phagephisher) has been designed as a simple pipeline to extract relevant information from complex and mixed datasets, and will improve the examination of bacteriophages, viruses, and virally related sequences, in a range of environments. Key aspects of the software include speed and ease of use; PhagePhisher can be used with limited operator knowledge of bioinformatics on a standard workstation. As a proof-of-concept, PhagePhisher was successfully implemented with bacteria–virus mixed samples of varying complexity. Furthermore, viral signals within microbial metagenomic datasets were easily and quickly identified by PhagePhisher, including those from prophages as well as lysogenic phages, an important and often neglected aspect of examining phage populations in the environment. PhagePhisher resolves viral-related sequences which may be obscured by or imbedded in bacterial genomes.

Linking Virus Genomes with Host Taxonomy

Environmental genomics can describe all forms of organisms—cellular and viral—present in a community. The analysis of such eco-systems biology data relies heavily on reference databases, e.g., taxonomy or gene function databases. Reference databases of symbiosis sensu lato, although essential for the analysis of organism interaction networks, are lacking. By mining existing databases and literature, we here provide a comprehensive and manually curated database of taxonomic links between viruses and their cellular hosts.

Computational approaches to predict bacteriophage-host relationships

Metagenomics has changed the face of virus discovery by enabling the accurate identification of viral genome sequences without requiring isolation of the viruses. As a result, metagenomic virus discovery leaves the first and most fundamental question about any novel virus unanswered: What host does the virus infect? The diversity of the global virosphere and the volumes of data obtained in metagenomic sequencing projects demand computational tools for virus–host prediction. We focus on bacteriophages (phages, viruses that infect bacteria), the most abundant and diverse group of viruses found in environmental metagenomes. By analyzing 820 phages with annotated hosts, we review and assess the predictive power of in silico phage–host signals. Sequence homology approaches are the most effective at identifying known phage–host pairs. Compositional and abundance-based methods contain significant signal for phage–host classification, providing opportunities for analyzing the unknowns in viral metagenomes. Together, these computational approaches further our knowledge of the interactions between phages and their hosts. Importantly, we find that all reviewed signals significantly link phages to their hosts, illustrating how current knowledge and insights about the interaction mechanisms and ecology of coevolving phages and bacteria can be exploited to predict phage–host relationships, with potential relevance for medical and industrial applications.

New viruses infecting bacteria are increasingly being discovered in many environments through sequence-based explorations. To understand their role in microbial ecosystems, computational tools are indispensable to prioritize and guide experimental efforts. This review assesses and discusses a range of bioinformatic approaches to predict bacteriophage–host relationships when all that is known is their genome sequence.

Complete genome sequence of bacteriophage P8625, the first lytic phage that infects Verrucomicrobia

Bacteriophage P8625 is a lytic bacteriophage that infects the verrucomicrobial strain IMCC8625, a marine bacterium affiliated with Verrucomicrobia subdivision 4. Both the bacteriophage and the host bacterial strain were isolated from surface seawater samples collected off the east coast of Korea. The phage particle has an icosahedral capsid with a diameter of ~47 nm and a long tail of ~75 nm in length, showing the distinctive morphology of the Siphoviridae family. The complete genome sequence of phage P8625 is 32,894 bp long with 51.0 % G + C content. This is the first report of the complete genome sequence of a lytic phage that infects the Verrucomicrobia, for which the name "verrucophage" is proposed.

First Insights into the Viral Communities of the Deep-sea Anoxic Brines of the Red Sea

The deep-sea brines of the Red Sea include some of the most extreme and unique environments on Earth. They combine high salinities with increases in temperature, heavy metals, hydrostatic pressure, and anoxic conditions, creating unique settings for thriving populations of novel extremophiles. Despite a recent increase of studies focusing on these unusual biotopes, their viral communities remain unexplored. The current survey explores four metagenomic datasets obtained from different brine–seawater interface samples, focusing specifically on the diversity of their viral communities. Data analysis confirmed that the particle-attached viral communities present in the brine–seawater interfaces were diverse and generally dominated by Caudovirales, yet appearing distinct from sample to sample. With a level of caution, we report the unexpected finding of Phycodnaviridae, which infects algae and plants, and trace amounts of insect-infecting Iridoviridae. Results from Kebrit Deep revealed stratification in the viral communities present in the interface: the upper-interface was enriched with viruses associated with typical marine bacteria, while the lower-interface was enriched with haloviruses and halophages. These results provide first insights into the unexplored viral communities present in deep-sea brines of the Red Sea, representing one of the first steps for ongoing and future sampling efforts and studies.

Shedding new light on viral photosynthesis

Viruses infecting the environmentally important marine cyanobacteria Prochlorococcus and Synechococcus encode 'auxiliary metabolic genes' (AMGs) involved in the light and dark reactions of photosynthesis. Here, we discuss progress on the inventory of such AMGs in the ever-increasing number of viral genome sequences as well as in metagenomic datasets. We contextualise these gene acquisitions with reference to a hypothesised fitness gain to the phage. We also report new evidence with regard to the sequence and predicted structural properties of viral petE genes encoding the soluble electron carrier plastocyanin. Viral copies of PetE exhibit extensive modifications to the N-terminal signal peptide and possess several novel residues in a region responsible for interaction with redox partners. We also highlight potential knowledge gaps in this field and discuss future opportunities to discover novel phage-host interactions involved in the photosynthetic process.

Challenges of metagenomics and single-cell genomics approaches for exploring cyanobacterial diversity

Cyanobacteria have played a crucial role in the history of early earth and continue to be instrumental in shaping our planet, yet applications of cutting edge technology have not yet been widely used to explore cyanobacterial diversity. To provide adequate background, we briefly review current sequencing technologies and their innovative uses in genomics and metagenomics. Next, we focus on current cell capture technologies and the challenges of using them with cyanobacteria. We illustrate the utility in coupling breakthroughs in DNA amplification with cell capture platforms, with an example of microfluidic isolation and subsequent targeted amplicon sequencing from individual terrestrial thermophilic cyanobacteria. Single cells of thermophilic, unicellular Synechococcus sp. JA-2-3-B'a(2-13) (Syn OS-B') were sorted in a microfluidic device, lysed, and subjected to whole genome amplification by multiple displacement amplification. We amplified regions from specific CRISPR spacer arrays, which are known to be highly diverse, contain semi-palindromic repeats which form secondary structure, and can be difficult to amplify. Cell capture, lysis, and genome amplification on a microfluidic device have been optimized, setting a stage for further investigations of individual cyanobacterial cells isolated directly from natural populations.

Evolution of viruses and cells: do we need a fourth domain of life to explain the origin of eukaryotes?

The recent discovery of diverse very large viruses, such as the mimivirus, has fostered a profusion of hypotheses positing that these viruses define a new domain of life together with the three cellular ones (Archaea, Bacteria and Eucarya). It has also been speculated that they have played a key role in the origin of eukaryotes as donors of important genes or even as the structures at the origin of the nucleus. Thanks to the increasing availability of genome sequences for these giant viruses, those hypotheses are amenable to testing via comparative genomic and phylogenetic analyses. This task is made very difficult by the high evolutionary rate of viruses, which induces phylogenetic artefacts, such as long branch attraction, when inadequate methods are applied. It can be demonstrated that phylogenetic trees supporting viruses as a fourth domain of life are artefactual. In most cases, the presence of homologues of cellular genes in viruses is best explained by recurrent horizontal gene transfer from cellular hosts to their infecting viruses and not the opposite. Today, there is no solid evidence for the existence of a viral domain of life or for a significant implication of viruses in the origin of the cellular domains.

The rapidly expanding universe of giant viruses: Mimivirus, Pandoravirus, Pithovirus and Mollivirus

More than a century ago, the term 'virus' was introduced to describe infectious agents that are invisible by light microscopy and capable of passing through sterilizing filters. In addition to their extremely small size, most viruses have minimal genomes and gene contents, and rely almost entirely on host cell-encoded functions to multiply. Unexpectedly, four different families of eukaryotic 'giant viruses' have been discovered over the past 10 years with genome sizes, gene contents and particle dimensions overlapping with that of cellular microbes. Their ongoing analyses are challenging accepted ideas about the diversity, evolution and origin of DNA viruses.

Ocean plankton. Patterns and ecological drivers of ocean viral communities

Viruses influence ecosystems by modulating microbial population size, diversity, metabolic outputs, and gene flow. Here, we use quantitative double-stranded DNA (dsDNA) viral-fraction metagenomes (viromes) and whole viral community morphological data sets from 43 Tara Oceans expedition samples to assess viral community patterns and structure in the upper ocean. Protein cluster cataloging defined pelagic upper-ocean viral community pan and core gene sets and suggested that this sequence space is well-sampled. Analyses of viral protein clusters, populations, and morphology revealed biogeographic patterns whereby viral communities were passively transported on oceanic currents and locally structured by environmental conditions that affect host community structure. Together, these investigations establish a global ocean dsDNA viromic data set with analyses supporting the seed-bank hypothesis to explain how oceanic viral communities maintain high local diversity.

The diversity and host interactions of Propionibacterium acnes bacteriophages on human skin

The viral population, including bacteriophages, is an important component of the human microbiota, yet is poorly understood. We aim to determine whether bacteriophages modulate the composition of the bacterial populations, thus potentially playing a role in health or disease. We investigated the diversity and host interactions of the bacteriophages of Propionibacterium acnes, a major human skin commensal implicated in acne pathogenesis. By sequencing 48 P. acnes phages isolated from acne patients and healthy individuals and by analyzing the P. acnes phage populations in healthy skin metagenomes, we revealed that P. acnes phage populations in the skin microbial community are often dominated by one strain. We also found phage strains shared among both related and unrelated individuals, suggesting that a pool of common phages exists in the human population and that transmission of phages may occur between individuals. To better understand the bacterium-phage interactions in the skin microbiota, we determined the outcomes of 74 genetically defined Propionibacterium strains challenged by 15 sequenced phages. Depending on the Propionibacterium lineage, phage infection can result in lysis, pseudolysogeny, or resistance. In type II P. acnes strains, we found that encoding matching clustered regularly interspaced short palindromic repeat spacers is insufficient to confer phage resistance. Overall, our findings suggest that the prey-predator relationship between bacteria and phages may have a role in modulating the composition of the microbiota. Our study also suggests that the microbiome structure of an individual may be an important factor in the design of phage-based therapy.

Biogeography of Viruses in the Sea

Viral ecology is a rapidly progressing area of research, as molecular methods have improved significantly for targeted research on specific populations and whole communities. To interpret and synthesize global viral diversity and distribution, it is feasible to assess whether macroecology concepts can apply to marine viruses. We review how viral and host life history and physical properties can influence viral distribution in light of biogeography and metacommunity ecology paradigms. We highlight analytical approaches that can be applied to emerging global data sets and meta-analyses to identify individual taxa with global influence and drivers of emergent properties that influence microbial community structure by drawing on examples across the spectrum of viral taxa, from RNA to ssDNA and dsDNA viruses.

Viral dark matter and virus-host interactions resolved from publicly available microbial genomes

The ecological importance of viruses is now widely recognized, yet our limited knowledge of viral sequence space and virus–host interactions precludes accurate prediction of their roles and impacts. In this study, we mined publicly available bacterial and archaeal genomic data sets to identify 12,498 high-confidence viral genomes linked to their microbial hosts. These data augment public data sets 10-fold, provide first viral sequences for 13 new bacterial phyla including ecologically abundant phyla, and help taxonomically identify 7–38% of ‘unknown’ sequence space in viromes. Genome- and network-based classification was largely consistent with accepted viral taxonomy and suggested that (i) 264 new viral genera were identified (doubling known genera) and (ii) cross-taxon genomic recombination is limited. Further analyses provided empirical data on extrachromosomal prophages and coinfection prevalences, as well as evaluation of in silico virus–host linkage predictions. Together these findings illustrate the value of mining viral signal from microbial genomes.

DOI:http://dx.doi.org/10.7554/eLife.08490.001

Viruses are infectious particles that can only multiply inside the cells of microbes and other organisms. Little is known about the genetic differences between virus particles (so-called ‘genetic diversity’), especially compared to what we know about the diversity of bacteria, archaea, and other single-celled microbes. This lack of knowledge hampers our understanding of the role viruses play in the evolution of microbial communities and their associated ecosystems.

Studying the genetics of the viruses in these communities is challenging. There is no single ‘marker’ gene that can be used to identify all viruses in environmental samples. Also, many of the fragments of viral genomes that have been identified have not yet been linked to their host microbes. Many viruses integrate their genome into the DNA of their host cell, and there are computational tools available that exploit this ability to identify viruses and link them to their host. However, other viruses can live and multiply inside cells without integrating their genome into the host's DNA.

Earlier in 2015, researchers developed a new computational tool called VirSorter that can predict virus genome sequences within the DNA extracted from microbes. VirSorter identifies viral genome sequences based on the presence of ‘hallmark’ genes that encode for components found in many virus particles, together with a reference database of genomes from many viruses.

Now, Roux et al.—including some of the researchers from the earlier work—use VirSorter to predict viral DNA from publicly available bacteria and archaea genome data. The study identifies over 12,000 viral genomes and links them to their microbial hosts. These data increase the number of viral genome sequences that are publically available by a factor of ten and identify the first viruses associated with 13 new types of bacteria, which include species that are abundant in particular environments.

It is possible for several different viruses to infect a single cell at the same time. Some viruses are known to be able to exchange DNA, and if this happens frequently in other viruses, it could have a big impact on how viruses evolve. Roux et al.'s findings suggest that although it is common for several different viruses to infect the same cell, it is relatively rare for these viruses to exchange genetic material.

Roux et al.'s findings demonstrate the value of searching publicly available microbial genome data for fragments of viral genomes. These new viral genomes will serve as a useful resource for researchers as they explore the communities of viruses and microbes in natural environments, the human body and in industrial processes.

DOI:http://dx.doi.org/10.7554/eLife.08490.002

Analysis of virus genomes from glacial environments reveals novel virus groups with unusual host interactions

Microbial communities in glacial ecosystems are diverse, active, and subjected to strong viral pressures and infection rates. In this study we analyse putative virus genomes assembled from three dsDNA viromes from cryoconite hole ecosystems of Svalbard and the Greenland Ice Sheet to assess the potential hosts and functional role viruses play in these habitats. We assembled 208 million reads from the virus-size fraction and developed a procedure to select genuine virus scaffolds from cellular contamination. Our curated virus library contained 546 scaffolds up to 230 Kb in length, 54 of which were circular virus consensus genomes. Analysis of virus marker genes revealed a wide range of viruses had been assembled, including bacteriophages, cyanophages, nucleocytoplasmic large DNA viruses and a virophage, with putative hosts identified as Cyanobacteria, Alphaproteobacteria, Gammaproteobacteria, Actinobacteria, Firmicutes, eukaryotic algae and amoebae. Whole genome comparisons revealed the majority of circular genome scaffolds (CGS) formed 12 novel groups, two of which contained multiple phage members with plasmid-like properties, including a group of phage-plasmids possessing plasmid-like partition genes and toxin-antitoxin addiction modules to ensure their replication and a satellite phage-plasmid group. Surprisingly we also assembled a phage that not only encoded plasmid partition genes, but a clustered regularly interspaced short palindromic repeat (CRISPR)/Cas adaptive bacterial immune system. One of the spacers was an exact match for another phage in our virome, indicating that in a novel use of the system, the lysogen was potentially capable of conferring immunity on its bacterial host against other phage. Together these results suggest that highly novel and diverse groups of viruses are present in glacial environments, some of which utilize very unusual life strategies and genes to control their replication and maintain a long-term relationship with their hosts.

Bayesian mixture analysis for metagenomic community profiling

Motivation: Deep sequencing of clinical samples is now an established tool for the detection of infectious pathogens, with direct medical applications. The large amount of data generated produces an opportunity to detect species even at very low levels, provided that computational tools can effectively profile the relevant metagenomic communities. Data interpretation is complicated by the fact that short sequencing reads can match multiple organisms and by the lack of completeness of existing databases, in particular for viral pathogens. Here we present metaMix, a Bayesian mixture model framework for resolving complex metagenomic mixtures. We show that the use of parallel Monte Carlo Markov chains for the exploration of the species space enables the identification of the set of species most likely to contribute to the mixture.

Results: We demonstrate the greater accuracy of metaMix compared with relevant methods, particularly for profiling complex communities consisting of several related species. We designed metaMix specifically for the analysis of deep transcriptome sequencing datasets, with a focus on viral pathogen detection; however, the principles are generally applicable to all types of metagenomic mixtures.

Availability and implementation: metaMix is implemented as a user friendly R package, freely available on CRAN: http://cran.r-project.org/web/packages/metaMix

Contact:sofia.morfopoulou.10@ucl.ac.uk

Supplementary information:Supplementary data are available at Bionformatics online.

Multidimensional metrics for estimating phage abundance, distribution, gene density, and sequence coverage in metagenomes

Phages are the most abundant biological entities on Earth and play major ecological roles, yet the current sequenced phage genomes do not adequately represent their diversity, and little is known about the abundance and distribution of these sequenced genomes in nature. Although the study of phage ecology has benefited tremendously from the emergence of metagenomic sequencing, a systematic survey of phage genes and genomes in various ecosystems is still lacking, and fundamental questions about phage biology, lifestyle, and ecology remain unanswered. To address these questions and improve comparative analysis of phages in different metagenomes, we screened a core set of publicly available metagenomic samples for sequences related to completely sequenced phages using the web tool, Phage Eco-Locator. We then adopted and deployed an array of mathematical and statistical metrics for a multidimensional estimation of the abundance and distribution of phage genes and genomes in various ecosystems. Experiments using those metrics individually showed their usefulness in emphasizing the pervasive, yet uneven, distribution of known phage sequences in environmental metagenomes. Using these metrics in combination allowed us to resolve phage genomes into clusters that correlated with their genotypes and taxonomic classes as well as their ecological properties. We propose adding this set of metrics to current metaviromic analysis pipelines, where they can provide insight regarding phage mosaicism, habitat specificity, and evolution.

Combining genomic sequencing methods to explore viral diversity and reveal potential virus-host interactions

Viral diversity and virus-host interactions in oxygen-starved regions of the ocean, also known as oxygen minimum zones (OMZs), remain relatively unexplored. Microbial community metabolism in OMZs alters nutrient and energy flow through marine food webs, resulting in biological nitrogen loss and greenhouse gas production. Thus, viruses infecting OMZ microbes have the potential to modulate community metabolism with resulting feedback on ecosystem function. Here, we describe viral communities inhabiting oxic surface (10 m) and oxygen-starved basin (200 m) waters of Saanich Inlet, a seasonally anoxic fjord on the coast of Vancouver Island, British Columbia using viral metagenomics and complete viral fosmid sequencing on samples collected between April 2007 and April 2010. Of 6459 open reading frames (ORFs) predicted across all 34 viral fosmids, 77.6% (n = 5010) had no homology to reference viral genomes. These fosmids recruited a higher proportion of viral metagenomic sequences from Saanich Inlet than from nearby northeastern subarctic Pacific Ocean (Line P) waters, indicating differences in the viral communities between coastal and open ocean locations. While functional annotations of fosmid ORFs were limited, recruitment to NCBI's non-redundant “nr” database and publicly available single-cell genomes identified putative viruses infecting marine thaumarchaeal and SUP05 proteobacteria to provide potential host linkages with relevance to coupled biogeochemical cycling processes in OMZ waters. Taken together, these results highlight the power of coupled analyses of multiple sequence data types, such as viral metagenomic and fosmid sequence data with prokaryotic single cell genomes, to chart viral diversity, elucidate genomic and ecological contexts for previously unclassifiable viral sequences, and identify novel host interactions in natural and engineered ecosystems.

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.

Polar freshwater cyanophage S-EIV1 represents a new widespread evolutionary lineage of phages

Cyanobacteria are often the dominant phototrophs in polar freshwater communities; yet, the phages that infect them remain unknown. Here, we present a genomic and morphological characterization of cyanophage S-EIV1 that was isolated from freshwaters on Ellesmere Island (Nunavut, High Arctic Canada), and which infects the polar Synechococcus sp., strain PCCC-A2c. S-EIV1 represents a newly discovered evolutionary lineage of bacteriophages whose representatives are widespread in aquatic systems. Among the 130 predicted open reading frames (ORFs) there is no recognizable similarity to genes that encode structural proteins other than the large terminase subunit and a distant viral morphogenesis protein, indicating that the genes encoding the structural proteins of S-EIV1 are distinct from other viruses. As well, only 19 predicted coding sequences on the 79 178 bp circularly permuted genome have homology with genes encoding proteins of known function. Although S-EIV1 is divergent from other sequenced phage isolates, it shares synteny with phage genes captured on a fosmid from the deep-chlorophyll maximum in the Mediterranean Sea, as well as with an incision element in the genome of Anabaena variabilis (ATCC 29413). Sequence recruitment of metagenomic data indicates that S-EIV1-like viruses are cosmopolitan and abundant in a wide range of aquatic systems, suggesting they have an important ecological role.

Diversity of Viruses Infecting the Green Microalga Ostreococcus lucimarinus

The functional diversity of eukaryotic viruses infecting a single host strain from seawater samples originating from distant marine locations is unknown. To estimate this diversity, we used lysis plaque assays to detect viruses that infect the widespread species Ostreococcus lucimarinus, which is found in coastal and mesotrophic systems, and O. tauri, which was isolated from coastal and lagoon sites from the northwest Mediterranean Sea. Detection of viral lytic activities against O. tauri was not observed using seawater from most sites, except those close to the area where the host strain was isolated. In contrast, the more cosmopolitan O. lucimarinus species recovered viruses from locations in the Atlantic and Pacific Oceans and the Mediterranean Sea. Six new O. lucimarinus viruses (OlVs) then were characterized and their genomes sequenced. Two subgroups of OlVs were distinguished based on their genetic distances and on the inversion of a central 32-kb-long DNA fragment, but overall their genomes displayed a high level of synteny. The two groups did not correspond to proximity of isolation sites, and the phylogenetic distance between these subgroups was higher than the distances observed among viruses infecting O. tauri. Our study demonstrates that viruses originating from very distant sites are able to infect the same algal host strain and can be more diverse than those infecting different species of the same genus. Finally, distinctive features and evolutionary distances between these different viral subgroups does not appear to be linked to biogeography of the viral isolates.

IMPORTANCE

Marine eukaryotic phytoplankton virus diversity has yet to be addressed, and more specifically, it is unclear whether diversity is connected to geographical distance and whether differential infection and lysis patterns exist among such viruses that infect the same host strain. Here, we assessed the genetic distance of geographically segregated viruses that infect the ubiquitous green microalga Ostreococcus. This study provides the first glimpse into the diversity of predicted gene functions in Ostreococcus viruses originating from distant sites and provides new insights into potential host distributions and restrictions in the world oceans.

Replicating phages in the epidermal mucosa of the eel (Anguilla anguilla)

In this work, we used the eel (Anguilla anguilla) as an animal model to test the hypothesis of Barr et al. (2013a,b) about the putative role of the epidermal mucosa as a phage enrichment layer. To this end, we analyzed the microbial content of the skin mucus of wild and farmed eels by using a metagenomic approach. We found a great abundance of replicating phage genomes (concatemers) in all the samples. They were assembled in four complete genomes of three Myovirus and one Podovirus. We also found evidences that ΦKZ and Podovirus phages could be part of the resident microbiota associated to the eel mucosal surface and persist on them over the time. Moreover, the viral abundance estimated by epiflorescent counts and by metagenomic recruitment from eel mucosa was higher than that of the surrounding water. Taken together, our results support the hypothesis that claims a possible role of phages in the animal mucus as agents controlling bacterial populations, including pathogenic species, providing a kind of innate immunity.

Rising to the challenge: accelerated pace of discovery transforms marine virology

Marine viruses have important roles in microbial mortality, gene transfer, metabolic reprogramming and biogeochemical cycling. In this Review, we discuss recent technological advances in marine virology including the use of near-quantitative, reproducible metagenomics for large-scale investigation of viral communities and the emergence of gene-based viral ecology. We also describe the reprogramming of microbially driven processes by viral metabolic genes, the identification of novel viruses using cultivation-dependent and cultivation-independent tools, and the potential for modelling studies to provide a framework for studying virus-host interactions. These transformative advances have set a rapid pace in exploring and predicting how marine viruses manipulate and respond to their environment.

Effects of actin-like proteins encoded by two Bacillus pumilus phages on unstable lysogeny, revealed by genomic analysis

We characterized two newly isolated myoviruses, Bp8p-C and Bp8p-T, infecting the ginger rhizome rot disease pathogen Bacillus pumilus GR8. The plaque of Bp8p-T exhibited a clear center with a turbid rim, suggesting that Bp8p-T could transform into latent phage. Lysogeny assays showed that both the two phages could form latent states, while Bp8p-T could form latent phage at a higher frequency and stability than Bp8p-C. The genomes of Bp8p-C and Bp8p-T were 151,417 and 151,419 bp, respectively; both encoded 212 putative proteins, and only differed by three nucleotides. Moreover, owing to this difference, Bp8p-C encoded a truncated, putative actin-like plasmid segregation protein Gp27-C. Functional analysis of protein Gp27 showed that Gp27-T encoded by Bp8p-T exhibited higher ATPase activity and assembly ability than Gp27-C. The results indicate that the difference in Gp27 affected the phage lysogenic ability. Structural proteome analysis of Bp8p-C virion resulted in the identification of 14 structural proteins, among which a pectin lyase-like protein, a putative poly-gamma-glutamate hydrolase, and three proteins with unknown function, were firstly identified as components of the phage virion. Both phages exhibited specific lytic ability to the host strain GR8. Bp8p-C showed better control effect on the pathogen in ginger rhizome slices than Bp8p-T, suggesting that Bp8p-C has a potential application in bio-control of ginger rhizome rot disease.

Emerging methods to study bacteriophage infection at the single-cell level

Bacteria and their viruses (phages) are abundant across diverse ecosystems and their interactions influence global biogeochemical cycles and incidence of disease. Problematically, both classical and metagenomic methods insufficiently assess the host specificity of phages and phage–host infection dynamics in nature. Here we review emerging methods to study phage–host interaction and infection dynamics with a focus on those that offer resolution at the single-cell level. These methods leverage ever-increasing sequence data to identify virus signals from single-cell amplified genome datasets or to produce primers/probes to target particular phage–bacteria pairs (digital PCR and phageFISH), even in complex communities. All three methods enable study of phage infection of uncultured bacteria from environmental samples, while the latter also discriminates between phage–host interaction outcomes (e.g., lytic, chronic, lysogenic) in model systems. Together these techniques enable quantitative, spatiotemporal studies of phage–bacteria interactions from environmental samples of any ecosystem, which will help elucidate and predict the ecological and evolutionary impacts of specific phage–host pairings in nature.

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

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

Deciphering diversity and ecological function from marine metagenomes

Metagenomic sequencing now represents a common, powerful approach for investigating diversity and functional relationships in marine ecosystems. High-throughput datasets generated from random fragments of environmental DNA can provide a less biased view of organismal abundance (versus PCR-based amplicon sequencing) and enable novel exploration of microbial genomes by recovering genome assemblies from uncultured species, identifying ecological functions, and reconstructing metabolic pathways. This review highlights the current state of knowledge in marine metagenomics, focusing on biological insights gained from recent environmental studies and detailing commonly employed methods for data collection and analysis.

The interaction of human microbial pathogens, particulate material and nutrients in estuarine environments and their impacts on recreational and shellfish waters

Anthropogenic activities have increased the load of faecal bacteria, pathogenic viruses and nutrients in rivers, estuaries and coastal areas through point and diffuse sources such as sewage discharges and agricultural runoff. These areas are used by humans for both commercial and recreational activities and are therefore protected by a range of European Directives. If water quality declines in these zones, significant economic losses can occur. Identifying the sources of pollution, however, is notoriously difficult due to the ephemeral nature of discharges, their diffuse source, and uncertainties associated with transport and transformation of the pollutants through the freshwater-marine interface. Further, significant interaction between nutrients, microorganisms and particulates can occur in the water column making prediction of the fate and potential infectivity of human pathogenic organisms difficult to ascertain. This interaction is most prevalent in estuarine environments due to the formation of flocs (suspended sediment) at the marine-freshwater interface. A range of physical, chemical and biological processes can induce the co-flocculation of microorganisms, organic matter and mineral particles resulting in pathogenic organisms becoming potentially protected from a range of biotic (e.g. predation) and abiotic stresses (e.g. UV, salinity). These flocs contain and retain macro- and micro- nutrients allowing the potential survival, growth and transfer of pathogenic organisms to commercially sensitive areas (e.g. beaches, shellfish harvesting waters). The flocs can either be transported directly to the coastal environment or can become deposited in the estuary forming cohesive sediments where pathogens can survive for long periods. Especially in response to storms, these sediments can be subsequently remobilised releasing pulses of potential pathogenic organisms back into the water column leading to contamination of marine waters long after the initial contamination event occurred. Further work, however, is still required to understand and predict the potential human infectivity of pathogenic organisms alongside the better design of early warning systems and surveillance measures for risk assessment purposes.

Ecology and evolution of viruses infecting uncultivated SUP05 bacteria as revealed by single-cell- and meta-genomics

Viruses modulate microbial communities and alter ecosystem functions. However, due to cultivation bottlenecks, specific virus-host interaction dynamics remain cryptic. In this study, we examined 127 single-cell amplified genomes (SAGs) from uncultivated SUP05 bacteria isolated from a model marine oxygen minimum zone (OMZ) to identify 69 viral contigs representing five new genera within dsDNA Caudovirales and ssDNA Microviridae. Infection frequencies suggest that ∼1/3 of SUP05 bacteria is viral-infected, with higher infection frequency where oxygen-deficiency was most severe. Observed Microviridae clonality suggests recovery of bloom-terminating viruses, while systematic co-infection between dsDNA and ssDNA viruses posits previously unrecognized cooperation modes. Analyses of 186 microbial and viral metagenomes revealed that SUP05 viruses persisted for years, but remained endemic to the OMZ. Finally, identification of virus-encoded dissimilatory sulfite reductase suggests SUP05 viruses reprogram their host's energy metabolism. Together, these results demonstrate closely coupled SUP05 virus-host co-evolutionary dynamics with the potential to modulate biogeochemical cycling in climate-critical and expanding OMZs.

Depth-stratified functional and taxonomic niche specialization in the 'core' and 'flexible' Pacific Ocean Virome

Microbes drive myriad ecosystem processes, and their viruses modulate microbial-driven processes through mortality, horizontal gene transfer, and metabolic reprogramming by viral-encoded auxiliary metabolic genes (AMGs). However, our knowledge of viral roles in the oceans is primarily limited to surface waters. Here we assess the depth distribution of protein clusters (PCs) in the first large-scale quantitative viral metagenomic data set that spans much of the pelagic depth continuum (the Pacific Ocean Virome; POV). This established 'core' (180 PCs; one-third new to science) and 'flexible' (423K PCs) community gene sets, including niche-defining genes in the latter (385 and 170 PCs are exclusive and core to the photic and aphotic zones, respectively). Taxonomic annotation suggested that tailed phages are ubiquitous, but not abundant (<5% of PCs) and revealed depth-related taxonomic patterns. Functional annotation, coupled with extensive analyses to document non-viral DNA contamination, uncovered 32 new AMGs (9 core, 20 photic and 3 aphotic) that introduce ways in which viruses manipulate infected host metabolism, and parallel depth-stratified host adaptations (for example, photic zone genes for iron-sulphur cluster modulation for phage production, and aphotic zone genes for high-pressure deep-sea survival). Finally, significant vertical flux of photic zone viruses to the deep sea was detected, which is critical for interpreting depth-related patterns in nature. Beyond the ecological advances outlined here, this catalog of viral core, flexible and niche-defining genes provides a resource for future investigation into the organization, function and evolution of microbial molecular networks to mechanistically understand and model viral roles in the biosphere.

Diverse small circular single-stranded DNA viruses identified in a freshwater pond on the McMurdo Ice Shelf (Antarctica)

Antarctica has some of the harshest environmental conditions for existence of life on Earth. In this pilot study we recovered eight diverse circular single-stranded DNA (ssDNA) viral genome sequences (1904-3120 nts) from benthic mats dominated by filamentous cyanobacteria in a freshwater pond on the McMurdo Ice Shelf sampled in 1988. All genomes contain two to three major open reading frames (ORFs) that are uni- or bi-directionally transcribed and all have an ORF encoding a replication-associated protein (Rep). In one genome, the second ORF has similarity to a capsid protein (CP) of Nepavirus which is most closely related to geminiviruses. Additionally, all genomes have two intergenic regions that contain putative stem loop structures, six genomes have NANTATTAC as the nonanucleotide motif, while one has CCTTATTAC, and another has a non-canonical stem loop. In the large intergenic region, we identified iterative sequences flanking the putative stem-loop elements which are a hallmark of most circular ssDNA viruses encoding rolling circle replication (RCR) initiators of the HUH endonuclease superfamily. The Reps encoded by ssDNA viral genomes recovered in this study shared <38% pairwise identity to all other Reps of known ssDNA viruses. A previous study on Lake Limnopolar (Livingston Island, South Shetland Islands), using next-generation sequencing identified circular ssDNA viruses and their putative Reps share <35% pairwise identity to those from the viral genomes removed in this study. It is evident from our pilot study that the global diversity of ssDNA viruses is grossly underestimated and there is limited knowledge on ssDNA viruses in Antarctica.

Tales from a thousand and one phages

The sequencing of marine metagenomic fosmids led to the discovery of several new complete phage genomes. Among the 21 major sequence groups, 10 totally novel groups of marine phages could be identified. Some of these represent the first phages infecting large marine prokaryotic phyla, such as the Verrucomicrobia and the recently described Ca. Actinomarinales. Coming from a single deep photic zone sample the diversity of phages found is astonishing, and the comparison with a metavirome from the same location indicates that only 2% of the real diversity was recovered. In addition to this large macro-diversity, rich micro-diversity was also found, affecting host-recognition modules, mirroring the variation of cell surface components in their host marine microbes.

Evidence for metaviromic islands in marine phages

Metagenomic islands (MGIs) have been defined as genomic regions in prokaryotic genomes that under-recruit from metagenomes where most of the same genome recruits at close to 100% identity over most of its length. The presence of MGIs in prokaryotes has been associated to the diversity of concurrent lineages that vary at this level to disperse the predatory pressure of phages that, reciprocally, maintain high clonal diversity in the population and improve ecosystem performance. This was proposed as a Constant-Diversity (C-D) model. Here we have investigated the regions of phage genomes under-recruiting in a metavirome constructed with a sample from the same habitat where they were retrieved. Some of the genes found to under-recruit are involved in host recognition as would be expected from the C-D model. Furthermore, the recruitment of intragenic regions known to be involved in molecular recognition also had a significant under-recruitment compared to the rest of the gene. However, other genes apparently disconnected from the recognition process under-recruited often, specifically the terminases involved in packaging of the phage genome in the capsid and a few others. In addition, some highly related phage genomes (at nucleotide sequence level) had no metaviromic islands (MVIs). We speculate that the latter might be generalist phages with broad infection range that do not require clone specific lineages.