key: cord-0907753-7m53i1h9
authors: Pope, Welkin H.; Jacobs-Sera, Deborah; Russell, Daniel A.; Rubin, Daniel H. F.; Kajee, Afsana; Msibi, Zama N. P.; Larsen, Michelle H.; Jacobs, William R.; Lawrence, Jeffrey G.; Hendrix, Roger W.; Hatfull, Graham F.
title: Genomics and Proteomics of Mycobacteriophage Patience, an Accidental Tourist in the Mycobacterium Neighborhood
date: 2014-12-02
journal: mBio
DOI: 10.1128/mbio.02145-14
sha: e100d3d45507da10018bc2b410c93072756fef47
doc_id: 907753
cord_uid: 7m53i1h9

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.

phages (average, 57.3% GC). The genome diversity and percent GC range of the mycobacteriophages support a model for phage genome evolution in which phages migrate much more rapidly across a diverse bacterial landscape than their genomes ameliorate toward that of any one host (4) . Thus, mycobacteriophages at the lower end of the percent GC spectrum are predicted to have infected lower-GC hosts in their recent evolutionary pasts (4) . We note that there are other phage-host systems with mismatched GC contents, such as T4 and the right arm of phage lambda (35.3% and 44.4%, respectively), relative to their Escherichia coli host (50.8%).

There is also substantial variation in tRNA content of the mycobacteriophages. Many have no tRNA genes, and some have only one or a small number, while others-such as the members of clusters C and M-have more than 20 (7, 8) . The rationale for carriage of these is unclear, and there is no obvious correlation between tRNA content and percent GC that might reflect tRNA acquisition to augment gene expression in newly acquired hosts. It has been noted that there is no close correlation between the tRNA specificities encoded by D29 and infrequently used codons (9) or between Bxz1 and phage codon preferences of putative highexpression genes (10) . Analysis of the codon usage of 32 mycobacteriophage genomes showed that there is variation in codon usage preferences (11) . In some phage genomes, the tRNAs may counteract host measures to protect themselves from infection by tRNA destruction (12) .

Here, we describe mycobacteriophage Patience, a newly isolated phage of M. smegmatis mc 2 155 that is a singleton with no close relatives and has a GC content of 50.3%, representing the extreme low end of the percent GC spectrum for mycobacteriophages. Of the predicted 109 Patience protein-coding genes, 61% are "orphams" with no close mycobacteriophage relatives in the Phamerator_285 database, and most of the 48 genes with homologues in other mycobacteriophages are most closely related to those in clusters H, R, and D, which also have below-average GC contents. However, Patience has a distinctly different codon usage profile from both its host and other mycobacteriophages and an abundance of codons that are rarely used in the host. Proteomic analysis using mass spectrometry provides evidence for expression of at least 83 Patience proteins, two of which are from cryptic open reading frames (ORFs) embedded within annotated genes. We propose that Patience is a relatively recent visitor to the Mycobacterium neighborhood, having evolved primarily in lower-GC hosts within the Actinomycetales, and is in the process of adapting to growth in its new high-GC genetic environment.

Isolation and genome sequencing of mycobacteriophage Patience. Mycobacteriophage Patience was isolated by direct plating of an environmental sample taken from near the Nelson Mandela School of Medicine at the University of KwaZulu-Natal (UKZN), Durban, South Africa, using M. smegmatis mc 2 155 as a host (3) . Isolation and purification of Patience were components of a 2-week workshop on mycobacterial genetics offered in July 2009; the genome was sequenced at the University of Pittsburgh and annotated in a second 2-week workshop at UKZN in July 2011. Patience forms normal-size hazy (~1-mm-diameter) plaques on M. smegmatis mc 2 155 at 37°C under standard conditions, although we have been unsuccessful in recovering stable lysogens. It can easily be propagated on solid media to titers greater than 10 10 PFU/ml. The genome is 70,506 bp long, circularly permuted, and presumably terminally redundant. For linear presentation, coordinate 1 is designated the beginning of an open reading frame upstream of the large terminase subunit consistent with the organization of the cluster H phages that contain homologues of the first open reading frame at their left ends (7) . The GenBank submission (JN412589) has been reported previously (3) . Genome annotation identified 109 putative open reading frames (ORFs) and one tRNA gene (Table 1) ; one additional ORF was annotated using mass spectrometry analysis (see below).

