Phage-host Interactions Research Articles

Genomic diversity of phages infecting the globally widespread genus Sulfurimonas

The widespread bacterial genus Sulfurimonas is metabolically versatile and occupies a key ecological niche in different habitats, but its interaction with bacteriophages remains unexplored. Here we systematically investigated the genetic diversity, taxonomy and interaction patterns of Sulfurimonas-associated phages based on sequenced microbial genomes and metagenomes. High-confidence phage contigs related to Sulfurimonas were retrieved from various ecosystems, clustered into 61 viral operational taxonomic units across three viral realms, including Duplodnaviria, Monodnaviria and Varidnaviria. Head-tail phages of Caudoviricetes were assigned to 19 genus-level viral clusters, the majority of which were distantly related to known viruses. Notably, diverse double jelly-roll viruses and inoviruses were also linked to Sulfurimonas, representing two commonly overlooked phage groups. Historical and current phage infections were revealed, implying viral impact on the evolution of host adaptive immunity. Additionally, phages carrying auxiliary metabolic genes might benefit hosts by compensating or augmenting sulfur metabolism. This study highlights the diversity and novelty of Sulfurimonas-associated phages with divergent tailless lineages, providing basis for further investigation of phage-host interactions within this genus.

Prediction of strain level phage-host interactions across the Escherichia genus using only genomic information.

Predicting bacteriophage infection of specific bacterial strains promises advancements in phage therapy and microbial ecology. Whether the dynamics of well-established phage-host model systems generalize to the wide diversity of microbes is currently unknown. Here we show that we could accurately predict the outcomes of phage-bacteria interactions at the strain level in natural isolates from the genus Escherichia using only genomic data (area under the receiver operating characteristic curve (AUROC) of 86%). We experimentally established a dataset of interactions between 403 diverse Escherichia strains and 96 phages. Most interactions are explained by adsorption factors as opposed to antiphage systems which play a marginal role. We trained predictive algorithms and pinpoint poorly predicted interactions to direct future research efforts. Finally, we established a pipeline to recommend tailored phage cocktails, demonstrating efficiency on 100 pathogenic E. coli isolates. This work provides quantitative insights into phage-host specificity and supports the use of predictive algorithms in phage therapy.

Identification of receptor-binding protein and host receptor of non-lytic dsRNA phage phiNY.

To date, complete genome sequences of 14 double-stranded RNA (dsRNA) phages are available, and studies have shown that the host range of dsRNA phages is limited. The hosts of most dsRNA phages belong to the genus Pseudomonas. However, the dsRNA phage phiNY, which has a non-lytic life cycle, was isolated from Microvirgula aerodenitrificans. Currently, the interaction between dsRNA phage phiNY and its host bacteria is unclear, which is not beneficial to a comprehensive understanding of dsRNA phage biology and the exploitation of dsRNA phage with non-lytic life cycle for biomedical applications and others. Phage adsorption is a crucial step through the interactions between receptor-binding protein (RBP) of the phage and its receptors to initiate the infection process, which dictates host range specificity. Thus, we identified the RBP and host receptor of phiNY. Through homology alignment, amino acid sequence similarity analysis, and the phylogenetic tree analysis, orf11, located in the M-segment of dsRNA phage phiNY, encodes a putative RBP. We further performed the whole-cell enzyme-linked immunosorbent assay (ELISA), western blotting assay, and indirect immunofluorescence assay and demonstrated that this orf11 is an RBP. Finally, using affinity chromatography, ELISA, and dynamic light scattering, we identified lipopolysaccharides (LPSs) on the surface of the host M. aerodenitrificans strain LH9 as host receptors involved in the adsorption of the dsRNA bacteriophage phiNY and observed the state of phiNY RBP after combining with LPS by atomic force microscopy. These results will guide future studies on phage-host interaction in a dsRNA phage with a non-lytic life cycle.IMPORTANCEThe interactions between the lytic dsRNA phages and their host receptors have been clarified in previous studies. However, the interaction between the dsRNA phage phiNY (which has a non-lytic life cycle) and its host receptors during the dsRNA phage adsorption process was unknown. Here, we found that phiNY uses the orf11 protein as a receptor-binding protein (RBP). In addition, we found that this orf11 recognizes lipopolysaccharide from the host bacterium Microvirgula aerodenitrificans strain LH9 as a specific receptor. These results suggest that phiNY, like lytic dsRNA phages, uses an RBP to bind to a similar host receptor (i.e., lipopolysaccharide). Determining the interaction between the dsRNA phage phiNY and its host receptors will help to elucidate the mechanisms underlying the phiNY non-lytic life cycle and enhance our understanding of its infection mechanism.

