Bacteriophage evolution differs by host, lifestyle and genome

Evolution of Superinfection Immunity in Cluster A Mycobacteriophages

Many aspects regarding superinfection, immunity, virulence, and the evolution of immune specificities are poorly understood due to the lack of large collections of isolated and sequenced phages with a spectrum of genetic diversity. Using a genetically diverse collection of Cluster A phages, we show that the classical and relatively straightforward patterns of homoimmunity, heteroimmunity, and virulence result from interactions between homotypic and heterotypic phages at the extreme edges of an evolutionary continuum of immune specificities. Genetic interactions between mesotypic phages result in more complex mesoimmunity phenotypes and virulence profiles. These results highlight that the evolution of immune specificities can be shaped by homotypic and mesotypic interactions and may be more dynamic than previously considered.

Temperate phages encode an immunity system to control lytic gene expression during lysogeny. This gene regulatory circuit consists of multiple interacting genetic elements, and although it is essential for controlling phage growth, it is subject to conflicting evolutionary pressures. During superinfection of a lysogen, the prophage’s circuit interacts with the superinfecting phage’s circuit and prevents lytic growth if the two circuits are closely related. The circuitry is advantageous since it provides the prophage with a defense mechanism, but the circuitry is also disadvantageous since it limits the phage’s host range during superinfection. Evolutionarily related phages have divergent, orthogonal immunity systems that no longer interact and are heteroimmune, but we do not understand how immunity systems evolve new specificities. Here, we use a group of Cluster A mycobacteriophages that exhibit a spectrum of genetic diversity to examine how immunity system evolution impacts superinfection immunity. We show that phages with mesotypic (i.e., genetically related but distinct) immunity systems exhibit asymmetric and incomplete superinfection phenotypes. They form complex immunity networks instead of well-defined immunity groups, and mutations conferring escape (i.e., virulence) from homotypic or mesotypic immunity have various escape specificities. Thus, virulence and the evolution of new immune specificities are shaped by interactions with homotypic and mesotypic immunity systems.

Marine DNA Viral Macro- and Microdiversity from Pole to Pole

Microbes drive most ecosystems and are modulated by viruses that impact their lifespan, gene flow, and metabolic outputs. However, ecosystem-level impacts of viral community diversity remain difficult to assess due to classification issues and few reference genomes. Here, we establish an ∼12-fold expanded global ocean DNA virome dataset of 195,728 viral populations, now including the Arctic Ocean, and validate that these populations form discrete genotypic clusters. Meta-community analyses revealed five ecological zones throughout the global ocean, including two distinct Arctic regions. Across the zones, local and global patterns and drivers in viral community diversity were established for both macrodiversity (inter-population diversity) and microdiversity (intra-population genetic variation). These patterns sometimes, but not always, paralleled those from macro-organisms and revealed temperate and tropical surface waters and the Arctic as biodiversity hotspots and mechanistic hypotheses to explain them. Such further understanding of ocean viruses is critical for broader inclusion in ecosystem models.

Staphylococci phages display vast genomic diversity and evolutionary relationships

Background

Bacteriophages are the most abundant and diverse entities in the biosphere, and this diversity is driven by constant predator–prey evolutionary dynamics and horizontal gene transfer. Phage genome sequences are under-sampled and therefore present an untapped and uncharacterized source of genetic diversity, typically characterized by highly mosaic genomes and no universal genes. To better understand the diversity and relationships among phages infecting human pathogens, we have analysed the complete genome sequences of 205 phages of Staphylococcus sp.

Results

These are predicted to encode 20,579 proteins, which can be sorted into 2139 phamilies (phams) of related sequences; 745 of these are orphams and possess only a single gene. Based on shared gene content, these phages were grouped into four clusters (A, B, C and D), 27 subclusters (A1-A2, B1-B17, C1-C6 and D1-D2) and one singleton. However, the genomes have mosaic architectures and individual genes with common ancestors are positioned in distinct genomic contexts in different clusters. The staphylococcal Cluster B siphoviridae are predicted to be temperate, and the integration cassettes are often closely-linked to genes implicated in bacterial virulence determinants. There are four unusual endolysin organization strategies found in Staphylococcus phage genomes, with endolysins predicted to be encoded as single genes, two genes spliced, two genes adjacent and as a single gene with inter-lytic-domain secondary translational start site. Comparison of the endolysins reveals multi-domain modularity, with conservation of the SH3 cell wall binding domain.

