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

A phage nucleus-associated RNA-binding protein is required for jumbo phage infection

Large-genome bacteriophages (jumbo phages) of the proposed family Chimalliviridae assemble a nucleus-like compartment bounded by a protein shell that protects the replicating phage genome from host-encoded restriction enzymes and DNA-targeting CRISPR-Cas nucleases. While the nuclear shell provides broad protection against host nucleases, it necessitates transport of mRNA out of the nucleus-like compartment for translation by host ribosomes, and transport of specific proteins into the nucleus-like compartment to support DNA replication and mRNA transcription. Here, we identify a conserved phage nuclear shell-associated protein that we term Chimallin C (ChmC), which adopts a nucleic acid-binding fold, binds RNA with high affinity in vitro, and binds phage mRNAs in infected cells. ChmC also forms phase-separated condensates with RNA in vitro. Targeted knockdown of ChmC using mRNA-targeting dCas13d results in accumulation of phage-encoded mRNAs in the phage nucleus, reduces phage protein production, and compromises virion assembly. Taken together, our data show that the conserved ChmC protein plays crucial roles in the viral life cycle, potentially by facilitating phage mRNA translocation through the nuclear shell to promote protein production and virion development.

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.

The enigmatic epitranscriptome of bacteriophages: putative RNA modifications in viral infections

RNA modifications play essential roles in modulating RNA function, stability, and fate across all kingdoms of life. The entirety of the RNA modifications within a cell is defined as the epitranscriptome. While eukaryotic RNA modifications are intensively studied, understanding bacterial RNA modifications remains limited, and knowledge about bacteriophage RNA modifications is almost nonexistent. In this review, we shed light on known mechanisms of bacterial RNA modifications and propose how this knowledge might be extended to bacteriophages. We build hypotheses on enzymes potentially responsible for regulating the epitranscriptome of bacteriophages and their host. This review highlights the exciting prospects of uncovering the unexplored field of bacteriophage epitranscriptomics and its potential role to shape bacteriophage-host interactions.

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.

Genome-Scale Analysis of the Structure and Function of RNA Pathways and Networks in Pseudomonas aeruginosa

In recent years, several genome-wide approaches based on RNA sequencing (RNA-seq) have been developed. These methods allow a comprehensive and dynamic view of the structure and function of the multi-layered RNA pathways and networks. Many of these approaches, including the promising one of single-cell transcriptome analysis, have been successfully applied to Pseudomonas aeruginosa. However, we are only at the beginning because only a few surrounding conditions have been considered. Here, we aim to illustrate the different types of approaches based on RNA-seq that will lead us in the future to a better understanding of the dynamics of RNA biology in P. aeruginosa.

A phage nucleus-associated RNA-binding protein is required for jumbo phage infection

Large-genome bacteriophages (jumbo phages) of the Chimalliviridae family assemble a nucleus-like compartment bounded by a protein shell that protects the replicating phage genome from host-encoded restriction enzymes and CRISPR/Cas nucleases. While the nuclear shell provides broad protection against host nucleases, it necessitates transport of mRNA out of the nucleus-like compartment for translation by host ribosomes, and transport of specific proteins into the nucleus-like compartment to support DNA replication and mRNA transcription. Here we identify a conserved phage nuclear shell-associated protein that we term Chimallin C (ChmC), which adopts a nucleic acid-binding fold, binds RNA with high affinity in vitro, and binds phage mRNAs in infected cells. ChmC also forms phase-separated condensates with RNA in vitro. Targeted knockdown of ChmC using mRNA-targeting dCas13d halts infections at an early stage. Taken together, our data suggest that the conserved ChmC protein acts as a chaperone for phage mRNAs, potentially stabilizing these mRNAs and driving their translocation through the nuclear shell to promote translation and infection progression.

Complete genome analysis of Tequatrovirus ufvareg1, a Tequatrovirus species inhibiting Escherichia coli O157:H7

Introduction

Bacteriophages infecting human pathogens have been considered potential biocontrol agents, and studying their genetic content is essential to their safe use in the food industry. Tequatrovirus ufvareg1 is a bacteriophage named UFV-AREG1, isolated from cowshed wastewater and previously tested for its ability to inhibit Escherichia coli O157:H7.

Methods

