Paenibacillus larvae and their phages; a community science approach to discovery and initial testing of prophylactic phage cocktails against American Foulbrood in New Zealand

Abstract

Background: American foulbrood (AFB) is a devastating disease of the European honey bee (Apis mellifera) and is found throughout the world. AFB is caused by the bacterium Paenibacillus larvae (P. larvae). Treatment with antibiotics is strictly forbidden in many regions, including New Zealand. Safe and natural prophylactic solutions to protect honey bees from AFB are needed. Bacteriophages are a well-studied alternative to antibiotics and have been shown to be effective against P. larvae in other countries.

Methods: We employed a community science approach to obtaining samples from around New Zealand to discover novel bacteriophages. Standard isolation approaches were employed for both bacteria and bacteriophages. Host range testing was performed by agar overlay spot tests, and cocktail formulation and in vitro testing were performed in 96-well plate assays, followed by sub-sampling and CFU visualization on agar plates.

Results: Herein, we describe the discovery and isolation of eight P. larvae bacterial isolates and 26 P. larvae bacteriophages that are novel and native to New Zealand. The phage genomes were sequenced and annotated, and their genomes were compared to extant sequenced P. larvae phage genomes. We test the host ranges of the bacteriophages and formulate cocktails to undertake in vitro testing on a set of representative bacterial strains. These results form the basis of a promising solution for protecting honey bees in New Zealand from AFB.

Keywords

Apis mellifera / phage / environmental sampling / community science

References

[1] Klein AM,Cane JH.Importance of pollinators in changing landscapes for world crops.Proc Biol Sci2007;274:303-13 PMCID:PMC1702377

[2] Newstrom-Lloyd LE. Pollination in New Zealand. In: Dymond JR, editor. Ecosystem services in New Zealand : conditions and trends. Lincoln: Manaaki Whenua Press; 2013. p. 408-431. Available from: http://www.mwpress.co.nz/__data/assets/pdf_file/0008/77057/2_11_Newstrom.pdf. [Last accessed on 26 Jul 2023]

[3] Ministry for Primary Industries. 2021 apiculture monitoring data. Available from: https://www.mpi.govt.nz/dmsdocument/48793-2021-Apiculture-monitoring-report-data. [Last accessed on 26 Jul 2023].

[4] Li G,Liu Z,Xu B.The wisdom of honeybee defenses against environmental stresses.Front Microbiol2018;9:722 PMCID:PMC5938604

[5] Genersch E.American Foulbrood in honeybees and its causative agent, Paenibacillus larvae.J Invertebr Pathol2010;103:S10-9

[6] Genersch E.Paenibacillus larvae and American Foulbrood - long since known and still surprising.J Verbr Lebensm2008;3:429-34

[7] Alippi AM,López AC,Aguilar OM.Molecular epidemiology of Paenibacillus larvae larvae and incidence of American Foulbrood in Argentinean honeys from Buenos Aires province.J Apic Res2004;43:135-43

[8] Rauch S,Hedtke K.Negative correlation between individual-insect-level virulence and colony-level virulence of Paenibacillus larvae, the etiological agent of American Foulbrood of honeybees.Appl Environ Microbiol2009;75:3344-7 PMCID:PMC2681656

[9] Alippi AM,Aguilar OM.Differentiation of Paenibacillus larvae subsp. larvae, the cause of American Foulbrood of honeybees, by using PCR and restriction fragment analysis of genes encoding 16S rRNA.Appl Environ Microbiol2022;68:3655-60 PMCID:PMC126810

[10] Lester P. Healthy bee, sick bee: the influence of parasites, pathogens, predators and pesticides on honey bees. Victoria University of Wellington Press; 2021. Available from: https://play.google.com/books/reader?id=fY4bEAAAQBAJ&pg=GBS.PP1&hl=en. [Last accessed on 26 Jul 2023]

[11] Kok DN.Save our bees: bacteriophages to protect honey bees against the pathogen causing American Foulbrood in New Zealand. N Z J Zool 2023;1-16.

