Inducible Prophage Mutant of Escherichia coli Can Lyse New Host and the Key Sites of Receptor Recognition Identification

The use of bacteriophages as therapeutic agents is hindered by their narrow and specific host range, and by a lack of the knowledge concerning the molecular mechanism of receptor recognition. Two P2-like coliphages, named P88 and pro147, were induced from Escherichia coli strains K88 and DE147, respectively. A comparison of the genomes of these two and other P2-like coliphages obtained from GenBank showed that the tail fiber protein genes, which are the key genes for receptor recognition in other myoviridae phages, showed more diversity than the conserved lysin, replicase, and terminase genes. Firstly, replacing hypervariable region 2 (HR2: amino acids 716–746) of the tail fiber protein of P88 with that of pro147 changed the host range of P88. Then, replacing six amino acids in HR2 with the corresponding residues from pro147 altered the host range only in these mutants with changes at position 730 (leucine) and 744 (glutamic acid). Thus, we predicted that these amino acids are vital to establish the host range of P88. This study provided a vector of lysogenic bacteria that could be used to change or expand the phage host range of P88. These results illustrated that, in P2-like phage P88, the tail fiber protein determined the receptor recognition. Amino acids 716–746 and the amino acids at positions 730 and 744 were important for receptor recognition.

Introduction

Escherichia coli is an important pathogen in intestinal and extraintestinal infections (Johnson and Russo, 2002; Kaper et al., 2004; Kim, 2012; Croxen et al., 2013). Antibiotics are commonly used to treat colibacillosis (infection with a bacteria called E. coli); however, antibiotic resistance among bacteria, especially that of multi-drug-resistant strains (Kim et al., 2007; Hasan et al., 2011; Liu et al., 2016) have led to the reconsideration of bacteriophages as alternative therapeutic agents (Brussow, 2012; Keen, 2012; Young and Gill, 2015; Gorski et al., 2016). Compared with traditional antibiotic therapy, phages can target specific bacteria and reduce damage to the normal flora of the host (Clark and March, 2006). Unfortunately, the main drawback of phages as therapeutic agents for bacterial infection is their narrow host range (Haq et al., 2012).

In general, phages can undergo either a lytic or lysogenic lifecycle, and the genome of a lysogenic phage is integrated into the host chromosome and replicated passively as the host genome is replicated. Lysogenic phages can exist as phage particles and as part of the host genome. When in the host chromosome, the phage genome can be modified by manipulating the bacterial genome using molecular biological techniques (Posfai et al., 1999; Datsenko and Wanner, 2000; Yu et al., 2008; Blank et al., 2011).

In the Myoviridae adsorption process, such as that of phage T4, binding to cell membrane surfaces during the infection processes has two steps: first, specific host recognition occurs through a reversible interaction of the tail fibers with lipopolysaccharides or with the outer membrane porin protein C (Yu and Mizushima, 1982; Bartual et al., 2010); second, the tail spikes (also called short tail fibers) extend and irreversibly bind to the lipopolysaccharides (Riede, 1987; Thomassen et al., 2003). The host selectivity of a phage is determined by its tail fiber specificity, meanwhile tail spike specificity is as important for the infection process and host recognition as tail fiber specificity (Kageyama et al., 2009). Tail fiber proteins of many phages have been identified as important regions for receptor recognition determination experimentally (Riede et al., 1985; Montag et al., 1990; Yoichi et al., 2005; Mahichi et al., 2009; Bartual et al., 2010). In P2 like phages, only the tail spike has been proven to be the important region for the receptor recognition (Kageyama et al., 2009), while the involvement of the tail fiber protein has not been proven experimentally. Thus, in the present study, we isolate P2-like phages and determined whether the tail fiber protein was the receptor recognition region. The different parts of the tail fiber of one phage were replaced with the corresponding parts of another phage, and then the host ranges of the phage mutants were analyzed. We aimed to identify the specific binding sequences that determine the phage host range and to alter or expand the phage host range in the vector of lysogenic bacteria.

Materials and Methods

Bacterial Strains, Plasmids, and Growth Conditions

Strains used for prophage induction and phage isolation were listed in Supplementary Tables S1 and S2, respectively. Strain K88 came from a swine source. In addition, the 101 avian pathogenic E. coli strains used in these experiments were isolated from the brains of ducks with clinical signs of septicemia and neurological symptoms at different times and in different areas in the east of China, as previously described (Ma et al., 2013) (Supplementary Table S2). All E. coli strains were grown in Luria-Bertani (LB) medium at 37°C. When necessary, the LB medium was supplemented with appropriate antibiotic: ampicillin (100 μg/mL), kanamycin (50 μg/mL), chloramphenicol (30 μg/mL), or tetracycline hydrochloride (5 μg/mL). The plasmids used in this study were listed in Supplementary Table S3. The Red- and I-SceI-expressing plasmid pWRG99 was a derivative of pKD46 (Datsenko and Wanner, 2000; Blank et al., 2011).

Prophage Induction from the Chromosome of 54 E. coli Strains

Bacteriophages were obtained following mitomycin C (Sigma, St. Louis, MO, USA) induction of 54 lysogenic strains (Supplementary Table S1). Each strain was cultured for 5 h in LB medium and then diluted 1:100 in 5 mL of fresh medium. When the optical density at 600 nm (OD600) of the cultures reached 0.2, mitomycin C was added to a final concentration of 500 ng/mL. The cultures were then incubated at 37°C for 10 h. The lysates were centrifuged for 10 min at 3,500 × g at 4°C, and the supernatants were collected for follow-up experiments. The double-layer agar plate method (Adams, 1959) was used to find the host of the induced bacteriophage from the 108 E. coli strains (Supplementary Table S2).