Conclusions

This study provides a high-resolution view of staphylococcal viral genetic diversity, and insights into their gene flux patterns within and across different phage groups (cluster and subclusters) providing insights into their evolution.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5647-8) contains supplementary material, which is available to authorized users.

Taxonomic assignment of uncultivated prokaryotic virus genomes is enabled by gene-sharing networks

Microbiomes from every environment contain a myriad of uncultivated archaeal and bacterial viruses, but studying these viruses is hampered by the lack of a universal, scalable taxonomic framework. We present vConTACT v.2.0, a network-based application utilizing whole genome gene-sharing profiles for virus taxonomy that integrates distance-based hierarchical clustering and confidence scores for all taxonomic predictions. We report near-identical (96%) replication of existing genus-level viral taxonomy assignments from the International Committee on Taxonomy of Viruses for National Center for Biotechnology Information virus RefSeq. Application of vConTACT v.2.0 to 1,364 previously unclassified viruses deposited in virus RefSeq as reference genomes produced automatic, high-confidence genus assignments for 820 of the 1,364. We applied vConTACT v.2.0 to analyze 15,280 Global Ocean Virome genome fragments and were able to provide taxonomic assignments for 31% of these data, which shows that our algorithm is scalable to very large metagenomic datasets. Our taxonomy tool can be automated and applied to metagenomes from any environment for virus classification.

The Third Annual Meeting of the European Virus Bioinformatics Center

The Third Annual Meeting of the European Virus Bioinformatics Center (EVBC) took place in Glasgow, United Kingdom, 28-29 March 2019. Virus bioinformatics has become central to virology research, and advances in bioinformatics have led to improved approaches to investigate viral infections and outbreaks, being successfully used to detect, control, and treat infections of humans and animals. This active field of research has attracted approximately 110 experts in virology and bioinformatics/computational biology from Europe and other parts of the world to attend the two-day meeting in Glasgow to increase scientific exchange between laboratory- and computer-based researchers. The meeting was held at the McIntyre Building of the University of Glasgow; a perfect location, as it was originally built to be a place for "rubbing your brains with those of other people", as Rector Stanley Baldwin described it. The goal of the meeting was to provide a meaningful and interactive scientific environment to promote discussion and collaboration and to inspire and suggest new research directions and questions. The meeting featured eight invited and twelve contributed talks, on the four main topics: (1) systems virology, (2) virus-host interactions and the virome, (3) virus classification and evolution and (4) epidemiology, surveillance and evolution. Further, the meeting featured 34 oral poster presentations, all of which focused on specific areas of virus bioinformatics. This report summarizes the main research findings and highlights presented at the meeting.

Rates of Mutation and Recombination in Siphoviridae Phage Genome Evolution over Three Decades

The evolution of asexual organisms is driven not only by the inheritance of genetic modification but also by the acquisition of foreign DNA. The contribution of vertical and horizontal processes to genome evolution depends on their rates per year and is quantified by the ratio of recombination to mutation. These rates have been estimated for bacteria; however, no estimates have been reported for phages. Here, we delineate the contribution of mutation and recombination to dsDNA phage genome evolution. We analyzed 34 isolates of the 936 group of Siphoviridae phages using a Lactococcus lactis strain from a single dairy over 29 years. We estimate a constant substitution rate of 1.9 × 10⁻⁴ substitutions per site per year due to mutation that is within the range of estimates for eukaryotic RNA and DNA viruses. The reconstruction of recombination events reveals a constant rate of five recombination events per year and 4.5 × 10⁻³ nucleotide alterations due to recombination per site per year. Thus, the recombination rate exceeds the substitution rate, resulting in a relative effect of recombination to mutation (r/m) of ∼24 that is homogenous over time. Especially in the early transcriptional region, we detect frequent gene loss and regain due to recombination with phages of the 936 group, demonstrating the role of the 936 group pangenome as a reservoir of genetic variation. The observed substitution rate homogeneity conforms to the neutral theory of evolution; hence, the neutral theory can be applied to phage genome evolution and also to genetic variation brought about by recombination.

