Characterization and complete genome of the virulent Myoviridae phage JD007 active against a variety of Staphylococcus aureus isolates from different hospitals in Shanghai, China

The implementation of phage therapy is re-emerging with the increase in widespread antibiotic-resistant bacteria. Staphylococcus phage JD007 was characterized and its complete genome sequence analysed. Staphylococcus phage JD007 was classified as belonging to the Myoviridae family based on its morphology, as observed by transmission electron microscopy. Its lytic activity was stable between pH 5–11 and below 42 °C; moreover, an absorbance curve showed that nearly 90% of the viral particles had adsorbed to its host after a 20 min co-incubation. The complete genome size is 141,836 bp, making JD007 one of the largest Staphylococcus phages of Myoviridae. No identifiable resistance or virulence genes were found in the JD007 genome. JD007 was able to lyse 95% of S. aureus isolates, including the prevalent ST239-MRSA and ST59-MRSA strains isolated from different hospitals in Shanghai, China, and inhibition assays showed that JD007 could inhibit S. aureus growth at a multiplicity of infection of 0.1. The results suggested that Staphylococcus phage JD007 can potentially be used in phage therapy or for the detection of S. aureus.


Background
Staphylococcus aureus is considered a commensal and major human pathogen responsible for a variety of acute and chronic diseases [1]. Due to the broad range of methicillin-resistant S. aureus (MRSA) strains and the emergence of vancomycin-resistant S. aureus (VRSA) strains [2], it is feasible that in the future, there will be a lack of antibiotics available to treat antibiotic-resistant infectious diseases. Accordingly, there is recent renewed interest in phage therapy. Indeed, phages that infect S. aureus have been used in clinical trials to treat chronic venous leg ulcers, and the results show no safety concerns with regard to the use of bacteriophage treatment [3]. Therefore, Staphylococcus phages can potentially be used for phage therapy.
Staphylococcus phages belonging to the Myoviridae family have been isolated worldwide [4][5][6][7][8]. The phages Romulus and Remus are able to infect approximately 70% of the tested S. aureus isolates and display lytic activity inside these hosts. Furthermore, both phages exhibited rapid initial adsorption and biofilm-degrading capacity [9]. In one study, phage ΦSA039 produced clear plaques on 13 of 15 Staphylococcus isolates (87%) [10], and phage Stau2 lysed 80% of the S. aureus isolates (164/205) obtained from hospitals in Taiwan, China [11]. Furthermore, phage ISP successfully infected 86% (31/ 34) of Staphylococcus isolates in another study, including relevant MRSA strains [12], and phage MSA6 infects a wide spectrum of staphylococcal strains originating from both humans and bovines [13]. In summary, these studies suggest that phages belonging to Myoviridae can kill a broad range of S. aureus.
In this paper, a Staphylococcus phage named JD007 was isolated from chicken faeces in Shanghai, China. Its morphology was assessed using transmission electron microscopy, its thermal and pH stability were evaluated, the inhibition assays and host adsorption rate were characterized, and the complete genome was sequenced and analysed. Furthermore, the host range of JD007 was characterized using prevalent strains of S. aureus isolated from different hospitals in Shanghai, China, to characterize JD007 in vitro bactericidal activity.

Staphylococcus phage JD007 morphology
A photograph of the bacteriophage JD007 was obtained by transmission electron microscopy, as shown in Fig. 1. Its head is approximately 60 nm in diameter, the length of its tail is approximately 110 nm, and the contractile sheath can be observed between the head and tail. Based on these morphological features, the phage belongs to Myoviridae. S. aureus phages belonging to Myoviridae have been isolated worldwide, including phages K, SA5, A5W, Sb-1, ISP, G1, GH15, SA11, Staphy1N, MSA6, 676Z, P4W, Fi200w, vB_SauM_Remus, vB_SauM_Romulus, Twort, phiIPLA-RODI, and phiIPLA-C1C [4][5][6][7][8]. Several studies have shown that these phages can prevent infections in animal models. JD007 was isolated from chicken faeces in Shanghai, China. These Myoviridae phages have genome sizes between 120 and 140 kb, exhibit a G + C content of 27.98 to 30.60%, and encode 183 to 217 open reading frames (ORFs).

