Next-Generation "-omics" Approaches Reveal a Massive Alteration of Host RNA Metabolism during Bacteriophage Infection of Pseudomonas aeruginosa

Profiling the interplay and coevolution of Microcystis aeruginosa and cyanosiphophage Mic1

The cyanosiphophage Mic1 specifically infects the bloom-forming Microcystis aeruginosa FACHB 1339 from Lake Chaohu, China. Previous genomic analysis showed that its 92,627 bp double-stranded DNA genome consists of 98 putative open reading frames, 63% of which are of unknown function. Here, we investigated the transcriptome dynamics of Mic1 and its host using RNA sequencing. In the early, middle, and late phases of the 10 h lytic cycle, the Mic1 genes are sequentially expressed and could be further temporally grouped into two distinct clusters in each phase. Notably, six early genes, including gp49 that encodes a TnpB-like transposase, immediately reach the highest transcriptional level in half an hour, representing a pioneer cluster that rapidly regulates and redirects host metabolism toward the phage. An in-depth analysis of the host transcriptomic profile in response to Mic1 infection revealed significant upregulation of a polyketide synthase pathway and a type III-B CRISPR system, accompanied by moderate downregulation of the photosynthesis and key metabolism pathways. The constant increase of phage transcripts and relatively low replacement rate over the host transcripts indicated that Mic1 utilizes a unique strategy to gradually take over a small portion of host metabolism pathways after infection. In addition, genomic analysis of a less-infective Mic1 and a Mic1-resistant host strain further confirmed their dynamic interplay and coevolution via the frequent horizontal gene transfer. These findings provide insights into the mutual benefit and symbiosis of the highly polymorphic cyanobacteria M. aeruginosa and cyanophages.

IMPORTANCE

The highly polymorphic Microcystis aeruginosa is one of the predominant bloom-forming cyanobacteria in eutrophic freshwater bodies and is infected by diverse and abundant cyanophages. The presence of a large number of defense systems in M. aeruginosa genome suggests a dynamic interplay and coevolution with the cyanophages. In this study, we investigated the temporal gene expression pattern of Mic1 after infection and the corresponding transcriptional responses of its host. Moreover, the identification of a less-infective Mic1 and a Mic1-resistant host strain provided the evolved genes in the phage-host coevolution during the multiple-generation cultivation in the laboratory. Our findings enrich the knowledge on the interplay and coevolution of M. aeruginosa and its cyanophages and lay the foundation for the future application of cyanophage as a potential eco-friendly and bio-safe agent in controlling the succession of harmful cyanobacterial blooms.

Transcriptional dynamics during Rhodococcus erythropolis infection with phage WC1

BACKGROUND

Belonging to the Actinobacteria phylum, members of the Rhodococcus genus thrive in soil, water, and even intracellularly. While most species are non-pathogenic, several cause respiratory disease in animals and, more rarely, in humans. Over 100 phages that infect Rhodococcus species have been isolated but despite their importance for Rhodococcus ecology and biotechnology applications, little is known regarding the molecular genetic interactions between phage and host during infection. To address this need, we report RNA-Seq analysis of a novel Rhodococcus erythopolis phage, WC1, analyzing both the phage and host transcriptome at various stages throughout the infection process.

RESULTS

By five minutes post-infection WC1 showed upregulation of a CAS-4 family exonuclease, putative immunity repressor, an anti-restriction protein, while the host showed strong upregulation of DNA replication, SOS repair, and ribosomal protein genes. By 30 min post-infection, WC1 DNA synthesis genes were strongly upregulated while the host showed increased expression of transcriptional and translational machinery and downregulation of genes involved in carbon, energy, and lipid metabolism pathways. By 60 min WC1 strongly upregulated structural genes while the host showed a dramatic disruption of metal ion homeostasis. There was significant expression of both host and phage non-coding genes at all time points. While host gene expression declined over the course of infection, our results indicate that phage may exert more selective control, preserving the host's regulatory mechanisms to create an environment conducive for virion production.

CONCLUSIONS

The Rhodococcus genus is well recognized for its ability to synthesize valuable compounds, particularly steroids, as well as its capacity to degrade a wide range of harmful environmental pollutants. A detailed understanding of these phage-host interactions and gene expression is not only essential for understanding the ecology of this important genus, but will also facilitate development of phage-mediated strategies for bioremediation as well as biocontrol in industrial processes and biomedical applications. Given the current lack of detailed global gene expression studies on any Rhodococcus species, our study addresses a pressing need to identify tools and genes, such as F6 and rpf, that can enhance the capacity of Rhodococcus species for bioremediation, biosynthesis and pathogen control.

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.

Exploring the transcriptional landscape of phage-host interactions using novel high-throughput approaches

In the last decade, powerful high-throughput sequencing approaches have emerged to analyse microbial transcriptomes at a global scale. However, to date, applications of these approaches to microbial viruses such as phages remain scarce. Tailoring these techniques to virus-infected bacteria promises to obtain a detailed picture of the underexplored RNA biology and molecular processes during infection. In addition, transcriptome study of stress and perturbations induced by phages in their infected bacterial hosts is likely to reveal new fundamental mechanisms of bacterial metabolism and gene regulation. Here, we provide references and blueprints to implement emerging transcriptomic approaches towards addressing transcriptome architecture, RNA-RNA and RNA-protein interactions, RNA modifications, structures and heterogeneity of transcription profiles in infected cells that will provide guides for future directions in phage-centric therapeutic applications and microbial synthetic biology.

Coordinated transcriptional response to environmental stress by a Synechococcus virus

Viruses are a major control on populations of microbes. Often, their virulence is examined in controlled laboratory conditions. Yet, in nature, environmental conditions lead to changes in host physiology and fitness that may impart both costs and benefits on viral success. Phosphorus (P) is a major abiotic control on the marine cyanobacterium Synechococcus. Some viruses infecting Synechococcus have acquired, from their host, a gene encoding a P substrate binding protein (PstS), thought to improve virus replication under phosphate starvation. Yet, pstS is uncommon among cyanobacterial viruses. Thus, we asked how infections with viruses lacking PstS are affected by P scarcity. We show that the production of infectious virus particles of such viruses is reduced in low P conditions. However, this reduction in progeny is not caused by impaired phage genome replication, thought to be a major sink for cellular phosphate. Instead, transcriptomic analysis showed that under low P conditions, a PstS-lacking cyanophage increased the expression of a specific gene set that included mazG, hli2, and gp43 encoding a pyrophosphatase, a high-light inducible protein and DNA polymerase, respectively. Moreover, several of the upregulated genes were controlled by the host's phoBR two-component system. We hypothesize that recycling and polymerization of nucleotides liberates free phosphate and thus allows viral morphogenesis, albeit at lower rates than when phosphate is replete or when phages encode pstS. Altogether, our data show how phage genomes, lacking obvious P-stress-related genes, have evolved to exploit their host's environmental sensing mechanisms to coordinate their own gene expression in response to resource limitation.

