Mycobacteriophages: genes and genomes

Whole genome comparison of a large collection of mycobacteriophages reveals a continuum of phage genetic diversity

The bacteriophage population is large, dynamic, ancient, and genetically diverse. Limited genomic information shows that phage genomes are mosaic, and the genetic architecture of phage populations remains ill-defined. To understand the population structure of phages infecting a single host strain, we isolated, sequenced, and compared 627 phages of Mycobacterium smegmatis. Their genetic diversity is considerable, and there are 28 distinct genomic types (clusters) with related nucleotide sequences. However, amino acid sequence comparisons show pervasive genomic mosaicism, and quantification of inter-cluster and intra-cluster relatedness reveals a continuum of genetic diversity, albeit with uneven representation of different phages. Furthermore, rarefaction analysis shows that the mycobacteriophage population is not closed, and there is a constant influx of genes from other sources. Phage isolation and analysis was performed by a large consortium of academic institutions, illustrating the substantial benefits of a disseminated, structured program involving large numbers of freshman undergraduates in scientific discovery.

DOI:http://dx.doi.org/10.7554/eLife.06416.001

Viruses are unable to replicate independently. To generate copies of itself, a virus must instead invade a target cell and commandeer that cell's replication machinery. Different viruses are able to invade different types of cell, and a group of viruses known as bacteriophages (or phages for short) replicate within bacteria. The enormous number and diversity of phages in the world means that they play an important role in virtually every ecosystem.

Despite their importance, relatively little is known about how different phage populations are related to each other and how they evolved. Many phages contain their genetic information in the form of strands of DNA. Using genetic sequencing to find out where and how different genes are encoded in the DNA can reveal information about how different viruses are related to each other. These relationships are particularly complicated in phages, as they can exchange genes with other viruses and microbes.

Previous studies comparing the genomes—the complete DNA sequence—of reasonably small numbers of phages that infect the Mycobacterium group of bacteria have found that the phages can be sorted into ‘clusters’ based on similarities in their genes and where these are encoded in their DNA. However, the number of phages investigated so far has been too small to conclude how different clusters are related. Are the clusters separate, or do they form a ‘continuum’ with different genes and DNA sequences shared between different clusters?

Here, Pope, Bowman, Russell et al. compare the individual genomes of 627 bacteriophages that infect the bacterial species Mycobacterium smegmatis. This is by far the largest number of phage genomes analyzed from a single host species. The large number of genomes analyzed allowed a much clearer understanding of the complexity and diversity of these phages to be obtained. The isolation, sequencing and analysis of the hundreds of M. smegmatis bacteriophage genomes was performed by an integrated research and education program, called the Science Education Alliance Phage Hunters Advancing Genomics and Evolutionary Science (SEA-PHAGES) program. This enabled thousands of undergraduate students from different institutions to contribute to the phage discovery and sequencing project, and co-author the report. SEA-PHAGES therefore shows that it is possible to successfully incorporate genuine scientific research into an undergraduate course, and that doing so can benefit both the students and researchers involved.

The results show that while the genomes could be categorized into 28 clusters, the genomes are not completely unrelated. Instead, a spread of diversity is seen, as genes and groups of genes are shared between different clusters. Pope, Bowman, Russell et al. further reveal that the phage population is in a constant state of change, and continuously acquires genes from other microorganisms and viruses.

DOI:http://dx.doi.org/10.7554/eLife.06416.002

Genomic and proteomic features of mycobacteriophage SWU1 isolated from China soil

Mycobacteriophage SWU1 is a newly isolated phage from soil sample collected in Sichuan province, China using Mycobacterium smegmatis mc(2)155 as host. Plaque, phage morphology and one-step growth curve were characterized. The complete genomic sequence of phage SWU1 was determined by shotgun sequencing. The ends of SWU1 were determined. Structural proteins of SWU1 were analyzed by NanoLC-ESI-MS/MS. Seven ORFs were identified as structural protein encoded by SWU1 genome. The genetic basis underlying the SWU1 plaque was explored using comparative genomics. Prophages homologous to SWU1 were identified in two pathogens, Segniliparus rugosus ATCC BAA-974 and Mycobacterium rhodesiae JS60. Genus Segniliparus is a member of the order Corynebacteriales. To our knowledge, this is the first report of Mycobacterium prophages in different genera.

Genomics and proteomics of mycobacteriophage patience, an accidental tourist in the Mycobacterium neighborhood

Newly emerging human viruses such as Ebola virus, severe acute respiratory syndrome (SARS) virus, and HIV likely originate within an extant population of viruses in nonhuman hosts and acquire the ability to infect and cause disease in humans. Although several mechanisms preventing viral infection of particular hosts have been described, the mechanisms and constraints on viral host expansion are ill defined. We describe here mycobacteriophage Patience, a newly isolated phage recovered using Mycobacterium smegmatis mc(2)155 as a host. Patience has genomic features distinct from its M. smegmatis host, including a much lower GC content (50.3% versus 67.4%) and an abundance of codons that are rarely used in M. smegmatis. Nonetheless, it propagates well in M. smegmatis, and we demonstrate the use of mass spectrometry to show expression of over 75% of the predicted proteins, to identify new genes, to refine the genome annotation, and to estimate protein abundance. We propose that Patience evolved primarily among lower-GC hosts and that the disparities between its genomic profile and that of M. smegmatis presented only a minimal barrier to host expansion. Rapid adaptions to its new host include recent acquisition of higher-GC genes, expression of out-of-frame proteins within predicted genes, and codon selection among highly expressed genes toward the translational apparatus of its new host.

IMPORTANCE