Virion morphology and virion PAGE analysis. Mycobacteriophage Patience is morphologically a member of the Siphoviridae with an isometric head and a long flexible tail (Fig. 1A) . The diameter of the heads averages 58 nm, and the tails average 358 nm in length, among the longest of the mycobacteriophages described to date, similar to those of the cluster M phages (348 nm [8] ); other phages with notably long tails are those in cluster H (approximately 290 nm) and cluster R (approximately 288 nm) (13) . SDS-PAGE analysis of Patience particles shows three abundant proteins and at least 12 other lower-abundance proteins between the sizes of 15 and 150 kDa (Fig. 1B) .

Relationship of Patience to other mycobacteriophages. The Patience genome is not closely related to other phage genomes, although dot plot analysis shows weak similarity to cluster H phages (Fig. 2) . The related segments span only 14%, 8%, and 6% of genome lengths with Barnyard, Predator, and Konstantine genomes, respectively (Fig. 2) , and the closest matching segment of Patience and Barnyard is a 1,091-bp region (73% nucleotide identity) corresponding to the capsid subunit genes. Alignment of the genome maps of Patience, Barnyard, Konstantine, and Predator in Phamerator shows their architectural relationships (Fig. 3) , and although the nucleotide sequence similarity is minimal, 37% of Patience genes have homologues in cluster H phages (Fig. 3) . However, more than half of the Patience ORFs are orphams (genes with no close mycobacteriophage relatives [ Fig. 3 ; Table 1 ]).

The overall GC content of the Patience genome is 50.3% and is maintained almost throughout the genome (Fig. 4A ). There are two notable departures suggesting recent acquisitions by horizontal exchange. One is within Patience gene 37, where the 3= end has 81% nucleotide identity to Rosebush gene 32 (subcluster B2) and the elevated percent GC reflects that of the B2 phages (Fig. 4B ). Patience gp37 is related to tail fibers of many other mycobacteriophages that are implicated in host range determination (4). The second example is gene 47 (62% GC [ Fig. 4C ]), although it is an orpham with no close mycobacteriophage relatives and is of unknown function.

Patience genome organization. Genome annotation of Patience indicates that all open reading frames and one tRNA gene are transcribed in the same direction ( Fig. 5 ; Table 1 ). The overall coding capacity is 94.9%, and there are no intergenic spaces greater than 250 bp. The recognizable virion structure and assembly functions lie within the gene 4 to 39 region, and the siphoviral syntenic arrangement of terminase, portal, protease, capsid, headtail connector proteins, major tail subunit, tail assembly chaperones, tape measure, and minor tail protein genes is observed, spanning 32 kbp of the genome (14) . This is atypically long not only because of the long tape measure gene corresponding to the long phage tail but because of a dozen genes of mostly unknown function located between the terminase large subunit (gene 6) and portal (gene 8) genes, between the portal (gene 8) and protease (gene 19) genes, and between the MuF-like gene (gene 20) and the putative scaffolding gene (gene 22) (Fig. 5) ; most of these are orphams (i.e., they do not have a close mycobacteriophage rela-tive), with the exceptions of genes 10 and 18 (Fig. 5 ). Gene 7 encodes a putative Endo VII protein and has weak similarity (Ͻ30% identity) to genes in phages Konstantine and Predator Table S1 in the supplemental material.

(subcluster H1), as well as Dori (singleton). The gene upstream of the terminase large subunit gene encodes a protein with similarities to putative HNH homing nucleases ( Fig. 5 ) but may function as part of the DNA packaging machinery as described for HK97 gp74 (15) . The four leftmost genes (genes 1 to 4) are of unknown function. The lysis cassette is located downstream of the virion structural genes and contains an endolysin gene (gene 41, lysin A) and a putative mycolylarabinogalactan esterase (gene 42, lysin B) (16, 17) . Genes 43 and 44 code for proteins containing two and four predicted transmembrane domains, respectively, that may act as holins or lysis chaperones (18) . The lysis cassette is not closely related to that of the cluster H phages, and lysin A is most closely related to the corresponding genes of cluster M phage PegLeg gp35 and Bongo gp35 (35% identity), with an Org-U domain organization (19) . In contrast, lysin B is most closely related to phages in subcluster B3 (66% identity to Gadget gp48), illustrating the mosaic nature of the lysis cassette (19) . Patience gp43 has no homologues, but gp44 has distant relatives in the cluster D phages (Ͻ30% identity), where the genes are located between lysin A and lysin B.

