Mycobacteriophages BPs, Angel and Halo: comparative genomics reveals a novel class of ultra-small mobile genetic elements

Bacteriophage infection and killing of intracellular Mycobacterium abscessus

IMPORTANCE

As we rapidly approach a post-antibiotic era, bacteriophage (phage) therapy may offer a solution for treating drug-resistant bacteria. Mycobacterium abscessus is an emerging, multidrug-resistant pathogen that causes disease in people with cystic fibrosis, chronic obstructive pulmonary disease, and other underlying lung diseases. M. abscessus can survive inside host cells, a niche that can limit access to antibiotics. As current treatment options for M. abscessus infections often fail, there is an urgent need for alternative therapies. Phage therapy is being used to treat M. abscessus infections as an option of last resort. However, little is known about the ability of phages to kill bacteria in the host environment and specifically in an intracellular environment. Here, we demonstrate the ability of phages to enter mammalian cells and to infect and kill intracellular M. abscessus. These findings support the use of phages to treat intracellular bacterial pathogens.

Therapeutically useful mycobacteriophages BPs and Muddy require trehalose polyphleates

Mycobacteriophages show promise as therapeutic agents for non-tuberculous mycobacterium infections. However, little is known about phage recognition of Mycobacterium cell surfaces or mechanisms of phage resistance. We show here that trehalose polyphleates (TPPs)-high-molecular-weight, surface-exposed glycolipids found in some mycobacterial species-are required for infection of Mycobacterium abscessus and Mycobacterium smegmatis by clinically useful phages BPs and Muddy. TPP loss leads to defects in adsorption and infection and confers resistance. Transposon mutagenesis shows that TPP disruption is the primary mechanism for phage resistance. Spontaneous phage resistance occurs through TPP loss by mutation, and some M. abscessus clinical isolates are naturally phage-insensitive due to TPP synthesis gene mutations. Both BPs and Muddy become TPP-independent through single amino acid substitutions in their tail spike proteins, and M. abscessus mutants resistant to TPP-independent phages reveal additional resistance mechanisms. Clinical use of BPs and Muddy TPP-independent mutants should preempt phage resistance caused by TPP loss.

Virion glycosylation influences mycobacteriophage immune recognition

Glycosylation of eukaryotic virus particles is common and influences their uptake, trafficking, and immune recognition. In contrast, glycosylation of bacteriophage particles has not been reported; phage virions typically do not enter the cytoplasm upon infection, and they do not generally inhabit eukaryotic systems. We show here that several genomically distinct phages of Mycobacteria are modified with glycans attached to the C terminus of capsid and tail tube protein subunits. These O-linked glycans influence antibody production and recognition, shielding viral particles from antibody binding and reducing production of neutralizing antibodies. Glycosylation is mediated by phage-encoded glycosyltransferases, and genomic analysis suggests that they are relatively common among mycobacteriophages. Putative glycosyltransferases are also encoded by some Gordonia and Streptomyces phages, but there is little evidence of glycosylation among the broader phage population. The immune response to glycosylated phage virions in mice suggests that glycosylation may be an advantageous property for phage therapy of Mycobacterium infections.

Mycobacterium trehalose polyphleates are required for infection by therapeutically useful mycobacteriophages BPs and Muddy

Mycobacteriophages are good model systems for understanding their bacterial hosts and show promise as therapeutic agents for nontuberculous mycobacterium infections. However, little is known about phage recognition of Mycobacterium cell surfaces, or mechanisms of phage resistance. We show here that surface-exposed trehalose polyphleates (TPPs) are required for infection of Mycobacterium abscessus and Mycobacterium smegmatis by clinically useful phages BPs and Muddy, and that TPP loss leads to defects in adsorption, infection, and confers resistance. Transposon mutagenesis indicates that TPP loss is the primary mechanism for phage resistance. Spontaneous phage resistance occurs through TPP loss, and some M. abscessus clinical isolates are phage-insensitive due to TPP absence. Both BPs and Muddy become TPP-independent through single amino acid substitutions in their tail spike proteins, and M. abscessus mutants resistant to TPP-independent phages reveal additional resistance mechanisms. Clinical use of BPs and Muddy TPP-independent mutants should preempt phage resistance caused by TPP loss.

Mycobacterial nucleoid-associated protein Lsr2 is required for productive mycobacteriophage infection

Mycobacteriophages are a diverse group of viruses infecting Mycobacterium with substantial therapeutic potential. However, as this potential becomes realized, the molecular details of phage infection and mechanisms of resistance remain ill-defined. Here we use live-cell fluorescence microscopy to visualize the spatiotemporal dynamics of mycobacteriophage infection in single cells and populations, showing that infection is dependent on the host nucleoid-associated Lsr2 protein. Mycobacteriophages preferentially adsorb at Mycobacterium smegmatis sites of new cell wall synthesis and following DNA injection, Lsr2 reorganizes away from host replication foci to establish zones of phage DNA replication (ZOPR). Cells lacking Lsr2 proceed through to cell lysis when infected but fail to generate consecutive phage bursts that trigger epidemic spread of phage particles to neighbouring cells. Many mycobacteriophages code for their own Lsr2-related proteins, and although their roles are unknown, they do not rescue the loss of host Lsr2.