The mycobacteriophage Patience genome has a notably lower GC content (50.3%) than its Mycobacterium smegmatis host (67.4%) and has markedly different codon usage biases. The viral genome has an abundance of codons that are rare in the host and are decoded by wobble tRNA pairing, although the phage grows well and expression of most of the genes is detected by mass spectrometry. Patience thus has the genomic profile of a virus that evolved primarily in one type of host genetic landscape (moderate-GC bacteria) but has found its way into a distinctly different high-GC environment. Although Patience genes are ill matched to the host expression apparatus, this is of little functional consequence and has not evidently imposed a barrier to migration across the microbial landscape. Interestingly, comparison of expression levels and codon usage profiles reveals evidence of codon selection as the genome evolves and adapts to its new environment.

Integration of genomic and proteomic analyses in the classification of the Siphoviridae family

Using a variety of genomic (BLASTN, ClustalW) and proteomic (Phage Proteomic Tree, CoreGenes) tools we have tackled the taxonomic status of members of the largest bacteriophage family, the Siphoviridae. In all over 400 phages were examined and we were able to propose 39 new genera, comprising 216 phage species, and add 62 species to two previously defined genera (Phic3unalikevirus; L5likevirus) grouping, in total, 390 fully sequenced phage isolates. Many of the remainders are orphans which the Bacterial and Archaeal Viruses Subcommittee of the International Committee on Taxonomy of Viruses (ICTV) chooses not to ascribe genus status at the time being.

Recombineering: A powerful tool for modification of bacteriophage genomes

Recombineering, a recently developed technique for efficient genetic manipulation of bacteria, is facilitated by phage-derived recombination proteins and has the advantage of using DNA substrates with short regions of homology. This system was first developed in E. coli but has since been adapted for use in other bacteria. It is now widely used in a number of different systems for a variety of purposes, and the construction of chromosomal gene knockouts, deletions, insertions, point mutations, as well as in vivo cloning, mutagenesis of bacterial artificial chromosomes and phasmids, and the construction of genomic libraries has been reported. However, these methods also can be effectively applied to the genetic modification of bacteriophage genomes, in both their prophage and lytically growing states. The ever-growing collection of fully sequenced bacteriophages raises more questions than they answer, including the unknown functions of vast numbers of genes with no known homologs and of unknown function. Recombineering of phage genomes is central to addressing these questions, enabling the simple construction of mutants, determination of gene essentiality, and elucidation of gene function. In turn, advances in our understanding of phage genomics should present similar recombineering tools for dissecting a multitude of other genetically naïve bacterial systems.

Factors affecting phage D29 infection: a tool to investigate different growth states of mycobacteria

Bacteriophages D29 and TM4 are able to infect a wide range of mycobacteria, including pathogenic and non-pathogenic species. Successful phage infection of both fast- and slow-growing mycobacteria can be rapidly detected using the phage amplification assay. Using this method, the effect of oxygen limitation during culture of mycobacteria on the success of phage infection was studied. Both D29 and TM4 were able to infect cultures of M. smegmatis and Mycobacterium avium subspecies paratuberculosis (MAP) grown in liquid with aeration. However when cultures were grown under oxygen limiting conditions, only TM4 could productively infect the cells. Cell attachment assays showed that D29 could bind to the cells surface but did not complete the lytic cycle. The ability of D29 to productively infect the cells was rapidly recovered (within 1 day) when the cultures were returned to an aerobic environment and this recovery required de novo RNA synthesis. These results indicated that under oxygen limiting conditions the cells are entering a growth state which inhibits phage D29 replication, and this change in host cell biology which can be detected by using both phage D29 and TM4 in the phage amplification assay.

The N-terminal domain of the repressor of Staphylococcus aureus phage Φ11 possesses an unusual dimerization ability and DNA binding affinity

Bacteriophage Φ11 uses Staphylococcus aureus as its host and, like lambdoid phages, harbors three homologous operators in between its two divergently oriented repressor genes. None of the repressors of Φ11, however, showed binding to all three operators, even at high concentrations. To understand why the DNA binding mechanism of Φ11 repressors does not match that of lambdoid phage repressors, we studied the N-terminal domain of the Φ11 lysogenic repressor, as it harbors a putative helix-turn-helix motif. Our data revealed that the secondary and tertiary structures of the N-terminal domain were different from those of the full-length repressor. Nonetheless, the N-terminal domain was able to dimerize and bind to the operators similar to the intact repressor. In addition, the operator base specificity, binding stoichiometry, and binding mechanism of this domain were nearly identical to those of the whole repressor. The binding affinities of the repressor and its N-terminal domain were reduced to a similar extent when the temperature was increased to 42°C. Both proteins also adequately dislodged a RNA polymerase from a Φ11 DNA fragment carrying two operators and a promoter. Unlike the intact repressor, the binding of the N-terminal domain to two adjacent operator sites was not cooperative in nature. Taken together, we suggest that the dimerization and DNA binding abilities of the N-terminal domain of the Φ11 repressor are distinct from those of the DNA binding domains of other phage repressors.

Enigmatic distribution, evolution, and function of inteins

Inteins are mobile genetic elements capable of self-splicing post-translationally. They exist in all three domains of life including in viruses and bacteriophage, where they have a sporadic distribution even among very closely related species. In this review, we address this anomalous distribution from the point of view of the evolution of the host species as well as the intrinsic features of the inteins that contribute to their genetic mobility. We also discuss the incidence of inteins in functionally important sites of their host proteins. Finally, we describe instances of conditional protein splicing. These latter observations lead us to the hypothesis that some inteins have adapted to become sensors that play regulatory roles within their host protein, to the advantage of the organism in which they reside.

Molecular Genetics of Mycobacteriophages

Mycobacteriophages have provided numerous essential tools for mycobacterial genetics, including delivery systems for transposons, reporter genes, and allelic exchange substrates, and components for plasmid vectors and mutagenesis. Their genetically diverse genomes also reveal insights into the broader nature of the phage population and the evolutionary mechanisms that give rise to it. The substantial advances in our understanding of the biology of mycobacteriophages including a large collection of completely sequenced genomes indicates a rich potential for further contributions in tuberculosis genetics and beyond.

