Prophage Tracer: precisely tracing prophages in prokaryotic genomes using overlapping split-read alignment

V- and VL-Scores Uncover Viral Signatures and Origins of Protein Families

Viruses are key drivers of microbial diversity, nutrient cycling, and co-evolution in ecosystems, yet their study is hindered due to challenges in culturing. Traditional gene-centric methods, which focus on a few hallmark genes like for capsids, miss much of the viral genome, leaving key viral proteins and functions undiscovered. Here, we introduce two powerful annotation-free metrics, V-score and VL-score, designed to quantify the "virus-likeness" of protein families and genomes and create an open-access searchable database, 'V-Score-Search'. By applying V- and VL-scores to public databases (KEGG, Pfam, and eggNOG), we link 38-77% of protein families with viruses, a 9-16x increase over current estimates. These metrics outperform existing approaches, enabling precise detection of viral genomes, prophages, and host-derived auxiliary viral genes (AVGs) from fragmented sequences, and significantly improving genome binning. Remarkably, we identify up to 17x more AVGs, dominated by non-metabolic proteins of unknown function. This innovation unlocks new insights into virus signatures and host interactions, with wide-ranging implications from genomics to biotechnology.

Genomic characterization of SNW-1, a novel prophage of the deep-sea vent chemolithoautotroph Sulfurimonas indica NW79

The globally widespread genus Sulfurimonas are playing important roles in different habitats, including the deep-sea hydrothermal vents. However, phages infecting Sulfurimonas have never been isolated and characterized to date. In the present study, a novel prophage SNW-1 was identified from Sulfurimonas indica NW79. Whole genome sequencing resulted in a circular, double-stranded DNA molecule of 37,096 bp with a mol% G+C content of 37. The genome includes 64 putative open reading frames, 33 of which code for proteins with predicted functions. Presence of hallmark genes associated with Caudoviricetes and genes involved in lysis and lysogeny indicated that SNW-1 should be a temperate, tailed phage. Phylogenetic and comparative proteomic analyses suggested that Sulfurimonas phage SNW-1 was distinct from other double stranded DNA phages and might represent a new viral genus.

A nuclease domain fused to the Snf2 helicase confers antiphage defence in coral-associated Halomonas meridiana

The coral reef microbiome plays a vital role in the health and resilience of reefs. Previous studies have examined phage therapy for coral pathogens and for modifying the coral reef microbiome, but defence systems against coral-associated bacteria have received limited attention. Phage defence systems play a crucial role in helping bacteria fight phage infections. In this study, we characterized a new defence system, Hma (HmaA-HmaB-HmaC), in the coral-associated Halomonas meridiana derived from the scleractinian coral Galaxea fascicularis. The Swi2/Snf2 helicase HmaA with a C-terminal nuclease domain exhibits antiviral activity against Escherichia phage T4. Mutation analysis revealed the nickase activity of the nuclease domain (belonging to PDD/EXK superfamily) of HmaA is essential in phage defence. Additionally, HmaA homologues are present in ~1000 bacterial and archaeal genomes. The high frequency of HmaA helicase in Halomonas strains indicates the widespread presence of these phage defence systems, while the insertion of defence genes in the hma region confirms the existence of a defence gene insertion hotspot. These findings offer insights into the diversity of phage defence systems in coral-associated bacteria and these diverse defence systems can be further applied into designing probiotics with high-phage resistance.

VIBES: a workflow for annotating and visualizing viral sequences integrated into bacterial genomes

Bacteriophages are viruses that infect bacteria. Many bacteriophages integrate their genomes into the bacterial chromosome and become prophages. Prophages may substantially burden or benefit host bacteria fitness, acting in some cases as parasites and in others as mutualists. Some prophages have been demonstrated to increase host virulence. The increasing ease of bacterial genome sequencing provides an opportunity to deeply explore prophage prevalence and insertion sites. Here we present VIBES (Viral Integrations in Bacterial genomES), a workflow intended to automate prophage annotation in complete bacterial genome sequences. VIBES provides additional context to prophage annotations by annotating bacterial genes and viral proteins in user-provided bacterial and viral genomes. The VIBES pipeline is implemented as a Nextflow-driven workflow, providing a simple, unified interface for execution on local, cluster and cloud computing environments. For each step of the pipeline, a container including all necessary software dependencies is provided. VIBES produces results in simple tab-separated format and generates intuitive and interactive visualizations for data exploration. Despite VIBES's primary emphasis on prophage annotation, its generic alignment-based design allows it to be deployed as a general-purpose sequence similarity search manager. We demonstrate the utility of the VIBES prophage annotation workflow by searching for 178 Pf phage genomes across 1072 Pseudomonas spp. genomes.

