Deborah Jacobs-Sera

Exploring the Mycobacteriophage Metaproteome: Phage Genomics as an Educational Platform

Bacteriophages are the most abundant forms of life in the biosphere and carry genomes characterized by high genetic diversity and mosaic architectures. The complete sequences of 30 mycobacteriophage genomes show them collectively to encode 101 tRNAs, three tmRNAs, and 3,357 proteins belonging to 1,536 “phamilies” of related sequences, and a statistical analysis predicts that these represent approximately 50% of the total number of phamilies in the mycobacteriophage population. These phamilies contain 2.19 proteins on average; more than half (774) of them contain just a single protein sequence. Only six phamilies have representatives in more than half of the 30 genomes, and only three—encoding tape-measure proteins, lysins, and minor tail proteins—are present in all 30 phages, although these phamilies are themselves highly modular, such that no single amino acid sequence element is present in all 30 mycobacteriophage genomes. Of the 1,536 phamilies, only 230 (15%) have amino acid sequence similarity to previously reported proteins, reflecting the enormous genetic diversity of the entire phage population. The abundance and diversity of phages, the simplicity of phage isolation, and the relatively small size of phage genomes support bacteriophage isolation and comparative genomic analysis as a highly suitable platform for discovery-based education.

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.

Phamerator: a bioinformatic tool for comparative bacteriophage genomics

BackgroundBacteriophage genomes have mosaic architectures and are replete with small open reading frames of unknown function, presenting challenges in their annotation, comparative analysis, and representation.ResultsWe describe here a bioinformatic tool, Phamerator, that assorts protein-coding genes into phamilies of related sequences using pairwise comparisons to generate a database of gene relationships. This database is used to generate genome maps of multiple phages that incorporate nucleotide and amino acid sequence relationships, as well as genes containing conserved domains. Phamerator also generates phamily circle representations of gene phamilies, facilitating analysis of the different evolutionary histories of individual genes that migrate through phage populations by horizontal genetic exchange.ConclusionsPhamerator represents a useful tool for comparative genomic analysis and comparative representations of bacteriophage genomes.

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

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

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.

Comparative genomic analysis of mycobacteriophage Tweety: evolutionary insights and construction of compatible site-specific integration vectors for mycobacteria.

Mycobacteriophage Tweety is a newly isolated phage of Mycobacterium smegmatis. It has a viral morphology with an isometric head and a long flexible tail, and forms turbid plaques from which stable lysogens can be isolated. The Tweety genome is 58 692 bp in length, contains 109 protein-coding genes, and shows significant but interrupted nucleotide sequence similarity with the previously described mycobacteriophages Llij, PMC and Che8. However, overall the genome possesses mosaic architecture, with gene products being related to other mycobacteriophages such as Che9d, Omega and Corndog. A gene encoding an integrase of the tyrosine-recombinase family is located close to the centre of the genome, and a putative attP site has been identified within a short intergenic region immediately upstream of int. This Tweety attP–int cassette was used to construct a new set of integration-proficient plasmid vectors that efficiently transform both fast- and slow-growing mycobacteria through plasmid integration at a chromosomal locus containing a tRNALys gene. These vectors are maintained well in the absence of selection and are completely compatible with integration vectors derived from mycobacteriophage L5, enabling the simple construction of complex recombinants with genes integrated simultaneously at different chromosomal positions.

Structure of the archaeal head-tailed virus HSTV-1 completes the HK97 fold story

It has been proposed that viruses can be divided into a small number of structure-based viral lineages. One of these lineages is exemplified by bacterial virus Hong Kong 97 (HK97), which represents the head-tailed dsDNA bacteriophages. Seemingly similar viruses also infect archaea. Here we demonstrate using genomic analysis, electron cryomicroscopy, and image reconstruction that the major coat protein fold of newly isolated archaeal Haloarcula sinaiiensis tailed virus 1 has the canonical coat protein fold of HK97. Although it has been anticipated previously, this is physical evidence that bacterial and archaeal head-tailed viruses share a common architectural principle. The HK97-like fold has previously been recognized also in herpesviruses, and this study expands the HK97-like lineage to viruses from all three domains of life. This is only the second established lineage to include archaeal, bacterial, and eukaryotic viruses. Thus, our findings support the hypothesis that the last common universal ancestor of cellular organisms was infected by a number of different viruses.

