Bacteriophages (or phages) play major roles in the evolution of bacterial pathogens via horizontal gene transfer. Multiple phages are often integrated in a host chromosome as prophages, not only carrying various novel virulence-related genetic determinants into host bacteria but also providing various possibilities for prophage-prophage interactions in bacterial cells. In particular, Escherichia coli strains such as Shiga toxin (Stx)-producing E. coli (STEC) and enteropathogenic E. coli (EPEC) strains have acquired more than 10 prophages (up to 21 prophages), many of which encode type III secretion system (T3SS) effector gene clusters. In these strains, some prophages are present at a single locus in tandem, which is usually interpreted as the integration of phages that use the same attachment (att) sequence. Here, we present phages integrating into T3SS effector gene cluster-associated loci in prophages, which are widely distributed in STEC and EPEC. Some of the phages integrated into prophages are Stx-encoding phages (Stx phages) and have induced the duplication of Stx phages in a single cell. The identified attB sequences in prophage genomes are apparently derived from host chromosomes. In addition, two or three different attB sequences are present in some prophages, which results in the generation of prophage clusters in various complex configurations. These phages integrating into prophages represent a medically and biologically important type of inter-phage interaction that promotes the accumulation of T3SS effector genes in STEC and EPEC, the duplication of Stx phages in STEC, and the conversion of EPEC to STEC and that may be distributed in other types of E. coli strains as well as other prophage-rich bacterial species.

Temperate phages (prophages) are ubiquitous in nature and persist as dormant components of host cells (lysogenic stage) before activating and lysing the host (lytic stage). Actively replicating prophages contribute to central community processes, such as enabling bacterial virulence, manipulating biogeochemical cycling, and driving microbial community diversification. Recent advances in sequencing technology have allowed for the identification and characterization of diverse phages, yet no approaches currently exist for identifying if a prophage has activated. Here, we present PropagAtE (Prophage Activity Estimator), an automated software tool for estimating if a prophage is in the lytic or lysogenic stage of infection. PropagAtE uses statistical analyses of prophage-to-host read coverage ratios to decipher actively replicating prophages, irrespective of whether prophages were induced or spontaneously activated. We demonstrate that PropagAtE is fast, accurate, and sensitive, regardless of sequencing depth. Application of PropagAtE to prophages from 348 complex metagenomes from human gut, murine gut, and soil environments identified distinct spatial and temporal prophage activation signatures, with the highest proportion of active prophages in murine gut samples. In infants treated with antibiotics or infants without treatment, we identified active prophage populations correlated with specific treatment groups. Within time series samples from the human gut, 11 prophage populations, some encoding the sulfur metabolism gene cysH or a rhuM-like virulence factor, were consistently present over time but not active. Overall, PropagAtE will facilitate accurate representations of viruses in microbiomes by associating prophages with their active roles in shaping microbial communities in nature. IMPORTANCE Viruses that infect bacteria are key components of microbiomes and ecosystems. They can kill and manipulate microorganisms, drive planetary-scale processes and biogeochemical cycling, and influence the structures of entire food networks. Prophages are viruses that can exist in a dormant state within the genome of their host (lysogenic stage) before activating in order to replicate and kill the host (lytic stage). Recent advances have allowed for the identification of diverse viruses in nature, but no approaches exist for characterizing prophages and their stages of infection (prophage activity). We develop and benchmark an automated approach, PropagAtE, to identify the stages of infection of prophages from genomic data. We provide evidence that active prophages vary in identity and abundance across multiple environments and scales. Our approach will enable accurate and unbiased analyses of viruses in microbiomes and ecosystems.

Prophages (viral genomes integrated within a host bacterial genome) can confer various phenotypic traits to their hosts, such as enhanced pathogenicity. Here we analyse >1300 genomes of 70 different Streptococcus species and identify nearly 800 prophages and satellite prophages (prophages that do not encode their own structural components but rely on the bacterial host and another helper prophage for survival). We show that prophages and satellite prophages are widely distributed among streptococci in a structured manner, and constitute two distinct entities with little effective genetic exchange between them. Cross-species transmission of prophages is not uncommon. Furthermore, a satellite prophage is associated with virulence in a mouse model of Streptococcus pneumoniae infection. Our findings highlight the potential importance of prophages in streptococcal biology and pathogenesis.

