Genome Modularization Reveals Overlapped Gene Topology Is Necessary for Efficient Viral Reproduction

A computational model for bacteriophage ϕX174 gene expression

Bacteriophage ϕX174 has been widely used as a model organism to study fundamental processes in molecular biology. However, several aspects of ϕX174 gene regulation are not fully resolved. Here we construct a computational model for ϕX174 and use the model to study gene regulation during the phage infection cycle. We estimate the relative strengths of transcription regulatory elements (promoters and terminators) by fitting the model to transcriptomics data. We show that the specific arrangement of a promoter followed immediately by a terminator, which occurs naturally in the ϕX174 genome, poses a parameter identifiability problem for the model, since the activity of one element can be partially compensated for by the other. We also simulate ϕX174 gene expression with two additional, putative transcription regulatory elements that have been proposed in prior studies. We find that the activities of these putative elements are estimated to be weak, and that variation in ϕX174 transcript abundances can be adequately explained without them. Overall, our work demonstrates that ϕX174 gene regulation is well described by the canonical set of promoters and terminators widely used in the literature.

An Improved Method for Eliminating or Creating Intragenic Bacterial Promoters

Recent advances in genomic refactoring have been hindered by the ever-present complication of internal or cryptic transcriptional regulation. Typical approaches to these features have been to randomize or perform mass alterations to the gene sequences thought to contain the regulatory motifs; however, this approach can cause problems by altering translational speeds, introducing long distance DNA-DNA interaction effects, and inducing RNA toxicity. Previously, we developed a rational design approach named COdon Restrained Promoter SilEncing (CORPSE) which takes externally identified promoter sequences and uses position-specific scoring matrices as proxy promoter strengths to make minimal changes to promoter sequences to disable their activity. Additionally, through inverting our system we were also able to modify weak internal promoters to increase their activity. In this chapter, we augment our previous process with the biophysical model Promoter Calculator v1.0 developed by LaFleur et al. to combine promoter identification and activity prediction, with our algorithm to silently modify promoter sequences, to provide more robust promoter elimination and creation.

Characterization and Genome Study of a Newly Isolated Temperate Phage Belonging to a New Genus Targeting Alicyclobacillus acidoterrestris

The spoilage of juices by Alicyclobacillus spp. remains a serious problem in industry and leads to economic losses. Compounds such as guaiacol and halophenols, which are produced by Alicyclobacillus, create undesirable flavors and odors and, thus, decrease the quality of juices. The inactivation of Alicyclobacillus spp. constitutes a challenge because it is resistant to environmental factors, such as high temperatures, and active acidity. However, the use of bacteriophages seems to be a promising approach. In this study, we aimed to isolate and comprehensively characterize a novel bacteriophage targeting Alicyclobacillus spp. The Alicyclobacillus phage strain KKP 3916 was isolated from orchard soil against the Alicyclobacillus acidoterrestris strain KKP 3133. The bacterial host's range and the effect of phage addition at different rates of multiplicity of infections (MOIs) on the host's growth kinetics were determined using a Bioscreen C Pro growth analyzer. The Alicyclobacillus phage strain KKP 3916, retained its activity in a wide range of temperatures (from 4 °C to 30 °C) and active acidity values (pH from 3 to 11). At 70 °C, the activity of the phage decreased by 99.9%. In turn, at 80 °C, no activity against the bacterial host was observed. Thirty minutes of exposure to UV reduced the activity of the phages by almost 99.99%. Based on transmission-electron microscopy (TEM) and whole-genome sequencing (WGS) analyses, the Alicyclobacillus phage strain KKP 3916 was classified as a tailed bacteriophage. The genomic sequencing revealed that the newly isolated phage had linear double-stranded DNA (dsDNA) with sizes of 120 bp and 131 bp and 40.3% G+C content. Of the 204 predicted proteins, 134 were of unknown function, while the remainder were annotated as structural, replication, and lysis proteins. No genes associated with antibiotic resistance were found in the genome of the newly isolated phage. However, several regions, including four associated with integration into the bacterial host genome and excisionase, were identified, which indicates the temperate (lysogenic) life cycle of the bacteriophage. Due to the risk of its potential involvement in horizontal gene transfer, this phage is not an appropriate candidate for further research on its use in food biocontrol. To the best of our knowledge, this is the first article on the isolation and whole-genome analysis of the Alicyclobacillus-specific phage.

Transcriptional analysis in bacteriophage Fc02 of Pseudomonas aeruginosa revealed two overlapping genes with exclusion activity

Little is known about the gene expression program during the transition from lysogenic to lytic cycles of temperate bacteriophages in Pseudomonas aeruginosa. To investigate this issue, we developed a thermo-sensitive repressor mutant in a lysogen and analyzed the phage transcriptional program by strand-specific RNA-Seq before and after thermo-induction. As expected, the repressor gene located on the phage DNA forward strand is transcribed in the lysogen at the permissive temperature of 30°C. Upstream the repressor gene, we noticed the presence of two overlapped ORFs apparently in the same transcript. One ORF is a gene that encodes a protein of 7.9 kDa mediating the exclusion of various super-infecting phages. The other ORF, placed in an alternate reading frame with a possible AUG initiation codon at 25 nucleotide downstream of the AUG of the first gene, is expected to encode a 20.7 kDa polypeptide of yet an unknown function. Upon lifting repression at 40°C, the transcription of an operon which is involved in the lytic cycle is started from a promoter on the reverse phage DNA strand. The first gene in the operon is a homolog of the antirepresor ner, a common gene in the lysis-lysogeny regulation region of other phages. Interestingly, the next gene after ner is gene 10 that on the reverse strand overlaps the overlapped gene olg1 on the forward strand. Curiously, gene 10 expression also shows superinfection exclusion. Strand-specific RNA-Seq also has uncovered the transcription succession of gene modules expressed during the phage lytic stage. The conservation of overlapped genes with similar functions may be evolutionarily selected.

