Sequence characteristics of T4-like bacteriophage IME08 benome termini revealed by high throughput sequencing

The Beauty of Bacteriophage T4 Research: Lindsay W. Black and the T4 Head Assembly

Viruses are biochemically complex structures and mainly consist of folded proteins that contain nucleic acids. Bacteriophage T4 is one of most prominent examples, having a tail structure that contracts during the infection process. Intracellular phage multiplication leads to separate self-directed assembly reactions of proheads, tails and tail fibers. The proheads are packaged with concatemeric DNA produced by tandem replication reactions of the parental DNA molecule. Once DNA packaging is completed, the head is joined with the tail and six long fibers are attached. The mature particles are then released from the cell via lysis, another tightly regulated process. These processes have been studied in molecular detail leading to a fascinating view of the protein-folding dynamics that direct the structural interplay of assembled complexes. Lindsay W. Black dedicated his career to identifying and defining the molecular events required to form the T4 virion. He leaves us with rich insights into the astonishingly precise molecular clockwork that co-ordinates all of the players in T4 assembly, both viral and cellular. Here, we summarize Lindsay’s key research contributions that are certain to stimulate our future science for many years to come.

Analyzing Genome Termini of Bacteriophage Through High-Throughput Sequencing

High-throughput sequencing (HTS) is an effective tool for bacteriophage genome and its termini analysis. HTS technology parallelizes the sequencing process, producing thousands to millions of reads concurrently. Terminal information of a bacteriophage genome is important and basic knowledge for understanding the biology of the bacteriophage. We have created a high-occurrence reads as termini theory and developed practical methods to determine the bacteriophage genome termini, which is based on the large data of HTS. With this method, the termini of the bacteriophage genome can be efficiently and reliably identified as a by-product of bacteriophage genome sequencing, by solely analyzing the sequence statistics of the raw sequencing data (reads), without any further lab experiments.

Predicting genome terminus sequences of Bacillus cereus-group bacteriophage using next generation sequencing data

Background

Most tailed bacteriophages (phages) feature linear dsDNA genomes. Characterizing novel phages requires an understanding of complete genome sequences, including the definition of genome physical ends.

Result

We sequenced 48 Bacillus cereus phage isolates and analyzed Next-generation sequencing (NGS) data to resolve the genome configuration of these novel phages. Most assembled contigs featured reads that mapped to both contig ends and formed circularized contigs. Independent assemblies of 31 nearly identical I48-like Bacillus phage isolates allowed us to observe that the assembly programs tended to produce random cleavage on circularized contigs. However, currently available assemblers were not capable of reporting the underlying phage genome configuration from sequence data. To identify the genome configuration of sequenced phage in silico, a terminus prediction method was developed by means of ‘neighboring coverage ratios’ and ‘read edge frequencies’ from read alignment files. Termini were confirmed by primer walking and supported by phylogenetic inference of large DNA terminase protein sequences.

Conclusions

The Terminus package using phage NGS data along with the contig circularity could efficiently identify the proximal positions of phage genome terminus. Complete phage genome sequences allow a proposed characterization of the potential packaging mechanisms and more precise genome annotation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3744-0) contains supplementary material, which is available to authorized users.

Characterization of Pseudomonas aeruginosa Phage C11 and Identification of Host Genes Required for Virion Maturation

The underlying mechanisms of phage-host interactions largely remained to be elucidated. In this work, Pseudomonas aeruginosa phage C11 was first characterized as a Myoviridae virus having a linear dsDNA molecule of 94109 bp with 1173 bp identical terminal direct repeats (TDR). Then the mutants resistant to phage C11 were screened in a Tn5G transposon mutant library of P. aeruginosa PAK, including two mutants with decreased adsorption rates (DAR) and five mutants with wild-type adsorption rates (WAR). When the WAR mutants were incubated with phage C11, their growth rates were significantly inhibited; the replication of the phage genomic DNA was detected in all the WAR mutants with the real-time quantitative PCR analysis; and the synthesized phage genomic DNA was processed into monomers for packaging evidenced by the southern blot analysis. Moreover, with strain PAK as indicator, small quantities of phage C11 were synthesized in the WAR mutants. Taken together, these data suggested the identified genes of the WAR mutants are necessary for efficient synthesis of the infectious phage particles. Finally, the WAR mutants were detected sensitive to two other Pseudomonas phages closely related with C11, further implying the evolved diversity and complexity of the phage-host interactions in both sides.

The search for therapeutic bacteriophages uncovers one new subfamily and two new genera of Pseudomonas-infecting Myoviridae

