Inferring bacteriophage infection strategies from genome sequence: analysis of bacteriophage 7-11 and related phages

Isolation and characterization of novel plasmid-dependent phages infecting bacteria carrying diverse conjugative plasmids

IMPORTANCE

This work was undertaken because plasmid-dependent phages can reduce the prevalence of conjugative plasmids and can be leveraged to prevent the acquisition and dissemination of ARGs by bacteria. The two novel phages described in this study, Lu221 and Hi226, can infect Escherichia coli, Salmonella enterica, Kluyvera sp. and Enterobacter sp. carrying conjugative plasmids. This was verified with plasmids carrying resistance determinants and belonging to the most common plasmid families among Gram-negative pathogens. Therefore, the newly isolated phages could have the potential to help control the spread of ARGs and thus help combat the antimicrobial resistance crisis.

Genomic characterization of bacteriophage pSal-SNUABM-01, a novel elongated-head phage infecting Salmonella sp

Salmonellosis is a disease of critical concern for public health, and the use of bacteriophages is among the most promising approaches to combating Salmonella. As Salmonella has various serotypes and strains, and bacteriophages are virulent to specific hosts, it is important to isolate phages and evaluate interactions with their hosts. In the present study, a novel Salmonella-infecting bacteriophage, pSal-SNUABM-01, was isolated and characterized. Transmission electron microscopy revealed that the bacteriophage is a member of the family Podoviridae and possesses an elongated head and a short tail. The phage genome is circular and 89,500 bp in size. A total of 162 open reading frames were predicted, eight of which were tRNAs. Morphological and genomic analysis revealed that pSal-SNUABM-01 is closely related to phage 7-11. In phylogenetic analysis, pSal-SNUABM-01 and 7-11 did not cluster together with the members of any established genus, suggesting that these two phages comprise a novel genus. The results of this study enhance our understanding of the phylogeny of the family Podoviridae and might be applicable to the development of bacteriophage treatments against Salmonella infections.

Infection Kinetics and Phylogenetic Analysis of vB_EcoD_SU57, a Virulent T1-Like Drexlerviridae Coliphage

The morphology, infection kinetics, genome sequence and phylogenetic characterization of the previously isolated bacteriophage vB_EcoD_SU57 are presented. The phage vB_EcoD_SU57 was isolated on Escherichia coli strain ECOR57 from the E. coli reference collection and was shown to produce four mm clear plaques with halos. Infection kinetics, as assessed by one-step growth analyses, suggest that vB_EcoD_SU57 is a virulent phage with an adsorption rate of 8.5 × 10^–10 mL × min^–1, a latency period of 14 min, and a burst size of 13 PFU per bacterium. Transmission electron microscopy confirmed vB_EcoD_SU57 to be a phage that used to be classified as a Siphoviridae phage. Bioinformatics analyses showed that the genome was 46,150 base pairs long, contained 29 genes with predicted protein functions, and 51 open reading frames encoding proteins with unknown function, many of which were gathered in clusters. A putative tRNA gene was also identified. Phylogenetic analyses showed that vB_EcoD_SU57 is a Braunvirinae phage of the newly formed Drexlerviridae family and closely related to T1-like E. coli phages vB_EcoS_ACG-M12 (Guelphvirus) and Rtp (Rtpvirus) as well as the unclassified phages vB_EcoS_CEB_EC3a and ECH1.

From biophysics to 'omics and systems biology

Recent decades brought a revolution to biology, driven mainly by exponentially increasing amounts of data coming from "'omics" sciences. To handle these data, bioinformatics often has to combine biologically heterogeneous signals, for which methods from statistics and engineering (e.g. machine learning) are often used. While such an approach is sometimes necessary, it effectively treats the underlying biological processes as a black box. Similarly, systems biology deals with inherently complex systems, characterized by a large number of degrees of freedom, and interactions that are highly non-linear. To deal with this complexity, the underlying physical interactions are often (over)simplified, such as in Boolean modelling of network dynamics. In this review, we argue for the utility of applying a biophysical approach in bioinformatics and systems biology, including discussion of two examples from our research which address sequence analysis and understanding intracellular gene expression dynamics.

Mix-and-matching as a promoter recognition mechanism by ECF σ factors

Background