Anaerobiosis, a neglected factor in phage-bacteria interactions

IMPORTANCE

Many parameters affect phage-bacteria interaction. Some of these parameters depend on the environment in which the bacteria are present. Anaerobiosis effect on phage infection in facultative anaerobic bacteria has not yet been studied. The absence of oxygen triggers metabolic changes in facultative bacteria and this affects phage infection and viral life cycle. Understanding how an anaerobic environment can alter the behavior of phages during infection is relevant for the phage therapy success.

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

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

In Situ Spatiotemporal SERS Measurements and Multivariate Analysis of Virally Infected Bacterial Biofilms Using Nanolaminated Plasmonic Crystals

In situ spatiotemporal biochemical characterization of the activity of living multicellular biofilms under external stimuli remains a significant challenge. Surface-enhanced Raman spectroscopy (SERS), combining the molecular fingerprint specificity of vibrational spectroscopy with the hotspot sensitivity of plasmonic nanostructures, has emerged as a promising noninvasive bioanalysis technique for living systems. However, most SERS devices do not allow reliable long-term spatiotemporal SERS measurements of multicellular systems because of challenges in producing spatially uniform and mechanically stable SERS hotspot arrays to interface with large cellular networks. Furthermore, very few studies have been conducted for multivariable analysis of spatiotemporal SERS datasets to extract spatially and temporally correlated biological information from multicellular systems. Here, we demonstrate in situ label-free spatiotemporal SERS measurements and multivariate analysis of Pseudomonas syringae biofilms during development and upon infection by bacteriophage virus Phi6 by employing nanolaminate plasmonic crystal SERS devices to interface mechanically stable, uniform, and spatially dense hotspot arrays with the P. syringae biofilms. We exploited unsupervised multivariate machine learning methods, including principal component analysis (PCA) and hierarchical cluster analysis (HCA), to resolve the spatiotemporal evolution and Phi6 dose-dependent changes of major Raman peaks originating from biochemical components in P. syringae biofilms, including cellular components, extracellular polymeric substances (EPS), metabolite molecules, and cell lysate-enriched extracellular media. We then employed supervised multivariate analysis using linear discriminant analysis (LDA) for the multiclass classification of Phi6 dose-dependent biofilm responses, demonstrating the potential for viral infection diagnosis. We envision extending the in situ spatiotemporal SERS method to monitor dynamic, heterogeneous interactions between viruses and bacterial networks for applications such as phage-based anti-biofilm therapy development and continuous pathogenic virus detection.

What Lies Beneath? Taking the Plunge into the Murky Waters of Phage Biology

The sequence revolution revealed that bacteria-infecting viruses, known as phages, are Earth's most abundant biological entities. Phages have far-reaching impacts on the form and function of microbial communities and play a fundamental role in ecological processes. However, even well into the sequencing revolution, we have only just begun to explore the murky waters around the phage biology iceberg. Many viral reads cannot be assigned to a culturable isolate, and reference databases are biased toward more easily collectible samples, which likely distorts our conclusions. This minireview points out alternatives to mapping reads to reference databases and highlights innovative bioinformatic and experimental approaches that can help us overcome some of the challenges in phage research and better decipher the impact of phages on microbial communities. Moving beyond the identification of novel phages, we highlight phage metabolomics as an important influencer of bacterial host cell physiology and hope to inspire the reader to consider the effects of phages on host metabolism and ecosystems at large. We encourage researchers to report unassigned/unknown sequencing reads and contigs and to continue developing alternative methods to investigate phages within sequence data.

The RNA-RNA interactome between a phage and its satellite virus reveals a small RNA that differentially regulates gene expression across both genomes

Satellite viruses are present across all domains of life, defined as subviral parasites that require infection by another virus for satellite progeny production. Phage satellites exhibit various regulatory mechanisms to manipulate phage gene expression to the benefit of the satellite, redirecting resources from the phage to the satellite, and often inhibiting phage progeny production. While small RNAs (sRNAs) are well documented as regulators of prokaryotic gene expression, they have not been shown to play a regulatory role in satellite-phage conflicts. Vibrio cholerae encodes the phage inducible chromosomal island-like element (PLE), a phage satellite, to defend itself against the lytic phage ICP1. Here, we use Hi-GRIL-seq to identify a complex RNA-RNA interactome between PLE and ICP1. Both inter- and intragenome RNA interactions were detected, headlined by the PLE sRNA, SviR. SviR is involved in regulating both PLE and ICP1 gene expression uniquely, decreasing ICP1 target translation and affecting PLE transcripts. The striking conservation of SviR across all known PLEs suggests the sRNA is deeply rooted in the PLE-ICP1 conflict and implicates sRNAs as unidentified regulators of gene expression in phage-satellite interactions.

"French Phage Network" Annual Conference-Seventh Meeting Report

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

Chemical-Biology and Metabolomics Studies in Phage-Host Interactions

For many years, several studies have explored the molecular mechanisms involved in the infection of bacteria by their specific phages to understand the main infection strategies and the host defense strategies. The modulation of the mechanisms involved in the infection, as well as the expression of key substances in the development of the different life cycles of phages, function as a natural source of strategies capable of promoting the control of different pathogens that are harmful to human and animal health. Therefore, this chapter aims to provide an overview of the mechanisms involved in virus-bacteria interaction to explore the main compounds produced or altered as a chemical survival strategy and the metabolism modulation when occurring a host-phage interaction. In this context, emphasis will be given to the chemistry of peptides/proteins and enzymes encoded by bacteriophages in the control of pathogenic bacteria and the use of secondary metabolites recently reported as active participants in the mechanisms of phage-bacteria interaction. Finally, metabolomics strategies developed to gain new insights into the metabolism involved in the phage-host interaction and the metabolomics workflow in host-phage interaction will be presented.

RNA-targeting CRISPR-Cas systems