Combining CRISPR/Cas9 and natural excision for the precise and complete removal of mobile genetic elements in bacteria

Horizontal gene transfer, facilitated by mobile genetic elements (MGEs), is an adaptive evolutionary process that contributes to the evolution of bacterial populations and infectious diseases. A variety of MGEs not only can integrate into the bacterial genome but also can survive or even replicate like plasmids in the cytoplasm, thus requiring precise and complete removal for studying their strategies in benefiting host cells. Existing methods for MGE removal, such as homologous recombination-based deletion and excisionase-based methods, have limitations in effectively eliminating certain MGEs. To overcome these limitations, we developed the Cas9-NE method, which combines the CRISPR/Cas9 system with the natural excision of MGEs. In this approach, a specialized single guide RNA (sgRNA) element is designed with a 20-nucleotide region that pairs with the MGE sequence. This sgRNA is expressed from a plasmid that also carries the Cas9 gene. By utilizing the Cas9-NE method, both the integrative and circular forms of MGEs can be precisely and completely eliminated through Cas9 cleavage, generating MGE-removed cells. We have successfully applied the Cas9-NE method to remove four representative MGEs, including plasmids, prophages, and genomic islands, from Vibrio strains. This new approach not only enables various investigations on MGEs but also has significant implications for the rapid generation of strains for commercial purposes.IMPORTANCEMobile genetic elements (MGEs) are of utmost importance for bacterial adaptation and pathogenicity, existing in various forms and multiple copies within bacterial cells. Integrated MGEs play dual roles in bacterial hosts, enhancing the fitness of the host by delivering cargo genes and potentially modifying the bacterial genome through the integration/excision process. This process can lead to alterations in promoters or coding sequences or even gene disruptions at integration sites, influencing the physiological functions of host bacteria. Here, we developed a new approach called Cas9-NE, allowing them to maintain the natural sequence changes associated with MGE excision. Cas9-NE allows the one-step removal of integrated and circular MGEs, addressing the challenge of eliminating various MGE forms efficiently. This approach simplifies MGE elimination in bacteria, expediting research on MGEs.

Active prophages in coral-associated Halomonas capable of lateral transduction

Temperate phages can interact with bacterial hosts through lytic and lysogenic cycles via different mechanisms. Lysogeny has been identified as the major form of bacteria-phage interaction in the coral-associated microbiome. However, the lysogenic-to-lytic switch of temperate phages in ecologically important coral-associated bacteria and its ecological impact have not been extensively investigated. By studying the prophages in coral-associated Halomonas meridiana, we found that two prophages, Phm1 and Phm3, are inducible by the DNA-damaging agent mitomycin C and that Phm3 is spontaneously activated under normal cultivation conditions. Furthermore, Phm3 undergoes an atypical lytic pathway that can amplify and package adjacent host DNA, potentially resulting in lateral transduction. The induction of Phm3 triggered a process of cell lysis accompanied by the formation of outer membrane vesicles (OMVs) and Phm3 attached to OMVs. This unique cell-lysis process was controlled by a four-gene lytic module within Phm3. Further analysis of the Tara Ocean dataset revealed that Phm3 represents a new group of temperate phages that are widely distributed and transcriptionally active in the ocean. Therefore, the combination of lateral transduction mediated by temperate phages and OMV transmission offers a versatile strategy for host-phage coevolution in marine ecosystems.

A universe of human gut-derived bacterial prophages: unveiling the hidden viral players in intestinal microecology

Prophages, which are an existing form of temperate phages that integrate into host genomes, have been found extensively present in diverse bacterial species. The human gut microbiome, characterized by its complexity, dynamism, and interconnectivity among multiple species, remains inadequately understood in terms of the global landscape of bacterial prophages. Here, we analyzed 43,942 human gut-derived bacterial genomes (439 species of 12 phyla) and identified 105,613 prophage regions in ~ 92% of them. 16254 complete prophages were distributed in ~ 24% of bacteria, indicating an extremely uneven prophage distribution across various species within the human gut. Among all identified prophages, ~4% possessed cross-genera (2-20 genera) integration capacity, while ~ 17% displayed broad infection host ranges (targeting 2-35 genera). Functional gene annotation revealed that antibiotic-resistance genes and toxin-related genes were detected in ~ 2.5% and ~ 5.8% of all prophages, respectively. Furthermore, through sequence alignments between our obtained prophages and publicly available human gut phageome contigs, we have observed that ~ 72% of non-redundant prophages are previously unreported genomes, and they illuminate ~ 6.5-9.5% of the individual intestinal "viral dark matter". Our study represents the first comprehensive depiction of human gut-derived prophages, provides a substantial collection of reference sequences for forthcoming human gut phageome-related investigations, and potentially enables better risk assessment of prophage dissemination.

