Prophage-mediated defence against viral attack and viral counter-defence

(Ph)ighting Phages: How Bacteria Resist Their Parasites

Bacteria are under constant attack from bacteriophages (phages), bacterial parasites that are the most abundant biological entity on earth. To resist phage infection, bacteria have evolved an impressive arsenal of anti-phage systems. Recent advances have significantly broadened and deepened our understanding of how bacteria battle phages, spearheaded by new systems like CRISPR-Cas. This review aims to summarize bacterial anti-phage mechanisms, with an emphasis on the most recent developments in the field.

Complete Genome Sequence of Microbacterium foliorum NRRL B-24224, a Host for Bacteriophage Discovery

We report the complete annotated genome sequence of Microbacterium foliorum NRRL B-24224, a type strain isolated from the phyllosphere of grasses and a commonly used host for bacteriophage discovery. The genome contains no identifiable prophage or CRISPR or restriction-modification system, which suggests that it may continue to be a fruitful host for phage discovery.

We report the complete annotated genome sequence of Microbacterium foliorum NRRL B-24224, a type strain isolated from the phyllosphere of grasses and a commonly used host for bacteriophage discovery. The genome contains no identifiable prophage or CRISPR or restriction-modification system, which suggests that it may continue to be a fruitful host for phage discovery.

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.

Mycobacteriophage ZoeJ: A broad host-range close relative of mycobacteriophage TM4

A collection of over 1600 sequenced bacteriophages isolated on a single host strain, Mycobacterium smegmatis mc²155, can be grouped into over two dozen types that have little or no nucleotide sequence similarity to each other. One group, Cluster K, can be divided into several subclusters, and the well-characterized and much exploited phage TM4 lies in Subcluster K2. Many of the Cluster K phages have broad host ranges and infect both fast- and slow-growing mycobacterial strains. Here we describe phage ZoeJ, a new Subcluster K2 member, which infects a broad spectrum of mycobacterial hosts including M. smegmatis, Mycobacterium tuberculosis, and Mycobacterium avium. ZoeJ has extensive sequence similarity to TM4, and comparative analysis reveals the precise deletion conferring the lytic phenotype of TM4. The ZoeJ immunity repressor was identified as gene 45, which is prophage-expressed, is required for lysogeny, and is sufficient to confer superinfection immunity to ZoeJ. ZoeJ gp45 also confers immunity to Subcluster K2 phage Milly, and Subcluster K1 phages Adephagia and CrimD, but surprisingly not to TM4. RNAseq analysis reveals the temporal pattern of early and late gene expressions in ZoeJ lytic growth and suggests a role for the ESAS motifs for gene regulation.

Enterococcus faecalis Countermeasures Defeat a Virulent Picovirinae Bacteriophage

Enterococcus faecalis is an opportunistic pathogen that has emerged as a major cause of nosocomial infections worldwide. Many clinical strains are indeed resistant to last resort antibiotics and there is consequently a reawakening of interest in exploiting virulent phages to combat them. However, little is still known about phage receptors and phage resistance mechanisms in enterococci. We made use of a prophageless derivative of the well-known clinical strain E. faecalis V583 to isolate a virulent phage belonging to the Picovirinae subfamily and to the P68 genus that we named Idefix. Interestingly, most isolates of E. faecalis tested—including V583—were resistant to this phage and we investigated more deeply into phage resistance mechanisms. We found that E. faecalis V583 prophage 6 was particularly efficient in resisting Idefix infection thanks to a new abortive infection (Abi) mechanism, which we designated Abiα. It corresponded to the Pfam domain family with unknown function DUF4393 and conferred a typical Abi phenotype by causing a premature lysis of infected E. faecalis. The abiα gene is widespread among prophages of enterococci and other Gram-positive bacteria. Furthermore, we identified two genes involved in the synthesis of the side chains of the surface rhamnopolysaccharide that are important for Idefix adsorption. Interestingly, mutants in these genes arose at a frequency of ~10⁻⁴ resistant mutants per generation, conferring a supplemental bacterial line of defense against Idefix.

