Integrative omics analysis of Pseudomonas aeruginosa virus PA5oct highlights the molecular complexity of jumbo phages

Nucleus-forming jumbophage PhiKZ therapeutically outcompetes non-nucleus-forming jumbophage Callisto

With the recent resurgence of phage therapy in modern medicine, jumbophages are currently under the spotlight due to their numerous advantages as anti-infective agents. However, most significant discoveries to date have primarily focused on nucleus-forming jumbophages, not their non-nucleus-forming counterparts. In this study, we compare the biological characteristics exhibited by two genetically diverse jumbophages: 1) the well-studied nucleus-forming jumbophage, PhiKZ; and 2) the newly discovered non-nucleus-forming jumbophage, Callisto. Single-cell infection studies further show that Callisto possesses different replication machinery, resulting in a delay in phage maturation compared to that of PhiKZ. The therapeutic potency of both phages was examined in vitro and in vivo, demonstrating that PhiKZ holds certain superior characteristics over Callisto. This research sheds light on the importance of the subcellular infection machinery and the organized progeny maturation process, which could potentially provide valuable insight in the future development of jumbophage-based therapeutics.

Refining the transcriptional landscapes for distinct clades of virulent phages infecting Pseudomonas aeruginosa

The introduction of high-throughput sequencing has resulted in a surge of available bacteriophage genomes, unveiling their tremendous genomic diversity. However, our current understanding of the complex transcriptional mechanisms that dictate their gene expression during infection is limited to a handful of model phages. Here, we applied ONT-cappable-seq to reveal the transcriptional architecture of six different clades of virulent phages infecting Pseudomonas aeruginosa. This long-read microbial transcriptomics approach is tailored to globally map transcription start and termination sites, transcription units, and putative RNA-based regulators on dense phage genomes. Specifically, the full-length transcriptomes of LUZ19, LUZ24, 14-1, YuA, PAK_P3, and giant phage phiKZ during early, middle, and late infection were collectively charted. Beyond pinpointing traditional promoter and terminator elements and transcription units, these transcriptional profiles provide insights in transcriptional attenuation and splicing events and allow straightforward validation of Group I intron activity. In addition, ONT-cappable-seq data can guide genome-wide discovery of novel regulatory element candidates, including noncoding RNAs and riboswitches. This work substantially expands the number of annotated phage-encoded transcriptional elements identified to date, shedding light on the intricate and diverse gene expression regulation mechanisms in Pseudomonas phages, which can ultimately be sourced as tools for biotechnological applications in phage and bacterial engineering.

Phage-Host Interaction Analysis Using Flow Cytometry

Phage-host interactions are commonly evaluated by culture-based methods. However, these techniques are very laborious and time-consuming. Therefore, other time-efficient, not labor-intensive, and cost-effective methods have been developed.This chapter describes the methodology used to assess the susceptibility of planktonic cultures of bacteria to phage infection and to study their interactions over time by flow cytometry.

Obtaining Detailed Phage Transcriptomes Using ONT-Cappable-Seq

Detailed transcription maps of bacteriophages are not usually explored, limiting our understanding of molecular phage biology and restricting their exploitation and engineering. The ONT-cappable-seq method described here brings phage transcriptomics to the accessible nanopore sequencing platform and provides an affordable and more detailed overview of transcriptional features compared to traditional RNA-seq experiments. With ONT-cappable-seq, primary transcripts are specifically capped, enriched, and prepared for long-read sequencing on the nanopore sequencing platform. This enables end-to-end sequencing of unprocessed transcripts covering both phage and host genome, thus providing insight on their operons. The subsequent analysis pipeline makes it possible to rapidly identify the most important transcriptional features such as transcription start and stop sites. The obtained data can thus provide a comprehensive overview of the transcription by your phage of interest.

Comparative genomics and proteomics analysis of phages infecting multi-drug resistant Escherichia coli O177 isolated from cattle faeces