Of the other 65 Patience ORFs, only seven can be assigned putative functions, most of which are involved in DNA metabolism (helicases, RecA, DNA polymerase III [Pol III] alpha subunit, The three most abundant proteins likely correspond to the major tail subunit (gp31), the capsid subunit (gp23), and gp15, with predicted molecular masses of 30.5 kDa, 35.3 kDa (after processing), and 32 kDa, respectively. The major tail subunit may migrate slower than its predicted mass, as observed in some other phages (37) . Numbers at right are molecular masses in kilodaltons. Table 1 ]). These all have homologues in a variety of mycobacteriophages as well as in phages of Gordonia and Corynebacterium. One (gene 103) encodes a MazG-like protein, a putative nucleoside triphosphate pyro-phosphohydrolase that is common in a variety of phage genomes (20) . Patience has a single tRNA gene, and the tRNA is predicted to be charged with glutamine and has the anticodon 5=-UUG [i.e., tRNA Gln (UUG)]. It is unclear how the Patience genes are tran- (35) and the database Mycobacteriophage_285, containing 285 complete genome sequences, and aligned by the left ends of the tape measure gene. Maps are displayed in three tiers with pairwise nucleotide sequence similarities displayed in spectrum coloring, with violet being the most similar and red the least similar (minimal BLASTN cutoff E value is 10 Ϫ4 ). Genes are shown as colored boxes, with colors reflecting pham assignments for each gene (gene members of a phamily have the same color). Genes shown in white are orphams and have no mycobacteriophage relatives with greater than 32.5% amino acid identity or a BLASTP E value lower than 10 Ϫ50 . Phamily assignments with the number of phamily members in parentheses are above each gene. scribed, and there are no strongly predicted SigA-like promoters, similar to other mycobacteriophages such as Giles (21), even though SigA-like promoters are active in other mycobacteriophages (22) (23) (24) . Table S1 in the supplemental material). Using stringent criteria for peptide identification, we identified 79 of the 109 previously annotated predicted gene products (82 different gene products, including previously unannotated genes) present in at least one of the samples, including 34 of the 35 putative virion structure and assembly proteins (genes 5 to 39), 27 of which are particle associated (Fig. 5) . Although gp9 was not detected, particle-associated peptides corresponding to a previously unannotated 37-codon gene between genes 8 and 9 were identified, which we designate gene 111 (Fig. 5) . Particles contain four additional products, gp1, gp4, gp79, and gp86, although gp86 is abundant in lytic growth and the few particle-associated peptides could be contaminants in the phage preparation (Table S1 ). The capsid subunit (gp23) appears to be proteolytically cleaved between residues 82 and 83 to generate a 35.3-kDa mature protein product; that and the gp31 major tail subunit (30.5 kDa) and gp15 (32 kDa) likely correspond to the three major species observed by SDS-PAGE (Fig. 1B) . It is noteworthy that-with the exception of 31-residue gp9 -all of the genes corresponding to insertions within the otherwise canonically syntenic virion structure and assembly operon are present in Patience particles, although only 3 spectra of gp21 were identified and could represent contaminants from the lysate (Fig. 5 ; Table S1).

Of the total of 81 Patience proteins identified in infected cells, only one-gp47-was not identified in the late-infected (2.5-hpostinfection) sample. Patience gp47 was identified in the earlyinfected sample and is presumably not expressed late in infection but also is turned over rapidly ( Fig. 5 ; see also Table S1 in the supplemental material). Many of the virion proteins (genes 1 to 40) are expressed at higher levels (using peptide spectral count as a surrogate for expression) in the late-infected sample, although there is also substantial expression of many of these at the early time. Of the 29 annotated proteins not identified in any of the samples, 11 are predicted to be smaller than 10 kDa and may have escaped detection if they are not expressed at high levels. Twelve are predicted to be membrane or wall associated and are likely excluded from the soluble fractions used for analysis. Peptide abundance in infected cells is expected to generally correlate with expression, although other factors such as protein size and trypsin cleavage efficiencies influence peptide detection.

