Complete Genome Sequences of Five Paenibacillus larvae Bacteriophages

Paenibacillus larvae is a pathogen of honeybees that causes American foulbrood (AFB). We isolated bacteriophages from soil containing bee debris collected near beehives in Utah. We announce five high-quality complete genome sequences, which represent the first completed genome sequences submitted to GenBank for any P. larvae bacteriophage.

pathogen that causes American foulbrood (AFB). AFB kills honeybee larvae (1), contributes to colony collapse disorder (2), and limits agricultural yields (3). Unfortunately, some P. larvae strains have become resistant to the antibiotics typically used for AFB treatment (4). Phage therapy is a potential treatment for AFB, yet few P. larvae-specific phages have been described (5)(6)(7), and a full-genome sequence for one has only recently become available (8).
Soil, honey, and larva samples were collected in the Utah, Salt Lake, and Davis counties of Utah. The samples were used for isolating new host strains or bacteriophages. P. larvae subsp. pulvifaciens (9) hosts were confirmed by 16S rRNA sequencing. The bacterial cultures were inoculated with soil samples to enrich for bacteriophages. Plated enrichment samples formed plaques from which bacteriophages were selected and purified by a minimum of three passages. High-titer lysates were filtered, incubated with 5 g/ml RNase and 10 g/ml DNase for 30 min at 37°C, and treated with 100 g/l proteinase K at 52°C for 1 h. Following phenol-chloroform extraction and ethanol precipitation, highquality DNA (10) was sequenced using 454 pyrosequencing.
No experiments were performed to determine the physical ends or the packing or replication strategies of the phage DNA. However, during manual finishing, overlapping contigs assembled the genome into an apparently circular genome. For the bacteriophages Jimmer1, Jimmer2, Davies, and Abouo, the first base of each genome was selected in the noncoding gap between the terminase gene and the prior gene. Since Emery did not have a terminase small subunit gene, the first base was selected in the first gap upstream of the large subunit terminase gene. The noncoding gap where the first base was selected for a genome contained multiple stops and lacked coding potential in any frame in all genomes.
Annotation was completed using DNA Master (http://cobamide2 .bio.pitt.edu). A coding potential map was generated using GeneMark 2.5p (14) for each phage based on Bacillus cereus strain ATCC 14579, the closest available relative to P. larvae. Our selection of gene calls emphasized the following criteria: GeneMark HHM and Glimmer autoannotation, BLAST alignment E values of Ͻ0.001, coding potential from GeneMark, start codon sequences, and Shine-Dalgarno (SD) scores of Ͼ200 nats using the Karlin position-specific scoring matrix (PSSM) for moderately to highly expressed genes.
All five phages were identified as myoviruses. Jimmer1 and Jimmer2 were isolated independently from the same soil sample and differ in their genomic sequences by only 80 bp of the 54,312-bp genomes (99.85% similarity). The differences between these related phages are real base pair changes because the sequencing fold coverages of these samples are statistically greater than the error rate of 454 sequencing (15). Nucleotide sequence accession numbers. The GenBank accession numbers for the five Paenibacillus larvae bacteriophages are listed in Table 1.