Mycobacteriophage CRB2 defines a new subcluster in mycobacteriophage classification

Mycobacteriophages are viruses -mostly temperates- that infect Mycobacterium smegmatis and sometimes Mycobacterium tuberculosis. Mycobacteriophages are grouped in clusters on the basis of the overall nucleotide sequence homology, being further divided in subclusters as more mycobacteriophage genomes are sequenced and annotated. As part of our on-going screening for novel isolates, we herein report the bioinformatics analysis of CRB2, a mycobacteriophage belonging into the Siphoviridae family that propagates at 30°C. CRB2 has a 72,217 bp genome with a 69.78% GC content that belongs to Cluster B; nucleotide comparison with other B cluster members positions CRB2 as the sole member of a new subcluster, B9, being mycobacteriophage Saguaro (belonging into subcluster B7) its closest relative. Sequencing and annotation of 14 mycobacteriophages isolated by our group has yielded six cluster A members, a singleton, four of the five members of subcluster B6, one of the three reported members of subcluster G4, and CRB2 which defines subcluster B9. Considering the massive mycobacteriophage search performed in USA and the relatively rarity of our phages, we propose that factors other than size of the sampling determine the variability of mycobacteriophage distribution, and thus a world-wide concerted mining would most likely bring extremely rare and yet undiscovered mycobacteriophages.

Mycobacteriophages

Mycobacteriophages are viruses that infect mycobacterial hosts. A large number of mycobacteriophages have been isolated and genomically characterized, providing insights into viral diversity and evolution, as well as fueling development of tools for mycobacterial genetics. Mycobacteriophages have intimate relationships with their hosts and provide insights into the genetics and physiology of the mycobacteria and tools for potential clinical applications such as drug development, diagnosis, vaccines, and potentially therapy.

Global In-Silico Scenario of tRNA Genes and Their Organization in Virus Genomes

Viruses are known to be highly dependent on the host translation machinery for their protein synthesis. However, tRNA genes are occasionally identified in such organisms, and in addition, few of them harbor tRNA gene clusters comprising dozens of genes. Recently, tRNA gene clusters have been shown to occur among the three domains of life. In such a scenario, the viruses could play a role in the dispersion of such structures among these organisms. Thus, in order to reveal the prevalence of tRNA genes as well as tRNA gene clusters in viruses, we performed an unbiased large-scale genome survey. Interestingly, tRNA genes were predicted in ssDNA (single-stranded DNA) and ssRNA (single-stranded RNA) viruses as well in many other dsDNA viruses of families from Caudovirales order. In the latter group, tRNA gene clusters composed of 15 to 37 tRNA genes were characterized, mainly in bacteriophages, enlarging the occurrence of such structures within viruses. These bacteriophages were from hosts that encompass five phyla and 34 genera. This in-silico study presents the current global scenario of tRNA genes and their organization in virus genomes, contributing and opening questions to be explored in further studies concerning the role of the translation apparatus in these organisms.

Contemporary Phage Biology: From Classic Models to New Insights

Bacteriophages, discovered about a century ago, have been pivotal as models for understanding the fundamental principles of molecular biology. While interest in phage biology declined after the phage "golden era," key recent developments, including advances in phage genomics, microscopy, and the discovery of the CRISPR-Cas anti-phage defense system, have sparked a renaissance in phage research in the past decade. This review highlights recently discovered unexpected complexities in phage biology, describes a new arsenal of phage genes that help them overcome bacterial defenses, and discusses advances toward documentation of the phage biodiversity on a global scale.

Minimum Information about an Uncultivated Virus Genome (MIUViG)

We present an extension of the Minimum Information about any (x) Sequence (MIxS) standard for reporting sequences of uncultivated virus genomes. Minimum Information about an Uncultivated Virus Genome (MIUViG) standards were developed within the Genomic Standards Consortium framework and include virus origin, genome quality, genome annotation, taxonomic classification, biogeographic distribution and in silico host prediction. Community-wide adoption of MIUViG standards, which complement the Minimum Information about a Single Amplified Genome (MISAG) and Metagenome-Assembled Genome (MIMAG) standards for uncultivated bacteria and archaea, will improve the reporting of uncultivated virus genomes in public databases. In turn, this should enable more robust comparative studies and a systematic exploration of the global virosphere.

