Skip to main content
Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
Microbiol Resour Announc. 2023 Nov; 12(11): e00852-23.
Published online 2023 Oct 25. doi: 10.1128/MRA.00852-23
PMCID: PMC10652884
PMID: 37877714

Isolation and genome sequence annotation of KillerTomato, a newly discovered cluster EE bacteriophage infecting Microbacterium

Omer Almahie, Data curation, Formal analysis, Investigation, Visualization, Writing – original draft, 1 Breanna Danklefsen, Data curation, Formal analysis, Investigation, Validation, Writing – original draft, 1 R. Scott Graves, Data curation, Investigation, Validation, Writing – original draft, 1 Sariah M. Hepworth, Data curation, Investigation, Validation, Writing – original draft, 1 Rylee T. Mathison, Data curation, Investigation, Validation, Writing – original draft, 1 Emily A. Sullivan, Data curation, Investigation, Validation, Visualization, Writing – original draft, 1 Paige R. Wood, Data curation, Investigation, Validation, Writing – original draft, 1 Skylar Bartholomew, Data curation, Investigation, Validation, Writing – original draft, 1 Reganne Brewster, Investigation, Validation, Visualization, Writing – original draft, 1 Dylan J. Hunt, Data curation, Investigation, Validation, 1 Carlos Serna, Data curation, Investigation, Validation, 1 Jack F. Shurley, Methodology, Supervision, Writing – review and editing, 1 Anna S. Grinath, Supervision, Validation, Writing – review and editing, 1 and Michael A. Thomas, Conceptualization, Data curation, Formal analysis, Funding acquisition, Methodology, Project administration, Supervision, Validation, Writing – review and editingcorresponding author 1
Simon Roux, Editor
Simon Roux, DOE Joint Genome Institute, Berkeley, California, USA ;

ABSTRACT

Discovered in Pocatello, Idaho, soil near a tomato garden, siphovirus KillerTomato infects Microbacterium foliorum NRRL B-24224. KillerTomato is a lytic cluster EE phage with a 17,442-bp genome and 68.6% GC content. Of 25 genes, 20 were assigned putative functions, including a putative tail assembly chaperone protein with a programmed frameshift and an endolysin.

KEYWORDS: bacteriophage, student-directed research, SEA-PHAGES

ANNOUNCEMENT

The discovery and annotation of bacteriophages are an important step in the development of phage therapy strategies for fighting bacterial infections (1). Here, we describe the characteristics of a lytic bacteriophage that infects Microbacterium foliorum.

KillerTomato was found in garden soil on 1 September 2022 on the Idaho State University campus in Pocatello, Idaho (42.867222°N, 112.428889° W), following established procedures (2, 3). Soil samples were washed with PYCa liquid media and filtered with a 0.22-μm filter. The filtrate was inoculated with Microbacterium foliorum NRRL B-24224. After incubation with shaking at 20°C for 48 h, the culture was filtered and the filtrate was plated on PYCa agar with M. foliorum and incubated at 20°C for 48 h. Plaques formed by KillerTomato were purified through three additional rounds of plating (Fig. 1a). Transmission electron micrographs of KillerTomato revealed a siphovirus with a long tail (Fig. 1b).

Using a Promega Wizard DNA clean up kit, KillerTomato DNA was extracted from a high-titer lysate and prepared for sequencing using a NEBNext Ultra II FS kit. DNA was sequenced using a shotgun sequencing approach on an Illumina MiSeq (v3 reagents). Sequencing generated 650,649 single-end 150-bp reads with 5,321-fold coverage. The raw reads were trimmed and assembled with Newbler v2.9 using default parameters, yielding a single contig; Consend v29 was used to check for completeness and accuracy and determine phage termini, both using default settings (4 6). KillerTomato has a genome length of 17,442 bp and a 3′ single-stranded sticky overhang of 9 bases (5′-CCCGCCCCA-3′). It has a GC content of 68.6%, which matches its host (68.7%) (7). KillerTomato belongs to cluster EE based on gene content and similarity to other EE phages in the Actinobacteriophage database (8).