Widespread and largely unknown prophage activity, diversity, and function in two genera of wheat phyllosphere bacteria

Environmental bacteria host an enormous number of prophages, but their diversity and natural functions remain largely elusive. Here, we investigate prophage activity and diversity in 63 Erwinia and Pseudomonas strains isolated from flag leaves of wheat grown in a single field. Introducing and validating Virion Induction Profiling Sequencing (VIP-Seq), we identify and quantify the activity of 120 spontaneously induced prophages, discovering that some phyllosphere bacteria produce more than 108 virions/mL in overnight cultures, with significant induction also observed in planta. Sequence analyses and plaque assays reveal E. aphidicola prophages contribute a majority of intraspecies genetic diversity and divide their bacterial hosts into antagonistic factions engaged in widespread microbial warfare, revealing the importance of prophage-mediated microdiversity. When comparing spontaneously active prophages with predicted prophages we also find insertion sequences are strongly correlated with non-active prophages. In conclusion, we discover widespread and largely unknown prophage diversity and function in phyllosphere bacteria.

VIBES: A Workflow for Annotating and Visualizing Viral Sequences Integrated into Bacterial Genomes

Bacteriophages are viruses that infect bacteria. Many bacteriophages integrate their genomes into the bacterial chromosome and become prophages. Prophages may substantially burden or benefit host bacteria fitness, acting in some cases as parasites and in others as mutualists, and have been demonstrated to increase host virulence. The increasing ease of bacterial genome sequencing provides an opportunity to deeply explore prophage prevalence and insertion sites. Here we present VIBES, a workflow intended to automate prophage annotation in complete bacterial genome sequences. VIBES provides additional context to prophage annotations by annotating bacterial genes and viral proteins in user-provided bacterial and viral genomes. The VIBES pipeline is implemented as a Nextflow-driven workflow, providing a simple, unified interface for execution on local, cluster, and cloud computing environments. For each step of the pipeline, a container including all necessary software dependencies is provided. VIBES produces results in simple tab separated format and generates intuitive and interactive visualizations for data exploration. Despite VIBES' primary emphasis on prophage annotation, its generic alignment-based design allows it to be deployed as a general-purpose sequence similarity search manager. We demonstrate the utility of the VIBES prophage annotation workflow by searching for 178 Pf phage genomes across 1,072 Pseudomonas spp. genomes. VIBES software is available at https://github.com/TravisWheelerLab/VIBES.

A systematic analysis of marine lysogens and proviruses

Viruses are ubiquitous in the oceans, exhibiting high abundance and diversity. Here, we systematically analyze existing genomic sequences of marine prokaryotes to compile a Marine Prokaryotic Genome Dataset (MPGD, consisting of over 12,000 bacterial and archaeal genomes) and a Marine Temperate Viral Genome Dataset (MTVGD). At least 40% of the MPGD genomes contain one or more proviral sequences, indicating that they are lysogens. The MTVGD includes over 12,900 viral contigs or putative proviruses, clustered into 10,897 viral genera. We show that lysogens and proviruses are abundant in marine ecosystems, particularly in the deep sea, and marine lysogens differ from non-lysogens in multiple genomic features and growth properties. We reveal several virus-host interaction networks of potential ecological relevance, and identify proviruses that appear to be able to infect (or to be transferred between) different bacterial classes and phyla. Auxiliary metabolic genes in the MTVGD are enriched in functions related to carbohydrate metabolism. Finally, we experimentally demonstrate the impact of a prophage on the transcriptome of a representative marine Shewanella bacterium. Our work contributes to a better understanding of the ecology of marine prokaryotes and their viruses.

A multidrug-resistant Salmonella enterica Typhimurium DT104 complex lineage circulating among humans and cattle in the USA lost the ability to produce pertussis-like toxin ArtAB

Salmonella enterica subsp. enterica serotype Typhimurium definitive type 104 (DT104) can infect both humans and animals and is often multidrug-resistant (MDR). Previous studies have indicated that, unlike most S . Typhimurium, the overwhelming majority of DT104 strains produce pertussis-like toxin ArtAB via prophage-encoded genes artAB . However, DT104 that lack artAB have been described on occasion. Here, we identify an MDR DT104 complex lineage circulating among humans and cattle in the USA, which lacks artAB (i.e. the ‘U.S. artAB -negative major clade’; n =42 genomes). Unlike most other bovine- and human-associated DT104 complex strains from the USA (n =230 total genomes), which harbour artAB on prophage Gifsy-1 (n =177), members of the U.S. artAB -negative major clade lack Gifsy-1, as well as anti-inflammatory effector gogB . The U.S. artAB -negative major clade encompasses human- and cattle-associated strains isolated from ≥11 USA states over a 20-year period. The clade was predicted to have lost artAB , Gifsy-1 and gogB circa 1985–1987 (95 % highest posterior density interval 1979.0–1992.1). When compared to DT104 genomes from other regions of the world (n =752 total genomes), several additional, sporadic artAB , Gifsy-1 and/or gogB loss events among clades encompassing five or fewer genomes were observed. Using phenotypic assays that simulate conditions encountered during human and/or bovine digestion, members of the U.S. artAB -negative major clade did not differ from closely related Gifsy-1/artAB /gogB -harbouring U.S. DT104 complex strains (ANOVA raw P >0.05); thus, future research is needed to elucidate the roles that artAB , gogB and Gifsy-1 play in DT104 virulence in humans and animals.