Creating De Novo Overlapped Genes

Future applications of synthetic biology will rely on deploying engineered cells outside of lab environments for long periods of time. Currently, a significant roadblock to this application is the potential for deactivating mutations in engineered genes. A recently developed method to protect engineered coding sequences from mutation is called Constraining Adaptive Mutations using Engineered Overlapping Sequences (CAMEOS). In this chapter we provide a workflow for utilizing CAMEOS to create synthetic overlaps between two genes, one essential (infA) and one non-essential (aroB), to protect the non-essential gene from mutation and loss of protein function. In this workflow we detail the methods to collect large numbers of related protein sequences, produce multiple sequence alignments (MSAs), use the MSAs to generate hidden Markov models and Markov random field models, and finally generate a library of overlapping coding sequences through CAMEOS scripts. To assist practitioners with basic coding skills to try out the CAMEOS method, we have created a virtual machine containing all the required packages already installed that can be downloaded and run locally.

Computational Tools for the Analysis of Uncultivated Phage Genomes

Over a century of bacteriophage research has uncovered a plethora of fundamental aspects of their biology, ecology, and evolution. Furthermore, the introduction of community-level studies through metagenomics has revealed unprecedented insights on the impact that phages have on a range of ecological and physiological processes. It was not until the introduction of viral metagenomics that we began to grasp the astonishing breadth of genetic diversity encompassed by phage genomes. Novel phage genomes have been reported from a diverse range of biomes at an increasing rate, which has prompted the development of computational tools that support the multilevel characterization of these novel phages based solely on their genome sequences. The impact of these technologies has been so large that, together with MAGs (Metagenomic Assembled Genomes), we now have UViGs (Uncultivated Viral Genomes), which are now officially recognized by the International Committee for the Taxonomy of Viruses (ICTV), and new taxonomic groups can now be created based exclusively on genomic sequence information. Even though the available tools have immensely contributed to our knowledge of phage diversity and ecology, the ongoing surge in software programs makes it challenging to keep up with them and the purpose each one is designed for. Therefore, in this review, we describe a comprehensive set of currently available computational tools designed for the characterization of phage genome sequences, focusing on five specific analyses: (i) assembly and identification of phage and prophage sequences, (ii) phage genome annotation, (iii) phage taxonomic classification, (iv) phage-host interaction analysis, and (v) phage microdiversity.

Gene Overlapping as a Modulator of Begomovirus Evolution

In RNA viruses, which have high mutation—and fast evolutionary— rates, gene overlapping (i.e., genomic regions that encode more than one protein) is a major factor controlling mutational load and therefore the virus evolvability. Although DNA viruses use host high-fidelity polymerases for their replication, and therefore should have lower mutation rates, it has been shown that some of them have evolutionary rates comparable to those of RNA viruses. Notably, these viruses have large proportions of their genes with at least one overlapping instance. Hence, gene overlapping could be a modulator of virus evolution beyond the RNA world. To test this hypothesis, we use the genus Begomovirus of plant viruses as a model. Through comparative genomic approaches, we show that terminal gene overlapping decreases the rate of virus evolution, which is associated with lower frequency of both synonymous and nonsynonymous mutations. In contrast, terminal overlapping has little effect on the pace of virus evolution. Overall, our analyses support a role for gene overlapping in the evolution of begomoviruses and provide novel information on the factors that shape their genetic diversity.

Overlapping genes in natural and engineered genomes

Modern genome-scale methods that identify new genes, such as proteogenomics and ribosome profiling, have revealed, to the surprise of many, that overlap in genes, open reading frames and even coding sequences is widespread and functionally integrated into prokaryotic, eukaryotic and viral genomes. In parallel, the constraints that overlapping regions place on genome sequences and their evolution can be harnessed in bioengineering to build more robust synthetic strains and constructs. With a focus on overlapping protein-coding and RNA-coding genes, this Review examines their discovery, topology and biogenesis in the context of their genome biology. We highlight exciting new uses for sequence overlap to control translation, compress synthetic genetic constructs, and protect against mutation.

Plaque Size Tool: An automated plaque analysis tool for simplifying and standardising bacteriophage plaque morphology measurements

Bacteriophage plaque size measurement is essential for phage characterisation, but manual size estimation requires a considerable amount of time and effort. In order to ease the work of phage researchers, we have developed an automated command-line application called Plaque Size Tool (PST) that can detect plaques of different morphology on the images of Petri dishes and measure plaque area and diameter. Plaque size measurements using PST showed no difference to those obtained with manual plaque size measurement in Fiji, indicating future results using PST are backwards compatible with prior measurements in the literature. PST can be applied to a range of lytic bacteriophages producing oval-shaped plaques, including bull's-eye and turbid morphology. The application can also be used for titer calculation if most of the plaques are stand-alone. As laboratory automation becomes more commonplace, standardised and flexible open-source analytical tools like PST will be important parts of biofoundry and cloud lab bacteriophage workflows.

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.