Prophage-mediated dynamics of 'Candidatus Liberibacter asiaticus' populations, the destructive bacterial pathogens of citrus huanglongbing

Prophages are highly dynamic components in the bacterial genome and play an important role in intraspecies variations. There are at least two prophages in the chromosomes of Candidatus Liberibacter asiaticus' (Las) Floridian isolates. Las is both unculturable and the most prevalent species of Liberibacter pathogens that cause huanglongbing (HLB), a worldwide destructive disease of citrus. In this study, seven new prophage variants resulting from two hyper-variable regions were identified by screening clone libraries of infected citrus, periwinkle and psyllids. Among them, Types A and B share highly conserved sequences and localize within the two prophages, FP1 and FP2, respectively. Although Types B and C were abundant in all three libraries, Type A was much more abundant in the libraries from the Las-infected psyllids than from the Las-infected plants, and Type D was only identified in libraries from the infected host plants but not from the infected psyllids. Sequence analysis of these variants revealed that the variations may result from recombination and rearrangement events. Conventional PCR results using type-specific molecular markers indicated that A, B, C and D are the four most abundant types in Las-infected citrus and periwinkle. However, only three types, A, B and C are abundant in Las-infected psyllids. Typing results for Las-infected citrus field samples indicated that mixed populations of Las bacteria present in Floridian isolates, but only the Type D population was correlated with the blotchy mottle symptom. Extended cloning and sequencing of the Type D region revealed a third prophage/phage in the Las genome, which may derive from the recombination of FP1 and FP2. Dramatic variations in these prophage regions were also found among the global Las isolates. These results are the first to demonstrate the prophage/phage-mediated dynamics of Las populations in plant and insect hosts, and their correlation with insect transmission and disease development.

Cluster M mycobacteriophages Bongo, PegLeg, and Rey with unusually large repertoires of tRNA isotypes

Genomic analysis of a large set of phages infecting the common host Mycobacterium smegmatis mc(2)155 shows that they span considerable genetic diversity. There are more than 20 distinct types that lack nucleotide similarity with each other, and there is considerable diversity within most of the groups. Three newly isolated temperate mycobacteriophages, Bongo, PegLeg, and Rey, constitute a new group (cluster M), with the closely related phages Bongo and PegLeg forming subcluster M1 and the more distantly related Rey forming subcluster M2. The cluster M mycobacteriophages have siphoviral morphologies with unusually long tails, are homoimmune, and have larger than average genomes (80.2 to 83.7 kbp). They exhibit a variety of features not previously described in other mycobacteriophages, including noncanonical genome architectures and several unusual sets of conserved repeated sequences suggesting novel regulatory systems for both transcription and translation. In addition to containing transfer-messenger RNA and RtcB-like RNA ligase genes, their genomes encode 21 to 24 tRNA genes encompassing complete or nearly complete sets of isotypes. We predict that these tRNAs are used in late lytic growth, likely compensating for the degradation or inadequacy of host tRNAs. They may represent a complete set of tRNAs necessary for late lytic growth, especially when taken together with the apparent lack of codons in the same late genes that correspond to tRNAs that the genomes of the phages do not obviously encode.

IMPORTANCE

The bacteriophage population is vast, dynamic, and old and plays a central role in bacterial pathogenicity. We know surprisingly little about the genetic diversity of the phage population, although metagenomic and phage genome sequencing indicates that it is great. Probing the depth of genetic diversity of phages of a common host, Mycobacterium smegmatis, provides a higher resolution of the phage population and how it has evolved. Three new phages constituting a new cluster M further expand the diversity of the mycobacteriophages and introduce novel features. As such, they provide insights into phage genome architecture, virion structure, and gene regulation at the transcriptional and translational levels.

Data-independent acquisition (MSE) with ion mobility provides a systematic method for analysis of a bacteriophage structural proteome

In this work, a method was developed to study the structural proteome of mycobacteriophage Marvin, a recent isolate from soil with 107 predicted coding sequences. This prototype method was applied for semi-quantitative analysis of the composition of this mycobacteriophage virion using ion mobility spectrometry and data-independent acquisition (MS(E)-IMS). MS(E)-IMS was compared to a more conventional proteomics technique employing mass spectrometry with a data-dependent acquisition strategy. MS(E)-IMS provided broad coverage of the virion proteome and high sequence coverage for individual proteins. This shotgun method does not depend on the limited sensitivity of visualization of protein bands by staining reagents inherent in gel-based methods. The method is comprehensive, provides high sequence coverage and is proposed as a particularly efficient method for the study of bacteriophage proteomes.

Twelve previously unknown phage genera are ubiquitous in global oceans

Viruses are fundamental to ecosystems ranging from oceans to humans, yet our ability to study them is bottlenecked by the lack of ecologically relevant isolates, resulting in "unknowns" dominating culture-independent surveys. Here we present genomes from 31 phages infecting multiple strains of the aquatic bacterium Cellulophaga baltica (Bacteroidetes) to provide data for an underrepresented and environmentally abundant bacterial lineage. Comparative genomics delineated 12 phage groups that (i) each represent a new genus, and (ii) represent one novel and four well-known viral families. This diversity contrasts the few well-studied marine phage systems, but parallels the diversity of phages infecting human-associated bacteria. Although all 12 Cellulophaga phages represent new genera, the podoviruses and icosahedral, nontailed ssDNA phages were exceptional, with genomes up to twice as large as those previously observed for each phage type. Structural novelty was also substantial, requiring experimental phage proteomics to identify 83% of the structural proteins. The presence of uncommon nucleotide metabolism genes in four genera likely underscores the importance of scavenging nutrient-rich molecules as previously seen for phages in marine environments. Metagenomic recruitment analyses suggest that these particular Cellulophaga phages are rare and may represent a first glimpse into the phage side of the rare biosphere. However, these analyses also revealed that these phage genera are widespread, occurring in 94% of 137 investigated metagenomes. Together, this diverse and novel collection of phages identifies a small but ubiquitous fraction of unknown marine viral diversity and provides numerous environmentally relevant phage-host systems for experimental hypothesis testing.