First, we note that peptides at or near the N terminus confirm the previously annotated translation start site for 54 proteins but also identify seven (gp4, gp17, gp29, gp47, gp53, gp89, and gp101) for which revisions of the annotated start sites are supported. For 16 proteins, the peptide coverage is insufficient to be informative about start site usage. Interestingly, seven proteins, all of them particle associated (gp1, gp20, gp21, gp26, gp31, gp37, and gp79), are acetylated at their N terminus (following methionine loss), although the functional significance-if any-is not known. Six are acetylated at an N-terminal threonine following methionine removal (the seventh is a serine acetylation), although not all proteins with N-terminal threonine residues are acetylated. gp21 is unusual in that the N-terminal-most peptides start at residue 35 of the annotated product following a glycine in the Ϫ1 position, indicating that they were not generated either by tryptic digestion or by translation initiation. Other processes such as posttranslational processing or intron splicing prior to translation may be involved. We also note that gp24 is present in the particles but that the peptide coverage (62 total spectra) is restricted to the N-terminal 50% of the protein, whereas in infected cells, the spectra reflect 100% coverage of the protein (149 total spectra). This could be explained by assembly-associated protein processing. Surprisingly, we identified several peptides corresponding to translation of regions wholly embedded within annotated genes. Two of these were identified only with peptides for which we have lower confidence and will not be considered further. The evidence supporting translation of the other two is, however, quite strong. The first of these is a 372-bp open reading frame transcribed on the same strand, in a different reading frame, within gene 91, the DNA polymerase III catalytic subunit (see Fig. S1A and S2A in the supplemental material). A total of 17 instances of seven different peptides were identified, and the spectra and fragmentation tables strongly support the peptide assignments (see Text S1 in the supplemental material). Moreover, 10 peptides correspond to the extreme N terminus with the methionine present, and there is a strong ribosome binding site upstream (Fig. S2) . Although the protein is seemingly expressed, the functional relevance is unclear. The predicted product has no close database relatives, and the evidence for conservation is ambiguous (Fig. S1A ). Barnyard contains a similar open reading frame within its polymerase, and the products share 30% amino acid identity, whereas the corresponding segments of the polymerases share only 53% identity. This region of the polymerase is more distantly related in Konstantine (40% identity), and the second open reading frame is not conserved (Fig. S1A) .

A similar scenario is seen within Patience gene 20, where two peptides corresponding to a 171-bp open reading frame on the same strand were identified with high confidence (see Fig. S1B in the supplemental material). Peptides corresponding to the predicted N terminus were not identified, although a start codon with a strong ribosome binding site is present (Fig. S2B) . The corresponding segment of Patience gp20 is conserved in Barnyard gp18 and Konstantine gp17, but the second open reading frame is not. Presumably, either expression of embedded out-of-frame genes is more common than anticipated or these reflect expression artifacts resulting from the use of noncognate host transcription and translation systems.

Codon selection and adaptation for mycobacterial growth. Because M. smegmatis has a relatively high GC content (67.4%), it is not surprising that these mutational biases (reflected in base composition) result in frequent usage of GC base pairs in thirdcodon positions (Fig. 6A ; see also Fig. S3 and S4 in the supplemental material). Each of the most commonly used of the synonymous codons has GC in the third position, and M. smegmatis mc 2 155 carries tRNAs with anticodons corresponding to all codons with GC in the third position, except for CGC (Fig. S3) . All of the NNU codons are present infrequently and (with the exception of CGU) are decoded by wobble pairing in the codon third position (Fig. S3  and S4) ; M. tuberculosis H37Rv has an almost identical tRNA profile.

Mycobacteriophages such as Twister (65% GC) and KayaCho (70% GC), with nucleotide compositions similar to that of M. smegmatis, have similar synonymous codon distributions (see Fig. S3 and S4 in the supplemental material). In contrast, Patience codon usage patterns are distinctly different, with notably different distributions of codons with respect to their third-position content relative to M. smegmatis ( Fig. 6A ; see also Fig. S3 and S4 ). There are a total of nine switches to a different most commonly used synonymous codon from M. smegmatis, seven of which are NNU codons that are rarely used in the host (five with normalized Fig. S3] ). In total, NNU codons (except CGU, for which there is a cognate host tRNA) represent over 25% of all Patience codons; in contrast, these are only 5% of M. smegmatis codons. This represents an extreme end of a trend between the usage of NNU codons and the overall GC content of the phage, where usage of NNU increases as percent GC decreases (Fig. S5 ).