Complete Genome Sequences of Cluster F1 Mycobacteriophages Akhila and MilanaBonita

Akhila and MilanaBonita are mycobacteriophages that were isolated from soil in New York using Mycobacterium smegmatis. Both phages have genomes that are 56,251 bp long and contain 99 genes; the genomes differ by only 1 nucleotide. Based on gene content similarity to phages in the Actinobacteriophage Database, both phages are assigned to cluster F1.

Phage Therapy for Nontuberculous Mycobacteria: Challenges and Opportunities

Non-tuberculous mycobacterium (NTM) infections are often clinically challenging, with lengthy antibiotic regimens that fail to resolve the infections with few good outcomes remaining. Mycobacteriophages-viruses that infect Mycobacterium hosts-show promise as therapeutic agents for NTM infections and have been used in 20 compassionate use cases. Favorable outcomes were observed in many but not all cases, although the phages show exceptional safety profiles and no evidence of phage resistance was observed, even when only a single phage was administered. Phage-specific antibodies are commonly present following intravenous administration and are often neutralizing for the phage in vitro. However, phage neutralization does not consistently correlate with poor treatment outcomes and may not be a therapeutic limitation in all patients, even when immunocompetent. Currently, the therapeutic potential of phages is substantially limited by the great variation in phage susceptibility and a relatively small repertoire of therapeutically useful phages. As many as 45% of clinical isolates can have a smooth colony morphotype, and phages that both efficiently infect and kill these strains have yet to be described. In contrast, ~ 75% of rough strains are susceptible to and killed by one or more phages and therapeutic options can be considered on a compassionate use basis. Although therapies must currently be personalized, elucidating the determinants of phage host specificity, expanding the useful phage repertoire, and identifying the key determinants of clinical outcomes will reveal their full therapeutic potential.

Characterization of the cluster MabR prophages of Mycobacterium abscessus and Mycobacterium chelonae

Mycobacterium abscessus is an emerging pathogen of concern in cystic fibrosis and immunocompromised patients and is considered one of the most drug-resistant mycobacteria. The majority of clinical Mycobacterium abscessus isolates carry 1 or more prophages that are hypothesized to contribute to virulence and bacterial fitness. The prophage McProf was identified in the genome of the Bergey strain of Mycobacterium chelonae and is distinct from previously described prophages of Mycobacterium abscessus. The McProf genome increases intrinsic antibiotic resistance of Mycobacterium chelonae and drives expression of the intrinsic antibiotic resistance gene, whiB7, when superinfected by a second phage. The prevalence of McProf-like genomes was determined in sequenced mycobacterial genomes. Related prophage genomes were identified in the genomes of 25 clinical isolates of Mycobacterium abscessus and assigned to the novel cluster, MabR. They share less than 10% gene content with previously described prophages; however, they share features typical of prophages, including polymorphic toxin–immunity systems.

A Review on Mycobacteriophages: From Classification to Applications

Mycobacterial infections are a group of life-threatening conditions triggered by fast- or slow-growing mycobacteria. Some mycobacteria, such as Mycobacterium tuberculosis, promote the deaths of millions of lives throughout the world annually. The control of mycobacterial infections is influenced by the challenges faced in the diagnosis of these bacteria and the capability of these pathogens to develop resistance against common antibiotics. Detection of mycobacterial infections is always demanding due to the intracellular nature of these pathogens that, along with the lipid-enriched structure of the cell wall, complicates the access to the internal contents of mycobacterial cells. Moreover, recent studies depicted that more than 20% of M. tuberculosis (Mtb) infections are multi-drug resistant (MDR), and only 50% of positive MDR-Mtb cases are responsive to standard treatments. Similarly, the susceptibility of nontuberculosis mycobacteria (NTM) to first-line tuberculosis antibiotics has also declined in recent years. Exploiting mycobacteriophages as viruses that infect mycobacteria has significantly accelerated the diagnosis and treatment of mycobacterial infections. This is because mycobacteriophages, regardless of their cycle type (temperate/lytic), can tackle barriers in the mycobacterial cell wall and make the infected bacteria replicate phage DNA along with their DNA. Although the infectivity of the majority of discovered mycobacteriophages has been evaluated in non-pathogenic M. smegmatis, more research is still ongoing to find mycobacteriophages specific to pathogenic mycobacteria, such as phage DS6A, which has been shown to be able to infect members of the M. tuberculosis complex. Accordingly, this review aimed to introduce some potential mycobacteriophages in the research, specifically those that are infective to the three troublesome mycobacteria, M. tuberculosis, M. avium subsp. paratuberculosis (MAP), and M. abscessus, highlighting their theranostic applications in medicine.

Mycobacteriophages: From Petri dish to patient