Prophages, which enables bacterial hosts to acquire novel traits, and increase genetic variation and evolutionary innovation, are considered to be one of the greatest drivers of bacterial diversity and evolution. Stenotrophomonas maltophilia is widely distributed and one of the most important multidrug resistant bacteria in hospitals. However, the distribution and genetic diversity of S. maltophilia prophages have not been elucidated. In this study, putative prophages were predicted in S. maltophilia genomes by using virus prediction tools, and the genetic diversity and phylogeny of S. maltophilia and the prophages they harbor were further analyzed. A total of 356 prophage regions were predicted from 88 S. maltophilia genomes. Among them, 144 were intact prophages, but 77.09% of the intact prophages did not match any known phage sequences in the public database. The number of prophage carried by S. maltophilia is related to its host habitat and is an important factor affecting the size of the host genome, but it is not related to the genetic diversity of the prophage. The prediction of auxiliary genes encoded by prophage showed that antibiotic resistance genes was not predicted for any of the prophages except for one questionable prophage, while 53 virulence genes and 169 carbohydrate active enzymes were predicted from 11.24 and 44.1% prophages, respectively. Most of the prophages (72.29%) mediated horizontal gene transfer of S. maltophilia genome, but only involved in 6.25% of the horizontal gene transfer events. In addition, CRISPR prediction indicated 97.75% S. maltophilia strains contained the CRISPR-Cas system containing 818 spacer sequences. However, these spacer sequences did not match any known S. maltophilia phages, and only a few S. maltophilia prophages. Comparative genomic analysis revealed a highly conserved and syntenic organization with genomic rearrangement between the prophages and the known related S. maltophilia phages. Our results indicate a high prevalence and genetic diversity of prophages in the genome of S. maltophilia, as well as the presence of a large number of uncharacterized phages. It provides an important complement to understanding the diversity and biological characteristics of phages, as well as the interactions and evolution between bacteria and phages.

Phages dominate every ecosystem on the planet. While virulent phages sculpt the microbiome by killing their bacterial hosts, temperate phages provide unique growth advantages to their hosts through lysogenic conversion. Many prophages benefit their host, and prophages are responsible for genotypic and phenotypic differences that separate individual microbial strains. However, the microbes also endure a cost to maintain those phages: additional DNA to replicate and proteins to transcribe and translate. We have never quantified those benefits and costs. Here, we analysed over two and a half million prophages from over half a million bacterial genome assemblies. Analysis of the whole dataset and a representative subset of taxonomically diverse bacterial genomes demonstrated that the normalised prophage density was uniform across all bacterial genomes above 2 Mbp. We identified a constant carrying capacity of phage DNA per bacterial DNA. We estimated that each prophage provides cellular services equivalent to approximately 2.4 % of the cell's energy or 0.9 ATP per bp per hour. We demonstrate analytical, taxonomic, geographic, and temporal disparities in identifying prophages in bacterial genomes that provide novel targets for identifying new phages. We anticipate that the benefits bacteria accrue from the presence of prophages balance the energetics involved in supporting prophages. Furthermore, our data will provide a new framework for identifying phages in environmental datasets, diverse bacterial phyla, and from different locations.