Although codon usage shows that Patience spent a significant portion of its evolutionary past in hosts with moderate GC contents, it is not clear what its current host range is. If Patience currently exploits more GC-rich hosts, then it would experience not only mutational biases, which would increase its GC content, but selective pressure on its most highly expressed genes to use tRNA pools poised to translate GC-rich genes. This selective pressure would be reflected in a preferred usage of GC-rich codons within highly expressed genes. Patience does experience codon selection, which is shown by the robust positive correlation between codon selection (adapative codon enrichment [ACE z ] [25] ) (see Materials and Methods) and level of gene expression, using total numbers of peptides as identified by mass spectrometry under stringent conditions as a surrogate for expression (Fig. 6B) . Codon selection is measured using ␦ values, or the ratio of codon frequencies in genes experiencing codon selection normalized to their frequencies in genes lacking codon selection. Preferred codons are those with ␦ values greater than 1, indicating more frequent use in genes experiencing codon selection, frequently those expressed to greater levels. Patience favors the use of GCrich codons in genes experiencing codon selection (Table 2) , with 24 of 33 preferred codons bearing GC base pairs in their third positions, whereas 21 of 26 nonpreferred codons end with AT base pairs. Moreover, the patterns of codon preference in Patience (which codons are preferred and which are not) are strongly correlated with that of M. smegmatis (Fig. 6C) , suggesting that this host imposes codon selection congruent with that currently being experienced by Patience. In contrast, Patience codon selection is not congruent with that imposed by hosts with more moderate GC content, such as Escherichia coli (Fig. 6C, inset) , although such surveys can never be conclusive.

Patience represents an intriguing example of a virus that has successfully entered the mycobacterial genetic neighborhood in its relatively recent evolutionary history. Its overall GC content and codon usage profiles are distinctly different from those of its mycobacterial host, suggesting that it primarily evolved in a moderate-GC (~50%) environment. Growth in high-GC bacteria may have required multiple events, including acquisition of part of a tail gene (gene 47) by lateral gene transfer. Nonetheless, the mismatch between viral and host genomic profiles does not appear to have been a substantial impediment to host range expan-sion, although the highly expressed viral genes are under codon selection for more efficient translation by the host apparatus. Interestingly, although Patience conceivably could have responded by acquisition of a tRNA repertoire to facilitate phage gene expression, this has not occurred; the only phage-carried tRNA is tRNA-Gln (UUG), and CAA is not a rare codon in the Patience genome (see Fig. S3 in the supplemental material).

Patience is the first phage to our knowledge for which the proteomic profile in infected cells has been examined by mass spectrometry. The approach is highly informative, providing strong evidence that many of the annotated reading frames are expressed-including 29 that are shorter than 120 codons-and providing support for many of the translational start sites, as well as revisions of start sites of several genes. Moreover, a previously unannotated gene was identified (gene 111) containing only 37 codons, and the revision of start codons indicates that two pairs of genes have significant overlaps (Ͼ60 bp). Phage genome annotation is generally more error-prone than that of other genomes because of the abundance of small open reading frames and relatively small gene size (average mycobacteriophage gene length is 640 bp). HPLC-MS/MS adds confidence to genome annotation, especially for phages such as Patience, whose coding potential does not closely match its bacterial host. It is surprising to find that at least two reading frames embedded out of frame within annotated genes are also expressed. These ORFs are generally not conserved and may not express functional products, but an intriguing possibility is that this expression is a consequence of movement into higher-GC hosts, presenting a small reservoir of new products available for selection and further adaption at little evolutionary cost. We note that ribosomal profiling of phage suggests that previously unannotated genes are expressed from its genome (26) .

