Advances in phage-host interaction prediction: in silico method enhances the development of phage therapies.

Phages can specifically recognize and kill bacteria, which lead to important application value of bacteriophage in bacterial identification and typing, livestock aquaculture and treatment of human bacterial infection. Considering the variety of human-infected bacteria and the continuous discovery of numerous pathogenic bacteria, screening suitable therapeutic phages that are capable of infecting pathogens from massive phage databases has been a principal step in phage therapy design. Experimental methods to identify phage-host interaction (PHI) are time-consuming and expensive; high-throughput computational method to predict PHI is therefore a potential substitute. Here, we systemically review bioinformatic methods for predicting PHI, introduce reference databases and in silico models applied in these methods and highlight the strengths and challenges of current tools. Finally, we discuss the application scope and future research direction of computational prediction methods, which contribute to the performance improvement of prediction models and the development of personalized phage therapy.

Phage tailspike modularity and horizontal gene transfer reveals specificity towards E. coli O-antigen serogroups.

BACKGROUND: The interaction between bacteriophages and their hosts is intricate and highly specific. Receptor-binding proteins (RBPs) of phages such as tail fibers and tailspikes initiate the infection process. These RBPs bind to diverse outer membrane structures, including the O-antigen, which is a serogroup-specific sugar-based component of the outer lipopolysaccharide layer of Gram-negative bacteria. Among the most virulent Escherichia coli strains is the Shiga toxin-producing E. coli (STEC) pathotype dominated by a subset of O-antigen serogroups. METHODS: Extensive phylogenetic and structural analyses were used to identify and validate specificity correlations between phage RBP subtypes and STEC O-antigen serogroups, relying on the principle of horizontal gene transfer as main driver for RBP evolution. RESULTS: We identified O-antigen specific RBP subtypes for seven out of nine most prevalent STEC serogroups (O26, O45, O103, O104, O111, O145 and O157) and seven additional E. coli serogroups (O2, O8, O16, O18, 4s/O22, O77 and O78). Eight phage genera (Gamaleya-, Justusliebig-, Kaguna-, Kayfuna-, Kutter-, Lederberg-, Nouzilly- and Uetakeviruses) emerged for their high proportion of serogroup-specific RBPs. Additionally, we reveal sequence motifs in the RBP region, potentially serving as recombination hotspots between lytic phages. CONCLUSION: The results contribute to a better understanding of mosaicism of phage RBPs, but also demonstrate a method to identify and validate new RBP subtypes for current and future emerging serogroups.

Understanding Bacteriophage Tail Fiber Interaction with Host Surface Receptor: The Key "Blueprint" for Reprogramming Phage Host Range.

Bacteriophages (phages), as natural antibacterial agents, are being rediscovered because of the growing threat of multi- and pan-drug-resistant bacterial pathogens globally. However, with an estimated 1031 phages on the planet, finding the right phage to recognize a specific bacterial host is like looking for a needle in a trillion haystacks. The host range of a phage is primarily determined by phage tail fibers (or spikes), which initially mediate reversible and specific recognition and adsorption by susceptible bacteria. Recent significant advances at single-molecule and atomic levels have begun to unravel the structural organization of tail fibers and underlying mechanisms of phage-host interactions. Here, we discuss the molecular mechanisms and models of the tail fibers of the well-characterized T4 phage's interaction with host surface receptors. Structure-function knowledge of tail fibers will pave the way for reprogramming phage host range and will bring future benefits through more-effective phage therapy in medicine. Furthermore, the design strategies of tail fiber engineering are briefly summarized, including machine-learning-assisted engineering inspired by the increasingly enormous amount of phage genetic information.

Introducing differential RNA-seq mapping to track the early infection phase for Pseudomonas phage ɸKZ.

As part of the ongoing renaissance of phage biology, more phage genomes are becoming available through DNA sequencing. However, our understanding of the transcriptome architecture that allows these genomes to be expressed during host infection is generally poor. Transcription start sites (TSSs) and operons have been mapped for very few phages, and an annotated global RNA map of a phage - alone or together with its infected host - is not available at all. Here, we applied differential RNA-seq (dRNA-seq) to study the early, host takeover phase of infection by assessing the transcriptome structure of Pseudomonas aeruginosa jumbo phage ɸKZ, a model phage for viral genetics and structural research. This map substantially expands the number of early expressed viral genes, defining TSSs that are active ten minutes after ɸKZ infection. Simultaneously, we record gene expression changes in the host transcriptome during this critical metabolism conversion. In addition to previously reported upregulation of genes associated with amino acid metabolism, we observe strong activation of genes with functions in biofilm formation (cdrAB) and iron storage (bfrB), as well as an activation of the antitoxin ParD. Conversely, ɸKZ infection rapidly down-regulates complexes IV and V of oxidative phosphorylation (atpCDGHF and cyoABCDE). Taken together, our data provide new insights into the transcriptional organization and infection process of the giant bacteriophage ɸKZ and adds a framework for the genome-wide transcriptomic analysis of phage-host interactions.