Generation of affinity-tagged fluoromycobacteriophages by mixed assembly of phage capsids

Addition of affinity tags to bacteriophage particles facilitates a variety of applications, including vaccine construction and diagnosis of bacterial infections. Addition of tags to phage capsids is desirable, as modification of the tails can lead to poor adsorption and loss of infectivity. Although tags can readily be included as fusions to head decoration proteins, many phages do not have decoration proteins as virion components. The addition of a small (10-amino-acid) Strep-tag II (STAG II) to the mycobacteriophage TM4 capsid subunit, gp9, was not tolerated as a genetically homogenous recombinant phage but could be incorporated into the head by growth of wild-type phage on a host expressing the capsid-STAG fusion. Particles with capsids composed of wild-type and STAG-tagged subunit mixtures could be grown to high titers, showed good infectivities, and could be used to isolate phage-bacterium complexes. Preparation of a STAG-labeled fluoromycobacteriophage enabled capture of bacterial complexes and identification of infected bacteria by fluorescence.

Generation of affinity-tagged fluoromycobacteriophages by mixed assembly of phage capsids

Addition of affinity tags to bacteriophage particles facilitates a variety of applications, including vaccine construction and diagnosis of bacterial infections. Addition of tags to phage capsids is desirable, as modification of the tails can lead to poor adsorption and loss of infectivity. Although tags can readily be included as fusions to head decoration proteins, many phages do not have decoration proteins as virion components. The addition of a small (10-amino-acid) Strep-tag II (STAG II) to the mycobacteriophage TM4 capsid subunit, gp9, was not tolerated as a genetically homogenous recombinant phage but could be incorporated into the head by growth of wild-type phage on a host expressing the capsid-STAG fusion. Particles with capsids composed of wild-type and STAG-tagged subunit mixtures could be grown to high titers, showed good infectivities, and could be used to isolate phage-bacterium complexes. Preparation of a STAG-labeled fluoromycobacteriophage enabled capture of bacterial complexes and identification of infected bacteria by fluorescence.

Cluster J mycobacteriophages: intron splicing in capsid and tail genes

Bacteriophages isolated on Mycobacterium smegmatis mc²155 represent many distinct genomes sharing little or no DNA sequence similarity. The genomes are architecturally mosaic and are replete with genes of unknown function. A new group of genomes sharing substantial nucleotide sequences constitute Cluster J. The six mycobacteriophages forming Cluster J are morphologically members of the Siphoviridae, but have unusually long genomes ranging from 106.3 to 117 kbp. Reconstruction of the capsid by cryo-electron microscopy of mycobacteriophage BAKA reveals an icosahedral structure with a triangulation number of 13. All six phages are temperate and homoimmune, and prophage establishment involves integration into a tRNA-Leu gene not previously identified as a mycobacterial attB site for phage integration. The Cluster J genomes provide two examples of intron splicing within the virion structural genes, one in a major capsid subunit gene, and one in a tail gene. These genomes also contain numerous free-standing HNH homing endonuclease, and comparative analysis reveals how these could contribute to genome mosaicism. The unusual Cluster J genomes provide new insights into phage genome architecture, gene function, capsid structure, gene mobility, intron splicing, and evolution.

The secret lives of mycobacteriophages

The study of mycobacteriophages provides insights into viral diversity and evolution, as well as the genetics and physiology of their pathogenic hosts. Genomic characterization of 80 mycobacteriophages reveals a high degree of genetic diversity and an especially rich reservoir of interesting genes. These include a vast number of genes of unknown function that do not match known database entries and many genes whose functions can be predicted but which are not typically found as components of phage genomes. Thus many mysteries surround these genomes, such as why the genes are there, what do they do, how are they expressed and regulated, how do they influence the physiology of the host bacterium, and what forces of evolution directed them to their genomic homes? Although the genetic diversity and novelty of these phages is full of intrigue, it is a godsend for the mycobacterial geneticist, presenting an abundantly rich toolbox that can be exploited to devise new and effective ways for understanding the genetics and physiology of human tuberculosis. As the number of sequenced genomes continues to grow, their mysteries continue to thicken, and the time has come to learn more about the secret lives of mycobacteriophages.

Phage cluster relationships identified through single gene analysis

Background

Phylogenetic comparison of bacteriophages requires whole genome approaches such as dotplot analysis, genome pairwise maps, and gene content analysis. Currently mycobacteriophages, a highly studied phage group, are categorized into related clusters based on the comparative analysis of whole genome sequences. With the recent explosion of phage isolation, a simple method for phage cluster prediction would facilitate analysis of crude or complex samples without whole genome isolation and sequencing. The hypothesis of this study was that mycobacteriophage-cluster prediction is possible using comparison of a single, ubiquitous, semi-conserved gene. Tape Measure Protein (TMP) was selected to test the hypothesis because it is typically the longest gene in mycobacteriophage genomes and because regions within the TMP gene are conserved.

Results

A single gene, TMP, identified the known Mycobacteriophage clusters and subclusters using a Gepard dotplot comparison or a phylogenetic tree constructed from global alignment and maximum likelihood comparisons. Gepard analysis of 247 mycobacteriophage TMP sequences appropriately recovered 98.8% of the subcluster assignments that were made by whole-genome comparison. Subcluster-specific primers within TMP allow for PCR determination of the mycobacteriophage subcluster from DNA samples. Using the single-gene comparison approach for siphovirus coliphages, phage groupings by TMP comparison reflected relationships observed in a whole genome dotplot comparison and confirm the potential utility of this approach to another widely studied group of phages.