The Revisited Genome of Bacillus subtilis Bacteriophage SPP1

Bacillus subtilis bacteriophage SPP1 is a lytic siphovirus first described 50 years ago [1]. Its complete DNA sequence was reported in 1997 [2]. Here we present an updated annotation of the 44,016 bp SPP1 genome and its correlation to different steps of the viral multiplication process. Five early polycistronic transcriptional units encode phage DNA replication proteins and lysis functions together with less characterized, mostly non-^{essential, functions. Late transcription drives synthesis of proteins necessary for SPP1 viral particles assembly and for cell lysis, together with a short set of proteins of unknown function. The extensive genetic, biochemical and structural biology studies on the molecular mechanisms of SPP1 DNA replication and phage particle assembly rendered it a model system for tailed phages research. We propose SPP1} as the reference species for a new SPP1-like viruses genus of the Siphoviridae family.

Close Encounters of Three Kinds: Bacteriophages, Commensal Bacteria, and Host Immunity

Recent years have witnessed an explosion of interest in the human microbiota. Although commensal bacteria have dominated research efforts to date, mounting evidence suggests that endogenous viral populations (the 'virome') play key roles in basic human physiology. The most numerous constituents of the human virome are not eukaryotic viruses but rather bacteriophages, viruses that infect bacteria. Here, we review phages' interactions with their immediate (prokaryotic) and extended (eukaryotic) hosts and with each other, with a particular emphasis on the temperate phages and prophages which dominate the human virome. We also discuss key outstanding questions in this emerging field and emphasize the urgent need for functional studies in animal models to complement previous in vitro work and current computational approaches.

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.

Detection of Bacteriophage Particles Containing Antibiotic Resistance Genes in the Sputum of Cystic Fibrosis Patients

Cystic fibrosis (CF) is a chronic disease in which the bacterial colonization of the lung is linked to an excessive inflammatory response that leads to respiratory failure. The microbiology of CF is complex. Staphylococcus aureus is the first bacterium to colonize the lungs in 30% of pediatric CF patients, and 80% of adult patients develop a chronic Pseudomonas aeruginosa infection, but other microorganisms can also be found. The use of antibiotics is essential to treat the disease, but antibiotic performance is compromised by resistance mechanisms. Among various mechanisms of transfer of antibiotic resistance genes (ARGs), the recently been reported bacteriophages are the least explored in clinical settings. To determine the role of phages in CF as mobile genetic elements (MGEs) carrying ARGs, we evaluated their presence in 71 CF patients. 71 sputum samples taken from these patients were screened for eight ARGs (bla_TEM, bla_CTX-M-1-group, bla_CTX-M-9-group, bla_OXA-48, bla_VIM, mecA, qnrA, and qnrS) in the bacteriophage DNA fraction. The phages found were also purified and observed by electron microscopy. 32.4% of CF patients harbored ARGs in phage DNA. β-lactamase genes, particularly bla_VIM and bla_TEM, were the most prevalent and abundant, whereas mecA, qnrA, and qnrS were very rare. Siphoviridae phage particles capable of infecting P. aeruginosa and Klebsiella pneumoniae were detected in CF sputum. Phage particles harboring ARGs were found to be abundant in the lungs of both CF patients and healthy individuals and could contribute to the colonization of multiresistant strains.

Mycobacteriophage Fruitloop gp52 inactivates Wag31 (DivIVA) to prevent heterotypic superinfection