A Specific Sugar Moiety in the Lactococcus lactis Cell Wall Pellicle Is Required for Infection by CHPC971, a Member of the Rare 1706 Phage Species.

Lactococcus lactis is a Gram-positive bacterium widely used as a starter culture for the production of different dairy products, especially a large variety of cheeses. Infection of lactococcal starter cultures by bacteriophages is one of the major causes of fermentation failure and often leads to production halt. Lactococcal bacteriophages belonging to the c2, 936, and P335 species are the most commonly isolated in dairy plants and have been extensively investigated in the past three decades. Information regarding bacteriophages belonging to less commonly isolated species is, on the other hand, less extensive, although these phages can also contribute to starter culture infection. Here, we report the nucleotide sequence of the newly isolated L. lactis phage CHPC971, belonging to the rare 1706 species of lactococcal phages. We investigated the nature of the host receptor recognized by the phage and collected evidence that strongly suggests that it binds to a specific sugar moiety in the cell wall pellicle of its host. An in silico analysis of the genome of phage CHPC971 identified the hypothetical genes involved in receptor binding.IMPORTANCE Gathering information on how lactococcal bacteriophages recognize their host and proliferate in the dairy environment is of vital importance for the establishment of proper starter culture rotation plans and to avoid fermentation failure and consequent great economic losses for dairy industries. We provide strong evidence on the type of receptor recognized by a newly isolated 1706-type lactococcal bacteriophage, increasing knowledge of phage-host interactions relevant to dairying. This information can help to prevent phage infection events that, so far, are hard to predict and avoid.

Long-read based de novo assembly of low-complexity metagenome samples results in finished genomes and reveals insights into strain diversity and an active phage system.

BACKGROUND: Complete and contiguous genome assemblies greatly improve the quality of subsequent systems-wide functional profiling studies and the ability to gain novel biological insights. While a de novo genome assembly of an isolated bacterial strain is in most cases straightforward, more informative data about co-existing bacteria as well as synergistic and antagonistic effects can be obtained from a direct analysis of microbial communities. However, the complexity of metagenomic samples represents a major challenge. While third generation sequencing technologies have been suggested to enable finished metagenome-assembled genomes, to our knowledge, the complete genome assembly of all dominant strains in a microbiome sample has not been demonstrated. Natural whey starter cultures (NWCs) are used in cheese production and represent low-complexity microbiomes. Previous studies of Swiss Gruyère and selected Italian hard cheeses, mostly based on amplicon metagenomics, concurred that three species generally pre-dominate: Streptococcus thermophilus, Lactobacillus helveticus and Lactobacillus delbrueckii. RESULTS: Two NWCs from Swiss Gruyère producers were subjected to whole metagenome shotgun sequencing using the Pacific Biosciences Sequel and Illumina MiSeq platforms. In addition, longer Oxford Nanopore Technologies MinION reads had to be generated for one to resolve repeat regions. Thereby, we achieved the complete assembly of all dominant bacterial genomes from these low-complexity NWCs, which was corroborated by a 16S rRNA amplicon survey. Moreover, two distinct L. helveticus strains were successfully co-assembled from the same sample. Besides bacterial chromosomes, we could also assemble several bacterial plasmids and phages and a corresponding prophage. Biologically relevant insights were uncovered by linking the plasmids and phages to their respective host genomes using DNA methylation motifs on the plasmids and by matching prokaryotic CRISPR spacers with the corresponding protospacers on the phages. These results could only be achieved by employing long-read sequencing data able to span intragenomic as well as intergenomic repeats. CONCLUSIONS: Here, we demonstrate the feasibility of complete de novo genome assembly of all dominant strains from low-complexity NWCs based on whole metagenomics shotgun sequencing data. This allowed to gain novel biological insights and is a fundamental basis for subsequent systems-wide omics analyses, functional profiling and phenotype to genotype analysis of specific microbial communities.

Viral interference of the bacterial RNA metabolism machinery.