Transcription initiation is in bacteria exhibited by different σ factors, most of which fall within σ⁷⁰ family. This family is diverse, ranging from the housekeeping Group I (RpoDs), to Group IV (ECF) σ factors, that transcribe smaller regulons under more stringent conditions. RpoDs employ a kinetic mix-and-match mechanism, where promoter elements complement each other binding strengths in achieving sufficient transcription activity. On the other hand, it is assumed that ECF σs, which are the most distant from the housekeeping σ factors, cannot exhibit mix-and-matching. However, mix-and-matching for ECF σ factors was not quantitatively checked before, and recent results show a much larger flexibility in the promoter recognition by the members of this group.

Results

To this end, we quantitatively investigate mix-and-matching in two canonical ECF σ family members (σ^E and σ^W), for which we use a biophysics based model of transcription initiation. For σ^E, we perform a separate analysis for in-vitro active and in-vitro inactive promoters, which allows us investigating how mix-and-matching depends on the external factors that may control transcription activity in the in-vitro inactive set. We show that the promoter elements of canonical ECF σs significantly complement each other strengths, where such mix-and-matching is in the in-vitro active set even stronger compared to the correlations observed for the housekeeping σs. This complementation however significantly decreases for the in-vitro inactive set, which we propose is due to mix-and-matching with regulatory sequences outside of the canonical promoter elements. In line with this proposition, we show that a conserved spacer element, which appears in the in-vitro inactive promoter set, significantly increases the promoter element complementation. While RpoD promoter elements mix-and-match to achieve sufficient total transcription activity, for σ^E they complement each other to achieve sufficiently strong total binding affinity, which we relate to differences in physiological responses between the two groups of σ factors.

Conclusion

Despite a common notion that smaller σ factor specificity leads to a larger mix-and-matching, we here obtain a larger promoter element complementation for σ^E compared to RpoDs. Finally, to explain this finding, we propose a simple model which relates the size of σ factor regulon with the extent of mix-and-matching, based on an assumption of a selection pressure on promoters that are near the non-specific binding boundary to remain functional.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0865-z) contains supplementary material, which is available to authorized users.

Genomic and Transcriptional Mapping of PaMx41, Archetype of a New Lineage of Bacteriophages Infecting Pseudomonas aeruginosa

Previously, a collection of virulent phages infecting Pseudomonas aeruginosa was isolated from open water reservoirs and residual waters. Here, we described the comparative genomics of a set of five related phages from the collection, the physical structure of the genome, the structural proteomics of the virion, and the transcriptional program of archetypal phage PaMx41. The phage genomes were closely associated with each other and with those of two other P. aeruginosa phages, 119X and PaP2, which were previously filed in the databases. Overall, the genomes were approximately 43 kb, harboring 53 conserved open reading frames (ORFs) and three short ORFs in indel regions and containing 45% GC content. The genome of PaMx41 was further characterized as a linear, terminally redundant DNA molecule. A total of 16 ORFs were associated with putative functions, including nucleic acid metabolism, morphogenesis, and lysis, and eight virion proteins were identified through mass spectrometry. However, the coding sequences without assigned functions represent 70% of the ORFs. The PaMx41 transcription program was organized in early, middle, and late expressed genomic modules, which correlated with regions containing functionally related genes. The high genomic conservation among these distantly isolated phages suggests that these viruses undergo selective pressure to remain unchanged. The 119X lineage represents a unique set of phages that corresponds to a novel phage group. The features recognized in the genomes and the broad host range of clinical strains suggest that these phages are candidates for therapy applications.

IMPORTANCE

Pseudomonas aeruginosa is an opportunistic pathogen that causes stubborn nosocomial infections that are frequently resistant to multiple antibiotics. Bacterial viruses (bacteriophages or phages) represent a natural mechanism for pathogenic bacterial control. Here, a group of virulent phages, previously shown to infect a broad range of clinical P. aeruginosa strains, was characterized at the genomic and molecular levels. These phages belong to a unique and tightly related group. In addition, we conducted a transcriptional study of an archetypal phage of this group to characterize the role of many unknown coding sequences based on expression temporalities. These results contribute to our knowledge of 119X-like phages and, in general, provide information concerning P. aeruginosa podophage diversity and lytic cycles.

Promoter Recognition by Extracytoplasmic Function σ Factors: Analyzing DNA and Protein Interaction Motifs