[12] Goodwin M.American Foulbrood control: the New Zealand approach.Bee World2005;86:44-5

[13] The management agency national american foulbrood pest management plan. Available from: https://afb.org.nz/wp-content/uploads/2018/07/BRIEFING-DOCUMENT-MPI-Government-01112017.pdf. [Last accessed on 31 Jul 2023]

[14] Biosecurity (National American Foulbrood Pest Management Plan) Order 1998. 1998. p. 1-21. Available from: https://afb.org.nz/wp-content/uploads/2018/10/Biosecurity-National-American-Foulbrood-Pest-Management-Plan-Order-1998.pdf. [Last accessed on 25 Jul 2023]

[15] Mushegian AR.Are there 10³¹ virus particles on earth, or more, or fewer?.J Bacteriol2020;202:e0052-20 PMCID:PMC7148134

[16] Hendrix RW,Burns RN,Hatfull GF.Evolutionary relationships among diverse bacteriophages and prophages: all the world’s a phage.Proc Natl Acad Sci USA.1999;96:2192-2197 PMCID:PMC26759

[17] Yost DG,Amy PS.Experimental bacteriophage treatment of honeybees (Apis mellif era) infected with Paenibacillus larvae, the causative agent of American Foulbrood Disease.Bacteriophage2016;6:e1122698 PMCID:PMC4836486

[18] Brady TS,Hilton JA,Stephenson MB.Bacteriophages as an alternative to conventional antibiotic use for the prevention or treatment of Paenibacillus larvae in honeybee hives.J Invertebr Pathol2017;150:94-100

[19] Stamereilers C,Walker JK.Genomic analysis of 48 Paenibacillus larvae bacteriophages.Viruses2018;10:377 PMCID:PMC6070908

[20] Graham SAM. American foulbrood and its causative agent, Paenibacillus larvae, in new zealand’s registered hives and apiaries. Available from: https://openaccess.wgtn.ac.nz/articles/thesis/American_foulbrood_and_its_causative_agent_Paenibacillus_larvae_in_New_Zealand_s_registered_hives_and_apiaries/17013008. [Last accessed on 26 Jul 2023]

[21] Dingman DW.Medium promoting sporulation of Bacillus larvae and metabolism of medium components.Appl Environ Microbiol1983;46:860-9 PMCID:PMC239480

[22] de Graaf DC,Antúnez K.Standard methods for American Foulbrood research.J Apic Res2013;52:1-28

[23] Dobbelaere W,Peeters JE.Development of a fast and reliable diagnostic method for American Foulbrood disease (Paenibacillus larvae subsp. larvae) using a 16S rRNA gene based PCR.Apidologie2001;32:363-70

[24] Bankevich A,Antipov D.SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing.J Comput Biol2012;19:455-77 PMCID:PMC3342519

[25] Prjibelski A,Meleshko D,Korobeynikov A.Using SPAdes De Novo assembler.Curr Protoc Bioinformatics2020;70:e102

[26] Overbeek R,Pusch GD.The SEED and the Rapid Annotation of microbial genomes using Subsystems Technology (RAST).Nucleic Acids Res2014;42:D206-14 PMCID:PMC3965101

[27] Brettin T,Disz T.RASTtk: a modular and extensible implementation of the RAST algorithm for building custom annotation pipelines and annotating batches of genomes.Sci Rep2015;5:8365 PMCID:PMC4322359

[28] Aziz RK,Best AA.The RAST Server: rapid annotations using subsystems technology.BMC Genomics2008;9:75 PMCID:PMC2265698

[29] Seemann T.Prokka: rapid prokaryotic genome annotation.Bioinformatics2014;30:2068-9

[30] Kok DN,Takeuchi N,Hendrickson HL.In Vitro evolution to increase the titers of difficult bacteriophages: RAMP-UP protocol.Phage2023;4:68-81 PMCID:PMC10282794

[31] Santos MA.An improved method for the small scale preparation of bacteriophage DNA based on phage precipitation by zinc chloride.Nucleic Acids Res1991;19:5442 PMCID:PMC328918