CRISPR-Cas is a widespread adaptive immune system in bacteria and archaea that protects against viral infection by targeting specific invading nucleic acid sequences. Whereas some CRISPR-Cas systems sense and cleave viral DNA, type III and type VI CRISPR-Cas systems sense RNA that results from viral transcription and perhaps invasion by RNA viruses. The sequence-specific detection of viral RNA evokes a cell-wide response that typically involves global damage to halt the infection. How can one make sense of an immune strategy that encompasses broad, collateral effects rather than specific, targeted destruction? In this Review, we summarize the current understanding of RNA-targeting CRISPR-Cas systems. We detail the composition and properties of type III and type VI systems, outline the cellular defence processes that are instigated upon viral RNA sensing and describe the biological rationale behind the broad RNA-activated immune responses as an effective strategy to combat viral infection.

Protist impacts on marine cyanovirocell metabolism

The fate of oceanic carbon and nutrients depends on interactions between viruses, prokaryotes, and unicellular eukaryotes (protists) in a highly interconnected planktonic food web. To date, few controlled mechanistic studies of these interactions exist, and where they do, they are largely pairwise, focusing either on viral infection (i.e., virocells) or protist predation. Here we studied population-level responses of Synechococcus cyanobacterial virocells (i.e., cyanovirocells) to the protist Oxyrrhis marina using transcriptomics, endo- and exo-metabolomics, photosynthetic efficiency measurements, and microscopy. Protist presence had no measurable impact on Synechococcus transcripts or endometabolites. The cyanovirocells alone had a smaller intracellular transcriptional and metabolic response than cyanovirocells co-cultured with protists, displaying known patterns of virus-mediated metabolic reprogramming while releasing diverse exometabolites during infection. When protists were added, several exometabolites disappeared, suggesting microbial consumption. In addition, the intracellular cyanovirocell impact was largest, with 4.5- and 10-fold more host transcripts and endometabolites, respectively, responding to protists, especially those involved in resource and energy production. Physiologically, photosynthetic efficiency also increased, and together with the transcriptomics and metabolomics findings suggest that cyanovirocell metabolic demand is highest when protists are present. These data illustrate cyanovirocell responses to protist presence that are not yet considered when linking microbial physiology to global-scale biogeochemical processes.

Insights into the global freshwater virome

Viruses are by far the most abundant life forms on this planet. Yet, the full viral diversity remains mostly unknown, especially in environments like freshwater. Therefore, we aimed to study freshwater viruses in a global context. To this end, we downloaded 380 publicly available viral metagenomes (>1 TB). More than 60% of these metagenomes were discarded based on their levels of cellular contamination assessed by ribosomal DNA content. For the remaining metagenomes, assembled contigs were decontaminated using two consecutive steps, eventually yielding 273,365 viral contigs longer than 1,000 bp. Long enough contigs (≥ 10 kb) were clustered to identify novel genomes/genome fragments. We could recover 549 complete circular and high-quality draft genomes, out of which 10 were recognized as being novel. Functional annotation of these genomes showed that most of the annotated coding sequences are DNA metabolic genes or phage structural genes. On the other hand, taxonomic analysis of viral contigs showed that most of the assigned contigs belonged to the order Caudovirales, particularly the families of Siphoviridae, Myoviridae, and Podoviridae. The recovered viral contigs contained several auxiliary metabolic genes belonging to several metabolic pathways, especially carbohydrate and amino acid metabolism in addition to photosynthesis as well as hydrocarbon degradation and antibiotic resistance. Overall, we present here a set of prudently chosen viral contigs, which should not only help better understanding of freshwater viruses but also be a valuable resource for future virome studies.

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.

The function of small RNA in Pseudomonas aeruginosa

Pseudomonas aeruginosa, the main conditional pathogen causing nosocomial infection, is a gram-negative bacterium with the largest genome among the known bacteria. The main reasons why Pseudomonas aeruginosa is prone to drug-resistant strains in clinic are: the drug-resistant genes in its genome and the drug resistance easily induced by single antibiotic treatment. With the development of high-throughput sequencing technology and bioinformatics, the functions of various small RNAs (sRNA) in Pseudomonas aeruginosa are being revealed. Different sRNAs regulate gene expression by binding to protein or mRNA to play an important role in the complex regulatory network. In this article, first, the importance and biological functions of different sRNAs in Pseudomonas aeruginosa are explored, and then the evidence and possibilities that sRNAs served as drug therapeutic targets are discussed, which may introduce new directions to develop novel disease treatment strategies.

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.

The gut environment regulates bacterial gene expression which modulates susceptibility to bacteriophage infection

Abundance and diversity of bacteria and their viral predators, bacteriophages (phages), in the digestive tract are associated with human health. Particularly intriguing is the long-term coexistence of these two antagonistic populations. We performed genome-wide RNA sequencing on a human enteroaggregative Escherichia coli isolate to identify genes differentially expressed between in vitro conditions and in murine intestines. We experimentally demonstrated that four of these differentially expressed genes modified the interactions between E. coli and three virulent phages by either increasing or decreasing its susceptibility/resistance pattern and also by interfering with biofilm formation. Therefore, the regulation of bacterial genes expression during the colonization of the digestive tract influences the coexistence of phages and bacteria, highlighting the intricacy of tripartite relationships between phages, bacteria, and the animal host in intestinal homeostasis.

A Phage Foundry Framework to Systematically Develop Viral Countermeasures to Combat Antibiotic-Resistant Bacterial Pathogens

At its current rate, the rise of antimicrobial-resistant (AMR) infections is predicted to paralyze our industries and healthcare facilities while becoming the leading global cause of loss of human life. With limited new antibiotics on the horizon, we need to invest in alternative solutions. Bacteriophages (phages)—viruses targeting bacteria—offer a powerful alternative approach to tackle bacterial infections. Despite recent advances in using phages to treat recalcitrant AMR infections, the field lacks systematic development of phage therapies scalable to different applications. We propose a Phage Foundry framework to establish metrics for phage characterization and to fill the knowledge and technological gaps in phage therapeutics. Coordinated investment in AMR surveillance, sampling, characterization, and data sharing procedures will enable rational exploitation of phages for treatments. A fully realized Phage Foundry will enhance the sharing of knowledge, technology, and viral reagents in an equitable manner and will accelerate the biobased economy.

Insights into Gene Transcriptional Regulation of Kayvirus Bacteriophages Obtained from Therapeutic Mixtures

