Abstract
Panton-Valentine leukocidin (PVL, encoded by lukSF-PV genes), a bi-component and pore-forming toxin, is carried by different staphylococcal bacteriophages. The prevalence of PVL in Staphylococcus aureus has been reported around the globe. However, the data on PVL-encoding phage types, lukSF-PV gene variation and chromosomal phage insertion sites for PVL-positive S. aureus are limited, especially in China. In order to obtain a more complete understanding of the molecular epidemiology of PVL-positive S. aureus, an integrated and modified PCR-based scheme was applied to detect the PVL-encoding phage types. Phage insertion locus and the lukSF-PV variant were determined by PCR and sequencing. Meanwhile, the genetic background was characterized by staphylococcal cassette chromosome mec (SCCmec) typing, staphylococcal protein A (spa) gene polymorphisms typing, pulsed-field gel electrophoresis (PFGE) typing, accessory gene regulator (agr) locus typing and multilocus sequence typing (MLST). Seventy eight (78/1175, 6.6%) isolates possessed the lukSF-PV genes and 59.0% (46/78) of PVL-positive strains belonged to CC59 lineage. Eight known different PVL-encoding phage types were detected, and Φ7247PVL/ΦST5967PVL (n = 13) and ΦPVL (n = 12) were the most prevalent among them. While 25 (25/78, 32.1%) isolates, belonging to ST30, and ST59 clones, were unable to be typed by the modified PCR-based scheme. Single nucleotide polymorphisms (SNPs) were identified at five locations in the lukSF-PV genes, two of which were non-synonymous. Maximum-likelihood tree analysis of attachment sites sequences detected six SNP profiles for attR and eight for attL, respectively. In conclusion, the PVL-positive S. aureus mainly harbored Φ7247PVL/ΦST5967PVL and ΦPVL in the regions studied. lukSF-PV gene sequences, PVL-encoding phages, and phage insertion locus generally varied with lineages. Moreover, PVL-positive clones that have emerged worldwide likely carry distinct phages.
Introduction
Staphylococcus aureus causes a spectrum of diseases from minor skin and soft tissue infections (SSTIs) to life-threatening conditions due to its potential to produce many toxins and efficiency at overcoming antibiotics (David and Daum, ; Uhlemann et al., ). Panton-Valentine leukocidin (PVL), found in some strains of S. aureus, plays a key role in the leukocytolysis and tissue necrosis (Shallcross et al., ). The genes (lukSF-PV) encoding PVL are composed of two co-transcribed open reading frames (lukS-PV and lukF-PV), and located on lysogenized bacteriophages integrated into S. aureus chromosome (Boyle-Vavra and Daum, ; Shallcross et al., ).
To date, all known S. aureus phages belong to the order Caudovirales, which can be separated into three major families (Podoviridae, Siphoviridae, and Myoviridae) depending on the tail morphology (Xia and Wolz, ). At least 10 PVL phages have been described and sequenced, and all of them belong to the Siphoviridae family characterized by double-stranded DNA and a long non-contractile tail (Zhang et al., ; Xia and Wolz, ). As the mobile genetic elements, phages can be mobilized and transferred among S. aureus strains. Hence, investigating the typing of PVL-encoding phages among S. aureus may help to better understand the evolution of these pathogens. The PCR-based scheme targeting a small part of the phage genome is a cursory means of predicting phages types inexpensively. However, there is no unified pattern for this scheme. To identify as many PVL-encoding phages as possible, this article summarized a strategy by integrating and modifying the expanded PCR-based scheme described previously (Ma et al., ; Boakes et al., ; Chen et al., ; Sanchini et al., ). Subsequently, by the usage of this strategy, fifteen-reaction PCR assay was carried out to detect 10 of the PVL-encoding phages (ΦPVL, Φ108PVL, Φtp310-1, ΦSa2958, ΦSa2mw, ΦSLT, ΦSa2USA, ΦTCH60, Φ7247PVL/ ΦST5967PVL, and ΦSa119) in S. aureus from China.
At least 22 single-nucleotide polymorphisms (SNPs) have been identified in the lukSF-PV genes (Wolter et al., ; O'Hara et al., ; Boakes et al., ; Li et al., ; Chen et al., ; Sanchini et al., ). And several PVL protein isoforms carrying non-synonymous mutations have recently been revealed (O'Hara et al., ; Li et al., ), which may implicate functional significance. As we all know, phages lysogenize into the bacterial chromosome through the integrative pathway. Boakes et al. () analyzed the sequence diversity at the insertion site for the different PVL-encoding phages and found two lineage-specific insertion sites within the S. aureus chromosome: Ins1 and Ins2. Of note, site-specific attachment sequences (attL and attR) are also conserved across lineages of PVL-carrying strains although some variations were found (Boakes et al., ; Chen et al., ).
Due to the absence of a uniform detection method, complete data on Chinese epidemiology of PVL-encoding phages are limited (Li et al., ; Hu et al., ). In addition, no details yet about chromosomal PVL-encoding phage insertion sites for S. aureus in China are reported. In the present study, we aimed to obtain a more complete description of the molecular epidemiology of PVL-positive S. aureus from China by detecting PVL-encoding phage types, analyzing PVL variant alleles and the chromosomal phage insertion junctions, and determining the genetic background. Moreover, the relationships among them were also elucidated.
Materials and methods
Bacterial isolates
During January 2010 to May 2015, a total of 1175 consecutive, non-duplicate clinical S. aureus isolates were collected from seven hospitals in China, namely, Shanghai General Hospital, Shanghai sixth People's Hospital, Tongren Hospital, Ruijin hospital, Shanghai People's Hospital of Putuo District, Zhejiang Xiaoshan Hospital, and The Central Hospital of Lishui City, Zhejiang province (Figure 1). All the isolates, including 924 methicillin-resistant S. aureus (MRSA) strains and 251 methicillin-susceptible S. aureus (MSSA) strains, were identified by VITEK Systems (BioMérieux, Marcy l′ Etoile, France), the susceptibility of cefoxitin (30 μg, Oxoid, Basingstoke, UK), and the presence of mecA and mecC genes (Bignardi et al., ; Ganesan et al., ; Clinical Laboratory Standards Institute, ). The presence of lukSF-PV genes was determined by PCR according to previously published method (Lina et al., ).
Figure 1
Single-nucleotide polymorphisms of the lukSF-PV genes
All isolates were cultured on blood agar and incubated overnight at 37°C. Genomic DNA was extracted by TIANamp Bacterial DNA Kit (TIANGEN BIOTECH Co., Ltd., Beijing, China). Three primer pairs were designed to amplify three fragments (654, 718, and 680 nucleotides in length, respectively) of lukSF-PV genes as described by Boakes et al. (). All the products of PCR were sequenced in both directions by Shanghai Sangon Biotech.
Characterization of PVL-positive isolates
PVL-positive S. aureus were characterized by staphylococcal cassette chromosome mec (SCCmec) typing (Zhang et al., ), staphylococcal protein A (spa) gene polymorphisms typing (Koreen et al., ), pulsed-field gel electrophoresis (PFGE) typing (Mulvey et al., ), accessory gene regulator (agr) locus typing (Lina et al., ) and multilocus sequence typing (MLST; Enright et al., ). The spa type for each isolate was obtained based on a website (http://www.ridom.de/spaserver/). Sequence types (STs) were determined by allelic profile according to the MLST database (http://saureus.mlst.net/). And clustering of related STs into clonal complexes (CCs) was analyzed using eBURST (http://www.mlst.net).
PVL-encoding phage typing
Integrating the expanded PCR-based scheme described by previous studies (Ma et al., ; Boakes et al., ; Chen et al., ; Sanchini et al., ) with a little modification, three sets of PCRs including 15 PCR reactions (Set 1: PCR-1, -4 -7, and -8, Set 2: PCR-2, -5, and -9 to -13, and Set 3: PCR-3, -6, -14, and -15), were carried out to detect 10 of the PVL-encoding phages (ΦPVL, Φ108PVL, Φtp310-1, ΦSa2958, ΦSa2mw, ΦSLT, ΦSa2USA, ΦTCH60, Φ7247PVL/ ΦST5967PVL, and ΦSa119) (Table 1). PCR-1 to -3 (step-1) targeting genes encoding icosahedral or elongated head shape allowed classification of elongated-head group, icosahedral-head group I and icosahedral-head group II phages; PCR-4 to -6 (step-2) were used to link these morphologically specific tail genes to the lukSF-PV genes with primer pairs commonly conserved among each of the groups; PCR-7 to -15 (step-3) were designed to detect 10 specific PVL-encoding phages. All the PVL-positive isolates were detected by PCR-1 to -3 (step-1). Only when a positive result was obtained from step-1, were step-2, and step-3 performed. Phages which were positive by icosahedral/elongated head classification assays (PCR-1 to -3) and characterization assays (PCR-7 to -15) but negative for linkage assays (PCR-4 and -6), were defined as “phage-like” depending on the existence of known individual phage type (Chen et al., ). The workflow and the detailed PCR-based scheme for PVL-encoding phage typing were summarized in Figure 2 and Table 1, respectively.
Table 1
| Step 1 | Step 2 | Step 3 | ||||||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| Primer use | Primer name | Target locus | Size of PCR product (bp) | References | PCR's order number | Primer use | Primer name | Target locus | Size of PCR product (bp) | References | PCR's order number | PVL phage(s) | Primer use | Primer (pair) name | Target locus | Size of PCR product (bp) | GenBank accession no. | References | PCR's order number | |
| Set 1 | Icosahedral-head group I | portal-1F | por | 569 | Ma et al., | PCR-1 | To link icosahedral-head group I specific tail genes to the PVL genes. | tail-ico-F | mtp | 8570 | Ma et al., | PCR-4 | ΦPVL | To Classify individual PVL phages by detecting the gene lineage between the integrase gene and the genes located downstream of the gene. | intF-2/PVL-aR | int/JP030 | 1411 | AB009866 | Ma et al., | PCR-7 |
| portal-1R | por | Ma et al., | 10,497 | Φ108PVL | intF-2/108-aR | int/ant | 4340 | AB243556 | Ma et al., | PCR-8 | ||||||||||
| tail-1F | mtp | 489 | Ma et al., | lukSR1 | lukS-PV | 8574 | Ma et al., | Φtp310-1* | intF-2/PVL-aR /108-aR* | int/JP030/ ant | 1411/4579* | EF462197 | Ma et al., ; Chen et al., | PCR-7, -8 | ||||||
| tail-1R | mtp | Ma et al., | ||||||||||||||||||
| Set 2 | Elongated-head group | portal-2F | por | 656 | Ma et al., | PCR-2 | To link elongated-head group specific tail genes to the PVL genes. | tailE-F2 | mtp | 9483 | Otter et al., | PCR-5 | ΦSa2958 | To Classify individual PVL phages by detecting the gene lineage between the integrase gene and the genes located downstream of the gene. | intF-2/2958-aR | int/JP004 | 2238 | AP009363 | Ma et al., | PCR-9 |
| portal-2R | por | Ma et al., | 9484 | ΦSa2mw | intF-2/MW2-aR | int/ cro | 4065 | BA000033 | Ma et al., | PCR-10 | ||||||||||
| tail-2F | mtp | 468 | Ma et al., | lukSR1 | lukS-PV | 9486 | Ma et al., | ΦSLT | intF-2/SLT-aR | int/ ssb | 8770 | AB045978 | Ma et al., | PCR-11 | ||||||
| tail-2R | mtp | Ma et al., | 9484 | ΦSa2USA | Sa2USA-F /Sa2USA-R2 | phiSLT ORF484-like/lukS-PV | 679 | CP000255 | Boakes et al., ; Sanchini et al., | PCR-12 | ||||||||||
| 9482 | ΦTCH60 | intF-2/ TCH60-aR | int/HMPREF0772-11656 | 2675 | NC-017342 | Ma et al., ; Chen et al., | PCR-13 | |||||||||||||
| Set 3 | Icosahedral-head group II | portal-3F | por | 535 | Sanchini et al., | PCR-3 | To link icosahedral-head group II specific tail genes to the PVL genes. | TAIL-5 | mtp | 9164 | Sanchini et al., | PCR-6 | Φ7247PVL/ΦST5967PVL | To Classify individual PVL phages by detecting the gene lineage between the integrase gene and the genes located downstream of the gene. | intF-2/repR | int/rep | 2965 | AP011956 /AP011955 | Ma et al., ; Chen et al., ; Sanchini et al., | PCR-14 |
| portal-3R | por | |||||||||||||||||||
| tail-3F | mtp | 842 | Sanchini et al., | lukSR1 | lukS-PV | 10.728 | Ma et al., | ΦSa119 | intF-2/SA119Ant1 | int/ant | 4918 | KJ596420 | Ma et al., ; Sanchini et al., | PCR-15 | ||||||
| tail-3R | mtp | |||||||||||||||||||
PCR-based scheme for PVL-phage typing.
Φtp310-1 harbors an icosahedral head and can be identified by the presence of PCR products specific both for ΦPVL and Φ108PVL (positive both for PCR-7 and -8).
Figure 2
In case of the fragments amplified >3000 nucleotides, the PCRs were carried out using PrimeSTAR® Max DNA Polymerase (TaKaRa, Dalian, China) adhering to the protocols recommended by manufacturer.
Phage insertion locus sequencing
PCR designed to detect the proximal and distal junctions of PVL-encoding phage insertion sites were performed by the method described previously (Boakes et al.,
Results
Frequency and characterization of S. aureus harboring lukSF-PV genes
Of 1175 S. aureus isolates, 78 (6.6%) isolates were lukSF-PV genes positive, including 62 MRSA and 16 MSSA. Among the 78 lukSF-PV genes positive strains, 66 (55 MRSA and 11 MSSA) were agr1, 9 (6 MRSA and 3 MSSA) were agr3, and 2 (MSSA) were agr4. One isolate was unable to be classified in any of the established agr groups, which may be an agr-negative variant associated with the extensive use of antibiotics (Paulander et al.,
Table 2
| Isolate ID | location | ST | CC | SCCmec | spa type | agr type | PFGE | attL cluster | attR cluster | Phage morphology | PVL-phage type | Isoform of PVL |
|---|---|---|---|---|---|---|---|---|---|---|---|---|
| MRSA (n = 62) | ||||||||||||
| 7 | Zhejinag | ST1 | CC1 | V | t114 | I | K | I | I | Icosahedral II | Φ7247PVL/ ΦST5967PVL | H2 |
| LS2183 | Zhejinag | ST59 | CC59 | NT | t437 | I | B1 | I | I | Icosahedral II | NT | H1 |
| 161 | Zhejinag | ST188 | CC1 | NT | t189 | I | M | NS | NS | Icosahedral II | NT | H2 |
| SH12 | Shanghai | ST88 | CC88 | NT | t12147 | III | A1 | IV | III | Icosahedral I | ΦPVL | H2 |
| 244 | Zhejinag | ST88 | CC88 | NT | t7637 | III | C7 | IV | III | Elongated | ΦTCH60 | H2 |
| 104 | Zhejinag | ST88 | CC88 | NT | t5269 | III | N | IV | III | Elongated | ΦTCH60* | R2 |
| 148 | Zhejinag | ST9 | CC9 | IVb | t309 | I | O | VI | II | Icosahedral I | ΦPVL | H2 |
| 203 | Zhejinag | ST1 | CC1 | IVa | t127 | III | H3 | I | NS | Icosahedral I | Φ108PVL* | H1 |
| SH9 | Shanghai | ST59 | CC59 | IVa | t441 | I | I1 | I | I | Icosahedral I | ΦPVL | H1 |
| 209 | Zhejinag | ST59 | CC59 | IVa | t437 | I | E3 | II | I | Elongated | ΦSa2958 | H1 |
| LS373 | Zhejinag | ST59 | CC59 | IVa | t1751 | I | B3 | I | I | Icosahedral II | NT | H2 |
| 60 | Zhejinag | ST59 | CC59 | IVa | t437 | I | D3 | NS | I | Icosahedral II | NT | H2 |
| 34 | Zhejinag | ST59 | CC59 | IVa | t437 | I | J | I | I | Icosahedral II | NT | H2 |
| 65 | Zhejinag | ST59 | CC59 | IVa | t437 | I | D1 | I | I | Icosahedral II | NT | H2 |
| 79 | Zhejinag | ST59 | CC59 | IVa | t437 | I | F2 | I | I | Icosahedral II | ΦSa119 | H2 |
| 186 | Zhejinag | ST59 | CC59 | IVa | t437 | I | I1 | I | I | Icosahedral II | ΦSa119 | H2 |
| 238 | Zhejinag | ST59 | CC59 | IVa | t437 | I | C3 | II | I | Icosahedral II | ΦSa119* | H2 |
| 212 | Zhejinag | ST59 | CC59 | IVa | t437 | I | E1 | II | I | Elongated | ΦSa2958 | H2 |
| 226 | Zhejinag | ST59 | CC59 | IVa | t437 | I | C3 | II | I | Elongated | ΦSa2958 | H2 |
| 246 | Zhejinag | ST338 | CC59 | IVa | t437 | I | C5 | II | III | Elongated | ΦTCH60 | H2 |
| 202 | Zhejinag | ST88 | CC88 | IVa | t7637 | III | H1 | IV | III | Elongated | ΦTCH60 | H2 |
| 213 | Zhejinag | ST59 | CC59 | IVa | t1451 | I | H1 | NS | NS | Icosahedral II | NT | H2 |
| LS1939 | Zhejinag | ST59 | CC59 | Iva | t437 | I | B3 | I | I | Icosahedral II | NT | H2 |
| SH6 | Shanghai | ST59 | CC59 | IVa | t441 | I | A3 | I | I | Icosahedral I | ΦPVL | H2 |
| 40 | Zhejinag | ST59 | CC59 | IVa | t437 | I | J | II | I | Icosahedral II | ΦSa119 | R1 |
| 218 | Zhejinag | ST59 | CC59 | IVa | t437 | NT | E2 | II | NS | Icosahedral II | NT | R2 |
| 206 | Zhejinag | ST59 | CC59 | IVa | t437 | I | H2 | II | I | Icosahedral I | ΦPVL | R2 |
| LS2032 | Zhejinag | ST59 | CC59 | III | t437 | I | A2 | I | I | Icosahedral II | NT | H1 |
| 130 | Zhejinag | ST59 | CC59 | III | t441 | I | F2 | I | I | Icosahedral II | ΦSa119 | H1 |
| 242 | Zhejinag | ST338 | CC59 | III | t441 | I | C2 | I | I | Icosahedral II | Φ7247PVL/ ΦST5967PVL | H1 |
| 256 | Zhejinag | ST59 | CC59 | III | t5983 | I | D2 | I | I | Icosahedral II | Φ7247PVL/ ΦST5967PVL | H1 |
| 211 | Zhejinag | ST22 | CC22 | III | t5983 | I | G1 | VI | II | Icosahedral I | Φtp310-1 | H1 |
| LS1268 | Zhejinag | ST59 | CC59 | III | t437 | I | B3 | I | I | Icosahedral II | NT | H1 |
| 233 | Zhejinag | ST30 | CC30 | III | t1749 | III | L | VII | V | Icosahedral II | NT | H2 |
| 219 | Zhejinag | ST338 | CC59 | III | t437 | I | C1 | I | I | Icosahedral II | NT | H2 |
| 82 | Zhejinag | ST59 | CC59 | III | t437 | I | F1 | I | I | Icosahedral II | NT | H2 |
| 253 | Zhejinag | ST59 | CC59 | III | t437 | I | D1 | I | I | Icosahedral II | NT | H2 |
| 8 | Zhejinag | ST59 | CC59 | III | t437 | I | D1 | I | NS | Icosahedral II | NT | H2 |
| 204 | Zhejinag | ST59 | CC59 | III | t437 | I | E3 | I | I | Icosahedral II | NT | H2 |
| 210 | Zhejinag | ST59 | CC59 | III | t437 | I | E4 | I | I | Icosahedral II | NT | H2 |
| 214 | Zhejinag | ST22 | CC22 | III | t5983 | I | G1 | IV | II | Icosahedral I | ΦPVL | H2 |
| 221 | Zhejinag | ST338 | CC59 | III | t437 | I | C1 | I | I | Icosahedral I | ΦPVL | H2 |
| 69 | Zhejinag | ST217 | CC22 | III | t309 | I | G3 | III | II | Icosahedral I | ΦPVL* | H2 |
| 195 | Zhejinag | ST59 | CC59 | III | t437 | I | P | I | I | Icosahedral II | ΦSa119 | H2 |
| 51 | Zhejinag | ST59 | CC59 | III | t441 | I | D2 | I | I | Elongated | ΦSa2958 | H2 |
| 220 | Zhejinag | ST338 | CC59 | III | t437 | I | C1 | I | I | Icosahedral II | Φ7247PVL/ ΦST5967PVL | H2 |
| 74 | Zhejinag | ST59 | CC59 | III | t437 | I | F1 | I | I | Icosahedral II | Φ7247PVL/ ΦST5967PVL | H2 |
| 170 | Zhejinag | ST59 | CC59 | III | t437 | I | F2 | I | I | Icosahedral II | Φ7247PVL/ ΦST5967PVL | H2 |
| 254 | Zhejinag | ST59 | CC59 | III | t437 | I | D1 | I | I | Icosahedral II | Φ7247PVL/ ΦST5967PVL | H2 |
| 237 | Zhejinag | ST59 | CC59 | III | t437 | I | C4 | I | I | Icosahedral II | Φ7247PVL/ ΦST5967PVL | H2 |
| 236 | Zhejinag | ST59 | CC59 | III | t437 | I | C4 | I | I | Icosahedral II | Φ7247PVL/ ΦST5967PVL | H2 |
| 108 | Zhejinag | ST59 | CC59 | III | t437 | I | F2 | I | I | Icosahedral II | Φ7247PVL/ ΦST5967PVL | H2 |
| 145 | Zhejinag | ST59 | CC59 | III | t437 | I | D2 | I | I | Icosahedral II | Φ7247PVL/ ΦST5967PVL | H2 |
| 217 | Zhejinag | ST59 | CC59 | III | t437 | I | E1 | I | I | Icosahedral II | Φ7247PVL/ ΦST5967PVL | H2 |
| 16 | Zhejinag | ST59 | CC59 | III | t437 | I | D3 | I | I | Icosahedral II | Φ7247PVL/ ΦST5967PVL | H2 |
| 205 | Zhejinag | ST22 | CC22 | III | t5983 | I | G1 | I | II | Icosahedral I | Φtp310-1 | H2 |
| 239 | Zhejinag | ST22 | CC22 | III | t5983 | I | G2 | I | II | Icosahedral I | Φtp310-1 | H2 |
| RJ117 | Shanghai | ST338 | CC59 | III | t437 | I | A4 | I | I | Icosahedral II | NT | H2 |
| LS2137 | Zhejinag | ST59 | CC59 | III | t2755 | I | A6 | I | I | Icosahedral II | NT | H2 |
| 3 | Zhejinag | ST59 | CC59 | III | t034 | I | Q | III | III | Elongated | ΦSa2USA | H2 |
| 222 | Zhejinag | ST149 | CC5 | II | t437 | I | C6 | I | I | Elongated | ΦSa2USA | H2 |
| 255 | Zhejinag | ST149 | CC5 | II | t437 | I | K | II | III | Icosahedral I | Φtp310-1 | H2 |
| MSSA (n = 16) | ||||||||||||
| LS126 | Zhejinag | ST30 | CC30 | – | t318 | I | B2 | I | I | Icosahedral II | NT | H1 |
| SH14 | Shanghai | ST398 | CC398 | – | t034 | I | A5 | III | III | Elongated | NT | H1 |
| LS1985 | Zhejinag | ST22 | CC22 | – | t309 | I | B1 | IV | II | Icosahedral I | ΦPVL | H1 |
| LS2074 | Zhejinag | ST25 | CC25 | – | t227 | I | B1 | NS | IV | Icosahedral I | ΦPVL | H1 |
| SH26 | Shanghai | ST25 | CC25 | – | t078 | I | A2 | NS | NS | Elongated | ΦSa2958 | H1 |
| SH135 | Shanghai | ST1301 | CC121 | – | t12145 | IV | A1 | V | VI | Icosahedral I | NT | H1 |
| LS1940 | Zhejinag | ST30 | CC30 | – | t318 | III | A7 | I | VI | Icosahedral II | NT | H1 |
| LS2078 | Zhejinag | ST172 | singleton | – | t078 | I | B1 | NS | VI | Icosahedral I | ΦPVL | H1 |
| LS1004 | Zhejinag | ST160 | CC121 | – | t15797 | IV | B3 | V | VI | Icosahedral I | NT | H2 |
| LS1966 | Zhejinag | ST88 | CC88 | – | t15796 | III | B2 | IV | III | Icosahedral I | ΦPVL | H2 |
| SH13 | Shanghai | ST88 | CC88 | – | t2310 | III | A3 | IV | III | Elongated | ΦTCH60 | H2 |
| SH25 | Shanghai | ST217 | CC22 | – | t309 | I | A4 | VIII | II | Elongated | ΦSa2958 | H2 |
| SH19 | Shanghai | ST217 | CC22 | – | t309 | I | A4 | VIII | II | Elongated | ΦSa2958 | H3 |
| SJ1775 | Shanghai | ST398 | CC398 | – | t034 | I | I2 | III | III | Elongated | ΦSa2USA | R1 |
| SH3 | Shanghai | ST398 | CC398 | – | t1255 | I | A3 | III | NS | Elongated | ΦSa2USA | R1 |
| LS1911 | Zhejinag | ST398 | CC398 | – | t034 | I | B3 | III | III | Elongated | ΦSa2USA | R1 |
Characteristics and PVL-phage typing of 78 PVL-positive S. aureus isolates.
Phages-like, positive by icosahedral/elongated head classification assays (PCR-1 to -3) and characterization assays (PCR-7 to -15) but negative for linkage assays (PCR-4 and -6). –, no SCCmec elements was detected (MSSA).
NT, non-typeable.
NS, could not be sequenced.
Typing of PVL-encoding phages
Applying the expanded PCR-based scheme described above, it was possible to identify 10 of known phages carrying lukSF-PV genes. Our results showed that 67.9% (53/78, 42 MRSA and 11 MSSA) of PVL-positive isolates could be divided into 8 phage types [ΦPVL (n = 12, including one ΦPVL-like), Φ108PVL-like (n = 1), Φtp310-1 (n = 4), ΦSa2958 (n = 7), ΦSa2USA (n = 5), ΦTCH60 (n = 5, including one ΦTCH60-like), Φ7247PVL/ΦST5967PVL (n = 13), and ΦSa119 (n = 6, including one ΦSa119-like)], which were almost equally split among elongated-head group (n = 17), icosahedral-head group I (n = 17), and II (n = 19). Phage typeability was 68.8% in MSSA (11 out of 16) and 67.7% in MRSA (42 out of 62). Unfortunately, 25 isolates (32.1%) could not be typed by the present scheme. Four MRSA isolates were considered to be “phage-like,” named ΦPVL-like, Φ108PVL-like, ΦTCH60-like and ΦSa119-like, according to the definition from materials and methods. The details of PVL-encoding phages types were shown in Table 2.
Variation of lukSF-PV genes
The amplification products of the lukSF-PV genes from 78 isolates were sequenced, and nucleotide variations were seen at five sites (position 527 and 663 located in the lukS locus and position 1022, 1396, 1729 located in the lukF locus) using the lukSF-PV genes of ΦSLT as a reference. Of 78 PVL-positive isolates, 71 (91.0%, 71/78) were of H variant as defined by O'Hara et al. (
DNA sequencing of phage/chromosome junctions
To investigate the chromosomal insertion site of phage in the present study, we performed PCR based on the known phage/chromosome junctions in the other strains (Boakes et al.,
Figure 3

ML analysis trees for (A) attR attachment site (71 isolates) and (B) attL attachment site (72 isolates). Among 78 PVL-positive isolates, the attR and attL attachment site sequences cannot be identified in 7 and 6 strains, respectively. Branching numbers represent bootstrap values.
Discussion
PVL-positive S. aureus, strongly associated with SSTI and necrotizing pneumonia, has gained great attention in recent years (David and Daum,
A previously published eight-reaction PCR assay was performed to identify five (ΦSa2958, ΦSa2mw, ΦPVL, Φ108PVL, and ΦSLT) of the PVL-encoding phages in S. aureus (Ma et al.,
Applying the PCR-based strategy, the distribution of PVL-encoding phages was investigated all over the world. Two morphologically distinct phages (ΦPVL and ΦSa2958) were found to be predominant among Japanese PVL-positive MRSA (Ma et al.,
These findings are similar with those reported in China previously (Li et al.,
The correlation between the presence of PVL in S. aureus infections and the clinical outcome has been controversial because of some conflicting data from epidemiological study or in vivo animal models (Hermos et al.,
Another interesting finding was that all the attR cluster I-harboring isolates possessed the same agr type (agr I) and were mostly observed within CC59 lineage. With the exception of one strain (strain 148), all of the attR sequences of cluster II strains were identified in CC22 lineage. The CC88 isolates were completely associated with attR cluster III and attL cluster IV. The main variations in attL were located in the phage binding sites (P1 and P2) and those in attR were located in the phage binding site P3 and chromosome binding sites B2. These indicated that the diversity of junction sequences was mainly due to differences in the phage DNA, which was different from the previous report that the main variations in attR were located in the 29-bp chromosome binding sites (B2) (Chen et al.,
There were two limitations in this study. First, the identification of PVL-encoding phages was on the basis of PCR assays targeting only a part of the phage genome. Since phages tend to have recombination events, positivity by PCR should be supposed to be a family of PVL-encoding phages with variable genomic portions, instead of being taken as a specific phage (Sanchini et al.,
Conclusion
This study characterized PVL-encoding phages, the chromosomal phage insertion sites, the polymorphism of lukSF-PV genes, and the genetic background of PVL-positive S. aureus clinical isolates from China, and found the existence of some correlation among them. Our findings may contribute to the understanding of the epidemiology and evolution of PVL-positive S. aureus, and add the evidence that PVL-positive strains disseminating worldwide likely carry distinct PVL phages.
Statements
Author contributions
HZ, FH performed the experiments; QL designed and conceived the study; HZ, QL, CH analyzed the clinical data and wrote this manuscript; QL, SJ, XX, YZ, BD, and FG collected the clinical samples. All authors read and approved the final manuscript.
Acknowledgments
The authors thank Jinwei Huang, Jianmin Ren, Chuanling Zhang, Lizhong Han, Zhengyin Zhang, and Jin Tang to provide the clinical S. aureus isolates from Zhejiang Province and Shanghai. This work was supported by grants from Natural Science Foundation, Science and Technology Commission of Shanghai (No. 12ZR1425000), the National Natural Science Foundation of China (No. 81371872), and Putuo District, Major Scientific and Technological Research Projects (No. KW15101).
Conflict of interest
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.
References
1
BignardiG. E.WoodfordN.ChapmanA.JohnsonA. P.SpellerD. C. (1996). Detection of the mec-A gene and phenotypic detection of resistance in Staphylococcus aureus isolates with borderline or low-level methicillin resistance. J. Antimicrob. Chemother.37, 53–63. 10.1093/jac/37.1.53
2
BoakesE.KearnsA. M.GannerM.PerryC.HillR. L.EllingtonM. J. (2011). Distinct bacteriophages encoding Panton-Valentine leukocidin (PVL) among international methicillin-resistant Staphylococcus aureus clones harboring PVL. J. Clin. Microbiol.49, 684–692. 10.1128/JCM.01917-10
3
Boyle-VavraS.DaumR. S. (2007). Community-acquired methicillin-resistant Staphylococcus aureus: the role of Panton-Valentine leukocidin. Lab. Invest.87, 3–9. 10.1038/labinvest.3700501
4
ChenL.ChavdaK. D.SolankiM.MediavillaJ. R.MathemaB.SchlievertP. M.et al. (2013). Genetic variation among Panton-Valentine leukocidin-encoding bacteriophages in Staphylococcus aureus clonal complex 30 strains. J. Clin. Microbiol.51, 914–919. 10.1128/JCM.03015-12
5
Clinical Laboratory Standards Institute (2014). Performance Standards for Antimicrobial Susceptibility Testing; Twenty-Second Informational Supplement (M100-S22), Vol. 34. Wayne, PA: Clinical and Laboratory Standards Institute.
6
DavidM. Z.DaumR. S. (2010). Community-associated methicillin-resistant Staphylococcus aureus: epidemiology and clinical consequences of an emerging epidemic. Clin. Microbiol. Rev.23, 616–687. 10.1128/CMR.00081-09
7
EnrightM. C.DayN. P.DaviesC. E.PeacockS. J.SprattB. G. (2000). Multilocus sequence typing for characterization of methicillin-resistant and methicillin-susceptible clones of Staphylococcus aureus. J. Clin. Microbiol.38, 1008–1015. 10.1128/JCM.43.9.4448-4454.2005
8
FanJ.ShuM.ZhangG.ZhouW.JiangY.ZhuY.et al. (2009). Biogeography and virulence of Staphylococcus aureus. PLoS ONE4:e6216. 10.1371/journal.pone.0006216
9
GanesanA.CrawfordK.MendeK.MurrayC. K.LloydB.EllisM.et al. (2013). Evaluation for a novel methicillin resistance (mecC) homologue in methicillin-resistant Staphylococcus aureus isolates obtained from injured military personnel. J. Clin. Microbiol.51, 3073–3075. 10.1128/JCM.01516-13
10
HermosC. R.YoongP.PierG. B. (2010). High levels of antibody to Panton-Valentine leukocidin are not associated with resistance to Staphylococcus aureus-associated skin and soft-tissue infection. Clin. Infect. Dis.51, 1138–1146. 10.1086/656742
11
HuQ.ChengH.YuanW.ZengF.ShangW.TangD.et al. (2015). Panton-Valentine leukocidin (PVL)-positive health care-associated methicillin-resistant Staphylococcus aureus isolates are associated with skin and soft tissue infections and colonized mainly by infective PVL-encoding bacteriophages. J. Clin. Microbiol.53, 67–72. 10.1128/JCM.01722-14
12
JohnJ. F.Jr.LindsayJ. A. (2008). Clones and drones: do variants of Panton-Valentine leukocidin extend the reach of community-associated methicillin-resistant Staphylococcus aureus?J. Infect. Dis.197, 175–178. 10.1086/524693
13
KoreenL.RamaswamyS. V.GravissE. A.NaidichS.MusserJ. M.KreiswirthB. N. (2004). spa typing method for discriminating among Staphylococcus aureus isolates: implications for use of a single marker to detect genetic micro- and macro variation. J. Clin. Microbiol.42, 792–799. 10.1128/JCM.42.2.792-799.2004
14
LiX.SunJ.WuD.WangL.YangY.WangC.et al. (2012). Panton-Valentine leukocidin gene sequence variation and phage in methicillin-resistant and methicillin-susceptible Staphylococcus aureus from children in mainland China. Microbiol. Immunol.56, 155–162. 10.1111/j.1348-0421.2011.00422.x
15
LiZ.StevensD. L.HamiltonS. M.ParimonT.MaY.KearnsA. M.et al. (2011). Fatal S. aureus hemorrhagic pneumonia: genetic analysis of a unique clinical isolate producing both PVL and TSST-1. PLoS ONE6:e27246. 10.1371/journal.pone.0027246
16
LinaG.BoutiteF.TristanA.BesM.EtienneJ.VandeneschF. (2003). Bacterial competition for human nasal cavity colonization: role of Staphylococcal agr alleles. Appl. Environ. Microbiol.69, 18–23. 10.1128/AEM.69.1.18-23.2003
17
LinaG.PiémontY.Godail-GamotF.BesM.PeterM. O.GauduchonV.et al. (1999). Involvement of Panton-Valentine leukocidin-producing Staphylococcus aureus in primary skin infections and pneumonia. Clin. Infect. Dis.29, 1128–1132. 10.1086/313461
18
MaX. X.ItoT.KondoY.ChoM.YoshizawaY.KanekoJ.et al. (2008). Two different Panton-Valentine leukocidin phage lineages predominate in Japan. J. Clin. Microbiol.46, 3246–3258. 10.1128/JCM.00136-08
19
MariemB. J.ItoT.ZhangM.JinJ.LiS.IIhemB. B.et al. (2013). Molecular characterization of methicillin-resistant Panton-valentine leukocidin positive Staphylococcus aureus clones disseminating in Tunisian hospitals and in the community. BMC Microbiol.13:2. 10.1186/1471-2180-13-2
20
MulveyM. R.ChuiL.IsmailJ.LouieL.MurphyC.ChangN.et al. (2001). Development of a Canadian standardized protocol for subtyping methicillin-resistant Staphylococcus aureus using pulsed-field gel electrophoresis. J. Clin. Microbiol.39, 3481–3485. 10.1128/JCM.39.10.3481-3485.2001
21
O'HaraF. P.GuexN.WordJ. M.MillerL. A.BeckerJ. A.WalshS. L.et al. (2008). A geographic variant of the Staphylococcus aureus Panton-Valentine leukocidin toxin and the origin of community-associated methicillin-resistant S. aureus USA300. J. Infect. Dis.197, 187–194. 10.1086/524684
22
OtterJ. A.KearnsA. M.FrenchG. L.EllingtonM. J. (2010). Panton-Valentine leukocidin-encoding bacteriophage and gene sequence variation in community-associated methicillin-resistant Staphylococcus aureus. Clin. Microbiol. Infect.16, 68–73. 10.1111/j.1469-0691.2009.02925.x
23
OttoM. (2011). A MRSA-terious enemy among us: end of the PVL controversy?Nat. Med.17, 169–170. 10.1038/nm0211-169
24
PaulanderW.NissenV. A.BækK. T.HaaberJ.FreesD.IngmerH. (2013). Antibiotic-mediated selection of quorum-sensing-negative Staphylococcus aureus. mBio3:e00459-12. 10.1128/mBio.00459-1213
25
SanchiniA.DelG. M.VillaL.AmmendoliaM. G.SupertiF.MonacoM.et al. (2014). Typing of Panton-Valentine leukocidin-encoding phages carried by methicillin-susceptible and methicillin-resistant Staphylococcus aureus from Italy. Clin. Microbiol. Infect.20, O840–O846. 10.1111/1469-0691.12679
26
ShallcrossL. J.FragaszyE.JohnsonA. M.HaywardA. C. (2013). The role of the Panton-Valentine leucocidin toxin in staphylococcal disease: a systematic review and meta-analysis. Lancet Infect. Dis.13, 43–54. 10.1016/S1473-3099(12)70238-4
27
UhlemannA. C.OttoM.LowyF. D.DeLeoF. R. (2014). Evolution of community- and healthcare-associated methicillin-resistant Staphylococcus aureus. Infect. Genet. Evol.21, 563–574. 10.1016/j.meegid.2013.04.030
28
WolterD. J.TenoverF. C.GoeringR. V. (2007). Allelic variation in genes encoding Panton-Valentine leukocidin from community-associated Staphylococcus aureus. Clin. Microbiol. Infect.13, 827–830. 10.1111/j.1469-0691.2007.01763.x
29
XiaG.WolzC. (2014). Phages of Staphylococcus aureus and their impact on host evolution. Infect. Genet. Evol.21, 593–601. 10.1016/j.meegid.2013.04.022
30
XiaoM.WangH.ZhaoY.MaoL. L.BrownM.YuY. S.et al. (2013). National surveillance of methicillin-resistant Staphylococcus aureus in China highlights a still-evolving epidemiology with 15 novel emerging multilocus sequence types. J. Clin. Microbiol.51, 3638–3644. 10.1128/JCM.01375-13
31
YuF.ChenZ.LiuC.ZhangX.LinX.ChiS.et al. (2008). Prevalence of Staphylococcus aureus carrying Panton-Valentine leukocidin genes among isolates from hospitalised patients in China. Clin. Microbiol. Infect.14, 381–384. 10.1111/j.1469-0691.2007.01927.x
32
ZhangK.McClureJ. A.ElsayedS.LouieT.ConlyJ. M. (2005). Novel multiplex PCR assay for characterization and concomitant subtyping of staphylococcal cassette chromosome mec types I to V in methicillin-resistant Staphylococcus aureus. J. Clin. Microbiol.43, 5026–5033. 10.1128/JCM.43.10.5026-5033.2005
33
ZhangM.ItoT.LiS.JinJ.TakeuchiF.LauderdaleT. L.et al. (2011). Identification of the third type of PVL phage in ST59 methicillin-resistant Staphylococcus aureus (MRSA) strains. FEMS Microbiol. Lett.323, 20–28. 10.1111/j.1574-6968.2011.02355.x
Summary
Keywords
Staphylococcus aureus, Panton-Valentine leukocidin, sequence variation, phage typing, genetic background
Citation
Zhao H, Hu F, Jin S, Xu X, Zou Y, Ding B, He C, Gong F and Liu Q (2016) Typing of Panton-Valentine Leukocidin-Encoding Phages and lukSF-PV Gene Sequence Variation in Staphylococcus aureus from China. Front. Microbiol. 7:1200. doi: 10.3389/fmicb.2016.01200
Received
27 May 2016
Accepted
19 July 2016
Published
03 August 2016
Volume
7 - 2016
Edited by
Jorge Blanco, University of Santiago de Compostela, Spain
Reviewed by
Ruud H. Deurenberg, University Medical Center Groningen, Netherlands; Frieder Schaumburg, University Hospital Münster, Germany
Updates

Check for updates
Copyright
© 2016 Zhao, Hu, Jin, Xu, Zou, Ding, He, Gong and Liu.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Qingzhong Liu jiaodamedicine@foxmail.com
This article was submitted to Infectious Diseases, a section of the journal Frontiers in Microbiology
†These authors have contributed equally to this work.
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.