Conclusions

TMP sequence comparison and PCR results support the hypothesis that a single gene can be used for distinguishing phage cluster and subcluster assignments. TMP single-gene analysis can quickly and accurately aid in mycobacteriophage classification.

Complete genome sequence analysis of Nocardia brasiliensis HUJEG-1 reveals a saprobic lifestyle and the genes needed for human pathogenesis

Nocardia brasiliensis is an important etiologic agent of mycetoma. These bacteria live as a saprobe in soil or organic material and enter the tissue via minor trauma. Mycetoma is characterized by tumefaction and the production of fistula and abscesses, with no spontaneous cure. By using mass sequencing, we determined the complete genomic nucleotide sequence of the bacteria. According to our data, the genome is a circular chromosome 9,436,348-bp long with 68% G+C content that encodes 8,414 proteins. We observed orthologs for virulence factors, a higher number of genes involved in lipid biosynthesis and catabolism, and gene clusters for the synthesis of bioactive compounds, such as antibiotics, terpenes, and polyketides. An in silico analysis of the sequence supports the conclusion that the bacteria acquired diverse genes by horizontal transfer from other soil bacteria, even from eukaryotic organisms. The genome composition reflects the evolution of bacteria via the acquisition of a large amount of DNA, which allows it to survive in new ecological niches, including humans.

The first structure of a mycobacteriophage, the Mycobacterium abscessus subsp. bolletii phage Araucaria

The unique characteristics of the waxy mycobacterial cell wall raise questions about specific structural features of their bacteriophages. No structure of any mycobacteriophage is available, although ∼3,500 have been described to date. To fill this gap, we embarked in a genomic and structural study of a bacteriophage from Mycobacterium abscessus subsp. bolletii, a member of the Mycobacterium abscessus group. This opportunistic pathogen is responsible for respiratory tract infections in patients with lung disorders, particularly cystic fibrosis. M. abscessus subsp. bolletii was isolated from respiratory tract specimens, and bacteriophages were observed in the cultures. We report here the genome annotation and characterization of the M. abscessus subsp. bolletii prophage Araucaria, as well as the first single-particle electron microscopy reconstruction of the whole virion. Araucaria belongs to Siphoviridae and possesses a 64-kb genome containing 89 open reading frames (ORFs), among which 27 could be annotated with certainty. Although its capsid and connector share close similarity with those of several phages from Gram-negative (Gram(-)) or Gram(+) bacteria, its most distinctive characteristic is the helical tail decorated by radial spikes, possibly host adhesion devices, according to which the phage name was chosen. Its host adsorption device, at the tail tip, assembles features observed in phages binding to protein receptors, such as phage SPP1. All together, these results suggest that Araucaria may infect its mycobacterial host using a mechanism involving adhesion to cell wall saccharides and protein, a feature that remains to be further explored.

Functional requirements for bacteriophage growth: gene essentiality and expression in mycobacteriophage Giles

Bacteriophages represent a majority of all life forms, and the vast, dynamic population with early origins is reflected in their enormous genetic diversity. A large number of bacteriophage genomes have been sequenced. They are replete with novel genes without known relatives. We know little about their functions, which genes are required for lytic growth, and how they are expressed. Furthermore, the diversity is such that even genes with required functions - such as virion proteins and repressors - cannot always be recognized. Here we describe a functional genomic dissection of mycobacteriophage Giles, in which the virion proteins are identified, genes required for lytic growth are determined, the repressor is identified, and the transcription patterns determined. We find that although all of the predicted phage genes are expressed either in lysogeny or in lytic growth, 45% of the predicted genes are non-essential for lytic growth. We also describe genes required for DNA replication, show that recombination is required for lytic growth, and that Giles encodes a novel repressor. RNAseq analysis reveals abundant expression of a small non-coding RNA in a lysogen and in late lytic growth, although it is non-essential for lytic growth and does not alter lysogeny.

Archaeal viruses, not archaeal phages: an archaeological dig

Viruses infect members of domains Bacteria, Eukarya, and Archaea. While those infecting domain Eukarya are nearly universally described as "Viruses", those of domain Bacteria, to a substantial extent, instead are called "Bacteriophages," or "Phages." Should the viruses of domain Archaea therefore be dubbed "Archaeal phages," "Archaeal viruses," or some other construct? Here we provide documentation of published, general descriptors of the viruses of domain Archaea. Though at first the term "Phage" or equivalent was used almost exclusively in the archaeal virus literature, there has been a nearly 30-year trend away from this usage, with some persistence of "Phage" to describe "Head-and-tail" archaeal viruses, "Halophage" to describe viruses of halophilic Archaea, use of "Prophage" rather than "Provirus," and so forth. We speculate on the root of the early 1980's transition from "Phage" to "Virus" to describe these infectious agents, consider the timing of introduction of "Archaeal virus" (which can be viewed as analogous to "Bacterial virus"), identify numerous proposed alternatives to "Archaeal virus," and also provide discussion of the general merits of the term, "Phage." Altogether we identify in excess of one dozen variations on how the viruses of domain Archaea are described, and document the timing of both their introduction and use.

Analysis of novel mycobacteriophages indicates the existence of different strategies for phage inheritance in mycobacteria