Mycobacteriophages-bacteriophages infecting Mycobacterium hosts-contribute substantially to our understanding of viral diversity and evolution, provide resources for advancing Mycobacterium genetics, are the basis of high-impact science education programs, and show considerable therapeutic potential. Over 10,000 individual mycobacteriophages have been isolated by high school and undergraduate students using the model organism Mycobacterium smegmatis mc2155 and 2,100 have been completely sequenced, giving a high-resolution view of the phages that infect a single common host strain. The phage genomes are revealed to be highly diverse and architecturally mosaic and are replete with genes of unknown function. Mycobacteriophages have provided many widely used tools for Mycobacterium genetics including integration-proficient vectors and recombineering systems, as well as systems for efficient delivery of reporter genes, transposons, and allelic exchange substrates. The genomic insights and engineering tools have facilitated exploration of phages for treatment of Mycobacterium infections, although their full therapeutic potential has yet to be realized.

Increased whiB7 expression and antibiotic resistance in Mycobacterium chelonae carrying two prophages

BACKGROUND

The global rise in the incidence of non-tuberculosis mycobacterial infections is of increasing concern due their high levels of intrinsic antibiotic resistance. Although integrated viral genomes, called prophage, are linked to increased antibiotic resistance in some bacterial species, we know little of their role in mycobacterial drug resistance.

RESULTS

We present here for the first time, evidence of increased antibiotic resistance and expression of intrinsic antibiotic resistance genes in a strain of Mycobacterium chelonae carrying prophage. Strains carrying the prophage McProf demonstrated increased resistance to amikacin. Resistance in these strains was further enhanced by exposure to sub-inhibitory concentrations of the antibiotic, acivicin, or by the presence of a second prophage, BPs. Increased expression of the virulence gene, whiB7, was observed in strains carrying both prophages, BPs and McProf, relative to strains carrying a single prophage or no prophages.

CONCLUSIONS

This study provides evidence that prophage alter expression of important mycobacterial intrinsic antibiotic resistance genes and additionally offers insight into the role prophage may play in mycobacterial adaptation to stress.

Toward a Phage Cocktail for Tuberculosis: Susceptibility and Tuberculocidal Action of Mycobacteriophages against Diverse Mycobacterium tuberculosis Strains

The global health burden of human tuberculosis (TB) and the widespread antibiotic resistance of its causative agent Mycobacterium tuberculosis warrant new strategies for TB control. The successful use of a bacteriophage cocktail to treat a Mycobacterium abscessus infection suggests that phages could play a role in tuberculosis therapy. To assemble a phage cocktail with optimal therapeutic potential for tuberculosis, we have explored mycobacteriophage diversity to identify phages that demonstrate tuberculocidal activity and determined the phage infection profiles for a diverse set of strains spanning the major lineages of human-adapted strains of the Mycobacterium tuberculosis complex. Using a combination of genome engineering and bacteriophage genetics, we have assembled a five-phage cocktail that minimizes the emergence of phage resistance and cross-resistance to multiple phages, and which efficiently kills the M. tuberculosis strains tested. Furthermore, these phages function without antagonizing antibiotic effectiveness, and infect both isoniazid-resistant and -sensitive strains.IMPORTANCE Tuberculosis kills 1.5 million people each year, and resistance to commonly used antibiotics contributes to treatment failures. The therapeutic potential of bacteriophages against Mycobacterium tuberculosis offers prospects for shortening antibiotic regimens, provides new tools for treating multiple drug-resistant (MDR)-TB and extensively drug-resistant (XDR)-TB infections, and protects newly developed antibiotics against rapidly emerging resistance to them. Identifying a suitable suite of phages active against diverse M. tuberculosis isolates circumvents many of the barriers to initiating clinical evaluation of phages as part of the arsenal of antituberculosis therapeutics.

Complete Genome Sequences of Mycobacteriophages HarryOW and Peeb

HarryOW and Peeb are Mycobacterium smegmatis mc² 155 Siphoviridae temperate phages with 52,935 and 41,876 base pairs in genome length, respectively. HarryOW belongs to the A1 subcluster and Peeb to the G1 subcluster. They were isolated and annotated by students from the SUNY Old Westbury Science and Technology Entry Program.

HarryOW and Peeb are Mycobacterium smegmatis mc² 155 Siphoviridae temperate phages with 52,935 and 41,876 base pairs in genome length, respectively. HarryOW belongs to the A1 subcluster and Peeb to the G1 subcluster. They were isolated and annotated by students from the SUNY Old Westbury Science and Technology Entry Program.

Phage Therapy for Mycobacterium Abscessus and Strategies to Improve Outcomes

Members of Mycobacterium abscessus complex are known for causing severe, chronic infections. Members of M. abscessus are a new "antibiotic nightmare" as one of the most resistant organisms to chemotherapeutic agents. Treatment of these infections is challenging due to the either intrinsic or acquired resistance of the M. abscessus complex to the available antibiotics. Recently, successful phage therapy with a cocktail of three phages (one natural lytic phage and two engineered phages) every 12 h for at least 32 weeks has been reported against a severe case of the disseminated M. abscessus subsp. massiliense infection, which underlines the high value of phages against drug-resistant superbugs. This report also highlighted the limitations of phage therapy, such as the absence of lytic phages with a broad host-range against all strains and subspecies of the M. abscessus complex and also the risk of phage resistant bacteria over treatment. Cutting-edge genomic technologies have facilitated the development of engineered phages for therapeutic purposes by introducing new desirable properties, changing host-range and arming the phages with additional killing genes. Here, we review the available literature and suggest new potential solutions based on the progress in phage engineering that can help to overcome the present limitations of M. abscessus treatment.