Bacteriophages (phages) of the genus Kayvirus of Staphylococcus aureus are promising agents for therapeutic applications. In this study, we isolated Kayvirus phages, SAM1 and SAM2, from the Fersisi commercial phage cocktail (George Eliava Institute, Tbilisi, Georgia), which exhibits high sequence homology with phage K (≥94%, BLASTn). We found that phages SAM1 and SAM2 infected 95% and 86% of 21 MRSA of differing sequence types (MLST, SCCmec type) obtained from the Irish National MRSA collection, respectively. We conducted differential transcriptomic analysis by RNA-Seq on phage SAM1 during host infection, showing differential expression of its genes at different points during host infection. This analysis also allowed the identification of potentially adverse outcomes in the application of these phages to target MRSA as therapy. The interaction of phage SAM1 on the host caused the upregulation of prophage genes. Additionally, phage infection was found to cause the slight upregulation of host genes implicated in virulence factors relating to hemolysins, immune evasion, and adhesion, but also the downregulation of genes associated with enterotoxins. The findings of this study give further insights into the biology of kayviruses and their use as therapeutics.

Nutrient driven transcriptional changes during phage infection in an aquatic Gammaproteobacterium

Phages modulate bacterial metabolism during infection by regulating gene expression, which influences aquatic nutrient cycling. However, the effects of shifting nutrient regimes are less understood. Here, we analyzed transcriptomes of an ecologically relevant Gammaproteobacterium and its lytic phage in high (HNM) and low (LNM) nutrient medium. Despite different infection characteristics, including reduced burst size and longer latent period in LNM, the phage had a fixed expression profile. Bacterial transcription was instead different depending on nutrient regime, with HNM bacteria focusing on growth while LNM bacteria focused on motility and membrane transport. Additionally, phage infection had a larger effect on bacterial gene expression in LNM compared to HNM, e.g. suppressing increased iron uptake and altering expression of phosphorus uptake genes. Overall, phage infection influenced host metabolism more in LNM, which was more similar to natural conditions, emphasizing the importance of considering natural conditions to understand phage and host ecology.

Comparative metabolomics revealed key pathways associated with the synergistic killing of multidrug-resistant Klebsiella pneumoniae by a bacteriophage-polymyxin combination

Resistance to the last-line polymyxins is emerging in multidrug-resistant Klebsiella pneumoniae and phage therapy is a promising alternative. However, phage monotherapy often rapidly causes resistance and few studies have examined antibiotic-phage combinations against K. pneumoniae. Here, we investigated the combination of polymyxin B with a novel phage pK8 against an mcr-1-carrying polymyxin-resistant clinical isolate Kp II-503 (polymyxin B MIC, 8 mg/L). The phage genome was sequenced and bacterial metabolomes were analysed at 4 and 24 h following the treatment with polymyxin B (16 mg/L), phage pK8 (10² PFU/mL) and their combination. Minimal metabolic changes across 24 h were observed with polymyxin B alone; whereas a significant inhibition of the citrate cycle, pentose phosphate pathway, amino acid and nucleotide metabolism occurred with the phage-polymyxin combination at both 4 and 24 h, but with phage alone only at 4 h. The development of resistance to phage alone was associated with enhanced membrane lipid and decreased amino acid biosynthesis in Kp II-503. Notably, cAMP, cGMP and cCMP were significantly enriched (3.1–6.6 log₂fold) by phage alone and the combination only at 4 h. This is the first systems pharmacology study to investigate the enhanced bacterial killing by polymyxin-phage combination and provides important mechanistic information on phage killing, resistance and antibiotic-phage combination in K. pneumoniae.

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.

Clostridioides difficile - phage relationship the RNA way

Clostridioides difficile (formerly Clostridium difficile)-associated diarrhea is currently the most frequently occurring nosocomial diarrhea worldwide. During its infection cycle this pathogen needs to survive in phage-rich gut communities. Recent data strongly suggest that regulatory RNAs control gene expression in C. difficile and many of these RNAs appear to modulate C. difficile-phage interactions. Of the 200 regulatory RNAs identified by deep sequencing and targeted approaches, many function as antitoxins within type I toxin-antitoxin modules and CRISPR RNAs for anti-phage defenses. In this review, we discuss recent insights into the role of RNAs in modulating interactions between C. difficile and phages in light of intriguing data in other prokaryotes.

Diversity, taxonomy, and evolution of archaeal viruses of the class Caudoviricetes

The archaeal tailed viruses (arTV), evolutionarily related to tailed double-stranded DNA (dsDNA) bacteriophages of the class Caudoviricetes, represent the most common isolates infecting halophilic archaea. Only a handful of these viruses have been genomically characterized, limiting our appreciation of their ecological impacts and evolution. Here, we present 37 new genomes of haloarchaeal tailed virus isolates, more than doubling the current number of sequenced arTVs. Analysis of all 63 available complete genomes of arTVs, which we propose to classify into 14 new families and 3 orders, suggests ancient divergence of archaeal and bacterial tailed viruses and points to an extensive sharing of genes involved in DNA metabolism and counterdefense mechanisms, illuminating common strategies of virus-host interactions with tailed bacteriophages. Coupling of the comparative genomics with the host range analysis on a broad panel of haloarchaeal species uncovered 4 distinct groups of viral tail fiber adhesins controlling the host range expansion. The survey of metagenomes using viral hallmark genes suggests that the global architecture of the arTV community is shaped through recurrent transfers between different biomes, including hypersaline, marine, and anoxic environments.

Differential transcription profiling of the phage LUZ19 infection process in different growth media

RNA sequencing of phage-infected bacterial cultures offers a snapshot of transcriptional events occurring during the infection process, providing insights into the phage transcriptional organization as well as the bacterial response. To better mimic real environmental contexts, we performed RNA-seq of Pseudomonas aeruginosa PAO1 cultures infected with phage LUZ19 in a mammalian cell culture medium to better simulate a phage therapy event and the data were compared to lysogeny broth medium. Regardless of the media, phage LUZ19 induces significant transcriptional changes in the bacterial host over time, particularly during early infection (t = 5 min) and gradually shuts down bacterial transcription. In a common response in both media, 56 P. aeruginosa PAO1 genes are differentially transcribed and clustered into several functional categories such as metabolism, translation and transcription. Our data allowed us to tease apart a medium-specific response during infection from the identified infection-associated responses. This reinforces the concept that phages overtake bacterial transcriptome in a strict manner to gain control of the bacterial machinery and reallocate resources for infection, in this case overcoming the nutritional limitations of the mammalian cell culture medium. From a phage therapy perspective, this study contributes towards a better understanding of phage-host interaction in human physiological conditions and demonstrates the versatility of phage LUZ19 to adapt to different environments.