The increasing prevalence of antimicrobial-resistant (AMR) pathogens has become a major global health concern. To address this challenge, innovative strategies such as bacteriophage therapy must be optimised. Genomic characterisation is a crucial step in identifying suitable phage candidates for combating AMR pathogens. The aim of this study was to characterise seven phages that infect the Escherichia coli O177 strain using a whole genome sequencing. The analysis of genome sequences revealed that these phages had linear dsDNA, with genome sizes spanning from 136, 483 to 166,791 bp and GC content varying from 35.39 to 43.63%. Taxonomically, the phages were classified under three different subfamilies (Stephanstirmvirinae, Tevenvirinae, and Vequintavirinae) and three genera (Phapecoctavirus, Tequatrovirus, and Vequintavirus) within the class Caudoviricetes. In silico PhageAI analysis predicted that all the phages were virulent, with confidence levels between 96.07 and 97.26%. The phage genomes contained between 66 and 82 ORFs, which encode hypothetical and putative functional proteins. In addition, the phage genomes contained core genes associated with molecular processes such as DNA replication, transcription modulation, nucleotide metabolism, phage structure (capsid and tail), and lysis. None of the genomes carried genes associated with undesirable traits such as integrase, antimicrobial resistance, virulence, and toxins. The study revealed high genome and proteome homology among E. coli O177 phages and other known Escherichia phages. The results suggest that the seven phages are new members of the genera Phapecoctavirus, Tequatrovirus, and Vequintavirus under the subfamilies Stephanstirmvirinae, Tevenvirinae, and Vequintavirinae, respectively.

Characterization and Comparative Genomic Analysis of Three Virulent E. coli Bacteriophages with the Potential to Reduce Antibiotic-Resistant Bacteria in the Environment

The emerging global crisis of antibiotic resistance demands new alternative antibacterial solutions. Although bacteriophages have been used to combat bacterial infections for over a century, a dramatic boost in phage studies has recently been observed. In the development of modern phage applications, a scientific rationale is strongly required and newly isolated phages need to be examined in detail. In this study, we present the full characterization of bacteriophages BF9, BF15, and BF17, with lytic activity against extended-spectrum β-lactamases (ESBLs)- and AmpC β-lactamases (AmpC)-producing Escherichia coli, the prevalence of which has increased significantly in livestock in recent decades, representing a great hazard to food safety and a public health risk. Comparative genomic and phylogenetic analysis indicated that BF9, BF15, and BF17 represent the genera Dhillonvirus, Tequatrovirus, and Asteriusvirus, respectively. All three phages significantly reduced in vitro growth of their bacterial host and retained the ability to lyse bacteria after preincubation at wide ranges of temperature (−20–40 °C) and pH (5–9). The results described herein indicate the lytic nature of BF9, BF15, and BF17, which, along with the absence of genes encoding toxins and bacterial virulence factors, represents an undoubted asset in terms of future phage application.

Bacteriophages inhibit and evade cGAS-like immune function in bacteria

A fundamental strategy of eukaryotic antiviral immunity involves the cGAS enzyme, which synthesizes 2',3'-cGAMP and activates the effector STING. Diverse bacteria contain cGAS-like enzymes that produce cyclic oligonucleotides and induce anti-phage activity, known as CBASS. However, this activity has only been demonstrated through heterologous expression. Whether bacteria harboring CBASS antagonize and co-evolve with phages is unknown. Here, we identified an endogenous cGAS-like enzyme in Pseudomonas aeruginosa that generates 3',3'-cGAMP during phage infection, signals to a phospholipase effector, and limits phage replication. In response, phages express an anti-CBASS protein ("Acb2") that forms a hexamer with three 3',3'-cGAMP molecules and reduces phospholipase activity. Acb2 also binds to molecules produced by other bacterial cGAS-like enzymes (3',3'-cUU/UA/UG/AA) and mammalian cGAS (2',3'-cGAMP), suggesting broad inhibition of cGAS-based immunity. Upon Acb2 deletion, CBASS blocks lytic phage replication and lysogenic induction, but rare phages evade CBASS through major capsid gene mutations. Altogether, we demonstrate endogenous CBASS anti-phage function and strategies of CBASS inhibition and evasion.

Interactions between Jumbo Phage SA1 and Staphylococcus: A Global Transcriptomic Analysis

Staphylococcus aureus (S. aureus) is an important zoonotic pathogen that poses a serious health concern to humans and cattle worldwide. Although it has been proven that lytic phages may successfully kill S. aureus, the interaction between the host and the phage has yet to be thoroughly investigated, which will likely limit the clinical application of phage. Here, RNA sequencing (RNA-seq) was used to examine the transcriptomics of jumbo phage SA1 and Staphylococcus JTB1-3 during a high multiplicity of infection (MOI) and RT-qPCR was used to confirm the results. The RNA-seq analysis revealed that phage SA1 took over the transcriptional resources of the host cells and that the genes were categorized as early, middle, and late, based on the expression levels during infection. A minor portion of the resources of the host was employed to enable phage replication after infection because only 35.73% (997/2790) of the host genes were identified as differentially expressed genes (DEGs). Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses showed that the phage infection mainly affected the nucleotide metabolism, protein metabolism, and energy-related metabolism of the host. Moreover, the expression of the host genes involved in anti-phage systems, virulence, and drug resistance significantly changed during infection. This research gives a fresh understanding of the relationship between jumbo phages and their Gram-positive bacteria hosts and provides a reference for studying phage treatment and antibiotics.