Mycobacteriophages

Mycobacteriophages are viruses that infect mycobacterial hosts. A large number of mycobacteriophages have been isolated and genomically characterized, providing insights into viral diversity and evolution, as well as fueling development of tools for mycobacterial genetics. Mycobacteriophages have intimate relationships with their hosts and provide insights into the genetics and physiology of the mycobacteria and tools for potential clinical applications such as drug development, diagnosis, vaccines, and potentially therapy.

A Parallel Adder Coordinates Mycobacterial Cell-Cycle Progression and Cell-Size Homeostasis in the Context of Asymmetric Growth and Organization

In model bacteria, such as E. coli and B. subtilis, regulation of cell-cycle progression and cellular organization achieves consistency in cell size, replication dynamics, and chromosome positioning [1-3]. Mycobacteria elongate and divide asymmetrically, giving rise to significant variation in cell size and elongation rate among closely related cells [4, 5]. Given the physical asymmetry of mycobacteria, the models that describe coordination of cellular organization and cell-cycle progression in model bacteria are not directly translatable [1, 2, 6-8]. Here, we used time-lapse microscopy and fluorescent reporters of DNA replication and chromosome positioning to examine the coordination of growth, division, and chromosome dynamics at a single-cell level in Mycobacterium smegmatis (M. smegmatis) and Mycobacterium bovis Bacillus Calmette-Guérin (BCG). By analyzing chromosome and replisome localization, we demonstrated that chromosome positioning is asymmetric and proportional to cell size. Furthermore, we found that cellular asymmetry is maintained throughout the cell cycle and is not established at division. Using measurements and stochastic modeling of mycobacterial cell size and cell-cycle timing in both slow and fast growth conditions, we found that well-studied models of cell-size control are insufficient to explain the mycobacterial cell cycle. Instead, we showed that mycobacterial cell-cycle progression is regulated by an unprecedented mechanism involving parallel adders (i.e., constant growth increments) that start at replication initiation. Together, these adders enable mycobacterial populations to regulate cell size, growth, and heterogeneity in the face of varying environments.

A widespread family of polymorphic toxins encoded by temperate phages

Background

Polymorphic toxins (PTs) are multi-domain bacterial exotoxins belonging to distinct families that share common features in terms of domain organization. PTs are found in all major bacterial clades, including many toxic effectors of type V and type VI secretion systems. PTs modulate the dynamics of microbial communities by killing or inhibiting the growth of bacterial competitors lacking protective immunity proteins.

Results

In this work, we identified a novel widespread family of PTs, named MuF toxins, which were exclusively encoded within temperate phages and their prophages. By analyzing the predicted proteomes of 1845 bacteriophages and 2464 bacterial genomes, we found that MuF-containing proteins were frequently part of the DNA packaging module of tailed phages. Interestingly, MuF toxins were abundant in the human gut microbiome.

Conclusions

Our results uncovered the presence of the MuF toxin family in the temperate phages of Firmicutes. The MuF toxin family is likely to play an important role in the ecology of the human microbiota where pathogens and commensal species belonging to the Firmicutes are abundant. We propose that MuF toxins could be delivered by phages into host bacteria and either influence the lysogeny decision or serve as bacterial weapons by inhibiting the growth of competing bacteria.

Electronic supplementary material

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

Genome Sequences of Mycobacteriophages Jane and Sneeze, New Members of Cluster G

Jane and Sneeze are newly isolated phages of Mycobacterium smegmatis mc²155 from Hillsborough, NJ, and Palo Verde, Costa Rica, respectively. Both are cluster G, subcluster G1 mycobacteriophages. Notable nucleotide differences exist between genomes in the right half, including the presence of mycobacteriophage mobile element 1 (MPME1) in Jane.

Evolutionary interpretations of mycobacteriophage biodiversity and host-range through the analysis of codon usage bias

In an genomics course sponsored by the Howard Hughes Medical Institute (HHMI), undergraduate students have isolated and sequenced the genomes of more than 1,150 mycobacteriophages, creating the largest database of sequenced bacteriophages able to infect a single host, Mycobacterium smegmatis, a soil bacterium. Genomic analysis indicates that these mycobacteriophages can be grouped into 26 clusters based on genetic similarity. These clusters span a continuum of genetic diversity, with extensive genomic mosaicism among phages in different clusters. However, little is known regarding the primary hosts of these mycobacteriophages in their natural habitats, nor of their broader host ranges. As such, it is possible that the primary host of many newly isolated mycobacteriophages is not M. smegmatis, but instead a range of closely related bacterial species. However, determining mycobacteriophage host range presents difficulties associated with mycobacterial cultivability, pathogenicity and growth. Another way to gain insight into mycobacteriophage host range and ecology is through bioinformatic analysis of their genomic sequences. To this end, we examined the correlations between the codon usage biases of 199 different mycobacteriophages and those of several fully sequenced mycobacterial species in order to gain insight into the natural host range of these mycobacteriophages. We find that UPGMA clustering tends to match, but not consistently, clustering by shared nucleotide sequence identify. In addition, analysis of GC content, tRNA usage and correlations between mycobacteriophage and mycobacterial codon usage bias suggests that the preferred host of many clustered mycobacteriophages is not M. smegmatis but other, as yet unknown, members of the mycobacteria complex or closely allied bacterial species.

