Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2020 Jun 18;15(6):e0234636.
doi: 10.1371/journal.pone.0234636. eCollection 2020.

Genomic diversity of bacteriophages infecting Microbacterium spp

Deborah Jacobs-Sera  1 , Lawrence A Abad  1 , Richard M Alvey  2 , Kirk R Anders  3 , Haley G Aull  1 , Suparna S Bhalla  4 , Lawrence S Blumer  5 , David W Bollivar  2 , J Alfred Bonilla  6 , Kristen A Butela  1 , Roy J Coomans  7 , Steven G Cresawn  8 , Tom D'Elia  9 , Arturo Diaz  10 , Ashley M Divens  1 , Nicholas P Edgington  11 , Gregory D Frederick  12 , Maria D Gainey  13 , Rebecca A Garlena  1 , Kenneth W Grant  14 , Susan M R Gurney  15 , Heather L Hendrickson  16 , Lee E Hughes  17 , Margaret A Kenna  18 , Karen K Klyczek  6 , Hari Kotturi  19 , Travis N Mavrich  1 , Angela L McKinney  20 , Evan C Merkhofer  4 , Jordan Moberg Parker  21 , Sally D Molloy  22 , Denise L Monti  23 , Dana A Pape-Zambito  24 , Richard S Pollenz  25 , Welkin H Pope  1 , Nathan S Reyna  26 , Claire A Rinehart  27 , Daniel A Russell  1 , Christopher D Shaffer  28 , Viknesh Sivanathan  29 , Ty H Stoner  1 , Joseph Stukey  30 , C Nicole Sunnen  31 , Sara S Tolsma  32 , Philippos K Tsourkas  33 , Jamie R Wallen  13 , Vassie C Ware  18 , Marcie H Warner  1 , Jacqueline M Washington  34 , Kristi M Westover  35 , JoAnn L Whitefleet-Smith  36 , Helen I Wiersma-Koch  9 , Daniel C Williams  37 , Kira M Zack  1 , Graham F Hatfull  1
Affiliations

Genomic diversity of bacteriophages infecting Microbacterium spp

Deborah Jacobs-Sera et al. PLoS One. .

Abstract

The bacteriophage population is vast, dynamic, old, and genetically diverse. The genomics of phages that infect bacterial hosts in the phylum Actinobacteria show them to not only be diverse but also pervasively mosaic, and replete with genes of unknown function. To further explore this broad group of bacteriophages, we describe here the isolation and genomic characterization of 116 phages that infect Microbacterium spp. Most of the phages are lytic, and can be grouped into twelve clusters according to their overall relatedness; seven of the phages are singletons with no close relatives. Genome sizes vary from 17.3 kbp to 97.7 kbp, and their G+C% content ranges from 51.4% to 71.4%, compared to ~67% for their Microbacterium hosts. The phages were isolated on five different Microbacterium species, but typically do not efficiently infect strains beyond the one on which they were isolated. These Microbacterium phages contain many novel features, including very large viral genes (13.5 kbp) and unusual fusions of structural proteins, including a fusion of VIP2 toxin and a MuF-like protein into a single gene. These phages and their genetic components such as integration systems, recombineering tools, and phage-mediated delivery systems, will be useful resources for advancing Microbacterium genetics.

PubMed Disclaimer

Conflict of interest statement

The authors have declared that no competing interests exist.

Figures