Bacteriophages engage in complex dynamic interactions with their bacterial hosts and with each other. Bacteria have numerous mechanisms to resist phage infection, and phages must co-evolve by overcoming bacterial resistance or by choosing an alternative host. Phages also compete with each other, both during lysogeny by prophage-mediated defense against viral attack and by superinfection exclusion during lytic replication. Phages are enormously diverse genetically and are replete with small genes of unknown function, many of which are not required for lytic growth, but which may modulate these bacteria-phage and phage-phage dynamics. Using cellular toxicity of phage gene overexpression as an assay, we identified the 93-residue protein gp52 encoded by Cluster F mycobacteriophage Fruitloop. The toxicity of Fruitloop gp52 overexpression results from interaction with and inactivation of Wag31 (DivIVA), an essential Mycobacterium smegmatis protein organizing cell wall biosynthesis at the growing cellular poles. Fruitloop gene 52 is expressed early in lytic growth and is not required for normal Fruitloop lytic replication but interferes with Subcluster B2 phages such as Hedgerow and Rosebush. We conclude that Hedgerow and Rosebush are Wag31-dependent phages and that Fruitloop gp52 confers heterotypic superinfection exclusion by inactivating Wag31.

Bacteriophages of the Urinary Microbiome

Bacterial viruses (bacteriophages) play a significant role in microbial community dynamics. Within the human gastrointestinal tract, for instance, associations among bacteriophages (phages), microbiota stability, and human health have been discovered. In contrast to the gastrointestinal tract, the phages associated with the urinary microbiota are largely unknown. Preliminary metagenomic surveys of the urinary virome indicate a rich diversity of novel lytic phage sequences at an abundance far outnumbering that of eukaryotic viruses. These surveys, however, exclude the lysogenic phages residing within the bacteria of the bladder. To characterize this phage population, we examined 181 genomes representative of the phylogenetic diversity of bacterial species within the female urinary microbiota and found 457 phage sequences, 226 of which were predicted with high confidence. Phages were prevalent within the bladder bacteria: 86% of the genomes examined contained at least one phage sequence. Most of these phages are novel, exhibiting no discernible sequence homology to sequences in public data repositories. The presence of phages with substantial sequence similarity within the microbiota of different women supports the existence of a core community of phages within the bladder. Furthermore, the observed variation between the phage populations of women with and without overactive bladder symptoms suggests that phages may contribute to urinary health. To complement our bioinformatic analyses, viable phages were cultivated from the bacterial isolates for characterization; a novel coliphage was isolated, which is obligately lytic in the laboratory strain Escherichia coli C. Sequencing of bacterial genomes facilitates a comprehensive cataloguing of the urinary virome and reveals phage-host interactions.IMPORTANCE Bacteriophages are abundant within the human body. However, while some niches have been well surveyed, the phage population within the urinary microbiome is largely unknown. Our study is the first survey of the lysogenic phage population within the urinary microbiota. Most notably, the abundance of prophage exceeds that of the bacteria. Furthermore, many of the prophage sequences identified exhibited no recognizable sequence homology to sequences in data repositories. This suggests a rich diversity of uncharacterized phage species present in the bladder. Additionally, we observed a variation in the abundances of phages between bacteria isolated from asymptomatic "healthy" individuals and those with urinary symptoms, thus suggesting that, like phages within the gut, phages within the bladder may contribute to urinary health.

An inclusive Research Education Community (iREC): Impact of the SEA-PHAGES program on research outcomes and student learning

The Science Education Alliance–Phage Hunters Advancing Genomics and Evolutionary Science program is an inclusive Research Education Community with centralized programmatic and scientific support, in which broad student engagement in authentic science is linked to increased accessibility to research experiences for students; increased persistence of these students in science, technology, engineering, and mathematics; and increased scientific productivity for students and faculty alike.

Engaging undergraduate students in scientific research promises substantial benefits, but it is not accessible to all students and is rarely implemented early in college education, when it will have the greatest impact. An inclusive Research Education Community (iREC) provides a centralized scientific and administrative infrastructure enabling engagement of large numbers of students at different types of institutions. The Science Education Alliance–Phage Hunters Advancing Genomics and Evolutionary Science (SEA-PHAGES) is an iREC that promotes engagement and continued involvement in science among beginning undergraduate students. The SEA-PHAGES students show strong gains correlated with persistence relative to those in traditional laboratory courses regardless of academic, ethnic, gender, and socioeconomic profiles. This persistent involvement in science is reflected in key measures, including project ownership, scientific community values, science identity, and scientific networking.