Mycobacteriophages have been essential in the development of mycobacterial genetics through their use in the construction of tools for genetic manipulation. Due to the simplicity of their isolation and variety of exploitable molecular features, we searched for and isolated 18 novel mycobacteriophages from environmental samples collected from several geographic locations. Characterization of these phages did not differ from most of the previously described ones in the predominant physical features (virion size in the 100–400 nm, genome size in the 50–70 kbp, morphological features compatible with those corresponding to the Siphoviridae family), however novel characteristics for propagation were noticed. Although all the mycobacteriophages propagated at 30°C, eight of them failed to propagate at 37°C. Since some of our phages yielded pinpoint plaques, we improved plaque detection by including sub-inhibitory concentrations of isoniazid or ampicillin-sulbactam in the culture medium. Thus, searches for novel mycobacteriophages at low temperature and in the presence of these drugs would allow for the isolation of novel members that would otherwise not be detected. Importantly, while eight phages lysogenized Mycobacterium smegmatis, four of them were also capable of lysogenizing Mycobacterium tuberculosis. Analysis of the complete genome sequence obtained for twelve mycobacteriophages (the remaining six rendered partial genomic sequences) allowed for the identification of a new singleton. Surprisingly, sequence analysis revealed the presence of parA or parA/parB genes in 7/18 phages including four that behaved as temperate in M. tuberculosis. In summary, we report here the isolation and preliminary characterization of mycobacteriophages that bring new information to the field.

Integration-dependent bacteriophage immunity provides insights into the evolution of genetic switches

Genetic switches are critical components of developmental circuits. Because temperate bacteriophages are vastly abundant and greatly diverse, they are rich resources for understanding the mechanisms and evolution of switches and the molecular control of genetic circuitry. Here, we describe a new class of small, compact, and simple switches that use site-specific recombination as the key decision point. The phage attachment site attP is located within the phage repressor gene such that chromosomal integration results in removal of a C-terminal tag that destabilizes the virally encoded form of the repressor. Integration thus not only confers prophage stability but also is a requirement for lysogenic establishment. The variety of these self-contained integration-dependent immunity systems in different genomic contexts suggests that these represent ancestral states in switch evolution from which more-complex switches have evolved. They also provide a powerful toolkit for building synthetic biological circuits.

Biology of a novel mycobacteriophage, SWU1, isolated from Chinese soil as revealed by genomic characteristics

Mycobacteriophage SWU1 is a newly isolated phage from a soil sample collected at Gongping village, Pingchang County, Sichuan Province, China, using Mycobacterium smegmatis mc(2)155 as a host. Plaques of SWU1 appear as a unique bull's-eye on an M. smegmatis lawn. In this paper, we report the complete genome sequence of SWU1 and some major findings from the analysis result.

Diversity in bacterial lysis systems: bacteriophages show the way

Bacteriophages have developed multiple host cell lysis strategies to promote release of descendant virions from infected bacteria. This review is focused on the lysis mechanisms employed by tailed double-stranded DNA bacteriophages, where new developments have recently emerged. These phages seem to use a least common denominator to induce lysis, the so-called holin-endolysin dyad. Endolysins are cell wall-degrading enzymes whereas holins form 'holes' in the cytoplasmic membrane at a precise scheduled time. The latter function was long viewed as essential to provide a pathway for endolysin escape to the cell wall. However, recent studies have shown that phages can also exploit the host cell secretion machinery to deliver endolysins to their target and subvert the bacterial autolytic arsenal to effectively accomplish lysis. In these systems the membrane-depolarizing holin function still seems to be essential to activate secreted endolysins. New lysis players have also been uncovered that promote degradation of particular bacterial cell envelopes, such as that of mycobacteria.

On the nature of mycobacteriophage diversity and host preference

The complete genome sequences of over 220 mycobacteriophages reveal them to be highly diverse, with numerous types sharing little or no nucleotide sequence identity with each other. We have determined the preferences of these phages for Mycobacterium tuberculosis and for other strains of Mycobacterium smegmatis, and find there is a correlation between genome type (cluster, subcluster, singleton) and host range. For many of the phages, expansion of host range occurs at relatively high frequencies, and we describe several examples in which host constraints occur at early stages of infection (adsorption or DNA injection), and phages have the ability to expand their host range through mutations in tail genes. We present a model in which phage diversity is a function of both the ability of phages to rapidly adapt to new hosts and the richness of the diversity of the bacterial population from which those phages are isolated.

Genetic engineering of Mycobacterium tuberculosis: a review

Genetic engineering has been used for decades to mutate and delete genes in the Mycobacterium tuberculosis genome with the translational goal of producing attenuated mutants with conserved susceptibility to antituberculous antibiotics. The development of plasmids and mycobacteriophages that can transfer DNA into the M. tuberculosis chromosome has effectively overcome M. tuberculosis slow growth rate and the capsule and mycolic acid wall, which limit DNA uptake. The use of genetic engineering techniques has shed light on many aspects of pathogenesis mechanisms, including cellular growth, mycolic acid biosynthesis, metabolism, drug resistance and virulence. Moreover, such research gave clues to the development of new vaccines or new drugs for routine clinical practice. The use of genetic engineering tools is mainly based on the underlying concept that altering or reducing the M. tuberculosis genome could decrease its virulence. A contrario, recent post-genomic analyses indicated that reduced bacterial genomes are often associated with increased bacterial virulence and that M. tuberculosis acquired genes by lateral genetic exchange during its evolution. Therefore, ancestors utilizing genetic engineering to add genes to the M. tuberculosis genome may lead to new vaccines and the availability of M. tuberculosis isolates with increased susceptibility to antituberculous antibiotics.

Bacteriophage based probes for pathogen detection

Rapid and specific detection of pathogenic bacteria is important for the proper treatment, containment and prevention of human, animal and plant diseases. Identifying unique biological probes to achieve a high degree of specificity and minimize false positives has therefore garnered much interest in recent years. Bacteriophages are obligate intracellular parasites that subvert bacterial cell resources for their own multiplication and production of disseminative new virions, which repeat the cycle by binding specifically to the host surface receptors and injecting genetic material into the bacterial cells. The precision of host recognition in phages is imparted by the receptor binding proteins (RBPs) that are often located in the tail-spike or tail fiber protein assemblies of the virions. Phage host recognition specificity has been traditionally exploited for bacterial typing using laborious and time consuming bacterial growth assays. At the same time this feature makes phage virions or RBPs an excellent choice for the development of probes capable of selectively capturing bacteria on solid surfaces with subsequent quick and automatic detection of the binding event. This review focuses on the description of pathogen detection approaches based on immobilized phage virions as well as pure recombinant RBPs. Specific advantages of RBP-based molecular probes are also discussed.