The recent evolutionary history of Patience supports a model in which the diversity of viruses of a given host is a function of rapid movement from one host to another coupled with a landscape of diverse but closely related hosts in which the viruses evolve (4) . The large collection of mycobacteriophages that infect the common host M. smegmatis mc 2 155 encompasses considerable diversity of sequence and GC content, and it seems likely that other groups of phages have entered the high-GC environment relatively recently, with Patience representing an extreme example. We note that there are bacterial strains within the order Actinomycetales, such as those of Corynebacterium pseudotuberculosis, Corynebacterium ulcerans, and Corynebacterium diphtheriae, which have moderate GC contents (52.5%, 53.4%, and 53.5%, respectively) and may be relatives of hosts that supported Patience growth in its earlier evolutionary history. Their codon usage profiles (see Fig. S6 in the supplemental material) more closely reflect those of the lower-GC mycobacteriophages such as Patience. Patience may share gene content with phages of such hosts, but to date, few have been characterized, other than some prophages such as phage Beta of C. diphtheriae (27, 28) . Deeper exploration into the phages of these hosts and other hosts is thus likely to be highly informative and provide insights into viral origins and viral evolution (29) .

DNA sequencing. Patience was isolated using standard methods as described previously (7, 30) , and the genome was sequenced using 454 technology at the University of Pittsburgh's Genomics and Proteomics Core Laboratories; a total of~45,000 reads were assembled to yield an average (34) , and Phamerator (35) . The Phamerator database used for genomic comparisons was Mycobacteriophage_285. Phams were built using BLASTP and/or ClustalW, with similarity cutoff E values of 10 Ϫ50 and 32.5% similarity or better as described elsewhere (35) . Codon usage and codon selection. Codon usage resulting from mutational biases was estimated using the collection of all genes from a genome. The fraction (f) of each codon for a given amino acid was calculated as the ratio of the codon count to the amino acid count. Relative synonymous codon usage (RSCU) normalizes codon frequencies so that the sum of RSCU for codons of each amino acid is equal to the number of synonymous codons for that amino acid.

To measure codon selection, a second codon usage table (f o ) was tabulated, limited to genes experiencing strong codon selection. For bacterial genomes, this table was constructed from homologues of Sharp's set of 40 genes whose products participate in translation (36) . For Patience, this table was constructed in three steps. First, a codon usage table was generated from 20% of the genome using the genes with the most extreme value of codon usage bias as determined by 2 (where expected codon usage is calculated from the nucleotide composition). Next, adaptive codon enrichment (ACE u ) values (25) are calculated for all genes as described previously (25) , using this table to represent codon frequencies under codon selection (f o ) and the frequencies of codons among all genes in the Patience genome to represent codon frequencies expected from mutational processes alone (f N (34) , and Phamerator (35) .

Electron microscopy. CsCl gradient-purified Patience particles were applied to glow-discharged Formvar-and carbon-coated copper grids (400 mesh) (Ted Pella). They were stained with 1% uranyl acetate and imaged with a Morgagni 268 transmission electron microscope fitted with a Hamamatsu Orca HR side-model digital camera and AMT540 software. SDS-PAGE. Patience particles were concentrated and purified via CsCl gradient and ultracentrifugation. The visible phage band was dialyzed against two changes of phage buffer; 500 l of the dialyzed CsCl band was pelleted by a 30-min spin at 14,000 rpm in a microcentrifuge. The pellet was resuspended in 75 l of 20 mM dithiothreitol (DTT), and then 2 l of 0.5 M EDTA and 1 l of 1 M MgSO 4 were added. The phage was disrupted by being heated to 75°C for 2 min and then sonicated on ice six times for 30 s to disrupt the DNA. The sample was then mixed with 25 l of 4ϫ SDS sample buffer and heated in a boiling bath for 3 min at 95°C. The sample was electrophoresed through a 12% polyacrylamide gel containing SDS and stained with Coomassie brilliant blue in methanol.