Global Transcriptomic Analysis of Bacteriophage-Host Interactions between a Kayvirus Therapeutic Phage and Staphylococcus aureus

Kayviruses are polyvalent broad host range staphylococcal phages with a potential to combat staphylococcal infections. However, the implementation of rational phage therapy in medicine requires a thorough understanding of the interactions between bacteriophages and pathogens at omics level. To evaluate the effect of a phage used in therapy on its host bacterium, we performed differential transcriptomic analysis by RNA-Seq from bacteriophage K of genus Kayvirus infecting two Staphylococcus aureus strains, prophage-less strain SH1000 and quadruple lysogenic strain Newman. The temporal transcriptional profile of phage K was comparable in both strains except for a few loci encoding hypothetical proteins. Stranded sequencing revealed transcription of phage noncoding RNAs that may play a role in the regulation of phage and host gene expression. The transcriptional response of S. aureus to phage K infection resembles a general stress response with differential expression of genes involved in a DNA damage response. The host transcriptional changes involved upregulation of nucleotide, amino acid and energy synthesis and transporter genes and downregulation of host transcription factors. The interaction of phage K with variable genetic elements of the host showed slight upregulation of gene expression of prophage integrases and antirepressors. The virulence genes involved in adhesion and immune evasion were only marginally affected, making phage K suitable for therapy. IMPORTANCE Bacterium Staphylococcus aureus is a common human and veterinary pathogen that causes mild to life-threatening infections. As strains of S. aureus are becoming increasingly resistant to multiple antibiotics, the need to search for new therapeutics is urgent. A promising alternative to antibiotic treatment of staphylococcal infections is a phage therapy using lytic phages from the genus Kayvirus. Here, we present a comprehensive view on the phage-bacterium interactions on transcriptomic level that improves the knowledge of molecular mechanisms underlying the Kayvirus lytic action. The results will ensure safer usage of the phage therapeutics and may also serve as a basis for the development of new antibacterial strategies.

Characterization and Genomic Analysis of a Novel Jumbo Bacteriophage vB_StaM_SA1 Infecting Staphylococcus aureus With Two Lysins

The development of new antimicrobial agents is critically needed due to the alarming increase in antibiotic resistance in bacterial pathogens. Phages have been widely considered as effective alternatives to antibiotics. A novel phage vB_StaM_SA1 (hereinafter as SA1) that can infect multiple Staphylococcus strains was isolated from untreated sewage of a pig farm, which belonged to Myoviridae family. At MOI of 0.1, the latent period of phage SA1 was 55 min, and the final titer reached about 10⁹ PFU/mL. The genome of phage SA1 was 260,727 bp, indicating that it can be classified as a jumbo phage. The genome of SA1 had 258 ORFs and a serine tRNA, while only 53 ORFs were annotated with functions. Phage SA1 contained a group of core genes that was characterized by multiple RNA polymerase subunits and also found in phiKZ-related jumbo phages. The phylogenetic tree showed that phage SA1 was a phiKZ-related phage and was closer to jumbo phages compared with Staphylococcus phages with small genome. Three proteins (lys4, lys210, and lys211) were predicted to be associated with lysins, and two proteins with lytic function were verified by recombinant expression and bacterial survival test. Both lys210 and lys211 possessed efficient bactericidal ability, and lys210 could lyse all test strains. The results show that phage SA1 and lys210/lys211 could be potentially used as antibiotic agents to treat Staphylococcus infection.

Phage-Host Interaction Analysis by Flow Cytometry Allows for Rapid and Efficient Screening of Phages