Mycobacteriophage endolysins: diverse and modular enzymes with multiple catalytic activities

The mycobacterial cell wall presents significant challenges to mycobacteriophages – viruses that infect mycobacterial hosts – because of its unusual structure containing a mycolic acid-rich mycobacterial outer membrane attached to an arabinogalactan layer that is in turn linked to the peptidoglycan. Although little is known about how mycobacteriophages circumvent these barriers during the process of infection, destroying it for lysis at the end of their lytic cycles requires an unusual set of functions. These include Lysin B proteins that cleave the linkage of mycolic acids to the arabinogalactan layer, chaperones required for endolysin delivery to peptidoglycan, holins that regulate lysis timing, and the endolysins (Lysin As) that hydrolyze peptidoglycan. Because mycobacterial peptidoglycan contains atypical features including 3→3 interpeptide linkages, it is not surprising that the mycobacteriophage endolysins also have non-canonical features. We present here a bioinformatic dissection of these lysins and show that they are highly diverse and extensively modular, with an impressive number of domain organizations. Most contain three domains with a novel N-terminal predicted peptidase, a centrally located amidase, muramidase, or transglycosylase, and a C-terminal putative cell wall binding domain.

Cangene Gold Medal Award Lecture — Genomic analysis and modification ofBurkholderia cepaciacomplex bacteriophages1This article is based on a presentation by Dr. Karlene Lynch at the 61st Annual Meeting of the Canadian Society of Microbiologists in St. John’s, Newfoundland and Labrador, on 21 June 2011. Dr. Lynch was the recipient of the 2011 Cangene Gold Medal as the Canadian Graduate Student Microbiologist of the Year, an annual award sponsored by Cangene Corporation intended to recognize excellence in graduate research

The Burkholderia cepacia complex (Bcc) is a group of 17 Gram-negative predominantly environmental bacterial species that cause potentially fatal opportunistic infections in cystic fibrosis (CF) patients. Although its prevalence in these individuals is lower than that of Staphylococcus aureus and Pseudomonas aeruginosa , the Bcc remains a serious problem in the CF community because of the pathogenicity, transmissibility, and inherent antibiotic resistance of these organisms. An alternative treatment for Bcc infections that is currently being developed is phage therapy, the clinical use of viruses that infect bacteria. To assess the suitability of individual phage isolates for therapeutic use, the complete genome sequences of a panel of Bcc‐specific phages were determined and analyzed. These sequences encode a broad range of proteins with a gradient of relatedness to phage and bacterial gene products from Burkholderia and other genera. The majority of these phages were found not to encode virulence factors, and despite their predominantly temperate nature, a proof-of-principle experiment has shown that they may be modified to a lytic form. Both the genomic characterization and subsequent engineering of Bcc‐specific phages are fundamental to the development of an effective phage therapy strategy for these bacteria.

Mycobacteriophage Marvin: a new singleton phage with an unusual genome organization

Mycobacteriophages represent a genetically diverse group of viruses that infect mycobacterial hosts. Although more than 80 genomes have been sequenced, these still poorly represent the likely diversity of the broader population of phages that can infect the host, Mycobacterium smegmatis mc(2)155. We describe here a newly discovered phage, Marvin, which is a singleton phage, having no previously identified close relatives. The 65,100-bp genome contains 107 predicted protein-coding genes arranged in a noncanonical genomic architecture in which a subset of the minor tail protein genes are displaced about 20 kbp from their typical location, situated among nonstructural genes anticipated to be expressed early in lytic growth. Marvin is not temperate, and stable lysogens cannot be recovered from infections, although the presence of a putative xis gene suggests that Marvin could be a relatively recent derivative of a temperate parent. The Marvin genome is replete with novel genes not present in other mycobacteriophage genomes, and although most are of unknown function, the presence of amidoligase and glutamine amidotransferase genes suggests intriguing possibilities for the interactions of Marvin with its mycobacterial hosts.

Genomics of bacterial and archaeal viruses: dynamics within the prokaryotic virosphere

Prokaryotes, bacteria and archaea, are the most abundant cellular organisms among those sharing the planet Earth with human beings (among others). However, numerous ecological studies have revealed that it is actually prokaryotic viruses that predominate on our planet and outnumber their hosts by at least an order of magnitude. An understanding of how this viral domain is organized and what are the mechanisms governing its evolution is therefore of great interest and importance. The vast majority of characterized prokaryotic viruses belong to the order Caudovirales, double-stranded DNA (dsDNA) bacteriophages with tails. Consequently, these viruses have been studied (and reviewed) extensively from both genomic and functional perspectives. However, albeit numerous, tailed phages represent only a minor fraction of the prokaryotic virus diversity. Therefore, the knowledge which has been generated for this viral system does not offer a comprehensive view of the prokaryotic virosphere. In this review, we discuss all families of bacterial and archaeal viruses that contain more than one characterized member and for which evolutionary conclusions can be attempted by use of comparative genomic analysis. We focus on the molecular mechanisms of their genome evolution as well as on the relationships between different viral groups and plasmids. It becomes clear that evolutionary mechanisms shaping the genomes of prokaryotic viruses vary between different families and depend on the type of the nucleic acid, characteristics of the virion structure, as well as the mode of the life cycle. We also point out that horizontal gene transfer is not equally prevalent in different virus families and is not uniformly unrestricted for diverse viral functions.