Different characteristics between mycobacteriophage Chy1 and D29, which were classified as cluster A2 mycobacteriophages

PURPOSE

The aim of this study was to isolate a novel mycobacteriophage and then explore its anti-tuberculosis (TB) potential.

MATERIALS AND METHODS

Phage was isolated from enriched soil sample. A total of 36 mycobacterial strains obtained from clinical specimens were subjected to investigate the host range of phage by the spot lysis assay. Biological characteristics were investigated through growth curve, host range and phage antimicrobial activity in vitro. Then, genome sequencing and further analysis were accomplished by using an ABI3730XL DNA sequencer and comparative genome, respectively.

RESULTS

A lytic mycobacteriophage (Chy1) was isolated and the plaque morphology was similar to D29. The genome of Chy1 was estimated to be about 47,198 base pair (bp) and strong similarity (97.4% identity) to D29, especially, the Chy1 gene 7 encoding holin which is considered as a clock controlling growth cycle of the corresponding phage, was identical (100% identity) to phage D29 gene 11, thus classifying Chy1 as a member of the cluster A2 family. However, to our surprise, Chy1 can infect a narrower range of host-mycobacterial strains than that of D29. The latent period of Chy1 was quite longer compared to D29. Moreover, Chy1 has a weaker ability to lyse Mycobacterium smegmatis compared to D29.

CONCLUSIONS

The sequence of Chy1 showed 97.4% homology with the genome sequence of D29, but there was a large difference in their biological characteristics. Overall, the results of this investigation indicate that Chy1 is not an ideal candidate for developing mycobacteriophage-based anti-TB therapies but for future researches to investigate the reason why biological characteristics of Chy1 and D29 were remarkably different.

Genome Sequences of Cluster G Mycobacteriophages Cambiare, FlagStaff, and MOOREtheMARYer

Mycobacteriophages Cambiare, FlagStaff, and MOOREtheMARYer are newly isolated phages of Mycobacterium smegmatis mc 2 155 recovered from soil samples in Pittsburgh, PA. All three genomes are closely related to cluster G mycobacteriophages but differ sufficiently in nucleotide sequence and gene content to warrant division of cluster G into several subclusters.

Mycobacteriophage-repressor-mediated immunity as a selectable genetic marker: Adephagia and BPs repressor selection

Mycobacteriophages provide an abundance of systems for use in mycobacterial genetics, including manipulation of Mycobacterium tuberculosis. Because of the dearth of antibiotic resistance cassettes and biosafety concerns in constructing recombinant virulent M. tuberculosis strains, we developed the use of mycobacteriophage-encoded repressor genes that can be selected in the presence of lytic versions of their cognate phages. The phage Adephagia repressor gene (43) was identified through its ability to confer immunity to Adephagia superinfection, together with the mapping of mutations in gene 43 that confer a clear-phage phenotype. Plasmid transformants containing either Adephagia 43 or the previously identified BPs repressor 33 can be readily selected following electroporation using engineered lytic derivatives of Adephagia and BPs, respectively. Selection is as efficient as antibiotic selection, can be used with either single-copy integration vectors or with extrachromosomal vectors, and works similarly in both Mycobacterium smegmatis and M. tuberculosis.

Comparative genomics of Cluster O mycobacteriophages

Mycobacteriophages--viruses of mycobacterial hosts--are genetically diverse but morphologically are all classified in the Caudovirales with double-stranded DNA and tails. We describe here a group of five closely related mycobacteriophages--Corndog, Catdawg, Dylan, Firecracker, and YungJamal--designated as Cluster O with long flexible tails but with unusual prolate capsids. Proteomic analysis of phage Corndog particles, Catdawg particles, and Corndog-infected cells confirms expression of half of the predicted gene products and indicates a non-canonical mechanism for translation of the Corndog tape measure protein. Bioinformatic analysis identifies 8-9 strongly predicted SigA promoters and all five Cluster O genomes contain more than 30 copies of a 17 bp repeat sequence with dyad symmetry located throughout the genomes. Comparison of the Cluster O phages provides insights into phage genome evolution including the processes of gene flux by horizontal genetic exchange.

Frequency and fitness consequences of bacteriophage φ6 host range mutations

Viruses readily mutate and gain the ability to infect novel hosts, but few data are available regarding the number of possible host range-expanding mutations allowing infection of any given novel host, and the fitness consequences of these mutations on original and novel hosts. To gain insight into the process of host range expansion, we isolated and sequenced 69 independent mutants of the dsRNA bacteriophage Φ6 able to infect the novel host, Pseudomonas pseudoalcaligenes. In total, we found at least 17 unique suites of mutations among these 69 mutants. We assayed fitness for 13 of 17 mutant genotypes on P. pseudoalcaligenes and the standard laboratory host, P. phaseolicola. Mutants exhibited significantly lower fitnesses on P. pseudoalcaligenes compared to P. phaseolicola. Furthermore, 12 of the 13 assayed mutants showed reduced fitness on P. phaseolicola compared to wildtype Φ6, confirming the prevalence of antagonistic pleiotropy during host range expansion. Further experiments revealed that the mechanistic basis of these fitness differences was likely variation in host attachment ability. In addition, using computational protein modeling, we show that host-range expanding mutations occurred in hotspots on the surface of the phage's host attachment protein opposite a putative hydrophobic anchoring domain.

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.