Deciphering the Genetic Architecture of Staphylococcus warneri Prophage vB_G30_01: A Comprehensive Molecular Analysis.

The current knowledge of Staphylococcus warneri phages is limited, with few genomes sequenced and characterized. In this study, a prophage, vB_G30_01, isolated from Staphylococcus warneri G30 was characterized and evaluated for its lysogenic host range. The phage was studied using transmission electron microscopy and a host range. The phage genome was sequenced and characterized in depth, including phylogenetic and taxonomic analyses. The linear dsDNA genome of vB_G30_01 contains 67 predicted open reading frames (ORFs), classifying it within Bronfenbrennervirinae. With a total of 10 ORFs involved in DNA replication-related and transcriptional regulator functions, vB_G30_01 may play a role in the genetics and transcription of a host. Additionally, vB_G30_01 possesses a complete set of genes related to host lysogeny and lysis, implying that vB_G30_01 may influence the survival and adaptation of its host. Furthermore, a comparative genomic analysis reveals that vB_G30_01 shares high genomic similarity with other Staphylococcus phages and is relatively closely related to those of Exiguobacterium and Bacillus, which, in combination with the cross-infection assay, suggests possible cross-species infection capabilities. This study enhances the understanding of Staphylococcus warneri prophages, providing insights into phage-host interactions and potential horizontal gene transfer.

Mechanistic modeling of the dynamics of phage attack during milk acidification in the cheese-making process

Bacteriophage attacks represent a major threat in the dairy industry. Here, an unstructured mechanistic model predicting the dynamics of milk acidification in case of phage attack was developed and experimentally validated. Multiple acidification experiments were run with different combinations of initial phage titers and bacterial concentrations and the resulting pH dynamics were recorded. The model could successfully predict the success or failure of milk acidification. Using the model, important biological parameters were deduced from simple, low-cost acidification measurements. These parameters included bacteria’s maximum growth and lysis rates, phages’ burst size, etc. Sensitivity analysis helped identify biologically relevant aspects of phage-host interactions. Growth and lysis kinetics were shown to have the most important impacts. This knowledge can be used to develop easy routine strategies to fight phage attack in the dairy industry. The model can be used to raise awareness amongst cheese makers on the importance of cleaning to avoid food and material waste.

Isolation and characterization of the new Streptomyces phages Kamino, Geonosis, Abafar, and Scarif infecting a broad range of host species.

Streptomyces, a multifaceted genus of soil-dwelling bacteria belonging to the phylum Actinomycetota, features intricate phage-host interactions shaped by its complex life cycle and the synthesis of a diverse array of specialized metabolites. Here, we describe the isolation and characterization of four novel Streptomyces phages infecting a variety of different host species. While phage Kamino, isolated on Streptomyces kasugaensis, is predicted to be temperate and encodes a serine integrase in its genome, phages Geonosis (isolated on Streptomyces griseus) and Abafar and Scarif, isolated on Streptomyces albidoflavus, are virulent phages. Phages Kamino and Geonosis were shown to amplify well in liquid culture leading to a pronounced culture collapse already at low titers. Determination of the host range by testing >40 different Streptomyces species identified phages Kamino, Abafar, and Scarif as broad host-range phages. Overall, the phages described in this study expand the publicly available portfolio of phages infecting Streptomyces and will be instrumental in advancing the mechanistic understanding of the intricate antiviral strategies employed by these multicellular bacteria.IMPORTANCEThe actinobacterial genus Streptomyces is characterized by multicellular, filamentous growth and the synthesis of a diverse range of bioactive molecules. These characteristics also play a role in shaping their interactions with the most abundant predator in the environment, bacteriophages-viruses infecting bacteria. In this study, we characterize four new phages infecting Streptomyces. Out of those, three phages feature a broad host range infecting up to 15 different species. The isolated phages were characterized with respect to plaque and virion morphology, host range, and amplification in liquid culture. In summary, the phages reported in this study contribute to the broader collection of publicly available phages infecting Streptomyces, playing a crucial role in advancing our mechanistic understanding of phage-host interactions of these multicellular bacteria.