Optimal Multiplicity of Infection (MOI) and One-Step Growth Curve of P88

Optimal MOI

Early log phase cells were infected with P88 at five different ratios (0.1, 0.01, 0.001, 0.0001, and 0.00001). After incubation for 3.5 h at 37°C, the phage lysates was centrifuged at 3,500 × g for 5 min. The supernatant was diluted to determine the phage titer (Lu et al., 2003).

One-Step Growth Curve

The phage and bacteria were mixed. Samples only containing absorbed phage were taken at 5 min intervals and then diluted to determine the plaque forming units (pfu) using the double-layer plate method (Chow et al., 1988; Ma and Lu, 2008; Hejnowicz et al., 2009).

Electron Microscopy

The phage filtrate was applied to a copper grid before negative staining with phosphotungstic acid (PTA, 2% w/v). Electron micrographs were observed using an H_7650 transmission electron microscope (TEM; Hitachi, Japan).

Phage DNA Extraction

Phage DNA extraction was performed as previously described (Ma and Lu, 2008) with some modifications. The phage pellet was digested with 2 μg/mL DNase and 10 μg/mL RNase for 30 min at 37°C. Phages were purified with NaCl-polyethylene glycol (PEG) 8000, and the DNA was isolated using SDS-Proteinase K. To analyze the quality of the DNA fragments, the molecules were separated by 0.8 % (w/v) agarose gel electrophoresis in TAE buffer (40 mM Tris-HCl, 500 mM sodium acetate, 50 mM EDTA, pH 7.2). The DNA was then suspended in ultrapure water.

Genome Sequencing

The lysogenic bacteria K88, host bacteria DE048, and phage P88 were analyzed by high-throughput sequencing to confirm if the phage P88 genome is located in the K88 genome. The phage P88 and pro147 genomes were sequenced by Roche 454 Sequencing, and reads were assembled into contigs using the Newbler (version 2.8). Strain K88 and the DE048 genomes were sequenced by Illumina MiSeq Sequencing. Reads were assembled into contigs using the Newbler.

Functional Analysis of Phage P88 and pro147

RAST (Aziz et al., 2008; Overbeek et al., 2014; Brettin et al., 2015) was used to predict the open reading frames (ORFs) of sequences and to analyze their corresponding functions. BLASTp was used to analyze the putative ORFs against the NCBI non-redundant proteins (NR) database. The phylogenetic trees of several selected genes were constructed with MEGA5.2 (Tamura et al., 2011) using the Neighbor-Joining algorithm. The structural model was generated using the Phyre¹ server¹ (Kelley et al., 2015).

Comparative Genome Analysis

Nucleic acid sequences of 10 P2-like phages [P88, pro147, P2 (NC_001895.1), Wphi, PsP3, 186, L-413C, Fels-2, fiAA91-ss, and P2 (KC618326.1)] were downloaded from GenBank. Comparisons of nucleic acid sequences were carried out with BLASTn. Comparisons of nucleic acid sequences of P88 (or pro147) with other P2-like phages were performed. Comparisons of nucleic acid sequences of P88 with pro147 were also performed. Phylogenetic trees based on the amino acid sequences (the lysins, the replicase, the tail fiber protein and the ATPase subunit of the terminase) were constructed with MEGA to analyze their evolutionary relationships.

Comparison of Tail Fiber Amino Acid Sequences

Amino acid sequences of the tail fiber proteins of 10 P2-like phages [P88, pro147, P2 (NC_001895.1), Wphi, PsP3, 186, L-413C, Fels-2, fiAA91-ss, and P2 (KC618326.1)] were downloaded from the NCBI database and comparisons of amino acid sequences were carried out with ClustalX (version 2.1). The phylogenetic tree of the tail fiber proteins showed that P88 had a close evolutionary relationship with pro147, P2 (NC_001895.1), and Wphi. The amino acid sequences of the tail fiber proteins of P88, pro147, P2 (NC_001895.1), and Wphi were then compared to analyze which segments of encoding DNA could be selected for gene manipulation. The amino acid sequences of the tail fiber proteins of P88 and pro147 were compared to select the replacement sections.

Construction of Gene Mutants to Analyze the Key Regions and Sites

The recombineering method of scarless mutagenesis was carried out according to the protocol described previously (Blank et al., 2011) with some modifications. The recombineering method of scarless mutagenesis was divided into two steps. In the first step, the target gene was replaced with a chloramphenicol resistance cassette (cat) gene and the I-SceI recognition site, which was PCR-amplified from plasmid pKD3 (Supplementary Table S3). The primers (Supplementary Table S3) were homologous to the flanking region of the corresponding target gene. The PCR products were then transformed by electroporation into K88 containing the lambda Red recombinase expression plasmid pKD46 (Datsenko and Wanner, 2000; Wang et al., 2011). After electroporation, samples were incubated at 30°C for 90 min in SOC broth (Datsenko and Wanner, 2000) and plated on LB agar with chloramphenicol to select for target gene mutants. Mutants were confirmed by PCR and sequencing using primers k1 and k2 (c1 and c2) (Supplementary Table S3). In the second step, the PCR products containing regions homologous to the flanking regions of the resistance gene for homologous recombination were transformed by electroporation into K88 containing the lambda Red recombinase expression and I-SceI endonuclease expression plasmid pWRG99 (Supplementary Table S3). After electroporation, samples were incubated at 30°C for 90 min in SOC broth and plated on LB agar with tetracycline hydrochloride. Selection of successful recombinants was mediated by the sequential expression of I-SceI endonuclease after Red-mediated recombination. If recombination occurred, the 200–500 bp dsDNA fragment harboring the same homologous region would replace the corresponding sequence of K88. If no recombination occurred, the unique I-SceI restriction site within the genome would still be present and the I-SceI endonuclease expressed from plasmid pWRG99 would kill unsuccessful recombinants by induction of double-strand breaks. Finally, pWRG99 was cured by streaking colonies and incubated at 37°C. Mutants were confirmed by PCR and sequencing using primers check1 and check2 (Supplementary Table S3).