HPLC-MS/MS. Five milliliters of exponentially growing M. smegmatis mc 2 155 (optical density at 600 nm [OD 600 ] of 0.4) in 7H9-ADC medium (30) was concentrated to a 500-l volume via low-speed centrifugation and infected with Patience at a multiplicity of infection (MOI) of 100. Phage particles were allowed to adsorb for 15 min, and then 4.5 ml of fresh 7H9 medium was added to the culture and incubated with shaking for 3 h at 37°C; the OD 600 was monitored throughout to follow cell growth and lysis. At 30 min and 150 min postadsorption, a 1-ml aliquot was removed from the culture, the cells were pelleted via centrifugation (1 min, 14,000 rpm in a microcentrifuge), and the supernatant was removed. The cell pellet was frozen at Ϫ80°C and then shipped overnight on wet ice to the University of California, Davis Proteomics Core (UCDPC) (http://proteomics.ucdavis.edu). There, the cells were lysed via a MagNA Lyser, the insoluble fraction was removed, and the soluble proteins were precipitated, digested with trypsin, and cleaned up using a MacroSpin column. The peptides were then separated using an Easy-LC II high-pressure liquid chromatography (HPLC) system and loaded into a Q Exactive Orbitrap mass spectrometer with a Proxeon nanospray source (Thermo) for tandem MS analysis. Detected spectra and fragmentation profiles were matched against a database comprised of a six-frame translation of the Patience genome, the annotated proteins of M. smegmatis mc 2 155, and UniProt using X! Tandem. Peptide matches were analyzed using Scaf-fold4. The "Relaxed" settings (as reported in Table S1 in the supplemental material) used a peptide false discovery rate (FDR) of 1% and a protein FDR of 5%; the "Stringent" settings used a peptide FDR of 0.1% and a protein FDR of 0.6%. Estimation of relative protein abundance was determined by normalizing the total number of spectra detected to the gene size.

Supplemental material for this article may be found at http://mbio.asm.org /lookup/suppl/doi:10.1128/mBio.02145-14/-/DCSupplemental.

Text S1, PDF file, 0.1 MB. Figure S1 , PDF file, 0.1 MB. Figure S2 , PDF file, 0.1 MB. Figure S3 , PDF file, 0.3 MB. Figure S4 , PDF file, 0.3 MB. Figure S5 , PDF file, 0.1 MB. Figure S6 , PDF file, 0.1 MB. Table S1 , PDF file, 0.1 MB.

The secret lives of mycobacteriophages

Science Education Alliance Phage Hunters Advancing Genomics and Evolutionary Science

Phage Hunters Integrating Research and Education (PHIRE) Program. 2013. Complete genome sequences of 63 mycobacteriophages

KwaZulu-Natal Research Institute for Tuberculosis and HIV Mycobacterial Genetics Course Students, Phage Hunters Integrating Research and Education Program

On the nature of mycobacteriophage diversity and host preference

Exploring the mycobacteriophage metaproteome: phage genomics as an educational platform

Deciphering the biology of Mycobacterium tuberculosis from the complete genome sequence

Origins of highly mosaic mycobacteriophage genomes

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

Functional role of mycobacteriophage transfer RNAs

Synonymous codon usage analysis of the mycobacteriophage Bxz1 and its plating bacteria M. smegmatis: identification of highly and lowly expressed genes of Bxz1 and the possible function of its tRNA species

Synonymous codon usage analysis of thirty two mycobacteriophage genomes

Anticodon nucleases

Mycobacteriophages: genes and genomes

Comparative phage genomics and the evolution of Siphoviridae: insights from dairy phages

HNH proteins are a widespread component of phage DNA packaging machines

Mycobacteriophage lysin B is a novel mycolylarabinogalactan esterase

The lytic cassette of mycobacteriophage Ms6 encodes an enzyme with lipolytic activity

The mycobacteriophage Ms6 encodes a chaperone-like protein involved in the endolysin delivery to the peptidoglycan

Mycobacteriophage endolysins: diverse and modular enzymes with multiple catalytic activities

Evidence for the intense exchange of MazG in marine cyanophages by horizontal gene transfer

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

Transcriptional silencing by the mycobacteriophage L5 repressor

Transcriptional regulation of repressor synthesis in mycobacteriophage L5

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

Quantification of codon selection for comparative bacterial genomics

High-resolution view of bacteriophage lambda gene expression by ribosome profiling

Characterization of bacteriophages from tox-containing, non-toxigenic isolates of Corynebacterium diphtheriae

A beta-related corynebacteriophage which lacks a tox allele but can acquire it by recombination with converting phage

Improved bacteriophage genome data is necessary for integrating viral and bacterial ecology

Gepard: a rapid and sensitive tool for creating dotplots on genome scale

ARAGORN, a program to detect tRNA genes and tmRNA genes in nucleotide sequences

tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence

The HHpred interactive server for protein homology detection and structure prediction

Phamerator: a bioinformatic tool for comparative bacteriophage genomics

Variation in the strength of selected codon usage bias among bacteria

DNA sequence, structure and gene expression of mycobacteriophage L5: a phage system for mycobacterial genetics