Biophysical stability infection parameters
The stability of phage JD007 was also evaluated. Lytic activity was stable between pH5 and 11 ( Fig. 2a), whereas activity was completely lost when kept below pH 4 for 2 h at room temperature. As shown in Fig. 2b, the activity of phage JD007 was stable between 37 and 42°C, however, when kept at 50°C for 1 h, activity decreased by 25% compared to 37°C. Furthermore, JD007 lost nearly all activity at a temperature above 60°C. Sangeeta et al. also have reported a lytic phage of Salmonella enterica serovar Paratyphi B that could keep viral activity stable between pH 4 and 9 or between 4 and 40°C [14]. Based on thermal and acid stability tests, phage JD007 remained stable at room temperature and between pH 5 and 11.
The adsorption rate of Staphylococcus phage JD007 We also sought to characterize the adsorption rate of phage JD007 to its host cell, as shown in Fig. 2d. Phage JD007 adsorbed to its host cell S. aureus Sa60 80, 90, and 97% of the time when incubated at 37°C for 15, 20, and 30 min, respectively. To identify the in vitro bactericidal activity of phage JD007, S. aureus Sa60 was cultured at an OD 600 of 0.2 and then infected with phage JD007 at a multiplicity of infection (MOI) of 0, 0.1 and 1 at 37°C with shaking at 100 rpm. As shown in Fig. 2c, phage JD007 inhibited the growth of S. aureus after co-culture for 150 min, whereas S. aureus Sa60 was able to grow normally without phage JD007. Our results suggest that Staphylococcus phage JD007 can be used for inhibition the growth of S. aureus. Gutierrez reported that Staphylococcus phages phiIPLA-RODI and phiIPLA-C1C had the same lytic activity on S. aureus IPLA16, since no viable bacteria were detected after 8 h of incubation and a considerable decrease of the bacterial population was already achieved after 6 h of treatment [15].

The host range of Staphylococcus phage JD007
A total of 175 S. aureus isolates from different hospitals in Shanghai, China, were chosen to identify the host range of Staphylococcus phage JD007. Our results showed that phage JD007 could kill 95% (166/175) of S. aureus isolates of different types. As shown in Table 1, 41 of the 175 isolates representing different types were identified and characterized through MLS and Spa typing. Staphylococcus phage JD007 formed an inhibition zone in a double-layer plate, some of which were clear and some faint. Staphylococcus phage JD007 formed clear inhibition zones using ST239-t030 that were isolated from different hospitals in Shanghai, China. Importantly, phage JD007 killed the prevalent strains of ST239-t030 that account for an increase in the incidence of infectious diseases [16]. Phage JD007 also killed other types of isolates, including ST1-t127, ST1128-t164, ST188-t189, and ST239-t459, with high efficiency. Furthermore, phage JD007 killed the majority of ST239-t030 and ST59-MRSA strains, with different efficiencies, though several strains of ST239-t030 were absolutely not susceptible to phage JD007. We suggest that JD007 has a wide host range with different lytic efficiencies. Furthermore, no evidence for a significant correlation between phage activity and antimicrobial susceptibility, ST, or Spa type was found because phage JD007 killed different strain types with different antimicrobial susceptibilities at variable efficiencies. b The x-axis represents different temperatures, and the y-axis represents the relative titres of Staphylococcus phage JD007 compared to that at 37°C. % stability = (N/N 0 ) × 100, where N is the number of viable phages after 1 h of incubation and N 0 is the initial number of phages. c Inhibition assays of Staphylococcus phage JD007. The x-axis represents the co-culture time of Staphylococcus phage JD007 and Sa60, and the y-axis represents the change in OD 600 in the mixture of Staphylococcus phage JD007 infecting its host at different MOIs. The red, blue and green circles represent MOI = 0, 0.1 and 1, respectively. d The adsorption rate of Staphylococcus phage JD007 to its host S. aureus Sa60, with the x-axis representing the incubation time of Staphylococcus phage JD007 and its host, and the y-axis represents the percent of titres of phage that had not been absorbed to the host. The error bars represent the s.d