Fig 1
Fig 1. Geographical distribution of Microbacterium phages.
Geographic distribution of the isolation sites of Microbacterium phages with completely sequenced genomes.
Fig 2
Fig 2. Microbacterium phage virion morphologies.
Representative virion particles of each Microbacterium phage cluster show a predominance of siphoviridae morphologies; Cluster EK and singleton Burro are podoviridae. All have isometric capsids with the exception of Count, which has a prolate capsid.
Fig 3
Fig 3. Heatmap of Microbacterium phage average nucleotide identities (ANIs).
Pairwise average nucleotide identities were calculated for 116 Microbacterium phages using DNA Master with default settings. The heatmap was generated using R and the ‘heatmap2’ function, which determines distances between each genome and calculates the optimal genome order for representation, using distance parameter and clustering methods of ‘maximum’ and ‘single’, respectively. Genome clusters are shown on the axes, colored according to cluster, with singletons show in black; subclusters are indicated with numbers alongside their cluster designations. Phage vB_MosX-ISF9 genome is not included in these analyses [30].
Fig 4
Fig 4. Relationships among Microbacterium phages and their bacterial hosts.
A. A network phylogeny of Microbacterium phages. The predicted proteins of all 116 Microbacterium phage were sorted into 1975 phamilies (phams) according to shared amino acid sequence similarities using Phamerator [34]. Each genome was then assigned a value reflecting the presence or absence of a pham member; the genomes were compared and displayed using Splitstree [35]. The clusters and subclusters derived from dotplot and ANI comparisons are indicated with larger colored circles. Smaller colored circles at the nodes indicate the bacterium host used for isolation, as noted in the key. Phages isolated on M. foliorum NRRL B-24224 have no small colored circle at the node. The scale bar indicates 0.01 substitutions/site. B. Evolutionary relationships of bacterial host taxa, using MEGA7 [–38]. The number of phages isolated on each host is shown in parentheses.
Fig 5
Fig 5. Pairwise alignment of Microbacterium phage Cluster EA genomes.
A representative genome from each EA subcluster phage is shown. The sole Subcluster EA7 phage Theresita is shown at the bottom as it is substantially different from the others in the right halves of the genomes. Pairwise nucleotide sequence similarities are displayed with spectrum-coloring between genomes, with violet representing greatest similarity and red the least similar, above a threshold E value of 10−5. Genes are represented as boxes above or below the genomes reflecting rightwards- and leftwards-transcription respectively. Genes are colored according to their phamily designations using Phamerator [34] and database Actinobacteriophage_2422. White boxes represent ‘orphams’, genes with no close relatives in this dataset.
Fig 6
Fig 6. Genome organization of Microbacterium Subcluster EA1 phage TeddyBear.
A. The genome of Microbacterium phage TeddyBear is shown with predicted genes represented as boxes above or below the genome reflecting rightwards- and leftwards-transcription respectively. Genes are colored according to their phamily designations using Phamerator [34] and database Actinobacteriophage_2422. The phamily numbers shown above each gene with the number of phamily members in parentheses. B. Pairwise alignment of Microbacterium phages Golden, Theresita, Goodman, and Johann genomes. See Fig 5 for map details. Theresita shares 45.2% and 23.6% average gene content with Golden and Goodman, respectively.
Fig 7
Fig 7. Genome organization of Microbacterium Cluster EB phage Dismas.
See Fig 6A for details.
Fig 8
Fig 8. Genome organization of Microbacterium cluster EC phage Quahwah.
See Fig 6A for details. Vertical arrows show the positions of 18 bp repeated motifs corresponding to the consensus 5’-TAGaCTCTaGGTgTaAgC, where upper case letters indicate complete conservation, and lower case letters are those appearing in at least 9 of the 12 instances. Each instance of the motif has no more than eight departures from the 18 bp consensus sequence and they are numbered as listed in S10A Fig.
Fig 9
Fig 9. Genome organization of Microbacterium Subcluster ED1 phage Hortus1.
See Fig 6A for details.
Fig 10
Fig 10. Genome organization of Microbacterium Cluster EE phage BurtonThePup.
See Fig 6A for details.
Fig 11
Fig 11. Genome organization of Microbacterium cluster EF phage AnnaSerena, See Fig 6A for details.
The small red vertical arrows indicate the locations containing both conserved SAS and ESAS sequences; black arrows indicate the locations with an ESAS motif that lack an associated SAS. The SAS motifs are numbered as listed in S10B and S10C Fig.
Fig 12
Fig 12. Genome organization of Microbacterium cluster EG phage hyperion.
See Fig 6A for details. Vertical arrows indicate the positions of conserved short inverted repeat sequences and numbered as listed in S10D Fig.
Fig 13
Fig 13. Genome organization of Microbacterium cluster EH phage percival.
See Fig 6A for details.
Fig 14
Fig 14. Genome organization of Microbacterium cluster EI phage MementoMori.
See Fig 6A for details. Vertical arrows indicate the locations of 30 bp repeated sequences in intergenic regions, numbered as shown in S10E.
Fig 15
Fig 15. Genome organization of Microbacterium Cluster EJ phage Goodman.
See Fig 6A for details.
Fig 16
Fig 16. Microbacterium cluster EK phages and singleton burro.
A. Genome organization of phages ArMaWen, Akoni and Burro, clusters EK1, EK2, and a singleton, respectively. See Fig 5 for details. B. SDS-PAGE of Burro virions. Lanes 1 and 2 correspond to Burro virions purified through one and two CsCl density gradients, respectively. M, marker with protein sizes indicated in kDa.
Fig 17
Fig 17. Genome organization of Microbacterium Cluster EL phages, Camille and Count.
See Fig 5 for details.
Fig 18
Fig 18. Genome organization of Microbacterium Singleton phage Appa.
See Fig 6A for details.
Fig 19
Fig 19. Genome organization of Microbacterium Singleton phage Hendrix.
See Fig 6A for details.
Fig 20
Fig 20. Genome organization of Microbacterium Singleton phage Triscuit.
See Fig 6A for details.
Fig 21
Fig 21. Genome organization of Microbacterium Singleton phage ValentiniPuff.
See Fig 6A for details.
Fig 22
Fig 22. Genome organization of Microbacterium Singleton phage Zeta1847.
See Fig 6A for details.
Fig 23
Fig 23. Comparison of phage genomic diversity by host genera.
A. Genomic similarity plot comparing gene content dissimilarity (gcd) to whole genome nucleotide distance, as previously reported [49]. Each data point reflects a genome comparison involving one phage that infects the indicated host genus and (orange) a second phage that infects the indicated host genus or (black) one phage that infects a different host genus. n = number of phages that infect the indicated host genus. B. Box plot comparing the MaxGCDGap of phages from different host genera. Each data point is a phage genome, the box depicts the middle 50% of data, and the black bar represents the median. Number of phages per host genus as in panel A. C. Three representative genome networks highlight genomic relationships of Microbacterium phages to phages of other host genera. Each node represents a phage genome and is colored according to host genus. Edges between nodes represent pairs of phages that exhibit ‘intra-cluster’ genomic similarities (as measured by gene content dissimilarity and nucleotide distance from panel A). Selected phage names or cluster designations are highlighted for reference.
See this image and copyright information in PMC