Comparative Genomic Analysis of Cold-Water Coral-Derived Sulfitobacter faviae: Insights into Their Habitat Adaptation and Metabolism

Sulfitobacter is one of the major sulfite-oxidizing alphaproteobacterial groups and is often associated with marine algae and corals. Their association with the eukaryotic host cell may have important ecological contexts due to their complex lifestyle and metabolism. However, the role of Sulfitobacter in cold-water corals remains largely unexplored. In this study, we explored the metabolism and mobile genetic elements (MGEs) in two closely related Sulfitobacter faviae strains isolated from cold-water black corals at a depth of ~1000 m by comparative genomic analysis. The two strains shared high sequence similarity in chromosomes, including two megaplasmids and two prophages, while both contained several distinct MGEs, including prophages and megaplasmids. Additionally, several toxin-antitoxin systems and other types of antiphage elements were also identified in both strains, potentially helping Sulfitobacter faviae overcome the threat of diverse lytic phages. Furthermore, the two strains shared similar secondary metabolite biosynthetic gene clusters and genes involved in dimethylsulfoniopropionate (DMSP) degradation pathways. Our results provide insight into the adaptive strategy of Sulfitobacter strains to thrive in ecological niches such as cold-water corals at the genomic level.

Advances in the field of phage-based therapy with special emphasis on computational resources

In the current era, one of the major challenges is to manage the treatment of drug/antibiotic-resistant strains of bacteria. Phage therapy, a century-old technique, may serve as an alternative to antibiotics in treating bacterial infections caused by drug-resistant strains of bacteria. In this review, a systematic attempt has been made to summarize phage-based therapy in depth. This review has been divided into the following two sections: general information and computer-aided phage therapy (CAPT). In the case of general information, we cover the history of phage therapy, the mechanism of action, the status of phage-based products (approved and clinical trials) and the challenges. This review emphasizes CAPT, where we have covered primary phage-associated resources, phage prediction methods and pipelines. This review covers a wide range of databases and resources, including viral genomes and proteins, phage receptors, host genomes of phages, phage-host interactions and lytic proteins. In the post-genomic era, identifying the most suitable phage for lysing a drug-resistant strain of bacterium is crucial for developing alternate treatments for drug-resistant bacteria and this remains a challenging problem. Thus, we compile all phage-associated prediction methods that include the prediction of phages for a bacterial strain, the host for a phage and the identification of interacting phage-host pairs. Most of these methods have been developed using machine learning and deep learning techniques. This review also discussed recent advances in the field of CAPT, where we briefly describe computational tools available for predicting phage virions, the life cycle of phages and prophage identification. Finally, we describe phage-based therapy's advantages, challenges and opportunities.

The coral pathogen Vibrio coralliilyticus kills non-pathogenic holobiont competitors by triggering prophage induction

The coral reef microbiome is central to reef health and resilience. Competitive interactions between opportunistic coral pathogens and other commensal microbes affect the health of coral. Despite great advances over the years in sequencing-based microbial profiling of healthy and diseased coral, the molecular mechanism underlying colonization competition has been much less explored. In this study, by examining the culturable bacteria inhabiting the gastric cavity of healthy Galaxea fascicularis, a scleractinian coral, we found that temperate phages played a major role in mediating colonization competition in the coral microbiota. Specifically, the non-toxigenic Vibrio sp. inhabiting the healthy coral had a much higher colonization capacity than the coral pathogen Vibrio coralliilyticus, yet this advantage was diminished by the latter killing the former. Pathogen-encoded LodAB, which produces hydrogen peroxide, triggers the lytic cycle of prophage in the non-toxicogenic Vibrio sp. Importantly, V. coralliilyticus could outcompete other coral symbiotic bacteria (for example, Endozoicomonas sp.) through LodAB-dependent prophage induction. Overall, we reveal that LodAB can be used by pathogens as an important weapon to gain a competitive advantage over lysogenic competitors when colonizing corals.