In a recent publication, we reported a unique interaction between a protein encoded by the giant myovirus phiKZ and the Pseudomonas aeruginosa RNA degradosome. Crystallography, site-directed mutagenesis and interactomics approaches revealed this 'degradosome interacting protein' or Dip, to adopt an 'open-claw' dimeric structure that presents acidic patches on its outer surface which hijack 2 conserved RNA binding sites on the scaffold domain of the RNase E component of the RNA degradosome. This interaction prevents substrate RNAs from being bound and degraded by the RNA degradosome during the virus infection cycle. In this commentary, we provide a perspective into the biological role of Dip, its structural analysis and its mysterious evolutionary origin, and we suggest some therapeutic and biotechnological applications of this distinctive viral protein.

A phage protein that inhibits the bacterial ATPase required for type IV pilus assembly.

Type IV pili (TFPs) are required for bacterial twitching motility and for phage infection in the opportunistic human pathogen Pseudomonas aeruginosa. Here we describe a phage-encoded protein, D3112 protein gp05 (hereafter referred to as Tip, representing twitching inhibitory protein), whose expression is necessary and sufficient to mediate the inhibition of twitching motility. Tip interacts with and blocks the activity of bacterial-encoded PilB, the TFP assembly/extension ATPase, at an internal 40-aa region unique to PilB. Tip expression results in the loss of surface piliation. Based on these observations and the fact that many P. aeruginosa phages require TFPs for infection, Tip-mediated twitching inhibition may represent a generalized strategy for superinfection exclusion. Moreover, because TFPs are required for full virulence, PilB may be an attractive target for the development of novel antiinfectives.

DNA motifs determining the efficiency of adaptation into the Escherichia coli CRISPR array.