Cluster K mycobacteriophages: insights into the evolutionary origins of mycobacteriophage TM4

Five newly isolated mycobacteriophages –Angelica, CrimD, Adephagia, Anaya, and Pixie – have similar genomic architectures to mycobacteriophage TM4, a previously characterized phage that is widely used in mycobacterial genetics. The nucleotide sequence similarities warrant grouping these into Cluster K, with subdivision into three subclusters: K1, K2, and K3. Although the overall genome architectures of these phages are similar, TM4 appears to have lost at least two segments of its genome, a central region containing the integration apparatus, and a segment at the right end. This suggests that TM4 is a recent derivative of a temperate parent, resolving a long-standing conundrum about its biology, in that it was reportedly recovered from a lysogenic strain of Mycobacterium avium, but it is not capable of forming lysogens in any mycobacterial host. Like TM4, all of the Cluster K phages infect both fast- and slow-growing mycobacteria, and all of them – with the exception of TM4 – form stable lysogens in both Mycobacterium smegmatis and Mycobacterium tuberculosis; immunity assays show that all five of these phages share the same immune specificity. TM4 infects these lysogens suggesting that it was either derived from a heteroimmune temperate parent or that it has acquired a virulent phenotype. We have also characterized a widely-used conditionally replicating derivative of TM4 and identified mutations conferring the temperature-sensitive phenotype. All of the Cluster K phages contain a series of well conserved 13 bp repeats associated with the translation initiation sites of a subset of the genes; approximately one half of these contain an additional sequence feature composed of imperfectly conserved 17 bp inverted repeats separated by a variable spacer. The K1 phages integrate into the host tmRNA and the Cluster K phages represent potential new tools for the genetics of M. tuberculosis and related species.

Phamerator: a bioinformatic tool for comparative bacteriophage genomics

BACKGROUND

Bacteriophage genomes have mosaic architectures and are replete with small open reading frames of unknown function, presenting challenges in their annotation, comparative analysis, and representation.

RESULTS

We describe here a bioinformatic tool, Phamerator, that assorts protein-coding genes into phamilies of related sequences using pairwise comparisons to generate a database of gene relationships. This database is used to generate genome maps of multiple phages that incorporate nucleotide and amino acid sequence relationships, as well as genes containing conserved domains. Phamerator also generates phamily circle representations of gene phamilies, facilitating analysis of the different evolutionary histories of individual genes that migrate through phage populations by horizontal genetic exchange.

CONCLUSIONS

Phamerator represents a useful tool for comparative genomic analysis and comparative representations of bacteriophage genomes.

A second endolysin gene is fully embedded in-frame with the lysA gene of mycobacteriophage Ms6

Mycobacteriophages are dsDNA viruses that infect mycobacterial hosts. The mycobacteriophage Ms6 accomplishes lysis by producing two cell wall hydrolytic enzymes, Lysin A (LysA) that possesses a central peptidoglycan recognition protein (PGRP) super-family conserved domain with the amidase catalytic site, that cleaves the amide bond between the N-acetylmuramic acid and L-alanine residues in the oligopeptide crosslinking chains of the peptidoglycan and Lysin B (LysB) a mycolylarabinogalactan esterase that hydrolyzes the mycolic acids from the mycolyl-arabinogalactan-peptidoglycan complex. Examination of the endolysin (lysA) DNA sequence revealed the existence of an embedded gene (lysA₂₄₁) encoded in the same reading frame and preceded by a consensus ribosome-binding site. In the present work we show that, even though lysA is essential for Ms6 viability, phage mutants that express only the longer (Lysin₃₈₄) or the shorter (Lysin₂₄₁) endolysin are viable, but defective in the normal timing, progression and completion of host cell lysis. In addition, both endolysins have peptidoglycan hydrolase activity and demonstrated broad growth inhibition activity against various Gram-positive bacteria and mycobacteria.

Evolution of mosaically related tailed bacteriophage genomes seen through the lens of phage P22 virion assembly

The mosaic composition of the genomes of dsDNA tailed bacteriophages (Caudovirales) is well known. Observations of this mosaicism have generally come from comparisons of small numbers of often rather distantly related phages, and little is known about the frequency or detailed nature of the processes that generate this kind of diversity. Here we review and examine the mosaicism within fifty-seven clusters of virion assembly genes from bacteriophage P22 and its "close" relatives. We compare these orthologous gene clusters, discuss their surprising diversity and document horizontal exchange of genetic information between subgroups of the P22-like phages as well as between these phages and other phage types. We also point out apparent restrictions in the locations of mosaic sequence boundaries in this gene cluster. The relatively large sample size and the fact that phage P22 virion structure and assembly are exceptionally well understood make the conclusions especially informative and convincing.

Evolutionarily conserved orthologous families in phages are relatively rare in their prokaryotic hosts

We have identified conserved orthologs in completely sequenced genomes of double-strand DNA phages and arranged them into evolutionary families (phage orthologous groups [POGs]). Using this resource to analyze the collection of known phage genomes, we find that most orthologs are unique in their genomes (having no diverged duplicates [paralogs]), and while many proteins contain multiple domains, the evolutionary recombination of these domains does not appear to be a major factor in evolution of these orthologous families. The number of POGs has been rapidly increasing over the past decade, the percentage of genes in phage genomes that have orthologs in other phages has also been increasing, and the percentage of unknown "ORFans" is decreasing as more proteins find homologs and establish a family. Other properties of phage genomes have remained relatively stable over time, most notably the high fraction of genes that are never or only rarely observed in their cellular hosts. This suggests that despite the renowned ability of phages to transduce cellular genes, these cellular "hitchhiker" genes do not dominate the phage genomic landscape, and a large fraction of the genes in phage genomes maintain an evolutionary trajectory that is distinct from that of the host genes.