Expression and evolutionary patterns of mycobacteriophage D29 and its temperate close relatives

Background

Mycobacteriophages are viruses that infect Mycobacterium hosts. A large collection of phages known to infect the same bacterial host strain – Mycobacterium smegmatis mc²155 – exhibit substantial diversity and characteristically mosaic architectures. The well-studied lytic mycobacteriophage D29 appears to be a deletion derivative of a putative temperate parent, although its parent has yet to be identified.

Results

Here we describe three newly-isolated temperate phages – Kerberos, Pomar16 and StarStuff – that are related to D29, and are predicted to be very close relatives of its putative temperate parent, revealing the repressor and additional genes that are lost in D29. Transcriptional profiles show the patterns of both lysogenic and lytic gene expression and identify highly-expressed, abundant, stable, small non-coding transcripts made from the P_left early lytic promoter, and which are toxic to M. smegmatis.

Conclusions

Comparative genomics of phages D29, Kerberos, Pomar16 and StarStuff provide insights into bacteriophage evolution, and comparative transcriptomics identifies the pattern of lysogenic and lytic expression with unusual features including highly expressed, small, non-coding RNAs.

Electronic supplementary material

The online version of this article (10.1186/s12866-017-1131-2) contains supplementary material, which is available to authorized users.

Virus-host protein-protein interactions of mycobacteriophage Giles

Mycobacteriophage are viruses that infect mycobacteria. More than 1,400 mycobacteriophage genomes have been sequenced, coding for over one hundred thousand proteins of unknown functions. Here we investigate mycobacteriophage Giles-host protein-protein interactions (PPIs) using yeast two-hybrid screening (Y2H). A total of 25 reproducible PPIs were found for a selected set of 10 Giles proteins, including a putative virion assembly protein (gp17), the phage integrase (gp29), the endolysin (gp31), the phage repressor (gp47), and six proteins of unknown function (gp34, gp35, gp54, gp56, gp64, and gp65). We note that overexpression of the proteins is toxic to M. smegmatis, although whether this toxicity and the associated changes in cellular morphology are related to the putative interactions revealed in the Y2H screen is unclear.

Genome Sequences of Chancellor, Mitti, and Wintermute, Three Subcluster K4 Phages Isolated Using Mycobacterium smegmatis mc 2 155

Mycobacteriophages Chancellor, Mitti, and Wintermute infect Mycobacterium smegmatis mc²155 and are closely related to phages Cheetobro and Fionnbharth in subcluster K4. Genome sizes range from 57,697 bp to 58,046 bp. Phages are predicted to be temperate and to infect the pathogen Mycobacterium tuberculosis.

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.

Small Alarmone Synthetases as novel bacterial RNA-binding proteins

The alarmone nucleotides guanosine pentaphosphate (pppGpp) and tetraphosphate (ppGpp), collectively referred to as (p)ppGpp, are key regulators of bacterial growth, stress adaptation, antibiotic tolerance and pathogenicity. We have recently shown that the Small Alarmone Synthetase (SAS) RelQ from the Gram-positive pathogen Enterococcus faecalis has an RNA-binding activity (Beljantseva et al. 2017). RelQ's activities as an enzyme and as an RNA-binding protein are mutually incompatible: binding of single-stranded RNA potently inhibits (p)ppGpp synthesis in a sequence-specific manner, and RelQ's enzymatic activity destabilizes the RNA:RelQ complex. RelQ's allosteric regulator, pppGpp, destabilizes RNA binding and activates RelQ's enzymatic activity. Since SAS enzymes are widely distributed in bacteria, and, as has been discovered recently, are also mobilized by phages (Dedrick et al. 2017), RNA binding to SASs could be a widespread mechanism. The initial discovery raises numerous questions regarding RNA-binding function of the SAS enzymes: What is the molecular mechanism underlying the incompatibility of RNA:SAS complex formation with pppGpp binding and (p)ppGpp synthesis? What are the RNA targets in living cells? What is the regulatory output of the system - (p)ppGpp synthesis, modulation of RNA structure and function, or both?