In a previous study, six virulent bacteriophages PAK_P1, PAK_P2, PAK_P3, PAK_P4, PAK_P5 and CHA_P1 were evaluated for their in vivo efficacy in treating Pseudomonas aeruginosa infections using a mouse model of lung infection. Here, we show that their genomes are closely related to five other Pseudomonas phages and allow a subdivision into two clades, PAK_P1-like and KPP10-like viruses, based on differences in genome size, %GC and genomic contents, as well as number of tRNAs. These two clades are well delineated, with a mean of 86% and 92% of proteins considered homologous within individual clades, and 25% proteins considered homologous between the two clades. By ESI-MS/MS analysis we determined that their virions are composed of at least 25 different proteins and electron microscopy revealed a morphology identical to the hallmark Salmonella phage Felix O1. A search for additional bacteriophage homologs, using profiles of protein families defined from the analysis of the 11 genomes, identified 10 additional candidates infecting hosts from different species. By carrying out a phylogenetic analysis using these 21 genomes we were able to define a new subfamily of viruses, the Felixounavirinae within the Myoviridae family. The new Felixounavirinae subfamily includes three genera: Felixounalikevirus, PAK_P1likevirus and KPP10likevirus. Sequencing genomes of bacteriophages with therapeutic potential increases the quantity of genomic data on closely related bacteriophages, leading to establishment of new taxonomic clades and the development of strategies for analyzing viral genomes as presented in this article.

Characterization, sequencing and comparative genomic analysis of vB_AbaM-IME-AB2, a novel lytic bacteriophage that infects multidrug-resistant Acinetobacter baumannii clinical isolates

BACKGROUND

With the use of broad-spectrum antibiotics, immunosuppressive drugs, and glucocorticoids, multidrug-resistant Acinetobacter baumannii (MDR-AB) has become a major nosocomial pathogen species. The recent renaissance of bacteriophage therapy may provide new treatment strategies for combatting drug-resistant bacterial infections. In this study, we isolated a lytic bacteriophage vB_AbaM-IME-AB2 has a short latent period and a small burst size, which clear its host's suspension quickly, was selected for characterization and a complete genomic comparative study.

RESULTS

The isolated bacteriophage vB_AbaM-IME-AB2 has an icosahedral head and displays morphology resembling Myoviridae family. Gel separation assays showed that the phage particle contains at least nine protein bands with molecular weights ranging 15-100 kDa. vB_AbaM-IME-AB2 could adsorb its host cells in 9 min with an adsorption rate more than 99% and showed a short latent period (20 min) and a small burst size (62 pfu/cell). It could form clear plaques in the double-layer assay and clear its host's suspension in just 4 hours. Whole genome of vB_AbaM-IME-AB2 was sequenced and annotated and the results showed that its genome is a double-stranded DNA molecule consisting of 43,665 nucleotides. The genome has a G + C content of 37.5% and 82 putative coding sequences (CDSs). We compared the characteristics and complete genome sequence of all known Acinetobacter baumannii bacteriophages. There are only three that have been sequenced Acinetobacter baumannii phages AB1, AP22, and phiAC-1, which have a relatively high similarity and own a coverage of 65%, 50%, 8% respectively when compared with our phage vB_AbaM-IME-AB2. A nucleotide alignment of the four Acinetobacter baumannii phages showed that some CDSs are similar, with no significant rearrangements observed. Yet some sections of these strains of phage are nonhomologous.

CONCLUSION

vB_AbaM-IME-AB2 was a novel and unique A. baumannii bacteriophage. These findings suggest a common ancestry and microbial diversity and evolution. A clear understanding of its characteristics and genes is conducive to the treatment of multidrug-resistant A. baumannii in the future.

A novel Pseudomonas aeruginosa bacteriophage, Ab31, a chimera formed from temperate phage PAJU2 and P. putida lytic phage AF: characteristics and mechanism of bacterial resistance

A novel temperate bacteriophage of Pseudomonas aeruginosa, phage vB_PaeP_Tr60_Ab31 (alias Ab31) is described. Its genome is composed of structural genes related to those of lytic P. putida phage AF, and regulatory genes similar to those of temperate phage PAJU2. The virion structure resembles that of phage AF and other lytic Podoviridae (S. enterica Epsilon 15 and E. coli phiv10) with similar tail spikes. Ab31 was able to infect P. aeruginosa strain PA14 and two genetically related strains called Tr60 and Tr162, out of 35 diverse strains from cystic fibrosis patients. Analysis of resistant host variants revealed different phenotypes, including induction of pigment and alginate overproduction. Whole genome sequencing of resistant variants highlighted the existence of a large deletion of 234 kbp in two strains, encompassing a cluster of genes required for the production of CupA fimbriae. Stable lysogens formed by Ab31 in strain Tr60, permitted the identification of the insertion site. During colonization of the lung in cystic fibrosis patients, P. aeruginosa adapts by modifying its genome. We suggest that bacteriophages such as Ab31 may play an important role in this adaptation by selecting for bacterial characteristics that favor persistence of bacteria in the lung.