Recently, phages have become popular as an alternative to antibiotics. This increased demand for phage therapy needs rapid and efficient methods to screen phages infecting specific hosts. Existing methods are time-consuming, and for clinical purposes, novel, quick, and reliable screening methods are highly needed. Flow cytometry (FC) allows a quick differentiation and enumeration of bacterial cell populations and has been used to assess in vitro the activity of antimicrobial compounds. In this work, we propose FC as a rapid and reliable method to assess the susceptibility of a bacterial population to phage infection. For that, the interaction of phages vB_PaeM_CEB_DP1 and vB_PaeP_PE3 with Pseudomonas aeruginosa PAO1 was characterized by FC. Synchronous infection assays were performed, and samples were collected at different time points and stained with SYTO BC and PI before analysis. Part of the collected samples was used to characterize the expression of early, middle, and late genes by qPCR. Both FC and qPCR results were correlated with phage propagation assays. Results showed that SYTO BC median fluorescence intensity (MFI) values increased in the first 25 min of PE3 and DP1 infection. The increase of fluorescence is due to the expression of phage genes observed by qPCR. Since SYTO BC MFI values increase with gene expression, it allows the determination of host susceptibility to a phage in a short period of time, avoiding false positives caused by lysis from without. In conclusion, this method may allow for a quick and high-throughput real-time screening of different phages to a specific host, which can be crucial for a quick phage selection in clinical practice.

Shopping for phages? Unpacking design rules for therapeutic phage cocktails

In bacteriophage therapy, the combination of different phages into a single cocktail is of critical importance to overcome the narrow host range of single phage isolates. Today, the design of therapeutic cocktails is often akin to a black box and relies largely on intuition and (pre-)availability of isolates in local collections. Here we show that straightforward host range analysis can disclose design rules and we propose to apply/translate a data mining approach, historically applied in the field of marketing ('shopping cart analysis') to explore patterns in phage combinations. The technique is broadly applicable to host range datasets and can serve in combination with other molecular-based approaches to propose rationales for phage cocktail design.

Digital phagograms: predicting phage infectivity through a multilayer machine learning approach

Machine learning has been broadly implemented to investigate biological systems. In this regard, the field of phage biology has embraced machine learning to elucidate and predict phage-host interactions, based on receptor-binding proteins, (anti-)defense systems, prophage detection, and life cycle recognition. Here, we highlight the enormous potential of integrating information from omics data with insights from systems biology to better understand phage-host interactions. We conceptualize and discuss the potential of a multilayer model that mirrors the phage infection process, integrating adsorption, bacterial pan-immune components and hijacking of the bacterial metabolism to predict phage infectivity. In the future, this model can offer insights into the underlying mechanisms of the infection process, and digital phagograms can support phage cocktail design and phage engineering.

Type I CRISPR-Cas provides robust immunity but incomplete attenuation of phage-induced cellular stress

During infection, phages manipulate bacteria to redirect metabolism towards viral proliferation. To counteract phages, some bacteria employ CRISPR-Cas systems that provide adaptive immunity. While CRISPR-Cas mechanisms have been studied extensively, their effects on both the phage and the host during phage infection remains poorly understood. Here, we analysed the infection of Serratia by a siphovirus (JS26) and the transcriptomic response with, or without type I-E or I-F CRISPR-Cas immunity. In non-immune Serratia, phage infection altered bacterial metabolism by upregulating anaerobic respiration and amino acid biosynthesis genes, while flagella production was suppressed. Furthermore, phage proliferation required a late-expressed viral Cas4 homologue, which did not influence CRISPR adaptation. While type I-E and I-F immunity provided robust defence against phage infection, phage development still impacted the bacterial host. Moreover, DNA repair and SOS response pathways were upregulated during type I immunity. We also discovered that the type I-F system is controlled by a positive autoregulatory feedback loop that is activated upon phage targeting during type I-F immunity, leading to a controlled anti-phage response. Overall, our results provide new insight into phage-host dynamics and the impact of CRISPR immunity within the infected cell.

From Trees to Clouds: PhageClouds for Fast Comparison of ∼640,000 Phage Genomic Sequences and Host-Centric Visualization Using Genomic Network Graphs

Background: Fast and computationally efficient strategies are required to explore genomic relationships within an increasingly large and diverse phage sequence space. Here, we present PhageClouds, a novel approach using a graph database of phage genomic sequences and their intergenomic distances to explore the phage genomic sequence space. Methods: A total of 640,000 phage genomic sequences were retrieved from a variety of databases and public virome assemblies. Intergenomic distances were calculated with dashing, an alignment-free method suitable for handling massive data sets. These data were used to build a Neo4j^® graph database. Results: PhageClouds supported the search of related phages among all complete phage genomes from GenBank for a single query phage in just 10 s. Moreover, PhageClouds expanded the number of closely related phage sequences detected for both finished and draft phage genomes, in comparison with searches exclusively targeting phage entries from GenBank. Conclusions: PhageClouds is a novel resource that will facilitate the analysis of phage genomic sequences and the characterization of assembled phage genomes.