[32] Lazeroff M,Harris SL.Phage commander, an application for rapid gene identification in bacteriophage genomes using multiple programs.Phage2021;2:204-13 PMCID:PMC9041506

[33] Pope WH.Annotation of bacteriophage genome sequences using DNA master: an overview. In: Clokie MR, Kropinski AM, Lavigne R, editors. Bacteriophages. New York: Springer; 2018. p. 217-29.

[34] Salisbury A.A method for improving the accuracy and efficiency of bacteriophage genome annotation.Int J Mol Sci2019;20:3391 PMCID:PMC6678273

[35] Jolley KA,Maiden MCJ.Open-access bacterial population genomics: BIGSdb software, the PubMLST.org website and their applications.Wellcome Open Res2018;3:124 PMCID:PMC6192448

[36] Morrissey BJ,Poppinga L,Genersch E.Biogeography of Paenibacillus larvae, the causative agent of American Foulbrood, using a new multilocus sequence typing scheme.Environ Microbiol2015;17:1414-24 PMCID:PMC4405054

[37] Papić B,Kušar D.Analysis of the global population structure of Paenibacillus larvae and outbreak investigation of American Foulbrood using a stable wgMLST scheme.Front Vet Sci2021;8:582677 PMCID:PMC7952629

[38] Binney BM,Foxwell J.Genomic analysis of the population structure of Paenibacillus larvae in New Zealand.Front Microbiol2023;14:1161926 PMCID:PMC10157257

[39] Grissa I,Pourcel C.CRISPRFinder: a web tool to identify clustered regularly interspaced short palindromic repeats.Nucleic Acids Res2007;35:W52-7 PMCID:PMC1933234

[40] Tesson F,Mordret E.Systematic and quantitative view of the antiviral arsenal of prokaryotes.Nat Commun2022;13:2561 PMCID:PMC9090908

[41] Abby SS,Ménager H,Rocha EP.MacSyFinder: a program to mine genomes for molecular systems with an application to CRISPR-Cas systems.PLoS One2014;9:e110726 PMCID:PMC4201578

[42] Oliveira PH,Rocha EP.The interplay of restriction-modification systems with mobile genetic elements and their prokaryotic hosts.Nucleic Acids Res2014;42:10618-31 PMCID:PMC4176335

[43] Gao L,Böhning F.Diverse enzymatic activities mediate antiviral immunity in prokaryotes.Science2020;369:1077-84 PMCID:PMC7985843

[44] Bernheim A,Touchon M.Atypical organizations and epistatic interactions of CRISPRs and cas clusters in genomes and their mobile genetic elements.Nucleic Acids Res2020;48:748-60 PMCID:PMC7145637

[45] Doron S,Ofir G.Systematic discovery of antiphage defense systems in the microbial pangenome.Science2018;359:eaar4120 PMCID:PMC6387622

[46] Millman A,Leavitt A.An expanded arsenal of immune systems that protect bacteria from phages.Cell Host Microbe2022;30:1556-69

[47] Arndt D,Marcu A.PHASTER: a better, faster version of the PHAST phage search tool.Nucleic Acids Res2016;44:W16-21 PMCID:PMC4987931

[48] Zhou Y,Lynch KH,Wishart DS.PHAST: a fast phage search tool.Nucleic Acids Res2011;39:W347-52 PMCID:PMC3125810

[49] Tsourkas PK,Krohn A.Complete genome sequences of nine phages capable of infecting Paenibacillus larvae, the causative agent of American Foulbrood disease in honeybees.Genome Announc2015;3:e01120-15 PMCID:PMC4611677

[50] Tsourkas PK.Paenibacillus larvae bacteriophages: obscure past, promising future.Microb Genom2020;6:e000329 PMCID:PMC7067210

[51] Ebeling J,Genersch E.The buzz about ADP-ribosylation toxins from Paenibacillus larvae, the causative agent of American Foulbrood in honey bees.Toxins2021;13:151 PMCID:PMC7919650

[52] Fünfhaus A,Genersch E.Identification and characterization of two novel toxins expressed by the lethal honey bee pathogen Paenibacillus larvae, the causative agent of American Foulbrood.Environ Microbiol2013;15:2951-65