Molecular Characteristics of Novel Phage vB_ShiP-A7 Infecting Multidrug-Resistant Shigella flexneri and Escherichia coli, and Its Bactericidal Effect in vitro and in vivo

In recent years, increasing evidence has shown that bacteriophages (phages) can inhibit infection caused by multidrug-resistant (MDR) bacteria. Here, we isolated a new phage, named vB_ShiP-A7, using MDR Shigella flexneri as the host. vB_ShiP-A7 is a novel member of Podoviridae, with a latency period of approximately 35 min and a burst size of approximately 100 phage particles/cell. The adsorption rate constant of phage vB_ShiP-A7 to its host S. flexneri was 1.405 × 10^–8 mL/min. The vB_ShiP-A7 genome is a linear double-stranded DNA composed of 40,058 bp with 177 bp terminal repeats, encoding 43 putative open reading frames. Comparative genomic analysis demonstrated that the genome sequence of vB_ShiP-A7 is closely related to 15 different phages, which can infect different strains. Mass spectrometry analysis revealed that 12 known proteins and 6 hypothetical proteins exist in the particles of phage vB_ShiP-A7. Our results confirmed that the genome of vB_ShiP-A7 is free of lysogen-related genes, bacterial virulence genes, and antibiotic resistance genes. vB_ShiP-A7 can significantly disrupt the growth of some MDR clinical strains of S. flexneri and Escherichia coli in liquid culture and biofilms in vitro. In addition, vB_ShiP-A7 can reduce the load of S. flexneri by approximately 3–10 folds in an infection model of mice. Therefore, vB_ShiP-A7 is a stable novel phage with the potential to treat infections caused by MDR strains of S. flexneri and E. coli.

Hostile Takeover: How Viruses Reprogram Prokaryotic Metabolism

To reproduce, prokaryotic viruses must hijack the cellular machinery of their hosts and redirect it toward the production of viral particles. While takeover of the host replication and protein synthesis apparatus has long been considered an essential feature of infection, recent studies indicate that extensive reprogramming of host primary metabolism is a widespread phenomenon among prokaryotic viruses that is required to fulfill the biosynthetic needs of virion production. In this review we provide an overview of the most significant recent findings regarding virus-induced reprogramming of prokaryotic metabolism and suggest how quantitative systems biology approaches may be used to provide a holistic understanding of metabolic remodeling during lytic viral infection. Expected final online publication date for the Annual Review of Microbiology, Volume 75 is October 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Evolutionary biology and development model of medicines: A necessary 'pas de deux' for future successful bacteriophage therapy

The increase in frequency of multidrug-resistant bacteria worldwide is largely the result of the massive use of antibiotics in the second half of the 20th century. These relatively recent changes in human societies revealed the great evolutionary capacities of bacteria towards drug resistance. In this article, we hypothesize that the success of future antibacterial strategies lies in taking into account both these evolutionary processes and the way human activities influence them. Faced with the increasing prevalence of multidrug-resistant bacteria and the scarcity of new antibacterial chemical molecules, the use of bacteriophages is considered as a complementary and/or alternative therapy. After presenting the evolutionary capacities of bacteriophages and bacteria, we show how the development model currently envisaged (based on the classification of bacteriophages as medicinal products similar to antibacterial chemical molecules) ignores the evolutionary processes inherent in bacteriophage therapy. This categorization imposes to bacteriophage therapy a specific conception of what a treatment and a therapeutic scheme should be as well as its mode of production and prescription. We argue that a new development model is needed that would allow the use of therapeutic bacteriophages fully adapted (after in vitro 'bacteriophage training') to the aetiologic bacteria and/or aimed at rendering bacteria either avirulent or antibiotic-susceptible ('bacteriophage steering'). To not repeat the mistakes made with antibiotics, we must now think about and learn from the ways in which the materialities of microbes (e.g. evolutionary capacities of both bacteriophages and bacteria) are intertwined with those of societies.

A phage protein-derived antipathogenic peptide that targets type IV pilus assembly

Phage-inspired antibacterial discovery is a new approach that recruits phages in search for antibacterials with new molecular targets, in that phages are the biological entities well adapted to hijack host bacterial physiology in favor of their own thrive. We previously observed that phage-mediated twitching motility inhibition was effective to control the acute infections caused by Pseudomonas aeruginosa and that the motility inhibition was attributed to the delocalization of PilB, the type IV pilus (TFP) assembly ATPase by binding of the 136-amino acid (aa) phage protein, Tip. Here, we created a series of truncated and point-mutant Tip proteins to identify the critical residues in the Tip bioactivity: N-terminal 80-aa residues were dispensable for the Tip activity; we identified that Asp82, Leu84, and Arg85 are crucial in the Tip function. Furthermore, a synthetic 15-aa peptide (P1) that corresponds to Leu73 to Ala87 is shown to suffice for PilB delocalization, twitching inhibition, and virulence attenuation upon exogenous administration. The transgenic flies expressing the 15-aa peptide were resistant to P. aeruginosa infections as well. Taken together, this proof-of-concept study reveals a new antipathogenic peptide hit targeting bacterial motility and provides an insight into antibacterial discovery targeting TFP assembly.

Individuals at risk for rheumatoid arthritis harbor differential intestinal bacteriophage communities with distinct metabolic potential

Rheumatoid arthritis (RA) is an autoimmune disease characterized in seropositive individuals by the presence of anti-cyclic citrullinated protein (CCP) antibodies. RA is linked to the intestinal microbiota, yet the association of microbes with CCP serology and their contribution to RA is unclear. We describe intestinal phage communities of individuals at risk for developing RA, with or without anti-CCP antibodies, whose first-degree relatives have been diagnosed with RA. We show that at-risk individuals harbor intestinal phage compositions that diverge based on CCP serology, are dominated by Streptococcaceae, Bacteroidaceae, and Lachnospiraceae phages, and may originate from disparate ecosystems. These phages encode unique repertoires of auxiliary metabolic genes, which associate with anti-CCP status, suggesting that these phages directly influence the metabolic and immunomodulatory capability of the microbiota. This work sets the stage for the use of phages as preclinical biomarkers and provides insight into a possible microbial-based causation of RA disease development.