Biology and Taxonomy of crAss-like Bacteriophages, the Most Abundant Virus in the Human Gut

CrAssphages represent the most abundant virus in the human gut microbiota, but the lack of available genome sequences for comparison has kept them enigmatic. Recently, sequence-based classification of distantly related crAss-like phages from multiple environments was reported, leading to a proposed familial-level taxonomic group. Here, we assembled the metagenomic sequencing reads from 702 human fecal virome/phageome samples and analyzed 99 complete circular crAss-like phage genomes and 150 contigs ≥70 kb. In silico comparative genomics and taxonomic analysis enabled a classification scheme of crAss-like phages from human fecal microbiomes into four candidate subfamilies composed of ten candidate genera. Laboratory analysis was performed on fecal samples from an individual harboring seven distinct crAss-like phages. We achieved crAss-like phage propagation in ex vivo human fecal fermentations and visualized short-tailed podoviruses by electron microscopy. Mass spectrometry of a crAss-like phage capsid protein could be linked to metagenomic sequencing data, confirming crAss-like phage structural annotations.

Evaluation of the genomic diversity of viruses infecting bacteria, archaea and eukaryotes using a common bioinformatic platform: steps towards a unified taxonomy

Genome Relationship Applied to Virus Taxonomy (GRAViTy) is a genetics-based tool that computes sequence relatedness between viruses. Composite generalized Jaccard (CGJ) distances combine measures of homology between encoded viral genes and similarities in genome organizational features (gene orders and orientations). This scoring framework effectively recapitulates the current, largely morphology and phenotypic-based, family-level classification of eukaryotic viruses. Eukaryotic virus families typically formed monophyletic groups with consistent CGJ distance cut-off dividing between and within family divergence ranges. In the current study, a parallel analysis of prokaryotic virus families revealed quite different sequence relationships, particularly those of tailed phage families (Siphoviridae, Myoviridae and Podoviridae), where members of the same family were generally far more divergent and often not detectably homologous to each other. Analysis of the 20 currently classified prokaryotic virus families indeed split them into 70 separate clusters of tailed phages genetically equivalent to family-level assignments of eukaryotic viruses. It further divided several bacterial (Sphaerolipoviridae, Tectiviridae) and archaeal (Lipothrixviridae) families. We also found that the subfamily-level groupings of tailed phages were generally more consistent with the family assignments of eukaryotic viruses, and this supports ongoing reclassifications, including Spounavirinae and Vi1virus taxa as new virus families. The current study applied a common benchmark with which to compare taxonomies of eukaryotic and prokaryotic viruses. The findings support the planned shift away from traditional morphology-based classifications of prokaryotic viruses towards a genome-based taxonomy. They demonstrate the feasibility of a unified taxonomy of viruses into which the vast body of metagenomic viral sequences may be consistently assigned.

Characterization and induction of prophages in human gut-associated Bifidobacterium hosts

In the current report, we describe the identification of three genetically distinct groups of prophages integrated into three different chromosomal sites of human gut-associated Bifidobacterium breve and Bifidobacterium longum strains. These bifidobacterial prophages are distantly related to temperate actinobacteriophages of several hosts. Some prophages, integrated within the dnaJ₂ gene, are competent for induction, excision, replication, assembly and lysis, suggesting that they are fully functional and can generate infectious particles, even though permissive hosts have not yet been identified. Interestingly, several of these phages harbor a putative phase variation shufflon (the Rin system) that generates variation of the tail-associated receptor binding protein (RBP). Unlike the analogous coliphage-associated shufflon Min, or simpler Cin and Gin inversion systems, Rin is predicted to use a tyrosine recombinase to promote inversion, the first reported phage-encoded tyrosine-family DNA invertase. The identification of bifidobacterial prophages with RBP diversification systems that are competent for assembly and lysis, yet fail to propagate lytically under laboratory conditions, suggests dynamic evolution of bifidobacteria and their phages in the human gut.

Novel phages of healthy skin metaviromes from South Africa

