Welkin H. Pope

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.

A Broadly Implementable Research Course in Phage Discovery and Genomics for First-Year Undergraduate Students

ABSTRACT 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. 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.

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

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

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.

Prophage-mediated defence against viral attack and viral counter-defence

Temperate phages are common, and prophages are abundant residents of sequenced bacterial genomes. Mycobacteriophages are viruses that infect mycobacterial hosts including Mycobacterium tuberculosis and Mycobacterium smegmatis, encompass substantial genetic diversity and are commonly temperate. Characterization of ten Cluster N temperate mycobacteriophages revealed at least five distinct prophage-expressed viral defence systems that interfere with the infection of lytic and temperate phages that are either closely related (homotypic defence) or unrelated (heterotypic defence) to the prophage. Target specificity is unpredictable, ranging from a single target phage to one-third of those tested. The defence systems include a single-subunit restriction system, a heterotypic exclusion system and a predicted (p)ppGpp synthetase, which blocks lytic phage growth, promotes bacterial survival and enables efficient lysogeny. The predicted (p)ppGpp synthetase coded by the Phrann prophage defends against phage Tweety infection, but Tweety codes for a tetrapeptide repeat protein, gp54, which acts as a highly effective counter-defence system. Prophage-mediated viral defence offers an efficient mechanism for bacterial success in host–virus dynamics, and counter-defence promotes phage co-evolution.

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

ABSTRACT Genomic analysis of a large set of phages infecting the common host Mycobacterium smegmatis mc2155 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.

Genomics and Proteomics of Mycobacteriophage Patience, an Accidental Tourist in the Mycobacterium Neighborhood

ABSTRACT 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 mc2155 as a host. Patience has genomic features distinct from its M. smegmatis host, including a much lower GC content (50.3% versus 67.4%) and an abundance of codons that are rarely used in M. smegmatis. Nonetheless, it propagates well in M. smegmatis, and we demonstrate the use of mass spectrometry to show expression of over 75% of the predicted proteins, to identify new genes, to refine the genome annotation, and to estimate protein abundance. We propose that Patience evolved primarily among lower-GC hosts and that the disparities between its genomic profile and that of M. smegmatis presented only a minimal barrier to host expansion. Rapid adaptions to its new host include recent acquisition of higher-GC genes, expression of out-of-frame proteins within predicted genes, and codon selection among highly expressed genes toward the translational apparatus of its new host. IMPORTANCE The mycobacteriophage Patience genome has a notably lower GC content (50.3%) than its Mycobacterium smegmatis host (67.4%) and has markedly different codon usage biases. The viral genome has an abundance of codons that are rare in the host and are decoded by wobble tRNA pairing, although the phage grows well and expression of most of the genes is detected by mass spectrometry. Patience thus has the genomic profile of a virus that evolved primarily in one type of host genetic landscape (moderate-GC bacteria) but has found its way into a distinctly different high-GC environment. Although Patience genes are ill matched to the host expression apparatus, this is of little functional consequence and has not evidently imposed a barrier to migration across the microbial landscape. Interestingly, comparison of expression levels and codon usage profiles reveals evidence of codon selection as the genome evolves and adapts to its new environment. The mycobacteriophage Patience genome has a notably lower GC content (50.3%) than its Mycobacterium smegmatis host (67.4%) and has markedly different codon usage biases. The viral genome has an abundance of codons that are rare in the host and are decoded by wobble tRNA pairing, although the phage grows well and expression of most of the genes is detected by mass spectrometry. Patience thus has the genomic profile of a virus that evolved primarily in one type of host genetic landscape (moderate-GC bacteria) but has found its way into a distinctly different high-GC environment. Although Patience genes are ill matched to the host expression apparatus, this is of little functional consequence and has not evidently imposed a barrier to migration across the microbial landscape. Interestingly, comparison of expression levels and codon usage profiles reveals evidence of codon selection as the genome evolves and adapts to its new environment.

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

ABSTRACT 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 mc2155. 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.

Protein Folding Failure Sets High-Temperature Limit on Growth of Phage P22 in Salmonella enterica Serovar Typhimurium

ABSTRACT The high-temperature limit for growth of microorganisms differs greatly depending on their species and habitat. The importance of an organisms ability to manage thermal stress is reflected in the ubiquitous distribution of the heat shock chaperones. Although many chaperones function to reduce protein folding defects, it has been difficult to identify the specific protein folding pathways that set the high-temperature limit of growth for a given microorganism. We have investigated this for a simple system, phage P22 infection of Salmonella enterica serovar Typhimurium. Production of infectious particles exhibited a broad maximum of 150 phage per cell when host cells were grown at between 30 and 39°C in minimal medium. Production of infectious phage declined sharply in the range of 40 to 41°C, and at 42°C, production had fallen to less than 1% of the maximum rate. The host cells maintained optimal division rates at these temperatures. The decrease in phage infectivity was steeper than the loss of physical particles, suggesting that noninfectious particles were formed at higher temperatures. Sodium dodecyl sulfate-polyacrylamide gel electrophoresis revealed a decrease in the tailspike adhesins assembled on phage particles purified from cultures incubated at higher temperatures. The infectivity of these particles was restored by in vitro incubation with soluble tailspike trimers. Examination of tailspike folding and assembly in lysates of phage-infected cells confirmed that the fraction of polypeptide chains able to reach the native state in vivo decreased with increasing temperature, indicating a thermal folding defect rather than a particle assembly defect. Thus, we believe that the folding pathway of the tailspike adhesin sets the high-temperature limit for P22 formation in Salmonella serovar Typhimurium.

Cluster J Mycobacteriophages: Intron Splicing in Capsid and Tail Genes

Bacteriophages isolated on Mycobacterium smegmatis mc2155 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.