Marine phage genomics: the tip of the iceberg

The long and short of it: benchmarking viromics using Illumina, Nanopore and PacBio sequencing technologies

Viral metagenomics has fuelled a rapid change in our understanding of global viral diversity and ecology. Long-read sequencing and hybrid assembly approaches that combine long- and short-read technologies are now being widely implemented in bacterial genomics and metagenomics. However, the use of long-read sequencing to investigate viral communities is still in its infancy. While Nanopore and PacBio technologies have been applied to viral metagenomics, it is not known to what extent different technologies will impact the reconstruction of the viral community. Thus, we constructed a mock bacteriophage community of previously sequenced phage genomes and sequenced them using Illumina, Nanopore and PacBio sequencing technologies and tested a number of different assembly approaches. When using a single sequencing technology, Illumina assemblies were the best at recovering phage genomes. Nanopore- and PacBio-only assemblies performed poorly in comparison to Illumina in both genome recovery and error rates, which both varied with the assembler used. The best Nanopore assembly had errors that manifested as SNPs and INDELs at frequencies 41 and 157 % higher than found in Illumina only assemblies, respectively. While the best PacBio assemblies had SNPs at frequencies 12 and 78 % higher than found in Illumina-only assemblies, respectively. Despite high-read coverage, long-read-only assemblies recovered a maximum of one complete genome from any assembly, unless reads were down-sampled prior to assembly. Overall the best approach was assembly by a combination of Illumina and Nanopore reads, which reduced error rates to levels comparable with short-read-only assemblies. When using a single technology, Illumina only was the best approach. The differences in genome recovery and error rates between technology and assembler had downstream impacts on gene prediction, viral prediction, and subsequent estimates of diversity within a sample. These findings will provide a starting point for others in the choice of reads and assembly algorithms for the analysis of viromes.

Revealing viral diversity in the Napahai plateau wetland based on metagenomics

We focused on exploring the diversity of viruses in the Napahai plateau wetland, a unique ecosystem located in Yunnan, China. While viruses in marine environments have been extensively studied for their influence on microbial metabolism and biogeochemical cycles, little is known about their composition and function in plateau wetlands. Metagenomic analysis was employed to investigate the viral diversity and biogeochemical impacts in the Napahai wetland. It revealed that the Caudoviricetes and Malgrandaviricetes class level was the most abundant viral category based on phylogenetic analysis. Additionally, a gene-sharing network highlighted the presence of numerous unexplored viruses and demonstrated their unique characteristics and significant variation within the viral community of the Napahai wetland. Furthermore, the study identified the auxiliary metabolic genes (AMGs). AMGs provide phages with additional functions, such as protection against host degradation and involvement in metabolic pathways, such as the pentose phosphate pathway and DNA biosynthesis. The viruses in the Napahai wetland were found to influence carbon, nitrogen, sulfur, and amino acid metabolism, indirectly contributing to biogeochemical cycling through these AMGs. Overall, the research sheds light on the diverse and unique viral communities in the Napahai plateau wetland and emphasizes the significant roles of viruses in microbial ecology. The findings contribute to a deeper understanding of the characteristics and ecological functions of viral communities in plateau wetland ecosystems.

Distantly related Alteromonas bacteriophages share tail fibers exhibiting properties of transient chaperone caps

The host recognition modules encoding the injection machinery and receptor binding proteins (RBPs) of bacteriophages are predisposed to mutation and recombination to maintain infectivity towards co-evolving bacterial hosts. In this study, we reveal how Alteromonas mediterranea schitovirus A5 shares its host recognition module, including tail fiber and cognate chaperone, with phages from distantly related families including Alteromonas myovirus V22. While the V22 chaperone is essential for producing active tail fibers, here we demonstrate production of functional A5 tail fibers regardless of chaperone co-expression. AlphaFold-generated models of tail fiber and chaperone pairs from phages A5, V22, and other Alteromonas phages reveal how amino acid insertions within both A5-like proteins results in a knob domain duplication in the tail fiber and a chaperone β-hairpin "tentacle" extension. These structural modifications are linked to differences in chaperone dependency between the A5 and V22 tail fibers. Structural similarity between the chaperones and intramolecular chaperone domains of other phage RBPs suggests an additional function of these chaperones as transient fiber "caps". Finally, our identification of homologous host recognition modules from morphologically distinct phages implies that horizontal gene transfer and recombination events between unrelated phages may be a more common process than previously thought among Caudoviricetes phages.

Diversity, evolution and life strategies of CbK-like phages

Caulobacter phage CbK has been extensively studied as a model system in virology and bacteriology. Lysogeny-related genes have been found in each CbK-like isolate, suggesting a life strategy of both lytic and lysogenic cycles. However, whether CbK-related phages can enter lysogeny is still undetermined. This study identified new CbK-like sequences and expanded the collection of CbK-related phages. A common ancestry with a temperate lifestyle was predicted for the group, however, which subsequently evolved into two clades of different genome sizes and host associations. Through the examination of phage recombinase genes, alignment of attachment sites on the phage and bacterial genomes (attP-attB pairing), and the experimental validation, different lifestyles were found among the different members. A majority of clade II members retain a lysogenic lifestyle, whereas all clade I members have evolved into an obligate lytic lifestyle via a loss of the gene encoding Cre-like recombinase and the coupled attP fragment. We postulated that the loss of lysogeny may be a by-product of the increase in phage genome size, and vice versa. Clade I is likely to overcome the costs through maintaining more auxiliary metabolic genes (AMGs), particularly for those involved in protein metabolism, to strengthen host takeover and further benefit virion production.