Proteomic and Transcriptomic Analysis of Microviridae φX174 Infection Reveals Broad Upregulation of Host Escherichia coli Membrane Damage and Heat Shock Responses

A major part of the healthy human gut microbiome is the Microviridae bacteriophage, exemplified by the model φX174 phage, and their E. coli hosts. Although much has been learned from studying φX174 over the last half-century, until this work, the E. coli host response to infection has never been investigated in detail.

Measuring host-bacteriophage dynamics is an important approach to understanding bacterial survival functions and responses to infection. The model Microviridae bacteriophage φX174 is endemic to the human gut and has been studied for over 70 years, but the host response to infection has never been investigated in detail. To address this gap in our understanding of this important interaction within our microbiome, we have measured host Escherichia coli C proteomic and transcriptomic response to φX174 infection. We used mass spectrometry and RNA sequencing (RNA-seq) to identify and quantify all 11 φX174 proteins and over 1,700 E. coli proteins, enabling us to comprehensively map host pathways involved in φX174 infection. Most notably, we see significant host responses centered on membrane damage and remodeling, cellular chaperone and translocon activity, and lipoprotein processing, which we speculate is due to the peptidoglycan-disruptive effects of the φX174 lysis protein E on MraY activity. We also observe the massive upregulation of small heat shock proteins IbpA/B, along with other heat shock pathway chaperones, and speculate on how the specific characteristics of holdase protein activity may be beneficial for viral infections. Together, this study enables us to begin to understand the proteomic and transcriptomic host responses of E. coli to Microviridae infections and contributes insights to the activities of this important model host-phage interaction.

IMPORTANCE A major part of the healthy human gut microbiome is the Microviridae bacteriophage, exemplified by the model φX174 phage, and their E. coli hosts. Although much has been learned from studying φX174 over the last half-century, until this work, the E. coli host response to infection has never been investigated in detail. We reveal the proteomic and transcriptomic pathways differentially regulated during the φX174 infection cycle and uncover the details of a coordinated cellular response to membrane damage that results in increased lipoprotein processing and membrane trafficking, likely due to the phage antibiotic-like lysis protein. We also reveal that small heat shock proteins IbpA/B are massively upregulated during infection and that these holdase chaperones are highly conserved across the domains of life, indicating that reliance on them is likely widespread across viruses.

The Development of Bacteriophage Resistance in Vibrio alginolyticus Depends on a Complex Metabolic Adaptation Strategy

Lytic bacteriophages have been well documented to play a pivotal role in microbial ecology due to their complex interactions with bacterial species, especially in aquatic habitats. Although the use of phages as antimicrobial agents, known as phage therapy, in the aquatic environment has been increasing, recent research has revealed drawbacks due to the development of phage-resistant strains among Gram-negative species. Acquired phage resistance in marine Vibrios has been proven to be a very complicated process utilizing biochemical, metabolic, and molecular adaptation strategies. The results of our multi-omics approach, incorporating transcriptome and metabolome analyses of Vibrio alginolyticus phage-resistant strains, corroborate this prospect. Our results provide insights into phage-tolerant strains diminishing the expression of phage receptors ompF, lamB, and btuB. The same pattern was observed for genes encoding natural nutrient channels, such as rbsA, ptsG, tryP, livH, lysE, and hisp, meaning that the cell needs to readjust its biochemistry to achieve phage resistance. The results showed reprogramming of bacterial metabolism by transcript regulations in key-metabolic pathways, such as the tricarboxylic acid cycle (TCA) and lysine biosynthesis, as well as the content of intracellular metabolites belonging to processes that could also significantly affect the cell physiology. Finally, SNP analysis in resistant strains revealed no evidence of amino acid alterations in the studied putative bacterial phage receptors, but several SNPs were detected in genes involved in transcriptional regulation. This phenomenon appears to be a phage-specific, fine-tuned metabolic engineering, imposed by the different phage genera the bacteria have interacted with, updating the role of lytic phages in microbial marine ecology.

Characterization of a novel roseophage and the morphological and transcriptional response of the sponge symbiont Ruegeria AU67 to infection

Sponges have recently been recognized to contain complex communities of bacteriophages; however, little is known about how they interact with their bacterial hosts. Here, we isolated a novel phage, called Ruegeria phage Tedan, and characterized its impact on the bacterial sponge symbiont Ruegeria AU67 on a morphological and molecular level. Phage Tedan was structurally, genomically and phylogenetically characterized to be affiliated with the genus Xiamenvirus of the family Siphoviridae. Through microscopic observations and transcriptomic analysis, we show that phage Tedan upon infection induces a process leading to metabolic and morphological changes in its host. These changes would render Ruegeria AU67 better adapted to inhabit the sponge holobiont due to an improved utilization of ecologically relevant energy and carbon sources as well as a potential impediment of phagocytosis by the sponge through cellular enlargement. An increased survival or better growth of the bacterium in the sponge environment will likely benefit the phage reproduction. Our results point towards the possibility that phages from host-associated environments require, and have thus evolved, different strategies to interact with their host when compared to those phages from free-living or planktonic environments.

Bacterial Vivisection: How Fluorescence-Based Imaging Techniques Shed a Light on the Inner Workings of Bacteria

The rise in fluorescence-based imaging techniques over the past 3 decades has improved the ability of researchers to scrutinize live cell biology at increased spatial and temporal resolution. In microbiology, these real-time vivisections structurally changed the view on the bacterial cell away from the "watery bag of enzymes" paradigm toward the perspective that these organisms are as complex as their eukaryotic counterparts. Capitalizing on the enormous potential of (time-lapse) fluorescence microscopy and the ever-extending pallet of corresponding probes, initial breakthroughs were made in unraveling the localization of proteins and monitoring real-time gene expression. However, later it became clear that the potential of this technique extends much further, paving the way for a focus-shift from observing single events within bacterial cells or populations to obtaining a more global picture at the intra- and intercellular level. In this review, we outline the current state of the art in fluorescence-based vivisection of bacteria and provide an overview of important case studies to exemplify how to use or combine different strategies to gain detailed information on the cell's physiology. The manuscript therefore consists of two separate (but interconnected) parts that can be read and consulted individually. The first part focuses on the fluorescent probe pallet and provides a perspective on modern methodologies for microscopy using these tools. The second section of the review takes the reader on a tour through the bacterial cell from cytoplasm to outer shell, describing strategies and methods to highlight architectural features and overall dynamics within cells.

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.