Recent skin metagenomic studies have investigated the harbored viral diversity and its possible influence on healthy skin microbial populations, and tried to establish global patterns of skin-phage evolution. However, the detail associated with the phages that potentially play a role in skin health has not been investigated. While skin metagenome and -metavirome studies have indicated that the skin virome is highly site specific and shows marked interpersonal variation, they have not assessed the presence/absence of individual phages. Here, we took a semi-culture independent approach (metaviromic) to better understand the composition of phage communities on skin from South African study participants. Our data set adds over 130 new phage species of the skin to existing databases. We demonstrated that identical phages were present on different individuals and in different body sites, and we conducted a detailed analysis of the structural organization of these phages. We further found that a bacteriophage related to the Staphylococcus capitis phage Stb20 may be a common skin commensal virus potentially regulating its host and its activities on the skin.

Biological species in the viral world

Due to their dependence on cellular organisms for metabolism and replication, viruses are typically named and assigned to species according to their genome structure and the original host that they infect. But because viruses often infect multiple hosts and the numbers of distinct lineages within a host can be vast, their delineation into species is often dictated by arbitrary sequence thresholds, which are highly inconsistent across lineages. Here we apply an approach to determine the boundaries of viral species based on the detection of gene flow within populations, thereby defining viral species according to the biological species concept (BSC). Despite the potential for gene transfer between highly divergent genomes, viruses, like the cellular organisms they infect, assort into reproductively isolated groups and can be organized into biological species. This approach revealed that BSC-defined viral species are often congruent with the taxonomic partitioning based on shared gene contents and host tropism, and that bacteriophages can similarly be classified in biological species. These results open the possibility to use a single, universal definition of species that is applicable across cellular and acellular lifeforms.

Evolutionary History of Bacteriophages in the Genus Paraburkholderia

The genus Paraburkholderia encompasses mostly environmental isolates with diverse predicted lifestyles. Genome analyses have shown that bacteriophages form a considerable portion of some Paraburkholderia genomes. Here, we analyzed the evolutionary history of prophages across all Paraburkholderia spp. Specifically, we investigated to what extent the presence of prophages and their distribution affect the diversity/diversification of Paraburkholderia spp., as well as to what extent phages coevolved with their respective hosts. Particular attention was given to the presence of CRISPR-Cas arrays as a reflection of past interactions with phages. We thus analyzed 36 genomes of Paraburkholderia spp., including those of 11 new strains, next to those of three Burkholderia species. Most genomes were found to contain at least one full prophage sequence. The highest number was found in Paraburkholderia sp. strain MF2-27; the nine prophages found amount to up to 4% of its genome. Among all prophages, potential moron genes (e.g., DNA adenine methylase) were found that might be advantageous for host cell fitness. Co-phylogenetic analyses indicated the existence of complex evolutionary scenarios between the different Paraburkholderia hosts and their prophages, including short-term co-speciation, duplication, host-switching and phage loss events. Analysis of the CRISPR-Cas systems showed a record of diverse, potentially recent, phage infections. We conclude that, overall, different phages have interacted in diverse ways with their Paraburkholderia hosts over evolutionary time.

A Decade of Streptococcus thermophilus Phage Evolution in an Irish Dairy Plant

Phages of Streptococcus thermophilus present a major threat to the production of many fermented dairy products. To date, only a few studies have assessed the biodiversity of S. thermophilus phages in dairy fermentations. In order to develop strategies to limit phage predation in this important industrial environment, it is imperative that such studies are undertaken and that phage-host interactions of this species are better defined. The present study investigated the biodiversity and evolution of phages within an Irish dairy fermentation facility over an 11-year period. This resulted in the isolation of 17 genetically distinct phages, all of which belong to the so-called cos group. The evolution of phages within the factory appears to be influenced by phages from other dairy plants introduced into the factory for whey protein powder production. Modular exchange, primarily within the regions encoding lysogeny and replication functions, was the major observation among the phages isolated between 2006 and 2016. Furthermore, the genotype of the first isolate in 2006 was observed continuously across the following decade, highlighting the ability of these phages to prevail in the factory setting for extended periods of time. The proteins responsible for host recognition were analyzed, and carbohydrate-binding domains (CBDs) were identified in the distal tail (Dit), the baseplate proteins, and the Tail-associated lysin (Tal) variable regions (VR1 and VR2) of many isolates. This supports the notion that S. thermophilus phages recognize a carbohydrate receptor on the cell surface of their host.IMPORTANCE Dairy fermentations are consistently threatened by the presence of bacterial viruses (bacteriophages or phages), which may lead to a reduction in acidification rates or even complete loss of the fermentate. These phages may persist in factories for long periods of time. The objective of the current study was to monitor the progression of phages infecting the dairy bacterium Streptococcus thermophilus over a period of 11 years in an Irish dairy plant so as to understand how these phages evolve. A focused analysis of the genomic region that encodes host recognition functions highlighted that the associated proteins harbor a variety of carbohydrate-binding domains, which corroborates the notion that phages of S. thermophilus recognize carbohydrate receptors at the initial stages of the phage cycle.