Induction, Lytic Capacity Analysis, and Host Range Test of P88 Mutants

Each K88 mutant strain was cultured for 5 h in LB medium and then diluted 1:100 in 5 mL fresh medium. When the culture reached an OD600 of 0.2, mitomycin C was added and then was incubated at 30°C for 10 h. The lysates were centrifuged and used to test their lytic capacity by the double-layer agar plate method with DE048, the host of P88.

Phage P88 could not infect the bacteria if the phage suspension was dropped onto the LB agar plate spread with bacteria, so the host range analysis of each P88 mutant strain constructed by the scarless mutagenesis method was carried out by the double-layer agar plate method with DE048, MC1061, DH5a and BL21.

Results

Two Prophages of P88 and pro147 Were Induced from 54 E. coli Strains

Phage induction of 54 E. coli strains was investigated by the addition of mitomycin C. After induction, cell lysates were prepared and tested for their ability to lyse various E. coli hosts in plaque assays using the 108 strains listed in Supplementary Table S2. Out of the 54 E. coli strains used, the lysates of E. coli strain K88 and DE147 produced phages that could lyse the clinical isolates of avian pathogenic E. coli (APEC) strain DE048. These two phages were designated as P88 and pro147. After growth on DE048 in double agar LB plates at 37°C for 12 h, P88 formed clear and round plaques of approximately 0.06 cm in diameter (Figure 1A); pro147 formed clear and round plaques of approximately 0.13 cm in diameter (Figure 1B). Among the 108 E. coli strains tested, P88 specifically lysed DE048, while pro147 specifically lysed DE048, DH5a, BL21, and MC1061.

FIGURE 1

FIGURE 1. Fundamental characteristics of P2-like phages. (A) Plaques formed by phage P88. (B) Plaques formed by phage pro147. Growth of the phages (P88 and pro147) on Escherichia coli DE048 in double agar LB plates for 6 h. P88 formed clear and round plaques of approximately 0.06 cm in diameter. Pro147 formed clear and round plaques of approximately 0.13 cm in diameter. (C) One-step growth curve. Phage and bacteria were mixed, and samples containing only absorbed phage were taken at 5 min intervals and diluted. The double-layer plate method was used to determine the plaque forming units (pfu). One-step growth curve for P88 showed a latent period of about 20 min, a rise period of 25 min, and an average burst size of about 65. (D) Electron microscopy image of P88. (E) Electron microscopy image of pro147. Electron microscopy showed that the heads of P88 and pro147 were icosahedrons, and the tails was long and could shrink. For P88, the diameter of its capsid head was estimated at 50 nm and the tail length was estimated at 140 nm. For pro147, the diameter of its capsid head was estimated at 50 nm and the tail length was estimated at 130 nm.

Optimal MOI and One-Step Growth Curve of P88

After incubation for 3.5 h, P88 had the highest concentration (1.1 × 10¹⁰) at the ratio of 0.01 (pfu/cfu); therefore, the ratio of 0.01 (pfu/cfu) was considered as the optimal MOI. The one-step growth curve obtained for P88 showed a latent period of about 20 min, a rise period of 25 min, and an average burst size of about 65 (Figure 1C).

Electron Microscopy of P88 and pro147

Electron microscopy showed that the heads of P88 and pro147 were icosahedrons, and the tail was long and could shrink. In P88, the diameter of the capsid head was estimated at 55 nm and the tail length was estimated at 137 nm (Figure 1D). For pro147, the diameter of the capsid head was estimated at 54 nm and the tail length was estimated at 131 nm (Figure 1E). Phage P88 and pro147 belong to Myoviridae phages.

Genome Sequencing and Analysis

Genome Submission

P88 and pro147 genomes were submitted to GenBank with the accession numbers NC_026014.1 and NC_028896.1, respectively. The Whole Genome Shotgun project of DE048 and K88 were deposited at DDBJ/EMBL/GenBank under the accession numbers LBBM00000000 and LBBN00000000, respectively.

Overview of the P88 and pro147 Genomes

Comparative genome analysis showed that the P88 genome came from the E. coli K88 genome, and the mechanism of P88 prophage genome excision from E. coli K88 genome was displayed in Figure 2A. The P88 genome consisted of a double-stranded 35,814-bp DNA with 53 putative ORFs. The attP site of the P88 genome was determined as a 14-bp fragment (GCCACCCGAAGGTG). The pro147 genome consisted of a double-stranded 32,675-bp DNA with 44 putative ORFs. BLAST analysis of the P88 sequence at the NCBI database revealed that it had high similarity to prophages in E. coli, Shigella boydii, and Salmonella enterica, with more than 76 strains having coverage higher than 68% (ident 96%). However, P88 had low homology with reported P2-like phages and shared low similarity with P2, covering 5% (ident 90%).

FIGURE 2

FIGURE 2. Mechanism of P88 prophage genome excision from the E. coli K88 genome and the rationale of homologous recombination of P88. (A) Mechanism of P88 prophage genome excision from the E. coli K88 genome. Schematic representation of part of the lysogenic bacteria K88 genome and the location of the attL and attR regions. During excision, the P88 genome is circularized between the flanking host chromosomes. (B) The rationale of homologous recombination of P88. The phage genome is modified in the genome of lysogenic bacteria and then prophage production is induced. First, the target gene of P88 was selected for to manipulation, and then the segments of target gene were replaced with the corresponding segments of pro147 which were transformed by electroporation into K88. Finally, K88 mutant prophages were induced.