Investigating potential auxiliary anaerobic digestion activity of phage under polyvinyl chloride microplastic stress

Polyvinyl chloride (PVC) microplastics present in sewage were trapped in sludge, thereby hindering anaerobic digestion performance of waste active sludge (WAS). Phages regulate virocell metabolism by encoding auxiliary metabolic genes (AMGs) related to energy acquisition and material degradation, supporting hosts survive in harsh environments and play a crucial role in biogeochemical cycles. This study investigated the potential effects of phages on the recovery of WAS anaerobic digestion under PVC stress. We observed a significant alteration in the phage community induced by PVC microplastics. Phages encoded AMGs related to anaerobic digestion and cell growth probably alleviate PVC microplastics inhibition on WAS anaerobic digestion, and 54.2 % of hydrolysis-related GHs and 40.8 % of acidification-related AMGs were actively transcribed in the PVC-exposed group. Additionally, the degradation of chitin and peptidoglycan during hydrolysis and the conversion of glucose to pyruvate during acidification were more susceptible to phages. Prediction of phage-host relationship indicated that the phyla Pseudomonadota were predominantly targeted hosts by hydrolysis-related and acidification-related phages, and PVC toxicity had minimal impact on phage-host interaction. Our findings highlight the importance of phages in anaerobic digestion and provide a novel strategy for using phages in the functional recovery of microplastic-exposed sludge.

Exploring Burkholderia pseudomallei-specific bacteriophages: overcoming O-antigen specificity and adaptive mutation in phage tail fiber

IntroductionBurkholderia pseudomallei, a Gram-negative bacterium inhabiting soil and fresh water, is the causative agent of melioidosis, a formidable disease in the tropics. The emergence of antibiotic resistance and the extended duration of treatment, up to 20 weeks, have posed significant challenges in combatting melioidosis. As an alternative approach, bacteriophage therapy is being explored.MethodsTo identify the most promising bacteriophage for future therapeutic applications, we designed a screening process to address the barrier of phage specificity due to the O-antigen receptor diversity. By using two biosafe strains, Bp82 (O-antigen type A) and 576mn (O-antigen type B), to represent the major serotype A and B, we screened 145 phage samples collected from soil and water in southern Thailand.ResultsTen of them demonstrated the ability to overcome differences in O-antigen types, yielding positive plaques formed on culture of both bacterial strains. Subsequently, we isolated 22 bacteriophages from these samples, one was adaptively mutated during the screening process, named ΦPK23V1, which had the ability to infect up to 83.3% (115/138) of tested B. pseudomallei strains, spanning both serogroups. Employing a panel of surface polysaccharide antigen mutant strains, we explored the role of capsular polysaccharide (CPS) and O-antigens as essential components for phage infection. All isolated phages were classified into the P2-like myophage group. Additionally, our research revealed a point mutation in the phage tail fiber gene (gpH), expanding the host range of ΦPK23V1, even in the absence of CPS and O-antigens.DiscussionHowever, it was evident that ΦPK23V1 is a lysogenic phage, which cannot be readily applied for therapeutic use. This discovery sheds light on the receptor binding domain of P2-like bacteriophages in B. pseudomallei. Collectively, our study has identified bacteriophages with a broad host range within B. pseudomallei strains, enhancing our understanding of phage–host interactions and offering insights into the role of the phage tail fiber gene in host cell entry.

Quorum sensing positively regulates CPS-dependent Autographiviridae phage infection in Vibrio alginolyticus.