Safety assessment of Staphylococcus phage JD007 based on the genome sequence
The complete genome of Staphylococcus phage JD007 is 141,836 bp, and contains 217 ORFs [17]. The entire genome structure of Staphylococcus phage JD007 is shown in Fig. 3; the arrows represent the predicted ORFs consisting of genes involved in the bacteriophage structure and in DNA replication as well as other predicted functions. The proteins of helicase and major capsid protein were separately used to draw phylogenetic trees, as showed in Fig. 4a and Fig. 4b; phage JD007 and other Staphylococcus phages belong to Myoviridae were clustered in the same branch of tree, these results were constant with the results of its morphology. The phages chosen for phylogenic analysis are all belonged to Myoviridae, their host are different: phage A511 infects Listeria monocytogenes [18], while others infect S. aureus. Phage GH15 and JD007 were separately isolated from Jilin province and Shanghai, China, there are more than 1000 km between these two places, these two phages were clustered in the same subbranch, it meant their closest genetic relationship. All of the ORFs predicted by Glimmer were searched by blast in the database of toxin and antibiotic resistance genes, with no identifiable toxin or antibiotic resistance genes found in the genome until present. While, phage genomes are replete with genes of unknown function and it is possible that the genome of this phage contains toxin or resistance genes, they may be just of types that have not yet been characterized until present, so it is necessary to systematically study these genes' function for fully understanding safety of the phage for therapy. The complete genomic comparisons of phage JD007 with phage K, Twort and GH15 were showed in Fig. 5, its coverage of GH15 was 93% with identity of 97%, while its coverage of K was 85% with identity of 98%, and its coverage of Twort was 74% with identity of 8%. The results indicated the mosaic genomic structure among these bacterophages.
Bacteriophages can be considered vectors that transmit antibiotic resistance genes or virulence genes amongst bacteria. The prophage may contribute the pathogenic traits of Enterococcus faecalis [19]. Antibiotic resistance genes, such as blaTEM, qnrA, and blaCTX-M-1, as well as toxin genes, have been found in the genomes of bacteriophages isolated from human faecal samples [20,21]. Thus, the predicted ORFs of Staphylococcus phage JD007 were searched using blast against the database of virulence factors and antibiotic resistance genes, and no virulence or antibiotic resistance genes were found until present. This indicates that Staphylococcus phage JD007 will likely not act as a vector of dissemination of such undesired genes amongst bacteria when used for the prevention or control of S. aureus. (The "+" represent the activity of phage JD007 to its correspondence hosts, the more "+" the more clear inhibition zones were formed by the phage JD007)

Conclusion
Staphylococcus phage JD007 is stable at room temperature and under acidic conditions. It can quickly adsorb to its host cell and propagate in that host cell and has a wide host range. In addition, no antibiotic resistance or toxin genes were found in its genome until present. The phage kills 95% of S. aureus strains of different types, and most importantly, it can kill ST239-MRSA III-Spa t030 strains, which are prevalent and accounts for the increase in infectious disease incidence in China. JD007 may be potentially useful for phage therapy or the detection of S. aureus clinically.

Bacteria isolates and culture conditions
A total of 175 isolates of S. aureus were obtained from Ruijin Hospital, the Sixth People's Hospital of Shanghai, the Armed Police General Hospital, and the Centre Hospital of Changning District in Shanghai, China. The strains isolated from Ruijin Hospital were kindly provided by Qingtian Li, The use of these isolates in this paper was approved by these individuals, the ownership of the strains individually belonged to them, and we were authorized to use these isolates in this paper. The isolates were grown in liquid LB (Luria-Bertani) medium at 37°C, on solid LB medium (1.5% agar), or in LB soft agar overlays (0.7% agar). Phage JD007 was isolated from chicken faeces collected from the chicken slaughter facility in the Madang food market located on No.349-2# Madang Road in the Huangpu District of Shanghai, China [17]. S. aureus strain Sa60 isolated from Ruijin Hospital was used for phage JD007 amplification and the following experiments. Forty-one strains of the 175 total strains represent different types of S. aureus identified and characterized using MLS and Spa typing methods, as described previously [22]. Briefly, the strains were spa-typed via the online database (http://www.spaserver.ridom.de/). The sequence type (ST) was characterized by multi-locus sequence typing (MLST), and the products of seven housekeeping gene fragments were sequenced (Sangon Biotech, Shanghai) and compared to allele profiles from the database of S. aureus (http://saureus.mlst.net/). All of the bacteria were identified using the VITEK2 compact system, and antimicrobial susceptibility testing was performed by the disk diffusion method according to Clinical and Laboratory Standards Institute guidelines (M100S, 26 th edition) for penicillin (P) (10 units), cefoxitin (30 μg), gentamicin (CN) (10 μg), erythromycin (E) (15 μg), teicoplanin (30 μg), clindamycin (CC) (2 μg), rifampicin (5 μg), and linezolid (30 μg). The minimum inhibitory concentration (MIC) of vancomycin was determined using an E-test. S. aureus ATCC25923 and ATCC29213 were used as quality controls for the disk diffusion test and MIC detection, respectively.