PHIRE and TWiV: Experiences in Bringing Virology to New Audiences

Virology encompasses a broad spectrum of topics touching upon many aspects of our everyday lives. However, appreciation of this impact is too often restricted to those who have specialized training and participate in virology research. The Phage Hunters Integrating Research and Education (PHIRE) program and the This Week in Virology (TWiV) podcast seek to bring virology to new audiences through two different approaches—direct involvement of undergraduates in discovering and genomically characterizing bacteriophages (PHIRE) and clear, accessible, and free discussions among experts of all topics in virology (TWiV). Here we discuss these two high-impact programs, the audiences that they serve, their broader impacts, and their future potential.

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.

Mutational analysis of the mycobacteriophage BPs promoter PR reveals context-dependent sequences for mycobacterial gene expression

The PR promoter of mycobacteriophage BPs directs early lytic gene expression and is under the control of the BPs repressor, gp33. Reporter gene fusions showed that PR has modest activity in an extrachromosomal context but has activity that is barely detectable in an integrated context, even in the absence of its repressor. Mutational dissection of PR showed that it uses a canonical -10 hexamer recognized by SigA, and mutants with mutations to the sequence 5'-TATAMT had the greatest activities. It does not contain a 5'-TGN-extended -10 sequence, although mutants with mutations creating an extended -10 sequence had substantially increased promoter activity. Mutations in the -35 hexamer also influenced promoter activity but were strongly context dependent, and similar substitutions in the -35 hexamer differentially affected promoter activity, depending on the -10 and extended -10 motifs. This warrants caution in the construction of synthetic promoters or the bioinformatic prediction of promoter activity. Combinations of mutations throughout PR generated a calibrated series of promoters for expression of stably integrated recombinant genes in both Mycobacterium smegmatis and M. tuberculosis, with maximal promoter activity being more than 2-fold that of the strong hsp60 promoter.

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.

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.

Characterization and classification of Bo4 as a cluster G mycobacteriophage that can infect and lyse M. tuberculosis

Mycobacteriophage therapy is a potential alternative treatment for Mycobacterium tuberculosis infection. Here, we further characterized a mycobacteriophage, Bo4, and evaluated its ability to infect and kill M. tuberculosis. We first found that Bo4 can infect M. tuberculosis and Mycobacterium smegmatis. The observed clear plaques created by Bo4 infection indicated that Bo4 might be a lytic phage able to lyse mycobacterial strains, which was confirmed by phage antimicrobial activity. Bo4 formed clear zones in a medium with pH values of 7.4 or 5.0, suggesting the possibility that Bo4 could lyse mycobacteria, such as M. tuberculosis, in blood as well as in lysosomal macrophages. Further investigation into the Bo4 genome revealed that Bo4 had a dsDNA genome. Moreover, Bo4 contained ~39,318 bp comprised of 66.76 % G+C content. Complete genome sequencing showed high nucleotide identity with cluster G mycobacteriophages, thus classifying Bo4 as a member of the cluster G family. Additionally, annotation of the Bo4 genome indicated that it was a lytic bacteriophage and did not contain any harmful genes that increased mycobacterial virulence or decreased human immunity. Overall, the results of investigation indicate that the Bo4 possesses the potential to destroy M. tuberculosis, making it a potentially useful tool for diagnosing and treating tuberculosis.

A broadly implementable research course in phage discovery and genomics for first-year undergraduate students

Engaging large numbers of undergraduates in authentic scientific discovery is desirable but difficult to achieve. We have developed a general model in which faculty and teaching assistants from diverse academic institutions are trained to teach a research course for first-year undergraduate students focused on bacteriophage discovery and genomics. The course is situated within a broader scientific context aimed at understanding viral diversity, such that faculty and students are collaborators with established researchers in the field. The Howard Hughes Medical Institute (HHMI) Science Education Alliance Phage Hunters Advancing Genomics and Evolutionary Science (SEA-PHAGES) course has been widely implemented and has been taken by over 4,800 students at 73 institutions. We show here that this alliance-sourced model not only substantially advances the field of phage genomics but also stimulates students' interest in science, positively influences academic achievement, and enhances persistence in science, technology, engineering, and mathematics (STEM) disciplines. Broad application of this model by integrating other research areas with large numbers of early-career undergraduate students has the potential to be transformative in science education and research training.

IMPORTANCE

Engagement of undergraduate students in scientific research at early stages in their careers presents an opportunity to excite students about science, technology, engineering, and mathematics (STEM) disciplines and promote continued interests in these areas. Many excellent course-based undergraduate research experiences have been developed, but scaling these to a broader impact with larger numbers of students is challenging. The Howard Hughes Medical Institute (HHMI) Science Education Alliance Phage Hunting Advancing Genomics and Evolutionary Science (SEA-PHAGES) program takes advantage of the huge size and diversity of the bacteriophage population to engage students in discovery of new viruses, genome annotation, and comparative genomics, with strong impacts on bacteriophage research, increased persistence in STEM fields, and student self-identification with learning gains, motivation, attitude, and career aspirations.

