Genome Sequences of Six Paenibacillus larvae Siphoviridae Phages

GENOME ANNOUNCEMENT

Phages were isolated from infected comb presented to the North Carolina Department of Agriculture from across the state. Phage isolation and propagation used Paenibacillus larvae host strains ATCC 9545 and ATCC 25747 grown on brain heart infusion (BHI) plus 0.4% glucose and 1 µg/ml thiamine in broth or agar at 30°C. Swabs were taken from infected comb cells and incubated in growth medium with selective bacteria for 24 to 48 h. Plaques were identified on top agar lawns of P. larvae and reisolated from plate streaks at least three times. Electron micrographs showed that all six phages have the Siphoviridae morphotype.

Confluent lysis on P. larvae-seeded top agar was used to prepare lysates from which DNA was extracted using the Promega DNA cleanup system. Genomic DNA was sequenced in the North Carolina State University (NCSU) Genomic Sciences Laboratory using fragmentation, sizing, adapter addition, and flow cell usage, as described (Illumina, Inc., San Diego, CA). Lily and Rani were sequenced by HiSeq 100-bp reads, with the remaining four genomes being sequenced using MiSeq 2 × 300-bp reads. CLC Genomics Workbench releases 2013 and 2014 were used for assembly. The physical ends of the Lily and Rani genomes were determined by ligation of genomic DNA, PCR amplification across the joined ends, and Sanger sequencing. All use cos packaging and are likely temperate phages. Lily has a 12-bp 5′ overhanging terminus (GGTGCGCGTGAG), and Rani has a 9-bp 3′ overhanging terminus (CGACTGCCC). Diva, Redbud, Shelly, and Sitara have physical genome ends like those of Rani, based on matching end nucleotides and similar sequence assembly patterns.

The genomes were annotated by students using DNA Master (http://cobamide2.bio.pitt.edu) on the NC State Virtual Computing Lab. Coding regions were predicted using Glimmer (1) and GeneMark (2), and start codons were chosen based on DNA Master ribosomal binding site (RBS) parameters. The absence of tRNAs was predicted using ARAGORN (3) and tRNAscan (4). Protein function and start sites were corroborated using NCBI BLASTp (5).

All six genomes encode two terminase subunits oriented at the 5′ end of the sequence, followed by genes for the tail and head that are generally syntenic. Rani and Redbud are nearly identical, and Lily is the most divergent. All contain genes typically seen in Siphoviridae phages (major capsid, portal, tape measure, tail, holin, endolysin, etc.). Most of the phages share regions of similarity with phiIBB_Pl23, a P. larvae Siphoviridae phage (6), and negligible similarity with the Myoviridae phages isolated from Utah (7). That all of the reported phages from Utah are Myoviridae and the North Carolina phages are Siphoviridae may reflect that the North Carolina phages were uniquely isolated from American foulbrood (AFB) diseased hives or that the Utah host used is substantially different from that of the P. larvae strains used in this work. Phage HB10c2 (GenBank accession no. KP202972), isolated from an AFB-diseased hive in Germany, has substantial sequence similarity to the genomes reported here.