Emerging Aspects of Jumbo Bacteriophages

The bacteriophages have been explored at a huge scale as a model system for their applications in many biological-related fields. Jumbo phages with a large genome size from 200 to 500 kbp were not previously assigned a great value, and characterized by complex structures coupled with large virions with a wide variety of hosts. The origin of most of the jumbo phages was not well understood; however, many other prominent features have been discovered recently. In the current review, we strive to unearth the most advanced characteristics of jumbo phages, particularly their significance and structural organization that holds immense value to the viral life cycle. The unique characteristics of jumbo phages are the basis of variations in different types of phages concerning their organization at the genomic level, virion structure, evolution, and progeny propagation. The presence of tRNA and additional translation-related genes along with chaperonin genes mark the ability of these phages for being independent of host molecular machinery enabling them to have wide host options. A large number of jumbo phages have been isolated from various sources through advanced standard screening methods. The current review has summarized the available data on jumbo phages and discussed the genome orientation of jumbo phages, translational machinery, diversity and evolution of jumbo phages. In the studies conducted, jumbo phages possessed special additional genes that helps to reduce the dependence of jumbo phages on their hosts. Furthermore, their genomes might have evolved from smaller genome phages.

Unraveling Protein Interactions between the Temperate Virus Bam35 and Its Bacillus Host Using an Integrative Yeast Two Hybrid-High Throughput Sequencing Approach

Bacillus virus Bam35 is the model Betatectivirus and member of the family Tectiviridae, which is composed of tailless, icosahedral, and membrane-containing bacteriophages. Interest in these viruses has greatly increased in recent years as they are thought to be an evolutionary link between diverse groups of prokaryotic and eukaryotic viruses. Additionally, betatectiviruses infect bacteria of the Bacillus cereus group, which are known for their applications in industry and notorious since it contains many pathogens. Here, we present the first protein–protein interactions (PPIs) network for a tectivirus–host system by studying the Bam35–Bacillus thuringiensis model using a novel approach that integrates the traditional yeast two-hybrid system and high-throughput sequencing (Y2H-HTS). We generated and thoroughly analyzed a genomic library of Bam35′s host B. thuringiensis HER1410 and screened interactions with all the viral proteins using different combinations of bait–prey couples. Initial analysis of the raw data enabled the identification of over 4000 candidate interactions, which were sequentially filtered to produce 182 high-confidence interactions that were defined as part of the core virus–host interactome. Overall, host metabolism proteins and peptidases were particularly enriched within the detected interactions, distinguishing this host–phage system from the other reported host–phage PPIs. Our approach also suggested biological roles for several Bam35 proteins of unknown function, including the membrane structural protein P25, which may be a viral hub with a role in host membrane modification during viral particle morphogenesis. This work resulted in a better understanding of the Bam35–B. thuringiensis interaction at the molecular level and holds great potential for the generalization of the Y2H-HTS approach for other virus–host models.

A systematic review from basics to omics on bacteriophage applications in poultry production and processing

The growing human population is currently facing an unprecedented challenge on global food production and sustainability. Despite recognizing poultry as one of the most successful and rapidly growing food industries to address this challenge; poultry health and safety remain major issues that entail immediate attention. Bacterial diseases including colibacillosis, salmonellosis, and necrotic enteritis have become increasingly prevalent during poultry production. Likewise, outbreaks caused by consumption of undercooked poultry products contaminated with zoonotic bacterial pathogens such as Salmonella, Campylobacter and Listeria, are a serious public health concern. With antimicrobial resistance problem and restricted use of antibiotics in food producing animals, bacteriophages are increasingly recognized as an attractive natural antibacterial alternative. Bacteriophages have recently shown promising results to treat diseases in poultry, reduce contamination of carcasses, and enhance the safety of poultry products. Omics technologies have been successfully employed to accurately characterize bacteriophages and their genes/proteins important for interaction with bacterial hosts. In this review, the potential of using lytic bacteriophages to mitigate the risk of major poultry-associated bacterial pathogens are explored. This study also explores challenges associated with the adoption of this technology by industries. Furthermore, the impact of omics approaches on studying bacteriophages, their host interaction and applications is discussed.

Finding functional associations between prokaryotic virus orthologous groups: a proof of concept

Background