Bacteriophage amplification and purification
High-titre phage stocks were obtained through amplification in liquid LB medium containing 10 μM MgCl 2 and 5 μM CaCl 2 . First, phage JD007 that infected S. aureus Sa60 cells at an MOI of 0.1 was incubated at 37°C overnight. Visible lysis of the liquid culture was obtained, and the lysate was then incubated with chloroform (final concentration of 2%) for 30 min with gentle shaking to kill the remaining bacteria. Bacteria debris was removed by centrifugation at 6,500 rpm (Beckman, JA18.0, USA) for 15 min. The phages in the supernatant were enriched at 4°C overnight using polyethylene glycol (PEG) 8000 and precipitated (final concentration 10%w/v) at 8,500 rpm for 20 min (Beckman, JA18.0). The pellet was dissolved in TM buffer and vortexed. PEG8000 was removed after adding the same volume of chloroform after vortexing, and the solution was centrifuged at 4,000 × g for 10 min, the supernatant contained a high concentration of phage. CsCl was added at a concentration of 0.5 g per 1 mL, and the phages were purified by discontinuous centrifugation through a CsCl gradient (1.33, 1.45, 1.50, and 1.70 g/cm 3 ) in TM buffer in Ultra-Clear tubes (Beckman Coulter, Inc., Fullerton, CA) at 120,000 × g for 4 h. The band of enriched phages was removed using a syringe, and the same was dialyzed against TM buffer and stored at 4°C.

Electron microscopic imaging
Using the purified phage obtained above. Then, phage particles were collected by centrifugation at 33,000 × g for 1 h and washed twice in 0.1 × PBS (pH7.4) using a Beckman high-speed centrifuge and a JA-18.1 fixedangle rotor. Following deposition onto a carbon-coated copper grid and staining with 2% (wt/vol) potassium phosphotungstate (pH 7.0), the grids were observed using a Hitachi H7500 transmission electron microscope (TEM) operating at 80 kV.

Phage physical characterization
The stability of Staphylococcus phage JD007 The acid-base stability of phage JD007 was assessed. One hundred-fold dilutions of the initial phage titre (10 8 pfu/mL) were performed using TM buffer at different pH values ranging from 2 to 11, followed by a 2-h incubation at 37°C. Several samples were adequately diluted, and 50 μL was collected to assess the phage titre using the double-layer plate method. For the thermal stability assessment of phage JD007, a titre of 10 8 pfu/mL was incubated separately at different temperatures (37, 42, 50, 60, 70 and 80°C) for 1 h, The titre of JD007 at different temperatures was confirmed using the double-layer plate method with S. aureus Sa60 cultured to an OD 600nm ≈ 0.4.

Adsorption rate of Staphylococcus phage JD007 to Sa60
Phage JD007 and S. aureus Sa60 were incubated at an MOI of 0.01 for 1, 2, 3, 4, 5, 6, 10, 15, 20, 30, or 40 min; the mixtures were centrifuged at 16,000 × g for 30 s. The titre of phage in the supernatant was identified using the double-layer agar plate method. The titre of phages previously mixed with S. aureus Sa60 was considered time "0". The proportion of the amount of non-adsorbed phages to the amount of phages used for infection, based on three independent experiments, is shown, and standard deviations are indicated.
The inhibition assays of Staphylococcus phage JD007 to S. aureus Sa60 S. aureus Sa60 was cultured to an OD 600nm ≈ 0.2 and then infected by phage JD007 at an MOI of 0, 0.1, or 1. The mixtures were incubated at 37°Cwith shaking at 100 rpm, and the OD 600nm value at different time points was measured using a BioPhotometer plus.

The host range of Staphylococcus phage JD007
The host range was analysed by spotting serial dilutions of phage JD007 on a double-layer soft agar lawn of different S. aureus isolates obtained from different hospitals. Two microlitres of concentrated phage lysate (≈10 8 pfu/mL) and serial dilutions were plated onto LB medium plates overlaid with S. aureus (OD 600nm ≈ 0.4) mixed with 0.7% top agar (cultured for 30 min before adding), followed by overnight incubation. The strains with the same infectivity as the control strain (S. aureus Sa60) have an efficiency of plating of 1 (probably the equivalent of your "++++"). The inhibition zone of the spots formed by spot testing were seen as a common system for assessing the success of infection by the phage: ++ ++ complete clearing; +++ clearing throughout but with faintly hazy background; ++ substantial turbidity throughout the cleared zone; + a few individual plaques; 0 no clearing, − but you may see a spot where the pipette tip touched the agar [23]. In total, 175 isolates of S. aureus were used to identify the host range of phage JD007.
Safety assessment of Staphylococcus phage JD007 based on the genome sequence The complete genome sequence of JD007 has been reported, GenBank accession number is JX878671 [17]. It was sequenced using Roche 454 Sequencing Technology and assembled using the method of Newbler Metrics Results Software Release v. 2.7 the platform provided. All of the annotated genes were compared to an antibiotic resistance gene database (ARDB, http://ardb.cbcb.umd.edu/) and a virulence factor database (http:// www.mgc.ac.cn/VFs/main.htm). Genes with more than 70% coverage and 30% identity were retained. The comparison of complete genome sequences of phage JD007 with phage K [24], Twort [4] and GH15 [25] using Mauve20150226 [26]. The maximum likelihood trees were constructed using MEGA5 [27].