High-throughput mapping of the phage resistance landscape in E. coli

Bacteriophages (phages) are critical players in the dynamics and function of microbial communities and drive processes as diverse as global biogeochemical cycles and human health. Phages tend to be predators finely tuned to attack specific hosts, even down to the strain level, which in turn defend themselves using an array of mechanisms. However, to date, efforts to rapidly and comprehensively identify bacterial host factors important in phage infection and resistance have yet to be fully realized. Here, we globally map the host genetic determinants involved in resistance to 14 phylogenetically diverse double-stranded DNA phages using two model Escherichia coli strains (K-12 and BL21) with known sequence divergence to demonstrate strain-specific differences. Using genome-wide loss-of-function and gain-of-function genetic technologies, we are able to confirm previously described phage receptors as well as uncover a number of previously unknown host factors that confer resistance to one or more of these phages. We uncover differences in resistance factors that strongly align with the susceptibility of K-12 and BL21 to specific phage. We also identify both phage-specific mechanisms, such as the unexpected role of cyclic-di-GMP in host sensitivity to phage N4, and more generic defenses, such as the overproduction of colanic acid capsular polysaccharide that defends against a wide array of phages. Our results indicate that host responses to phages can occur via diverse cellular mechanisms. Our systematic and high-throughput genetic workflow to characterize phage-host interaction determinants can be extended to diverse bacteria to generate datasets that allow predictive models of how phage-mediated selection will shape bacterial phenotype and evolution. The results of this study and future efforts to map the phage resistance landscape will lead to new insights into the coevolution of hosts and their phage, which can ultimately be used to design better phage therapeutic treatments and tools for precision microbiome engineering.

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.

Transcriptomic determinants of the response of ST-111 Pseudomonas aeruginosa AG1 to ciprofloxacin identified by a top-down systems biology approach

Pseudomonas aeruginosa is an opportunistic pathogen that thrives in diverse environments and causes a variety of human infections. Pseudomonas aeruginosa AG1 (PaeAG1) is a high-risk sequence type 111 (ST-111) strain isolated from a Costa Rican hospital in 2010. PaeAG1 has both blaVIM-2 and blaIMP-18 genes encoding for metallo-β-lactamases, and it is resistant to β-lactams (including carbapenems), aminoglycosides, and fluoroquinolones. Ciprofloxacin (CIP) is an antibiotic commonly used to treat P. aeruginosa infections, and it is known to produce DNA damage, triggering a complex molecular response. In order to evaluate the effects of a sub-inhibitory CIP concentration on PaeAG1, growth curves using increasing CIP concentrations were compared. We then measured gene expression using RNA-Seq at three time points (0, 2.5 and 5 h) after CIP exposure to identify the transcriptomic determinants of the response (i.e. hub genes, gene clusters and enriched pathways). Changes in expression were determined using differential expression analysis and network analysis using a top–down systems biology approach. A hybrid model using database-based and co-expression analysis approaches was implemented to predict gene–gene interactions. We observed a reduction of the growth curve rate as the sub-inhibitory CIP concentrations were increased. In the transcriptomic analysis, we detected that over time CIP treatment resulted in the differential expression of 518 genes, showing a complex impact at the molecular level. The transcriptomic determinants were 14 hub genes, multiple gene clusters at different levels (associated to hub genes or as co-expression modules) and 15 enriched pathways. Down-regulation of genes implicated in several metabolism pathways, virulence elements and ribosomal activity was observed. In contrast, amino acid catabolism, RpoS factor, proteases, and phenazines genes were up-regulated. Remarkably, > 80 resident-phage genes were up-regulated after CIP treatment, which was validated at phenomic level using a phage plaque assay. Thus, reduction of the growth curve rate and increasing phage induction was evidenced as the CIP concentrations were increased. In summary, transcriptomic and network analyses, as well as the growth curves and phage plaque assays provide evidence that PaeAG1 presents a complex, concentration-dependent response to sub-inhibitory CIP exposure, showing pleiotropic effects at the systems level. Manipulation of these determinants, such as phage genes, could be used to gain more insights about the regulation of responses in PaeAG1 as well as the identification of possible therapeutic targets. To our knowledge, this is the first report of the transcriptomic analysis of CIP response in a ST-111 high-risk P. aeruginosa strain, in particular using a top-down systems biology approach.

Bacteriophages as sources of small non-coding RNA molecules

Bacteriophages play an essential role in the transferring of genes that contribute to the bacterial virulence and whose products are dangerous to human health. Interestingly, phages carrying virulence genes are mostly temperate and in contrast to lytic phages undergo both lysogenic and lytic cycles. Importantly, expression of the majority of phage genes and subsequent production of phage encoded proteins is suppressed during lysogeny. The expression of the majority of phage genes is tightly linked to lytic development. Among others, small non-coding RNAs (sRNAs) of phage origin are involved in the regulation of phage gene expression and thus play an important role in both phage and host development. In the case of bacteria, sRNAs affect processes such as virulence, colonization ability, motility and cell growth or death. In turn, in the case of phages, they play essential roles during the early stage of infection, maintaining the state of lysogeny and silencing the expression of late structural genes, thereby regulating the transition between phage life cycles. Interestingly, sRNAs have been identified in both lytic and temperate phages and they have been discussed in this work according to this classification. Particular attention was paid to viral sRNAs resembling eukaryotic microRNAs.

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

Pseudomonas virus vB_PaeM_PA5oct is proposed as a model jumbo bacteriophage to investigate phage-bacteria interactions and is a candidate for phage therapy applications. Combining hybrid sequencing, RNA-Seq and mass spectrometry allowed us to accurately annotate its 286,783 bp genome with 461 coding regions including four non-coding RNAs (ncRNAs) and 93 virion-associated proteins. PA5oct relies on the host RNA polymerase for the infection cycle and RNA-Seq revealed a gradual take-over of the total cell transcriptome from 21% in early infection to 93% in late infection. PA5oct is not organized into strictly contiguous regions of temporal transcription, but some genomic regions transcribed in early, middle and late phases of infection can be discriminated. Interestingly, we observe regions showing limited transcription activity throughout the infection cycle. We show that PA5oct upregulates specific bacterial operons during infection including operons pncA-pncB1-nadE involved in NAD biosynthesis, psl for exopolysaccharide biosynthesis and nap for periplasmic nitrate reductase production. We also observe a downregulation of T4P gene products suggesting mechanisms of superinfection exclusion. We used the proteome of PA5oct to position our isolate amongst other phages using a gene-sharing network. This integrative omics study illustrates the molecular diversity of jumbo viruses and raises new questions towards cellular regulation and phage-encoded hijacking mechanisms.