Improving high-throughput techniques for bacteriophage discovery in multi-well plates

Bacteriophages (also called phages) are viruses of bacteria that have numerous applications in medicine, agriculture, ecology, and molecular biology. With the increasing interest in phages for their many uses, it is now especially important to make phage discovery more efficient and economical. Using the host Mycobacterium smegmatis mc²155, which is a model organism for phage discovery research and is closely related to important pathogens of humans and other animals, we investigated three procedures that are an integral part of phage discovery: enrichment of environmental samples, phage isolation and detection (which can also be used for host range determination), and phage purification. Enrichment in 6-well plates was successful with most environmental samples, and enrichment in 24- and 96-well plates was successful with some environmental samples, demonstrating that larger sample volumes are preferred when possible, but smaller sample volumes may be acceptable if the starting concentration of phages is sufficiently high. Measuring absorbance in multi-well plates was at least as sensitive as the traditional plaque assay for the detection of phages. We also demonstrated a technique for the purification of single phage types from mixed cultures in liquid medium. Multi-well techniques can be used as alternatives or complementary approaches to traditional methods of phage discovery and characterization depending on the needs of the researcher in terms of time, available resources, host species, phage-bacteria matches, and specific goals. In the future, these techniques could be applied to the discovery of phages of aquatic mycobacteria and other hosts for which few phages have currently been isolated.

Comparative genomic analysis of dwarf Vibrio myoviruses defines a conserved gene cluster for successful phage infection

Tailed bacteriophages have been at the center of attention, not only for their ability to infect and kill pathogenic bacteria but also due to their peculiar and intriguing complex contractile tail structure. Tailed bacteriophages with contractile tails are known to have a Myoviridae morphotype and are members of the order Caudovirales. Large bacteriophages with a genome larger than 150 kbp have been studied for their ability to use multiple infection and lysis strategies to replicate more efficiently. On the other hand, smaller bacteriophages with fewer genes are represented in the GenBank database in greater numbers, and have several genes with unknown function. Isolation and molecular characterization of a newly reported bacteriophage named Athena1 revealed that it is a strongly lytic bacteriophage with a genome size of 39,826 bp. This prompted us to perform a comparative genomic analysis of Vibrio myoviruses with a genome size of no more than 50 kbp. The results revealed a pattern of genomic organization that includes sets of genes responsible for virion morphogenesis, replication/recombination of DNA, and lysis/lysogeny switching. By studying phylogenetic gene markers, we were able to draw conclusions about evolutionary events that shaped the genomic mosaicism of these phages, pinpointing the importance of a conserved organization of the genomic region encoding the baseplate protein for successful infection of Gram-negative bacteria. In addition, we propose the creation of new genera for dwarf Vibrio myoviruses. Comparative genomics of phages infecting aquatic bacteria could provide information that is useful for combating fish pathogens in aquaculture, using novel strategies.

INfrastructure for a PHAge REference Database: Identification of Large-Scale Biases in the Current Collection of Cultured Phage Genomes