The bacterial archetypal adaptive immune system, CRISPR-Cas, is thought to be repressed in the best-studied bacterium, Escherichia coli K-12. We show here that the E. coli CRISPR-Cas system is active and serves to inhibit its nine defective (i.e., cryptic) prophages. Specifically, compared to the wild-type strain, reducing the amounts of specific interfering RNAs (crRNA) decreases growth by 40%, increases cell death by 700%, and prevents persister cell resuscitation. Similar results were obtained by inactivating CRISPR-Cas by deleting the entire 13 spacer region (CRISPR array); hence, CRISPR-Cas serves to inhibit the remaining deleterious effects of these cryptic prophages, most likely through CRISPR array-derived crRNA binding to cryptic prophage mRNA rather than through cleavage of cryptic prophage DNA, i.e., self-targeting. Consistently, four of the 13 E. coli spacers contain complementary regions to the mRNA sequences of seven cryptic prophages, and inactivation of CRISPR-Cas increases the level of mRNA for lysis protein YdfD of cryptic prophage Qin and lysis protein RzoD of cryptic prophage DLP-12. In addition, lysis is clearly seen via transmission electron microscopy when the whole CRISPR-Cas array is deleted, and eliminating spacer #12, which encodes crRNA with complementary regions for DLP-12 (including rzoD), Rac, Qin (including ydfD), and CP4-57 cryptic prophages, also results in growth inhibition and cell lysis. Therefore, we report the novel results that (i) CRISPR-Cas is active in E. coli and (ii) CRISPR-Cas is used to tame cryptic prophages, likely through RNAi, i.e., unlike with active lysogens, active CRISPR-Cas and cryptic prophages may stably co-exist.

Mobile genetic elements (MGEs) drive bacterial evolution, alter gene availability within microbial communities, and facilitate adaptation to ecological niches. In natural systems, bacteria simultaneously possess or encounter multiple MGEs, yet their combined influences on microbial communities are poorly understood. Here, we investigate interactions among MGEs in the marine bacterium Sulfitobacter pontiacus. Two related strains, CB-D and CB-A, each harbor a single prophage. These prophages share high sequence identity with one another and an integration site within the host genome, yet these strains exhibit differences in "spontaneous" prophage induction (SPI) and consequent fitness. To better understand mechanisms underlying variation in SPI between these lysogens, we closed their genomes, which revealed that in addition to harboring different prophage genotypes, CB-A lacks two of the four large, low-copy-number plasmids possessed by CB-D. To assess the relative roles of plasmid content versus prophage genotype on host physiology, a panel of derivative strains varying in MGE content were generated. Characterization of these derivatives revealed a robust link between plasmid content and SPI, regardless of prophage genotype. Strains possessing all four plasmids had undetectable phage in cell-free lysates, while strains lacking either one plasmid (pSpoCB-1) or a combination of two plasmids (pSpoCB-2 and pSpoCB-4) produced high (>105 PFU/mL) phage titers. Homologous plasmid sequences were identified in related bacteria, and plasmid and phage genes were found to be widespread in Tara Oceans metagenomic data sets. This suggests that plasmid-dependent stabilization of prophages may be commonplace throughout the oceans. IMPORTANCE The consequences of prophage induction on the physiology of microbial populations are varied and include enhanced biofilm formation, conferral of virulence, and increased opportunity for horizontal gene transfer. These traits lead to competitive advantages for lysogenized bacteria and influence bacterial lifestyles in a variety of niches. However, biological controls of "spontaneous" prophage induction, the initiation of phage replication and phage-mediated cell lysis without an overt stressor, are not well understood. In this study, we observed a novel interaction between plasmids and prophages in the marine bacterium Sulfitobacter pontiacus. We found that loss of one or more distinct plasmids-which we show carry genes ubiquitous in the world's oceans-resulted in a marked increase in prophage induction within lysogenized strains. These results demonstrate cross talk between different mobile genetic elements and have implications for our understanding of the lysogenic-lytic switches of prophages found not only in marine environments, but throughout all ecosystems.