Genomic characterization of mycobacteriophage giles: Evidence for phage acquisition of host DNA by illegitimate recombination

A characteristic feature of bacteriophage genomes is that they are architecturally mosaic, with each individual genome representing a unique assemblage of individual exchangeable modules. Plausible mechanisms for generating mosaicism include homologous recombination at shared boundary sequences of module junctions, illegitimate recombination in a non-sequence-directed process, and site-specific recombination. Analysis of the novel mycobacteriophage Giles genome not only extends our current perspective on bacteriophage genetic diversity, with more than 60% of the genes unrelated to other mycobacteriophages, but offers novel insights into how mosaic genomes are created. In one example, the integration/excision cassette is atypically situated within the structural gene operon and could have moved there either by illegitimate recombination or more plausibly via integrase-mediated site-specific recombination. In a second example, a DNA segment has been recently acquired from the host bacterial chromosome by illegitimate recombination, providing further evidence that phage genomic mosaicism is generated by nontargeted recombination processes.

Propionibacterium acnes Bacteriophages Display Limited Genetic Diversity and Broad Killing Activity against Bacterial Skin Isolates

ABSTRACT Investigation of the human microbiome has revealed diverse and complex microbial communities at distinct anatomic sites. The microbiome of the human sebaceous follicle provides a tractable model in which to study its dominant bacterial inhabitant, Propionibacterium acnes, which is thought to contribute to the pathogenesis of the human disease acne. To explore the diversity of the bacteriophages that infect P. acnes, 11 P. acnes phages were isolated from the sebaceous follicles of donors with healthy skin or acne and their genomes were sequenced. Comparative genomic analysis of the P. acnes phage population, which spans a 30-year temporal period and a broad geographic range, reveals striking similarity in terms of genome length, percent GC content, nucleotide identity (>85%), and gene content. This was unexpected, given the far-ranging diversity observed in virtually all other phage populations. Although the P. acnes phages display a broad host range against clinical isolates of P. acnes, two bacterial isolates were resistant to many of these phages. Moreover, the patterns of phage resistance correlate closely with the presence of clustered regularly interspaced short palindromic repeat elements in the bacteria that target a specific subset of phages, conferring a system of prokaryotic innate immunity. The limited diversity of the P. acnes bacteriophages, which may relate to the unique evolutionary constraints imposed by the lipid-rich anaerobic environment in which their bacterial hosts reside, points to the potential utility of phage-based antimicrobial therapy for acne. IMPORTANCE Propionibacterium acnes is a dominant member of the skin microflora and has also been implicated in the pathogenesis of acne; however, little is known about the bacteriophages that coexist with and infect this bacterium. Here we present the novel genome sequences of 11 P. acnes phages, thereby substantially increasing the amount of available genomic information about this phage population. Surprisingly, we find that, unlike other well-studied bacteriophages, P. acnes phages are highly homogeneous and show a striking lack of genetic diversity, which is perhaps related to their unique and restricted habitat. They also share a broad ability to kill clinical isolates of P. acnes; phage resistance is not prevalent, but when detected, it appears to be conferred by chromosomally encoded immunity elements within the host genome. We believe that these phages display numerous features that would make them ideal candidates for the development of a phage-based therapy for acne. Propionibacterium acnes is a dominant member of the skin microflora and has also been implicated in the pathogenesis of acne; however, little is known about the bacteriophages that coexist with and infect this bacterium. Here we present the novel genome sequences of 11 P. acnes phages, thereby substantially increasing the amount of available genomic information for this phage population. Surprisingly, we find that, unlike other well-studied bacteriophages, P. acnes phages are highly homogeneous and show a striking lack of genetic diversity, which is perhaps related to their unique and restricted habitat. They also share a broad ability to kill clinical isolates of P. acnes; phage resistance is not prevalent, but when detected, it appears to be conferred by chromosomally encoded immunity elements within the host genome. We believe that these phages display numerous features that would make them ideal candidates for the development of a phage-based therapy for acne.