Background: With advances in sequencing technology and decreasing costs, the number of phage genomes that have been sequenced has increased markedly in the past decade. Materials and Methods: We developed an automated retrieval and analysis system for phage genomes (https://github.com/RyanCook94/inphared) to produce the INfrastructure for a PHAge REference Database (INPHARED) of phage genomes and associated metadata. Results: As of January 2021, 14,244 complete phage genomes have been sequenced. The INPHARED data set is dominated by phages that infect a small number of bacterial genera, with 75% of phages isolated on only 30 bacterial genera. There is further bias, with significantly more lytic phage genomes (∼70%) than temperate (∼30%) within our database. Collectively, this results in ∼54% of temperate phage genomes originating from just three host genera. With much debate on the carriage of antibiotic resistance genes and their potential safety in phage therapy, we searched for putative antibiotic resistance genes. Frequency of antibiotic resistance gene carriage was found to be higher in temperate phages than in lytic phages and again varied with host. Conclusions: Given the bias of currently sequenced phage genomes, we suggest to fully understand phage diversity, efforts should be made to isolate and sequence a larger number of phages, in particular temperate phages, from a greater diversity of hosts.

Identification of the tail assembly chaperone genes of T4-Like phages suggests a mechanism other than translational frameshifting for biogenesis of their encoded proteins

Tape measure (TM) proteins are essential for the formation of long-tailed phages. TM protein assembly into tails requires the action of tail assembly chaperones (TACs). TACs (e.g. gpG and gpT of E. coli phage lambda) are usually produced in a short (TAC-N) and long form (TAC-NC) with the latter comprised of TAC-N with an additional C-terminal domain (TAC-C). TAC-NC is generally synthesized through a ribosomal frameshifting mechanism. TAC encoding genes have never been identified in the intensively studied Escherichia coli phage T4, or any related phages. Here, we have bioinformatically identified putative TAC encoding genes in diverse T4-like phage genomes. The frameshifting mechanism for producing TAC-NC appears to be conserved in several T4-like phage groups. However, the group including phage T4 itself likely employs a different strategy whereby TAC-N and TAC-NC are encoded by separate genes (26 and 51 in phage T4).

Viruses infecting a warm water picoeukaryote shed light on spatial co-occurrence dynamics of marine viruses and their hosts

The marine picoeukaryote Bathycoccus prasinos has been considered a cosmopolitan alga, although recent studies indicate two ecotypes exist, Clade BI (B. prasinos) and Clade BII. Viruses that infect Bathycoccus Clade BI are known (BpVs), but not that infect BII. We isolated three dsDNA prasinoviruses from the Sargasso Sea against Clade BII isolate RCC716. The BII-Vs do not infect BI, and two (BII-V2 and BII-V3) have larger genomes (~210 kb) than BI-Viruses and BII-V1. BII-Vs share ~90% of their proteins, and between 65% to 83% of their proteins with sequenced BpVs. Phylogenomic reconstructions and PolB analyses establish close-relatedness of BII-V2 and BII-V3, yet BII-V2 has 10-fold higher infectivity and induces greater mortality on host isolate RCC716. BII-V1 is more distant, has a shorter latent period, and infects both available BII isolates, RCC716 and RCC715, while BII-V2 and BII-V3 do not exhibit productive infection of the latter in our experiments. Global metagenome analyses show Clade BI and BII algal relative abundances correlate positively with their respective viruses. The distributions delineate BI/BpVs as occupying lower temperature mesotrophic and coastal systems, whereas BII/BII-Vs occupy warmer temperature, higher salinity ecosystems. Accordingly, with molecular diagnostic support, we name Clade BII Bathycoccus calidus sp. nov. and propose that molecular diversity within this new species likely connects to the differentiated host-virus dynamics observed in our time course experiments. Overall, the tightly linked biogeography of Bathycoccus host and virus clades observed herein supports species-level host specificity, with strain-level variations in infection parameters.

Bacteriophage classification for assembled contigs using graph convolutional network

Motivation

Bacteriophages (aka phages), which mainly infect bacteria, play key roles in the biology of microbes. As the most abundant biological entities on the planet, the number of discovered phages is only the tip of the iceberg. Recently, many new phages have been revealed using high-throughput sequencing, particularly metagenomic sequencing. Compared to the fast accumulation of phage-like sequences, there is a serious lag in taxonomic classification of phages. High diversity, abundance and limited known phages pose great challenges for taxonomic analysis. In particular, alignment-based tools have difficulty in classifying fast accumulating contigs assembled from metagenomic data.

Results

In this work, we present a novel semi-supervised learning model, named PhaGCN, to conduct taxonomic classification for phage contigs. In this learning model, we construct a knowledge graph by combining the DNA sequence features learned by convolutional neural network and protein sequence similarity gained from gene-sharing network. Then we apply graph convolutional network to utilize both the labeled and unlabeled samples in training to enhance the learning ability. We tested PhaGCN on both simulated and real sequencing data. The results clearly show that our method competes favorably against available phage classification tools.

Availability and implementation

The source code of PhaGCN is available via: https://github.com/KennthShang/PhaGCN.

Potential virus-mediated nitrogen cycling in oxygen-depleted oceanic waters

Viruses play an important role in the ecology and biogeochemistry of marine ecosystems. Beyond mortality and gene transfer, viruses can reprogram microbial metabolism during infection by expressing auxiliary metabolic genes (AMGs) involved in photosynthesis, central carbon metabolism, and nutrient cycling. While previous studies have focused on AMG diversity in the sunlit and dark ocean, less is known about the role of viruses in shaping metabolic networks along redox gradients associated with marine oxygen minimum zones (OMZs). Here, we analyzed relatively quantitative viral metagenomic datasets that profiled the oxygen gradient across Eastern Tropical South Pacific (ETSP) OMZ waters, assessing whether OMZ viruses might impact nitrogen (N) cycling via AMGs. Identified viral genomes encoded six N-cycle AMGs associated with denitrification, nitrification, assimilatory nitrate reduction, and nitrite transport. The majority of these AMGs (80%) were identified in T4-like Myoviridae phages, predicted to infect Cyanobacteria and Proteobacteria, or in unclassified archaeal viruses predicted to infect Thaumarchaeota. Four AMGs were exclusive to anoxic waters and had distributions that paralleled homologous microbial genes. Together, these findings suggest viruses modulate N-cycling processes within the ETSP OMZ and may contribute to nitrogen loss throughout the global oceans thus providing a baseline for their inclusion in the ecosystem and geochemical models.

The genomic content and context of auxiliary metabolic genes in roseophages

Marine bacteriophages frequently possess auxiliary metabolic genes (AMGs) that accelerate host metabolism during phage infection. The significance of AMGs in phage infecting the ecologically important Roseobacter clade, found predominantly in marine environments, remains to be determined. Here, we analysed the distribution and genomic context of 180 AMGs, annotated into 20 types, across 50 roseophage genomes. Roseophages share seven high-frequency AMGs (trx, grx, RNR, thyX, DCD, phoH, and mazG)， most of them involved in the nucleotide biosynthesis pathway that represent conserved intra and inter operational taxonomic units (OTUs), and share ≥97% full-length DNA sequence similarity. Sporadic AMGs (dUTPase, lexA, degS, Que, NAPRT, AHL, pcnB, ctrA, RTX, RNR-nrdA, RNR-nrdE, wclP, and flgJ), present in only one or two OTUs, show high functional diversity. The roseophage AMG repertoire weakly correlates with environmental factors, while host range partially explains the sporadic AMG distribution. Locally co-linear blocks distribution index (LDI) analysis indicated that high-frequency roseopodovirus AMGs are restricted to particular genomic islands, possibly originating from limited historical acquisition events. Low-frequency roseopodovirus AMGs and all roseosiphovirus AMGs have high LDI values, implying multiple historical acquisition events. In summary, roseophages have acquired a range of AMGs through horizontal gene transfer, and the forces shaping the evolution of roseophages are described.

A plate of viruses: Viral metagenomics of supermarket chicken, pork and beef from Brazil

A viral metagenomics study was conducted in beef, pork, and chicken sold in supermarkets from Southern Brazil. From chicken, six distinct gyroviruses (GyV) were detected, including GyV3 and GyV6, which for the first time were detected in samples from avian species, plus a novel smacovirus species and two highly divergent circular Rep-encoding ssDNA (CRESS-DNA) viruses. From pork, genomes of numerous anelloviruses, porcine parvovirus 5 (PPV5) and 6 (PPV6), two new genomoviruses and two new CRESS-DNA viruses were found. Finally, two new CRESS-DNA genomes were recovered from beef. Although none of these viruses have history of transmission to humans, the findings reported here reveal that such agents are inevitably consumed in diets that include these types of meat.

Characterization and Genome Analysis of a Novel Mu-like Phage VW-6B Isolated from the Napahai Plateau Wetland of China

Although bacteriophages are more numerous and have smaller genomes than their bacterial hosts, relatively few have their genomes sequenced. Here, we isolated the Pseudomonas fluorescens bacteriophage from Napahai plateau wetland and performed de novo genome sequencing. Based on the previous biological characteristics and bioinformatics analysis, it was determined that VW-6B was a linear double-stranded DNA (dsDNA) phage with 35,306 bp, with 56.76% G+C content and 197 bp tandem repeats. The VW-6B genome contained 46 open-reading frames (ORFs), and no tRNA genes were found. Based on phage genome structure, sequence comparison, and collinear analysis, VW-6B should be classified into the family Siphoviridae and be considered as a member of a new species in the Mu-like phage. The newly isolated bacteriophage can specifically infect P. fluorescens, which further enriches the diversity of known bacteriophages and provides a basis for the subsequent research and application of bacteriophages.

Lysogeny in the oceans: Lessons from cultivated model systems and a reanalysis of its prevalence

In the oceans, viruses that infect bacteria (phages) influence a variety of microbially mediated processes that drive global biogeochemical cycles. The nature of their influence is dependent upon infection mode, be it lytic or lysogenic. Temperate phages are predicted to be prevalent in marine systems where they are expected to execute both types of infection modes. Understanding the range and outcomes of temperate phage-host interactions is fundamental for evaluating their ecological impact. Here, we (i) review phage-mediated rewiring of host metabolism, with a focus on marine systems, (ii) consider the range and nature of temperate phage-host interactions, and (iii) draw on studies of cultivated model systems to examine the consequences of lysogeny among several dominant marine bacterial lineages. We also readdress the prevalence of lysogeny among marine bacteria by probing a collection of 1239 publicly available bacterial genomes, representing cultured and uncultivated strains, for evidence of complete prophages. Our conservative analysis, anticipated to underestimate true prevalence, predicts 18% of the genomes examined contain at least one prophage, the majority (97%) were found within genomes of cultured isolates. These results highlight the need for cultivation of additional model systems to better capture the diversity of temperate phage-host interactions in the oceans.

Alteromonas Myovirus V22 Represents a New Genus of Marine Bacteriophages Requiring a Tail Fiber Chaperone for Host Recognition

Marine phages play a variety of critical roles in regulating the microbial composition of our oceans. Despite constituting the majority of genetic diversity within these environments, there are relatively few isolates with complete genome sequences or in-depth analyses of their host interaction mechanisms, such as characterization of their receptor binding proteins (RBPs). Here, we present the 92,760-bp genome of the Alteromonas-targeting phage V22. Genomic and morphological analyses identify V22 as a myovirus; however, due to a lack of sequence similarity to any other known myoviruses, we propose that V22 be classified as the type phage of a new Myoalterovirus genus within the Myoviridae family. V22 shows gene homology and synteny with two different subfamilies of phages infecting enterobacteria, specifically within the structural region of its genome. To improve our understanding of the V22 adsorption process, we identified putative RBPs (gp23, gp24, and gp26) and tested their ability to decorate the V22 propagation strain, Alteromonas mediterranea PT11, as recombinant green fluorescent protein (GFP)-tagged constructs. Only GFP-gp26 was capable of bacterial recognition and identified as the V22 RBP. Interestingly, production of functional GFP-gp26 required coexpression with the downstream protein gp27. GFP-gp26 could be expressed alone but was incapable of host recognition. By combining size-exclusion chromatography with fluorescence microscopy, we reveal how gp27 is not a component of the final RBP complex but instead is identified as a new type of phage-encoded intermolecular chaperone that is essential for maturation of the gp26 RBP.IMPORTANCE Host recognition by phage-encoded receptor binding proteins (RBPs) constitutes the first step in all phage infections and the most critical determinant of host specificity. By characterizing new types of RBPs and identifying their essential chaperones, we hope to expand the repertoire of known phage-host recognition machineries. Due to their genetic plasticity, studying RBPs and their associated chaperones can shed new light onto viral evolution affecting phage-host interactions, which is essential for fields such as phage therapy or biotechnology. In addition, since marine phages constitute one of the most important reservoirs of noncharacterized genetic diversity on the planet, their genomic and functional characterization may be of paramount importance for the discovery of novel genes with potential applications.

Dead or alive: sediment DNA archives as tools for tracking aquatic evolution and adaptation

DNA can be preserved in marine and freshwater sediments both in bulk sediment and in intact, viable resting stages. Here, we assess the potential for combined use of ancient, environmental, DNA and timeseries of resurrected long-term dormant organisms, to reconstruct trophic interactions and evolutionary adaptation to changing environments. These new methods, coupled with independent evidence of biotic and abiotic forcing factors, can provide a holistic view of past ecosystems beyond that offered by standard palaeoecology, help us assess implications of ecological and molecular change for contemporary ecosystem functioning and services, and improve our ability to predict adaptation to environmental stress.

Ellegaard et al. discuss the potential for using ancient environmental DNA (eDNA), combined with resurrection ecology, to analyse trophic interactions and evolutionary adaptation to changing environments. Their Review suggests that these techniques will improve our ability to predict genetic and phenotypic adaptation to environmental stress.

Game-Theoretical Modeling of Interviral Conflicts Mediated by Mini-CRISPR Arrays

All cellular organisms coevolve with multiple viruses, so that both virus-host and intervirus conflicts are major factors of evolution. Accordingly, hosts evolve multiple, elaborate defense systems and viruses respond by evolving means of antidefense. Although less thoroughly characterized, several dedicated mechanisms of intervirus competition have been described as well. Recently, the genomes of some bacterial and archaeal viruses have been shown to harbor CRISPR mini-arrays that typically contain a single spacer targeting a closely related virus. The involvement of mini-arrays in an intervirus conflict has been experimentally demonstrated for a pair of archaeal viruses. We model the evolution of virus-encoded CRISPR mini-arrays using a game theoretical approach. Analysis of the model reveals multiple equilibria that include mutual targeting, unidirectional targeting, no targeting, cyclic polymorphism, and bistability. The choice between these evolutionary regimes depends on the model parameters including the coinfection frequency, differential productivity of the conflicting viruses, and the fitness cost of mini-arrays. At high coinfection frequencies, the model becomes a version of the Prisoner’s dilemma in which defection, i.e., mutual targeting between the competing viruses, is the winning strategy.

Phage Therapy with a Focus on the Human Microbiota

Bacteriophages are viruses that infect bacteria. After their discovery in the early 1900s, bacteriophages were a primary cure against infectious disease for almost 25 years, before being completely overshadowed by antibiotics. With the rise of antibiotic resistance, bacteriophages are being explored again for their antibacterial activity. One of the critical apprehensions regarding bacteriophage therapy, however, is the possibility of genome evolution, development of phage resistance, and subsequent perturbations to our microbiota. Through this review, we set out to explore the principles supporting the use of bacteriophages as a therapeutic agent, discuss the human gut microbiome in relation to the utilization of phage therapy, and the co-evolutionary arms race between host bacteria and phage in the context of the human microbiota.

Bacteriophages in Natural and Artificial Environments

Bacteriophages (phages) are biological entities that have attracted a great deal of attention in recent years. They have been reported as the most abundant biological entities on the planet and their ability to impact the composition of bacterial communities is of great interest. In this review, we aim to explore where phages exist in natural and artificial environments and how they impact communities. The natural environment in this review will focus on the human body, soils, and the marine environment. In these naturally occurring environments there is an abundance of phages suggesting a role in the maintenance of bacterial community homeostasis. The artificial environment focuses on wastewater treatment plants, industrial processes, followed by pharmaceutical formulations. As in natural environments, the existence of bacteria in manmade wastewater treatment plants and industrial processes inevitably attracts phages. The presence of phages in these environments can inhibit the bacteria required for efficient water treatment or food production. Alternatively, they can have a positive impact by eliminating recalcitrant organisms. Finally, we conclude by describing how phages can be manipulated or formulated into pharmaceutical products in the laboratory for use in natural or artificial environments.

Riding the wave of genomics to investigate aquatic coliphage diversity and activity

Bacteriophages infecting Escherichia coli (coliphages) have been used as a proxy for faecal matter and water quality from a variety of environments. However, the diversity of coliphages that is present in seawater remains largely unknown, with previous studies largely focusing on morphological diversity. Here, we isolated and characterized coliphages from three coastal locations in the United Kingdom and Poland. Comparative genomics and phylogenetic analysis of phage isolates facilitated the identification of putative new species within the genera Rb69virus and T5virus and a putative new genus within the subfamily Tunavirinae. Furthermore, genomic and proteomic analysis combined with host range analysis allowed the identification of a putative tail fibre that is likely responsible for the observed differences in host range of phages vB_Eco_mar003J3 and vB_Eco_mar004NP2.

Dynamic marine viral infections and major contribution to photosynthetic processes shown by spatiotemporal picoplankton metatranscriptomes

Viruses provide top-down control on microbial communities, yet their direct study in natural environments was hindered by culture limitations. The advance of bioinformatics enables cultivation-independent study of viruses. Many studies assemble new viral genomes and study viral diversity using marker genes from free viruses. Here we use cellular metatranscriptomics to study active community-wide viral infections. Recruitment to viral contigs allows tracking infection dynamics over time and space. Our assemblies represent viral populations, but appear biased towards low diversity viral taxa. Tracking relatives of published T4-like cyanophages and pelagiphages reveals high genomic continuity. We determine potential hosts by matching dynamics of infection with abundance of particular microbial taxa. Finally, we quantify the relative contribution of cyanobacteria and viruses to photosystem-II psbA (reaction center) expression in our study sites. We show sometimes >50% of all cyanobacterial+viral psbA expression is of viral origin, highlighting the contribution of viruses to photosynthesis and oxygen production.

Here, Sieradzki et al. use metatranscriptomics to study active community-wide viral infections at three coastal California sites throughout a year, identify potential viral hosts, and show that viruses can contribute a substantial amount to photosystem-II psbA expression.

Extremophile deep-sea viral communities from hydrothermal vents: Structural and functional analysis

Ten publicly available metagenomic data sets from hydrothermal vents were analyzed to determine the taxonomic structure of the viral communities present, as well as their potential metabolic functions. The type of natural selection on two auxiliary metabolic genes was also analyzed. The structure of the virome in the hydrothermal vents was quite different in comparison with the viruses present in sediments, with specific populations being present in greater abundance in the plume samples when compared with the sediment samples. ssDNA genomes such as Circoviridae and Microviridae were predominantly present in the sediment samples, with Caudovirales which are dsDNA being present in the vent samples. Genes potentially encoding enzymes that participate in carbon, nitrogen and sulfur metabolic pathways were found in greater abundance, than those involved in the oxygen cycle, in the hydrothermal vents. Functional profiling of the viromes, resulted in the discovery of genes encoding proteins involved in bacteriophage capsids, DNA synthesis, nucleotide synthesis, DNA repair, as well as viral auxiliary metabolic genes such as cytitidyltransferase and ribonucleotide reductase. These auxiliary metabolic genes participate in the synthesis of phospholipids and nucleotides respectively and are likely to contribute to enhancing the fitness of their bacterial hosts within the hydrothermal vent communities. Finally, evolutionary analysis suggested that these auxiliary metabolic genes are highly conserved and evolve under purifying selection, and are thus maintained in their genome.

Cylindrospermopsis raciborskii Virus and host: genomic characterization and ecological relevance

Cylindrospermopsis (Raphidiopsis) raciborskii is an invasive, filamentous, nitrogen-fixing cyanobacterium that forms frequent blooms in freshwater habitats. While viruses play key roles in regulating the abundance, production and diversity of their hosts in aquatic ecosystems, the role(s) of viruses in the ecology of C. raciborskii is almost unexplored. Progress in this field has been hindered by the absence of a characterized virus-host system in C. raciborskii. To bridge this gap, we sequenced the genome of CrV-01T, a previously isolated cyanosiphovirus, and its host, C. raciborskii strain Cr2010. Analyses suggest that CrV-01T represents a distinct clade of siphoviruses infecting, and perhaps lysogenizing, filamentous cyanobacteria. Its genome contains unique features that include an intact CRISPR array and a 12 kb inverted duplication. Evidence suggests CrV-01T recently gained the ability to infect Cr2010 and recently lost the ability to form lysogens. The cyanobacterial host contains a CRISPR-Cas system with CRISPR spacers matching protospacers within the inverted duplication of the CrV-01T genome. Examination of metagenomes demonstrates that viruses with high genetic identity to CrV-01T, but lacking the inverted duplication, are present in C. raciborskii blooms in Australia. The unique genomic features of the CrV/Cr2010 system offers opportunities to investigate in more detail virus-host interactions in an ecologically important bloom-forming cyanobacterium.

Comparative Genomics of Bacteriophage of the Genus Seuratvirus

Despite being more abundant and having smaller genomes than their bacterial host, relatively few bacteriophages have had their genomes sequenced. Here, we isolated 14 bacteriophages from cattle slurry and performed de novo genome sequencing, assembly, and annotation. The commonly used marker genes polB and terL showed these bacteriophages to be closely related to members of the genus Seuratvirus. We performed a core-gene analysis using the 14 new and four closely related genomes. A total of 58 core genes were identified, the majority of which has no known function. These genes were used to construct a core-gene phylogeny, the results of which confirmed the new isolates to be part of the genus Seuratvirus and expanded the number of species within this genus to four. All bacteriophages within the genus contained the genes queCDE encoding enzymes involved in queuosine biosynthesis. We suggest these genes are carried as a mechanism to modify DNA in order to protect these bacteriophages against host endonucleases.

The smallest ssDNA phage infecting a marine bacterium

In the marine environment, only a few lytic single-stranded DNA (ssDNA) phages have been isolated and characterized, despite the fact that diverse ssDNA bacteriophages have been discovered via metagenomic studies. In this study, we isolated and characterized a new ssDNA phage, vB_RpoMi-Mini, which infects a marine bacterium Ruegeria pomeroyi DSS-3. With a genome size of 4248 bp and only four putative open reading frames (ORF), vB_RpoMi-Mini becomes the smallest ssDNA phage among the known ssDNA phage isolates and represents the DNA bacteriophage with the least number of ORFs. Genome-wide analysis reveals that bacteriophage Mini is distantly related to the known ssDNA phages and belongs to an unclassified ssDNA phage within the Microviridae family. The presence of peptidase in vB_RpoMi-Mini genome further implies that horizontal gene transfer could be an important driving force in the evolution of ssDNA phages. Bacteriophage Mini seems to have lost the spike protein commonly seen in ssDNA phages, suggesting that ssDNA phage can be more diverse than previously thought. Metagenomic analysis indicates that Mini-like phages are widely distributed in the environments. The discovery of vB_RpoMi-Mini expands our understanding of ssDNA phages in nature, and also indicates our dearth of knowledge regarding of ssDNA phages.

Phage spanins: diversity, topological dynamics and gene convergence

BACKGROUND

Spanins are phage lysis proteins required to disrupt the outer membrane. Phages employ either two-component spanins or unimolecular spanins in this final step of Gram-negative host lysis. Two-component spanins like Rz-Rz1 from phage lambda consist of an integral inner membrane protein: i-spanin, and an outer membrane lipoprotein: o-spanin, that form a complex spanning the periplasm. Two-component spanins exist in three different genetic architectures; embedded, overlapped and separated. In contrast, the unimolecular spanins, like gp11 from phage T1, have an N-terminal lipoylation signal sequence and a C-terminal transmembrane domain to account for the topology requirements. Our proposed model for spanin function, for both spanin types, follows a common theme of the outer membrane getting fused with the inner membrane, effecting the release of progeny virions.

RESULTS

Here we present a SpaninDataBase which consists of 528 two-component spanins and 58 unimolecular spanins identified in this analysis. Primary analysis revealed significant differences in the secondary structure predictions for the periplasmic domains of the two-component and unimolecular spanin types, as well as within the three different genetic architectures of the two-component spanins. Using a threshold of 40% sequence identity over 40% sequence length, we were able to group the spanins into 143 i-spanin, 125 o-spanin and 13 u-spanin families. More than 40% of these families from each type were singletons, underlining the extreme diversity of this class of lysis proteins. Multiple sequence alignments of periplasmic domains demonstrated conserved secondary structure patterns and domain organization within family members. Furthermore, analysis of families with members from different architecture allowed us to interpret the evolutionary dynamics of spanin gene arrangement. Also, the potential universal role of intermolecular disulfide bonds in two-component spanin function was substantiated through bioinformatic and genetic approaches. Additionally, a novel lipobox motif, AWAC, was identified and experimentally verified.

CONCLUSIONS

The findings from this bioinformatic approach gave us instructive insights into spanin function, evolution, domain organization and provide a platform for future spanin annotation, as well as biochemical and genetic experiments. They also establish that spanins, like viral membrane fusion proteins, adopt different strategies to achieve fusion of the inner and outer membranes.

Differing responses of red abalone (Haliotis rufescens) and white abalone (H. sorenseni) to infection with phage-associated Candidatus Xenohaliotis californiensis

The Rickettsiales-like prokaryote and causative agent of Withering Syndrome (WS)-Candidatus Xenohaliotis californiensis (Ca. Xc)-decimated black abalone populations along the Pacific coast of North America. White abalone-Haliotis sorenseni-are also susceptible to WS and have become nearly extinct in the wild due to overfishing in the 1970s. Candidatus Xenohaliotis californiensis proliferates within epithelial cells of the abalone gastrointestinal tract and causes clinical signs of starvation. In 2012, evidence of a putative bacteriophage associated with Ca. Xc in red abalone-Haliotis rufescens-was described. Recently, histologic examination of animals with Ca. Xc infection in California abalone populations universally appear to have the phage-containing inclusions. In this study, we investigated the current virulence of Ca. Xc in red abalone and white abalone at different environmental temperatures. Using a comparative experimental design, we observed differences over time between the two abalone species in mortality, body condition, and bacterial load by quantitative real time PCR (qPCR). By day 251, all white abalone exposed to the current variant of Ca. Xc held in the warm water (18.5 °C) treatment died, while red abalone exposed to the same conditions had a mortality rate of only 10%, despite a relatively heavy bacterial burden as determined by qPCR of posterior esophagus tissue and histological assessment at the termination of the experiment. These data support the current status of Ca. Xc as less virulent in red abalone, and may provide correlative evidence of a protective phage interaction. However, white abalone appear to remain highly susceptible to this disease. These findings have important implications for implementation of a white abalone recovery program, particularly with respect to the thermal regimes of locations where captively-reared individuals will be outplanted.

Pseudomonas PB1-Like Phages: Whole Genomes from Metagenomes Offer Insight into an Abundant Group of Bacteriophages

Despite the abundance, ubiquity and impact of environmental viruses, their inherent genomic plasticity and extreme diversity pose significant challenges for the examination of bacteriophages on Earth. Viral metagenomic studies have offered insight into broader aspects of phage ecology and repeatedly uncover genes to which we are currently unable to assign function. A combined effort of phage isolation and metagenomic survey of Chicago’s nearshore waters of Lake Michigan revealed the presence of Pbunaviruses, relatives of the Pseudomonas phage PB1. This prompted our expansive investigation of PB1-like phages. Genomic signatures of PB1-like phages and Pbunaviruses were identified, permitting the unambiguous distinction between the presence/absence of these phages in soils, freshwater and wastewater samples, as well as publicly available viral metagenomic datasets. This bioinformatic analysis led to the de novo assembly of nine novel PB1-like phage genomes from a metagenomic survey of samples collected from Lake Michigan. While this study finds that Pbunaviruses are abundant in various environments of Northern Illinois, genomic variation also exists to a considerable extent within individual communities.

Multiple mechanisms drive phage infection efficiency in nearly identical hosts

Phage-host interactions are critical to ecology, evolution, and biotechnology. Central to those is infection efficiency, which remains poorly understood, particularly in nature. Here we apply genome-wide transcriptomics and proteomics to investigate infection efficiency in nature's own experiment: two nearly identical (genetically and physiologically) Bacteroidetes bacterial strains (host18 and host38) that are genetically intractable, but environmentally important, where phage infection efficiency varies. On host18, specialist phage phi18:3 infects efficiently, whereas generalist phi38:1 infects inefficiently. On host38, only phi38:1 infects, and efficiently. Overall, phi18:3 globally repressed host18's transcriptome and proteome, expressed genes that likely evaded host restriction/modification (R/M) defenses and controlled its metabolism, and synchronized phage transcription with translation. In contrast, phi38:1 failed to repress host18's transcriptome and proteome, did not evade host R/M defenses or express genes for metabolism control, did not synchronize transcripts with proteins and its protein abundances were likely targeted by host proteases. However, on host38, phi38:1 globally repressed host transcriptome and proteome, synchronized phage transcription with translation, and infected host38 efficiently. Together these findings reveal multiple infection inefficiencies. While this contrasts the single mechanisms often revealed in laboratory mutant studies, it likely better reflects the phage-host interaction dynamics that occur in nature.

Namib Desert Soil Microbial Community Diversity, Assembly, and Function Along a Natural Xeric Gradient

The hyperarid Namib desert is a coastal desert in southwestern Africa and one of the oldest and driest deserts on the planet. It is characterized by a west/east increasing precipitation gradient and by regular coastal fog events (extending up to 75 km inland) that can also provide soil moisture. In this study, we evaluated the role of this natural aridity and xeric gradient on edaphic microbial community structure and function in the Namib desert. A total of 80 individual soil samples were collected at 10-km intervals along a 190-km transect from the fog-dominated western coastal region to the eastern desert boundary. Seventeen physicochemical parameters were measured for each soil sample. Soil parameters reflected the three a priori defined climatic/xeric zones along the transect ("fog," "low rain," and "high rain"). Microbial community structures were characterized by terminal restriction fragment length polymorphism fingerprinting and shotgun metaviromics, and their functional capacities were determined by extracellular enzyme activity assays. Both microbial community structures and activities differed significantly between the three xeric zones. The deep sequencing of surface soil metavirome libraries also showed shifts in viral composition along the xeric transect. While bacterial community assembly was influenced by soil chemistry and stochasticity along the transect, variations in community "function" were apparently tuned by xeric stress.

New insight into the RNA aquatic virosphere via viromics

RNA viruses that infect microbes are now recognized as an active, persistent and important component of the aquatic microbial community. While some information about the diversity and dynamics of the RNA virioplankton has been derived from culture-based and single gene approaches, research based on viromic and metatransciptomic methods has generated unprecedented insight into this relatively understudied class of microbes. Here, the relevant literature is summarized and discussed, including viromic studies of extracellular aquatic RNA viral assemblages, and transcriptomic studies of active and associated RNA viruses from aquatic environments followed by commentary on the present challenges and future directions of this field of research.

Viruses in Soil Ecosystems: An Unknown Quantity Within an Unexplored Territory

Viral abundance in soils can range from below detection limits in hot deserts to over 1 billion per gram in wetlands. Abundance appears to be strongly influenced by water availability and temperature, but a lack of informational standards creates difficulties for cross-study analysis. Soil viral diversity is severely underestimated and undersampled, although current measures of viral richness are higher for soils than for aquatic ecosystems. Both morphometric and metagenomic analyses have raised questions about the prevalence of nontailed, ssDNA viruses in soils. Soil is complex and critically important to terrestrial biodiversity and human civilization, but impacts of viral activities on soil ecosystem services are poorly understood. While information from aquatic systems and medical microbiology suggests the potential for viral influences on nutrient cycles, food web interactions, gene transfer, and other key processes in soils, very few empirical data are available. To understand the soil virome, much work remains.

Bacteriophage evolution differs by host, lifestyle and genome

Bacteriophages play key roles in microbial evolution^1,2, marine nutrient cycling³ and human disease⁴. Phages are genetically diverse, and their genome architectures are characteristically mosaic, driven by horizontal gene transfer with other phages and host genomes⁵. As a consequence, phage evolution is complex and their genomes are composed of genes with distinct and varied evolutionary histories^6,7. However, there are conflicting perspectives on the roles of mosaicism and the extent to which it generates a spectrum of genome diversity⁸ or genetically discrete populations^9,10. Here, we show that bacteriophages evolve within two general evolutionary modes that differ in the extent of horizontal gene transfer by an order of magnitude. Temperate phages distribute into high and low gene flux modes, whereas lytic phages share only the lower gene flux mode. The evolutionary modes are also a function of the bacterial host and different proportions of temperate and lytic phages are distributed in either mode depending on the host phylum. Groups of genetically related phages fall into either the high or low gene flux modes, suggesting there are genetic as well as ecological drivers of horizontal gene transfer rates. Consequently, genome mosaicism varies depending on the host, lifestyle and genetic constitution of phages.