A broadly implementable research course in phage discovery and genomics for first-year undergraduate students

Engaging large numbers of undergraduates in authentic scientific discovery is desirable but difficult to achieve. We have developed a general model in which faculty and teaching assistants from diverse academic institutions are trained to teach a research course for first-year undergraduate students focused on bacteriophage discovery and genomics. The course is situated within a broader scientific context aimed at understanding viral diversity, such that faculty and students are collaborators with established researchers in the field. The Howard Hughes Medical Institute (HHMI) Science Education Alliance Phage Hunters Advancing Genomics and Evolutionary Science (SEA-PHAGES) course has been widely implemented and has been taken by over 4,800 students at 73 institutions. We show here that this alliance-sourced model not only substantially advances the field of phage genomics but also stimulates students' interest in science, positively influences academic achievement, and enhances persistence in science, technology, engineering, and mathematics (STEM) disciplines. Broad application of this model by integrating other research areas with large numbers of early-career undergraduate students has the potential to be transformative in science education and research training.

IMPORTANCE

Engagement of undergraduate students in scientific research at early stages in their careers presents an opportunity to excite students about science, technology, engineering, and mathematics (STEM) disciplines and promote continued interests in these areas. Many excellent course-based undergraduate research experiences have been developed, but scaling these to a broader impact with larger numbers of students is challenging. The Howard Hughes Medical Institute (HHMI) Science Education Alliance Phage Hunting Advancing Genomics and Evolutionary Science (SEA-PHAGES) program takes advantage of the huge size and diversity of the bacteriophage population to engage students in discovery of new viruses, genome annotation, and comparative genomics, with strong impacts on bacteriophage research, increased persistence in STEM fields, and student self-identification with learning gains, motivation, attitude, and career aspirations.

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.

Evolutionary relationships among actinophages and a putative adaptation for growth in Streptomyces spp

The genome sequences of eight Streptomyces phages are presented, four of which were isolated for this study. Phages R4, TG1, Hau3, and SV1 were isolated previously and have been exploited as tools for understanding and genetically manipulating Streptomyces spp. We also extracted five apparently intact prophages from recent Streptomyces spp. genome projects and, together with six phage genomes in the database, we analyzed all 19 Streptomyces phage genomes with a view to understanding their relationships to each other and to other actinophages, particularly the mycobacteriophages. Fifteen of the Streptomyces phages group into four clusters of related genomes. Although the R4-like phages do not share nucleotide sequence similarity with other phages, they clearly have common ancestry with cluster A mycobacteriophages, sharing many protein homologues, common gene syntenies, and similar repressor-stoperator regulatory systems. The R4-like phage Hau3 and the prophage StrepC.1 (from Streptomyces sp. strain C) appear to have hijacked a unique adaptation of the streptomycetes, i.e., use of the rare UUA codon, to control translation of the essential phage protein, the terminase. The Streptomyces venezuelae generalized transducing phage SV1 was used to predict the presence of other generalized transducing phages for different Streptomyces species.

Degeneration of a CRISPR/Cas system and its regulatory target during the evolution of a pathogen

CRISPR/Cas systems are bacterial RNA-guided endonuclease machineries that target foreign nucleic acids. Recently, we demonstrated that the Cas protein Cas9 controls gene expression and virulence in Francisella novicida by altering the stability of the mRNA for an immunostimulatory bacterial lipoprotein (BLP). Genomic analyses, however, revealed that Francisella species with increased virulence harbor degenerated CRISPR/Cas systems. We hypothesize that CRISPR/Cas degeneration removed a barrier against genome alterations, which resulted in enhanced virulence. Importantly, the BLP locus was also lost; likely a necessary adaptation in the absence of Cas9-mediated repression. CRISPR/Cas systems likely play regulatory roles in numerous bacteria, and these data suggest additional genomic changes may be required to maintain fitness after CRISPR/Cas loss in such bacteria, having important evolutionary implications.

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.

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.

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.

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.

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.

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.

Evolutionary genomics of a temperate bacteriophage in an obligate intracellular bacteria (Wolbachia)

Genome evolution of bacteria is usually influenced by ecology, such that bacteria with a free-living stage have large genomes and high rates of horizontal gene transfer, while obligate intracellular bacteria have small genomes with typically low amounts of gene exchange. However, recent studies indicate that obligate intracellular species that host-switch frequently harbor agents of horizontal transfer such as mobile elements. For example, the temperate double-stranded DNA bacteriophage WO in Wolbachia persistently transfers between bacterial coinfections in the same host. Here we show that despite the phage's rampant mobility between coinfections, the prophage's genome displays features of constraint related to its intracellular niche. First, there is always at least one intact prophage WO and usually several degenerate, independently-acquired WO prophages in each Wolbachia genome. Second, while the prophage genomes are modular in composition with genes of similar function grouping together, the modules are generally not interchangeable with other unrelated phages and thus do not evolve by the Modular Theory. Third, there is an unusual core genome that strictly consists of head and baseplate genes; other gene modules are frequently deleted. Fourth, the prophage recombinases are diverse and there is no conserved integration sequence. Finally, the molecular evolutionary forces acting on prophage WO are point mutation, intragenic recombination, deletion, and purifying selection. Taken together, these analyses indicate that while lateral transfer of phage WO is pervasive between Wolbachia with occasional new gene uptake, constraints of the intracellular niche obstruct extensive mixture between WO and the global phage population. Although the Modular Theory has long been considered the paradigm of temperate bacteriophage evolution in free-living bacteria, it appears irrelevant in phages of obligate intracellular bacteria.

