key: cord-0824454-ifok0gm9
authors: Chen, Yibao; Yang, Lan; Sun, Erchao; Song, Jiaoyang; Wu, Bin
title: Characterisation of a newly detected bacteriophage infecting Bordetella bronchiseptica in swine
date: 2018-09-18
journal: Arch Virol
DOI: 10.1007/s00705-018-4034-0
sha: 2f091e6686b69583cabbb96a3cf80e41a6431578
doc_id: 824454
cord_uid: ifok0gm9

A novel virulent bacteriophage, vB_BbrM_PHB04, infecting Bordetella bronchiseptica was isolated from wastewater collected at a swine farm in China. Phage vB_BbrM_PHB04 exhibited growth over a wide range of temperature and pH conditions and showed different efficiency of plating values and lytic spectra within the same strains at 25 °C and 37 °C. High-throughput sequencing revealed that vB_BbrM_PHB04 has a linear double-stranded DNA genome with 124 putative open reading frames. Overall, the genome of vB_BbrM_PHB04 showed very low similarity (the highest nucleotide identity 82%, 1% coverage) to other phage sequences in the GenBank database. Phylogenetic analysis indicated that vB_BbrM_PHB04 is a new member of the family Myoviridae. In addition, polymerase chain reaction-based detection of phage genes in phage-resistant B. bronchiseptica variants revealed no evidence of lysogenic activity of phage vB_BbrM_PHB04. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (10.1007/s00705-018-4034-0) contains supplementary material, which is available to authorized users.

Bordetella bronchiseptica is a widespread Gram-negative pathogen that can be isolated from many different mammalian species. It is one of the primary causative agents of infectious respiratory diseases such as swine pneumonia and atrophic rhinitis [1] . In addition, B. bronchiseptica is often isolated in combination with other pathogens from swine suffering from pneumonia. Previous studies have shown that co-infection of viral and bacterial pathogens, including porcine reproductive and respiratory syndrome virus, porcine respiratory coronavirus, Haemophilus parasuis, Pasteurella multocida, and Streptococcus suis, exacerbates the severity of respiratory disease in pigs [2] . Importantly, human B. bronchiseptica infections have also been reported [3, 4] .

B. bronchiseptica virulence factors are regulated by the complex Bordetella virulence gene product BvgAS [5] . BvgAS is a two-component signal transduction system that regulates the expression of multiple virulence factors in response to changes in environmental conditions. Different virulence factors are expressed in the three different B. bronchiseptica phases: the bvg + , bvg i , and bvg − phases. In the virulent bvg + phase, which occurs at 37 °C in the absence of MgSO 4 and nicotinic acid, B. bronchiseptica can rapidly colonize the respiratory tract of its host through the expression of adhesins, autotransporters, and toxins. The avirulent bvg − phase occurs at growth temperatures of ~25 °C or in the presence of MgSO 4 or nicotinic acid. Finally, the intermediate (bvg i ) phase is responsible for biofilm formation by B. bronchiseptica [5, 6] .

With the increasing prevalence of multidrug-resistant bacteria, several reports have indicated the emergence of B. bronchiseptica isolates with decreased susceptibility to β-lactam antibiotics such as penicillins and cephalosporins [7] . To combat multidrug-resistant bacteria, phages offer a potential alternative to antibiotics. In the past, temperate B. bronchiseptica phages have been induced in clinical bacterial strains, but their phenotypic characteristics have rarely been described. Petrovic et al [8] . isolated 29 lytic B. bronchiseptica phages from the environment and described their phenotypic characteristics; however, the genomes of these phages were only partially sequenced to determine whether they were lytic or temperate phages. Here, we isolated a novel lytic phage from wastewater collected from a swine farm in China and studied its phenotypic and genetic characteristics. To our knowledge, the complete genomes of only three temperate B. bronchiseptica phages have been sequenced [9] . 