Phage resistance is a direct threat to phage therapy, and understanding phage-host interactions, especially how bacteria block phage infection, is essential for developing successful phage therapy. In the present study, we demonstrate for the first time that Vibrio alginolyticus uses quorum sensing (QS) to promote capsular polysaccharide (CPS)-specific phage infection by upregulating ugd expression, which is necessary for the synthesis of Autographiviridae phage receptor CPS. Although increased CPS-specific phage susceptibility is a novel trade-off mediated by QS, it results in the upregulation of virulence factors, promoting biofilm development and enhanced capsular polysaccharide production in V. alginolyticus. This suggests that inhibiting QS may improve the effectiveness of antibiotic treatment, but it may also reduce the efficacy of phage therapy.

Characterization and comparative genomic analysis of a marine Bacillus phage reveal a novel viral genus.

Although recent metagenomics research has obtained a wealth of phage genetic information, much of it is considered "dark matter" because of the lack of similarity with known sequences in the database. Therefore, the isolation and characterization of novel phages will help to interpret the vast unknown viral metagenome data and improve our understanding of phage diversity and phage-host interactions. Bacillus pumilus shows high economic relevance due to its wide applications in biotechnology, industry, biopharma, and environmental sectors. Since phages influence the abundance, metabolism, evolution, fitness, and ecological functions of bacteria through complex interactions, the significance of isolation and characterization of novel phages infecting B. pumilus is apparent. In this study, we isolated and characterized a B. pumilus phage belonging to a novel viral genus, which provides essential knowledge for phage biology as well as the industrial application of B. pumilus.

DepoScope: Accurate phage depolymerase annotation and domain delineation using large language models.

Bacteriophages (phages) are viruses that infect bacteria. Many of them produce specific enzymes called depolymerases to break down external polysaccharide structures. Accurate annotation and domain identification of these depolymerases are challenging due to their inherent sequence diversity. Hence, we present DepoScope, a machine learning tool that combines a fine-tuned ESM-2 model with a convolutional neural network to identify depolymerase sequences and their enzymatic domains precisely. To accomplish this, we curated a dataset from the INPHARED phage genome database, created a polysaccharide-degrading domain database, and applied sequential filters to construct a high-quality dataset, which is subsequently used to train DepoScope. Our work is the first approach that combines sequence-level predictions with amino-acid-level predictions for accurate depolymerase detection and functional domain identification. In that way, we believe that DepoScope can greatly enhance our understanding of phage-host interactions at the level of depolymerases.

Serratia marcescens ATCC 274 increases production of the red pigment prodigiosin in response to Chi phage infection

Serratia marcescens is an opportunistic human pathogen that produces a vibrant red pigment called prodigiosin. Prodigiosin has implications in virulence of S. marcescens and promising clinical applications. We discovered that addition of the virulent flagellotropic bacteriophage χ (Chi) to a culture of S. marcescens stimulates a greater than fivefold overproduction of prodigiosin. Active phage infection is required for the effect, as a χ-resistant strain lacking flagella does not respond to phage presence. Via a reporter fusion assay, we have determined that the addition of a χ-induced S. marcescens cell lysate to an uninfected culture causes a threefold increase in transcription of the pig operon, containing genes essential for pigment biosynthesis. Replacement of the pig promoter with a constitutive promoter abolished the pigmentation increase, indicating that regulatory elements present in the pig promoter likely mediate the phenomenon. We hypothesize that S. marcescens detects the threat of phage-mediated cell death and reacts by producing prodigiosin as a stress response. Our findings are of clinical significance for two main reasons: (i) elucidating complex phage-host interactions is crucial for development of therapeutic phage treatments, and (ii) overproduction of prodigiosin in response to phage could be exploited for its biosynthesis and use as a pharmaceutical.

DeePhafier: a phage lifestyle classifier using a multilayer self-attention neural network combining protein information.

Bacteriophages are the viruses that infect bacterial cells. They are the most diverse biological entities on earth and play important roles in microbiome. According to the phage lifestyle, phages can be divided into the virulent phages and the temperate phages. Classifying virulent and temperate phages is crucial for further understanding of the phage-host interactions. Although there are several methods designed for phage lifestyle classification, they merely either consider sequence features or gene features, leading to low accuracy. A new computational method, DeePhafier, is proposed to improve classification performance on phage lifestyle. Built by several multilayer self-attention neural networks, a global self-attention neural network, and being combined by protein features of the Position Specific Scoring Matrix matrix, DeePhafier improves the classification accuracy and outperforms two benchmark methods. The accuracy of DeePhafier on five-fold cross-validation is as high as 87.54% for sequences with length >2000bp.

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.