BLAST analysis of the pro147 sequence at the NCBI database revealed that it had high similarity to the prophages in E. coli, S. sonnei, Klebsiella pneumonia, and S. enterica, with more than 68 strains having coverage higher than 68% (ident 96%). Phage pro147 had high similarity to the reported P2-like phages [P2: coverage 75% (ident 97%); Wphi: coverage 75% (ident 98%); L-413C: coverage 73% (ident 96%); fiAA91-ss: coverage 73% (ident 98%)]. Although both P88 and pro147 could infect E. coli strain DE048, the P88 genome exhibited low similarity with pro147, covering 6% (ident 92%).

Functional Module Analysis

According to the putative functions of the ORFs predicted by RAST and BLASTp analyses, the P88 genome could be organized into six major functional modules involved in morphogenesis, replication, regulation, packaging, lysogeny, and lysis. Functional module and comparative genome analysis demonstrated that P88 and six other P2-like phages [P2 (NC_001895.1), Wphi, PsP3, 186, L-413C, and pro147] had similar functional modules (Figure 3).

FIGURE 3

FIGURE 3. Comparative genome alignments of the seven P2-like phages with similar modules. Modular organization of seven P2-like phage genomes [P88, pro147, P2 (NC_001895.1), Wphi, PsP3, 186 and L-413C]. Genes are grouped into six functional modules associated with the phage life cycle: lysogeny, replication, regulation, packaging, morphogenesis, and lysis. Sets of genes with similar function are grouped using the same color. The bottom line provides a base pair scale of the genomes.

Evolutionary Relationship Analysis

Phylogenetic trees based on the amino acid sequences (the lysins, the replicase, the tail fiber protein, and the ATPase subunit of the terminase) of the 10 P2 like phages [P88, pro147, P2 (NC_001895.1), Wphi, PsP3, 186, L-413C, Fels-2, fiAA91-ss, and P2 (KC618326.1)] were constructed. These 10 prophage genomes, which were P2-like phages whose hosts belong to the Enterobacteriaceae, were obtained from the NCBI database (Table 1).

TABLE 1

TABLE 1. Properties of P2-like phages whose hosts belong to the Enterobacteriaceae.

The phylogenetic trees based on the amino acid sequences of lysin, replicase and terminase (ATPase subunit) showed that phage P88 appeared to be phylogenetically distinct from the other P2-like phages, while, pro147 had close evolutionary relationships with P2 (NC_001895.1), Wphi, L-413C, and fiAA91-ss (Figures 4A,B,D). In the phylogenetic tree of the tail fiber protein, P88 had a close evolutionary relationship with pro147, P2 (NC_001895.1) and Wphi (Figure 4C). These results showed that the genes of tail fiber proteins were more diverse than the conserved genes of lysin, replicase and terminase.

FIGURE 4

FIGURE 4. Phylogenetic trees (the lysins, the replicase, the tail fiber protein, and the ATPase of the terminase) of the 10 P2-like phages whose hosts belong to the Enterobacteriaceae. Neighbor-joining tree analysis and bootstrap analysis (500 bootstrap Replications) based on the alignment of the amino acid sequence of the lysins (A), the replicase (B), the tail fiber protein (C), and the ATPase subunits (D) of the 10 P2-like phages [P88, pro147, P2 (NC_001895.1), Wphi, PsP3, 186, L-413C, Fels-2, fiAA91-ss, and P2 (KC618326.1)]. The numbers at the nodes indicate the bootstrap probabilities of that particular branch.

Comparison of Tail Fiber Amino Acid Sequences

The P88 genome exhibited low similarity with pro147, covering 6% (ident 92%); however, the tail fiber protein sequence of P88 shared high similarity with that of pro147, covering 100% (ident 70%). In addition, the P88 tail fiber protein shared high similarity with that of P2 (NC_001895.1), covering 100% (ident 65%) and Wphi, covering 100% (ident 62%). Alignment of the P88 tail fiber protein with those of pro147, P2 and Wphi showed that there were two hypervariable regions of 178aa (517–695 in P88) and 31 aa (716–746 in P88) in the C-terminal region (Figure 5A). Alignment of the P88 tail fiber protein with that of pro147 showed that there were two segments of hypervariable regions of 120aa (576–695 in P88) and 31 aa (716–746 in P88) in the C-terminal region (Figure 7A).

FIGURE 5