The conventional double-layer agar method was used to examine phage activity [10] . A sewage water sample (20 ml) was collected in Wuhan, China, and then sterilized by filtration through a 0.22-μm membrane. A 5-ml aliquot of the filtrate was mixed with 10 ml of an exponential phase B. bronchiseptica 4715 culture, and the mixture was incubated at 37 °C for 8 h. The culture was then centrifuged at 12,000 × g for 10 min, and the resulting supernatant was filtered using a 0.22-μm membrane. A 300-μl aliquot of the filtered supernatant was then mixed with 500 μl of B. bronchiseptica 4715, added to 6 ml of molten soft TSB (0.75% agar), supplemented with 10% sterile, defibrinated sheep blood and then poured onto the surface of a prepared TSA plate supplemented with 10% (v/v) sterile, defibrinated sheep blood. Plaques appearing after an overnight incubation at 37 °C were enumerated [10] .

A single plaque was picked aseptically and re-suspended at 180 rmp in 6 ml of SM buffer (5.8 g of NaCl, 2.0 g of MgSO 4 ·7H 2 O, 50 ml of Tris-HCl [pH7.4], 5.0 ml of 2% gelatin) for 3 h. The phage-containing SM buffer was collected and centrifuged at 12,000 × g for 30 s, and then the supernatant was filtered through a 0.22-μm membrane. This phage preparation was subjected to serial 10-fold dilutions in sterile SM buffer. Phage isolation by the double-layer agar method was repeated four more times, and the isolates were stored at 4 °C.

The phage filtrate was stained negatively with 2% uranyl acetate after addition of a drop of a phage suspension onto a grid surface, and then the excess stain was removed immediately. The morphology of the phages was examined under a 100-kV transmission electron microscope (HITACHI H-7650, Japan).

To examine these properties of the phage, we used a known method [10] . Briefly, stability of the purified phage at various temperatures (4, 20, 40, 50, 60, 70, and 80 °C) was assessed by incubation for 1 h. To evaluate the stability of the phage at different pH levels (3, 5, 7, 9, and 11) , the purified phage particles were incubated at 37 °C for 1 h. Each assay was performed in triplicate, and then the samples were titered by the double-layer agar plate method. 

To analyze the infectivity and replication ability of phages, a one-step growth curve test was performed as previous described [10] . Briefly, phage vB_BbrM_PHB04 was mixed with the host strain in exponential phase at a multiplicity of infection (MOI) of 0.01 and incubated at 37 °C for 5 min. The mixtures were then centrifuged at 12,000 × g for 30 s to remove any unabsorbed free phage. The precipitate was washed with 37 °C TSB, then the suspension was transferred to 30 ml of TSB supplemented with 10% (v/v) sterile, defibrinated sheep blood followed by incubation at 25 °C and 37 °C. This time point was defined as t = 0 s. At 10-min intervals, a 0.5-ml sample was collected for a total of 200 min. The titration of phage particles was conducted by the double-layer agar method. The experiment was repeated three times.

The bacterial strain susceptibility levels were detected using spotting methods as previously described [10] . 