Bacteria of the genus Desulfovibrio belong to the group of Sulphate Reducing Bacteria (SRB). SRB generate significant liabilities in the petroleum industry, mainly due to their ability to microbiologically induce corrosion, biofilm formation and H2S production. Bacteriophages are an alternative control method for SRB, whose information for this group of bacteria however, is scarce. The present study developed a workflow for the identification of complete prophages in Desulfovibrio. Poly-lysogenesis was shown to be common in Desulfovibrio. In the 47 genomes analyzed 53 complete prophages were identified. These were classified within the order Caudovirales, with 69.82% belonging to the Myoviridade family. More than half the prophages identified have genes coding for lysozyme or holin. Four of the analyzed bacterial genomes present prophages with identity above 50% in the same strain, whose comparative analysis demonstrated the existence of colinearity between the sequences. Of the 17 closed bacterial genomes analyzed, 6 have the CRISPR-Cas system classified as inactive. The identification of bacterial poly-lysogeny, the proximity between the complete prophages and the possible inactivity of the CRISPR-Cas in closed bacterial genomes analyzed allowed the choice of poly-lysogenic strains with prophages belonging to the Myoviridae family for the isolation of prophages and testing of related strains for subsequent studies.

This report describes the morphological characterization and genome analysis of an induced prophage (PLg-TB25) from a dairy strain of Lactococcus garvieae. The phage belongs to the Siphoviridae family and its morphology is typical of other lactococcal phages. A general analysis of its genome did not reveal similarities with other lactococcal phage genomes, confirming its novelty. However, similarities were found between genes of its morphogenesis cluster and genes of Gram-positive bacteria, suggesting that this phage genome resulted from recombination events that took place in a heterogeneous microbial environment. An in silico search for other prophages in 16 L. garvieae genomes available in public databases, uncovered eight seemingly complete prophages in strains isolated from dairy and fish niches. Genome analyses of these prophages revealed three novel L. garvieae phages. The remaining prophages had homology to phages of Lactococcus lactis (P335 group) suggesting a close relationship between these lactococcal species. The similarity in GC content of L. garvieae prophages to the genomes of L. lactis phages further supports the hypothesis that these phages likely originated from the same ancestor.

Identifying active prophages is critical for studying coevolution of phage and bacteria, investigating phage physiology and biochemistry, and engineering designer phages for diverse applications. We present Prophage Hunter, a tool aimed at hunting for active prophages from whole genome assembly of bacteria. Combining sequence similarity-based matching and genetic features-based machine learning classification, we developed a novel scoring system that exhibits higher accuracy than current tools in predicting active prophages on the validation datasets. The option of skipping similarity matching is also available so that there's higher chance for novel phages to be discovered. Prophage Hunter provides a one-stop web service to extract prophage genomes from bacterial genomes, evaluate the activity of the prophages, identify phylogenetically related phages, and annotate the function of phage proteins. Prophage Hunter is freely available at https://pro-hunter.bgi.com/.

We previously showed that zonula occludens toxin (Zot) encoded by Campylobacter concisus zot (808T) gene has the potential to initiate inflammatory bowel disease. This Zot protein caused prolonged intestinal epithelial barrier damage, induced intestinal epithelial and macrophage production of tumor necrosis factor-α and enhanced the responses of macrophages to other microbes. In order to understand the potential virulence of Zot proteins in other Campylobacter species, in this study we examined their presence, similarities, motifs and prophages.

DNA of viral origin represents a ubiquitous element of bacterial genomes. Its integration into host regulatory circuits is a pivotal driver of microbial evolution but requires the stringent regulation of phage gene activity. In this study, we describe the nucleoid-associated protein CgpS, which represents an essential protein functioning as a xenogeneic silencer in the Gram-positive Corynebacterium glutamicum CgpS is encoded by the cryptic prophage CGP3 of the C. glutamicum strain ATCC 13032 and was first identified by DNA affinity chromatography using an early phage promoter of CGP3. Genome-wide profiling of CgpS binding using chromatin affinity purification and sequencing (ChAP-Seq) revealed its association with AT-rich DNA elements, including the entire CGP3 prophage region (187 kbp), as well as several other elements acquired by horizontal gene transfer. Countersilencing of CgpS resulted in a significantly increased induction frequency of the CGP3 prophage. In contrast, a strain lacking the CGP3 prophage was not affected and displayed stable growth. In a bioinformatics approach, cgpS orthologs were identified primarily in actinobacterial genomes as well as several phage and prophage genomes. Sequence analysis of 618 orthologous proteins revealed a strong conservation of the secondary structure, supporting an ancient function of these xenogeneic silencers in phage-host interaction.