The lytic Myoviridae of Enterobacteriaceae form tight recombining assemblages separated by discontinuities in genome average nucleotide identity and lateral gene flow

In Bacteria, a working consensus of species circumscription may have been reached and one of the most prominent criteria is high average nucleotide identity (ANI). ANI in effect groups strains that may recombine more or less frequently, depending on their biology, as opposed to rare interspecies gene transfer. For bacteriophages, which show various lifestyles, the nature of the fundamental natural unit, if it exists in a biological sense, is not well understood and defined. The approaches based on dot-plots are useful to group similar bacteriophages, yet are not quantitative and use arbitrarily set cut-offs. Here, we focus on lytic Myoviridae and test the ANI metric for group delineation. We show that ANI-based groups are in agreement with the International Committee on Taxonomy of Viruses (ICTV) classification and already established dot-plot groups, which are occasionally further refined owing to higher resolution of ANI analysis. Furthermore, these groups are separated among themselves by clear ANI discontinuities. Their members readily exchange core genes with each other while they do not with bacteriophages of other ANI groups, not even with the most similar. Thus, ANI-delineated groups may represent the natural units in lytic Myoviridae evolution with features that resemble those encountered in bacterial species.

Ecological and Evolutionary Benefits of Temperate Phage: What Does or Doesn't Kill You Makes You Stronger

Infection by a temperate phage can lead to death of the bacterial cell, but sometimes these phages integrate into the bacterial chromosome, offering the potential for a more long-lasting relationship to be established. Here we define three major ecological and evolutionary benefits of temperate phage for bacteria: as agents of horizontal gene transfer (HGT), as sources of genetic variation for evolutionary innovation, and as weapons of bacterial competition. We suggest that a coevolutionary perspective is required to understand the roles of temperate phages in bacterial populations.

Bacteriophages of Gordonia spp. Display a Spectrum of Diversity and Genetic Relationships

The global bacteriophage population is large, dynamic, old, and highly diverse genetically. Many phages are tailed and contain double-stranded DNA, but these remain poorly characterized genomically. A collection of over 1,000 phages infecting Mycobacterium smegmatis reveals the diversity of phages of a common bacterial host, but their relationships to phages of phylogenetically proximal hosts are not known. Comparative sequence analysis of 79 phages isolated on Gordonia shows these also to be diverse and that the phages can be grouped into 14 clusters of related genomes, with an additional 14 phages that are “singletons” with no closely related genomes. One group of six phages is closely related to Cluster A mycobacteriophages, but the other Gordonia phages are distant relatives and share only 10% of their genes with the mycobacteriophages. The Gordonia phage genomes vary in genome length (17.1 to 103.4 kb), percentage of GC content (47 to 68.8%), and genome architecture and contain a variety of features not seen in other phage genomes. Like the mycobacteriophages, the highly mosaic Gordonia phages demonstrate a spectrum of genetic relationships. We show this is a general property of bacteriophages and suggest that any barriers to genetic exchange are soft and readily violable.

Despite the numerical dominance of bacteriophages in the biosphere, there is a dearth of complete genomic sequences. Current genomic information reveals that phages are highly diverse genomically and have mosaic architectures formed by extensive horizontal genetic exchange. Comparative analysis of 79 phages of Gordonia shows them to not only be highly diverse, but to present a spectrum of relatedness. Most are distantly related to phages of the phylogenetically proximal host Mycobacterium smegmatis, although one group of Gordonia phages is more closely related to mycobacteriophages than to the other Gordonia phages. Phage genome sequence space remains largely unexplored, but further isolation and genomic comparison of phages targeted at related groups of hosts promise to reveal pathways of bacteriophage evolution.