T. ufvareg1 was previously isolated using E. coli O157:H7 (ATCC 43895) as a bacterial host. The same strain was used for bacteriophage propagation and the one-step growth curve. The genome of the T. ufvareg1 was sequenced using 305 Illumina HiSeq, and the genome comparison was calculated by VIRIDIC and VIPTree.

Results

Here, we characterize its genome and compare it to other Tequatrovirus. T. ufvareg1 virions have an icosahedral head (114 x 86 nm) and a contracted tail (117 x 23 nm), with a latent period of 25 min, and an average burst size was 18 phage particles per infected E. coli cell. The genome of the bacteriophage T. ufvareg1 contains 268 coding DNA sequences (CDS) and ten tRNA genes distributed in both negative and positive strains. T. ufvareg1 genome also contains 40 promoters on its regulatory regions and two rho-independent terminators. T. ufvareg1 shares an average intergenomic similarity (VIRIDC) of 88.77% and an average genomic similarity score (VipTree) of 88.91% with eight four reference genomes for Tequatrovirus available in the NCBI RefSeq database. The pan-genomic analysis confirmed the high conservation of Tequatrovirus genomes. Among all CDS annotated in the T. ufvareg1 genome, there are 123 core genes, 38 softcore genes, 94 shell genes, and 13 cloud genes. None of 268 CDS was classified as being exclusive of T. ufvareg1.

Conclusion

The results in this paper, combined with other previously published findings, indicate that T. ufvareg1 bacteriophage is a potential candidate for food protection against E. coli O157:H7 in foods.

Pf4 Phage Variant Infection Reduces Virulence-Associated Traits in Pseudomonas aeruginosa

Pf4 is a filamentous bacteriophage integrated as a prophage into the genome of Pseudomonas aeruginosa PAO1. Pf4 virions can be produced without killing P. aeruginosa. However, cell lysis can occur during superinfection when Pf virions successfully infect a host lysogenized by a Pf superinfective variant. We have previously shown that infection of P. aeruginosa PAO1 with a superinfective Pf4 variant abolished twitching motility and altered biofilm architecture. More precisely, most of the cells embedded into the biofilm were showing a filamentous morphology, suggesting the activation of the cell envelope stress response involving both AlgU and SigX extracytoplasmic function sigma factors. Here, we show that Pf4 variant infection results in a drastic dysregulation of 3,360 genes representing about 58% of P. aeruginosa genome; of these, 70% of the virulence factors encoding genes show a dysregulation. Accordingly, Pf4 variant infection (termed Pf4*) causes in vivo reduction of P. aeruginosa virulence and decreased production of N-acyl-homoserine lactones and 2-alkyl-4-quinolones quorum-sensing molecules and related virulence factors, such as pyocyanin, elastase, and pyoverdine. In addition, the expression of genes involved in metabolism, including energy generation and iron homeostasis, was affected, suggesting further relationships between virulence and central metabolism. Altogether, these data show that Pf4 phage variant infection results in complex network dysregulation, leading to reducing acute virulence in P. aeruginosa. This study contributes to the comprehension of the bacterial response to filamentous phage infection. IMPORTANCE Filamentous bacteriophages can become superinfective and infect P. aeruginosa, even though they are inserted in the genome as lysogens. Despite this productive infection, growth of the host is only mildly affected, allowing the study of the interaction between the phage and the host, which is not possible in the case of lytic phages killing rapidly their host. Here, we demonstrate by transcriptome and phenotypic analysis that the infection by a superinfective filamentous phage variant causes a massive disruption in gene expression, including those coding for virulence factors and metabolic pathways.

SAPPHIRE.CNN: Implementation of dRNA-seq-driven, species-specific promoter prediction using convolutional neural networks

Data availability is a consistent bottleneck for the development of bacterial species-specific promoter prediction software. In this work we leverage genome-wide promoter datasets generated with dRNA-seq in the Gram-negative bacteria Pseudomonas aeruginosa and Salmonella enterica for promoter prediction. Convolutional neural networks are presented as an optimal architecture for model training and are further modified and tailored for promoter prediction. The resulting predictors reach high binary accuracies (95% and 94.9%) on test sets and outperform each other when predicting promoters in their associated species. SAPPHIRE.CNN is available online and can also be downloaded to run locally. Our results indicate a dependency of binary promoter classification on an organism’s GC content and a decreased performance of our classifiers on genera they were not trained for, further supporting the need for dedicated, species-specific promoter classification tools.

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.