Prophages have been considered genetic units that have an intimate association with novel phenotypic properties of bacterial hosts, such as pathogenicity and genomic variation. Little is known about the genetic information of prophages in the genome of Streptococcus mutans, a major pathogen of human dental caries. In this study, we identified 35 prophage-like elements in S. mutans genomes and performed a comparative genomic analysis. Comparative genomic and phylogenetic analyses of prophage sequences revealed that the prophages could be classified into three main large clusters: Cluster A, Cluster B, and Cluster C. The S. mutans prophages in each cluster were compared. The genomic sequences of phismuN66-1, phismuNLML9-1, and phismu24-1 all shared similarities with the previously reported S. mutans phages M102, M102AD, and ϕAPCM01. The genomes were organized into seven major gene clusters according to the putative functions of the predicted open reading frames: packaging and structural modules, integrase, host lysis modules, DNA replication/recombination modules, transcriptional regulatory modules, other protein modules, and hypothetical protein modules. Moreover, an integrase gene was only identified in phismuNLML9-1 prophages.

Lactic acid bacteria (LAB) are used as starters in a wide variety of food fermentations. While the number of reports of phages infecting other LAB steadily increased over the years, information about phage associated with Latilactobacillus sakei, a frequently used meat starter, remains scarce.

The prevalence of methicillin-resistant Staphylococcus aureus (MRSA) type SCCmec IV or V is increasing in Taiwan. It has been suggested that the surface protein SasX is responsible for their transmission. However, the sasX gene was not detected in our SCCmec IV or V isolates. Since sasX was originally found in S. epidermidis and believed to be transferred to S. aureus by a prophage, studies were conducted to detect and type this prophage in our clinical isolates.

Staphylococcus aureus infections are of growing concern given the increased incidence of antibiotic resistant strains. Egypt, like several other countries, has seen alarming increases in methicillin-resistant S. aureus (MRSA) infections. This species can rapidly acquire genes associated with resistance, as well as virulence factors, through mobile genetic elements, including phages. Recently, we sequenced 56 S. aureus genomes from Alexandria Main University Hospital in Alexandria, Egypt, complementing 17 S. aureus genomes publicly available from other sites in Egypt. In the current study, we found that the majority (73.6%) of these strains contain intact prophages, including Biseptimaviruses, Phietaviruses, and Triaviruses. Further investigation of these prophages revealed evidence of horizontal exchange of the integrase for two of the prophages. These Egyptian S. aureus prophages are predicted to encode numerous virulence factors, including genes associated with immune evasion and toxins, including the Panton-Valentine leukocidin (PVL)-associated genes lukF-PV/lukS-PV. Thus, prophages are likely to be a major contributor to the virulence of S. aureus strains in circulation in Egypt.

Oenococcus oeni is the most exploited lactic acid bacterium in the wine industry and drives the malolactic fermentation of wines. Although prophage-like sequences have been identified in the species, many are not characterized, and a global view of their integration and distribution amongst strains is currently lacking. In this work, we analyzed the complete genomes of 231 strains for the occurrence of prophages, and analyzed their size and positions of insertion. Our data show the limited variation in the number of prophages in O. oeni genomes, and that six sites of insertion within the bacterial genome are being used for site-specific recombination. Prophage diversity patterns varied significantly for different host lineages, and environmental niches. Overall, the findings highlight the pervasive presence of prophages in the O. oeni species, their role as a major source of within-species bacterial diversity and drivers of horizontal gene transfer. Our data also have implications for enhanced understanding of the prophage recombination events which occurred during evolution of O. oeni, as well as the potential of prophages in influencing the fitness of these bacteria in their distinct niches.