The field of viromics has greatly benefited from recent developments in metagenomics, with significant efforts focusing on viral discovery. However, functional annotation of the increasing number of viral genomes is lagging behind. This is highlighted by the degree of annotation of the protein clusters in the prokaryotic Virus Orthologous Groups (pVOGs) database, with 83% of its current 9518 pVOGs having an unknown function.

Results

In this study we describe a machine learning approach to explore potential functional associations between pVOGs. We measure seven genomic features and use them as input to a Random Forest classifier to predict protein–protein interactions between pairs of pVOGs. After systematic evaluation of the model’s performance on 10 different datasets, we obtained a predictor with a mean accuracy of 0.77 and Area Under Receiving Operation Characteristic (AUROC) score of 0.83. Its application to a set of 2,133,027 pVOG-pVOG interactions allowed us to predict 267,265 putative interactions with a reported probability greater than 0.65. At an expected false discovery rate of 0.27, we placed 95.6% of the previously unannotated pVOGs in a functional context, by predicting their interaction with a pVOG that is functionally annotated.

Conclusions

We believe that this proof-of-concept methodology, wrapped in a reproducible and automated workflow, can represent a significant step towards obtaining a more complete picture of bacteriophage biology.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04343-w.

Genome-driven elucidation of phage-host interplay and impact of phage resistance evolution on bacterial fitness

When considering the interactions between bacteriophages and their host, the issue of phage-resistance emergence is a key element in understanding the ecological impact of phages on the bacterial population. It is also an essential parameter for the implementation of phage therapy to combat antibiotic-resistant pathogens. This study investigates the phenotypic and genetic responses of five Pseudomonas aeruginosa strains (PAO1, A5803, AA43, CHA, and PAK) to the infection by seven phages with distinct evolutionary backgrounds and recognised receptors (LPS/T4P). Emerging phage-insensitivity was generally accompanied by self and cross-resistance mechanisms. Significant differences were observed between the reference PAO1 responses compared to other clinical representatives. LPS-dependent phage infections in clinical strains selected for mutations in the "global regulatory" and "other" genes, rather than in the LPS-synthesis clusters detected in PAO1 clones. Reduced fitness, as proxied by the growth rate, was correlated with large deletion (20-500 kbp) and phage carrier state. Multi-phage resistance was significantly correlated with a reduced growth rate but only in the PAO1 population. In addition, we observed that the presence of prophages decreased the lytic phage maintenance seemingly protecting the host against carrier state and occasional lytic phage propagation, thus preventing a significant reduction in bacterial growth rate.

Characterization and Genomic Analysis of PALS2, a Novel Staphylococcus Jumbo Bacteriophage

Staphylococcus aureus is an important human pathogen that can be frequently encountered in clinical and food-processing surroundings. Among the various countermeasures, bacteriophages have been considered to be promising alternatives to antibiotics. In this study, the bacteriophage PALS2 was isolated from bird feces, and the genomic and biological characteristics of this phage were investigated. PALS2 was determined to belong to the Myoviridae family and exhibited extended host inhibition that persisted for up to 24 h with repeated bursts of 12 plaque-forming units/cell. The complete genome of PALS2 measured 268,746 base pairs (bp), indicating that PALS2 could be classified as a jumbo phage. The PALS2 genome contained 279 ORFs and 1 tRNA covering asparagine, and the majority of predicted PALS2 genes encoded hypothetical proteins. Additional genes involved in DNA replication and repair, nucleotide metabolism, and genes encoding multisubunit RNA polymerase were identified in the PALS2 genome, which is a common feature of typical jumbo phages. Comparative genomic analysis indicated that PALS2 is a phiKZ-related virus and is more similar to typical jumbo phages than to staphylococcal phages. Additionally, the effective antimicrobial activities of phage PALS2 suggest its possible use as a biocontrol agent in various clinical and food processing environments.

Molecular Determinants of Antibiotic Resistance in the Costa Rican Pseudomonas aeruginosa AG1 by a Multi-omics Approach: A Review of 10 Years of Study