Scrutinizing virus genome termini by high-throughput sequencing

Analysis of genomic terminal sequences has been a major step in studies on viral DNA replication and packaging mechanisms. However, traditional methods to study genome termini are challenging due to the time-consuming protocols and their inefficiency where critical details are lost easily. Recent advances in next generation sequencing (NGS) have enabled it to be a powerful tool to study genome termini. In this study, using NGS we sequenced one iridovirus genome and twenty phage genomes and confirmed for the first time that the high frequency sequences (HFSs) found in the NGS reads are indeed the terminal sequences of viral genomes. Further, we established a criterion to distinguish the type of termini and the viral packaging mode. We also obtained additional terminal details such as terminal repeats, multi-termini, asymmetric termini. With this approach, we were able to simultaneously detect details of the genome termini as well as obtain the complete sequence of bacteriophage genomes. Theoretically, this application can be further extended to analyze larger and more complicated genomes of plant and animal viruses. This study proposed a novel and efficient method for research on viral replication, packaging, terminase activity, transcription regulation, and metabolism of the host cell.

Characterization of the Salmonella bacteriophage vB_SenS-Ent1

The bacteriophage vB_SenS-Ent1 (Ent1) is a member of the family Siphoviridae of tailed bacteriophages and infects a broad range of serovars of the enteric pathogen Salmonella enterica. The virion particle is composed of an icosahedral head 64 nm in diameter and a flexible, non-contractile tail of 116 × 8.5 nm possessing terminal fibres. The adsorption rate constant at 37 °C is 6.73 × 10(-9) ml min(-1). Latent and eclipse periods are 25 and 20 min, respectively, and the burst size is 35 progeny particles per cell after 35 min at 37 °C. Sequencing revealed a circularly permuted, 42 391 bp dsDNA genome containing 58 ORFs organized into four major transcriptional units. Comparisons with the genome sequences of other bacteriophages revealed a high level of nucleotide sequence identity and shared orthologous proteins with the Salmonella phages SETP3, SE2 and KS7 (SS3e) and the Escherichia phages K1G, K1H, K1ind1 and K1ind3.

Structure, assembly, and DNA packaging of the bacteriophage T4 head

The bacteriophage T4 head is an elongated icosahedron packed with 172 kb of linear double-stranded DNA and numerous proteins. The capsid is built from three essential proteins: gp23*, which forms the hexagonal capsid lattice; gp24*, which forms pentamers at 11 of the 12 vertices; and gp20, which forms the unique dodecameric portal vertex through which DNA enters during packaging and exits during infection. Intensive work over more than half a century has led to a deep understanding of the phage T4 head. The atomic structure of gp24 has been determined. A structural model built for gp23 using its similarity to gp24 showed that the phage T4 major capsid protein has the same fold as numerous other icosahedral bacteriophages. However, phage T4 displays an unusual membrane and portal initiated assembly of a shape determining self-sufficient scaffolding core. Folding of gp23 requires the assistance of two chaperones, the Escherichia coli chaperone GroEL acting with the phage-coded gp23-specific cochaperone, gp31. The capsid also contains two nonessential outer capsid proteins, Hoc and Soc, which decorate the capsid surface. Through binding to adjacent gp23 subunits, Soc reinforces the capsid structure. Hoc and Soc have been used extensively in bipartite peptide display libraries and to display pathogen antigens, including those from human immunodeficiency virus (HIV), Neisseria meningitides, Bacillus anthracis, and foot and mouth disease virus. The structure of Ip1*, one of a number of multiple (>100) copy proteins packed and injected with DNA from the full head, shows it to be an inhibitor of one specific restriction endonuclease specifically targeting glycosylated hydroxymethyl cytosine DNA. Extensive mutagenesis, combined with atomic structures of the DNA packaging/terminase proteins gp16 and gp17, elucidated the ATPase and nuclease functional motifs involved in DNA translocation and headful DNA cutting. The cryoelectron microscopy structure of the T4 packaging machine showed a pentameric motor assembled with gp17 subunits on the portal vertex. Single molecule optical tweezers and fluorescence studies showed that the T4 motor packages DNA at the highest rate known and can package multiple segments. Förster resonance energy transfer-fluorescence correlation spectroscopy studies indicate that DNA gets compressed in the stalled motor and that the terminase-to-portal distance changes during translocation. Current evidence suggests a linear two-component (large terminase plus portal) translocation motor in which electrostatic forces generated by ATP hydrolysis drive DNA translocation by alternating the motor between tensed and relaxed states.