Burkholderia spp. are often resistant to antibiotics, and infections with these organisms are difficult to treat. A potential alternative treatment for Burkholderia spp. infections is bacteriophage (phage) therapy; however, it can be difficult to locate phages that target these bacteria. Prophages incorporated into the bacterial genome have been identified within Burkholderia spp. and may represent a source of useful phages for therapy. Here we investigate whether prophages within Burkholderia spp. clinical isolates can kill conspecific and heterospecific isolates. Thirty-two Burkholderia spp. isolates were induced for prophage release, and harvested prophages were tested for lytic activity against the same 32 isolates. Lytic phages were passaged and their host ranges were determined, resulting in four unique phages of prophage origin that showed different ranges of lytic activity. We also analyzed the prophage content of 35 Burkholderia spp. clinical isolate genomes, and identified several prophages present in the genomes of multiple isolates of the same species. Finally, we observed that B. cenocepacia isolates were more phage-susceptible than Burkholderia multivorans isolates. Overall, our findings suggest that prophages present within Burkholderia spp. genomes are a potentially useful starting point for the isolation and development of novel phages for use in phage therapy.

Klebsiella pneumoniae is an increasing threat to public health and represents one of the most concerning pathogens involved in life-threatening infections. The resistant and virulence determinants are coded by mobile genetic elements which can easily spread between bacteria populations and co-evolve with its genomic host. In this study, we present the full genomic sequences, insertion sites and phylogenetic analysis of 150 prophages found in 40 K. pneumoniae clinical isolates obtained from an outbreak in a Portuguese hospital. All strains harbored at least one prophage and we identified 104 intact prophages (69.3%). The prophage size ranges from 29.7 to 50.6 kbp, coding between 32 and 78 putative genes. The prophage GC content is 51.2%, lower than the average GC content of 57.1% in K. pneumoniae. Complete prophages were classified into three families in the order Caudolovirales: Myoviridae (59.6%), Siphoviridae (38.5%) and Podoviridae (1.9%). In addition, an alignment and phylogenetic analysis revealed nine distinct clusters. Evidence of recombination was detected within the genome of some prophages but, in most cases, proteins involved in viral structure, transcription, replication and regulation (lysogenic/lysis) were maintained. These results support the knowledge that prophages are diverse and widely disseminated in K. pneumoniae genomes, contributing to the evolution of this species and conferring additional phenotypes. Moreover, we identified K. pneumoniae prophages in a set of endolysin genes, which were found to code for proteins with lysozyme activity, cleaving the β-1,4 linkages between N-acetylmuramic acid and N-acetyl-D-glucosamine residues in the peptidoglycan network and thus representing genes with the potential for lysin phage therapy.

Burkholderia cenocepacia, is a Gram-negative opportunistic pathogen that belongs to Burkholderia cepacia complex (BCC) group. BCC representatives carry various pathogenicity factors and can infect humans and plants. Phages as bacterial viruses play a significant role in biodiversity and ecological balance in the environment. Specifically, horizontal gene transfer (HGT) and lysogenic conversion (temperate phages) influence microbial diversification and fitness. In this study, we describe the prevalence and gene content of prophages in 16 fully sequenced B. cenocepacia genomes stored in NCBI database. The analysis was conducted in silico by manual and automatic approaches. Sixty-three potential prophage regions were found and classified as intact, incomplete, questionable, and artifacts. The regions were investigated for the presence of known virulence factors, resulting in the location of sixteen potential pathogenicity mechanisms, including toxin⁻antitoxin systems (TA), Major Facilitator Superfamily (MFS) transporters and responsible for drug resistance. Investigation of the region's closest neighborhood highlighted three groups of genes with the highest occurrence-tRNA-Arg, dehydrogenase family proteins, and ABC transporter substrate-binding proteins. Searches for antiphage systems such as BacteRiophage EXclusion (BREX) and Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) in the analyzed strains suggested 10 sequence sets of CRISPR elements. Our results suggest that intact B. cenocepacia prophages may provide an evolutionary advantage to the bacterium, while domesticated prophages may help to maintain important genes.