Pseudomonas aeruginosa AG1 (PaeAG1) is a Costa Rican strain that was isolated in 2010 in a major Hospital. This strain has resistance to multiple antibiotics such as β-lactams (including carbapenems), aminoglycosides, and fluoroquinolones. PaeAG1 is considered critical (Priority 1) due to its resistance to carbapenems, and it was the first report of a P. aeruginosa isolate carrying both VIM-2 and IMP-18 genes encoding for metallo-β-lactamases (MBL) enzymes (both with carbapenemase activity). Owing to these traits, we have studied this model for 10 years using diverse approaches including multi-omics. In this review, we summarize the main points of the different steps that we have studied in PaeAG1: preliminary analyses of this strain at the genomic and phenomic levels revealed that this microorganism has particular features of antibiotic resistance. In the multi-omics approach, the genome assembly was the initial step to identify the genomic determinants of this strain, including virulence factors, antibiotic resistance genes, as well as a complex accessory genome. Second, a comparative genomic approach was implemented to define and update the phylogenetic relationship among complete P. aeruginosa genomes, the genomic island content in other strains, and the architecture of the two MBL-carrying integrons. Third, the proteomic profile of PaeAG1 was studied after exposure to antibiotics using 2-dimensional gel electrophoresis (2D-GE). Fourth, to study the central response to multiple perturbations in P. aeruginosa, i.e., the core perturbome, a machine learning approach was used. The analysis revealed biological functions and determinants that are shared by different disturbances. Finally, to evaluate the effects of ciprofloxacin (CIP) on PaeAG1, a growth curve comparison, differential expression analysis (RNA-Seq), and network analysis were performed. Using the results of the core perturbome (pathways that also were found in this perturbation with CIP), it was possible to identify the “exclusive” response and determinants of PaeAG1 after exposure to CIP. Altogether, after a decade of study using a multi-omics approach (at genomics, comparative genomics, perturbomics, transcriptomics, proteomics, and phenomics levels), we have provided new insights about the genomic and transcriptomic determinants associated with antibiotic resistance in PaeAG1. These results not only partially explain the high-risk condition of this strain that enables it to conquer nosocomial environments and its multi-resistance profile, but also this information may eventually be used as part of the strategies to fight this pathogen.

Intracellular Organization by Jumbo Bacteriophages

Since their discovery more than 100 years ago, the viruses that infect bacteria (bacteriophages) have been widely studied as model systems. Largely overlooked, however, have been "jumbo phages," with genome sizes ranging from 200 to 500 kbp. Jumbo phages generally have large virions with complex structures and a broad host spectrum. While the majority of jumbo phage genes are poorly functionally characterized, recent work has discovered many unique biological features, including a conserved tubulin homolog that coordinates a proteinaceous nucleus-like compartment that houses and segregates phage DNA. The tubulin spindle displays dynamic instability and centers the phage nucleus within the bacterial host during phage infection for optimal reproduction. The shell provides robust physical protection for the enclosed phage genomes against attack from DNA-targeting bacterial immune systems, thereby endowing jumbo phages with broad resistance. In this review, we focus on the current knowledge of the cytoskeletal elements and the specialized nuclear compartment derived from jumbo phages, and we highlight their importance in facilitating spatial and temporal organization over the viral life cycle. Additionally, we discuss the evolutionary relationships between jumbo phages and eukaryotic viruses, as well as the therapeutic potential and drawbacks of jumbo phages as antimicrobial agents in phage therapy.

Looking into the future of phage-based control of zoonotic pathogens in food and animal production

Using bacteriophages (phages) to control zoonotic pathogens in food and animal production is a realistic and promising antimicrobial approach. Recent studies have demonstrated their efficacy and safety, yet bringing phage products on the market remains a challenge. Here we summarize the procedure for advancing phage applications from the laboratory to simplified model systems and testing in pilot scale, to farms and food industries. We highlight the most important contributions concerning phages in food matrices and animal guts, and propose directions for future research required to understand interactions in such complex systems. Finally, we propose a holistic approach combining a data repository with modelling, multi-omic techniques and data analysis to modernize phage-based control of zoonotic pathogens.

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 Family of Viral Satellites Manipulates Invading Virus Gene Expression and Can Affect Cholera Toxin Mobilization