Similar articles

  • Bacteriophages of Gordonia spp. Display a Spectrum of Diversity and Genetic Relationships.
    Pope WH, Mavrich TN, Garlena RA, Guerrero-Bustamante CA, Jacobs-Sera D, Montgomery MT, Russell DA, Warner MH; Science Education Alliance-Phage Hunters Advancing Genomics and Evolutionary Science (SEA-PHAGES); Hatfull GF. Pope WH, et al. mBio. 2017 Aug 15;8(4):e01069-17. doi: 10.1128/mBio.01069-17. mBio. 2017. PMID: 28811342 Free PMC article.
  • Actinobacteriophages: Genomics, Dynamics, and Applications.
    Hatfull GF. Hatfull GF. Annu Rev Virol. 2020 Sep 29;7(1):37-61. doi: 10.1146/annurev-virology-122019-070009. Annu Rev Virol. 2020. PMID: 32991269 Free PMC article. Review.
  • Tales of diversity: Genomic and morphological characteristics of forty-six Arthrobacter phages.
    Klyczek KK, Bonilla JA, Jacobs-Sera D, Adair TL, Afram P, Allen KG, Archambault ML, Aziz RM, Bagnasco FG, Ball SL, Barrett NA, Benjamin RC, Blasi CJ, Borst K, Braun MA, Broomell H, Brown CB, Brynell ZS, Bue AB, Burke SO, Casazza W, Cautela JA, Chen K, Chimalakonda NS, Chudoff D, Connor JA, Cross TS, Curtis KN, Dahlke JA, Deaton BM, Degroote SJ, DeNigris DM, DeRuff KC, Dolan M, Dunbar D, Egan MS, Evans DR, Fahnestock AK, Farooq A, Finn G, Fratus CR, Gaffney BL, Garlena RA, Garrigan KE, Gibbon BC, Goedde MA, Guerrero Bustamante CA, Harrison M, Hartwell MC, Heckman EL, Huang J, Hughes LE, Hyduchak KM, Jacob AE, Kaku M, Karstens AW, Kenna MA, Khetarpal S, King RA, Kobokovich AL, Kolev H, Konde SA, Kriese E, Lamey ME, Lantz CN, Lapin JS, Lawson TO, Lee IY, Lee SM, Lee-Soety JY, Lehmann EM, London SC, Lopez AJ, Lynch KC, Mageeney CM, Martynyuk T, Mathew KJ, Mavrich TN, McDaniel CM, McDonald H, McManus CJ, Medrano JE, Mele FE, Menninger JE, Miller SN, Minick JE, Nabua CT, Napoli CK, Nkangabwa M, Oates EA, Ott CT, Pellerino SK, Pinamont WJ, Pirnie RT, Pizzorno MC, Plautz EJ, Pope WH, Pruett KM, Rickstrew G, Rimple PA, Rinehart CA, Robinson KM, Rose VA, Russell DA, Schick AM, Schlossman J,… See abstract for full author list ➔ Klyczek KK, et al. PLoS One. 2017 Jul 17;12(7):e0180517. doi: 10.1371/journal.pone.0180517. eCollection 2017. PLoS One. 2017. PMID: 28715480 Free PMC article.
  • Complete genomic sequences of Propionibacterium freudenreichii phages from Swiss cheese reveal greater diversity than Cutibacterium (formerly Propionibacterium) acnes phages.
    Cheng L, Marinelli LJ, Grosset N, Fitz-Gibbon ST, Bowman CA, Dang BQ, Russell DA, Jacobs-Sera D, Shi B, Pellegrini M, Miller JF, Gautier M, Hatfull GF, Modlin RL. Cheng L, et al. BMC Microbiol. 2018 Mar 1;18(1):19. doi: 10.1186/s12866-018-1159-y. BMC Microbiol. 2018. PMID: 29490612 Free PMC article.
  • Bacteriophages and their genomes.
    Hatfull GF, Hendrix RW. Hatfull GF, et al. Curr Opin Virol. 2011 Oct;1(4):298-303. doi: 10.1016/j.coviro.2011.06.009. Curr Opin Virol. 2011. PMID: 22034588 Free PMC article. Review.