Extracytoplasmic function (ECF) σ factors are the largest and the most diverse group of alternative σ factors, but their mechanisms of transcription are poorly studied. This subfamily is considered to exhibit a rigid promoter structure and an absence of mixing and matching; both -35 and -10 elements are considered necessary for initiating transcription. This paradigm, however, is based on very limited data, which bias the analysis of diverse ECF σ subgroups. Here we investigate DNA and protein recognition motifs involved in ECF σ factor transcription by a computational analysis of canonical ECF subfamily members, much less studied ECF σ subgroups, and the group outliers, obtained from recently sequenced bacteriophages. The analysis identifies an extended -10 element in promoters for phage ECF σ factors; a comparison with bacterial σ factors points to a putative 6-amino-acid motif just C-terminal of domain σ2, which is responsible for the interaction with the identified extension of the -10 element. Interestingly, a similar protein motif is found C-terminal of domain σ2 in canonical ECF σ factors, at a position where it is expected to interact with a conserved motif further upstream of the -10 element. Moreover, the phiEco32 ECF σ factor lacks a recognizable -35 element and σ4 domain, which we identify in a homologous phage, 7-11, indicating that the extended -10 element can compensate for the lack of -35 element interactions. Overall, the results reveal greater flexibility in promoter recognition by ECF σ factors than previously recognized and raise the possibility that mixing and matching also apply to this group, a notion that remains to be biochemically tested.

IMPORTANCE

ECF σ factors are the most numerous group of alternative σ factors but have been little studied. Their promoter recognition mechanisms are obscured by the large diversity within the ECF σ factor group and the limited similarity with the well-studied housekeeping σ factors. Here we extensively compare bacterial and bacteriophage ECF σ factors and their promoters in order to infer DNA and protein recognition motifs involved in transcription initiation. We predict a more flexible promoter structure than is recognized by the current paradigm, which assumes rigidness, and propose that ECF σ promoter elements may complement (mix and match with) each other's strengths. These results warrant the refocusing of research efforts from the well-studied housekeeping σ factors toward the physiologically highly important, but insufficiently understood, alternative σ factors.

Characterization of Five Novel Brevibacillus Bacteriophages and Genomic Comparison of Brevibacillus Phages

Brevibacillus laterosporus is a spore-forming bacterium that causes a secondary infection in beehives following European Foulbrood disease. To better understand the contributions of Brevibacillus bacteriophages to the evolution of their hosts, five novel phages (Jenst, Osiris, Powder, SecTim467, and Sundance) were isolated and characterized. When compared with the five Brevibacillus phages currently in NCBI, these phages were assigned to clusters based on whole genome and proteome synteny. Powder and Osiris, both myoviruses, were assigned to the previously described Jimmer-like cluster. SecTim467 and Jenst, both siphoviruses, formed a novel phage cluster. Sundance, a siphovirus, was assigned as a singleton phage along with the previously isolated singleton, Emery. In addition to characterizing the basic relationships between these phages, several genomic features were observed. A motif repeated throughout phages Jenst and SecTim467 was frequently upstream of genes predicted to function in DNA replication, nucleotide metabolism, and transcription, suggesting transcriptional co-regulation. In addition, paralogous gene pairs that encode a putative transcriptional regulator were identified in four Brevibacillus phages. These paralogs likely evolved to bind different DNA sequences due to variation at amino acid residues predicted to bind specific nucleotides. Finally, a putative transposable element was identified in SecTim467 and Sundance that carries genes homologous to those found in Brevibacillus chromosomes. Remnants of this transposable element were also identified in phage Jenst. These discoveries provide a greater understanding of the diversity of phages, their behavior, and their evolutionary relationships to one another and to their host. In addition, they provide a foundation with which further Brevibacillus phages can be compared.

Bioinformatics as a first-line approach for understanding bacteriophage transcription

Current approach to understanding bacteriophage transcription strategies during infection includes a combination of experimental and bioinformatics approaches, which is often time and resource consuming. Given the exponentially growing number of sequenced bacteriophage genomes, it becomes sensible asking to what extent one can understand bacteriophage transcription by using bioinformatics methods alone. We here argue that a suitable choice of computational methods may provide a highly efficient first-line approach for underst-anding bacteriophage transcription.