Extremely halophilic brine community manipulation shows higher robustness of microbiomes inhabiting human-driven solar saltern than naturally driven lake.

Viruses greatly influence succession and diversification of their hosts, yet the effects of viral infection on the ecological dynamics of natural microbial populations remain poorly understood, especially at finer scales of diversity. By manipulating the viral predation pressure by autochthonous and allochthonous viruses, we uncovered potential phage-host interaction, and their important role in structuring the prokaryote community at an ecotype level.

Expanding the Phage Galaxy: Isolation and Characterization of Five Novel Streptomyces Siphoviruses Ankus, Byblos, DekoNeimoidia, Mandalore, and Naboo.

Key features of the actinobacterial genus Streptomyces are multicellular, filamentous growth, and production of a broad portfolio of bioactive molecules. These characteristics appear to play an important role in phage-host interactions and are modulated by phages during infection. To accelerate research of such interactions and the investigation of novel immune systems in multicellular bacteria, phage isolation, sequencing, and characterization are needed. This is a prerequisite for establishing systematic collections that appropriately cover phage diversity for comparative analyses. As part of a public outreach program within the priority program SPP 2330, involving local schools, we describe the isolation and characterization of five novel Streptomyces siphoviruses infecting S. griseus, S. venezuelae, and S. olivaceus. All isolates are virulent members of two existing genera and, additionally, establish a new genus in the Stanwilliamsviridae family. In addition to an extensive set of tRNAs and proteins involved in phage replication, about 80% of phage genes encode hypothetical proteins, underlining the yet underexplored phage diversity and genomic dark matter still found in bacteriophages infecting actinobacteria. Taken together, phages Ankus, Byblos, DekoNeimoidia, Mandalore, and Naboo expand the phage diversity and contribute to ongoing research in the field of Streptomyces phage-host interactions.

Determinants in the phage life cycle: The dynamic nature of ssDNA phage FLiP and host interactions under varying environmental conditions and growth phases.

The influence of environmental factors on the interactions between phages and bacteria, particularly single-stranded DNA (ssDNA) phages, has been largely unexplored. In this study, we used Finnlakevirus FLiP, the first known ssDNA phage species with a lipid membrane, as our model phage. We examined the infectivity of FLiP with three Flavobacterium host strains, B330, B167 and B114. We discovered that FLiP infection is contingent on the host strain and conditions such as temperature and bacterial growth phase. FLiP can infect its hosts across a wide temperature range, but optimal phage replication varies with each host. We uncovered some unique aspects of phage infectivity: FLiP has limited infectivity in liquid-suspended cells, but it improves when cells are surface-attached. Moreover, FLiP infects stationary phase B167 and B114 cells more rapidly and efficiently than exponentially growing cells, a pattern not observed with the B330 host. We also present the first experimental evidence of endolysin function in ssDNA phages. The activity of FLiP's lytic enzymes was found to be condition-dependent. Our findings underscore the importance of studying phage ecology in contexts that are relevant to the environment, as both the host and the surrounding conditions can significantly alter the outcome of phage-host interactions.

Screening of the PA14NR Transposon Mutant Library Identifies Genes Involved in Resistance to Bacteriophage Infection in Pseudomomas aeruginosa.

Multidrug-resistant P. aeruginosa infections pose a serious public health threat due to the rise in antimicrobial resistance. Phage therapy has emerged as a promising alternative. However, P. aeruginosa has evolved various mechanisms to thwart phage attacks, making it crucial to decipher these resistance mechanisms to develop effective therapeutic strategies. In this study, we conducted a forward-genetic screen of the P. aeruginosa PA14 non-redundant transposon library (PA14NR) to identify dominant-negative mutants displaying phage-resistant phenotypes. Our screening process revealed 78 mutants capable of thriving in the presence of phages, with 23 of them carrying insertions in genes associated with membrane composition. Six mutants exhibited total resistance to phage infection. Transposon insertions were found in genes known to be linked to phage-resistance such as galU and a glycosyl transferase gene, as well as novel genes such as mexB, lasB, and two hypothetical proteins. Functional experiments demonstrated that these genes played pivotal roles in phage adsorption and biofilm formation, indicating that altering the bacterial membrane composition commonly leads to phage resistance in P. aeruginosa. Importantly, these mutants displayed phenotypic trade-offs, as their resistance to phages inversely affected antibiotic resistance and hindered biofilm formation, shedding light on the complex interplay between phage susceptibility and bacterial fitness. This study highlights the potential of transposon mutant libraries and forward-genetic screens in identifying key genes involved in phage-host interactions and resistance mechanisms. These findings support the development of innovative strategies for combating antibiotic-resistant pathogens.