See all similar articles

Cited by

See all "Cited by" articles

References

    1. Hatfull GF. Mycobacteriophages. Microbiol Spectr. 2018;6(5). Epub 2018/10/07. 10.1128/microbiolspec.GPP3-0026-2018 . - DOI - PMC - PubMed
    1. Pope WH, Mavrich TN, Garlena RA, Guerrero-Bustamante CA, Jacobs-Sera D, Montgomery MT, et al. Bacteriophages of Gordonia spp. Display a Spectrum of Diversity and Genetic Relationships. MBio. 2017;8(4). Epub 2017/08/16. 10.1128/mBio.01069-17 - DOI - PMC - PubMed
    1. Klyczek KK, Bonilla JA, Jacobs-Sera D, Adair TL, Afram P, Allen KG, et al. Tales of diversity: Genomic and morphological characteristics of forty-six Arthrobacter phages. PLoS One. 2017;12(7):e0180517 Epub 2017/07/18. 10.1371/journal.pone.0180517 - DOI - PMC - PubMed
    1. Marinelli LJ, Fitz-Gibbon S, Hayes C, Bowman C, Inkeles M, Loncaric A, et al. Propionibacterium acnes Bacteriophages Display Limited Genetic Diversity and Broad Killing Activity against Bacterial Skin Isolates. MBio. 2012;3(5). Epub 2012/09/28. 10.1128/mBio.00279-12 . - DOI - PMC - PubMed
    1. Russell DA, Hatfull GF. PhagesDB: the actinobacteriophage database. Bioinformatics. 2017;33(5):784–6. 10.1093/bioinformatics/btw711 . - DOI - PMC - PubMed

Publication types

-