Viral communities of the human gut: metagenomic analysis of composition and dynamics

Background

The numerically most abundant biological entities on Earth are viruses. Vast populations prey on the cellular microbiota in all habitats, including the human gut.

Main body

Here we review approaches for studying the human virome, and some recent results on movement of viral sequences between bacterial cells and eukaryotic hosts. We first overview biochemical and bioinformatic methods, emphasizing that specific choices in the methods used can have strong effects on the results obtained. We then review studies characterizing the virome of the healthy human gut, which reveal that most of the viruses detected are typically uncharacterized phage - the viral dark matter - and that viruses that infect human cells are encountered only rarely. We then review movement of phage between bacterial cells during antibiotic treatment. Here a radical proposal for extensive movement of antibiotic genes on phage has been challenged by a careful reanalysis of the metagenomic annotation methods used. We then review two recent studies of movement of whole phage communities between human individuals during fecal microbial transplantation, which emphasize the possible role of lysogeny in dispersal.

Short conclusion

Methods for studying the human gut virome are improving, yielding interesting data on movement of phage genes between cells and mammalian host organisms. However, viral populations are vast, and studies of their composition and function are just beginning.

The Discovery, Mechanisms, and Evolutionary Impact of Anti-CRISPRs

Bacteria and archaea use CRISPR-Cas adaptive immune systems to defend themselves from infection by bacteriophages (phages). These RNA-guided nucleases are powerful weapons in the fight against foreign DNA, such as phages and plasmids, as well as a revolutionary gene editing tool. Phages are not passive bystanders in their interactions with CRISPR-Cas systems, however; recent discoveries have described phage genes that inhibit CRISPR-Cas function. More than 20 protein families, previously of unknown function, have been ascribed anti-CRISPR function. Here, we discuss how these CRISPR-Cas inhibitors were discovered and their modes of action were elucidated. We also consider the potential impact of anti-CRISPRs on bacterial and phage evolution. Finally, we speculate about the future of this field.

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.

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.

Embracing the enemy: the diversification of microbial gene repertoires by phage-mediated horizontal gene transfer

Bacteriophages and archaeal viruses contribute, through lysogenic conversion or transduction, to the horizontal transfer of genetic material between microbial genomes. Recent genomics, metagenomics, and single cell studies have shown that lysogenic conversion is widespread and provides hosts with adaptive traits often associated with biotic interactions. The quantification of the evolutionary impact of transduction has lagged behind and requires further theoretical and experimental work. Nevertheless, recent studies suggested that generalized transduction plays a role in the transfer of antibiotic resistance genes and in the acquisition of novel genes during intra-specific bacterial competition. The characteristics of transduction and lysogenic conversion complement those of other mechanisms of transfer, and could play a key role in the spread of adaptive genes between communities.

Variability and host density independence in inductions-based estimates of environmental lysogeny

Temperate bacterial viruses (phages) may enter a symbiosis with their host cell, forming a unit called a lysogen. Infection and viral replication are disassociated in lysogens until an induction event such as DNA damage occurs, triggering viral-mediated lysis. The lysogen-lytic viral reproduction switch is central to viral ecology, with diverse ecosystem impacts. It has been argued that lysogeny is favoured in phages at low host densities. This paradigm is based on the fraction of chemically inducible cells (FCIC) lysogeny proxy determined using DNA-damaging mitomycin C inductions. Contrary to the established paradigm, a survey of 39 inductions publications found FCIC to be highly variable and pervasively insensitive to bacterial host density at global, within-environment and within-study levels. Attempts to determine the source(s) of variability highlighted the inherent complications in using the FCIC proxy in mixed communities, including dissociation between rates of lysogeny and FCIC values. Ultimately, FCIC studies do not provide robust measures of lysogeny or consistent evidence of either positive or negative host density dependence to the lytic-lysogenic switch. Other metrics are therefore needed to understand the drivers of the lytic-lysogenic decision in viral communities and to test models of the host density-dependent viral lytic-lysogenic switch.