Clustered regularly interspaced short palindromic repeats (CRISPR) and their associated proteins constitute a recently identified prokaryotic defense system against invading nucleic acids. DNA segments, termed protospacers, are integrated into the CRISPR array in a process called adaptation. Here, we establish a PCR-based assay that enables evaluating the adaptation efficiency of specific spacers into the type I-E Escherichia coli CRISPR array. Using this assay, we provide direct evidence that the protospacer adjacent motif along with the first base of the protospacer (5'-AAG) partially affect the efficiency of spacer acquisition. Remarkably, we identified a unique dinucleotide, 5'-AA, positioned at the 3' end of the spacer, that enhances efficiency of the spacer's acquisition. Insertion of this dinucleotide increased acquisition efficiency of two different spacers. DNA sequencing of newly adapted CRISPR arrays revealed that the position of the newly identified motif with respect to the 5'-AAG is important for affecting acquisition efficiency. Analysis of approximately 1 million spacers showed that this motif is overrepresented in frequently acquired spacers compared with those acquired rarely. Our results represent an example of a short nonprotospacer adjacent motif sequence that affects acquisition efficiency and suggest that other as yet unknown motifs affect acquisition efficiency in other CRISPR systems as well.

Phage SEP1 hijacks Staphylococcus epidermidis stationary cells' metabolism to replicate.

In nature, bacteria often survive in a stationary state with low metabolic activity. Phages use the metabolic machinery of the host cell to replicate, and, therefore, their efficacy against non-dividing cells is usually limited. Nevertheless, it was previously shown that the Staphylococcus epidermidis phage SEP1 has the remarkable capacity to actively replicate in stationary-phase cells, reducing their numbers. Here, we studied for the first time the transcriptomic profiles of both exponential and stationary cells infected with SEP1 phage using RNA-seq to gain a better understanding of this rare phenomenon. We showed that SEP1 successfully takes over the transcriptional apparatus of both exponential and stationary cells. Infection was, however, delayed in stationary cells, with genes within the gp142-gp154 module putatively implicated in host takeover. S. epidermidis responded to SEP1 infection by upregulating three genes involved in a DNA modification system, with this being observed already 5 min after infection in exponential cells and later in stationary cells. In stationary cells, a significant number of genes involved in translation and RNA metabolic and biosynthetic processes were upregulated after 15 and 30 min of SEP1 infection in comparison with the uninfected control, showing that SEP1 activates metabolic and biosynthetic pathways necessary to its successful replication.IMPORTANCEMost phage-host interaction studies are performed with exponentially growing cells. However, this cell state is not representative of what happens in natural environments. Additionally, most phages fail to replicate in stationary cells. The Staphylococcus epidermidis phage SEP1 is one of the few phages reported to date to be able to infect stationary cells. Here, we unveiled the interaction of SEP1 with its host in both exponential and stationary states of growth at the transcriptomic level. The findings of this study provide valuable insights for a better implementation of phage therapy since phages able to infect stationary cells could be more efficient in the treatment of recalcitrant infections.

Unique relationships between phages and endospore-forming hosts.

As part of their survival strategy under harsh environmental conditions, endospore-forming bacteria can trigger a sporulation developmental program. Although the regulatory cascades that precisely control the transformation of vegetative bacteria into mother cells and resilient spores have been described in detail, less is known about how bacteriophages that prey on endospore-formers exploit sporulation. Herein, we argue that phages infecting these bacteria have evolved several specific molecular mechanisms, not yet known in other bacteria, that manifest from the phage-driven alliance to negative effects on the host. We anticipate that the relationships between phages and endospore-formers outlined here will inspire studies on phage ecology and evolution, and could facilitate important advances in the development of phage therapies against pathogenic spore-formers.

Antiphage small molecules produced by bacteria - beyond protein-mediated defenses.

Bacterial populations face the constant threat of viral predation exerted by bacteriophages ('phages'). In response, bacteria have evolved a wide range of defense mechanisms against phage challenges. Yet the vast majority of antiphage defense systems described until now are mediated by proteins or RNA complexes acting at the single-cell level. Here, we review small molecule-based defense strategies against phage infection, with a focus on the antiphage molecules described recently. Importantly, inhibition of phage infection by excreted small molecules has the potential to protect entire bacterial communities, highlighting the ecological significance of these antiphage strategies. Considering the immense repertoire of bacterial metabolites, we envision that the list of antiphage small molecules will be further expanded in the future.

Transcriptional Landscapes of Herelleviridae Bacteriophages and Staphylococcus aureus during Phage Infection: An Overview.

The issue of antibiotic resistance in healthcare worldwide has led to a pressing need to explore and develop alternative approaches to combat infectious diseases. Among these methods, phage therapy has emerged as a potential solution to tackle this growing challenge. Virulent phages of the Herelleviridae family, known for their ability to cause lysis of Staphylococcus aureus, a clinically significant pathogen frequently associated with multidrug resistance, have proven to be one of the most effective viruses utilized in phage therapy. In order to utilize phages for therapeutic purposes effectively, a thorough investigation into their physiology and mechanisms of action on infected cells is essential. The use of omics technologies, particularly total RNA sequencing, is a promising approach for analyzing the interaction between phages and their hosts, allowing for the assessment of both the behavior of the phage during infection and the cell's response. This review aims to provide a comprehensive overview of the physiology of the Herelleviridae family, utilizing existing analyses of their total phage transcriptomes. Additionally, it sheds light on the changes that occur in the metabolism of S. aureus when infected with virulent bacteriophages, contributing to a deeper understanding of the phage-host interaction.

An ensemble method for prediction of phage-based therapy against bacterial infections.

Phage therapy is a viable alternative to antibiotics for treating microbial infections, particularly managing drug-resistant strains of bacteria. One of the major challenges in designing phage-based therapy is to identify the most appropriate potential phage candidate to treat bacterial infections. In this study, an attempt has been made to predict phage-host interactions with high accuracy to identify the potential bacteriophage that can be used for treating a bacterial infection. The developed models have been created using a training dataset containing 826 phage- host interactions, and have been evaluated on a validation dataset comprising 1,201 phage-host interactions. Firstly, alignment-based models have been developed using similarity between phage-phage (BLASTPhage), host-host (BLASTHost) and phage-CRISPR (CRISPRPred), where we achieved accuracy between 42.4-66.2% for BLASTPhage, 55-78.4% for BLASTHost, and 43.7-80.2% for CRISPRPred across five taxonomic levels. Secondly, alignment free models have been developed using machine learning techniques. Thirdly, hybrid models have been developed by integrating the alignment-free models and the similarity-scores where we achieved maximum performance of (60.6-93.5%). Finally, an ensemble model has been developed that combines the hybrid and alignment-based models. Our ensemble model achieved highest accuracy of 67.9, 80.6, 85.5, 90, and 93.5% at Genus, Family, Order, Class, and Phylum levels on validation dataset. In order to serve the scientific community, we have also developed a webserver named PhageTB and provided a standalone software package (https://webs.iiitd.edu.in/raghava/phagetb/) for the same.

"French Phage Network" Annual Conference-Seventh Meeting Report.

The French Phage Network (Phages.fr) has continuously grown since its foundation, eight years ago. The annual conference, held at the Institut Pasteur in Paris, attracted 164 participants from the 11th to the 13th of October 2022. Researchers from academic laboratories, hospitals and private companies shared their ongoing projects and breakthroughs in the very institute where Felix d'Hérelle developed phage therapy over a century ago. The conference was divided into four thematic sessions, each opened by a keynote lecture: "Interaction between phages, mobile genetic elements and bacterial immune system," "Ecology and evolution of phage-bacteria interactions," "Molecular interplay between phages and their hosts" and "Therapeutic and biotechnological applications of phages." A total of 32 talks and 33 posters were presented during the conference.

Engineering bacteriophages for enhanced host range and efficacy: insights from bacteriophage-bacteria interactions.

Bacteriophages, the most abundant organisms on earth, have the potential to address the rise of multidrug-resistant bacteria resulting from the overuse of antibiotics. However, their high specificity and limited host range can hinder their effectiveness. Phage engineering, through the use of gene editing techniques, offers a means to enhance the host range of bacteria, improve phage efficacy, and facilitate efficient cell-free production of phage drugs. To engineer phages effectively, it is necessary to understand the interaction between phages and host bacteria. Understanding the interaction between the receptor recognition protein of bacteriophages and host receptors can serve as a valuable guide for modifying or replacing these proteins, thereby altering the receptor range of the bacteriophage. Research and development focused on the CRISPR-Cas bacterial immune system against bacteriophage nucleic acids can provide the necessary tools to promote recombination and counter-selection in engineered bacteriophage programs. Additionally, studying the transcription and assembly functions of bacteriophages in host bacteria can facilitate the engineered assembly of bacteriophage genomes in non-host environments. This review highlights a comprehensive summary of phage engineering methods, including in-host and out-of-host engineering, and the use of high-throughput methods to understand their role. The main aim of these techniques is to harness the intricate interactions between bacteriophages and hosts to inform and guide the engineering of bacteriophages, particularly in the context of studying and manipulating the host range of bacteriophages. By employing advanced high-throughput methods to identify specific bacteriophage receptor recognition genes, and subsequently introducing modifications or performing gene swapping through in-host recombination or out-of-host synthesis, it becomes possible to strategically alter the host range of bacteriophages. This capability holds immense significance for leveraging bacteriophages as a promising therapeutic approach against antibiotic-resistant bacteria.

Host-phage interactions and modeling for therapy.

Phage are drivers of numerous ecological processes on the planet and have the potential to be developed into a therapy alternative to antibiotics. Phage at all points of their life cycle, from initiation of infection to their release, interact with their host in some manner. More importantly, to harness their antimicrobial potential it is vital to understand how phage interact with the eukaryotic environment in the context of applying phage for therapy. In this chapter, the various mechanisms of phage interplay with their hosts as part of their natural life cycle are discussed in depth for Gram-positive and negative bacteria. Further, the literature surrounding the various models utilized to develop phage as a therapeutic are examined, and how these models may improve our understanding of phage-host interactions and current progress in utilizing phage for therapy in the clinical environment.

Proteomic Study of the Interactions between Phages and the Bacterial Host Klebsiella pneumoniae.

Phages and bacteria have acquired resistance mechanisms for protection. In this context, the aims of the present study were to analyze the proteins isolated from 21 novel lytic phages of Klebsiella pneumoniae in search of defense mechanisms against bacteria and also to determine the infective capacity of the phages. A proteomic study was also conducted to investigate the defense mechanisms of two clinical isolates of K. pneumoniae infected by phages. For this purpose, the 21 lytic phages were sequenced and de novo assembled. The host range was determined in a collection of 47 clinical isolates of K. pneumoniae, revealing the variable infective capacity of the phages. Genome sequencing showed that all of the phages were lytic phages belonging to the order Caudovirales. Phage sequence analysis revealed that the proteins were organized in functional modules within the genome. Although most of the proteins have unknown functions, multiple proteins were associated with defense mechanisms against bacteria, including the restriction-modification system, the toxin-antitoxin system, evasion of DNA degradation, blocking of host restriction and modification, the orphan CRISPR-Cas system, and the anti-CRISPR system. Proteomic study of the phage-host interactions (i.e., between isolates K3574 and K3320, which have intact CRISPR-Cas systems, and phages vB_KpnS-VAC35 and vB_KpnM-VAC36, respectively) revealed the presence of several defense mechanisms against phage infection (prophage, defense/virulence/resistance, oxidative stress and plasmid proteins) in the bacteria, and of the Acr candidate (anti-CRISPR protein) in the phages. IMPORTANCE Researchers, including microbiologists and infectious disease specialists, require more knowledge about the interactions between phages and their bacterial hosts and about their defense mechanisms. In this study, we analyzed the molecular mechanisms of viral and bacterial defense in phages infecting clinical isolates of K. pneumoniae. Viral defense mechanisms included restriction-modification system evasion, the toxin-antitoxin (TA) system, DNA degradation evasion, blocking of host restriction and modification, and resistance to the abortive infection system, anti-CRISPR and CRISPR-Cas systems. Regarding bacterial defense mechanisms, proteomic analysis revealed expression of proteins involved in the prophage (FtsH protease modulator), plasmid (cupin phosphomannose isomerase protein), defense/virulence/resistance (porins, efflux pumps, lipopolysaccharide, pilus elements, quorum network proteins, TA systems, and methyltransferases), oxidative stress mechanisms, and Acr candidates (anti-CRISPR protein). The findings reveal some important molecular mechanisms involved in the phage-host bacterial interactions; however, further study in this field is required to improve the efficacy of phage therapy.

Host interactions of novel Crassvirales species belonging to multiple families infecting bacterial host, Bacteroides cellulosilyticus WH2.

Bacteroides, the prominent bacteria in the human gut, play a crucial role in degrading complex polysaccharides. Their abundance is influenced by phages belonging to the Crassvirales order. Despite identifying over 600 Crassvirales genomes computationally, only few have been successfully isolated. Continued efforts in isolation of more Crassvirales genomes can provide insights into phage-host-evolution and infection mechanisms. We focused on wastewater samples, as potential sources of phages infecting various Bacteroides hosts. Sequencing, assembly, and characterization of isolated phages revealed 14 complete genomes belonging to three novel Crassvirales species infecting Bacteroides cellulosilyticus WH2. These species, Kehishuvirus sp. 'tikkala' strain Bc01, Kolpuevirus sp. 'frurule' strain Bc03, and 'Rudgehvirus jaberico' strain Bc11, spanned two families, and three genera, displaying a broad range of virion productions. Upon testing all successfully cultured Crassvirales species and their respective bacterial hosts, we discovered that they do not exhibit co-evolutionary patterns with their bacterial hosts. Furthermore, we observed variations in gene similarity, with greater shared similarity observed within genera. However, despite belonging to different genera, the three novel species shared a unique structural gene that encodes the tail spike protein. When investigating the relationship between this gene and host interaction, we discovered evidence of purifying selection, indicating its functional importance. Moreover, our analysis demonstrated that this tail spike protein binds to the TonB-dependent receptors present on the bacterial host surface. Combining these observations, our findings provide insights into phage-host interactions and present three Crassvirales species as an ideal system for controlled infectivity experiments on one of the most dominant members of the human enteric virome.

PHISDetector: A Tool to Detect Diverse In Silico Phage-host Interaction Signals for Virome Studies.

Phage-microbe interactions are appealing systems to study coevolution, and have also been increasingly emphasized due to their roles in human health, disease, and the development of novel therapeutics. Phage-microbe interactions leave diverse signals in bacterial and phage genomic sequences, defined as phage-host interaction signals (PHISs), which include clustered regularly interspaced short palindromic repeats (CRISPR) targeting, prophage, and protein-protein interaction signals. In the present study, we developed a novel tool phage-host interaction signal detector (PHISDetector) to predict phage-host interactions by detecting and integrating diverse in silico PHISs, and scoring the probability of phage-host interactions using machine learning models based on PHIS features. We evaluated the performance of PHISDetector on multiple benchmark datasets and application cases. When tested on a dataset of 758 annotated phage-host pairs, PHISDetector yields the prediction accuracies of 0.51 and 0.73 at the species and genus levels, respectively, outperforming other phage-host prediction tools. When applied to on 125,842 metagenomic viral contigs (mVCs) derived from 3042 geographically diverse samples, a detection rate of 54.54% could be achieved. Furthermore, PHISDetector could predict infecting phages for 85.6% of 368 multidrug-resistant (MDR) bacteria and 30% of 454 human gut bacteria obtained from the National Institutes of Health (NIH) Human Microbiome Project (HMP). The PHISDetector can be run either as a web server (http://www.microbiome-bigdata.com/PHISDetector/) for general users to study individual inputs or as a stand-alone version (https://github.com/HIT-ImmunologyLab/PHISDetector) to process massive phage contigs from virome studies.