[53] Ackermann HW.Bacteriophage observations and evolution.Res Microbiol2003;154:245-51

[54] Ministry for Primary Industries. ApiWellbeing. Available from: https://www.mpi.govt.nz/biosecurity/how-to-find-report-and-prevent-pests-and-diseases/bee-biosecurity/apiwellbeing/. [Last accessed on 26 Jul 2023]

[55] King C.American foulbrood.Surveillance2020;47:42Available from: https://www.mpi.govt.nz/dmsdocument/43978. [Last accessed on 31 Jul 2023]

[56] Abedon ST,Wozniak DJ.Phage cocktail development for bacteriophage therapy: toward improving spectrum of activity breadth and depth.Pharmaceuticals2021;14:1019 PMCID:PMC8541335

[57] Stamereilers C,Tsourkas PK.Characterization of CRISPR spacer and protospacer sequences in Paenibacillus larvae and its bacteriophages.Viruses2021;13:459 PMCID:PMC7998209

[58] Ribeiro HG,Melo LDR.Analysis of intact prophages in genomes of Paenibacillus larvae: an important pathogen for bees.Front Microbiol2022;13:903861 PMCID:PMC9341999

[59] Hille F,Wong SP,Ressel S.The biology of CRISPR-Cas: backward and forward.Cell2018;172:1239-59

[60] Loenen WA.The other face of restriction: modification-dependent enzymes.Nucleic Acids Res2014;42:56-69 PMCID:PMC3874153

[61] Gochnauer T A.Some properties of a bacteriophage from Bacillus larvae.J Invertebr Pathol1970;15:149-156

[62] Merrill BD,Hilton JA.Complete genome sequences of 18 Paenibacillus larvae phages from the western United States.Microbiol Resour Announc2018;7:e00966-18 PMCID:PMC6256562

[63] Jończyk-Matysiak E,Popiela E.Isolation and characterization of phages active against Paenibacillus larvae causing American Foulbrood in honeybees in Poland.Viruses2021;13:1217 PMCID:PMC8310151

[64] Carson S,DeFoor W.Genome sequences of six Paenibacillus larvae siphoviridae phages.Genome Announc2015;3:e00101-15 PMCID:PMC4472882

[65] Walker JK,Berg JA.Complete genome sequences of Paenibacillus larvae phages BN12, Dragolir, Kiel007, Leyra, Likha, Pagassa, PBL1c, and Tadhana.Genome Announc2018;6:e01602-17 PMCID:PMC6003738

[66] Beims H,Bunk B.Paenibacillus larvae-directed bacteriophage HB10c2 and its application in American Foulbrood-affected honey bee larvae.Appl Environ Microbiol2015;81:5411-9 PMCID:PMC4510184

[67] Oliveira A,Kropinski AM.Complete genome sequence of the broad-host-range Paenibacillus larvae phage phiIBB_Pl23.Genome Announc2013;1:e00438-13 PMCID:PMC3764407

[68] Ribeiro HG,Oliveira H.Characterization of a new podovirus infecting Paenibacillus larvae.Sci Rep2019;9:20355 PMCID:PMC6937236

[69] Yost DG,LeBlanc L.Complete genome sequences of Paenibacillus larvae phages halcyone, heath, scottie, and unity from Las Vegas, Nevada.Microbiol Resour Announc2018;7:e00977-18 PMCID:PMC6256684

[70] Citizen Phage Library. Available from: https://www.citizenphage.com. [Last accessed on 26 Jul 2023]

[71] Forti F,Cafora M.Design of a broad-range bacteriophage cocktail that reduces Pseudomonas aeruginosa biofilms and treats acute infections in two animal models.Antimicrob Agents Chemother2018;62:e02573-17 PMCID:PMC5971607

[72] Niu YD,Du H.Efficacy of individual bacteriophages does not predict efficacy of bacteriophage cocktails for control of Escherichia coli O157.Front Microbiol2021;12:616712 PMCID:PMC7943454