Many viruses possess temporally unfolding gene expression patterns aimed at subverting host defenses, commandeering host metabolism, and ultimately producing a large number of progeny virions. High-throughput omics tools, such as RNA sequencing (RNA-seq), have dramatically enhanced the resolution of expression patterns during infection. Less studied have been viral satellites, mobile genomes that parasitize viruses. By performing RNA-seq on infection time courses, we have obtained the first time-resolved transcriptomes for bacteriophage satellites during lytic infection. Specifically, we have acquired transcriptomes for the lytic Vibrio cholerae phage ICP1 and all five known variants of ICP1's parasite, the phage inducible chromosomal island-like elements (PLEs). PLEs rely on ICP1 for both DNA replication and mobilization and abolish production of ICP1 progeny in infected cells. We investigated PLEs' impact on ICP1 gene expression and found that PLEs did not broadly restrict or reduce ICP1 gene expression. A major exception occurred in ICP1's capsid morphogenesis operon, which was downregulated by each of the PLE variants. Surprisingly, PLEs were also found to alter the gene expression of CTXΦ, the integrative phage that encodes cholera toxin and is necessary for virulence of toxigenic V. cholerae One PLE, PLE1, upregulated CTXΦ genes involved in replication and integration and boosted CTXΦ mobility following induction of the SOS response.IMPORTANCE Viral satellites are found in all domains of life and can have profound fitness effects on both the viruses they parasitize and the cells they reside in. In this study, we have acquired the first RNA sequencing (RNA-seq) transcriptomes of viral satellites outside plants, as well as the transcriptome of the phage ICP1, a predominant predator of pandemic Vibrio cholerae Capsid downregulation, previously observed in an unrelated phage satellite, is conserved among phage inducible chromosomal island-like elements (PLEs), suggesting that viral satellites are under strong selective pressure to reduce the capsid expression of their larger host viruses. Despite conserved manipulation of capsid expression, PLEs exhibit divergent effects on CTXΦ transcription and mobility. Our results demonstrate that PLEs can influence both their hosts' resistance to phage and the mobility of virulence-encoding elements, suggesting that PLEs can play a substantial role in shaping Vibrio cholerae evolution.

Transcriptomic Analysis Reveals the Dependency of Pseudomonas aeruginosa Genes for Double-Stranded RNA Bacteriophage phiYY Infection Cycle

Bacteriophage phiYY is currently the only double-stranded RNA (dsRNA) phage that infects Pseudomonas aeruginosa and is a potential candidate for phage therapy. Here we applied RNA-seq to investigate the lytic cycle of phiYY infecting P. aeruginosa strain PAO1r. About 12.45% (651/5,229) of the host genes were determined to be differentially expressed genes. Moreover, oxidative stress response genes katB and ahpB are upregulated 64- to 128-fold after phage infection, and the single deletion of each gene blocked phiYY infection, indicating that phiYY is extremely sensitive to oxidative stress. On the contrary, another upregulated gene PA0800 might constrain phage infection, because the deletion of PA0800 resulted in a 3.5-fold increase of the efficiency of plating. Our study highlights a complicated dsRNA phage-phage global interaction and raises new questions toward the host defense mechanisms against dsRNA phage and dsRNA phage-encoded hijacking mechanisms.

Rethinking Phage Ecology by Rooting it Within an Established Plant Framework

Despite the abundance and significance of bacteriophages to microbial ecosystems, no broad ecological frameworks exist within which to determine "bacteriophage types" that reflect their ecological strategies and ways in which they interact with bacterial cells. To address this, we repurposed the well-established Grime's triangular CSR framework, which classifies plants according to three axes: competitiveness (C), ability to tolerate stress (S), and capacity to cope with disturbance (R). This framework is distinguished from other accepted schemes, as it seeks to identify individual characteristics of plants to understand their biological strategies and roles within an ecosystem. Our repurposing of the CSR triangle is based on phage transcription and the observation that typically phages have three major distinguishable transcription phases: early, middle, and late. We hypothesize that the proportion of genes expressed in these phases reflects key information about the phage "ecological strategy," namely the C, S, and R strategies, allowing us to examine phages in a similar way to how plants are projected onto the triangle. In the "phage version" of this scheme, we suggest: (1) that some phages prioritize the early phase of transcription that shuts off host defense mechanisms, which reflects competitiveness; (2) other phages prioritize tuning resource management mechanisms in the cell such as nucleotide metabolism during their "mid" expression profile to tolerate stress; and (3) a further subset of phages (termed Ruderals) survive disturbance by investing significant resources into regeneration so they express a higher proportion of their genes during late infection. We examined 42 published phage transcriptomes and show that they fall into discrete CSR categories according to their expression profiles. We discuss these positions in the context of their biology, which is largely consistent with our predictions of specific phage characteristics. In this opinion article, we suggest a starting point to ascribe phages into different functional types and thus understand them in an ecological framework. We suggest that this may have far-reaching implications for the application of phages in therapy and their exploitation to manipulate bacterial communities. We invite further use of this framework via our online tool; www.PhageCSR.ml.