Relevance of the bacteriophage adherence to mucus model for Pseudomonas aeruginosa phages.

Pseudomonas aeruginosa infections are getting increasingly serious as antimicrobial resistance spreads. Phage therapy may be a solution to the problem, especially if improved by current advances on phage-host studies. As a mucosal pathogen, we hypothesize that P. aeruginosa and its phages are linked to the bacteriophage adherence to mucus (BAM) model. This means that phage-host interactions could be influenced by mucin presence, impacting the success of phage infections on the P. aeruginosa host and consequently leading to the protection of the metazoan host. By using a group of four different phages, we tested three important phenotypes associated with the BAM model: phage binding to mucin, phage growth in mucin-exposed hosts, and the influence of mucin on CRISPR immunity of the bacterium. Three of the tested phages significantly bound to mucin, while two had improved growth rates in mucin-exposed hosts. Improved phage growth was likely the result of phage exploitation of mucin-induced physiological changes in the host. We could not detect CRISPR activity in our system but identified two putative anti-CRISPR proteins coded by the phage. Overall, the differential responses seen for the phages tested show that the same bacterial species can be targeted by mucosal-associated phages or by phages not affected by mucus presence. In conclusion, the BAM model is relevant for phage-bacterium interactions in P. aeruginosa, opening new possibilities to improve phage therapy against this important pathogen by considering mucosal interaction dynamics.IMPORTANCESome bacteriophages are involved in a symbiotic relationship with animals, in which phages held in mucosal surfaces protect them from invading bacteria. Pseudomonas aeruginosa is one of the many bacterial pathogens threatening humankind during the current antimicrobial resistance crisis. Here, we have tested whether P. aeruginosa and its phages are affected by mucosal conditions. We discovered by using a collection of four phages that, indeed, mucosal interaction dynamics can be seen in this model. Three of the tested phages significantly bound to mucin, while two had improved growth rates in mucin-exposed hosts. These results link P. aeruginosa and its phages to the bacteriophage adherence to the mucus model and open opportunities to explore this to improve phage therapy, be it by exploiting the phenotypes detected or by actively selecting mucosal-adapted phages for treatment.

The trade-off of Vibrio parahaemolyticus between bacteriophage resistance and growth competitiveness.

Vibrio parahaemolyticus is a food-borne pathogen, which is often isolated from various seafood products. In this study, two kinds of bacteriophages was isolated from the offshore sediments samples. The anti-phage mutant strain were obtained after seventeen rounds of co-culture of Vibrio parahaemolyticus and mixed bacteriophage, multigroup sequencing was carried out on spontaneous the anti-phage mutant strain and the wild-type strain. We used the Sanger sequencing to verify the accuracy of the mutation sites. Biolog GEN III MicroPlates were used to evaluate the metabolic capacity of wild-type strains and the anti-phage mutant strain. In this study, we found that with flaG gene (slight homology to N terminus of multiple flagellins) mutated, making the bacteriophage unable to absorb to the cell surface of the host. And, the growth competitiveness of the anti-phage mutant strain is lower than the wild-type strain. These results indicated that the fitness cost, including loss of the growth competitiveness, constitutes a barrier to the prevalence of these defense mechanisms. And the selection pressure on different anti-phage strategies depends on the trade-off between mortality imposed by bacteriophages and fitness cost of the defense strategy under the given environmental conditions. In conclusion, this study provides valuable insights into the phage-host interaction and phage resistance in Vibrio parahaemolyticus. Our study provided knowledge for the evolutionary adaption of bacteria against the bacteriophage, which could add more information to understand the phage resistance mechanism before applying in the industry.