Temporal Transcriptional Responses of a Vibrio alginolyticus Strain to Podoviridae Phage HH109 Revealed by RNA-Seq

Vibrio alginolyticus is a common opportunistic pathogen that causes mass mortality in cultured marine animals. Phage HH109 lyses pathogenic V. alginolyticus strain E110 with high efficiency and thus serves as a useful model to understand the dynamic interplay of a phage and its host.

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.

Genomic Diversity of Bacteriophages Infecting the Genus Acinetobacter

The number of sequenced Acinetobacter phage genomes in the International Nucleotide Sequence Database Collaboration has increased significantly in recent years, from 37 in 2017 to a total of 139 as of January 2021 with genome sizes ranging from 31 to 378 kb. Here, we explored the genetic diversity of the Acinetobacter phages using comparative genomics approaches that included assessment of nucleotide similarity, shared gene content, single gene phylogeny, and the network-based classification tool vConTACT2. Phages infecting Acinetobacter sp. are genetically diverse and can be grouped into 8 clusters (subfamilies) and 46 sub-clusters (genera), of which 8 represent genomic singletons (additional genera). We propose the creation of five new subfamilies and suggest a reorganisation of the genus Obolenskvirus. These results provide an updated view of the viruses infecting Acinetobacter species, providing insights into their diversity.

Shaping the Bacterial Epitranscriptome-5'-Terminal and Internal RNA Modifications

All domains of life utilize a diverse set of modified ribonucleotides that can impact the sequence, structure, function, stability, and the fate of RNAs, as well as their interactions with other molecules. Today, more than 160 different RNA modifications are known that decorate the RNA at the 5'-terminus or internal RNA positions. The boost of next-generation sequencing technologies sets the foundation to identify and study the functional role of RNA modifications. The recent advances in the field of RNA modifications reveal a novel regulatory layer between RNA modifications and proteins, which is central to developing a novel concept called "epitranscriptomics." The majority of RNA modifications studies focus on the eukaryotic epitranscriptome. In contrast, RNA modifications in prokaryotes are poorly characterized. This review outlines the current knowledge of the prokaryotic epitranscriptome focusing on mRNA modifications. Here, it is described that several internal and 5'-terminal RNA modifications either present or likely present in prokaryotic mRNA. Thereby, the individual techniques to identify these epitranscriptomic modifications, their writers, readers and erasers, and their proposed functions are explored. Besides that, still unanswered questions in the field of prokaryotic epitranscriptomics are pointed out, and its future perspectives in the dawn of next-generation sequencing technologies are outlined.

A Grad-seq View of RNA and Protein Complexes in Pseudomonas aeruginosa under Standard and Bacteriophage Predation Conditions

The Gram-negative rod-shaped bacterium Pseudomonas aeruginosa is not only a major cause of nosocomial infections but also serves as a model species of bacterial RNA biology. While its transcriptome architecture and posttranscriptional regulation through the RNA-binding proteins Hfq, RsmA, and RsmN have been studied in detail, global information about stable RNA-protein complexes in this human pathogen is currently lacking. Here, we implement gradient profiling by sequencing (Grad-seq) in exponentially growing P. aeruginosa cells to comprehensively predict RNA and protein complexes, based on glycerol gradient sedimentation profiles of >73% of all transcripts and ∼40% of all proteins. As to benchmarking, our global profiles readily reported complexes of stable RNAs of P. aeruginosa, including 6S RNA with RNA polymerase and associated product RNAs (pRNAs). We observe specific clusters of noncoding RNAs, which correlate with Hfq and RsmA/N, and provide a first hint that P. aeruginosa expresses a ProQ-like FinO domain-containing RNA-binding protein. To understand how biological stress may perturb cellular RNA/protein complexes, we performed Grad-seq after infection by the bacteriophage ΦKZ. This model phage, which has a well-defined transcription profile during host takeover, displayed efficient translational utilization of phage mRNAs and tRNAs, as evident from their increased cosedimentation with ribosomal subunits. Additionally, Grad-seq experimentally determines previously overlooked phage-encoded noncoding RNAs. Taken together, the Pseudomonas protein and RNA complex data provided here will pave the way to a better understanding of RNA-protein interactions during viral predation of the bacterial cell.