Expanding the diversity of mycobacteriophages: insights into genome architecture and evolution

Mycobacteriophages are viruses that infect mycobacterial hosts such as Mycobacterium smegmatis and Mycobacterium tuberculosis. All mycobacteriophages characterized to date are dsDNA tailed phages, and have either siphoviral or myoviral morphotypes. However, their genetic diversity is considerable, and although sixty-two genomes have been sequenced and comparatively analyzed, these likely represent only a small portion of the diversity of the mycobacteriophage population at large. Here we report the isolation, sequencing and comparative genomic analysis of 18 new mycobacteriophages isolated from geographically distinct locations within the United States. Although no clear correlation between location and genome type can be discerned, these genomes expand our knowledge of mycobacteriophage diversity and enhance our understanding of the roles of mobile elements in viral evolution. Expansion of the number of mycobacteriophages grouped within Cluster A provides insights into the basis of immune specificity in these temperate phages, and we also describe a novel example of apparent immunity theft. The isolation and genomic analysis of bacteriophages by freshman college students provides an example of an authentic research experience for novice scientists.

Mycobacteriophages: genes and genomes

Viruses are powerful tools for investigating and manipulating their hosts, but the enormous size and amazing genetic diversity of the bacteriophage population have emerged as something of a surprise. In light of the evident importance of mycobacteria to human health--especially Mycobacterium tuberculosis, which causes tuberculosis--and the difficulties that have plagued their genetic manipulation, mycobacteriophages are especially appealing subjects for discovery, genomic characterization, and manipulation. With more than 70 complete genome sequences available, the mycobacteriophages have provided a wealth of information on the diversity of phages that infect a common bacterial host, revealed the pervasively mosaic nature of phage genome architectures, and identified a huge number of genes of unknown function. Mycobacteriophages have provided key tools for tuberculosis genetics, and new methods for simple construction of mycobacteriophage recombinants will facilitate postgenomic explorations into mycobacteriophage biology.

Comparative genomic analysis of 60 Mycobacteriophage genomes: genome clustering, gene acquisition, and gene size

Mycobacteriophages are viruses that infect mycobacterial hosts. Expansion of a collection of sequenced phage genomes to a total of 60-all infecting a common bacterial host-provides further insight into their diversity and evolution. Of the 60 phage genomes, 55 can be grouped into nine clusters according to their nucleotide sequence similarities, 5 of which can be further divided into subclusters; 5 genomes do not cluster with other phages. The sequence diversity between genomes within a cluster varies greatly; for example, the 6 genomes in Cluster D share more than 97.5% average nucleotide similarity with one another. In contrast, similarity between the 2 genomes in Cluster I is barely detectable by diagonal plot analysis. In total, 6858 predicted open-reading frames have been grouped into 1523 phamilies (phams) of related sequences, 46% of which possess only a single member. Only 18.8% of the phams have sequence similarity to non-mycobacteriophage database entries, and fewer than 10% of all phams can be assigned functions based on database searching or synteny. Genome clustering facilitates the identification of genes that are in greatest genetic flux and are more likely to have been exchanged horizontally in relatively recent evolutionary time. Although mycobacteriophage genes exhibit a smaller average size than genes of their host (205 residues compared with 315), phage genes in higher flux average only 100 amino acids, suggesting that the primary units of genetic exchange correspond to single protein domains.

Defects in glycopeptidolipid biosynthesis confer phage I3 resistance in Mycobacterium smegmatis

Mycobacteriophages have played an important role in the development of genetic tools and diagnostics for pathogenic mycobacteria, including Mycobacterium tuberculosis. However, despite the isolation of numerous phages that infect mycobacteria, the mechanisms of mycobacteriophage infection remain poorly understood, and knowledge about phage receptors is minimal. In an effort to identify the receptor for phage I3, we screened a library of Mycobacterium smegmatis transposon mutants for phage-resistant strains. All four phage I3-resistant mutants isolated were found to have transposon insertions in genes located in a cluster involved in the biosynthesis of the cell-wall-associated glycopeptidolipid (GPL), and consequently the mutants did not synthesize GPLs. The loss of GPLs correlated specifically with phage I3 resistance, as all mutants retained sensitivity to two other mycobacteriophages: D29 and Bxz1. In order to define the minimal receptor for phage I3, we then tested the phage sensitivity of previously described GPL-deficient mutants of M. smegmatis that accumulate biosynthesis intermediates of GPLs. The results indicated that, while the removal of most sugar residues from the fatty acyl tetrapeptide (FATP) core of GPL did not affect sensitivity to phage I3, a single methylated rhamnose, transferred by the rhamnosyltransferase Gtf2 to the FATP core, was critical for phage binding.