Isolation and Characterization of the Novel Phage JD032 and Global Transcriptomic Response during JD032 Infection of Clostridioides difficile Ribotype 078

C. difficile is one of the most clinically significant intestinal pathogens. Although phages have been shown to effectively control C. difficile infection, the host responses to phage predation have not been fully studied. In this study, we reported the isolation and characterization of a new phage, JD032, and analyzed the global transcriptomic changes in the hypervirulent RT078 C. difficile strain, TW11, during phage JD032 infection. We found that bacterial host mRNA was progressively replaced with phage transcripts, three temporal categories of JD032 gene expression, the extensive interplay between phage-bacterium, antiphage-like responses of the host and phage evasion, and decreased expression of sporulation- and virulence-related genes of the host after phage infection. These findings confirmed the complexity of interactions between C. difficile and phages and suggest that phages undergoing a lytic cycle may also cause different phenotypes in hosts, similar to prophages, which may inspire phage therapy for the control of C. difficile.

Insights into the interaction between phages and their bacterial hosts are crucial for the development of phage therapy. However, only one study has investigated global gene expression of Clostridioides (formerly Clostridium) difficile carrying prophage, and transcriptional reprogramming during lytic infection has not been studied. Here, we presented the isolation, propagation, and characterization of a newly discovered 35,109-bp phage, JD032, and investigated the global transcriptomes of both JD032 and C. difficile ribotype 078 (RT078) strain TW11 during JD032 infection. Transcriptome sequencing (RNA-seq) revealed the progressive replacement of bacterial host mRNA with phage transcripts. The expressed genes of JD032 were clustered into early, middle, and late temporal categories that were functionally similar. Specifically, a gene (JD032_orf016) involved in the lysis-lysogeny decision was identified as an early expression gene. Only 17.7% (668/3,781) of the host genes were differentially expressed, and more genes were downregulated than upregulated. The expression of genes involved in host macromolecular synthesis (DNA/RNA/proteins) was altered by JD032 at the level of transcription. In particular, the expression of the ropA operon was downregulated. Most noteworthy is that the gene expression of some antiphage systems, including CRISPR-Cas, restriction-modification, and toxin-antitoxin systems, was suppressed by JD032 during infection. In addition, bacterial sporulation, adhesion, and virulence factor genes were significantly downregulated. This study provides the first description of the interaction between anaerobic spore-forming bacteria and phages during lytic infection and highlights new aspects of C. difficile phage-host interactions.

IMPORTANCEC. difficile is one of the most clinically significant intestinal pathogens. Although phages have been shown to effectively control C. difficile infection, the host responses to phage predation have not been fully studied. In this study, we reported the isolation and characterization of a new phage, JD032, and analyzed the global transcriptomic changes in the hypervirulent RT078 C. difficile strain, TW11, during phage JD032 infection. We found that bacterial host mRNA was progressively replaced with phage transcripts, three temporal categories of JD032 gene expression, the extensive interplay between phage-bacterium, antiphage-like responses of the host and phage evasion, and decreased expression of sporulation- and virulence-related genes of the host after phage infection. These findings confirmed the complexity of interactions between C. difficile and phages and suggest that phages undergoing a lytic cycle may also cause different phenotypes in hosts, similar to prophages, which may inspire phage therapy for the control of C. difficile.

Prevalence of small base-pairing RNAs derived from diverse genomic loci

Small RNAs (sRNAs) that act by base-pairing have been shown to play important roles in fine-tuning the levels and translation of their target transcripts across a variety of model and pathogenic organisms. Work from many different groups in a wide range of bacterial species has provided evidence for the importance and complexity of sRNA regulatory networks, which allow bacteria to quickly respond to changes in their environment. However, despite the expansive literature, much remains to be learned about all aspects of sRNA-mediated regulation, particularly in bacteria beyond the well-characterized Escherichia coli and Salmonella enterica species. Here we discuss what is known, and what remains to be learned, about the identification of regulatory base-pairing RNAs produced from diverse genomic loci including how their expression is regulated. This article is part of a Special Issue entitled: RNA and gene control in bacteria edited by Dr. M. Guillier and F. Repoila.

Parallel Genomics Uncover Novel Enterococcal-Bacteriophage Interactions

We lack fundamental understanding of how phage infection influences bacterial gene expression and, consequently, how bacterial responses to phage infection affect the assembly of polymicrobial communities. Using parallel genomic approaches, we have discovered novel transcriptional regulators and metabolic genes that influence phage infection. The integration of whole-genome transcriptomic profiling during phage infection has revealed the differential regulation of genes important for group behaviors and polymicrobial interactions. Our work suggests that therapeutic phages could more broadly influence bacterial community composition outside their intended host targets.

Bacteriophages (phages) have been proposed as alternative therapeutics for the treatment of multidrug-resistant bacterial infections. However, there are major gaps in our understanding of the molecular events in bacterial cells that control how bacteria respond to phage predation. Using the model organism Enterococcus faecalis, we used two distinct genomic approaches, namely, transposon library screening and RNA sequencing, to investigate the interaction of E. faecalis with a virulent phage. We discovered that a transcription factor encoding a LytR family response regulator controls the expression of enterococcal polysaccharide antigen (epa) genes that are involved in phage infection and bacterial fitness. In addition, we discovered that DNA mismatch repair mutants rapidly evolve phage adsorption deficiencies, underpinning a molecular basis for epa mutation during phage infection. Transcriptomic profiling of phage-infected E. faecalis revealed broad transcriptional changes influencing viral replication and progeny burst size. We also demonstrate that phage infection alters the expression of bacterial genes associated with intra- and interbacterial interactions, including genes involved in quorum sensing and polymicrobial competition. Together, our results suggest that phage predation has the potential to influence complex microbial behavior and may dictate how bacteria respond to external environmental stimuli. These responses could have collateral effects (positive or negative) on microbial communities, such as the host microbiota, during phage therapy.