The genome of KillerTomato was auto-annotated with PECAAN (http://discover.kbrinsgd.org) and DNAmaster v5.23.6 (http://cobamide2.bio.pitt.edu) with manual inspection and confirmation. Genemark v2.5 (9) and Glimmer v3.0 (10) were used to analyze the start sites and coding potential of each gene and for initial auto-annotation. We used Phamerator v454 (11) for cluster assignment and synteny analysis. Starterator (http://phages.wustl.edu/starterator/) was applied to compare start sites manually annotated for other EE phages to KillerTomato. HHPred v2.0 (12) (searches of pdb, pfamA, and NCBI CD databases) and BLAST v2.13.0 (13) (searches of the PhagesDB and NCBI nr databases) were used to predict putative gene functions. No tRNAs were identified using Aragorn v1.2.41 (14) and tRNA-SE v2.0 (15). Using default parameters, the annotation process yielded 25 protein-coding genes, 20 of which could be assigned functions (Fig. 1C).

An external file that holds a picture, illustration, etc.
Object name is mra.00852-23.f001.jpg

Characterization of bacteriophage KillerTomato. (a) KillerTomato forms round, medium-sized plaques. (b) Siphovirus morphology of KillerTomato with a capsid diameter of approximately 44 nm and a tail length of approximately 108 nm. Image produced by a Zeiss EM900 TEM with an accelerating voltage of 80 kV and uranyl acetate negative staining; 40–50 virions were measured. (c) The genome of KillerTomato contains 25 putative protein-coding genes (portrayed by the boxes with numbers depicting the gene number) on both forward and reverse strands.

Consistent with other EE phages, KillerTomato has a small genome with several fused genes (16). All but three genes (gp20 to gp22) are transcribed on one strand. KillerTomato has a putative tail assembly chaperone gene with a programmed frameshift (gp10 and gp11). KillerTomato also appears to encode an endolysin (gp17), a key protein involved in hydrolyzing the peptidoglycan layer of its host during a productive infection and that offers a promising strategy for treating infections by antibiotic-resistant bacteria (17). Along with plaque morphology, comparison with other cluster EE phages, which lack immunity repressor and integrase functions, suggests KillerTomato is a lytic phage.

ACKNOWLEDGMENTS

This study was supported by the Idaho State University Department of Biological Sciences and the Howard Hughes Medical Institute Science Education Alliance Phage Hunters Advancing Genomics and Evolutionary Science (HHMI SEA-PHAGES) program.

The soil sample was collected by R.B.

We thank Billy Biederman, Graham Hatfull, Deborah Jacobs-Sera, Daniel Russell, and Vic Sivanathan for training and resources and Rebecca Garlena and Daniel Russell for phage DNA sequencing. Special thanks are due to Jason Werth and Lisa McDougall from ISU’s Molecular Research Core Facility (MRCF) for the assistance with transmission electron microscopy imaging.

DATA AVAILABILITY

GenBank and SRA accession numbers for KillerTomato are OR159653 and SRX19690840, respectively.

REFERENCES

1. Housby JN, Mann NH. 2009. Phage therapy. Drug Discovery Today 14:536–540. doi: 10.1016/j.drudis.2009.03.006 [PubMed] [CrossRef] [Google Scholar]
2. Poxleitner M, Pope W, Jacobs-Sera D, Sivanathan V, Hatfull GF.. 2018. Phage discovery guide. Howard Hughes Medical Institute, Chevy Chase, MD, USA. https://seaphagesphagediscoveryguide.helpdocsonline.com. [Google Scholar]
3. Adams BM, Adams JB, Brewster RL, Cutler MS, Davis AE, Gallegos AH, Hernandez JS, May LH, Montoya EG, Reagan AT, Shurley JF, Grinath AS, Thomas MA. 2022. Annotation of the complete genome sequences of bacteriophages Ssara and birdfeeder. Microbiol Resour Announc 11:e0078022. doi: 10.1128/mra.00780-22 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
4. Russell DA. 2018. Sequencing, assembling, and finishing complete bacteriophage genomes. Methods Mol Biol 1681:109–125. doi: 10.1007/978-1-4939-7343-9_9 [PubMed] [CrossRef] [Google Scholar]
5. Gordon D, Green P. 2013. Consed: a graphical editor for next-generation sequencing. Bioinformatics 29:2936–2937. doi: 10.1093/bioinformatics/btt515 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
6. Chung C-H, Walter MH, Yang L, Chen S-C, Winston V, Thomas MA. 2017. Predicting genome terminus sequences of Bacillus cereus-group bacteriophage using next generation sequencing data. BMC Genomics 18:350. doi: 10.1186/s12864-017-3744-0 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
7. Russell DA, Garlena RA, Hatfull GF, Maresca JA. 2019. Complete genome sequence of Microbacterium foliorum NRRL B-24224, a host for bacteriophage discovery. Microbiol Resour Announc 8:e01467-18. doi: 10.1128/MRA.01467-18 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
8. Russell DA, Hatfull GF. 2017. PhagesDB: the actinobacteriophage database. Bioinformatics 33:784–786. doi: 10.1093/bioinformatics/btw711 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
9. Besemer J, Borodovsky M. 2005. GeneMark: web software for gene finding in prokaryotes, eukaryotes and viruses. Nucleic Acids Res 33:W451–4. doi: 10.1093/nar/gki487 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
10. Delcher AL, Bratke KA, Powers EC, Salzberg SL. 2007. Identifying bacterial genes and endosymbiont DNA with glimmer. Bioinformatics 23:673–679. doi: 10.1093/bioinformatics/btm009 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
11. Cresawn SG, Bogel M, Day N, Jacobs-Sera D, Hendrix RW, Hatfull GF. 2011. Phamerator: a bioinformatic tool for comparative bacteriophage genomics. BMC Bioinformatics 12:395. doi: 10.1186/1471-2105-12-395 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
12. Söding J, Biegert A, Lupas AN. 2005. The HHpred interactive server for protein homology detection and structure prediction. Nucleic Acids Res 33:W244–8. doi: 10.1093/nar/gki408 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
13. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. 1990. Basic local alignment search tool. J Mol Biol 215:403–410. doi: 10.1016/S0022-2836(05)80360-2 [PubMed] [CrossRef] [Google Scholar]
14. Laslett D, Canback B. 2004. ARAGORN, a program to detect tRNA genes and tmRNA genes in nucleotide sequences. Nucleic Acids Research 32:11–16. doi: 10.1093/nar/gkh152 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
15. Lowe TM, Eddy SR. 1997. tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res 25:955–964. doi: 10.1093/nar/25.5.955 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
16. Jacobs-Sera D, Abad LA, Alvey RM, Anders KR, Aull HG, Bhalla SS, Blumer LS, Bollivar DW, Bonilla JA, Butela KA, Coomans RJ, Cresawn SG, D’Elia T, Diaz A, Divens AM, Edgington NP, Frederick GD, Gainey MD, Garlena RA, Grant KW, Gurney SMR, Hendrickson HL, Hughes LE, Kenna MA, Klyczek KK, Kotturi H, Mavrich TN, McKinney AL, Merkhofer EC, Moberg Parker J, Molloy SD, Monti DL, Pape-Zambito DA, Pollenz RS, Pope WH, Reyna NS, Rinehart CA, Russell DA, Shaffer CD, Sivanathan V, Stoner TH, Stukey J, Sunnen CN, Tolsma SS, Tsourkas PK, Wallen JR, Ware VC, Warner MH, Washington JM, Westover KM, Whitefleet-Smith JL, Wiersma-Koch HI, Williams DC, Zack KM, Hatfull GF. 2020. Genomic diversity of bacteriophages infecting Microbacterium spp. PLoS One 15:e0234636. doi: 10.1371/journal.pone.0234636 [PMC free article] [PubMed] [CrossRef] [Google Scholar]
17. Rahman MU, Wang W, Sun Q, Shah JA, Li C, Sun Y, Li Y, Zhang B, Chen W, Wang S. 2021. Endolysin, a promising solution against antimicrobial resistance. Antibiotics (Basel) 10:11. doi: 10.3390/antibiotics10111277 [PMC free article] [PubMed] [CrossRef] [Google Scholar]

Articles from Microbiology Resource Announcements are provided here courtesy of American Society for Microbiology (ASM)

-