The phages' genomic DNA was extracted using the phenolchloroform protocol [11] . Briefly, purified phages were lysed by the addition of proteinase K (100 mg/ml), SDS (10%, wt/ vol) and EDTA (0.5 mM, pH 8.0) and incubation at 56 °C in water for 2 h. The digested sample was mixed with an equal volume of phenol/chloroform/isoamyl alcohol (25:24:1) centrifuged at 12,000× g (4 °C) for 10 min to remove any debris. This step was repeated three times. The supernatant was mixed with isoamyl alcohol kept at -20 °C overnight. Then the air-dried precipitate was washed three times with cold 75% ethanol. Finally, phages' genomic DNA was dissolved in TE buffer (10 mM Tris-HCl, 1 mM EDTA [pH 8.0]) and was sequenced on an Illumina HiSeq 2500 sequencer with 2×100 bp read length. More than 6,200-fold coverage of phage genome was achieved by sequencing the cloned fragments. The short reads were assembled into the genome by means of SOAPdenovo (http://sourceforge.net/projects/ soapdenovo2/files/SOAPdenovo2/) [12] . Protein-encoding putative open reading frames (ORFs) were predicted using Glimmer [13, 14] . The final assembled sequence was searched against the current protein and nucleotide databases (http://www.ncbi.nlm.nih.gov/) by means of the basic local alignment search tool (BLAST). Protein BLAST (BLASTP) (http://www.ncbi.nlm.nih.gov/BLAST /) was used to identify putative homologies and proteins sharing similarities with the predicted phage proteins. Phylogenetic analysis of phage large terminase subunit sequences was performed using the ClustalW program in MEGA 6 [15] .

To test for possible lysogens, host B. bronchiseptica 4715 cells in exponential growth phase (approximately 10 7 colony-forming units/ml) were incubated with phage vB_ BbrM_PHB04 at a multiplicity of infection (MOI) of 100 in SM buffer and the mixture was incubated for 10 min to allow adsorption of phages [16] . The conventional double-layer agar method was used to pick 40 random phage-resistant colonies. Single colonies were cultured with TSB supplemented with 10% (v/v) sterile, defibrinated sheep blood at 37 °C for 24 h and examined for phage vB_BbrM_PHB04 sensitivity by spot testing, applying 10 μl (10 7 pfu/ml) phage to visualize individual plaques. Genomic DNA was extracted from each of the bacteria resistance to phage using a commercial extraction kit (Roche, Switzerland). The PCR primers used in this study are described in Table S1 . Primers were tested for specificity using the phage vB_BbrM_PHB04 genome as a positive control. This experiment was repeated three times.

Phage vB_BbrM_PHB04 was isolated from its host, B. bronchiseptica strain 4715, which was isolated from a sewage sample from a swine farm in China. Transmission electron microscopy of the purified and concentrated phage particles showed that phage vB_BbrM_PHB04 had an isometric head (86.3 ± 1.3 nm in diameter) and a contractile tail (72.4 ± 3.2 nm long and 23.3 ± 1.8 nm wide) (Fig 1) . Based on the current International Committee on Taxonomy of Viruses classification system, the phage was assigned to the family Myoviridae, order Caudovirales.

Phage vB_BbrM_PHB04 was stable at temperatures ranging from 4-50 °C, but a reduction in viral titers occurred at 60-70 °C, and titers were dramatically reduced at 80 °C (Fig 2A) . The phage was stable over a broad pH range (from 5.0-9.0), but the titer dropped by approximately 1.2 log at pH 3.0 and pH 11.0 (Fig 2B) . A one-step growth experiment was conducted to determine the duration of the latent period and burst size of phage vB_BbrM_ PHB04. The phage had a long latent period at both 25 °C and 37 °C (~130 min and ~140 min, respectively), with an average burst size of 51 and 40 phage particles, respectively, per infected cell after 180 min (Fig 3) . 

High-throughput sequencing of phage genomic DNA was performed using the Illumina HiSeq system according to the manufacturer's instructions. The assembly of the random library of sequences yielded a closed, circular genome by means of SOAPdenovo. The complete linear, doublestranded DNA genome of phage vB_BbrM_PHB04 consisted of 94,005 bp, with a GC content of 65.5% (Fig 4) . BLASTn analysis revealed that the genome of phage vB_ BbrM_PHB04 showed the highest nucleotide identity (82%, though only 1% coverage) to the genome sequence of Burkholderia phage strain KACC 11889. This low sequence similarity suggested that phage vB_BbrM_PHB04 is a novel species. The genome contained 124 putative open reading frames (ORFs), and tRNA prediction, conducted using tRNAscan-SE analysis, identified three tRNA encoding genes: tRNA-Met-CAT, tRNA-Asn-GTT, and tRNA-Pro-TGG (Fig 4) . The predicted ORFs encoding proteins with known functions were classified into three groups: those associated with morphogenesis and structure, those associated with DNA replication and metabolism, and those associated with lysis or lysogeny (Fig 4) . Proteins involved in structure and assembly were encoded by genes in the phage vB_BbrM_PHB04 genome. BLASTn analysis of these regions, including the baseplate (ORF3), major tail sheath (ORF12), major capsid protein (ORF20), and tail fiber protein (ORF119, ORF121, ORF122, and ORF124), showed low identity to the corresponding proteins in other phage (Fig 4; Table S1 ). Genes encoding DNA replication and regulation modules such as the helicase (ORF25), DNA primase (ORF98), ssDNA-binding protein (ORF101), and RecD-like DNA helicase (ORF112) were located in the middle of the phage vB_BbrM_PHB04 genome (Fig 4) . However, all of these sequences showed only low levels of identity to genes encoding phage proteins of similar function. Additionally, we identified 96 genes of unknown function. With respect to lysis or lysogeny modules, ORF117 encoded a protein with 56% amino acid sequence identity to the endolysin protein of Pseudomonas phage LkO4 (Fig 4; Table S2 ). Holin, a protein that forms large holes in the cytoplasmic membrane of phage-infected bacteria is encoded by ORF118, showed 42% amino acid sequence identity to the holin protein from Pseudomonas phage LkO4 (Fig 4; Table S2 ). In addition, the protein encoded by ORF8 contained domains from multiple families, including LT_GEWL, SLT, MitE, PRK11619, and DUF3584 (Fig 4; Table S2 ). Proteins LT_GEWL and SLT belong to the lysozyme-like superfamily, all members of which participate in the hydrolysis of β-1,4-linked polysaccharides.

The amino acid sequences of the terminase large subunit proteins (encoded by ORF55) from phage vB_BbrM_ PHB04 and other phages belonging to the Myoviridae family were aligned using MEGA6 software. A neighbor-joining tree generated from these sequences revealed that phage vB_BbrM_PHB04 clustered with phages such as Bacillus phage SPO1, Bacillus phage TsarBomba, and Staphylococcus phage K in the subfamily Spounavirinae (Fig 5) . This clear phylogenetic distinction from other phages and the novel characteristics of the genome indicate that phage vB_ BbrM_PHB04 is a new member of the family Myoviridae, order Caudovirales.

The lysogenic activity of phage genes in phage-resistant bacteria was examined by PCR. In each round of screening, 40 random phage-resistant colonies were picked for further analysis and examined for vB_BbrM_PHB04 sensitivity by spot testing. The PCR results showed that the phage major capsid gene (ORF20), the terminase gene (ORF55), the helicase gene (ORF102), and the integrase gene (ORF113) were all absent from the genomes of phage-resistant colonies, suggesting that phage vB_BbrM_PHB04 does not have lysogenic activity in host bacteria (the results are not shown).

In this study, we initially screened 17 domestic sewage samples for the presence of B. bronchiseptica phage, but failed to identify the phage in any of the samples (the results are not shown). Petrovic et al [8] . also reported a failed attempt to isolate B. bronchiseptica phage from samples collected from swine water troughs and farm wastewater. Interestingly, although we eventually isolated phage vB_BbrM_PHB04 from swine wastewater, the farm in question had no reported outbreaks or even suspected cases of bronchopneumonia or atrophic rhinitis in recent years. This may therefore indicate latent B. bronchiseptica infection of pigs.

Phage vB_BbrM_PHB04 was assigned to the family Myoviridae based on its morphological characteristics. In addition, phylogenetic analysis of the phage-encoded large terminase subunit protein revealed a close phylogenetic relationship between phage vB_BbrM_PHB04 and phages belonging to the Spounavirinae subfamily. Liu et al [9] . reported on three temperate phages isolated from B. bronchiseptica that were assigned to the family Podoviridae, while Petrovic et al [8] . isolated 29 B. bronchiseptica phages assigned to the family Siphoviridae. Based on its morphology and genome features, phage vB_BbrM_PHB04 Phylogenetic tree analysis of the phage-encoded large terminase subunit (ORF55) of phage vB_BbrM_PHB04. The amino acid sequences of the large terminase subunits were compared using MEGA6 in the analysis, and phylogenetic tree was generated using the neighbour-joining method and 1000 bootstrap replicates was confirmed as a new member of the family Myoviridae, order Caudovirales.

One-step growth curve analysis revealed a significant difference in the duration of the latent period and burst size of phage vB_BbrM_PHB04 at 25 °C vs. 37 °C, which agrees with the findings of Petrovic et al [8] ., who determined that the latent period and burst size of phage vB_BbrS_CN1 is dependent on temperature. In addition, the EOP values of phage vB_BbrM_PHB04 showed slight differences in phage host range at the two temperatures. B. bronchiseptica displays phenotypic variation in response to temperature change, and this phase variation may have an impact on phage receptor expression. Moreover, phage vB_BbrM_ PHB04 had higher EOP values in mutant B. bronchiseptica strain 3814ΔbvgA than in wild-type strain 3814. We hypothesize that the loss of the virulence gene regulatory protein in the mutant strain resulted in enhanced vB_BbrM_PHB04 adsorption and subsequent infection.

The genome sequence of phage vB_BbrM_PHB04 showed the highest degree of nucleotide sequence identity to the genome sequence of Burkholderia phage strain KACC 11889. Functional predictions were possible for 27 of the predicted ORFs, while the remaining ORFs exhibited limited amino acid sequence identity to other hypothetical proteins, indicating that the phage vB_BbrM_PHB04 genome could be unique. In general, genome recombination occurs frequently between bacteria and phages and is conducive to phage adaptation to the host bacterium. The protein encoded by ORF113 showed 45% amino acid sequence identity to the putative integrase family protein of Pseudomonas phage PPpW-3. We know that temperate phage genomes can be integrated into bacterial genomes or plasmids and that prophages confer immunity that protects host bacteria against infection by other phages. In addition, the integrase gene is usually associated with lysogenic activity. To determine whether phage vB_BbrM_PHB04 is temperate, we designed four pairs of primers targeting essential phage genes. The results showed that none of the 40 phage-resistant bacterial clones were positive for the essential phage genes, suggesting that phage vB_BbrM_PHB04 does not have lysogenic activity.

Although no known lysogenic factors such as an excisionase, or anti-repressor proteins were identified in the genome of phage vB_BbrM_PHB04, the large number of predicted proteins with no known function may carry out lysogenrelated functions. Phages containing lysogen-related genes, but not being temperate themselves, have been described previously. For example, while phage UAB_Phi20 contains all the elements needed for a lysogenic cycle, phage PPpW-3 only contains an integrase gene; however, both are lytic phages [16, 17] .

We did not consider the spontaneous release of temperate phage or prophage induction from the host strain for life cycle of phage vB_BbrM_PHB04 assay. Moreover, an endolysin-encoding gene was identified in the phage vB_ BbrM_PHB04 genome. Future work should therefore be conducted to express the endolysin, and to examine potential applications for this enzyme in the treatment of multidrugresistant B. bronchiseptica.

Biology of Bordetella bronchiseptica

Prior infection with Bordetella bronchiseptica increases nasal colonization by Haemophilus parasuis in swine

Cavitary pneumonia in an AIDS patient caused by an unusual Bordetella bronchiseptica variant producing reduced amounts of pertactin and other major antigens

Bordetella bronchiseptica in the immunosuppressed population-A case series and review

The bvgAS locus negatively controls motility and synthesis of flagella in Bordetella bronchiseptica

Evaluation of the role of the Bvg intermediate phase in Bordetella pertussis during experimental respiratory infection

Antimicrobial susceptibility of Bordetella bronchiseptica isolates from swine and companion animals and detection of resistance genes

The first Siphoviridae family bacteriophages infecting Bordetella bronchiseptica isolated from environment

Reverse transcriptase-mediated tropism switching in Bordetella bacteriophage

Complete genome sequence of a novel T7-like bacteriophage from a Pasteurella multocida capsular type A isolate

A novel roseosiphophage isolated from the oligotrophic South China Sea

De novo assembly of human genomes with massively parallel short read sequencing

Microbial gene identification using interpolated Markov models

Identifying bacterial genes and endosymbiont DNA with Glimmer

MEGA6: molecular evolutionary genetics analysis version 6.0

Complete genome sequence analysis of two Pseudomonas plecoglossicida phages, potential therapeutic agents

Genomics of three new bacteriophages useful in the biocontrol of Salmonella

The complete genome sequence of phage vB_ BbrM_PHB04 with annotation has been deposited in GenBank under accession number MF663786.Funding This study was funded by Research projects of agricultural public welfare industry of China (201403054), and National Science supported planning projects of Hubei Province of China (2014BBB010).

Ethical approval This article does not contain any studies with human participants or animals performed by any of the authors. I would like to declare on behalf of my co-authors that we have no conflict of interest.