FIGURE 5. Alignment of the tail fiber structural proteins and structure prediction and modeling of the tail fiber proteins of P88 and pro147. (A) Alignment of the P88 tail fiber protein with similar regions (tail fiber proteins) from phages pro147, P2 (NC_001895.1) and Wphi. The tail fiber protein sequence of P88 shared high similarity with that of pro147, covering 100% (ident 70%); P2 (NC_001895.1), covering 100% (ident 65%), and Wphi, covering 100% (ident 62%). There were two hypervariable regions of 178aa (517–695 in P88) and 31 aa (716–746 in P88), respectively, in the C-terminal region. (B) Structure prediction and modeling of the tail fiber proteins of P88 and pro147. Sequence alignment of phage P88 tail fiber (or pro147) with its homologs from E. coli phage T4 gp34 were processed. (B1) The model crystal structure of amino acids 937–1127 of T4 gp34 (PDB code: c4uxeB.1). (B2) Homology modeling of P88 tail fiber based on the crystal structure of amino acids 937–1127 of T4 gp34 (PDB code: c4uxeB.1). The model was generated using the Phyre2 server (http://www.sbg.bio.ic.ac.uk/phyre2/html/). (B3) The model crystal structure of amino acids 984–1108 of T4 gp34 (PDB code: c4uxeB.1). (B4) Homology modeling of pro147 tail fiber based on the crystal structure of amino acids 984–1108 of T4 gp34 (PDB code: c4uxeB.1). The model was generated using the Phyre2 server (http://www.sbg.bio.ic.ac.uk/phyre2/html/).

Structure prediction and homology modeling of the tail fiber proteins of P88 and pro147 using the Phyre2 server² showed that the C-terminal domain (amino acids 426–612) of the P88 tail fiber (Figure 5B2) had the same folding topology as that of the domain (amino acids 937–1127) of T4 phage gp34 (Figure 5B1), with confidence of 97.17%; the C-terminal domain (amino acids 423–546) of the pro147 tail fiber (Figure 5B4) had the same folding topology as that of the domain (amino acids 984–1108) of T4 phage gp34 (Figure 5B3), with confidence of 97.17%. Despite the high degree of structural similarity, the tail fiber domains of P88 (or pro147) and T4 phage gp34 shared low amino acid identity over the relevant sequences with identity of 18 and 21%, respectively (Figure 5B).

Gene Mutants Were Constructed to Analyze the Key Regions and Sites

Construction of gene mutants was carried out by modifying the phage genome in the genome of lysogenic bacteria and then inducing the production of the prophage (Figure 2B). The scarless mutagenesis method (Figure 6) was used to produce mutants of P88. We modified the tail fiber protein of P88 by replacing different parts of it with the corresponding parts of pro147. After tail fiber comparison of P88 and pro147, hypervariable regions of amino acids 576–695 (HR1) and 716–746 (HR2) (Figure 7A) in the C-terminus of the P88 tail fiber protein were replaced with the corresponding parts of pro147. After replacement of HR2, the phage mutant P88HR2 displayed a different host range to that of P88; therefore the 12 different amino acids in the region of HR2 [amino acid sets SA (716), SB (718), SC (719), SD (721), SE (722), SF (725), SG (729), SH (730), SI (734), SJ (736), FC (744), and FB (746)] (Figure 7B) and amino acid set FA (757) (Figure 7B) of P88 were replaced with the corresponding amino acids of pro147. Finally, nine recombinant bacteria of K88HR1 (576–695), K88HR2 (716–746), K88SD (721), K88SE (722), K88SH (730), K88SI (734), K88FA (757), K88FB (746) and K88FC (744) were constructed (Table 2).

FIGURE 6

FIGURE 6. Rationale of the recombineering method of scarless mutagenesis. (A) chloramphenicol resistance cassette (cat) and an I-SceI recognition site were PCR amplified from plasmid pKD3 using primers with 40 bp homology extensions (red and green rectangles) and were integrated within a target gene “orfA” via Red-mediated recombination. The plasmid pWRG99 (containing the lambda Red recombinase and I-SceI endonuclease) and a 200–500 bp dsDNA homologous recombination fragment were transformed by electroporation sequentially. (A) Site-directed mutagenesis of the region adjacent to the I-SceI recognition site and the resistance cassette. With the help of Red-mediated recombination, homologous recombination occurred between the 200 and 500bp dsDNA fragment harboring the same homologous regions (red and green rectangles) and the corresponding sequence of K88. (B) The bacteria would be killed if homologous recombination did not occur. The I-SceI endonuclease expressed from plasmid pWRG99 would kill unsuccessful recombinants by induction of double-strand breaks.

FIGURE 7

FIGURE 7. Protein fragments and amino acids of tail fiber used for exchange mutagenesis. (A) Protein fragments of the tail fiber for exchange mutagenesis. Alignment of P88 tail fiber protein with that from phages pro147 showed that there were two hypervariable regions of 120aa (576–695 in P88) and 31 aa (716–746 in P88) in the C-terminal region. Hypervariable region of 120aa (HR1) and the region of 31 aa (HR2) of P88 were replaced with the corresponding parts of pro147. (B) Amino acids of the tail fiber for exchange mutagenesis. Twelve different amino acids in the HR2 region (amino acid sets 716, 718, 719, 721, 722, 725, 729, 730, 734, 736, 744, and 746) and amino acid 757 of P88 were replaced with the corresponding amino acids of pro147.

TABLE 2

TABLE 2. Results of the mutation experiment using scarless mutagenesis.

Induction and Lytic Capacity Analysis of P88 Mutants

The mutant bacteria were then induced by mitomycin C to evaluate their lysis capacity. Nine mutants of P88 were investigated. Eight mutants of K88HR1(576–695), K88HR2 (716–746), K88SD (721), K88SE (722), K88SH (730), K88SI (734), K88FA (757), and K88FB (746) could be induced to produce phages (P88HR1, P88HR2, P88SD, P88SE, P88SH, P88SI, P88FA, and P88FB) those were able to infect DE048 (Table 2; Figure 8), the host of P88; however, the bacteria mutant K88FC (744: glutamic acid to proline) could not be induced to produce phage that infected DE048 (Table 2). Meanwhile, K88FC could not be induced to produce phage to infect MC1061, DH5a, and BL21. The phage particles induced from the mutants that could lyse the E. coli strain DE048 were enriched for the subsequent experiment of host range determination and electron microscopy.

FIGURE 8

FIGURE 8. Host range of P88, pro147 and P88 mutants (P88HR1, P88HR2, P88SD, P88SE, P88SH, P88SI, P88FA, and P88FB). The host ranges of P88, pro147, and P88 mutants (P88HR1, P88HR2, P88SD, P88SE, P88SH, P88SI, P88FA, and P88FB) were displayed. On the lawn of DE048, all the phages could form clear plaques; on the lawn of MC1061 and DH5a, phages pro147, P88HR2, and P88SH could form clear plaques; on the lawn of BL21, phages pro147, and P88HR2 could form clear plaques. Phages pro147 and P88HR2 had the same host ranges.

The Host Range Test of P88 Mutants

The host ranges of P88, pro147 and P88 mutants (P88HR1, P88HR2, P88SD, P88SE, P88SH, P88SI, P88FA, and P88FB) were analyzed. After incubation at 37°C for 10 h, on the lawn of DE048, all the phages could form clear plaques; on the lawn of MC1061 and DH5a, phages pro147, P88HR2 and P88SH (730: leucine to isoleucine) could form clear plaques; on the lawn of BL21, phages pro147 and P88HR2 could form clear plaques; pro147 and P88HR2 had the same host ranges (Figure 8). P88, pro147 and all the P88 mutants could not infect DE048 when incubated at 28°C.

Electron micrographs showed that compared with those of the original P88, the head and the tail of the phage mutants with altered host ranges (P88HR2 and P88SH) and with no change in host range (P88HR1, P88SE, and P88FB) did not have visible differences (Figure 8).

Discussion

We screened two P2-like coliphages, named P88 and pro147, and identified their corresponding host bacteria. Analysis of the key regions and sites that might determine receptor recognition showed that the host range of P88 was altered after replacing the hypervariable regions HR2 (amino acids 716–746) of the tail fiber protein of P88 with that of pro147; the host ranges of P88 with replacement of amino acid 730 (leucine) and 744 (glutamic acid) were changed after replacing six amino acids in HR2 with the corresponding amino acids of pro147. These results showed that, in P88, the tail fiber determined the receptor recognition; the region of amino acids 716–746 was important for receptor recognition; amino acids 730 (leucine) and 744 (glutamic acid) were the key amino acids of P88 that determine the host range.

Replacing regions of the tail fiber protein in P88 with those of pro147 changed the host range, which indicated that the corresponding regions of pro147 would determine its receptor recognition. Furthermore, the high homology of the tail fiber sequences in P2-like phages P88, pro147, P2 (NC_001895.1), and Wphi (005056.1) caused us to speculate that the corresponding parts in P2 and Wphi would also determine their receptor recognition specificity.

For temperate phages, when the phage genome is modified by manipulating the bacterial genome, the mutant strains constructed could be confirmed by PCR and the prophage could then be introduced and characterized. If the phage genome was modified in the host cell when it was separated from the bacterial chromosome by homologous recombination with a homologous fragment in the form of plasmid or DNA fragment, like lytic phages (Yoichi et al., 2005; Mahichi et al., 2009; Thomason et al., 2009), a plaque assay was used to screen the mutant phages; however, the recombinant phage would not be selected without a sensitive host strain. Hence, to construct mutants of temperate phages, modifying the phage genome by manipulating the bacterial genome was preferred.

The lysogeny module is required to form stable lysogens (Court et al., 2007). It is thought that a phage should be lytic and not lysogenic for use as an antibacterial agent (Chan et al., 2013). Here, the phage modification experiment was conducted on lysogenic bacteria (phage in the lysogenic state). In other words, the induced phage was temperate and able to perform lysogenic infections. Therefore, to disrupt lysogenic infection, the genes of lysogeny module could be deleted without disrupting the lytic ability of the phage (Zhang et al., 2013). Thus, compared with wild temperate phage strains, the novel modified phages are expected to minimize potential hazards to be integrated into the bacteria in normal flora if used as antibacterial agents.

In the next step of our research, we could create a homologous sequence library that contains sufficient allelic sequences to modify specific binding sequences by gene synthesis and homologous sequences enriched from natural samples. The crystal structure of the tail fiber protein with the key sites for receptor binding would be analyzed to understand how P2-like phages recognize their host.

Conclusion

This study located the specific binding sequences that determined the phage host range of P88 and increased our understanding of the function of the tail fibers in P2-like phages. The result revealed that lysogenic bacteria could be used as vectors to change or broaden a phage’s host range. It might also be possible to modify the genomes to analyze gene function or change/expand the host range of other bacteriophages using this method.

Author Contributions

WZ and MC conceived the study and designed experiments. MC, HY, and CL performed the data analysis. MC, LZ, and SX carried out experiments. MC and LZ analyzed experimental results. MC and WZ wrote the manuscript.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

This work was supported by grants from National Natural Science Foundation of China (31322054), Priority Academic Program Development of Jiangsu Higher Education Institutions and the Special Fund for Public Welfare Industry of Chinese Ministry of Agriculture (201303041).

Supplementary Material

The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fmicb.2017.00147/full#supplementary-material

Footnotes

1. ^http://www.sbg.bio.ic.ac.uk/phyre2/html/

References

Adams, M. H. (1959). Bacteriophages. New York: Interscience.

Google Scholar

Aziz, R. K., Bartels, D., Best, A. A., Dejongh, M., Disz, T., Edwards, R. A., et al. (2008). The RAST server: rapid annotations using subsystems technology. BMC Genomics 9:75 doi: 10.1186/1471-2164-9-75

PubMed Abstract | CrossRef Full Text | Google Scholar

Bartual, S. G., Otero, J. M., Garcia-Doval, C., Llamas-Saiz, A. L., Kahn, R., Fox, G. C., et al. (2010). Structure of the bacteriophage T4 long tail fiber receptor-binding tip. Proc. Natl. Acad. Sci. U.S.A. 107, 20287–20292. doi: 10.1073/pnas.1011218107

PubMed Abstract | CrossRef Full Text | Google Scholar

Blank, K., Hensel, M., and Gerlach, R. G. (2011). Rapid and highly efficient method for scarless mutagenesis within the Salmonella enterica chromosome. PLoS ONE 6:e15763. doi: 10.1371/journal.pone.0015763

PubMed Abstract | CrossRef Full Text | Google Scholar

Brettin, T., Davis, J. J., Disz, T., Edwards, R. A., Gerdes, S., Olsen, G. J., et al. (2015). RASTtk: a modular and extensible implementation of the RAST algorithm for building custom annotation pipelines and annotating batches of genomes. Sci. Rep. 5:8365 doi: 10.1038/srep08365

PubMed Abstract | CrossRef Full Text | Google Scholar

Brussow, H. (2012). What is needed for phage therapy to become a reality in Western medicine? Virology 434, 138–142. doi: 10.1016/j.virol.2012.09.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Chan, B. K., Abedon, S. T., and Loc-Carrillo, C. (2013). Phage cocktails and the future of phage therapy. Future Microbiol. 8, 769–783. doi: 10.2217/fmb.13.47

PubMed Abstract | CrossRef Full Text | Google Scholar

Chow, J. J., Batt, C. A., and Sinskey, A. J. (1988). Characterization of Lactobacillus bulgaricus bacteriophage ch2. Appl. Environ. Microbiol. 54, 1138–1142.

PubMed Abstract | Google Scholar

Clark, J. R., and March, J. B. (2006). Bacteriophages and biotechnology: vaccines, gene therapy and antibacterials. Trends Biotechnol. 24, 212–218. doi: 10.1016/j.tibtech.2006.03.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Court, D. L., Oppertheim, A. B., and Adhya, S. L. (2007). A new look at bacteriophage lambda genetic networks. J. Bacteriol. 189, 298–304. doi: 10.1128/JB.01215-06

PubMed Abstract | CrossRef Full Text | Google Scholar

Croxen, M. A., Law, R. J., Scholz, R., Keeney, K. M., Wlodarska, M., and Finlay, B. B. (2013). Recent advances in understanding enteric pathogenic Escherichia coli. Clin. Microbiol. Rev. 26, 822–880. doi: 10.1128/CMR.00022-13

PubMed Abstract | CrossRef Full Text | Google Scholar

Datsenko, K. A., and Wanner, B. L. (2000). One-step inactivation of chromosomal genes in Escherichia coli K-12 using PCR products. Proc. Natl. Acad. Sci. U.S.A. 97, 6640–6645. doi: 10.1073/pnas.120163297

PubMed Abstract | CrossRef Full Text | Google Scholar

Gorski, A., Miedzybrodzki, R., Weber-Dabrowska, B., Fortuna, W., Letkiewicz, S., Rogoz, P., et al. (2016). Phage therapy: combating infections with potential for evolving from merely a treatment for complications to targeting diseases. Front. Microbiol. 7:1515. doi: 10.3389/fmicb.2016.01515

PubMed Abstract | CrossRef Full Text | Google Scholar

Haq, I. U., Chaudhry, W. N., Akhtar, M. N., Andleeb, S., and Qadri, I. (2012). Bacteriophages and their implications on future biotechnology: a review. Virol. J. 9:9 doi: 10.1186/1743-422X-9-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Hasan, B., Faruque, R., Drobni, M., Waldenstrom, J., Sadique, A., Ahmed, K. U., et al. (2011). High prevalence of antibiotic resistance in pathogenic Escherichia coli from large- and small-scale poultry farms in Bangladesh. Avian Dis. 55, 689–692. doi: 10.1637/9686-021411-Reg.1

PubMed Abstract | CrossRef Full Text | Google Scholar

Hejnowicz, M. S., Golebiewski, M., and Bardowski, J. (2009). Analysis of the complete genome sequence of the lactococcal bacteriophage bIBB29. Int. J. Food Microbiol. 131, 52–61. doi: 10.1016/j.ijfoodmicro.2008.06.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Johnson, J. R., and Russo, T. A. (2002). Uropathogenic Escherichia coli as agents of diverse non-urinary tract extraintestinal infections. J. Infect. Dis. 186, 859–864. doi: 10.1086/342490

PubMed Abstract | CrossRef Full Text | Google Scholar

Kageyama, Y., Murayama, M., Onodera, T., Yamada, S., Fukada, H., Kudou, M., et al. (2009). Observation of the membrane binding activity and domain structure of gpV, which comprises the tail spike of bacteriophage P2. Biochemistry 48, 10129–10135. doi: 10.1021/bi900928n

PubMed Abstract | CrossRef Full Text | Google Scholar

Kaper, J. B., Nataro, J. P., and Mobley, H. L. T. (2004). Pathogenic Escherichia coli. Nat. Rev. Microbiol. 2, 123–140. doi: 10.1038/nrmicro818

PubMed Abstract | CrossRef Full Text | Google Scholar

Keen, E. C. (2012). Phage therapy: concept to cure. Front. Microbiol. 3:238 doi: 10.3389/fmicb.2012.00238

PubMed Abstract | CrossRef Full Text | Google Scholar

Kelley, L. A., Mezulis, S., Yates, C. M., Wass, M. N., and Sternberg, M. J. (2015). The Phyre2 web portal for protein modeling, prediction and analysis. Nat. Protoc. 10, 845–858. doi: 10.1038/nprot.2015.053

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, K. S. (2012). Current concepts on the pathogenesis of Escherichia coli meningitis: implications for therapy and prevention. Curr. Opin. Infect. Dis. 25, 273–278. doi: 10.1097/QCO.0b013e3283521eb0

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, T. E., Jeong, Y. W., Cho, S. H., Kim, S. J., and Kwon, H. J. (2007). Chronological study of antibiotic resistances and their relevant genes in Korean avian pathogenic Escherichia coli isolates. J. Clin. Microbiol. 45, 3309–3315. doi: 10.1128/JCM.01922-06

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Y. Y., Wang, Y., Walsh, T. R., Yi, L. X., Zhang, R., Spencer, J., et al. (2016). Emergence of plasmid-mediated colistin resistance mechanism MCR-1 in animals and human beings in China: a microbiological and molecular biological study. Lancet Infect. Dis. 16, 161–168. doi: 10.1016/S1473-3099(15)00424-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Lu, Z., Breidt, F., Fleming, H. P., Altermann, E., and Klaenhammer, T. R. (2003). Isolation and characterization of a Lactobacillus plantarum bacteriophage, Phi JL-1, from a cucumber fermentation. Int. J. Food Microbiol. 84, 225–235.

Google Scholar

Ma, J., Sun, M., Bao, Y., Pan, Z., Zhang, W., Lu, C., et al. (2013). Genetic diversity and features analysis of type VI secretion systems loci in avian pathogenic Escherichia coli by wide genomic scanning. Infect. Genet. Evol. 20, 454–464. doi: 10.1016/j.meegid.2013.09.031

PubMed Abstract | CrossRef Full Text | Google Scholar

Ma, Y. L., and Lu, C. P. (2008). Isolation and identification of a bacteriophage capable of infecting Streptococcus suis type 2 strains. Vet. Microbiol. 132, 340–347. doi: 10.1016/j.vetmic.2008.05.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Mahichi, F., Synnott, A. J., Yamamichi, K., Osada, T., and Tanji, Y. (2009). Site-specific recombination of T2 phage using IP008 long tail fiber genes provides a targeted method for expanding host range while retaining lytic activity. FEMS Microbiol. Lett. 295, 211–217. doi: 10.1111/j.1574-6968.2009.01588.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Montag, D., Hashemolhosseini, S., and Henning, U. (1990). Receptor-recognizing proteins of T-even type bacteriophages. The receptor-recognizing area of proteins 37 of phages T4 TuIa and TuIb. J. Mol. Biol. 216, 327–334.

Google Scholar

Overbeek, R., Olson, R., Pusch, G. D., Olsen, G. J., Davis, J. J., Disz, T., et al. (2014). The SEED and the rapid annotation of microbial genomes using subsystems technology (RAST). Nucleic Acids Res. 42, D206–D214. doi: 10.1093/nar/gkt1226

PubMed Abstract | CrossRef Full Text | Google Scholar

Posfai, G., Kolisnychenko, V., Bereczki, Z., and Blattner, F. R. (1999). Markerless gene replacement in Escherichia coli stimulated by a double-strand break in the chromosome. Nucleic Acids Res. 27, 4409–4415. doi: 10.1093/nar/27.22.4409

PubMed Abstract | CrossRef Full Text | Google Scholar

Riede, I. (1987). Receptor specificity of the short tail fibres (gp12) of T-even type Escherichia coli phages. Mol. Gen. Genet. 206, 110–115. doi: 10.1007/BF00326544

PubMed Abstract | CrossRef Full Text | Google Scholar

Riede, I., Degen, M., and Henning, U. (1985). The receptor specificity of bacteriophages can be determined by a tail fiber modifying protein. EMBO J. 4, 2343–2346.

PubMed Abstract | Google Scholar

Tamura, K., Peterson, D., Peterson, N., Stecher, G., Nei, M., and Kumar, S. (2011). MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol. Biol. Evol. 28, 2731–2739. doi: 10.1093/molbev/msr121

PubMed Abstract | CrossRef Full Text | Google Scholar

Thomason, L. C., Oppenheim, A. B., and Court, D. L. (2009). Modifying bacteriophage lambda with recombineering. Methods Mol. Biol. 501, 239–251. doi: 10.1007/978-1-60327-164-6_21

CrossRef Full Text | Google Scholar

Thomassen, E., Gielen, G., Schutz, M., Schoehn, G., Abrahams, J. P., Miller, S., et al. (2003). The structure of the receptor-binding domain of the bacteriophage T4 short tail fibre reveals a knitted trimeric metal-binding fold. J. Mol. Biol. 331, 361–373. doi: 10.1016/S0022-2836(03)00755-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, S., Niu, C., Shi, Z., Xia, Y., Yaqoob, M., Dai, J., et al. (2011). Effects of ibeA deletion on virulence and biofilm formation of avian pathogenic Escherichia coli. Infect. Immunol. 79, 279–287. doi: 10.1128/IAI.00821-10

PubMed Abstract | CrossRef Full Text | Google Scholar

Yoichi, M., Abe, M., Miyanaga, K., Unno, H., and Tanji, Y. (2005). Alteration of tail fiber protein gp38 enables T2 phage to infect Escherichia coli O157:H7. J. Biotechnol. 115, 101–107. doi: 10.1016/j.jbiotec.2004.08.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Young, R., and Gill, J. J. (2015). Phage therapy redux - What is to be done? Science 350, 1163–1164. doi: 10.1126/science.aad6791

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, B. J., Kang, K. H., Lee, J. H., Sung, B. H., Kim, M. S., and Kim, S. C. (2008). Rapid and efficient construction of markerless deletions in the Escherichia coli genome. Nucleic Acids Res. 36:e84 doi: 10.1093/nar/gkn359

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, F., and Mizushima, S. (1982). Roles of lipopolysaccharide and outer membrane protein OmpC of Escherichia coli K-12 in the receptor function for bacteriophage T4. J. Bacteriol. 151, 718–722.

PubMed Abstract | Google Scholar

Zhang, H., Fouts, D. E., Depew, J., and Stevens, R. H. (2013). Genetic modifications to temperate Enterococcus faecalis phage Ef11 that abolish the establishment of lysogeny and sensitivity to repressor, and increase host range and productivity of lytic infection. Microbiology 159, 1023–1035. doi: 10.1099/mic.0.067116-0

PubMed Abstract | CrossRef Full Text | Google Scholar