Genetic Variation among Panton-Valentine Leukocidin-Encoding Bacteriophages in Staphylococcus aureus Clonal Complex 30 Strains

Abstract

Clonal complex 30 (CC30), one of the major Staphylococcus aureus lineages, has caused extensive hospital-acquired and community-acquired infections worldwide. Recent comparative genomics studies have demonstrated that three CC30 clones-phage type 80/81, Southwest Pacific (SWP), and contemporary EMRSA-16 associated (Con) strains-shared a recent common ancestor more than 100 years ago. Panton-Valentine leukocidin (PVL), a bacteriophage encoded toxin that has been epidemiologically linked with community-associated methicillin-resistant S. aureus (CA-MRSA), has frequently been identified in CC30 clones, although the pvl gene variation and distribution of PVL-encoding phages are poorly understood. We determined here the distribution of PVL phages, PVL gene sequences, and chromosomal phage insertion sites in 52 S. aureus CC30 PVL-harboring isolates, collected from four continents over a 75-year period. Our results indicate that PVL phages with icosahedral heads, including Φ108PVL and ΦPVL, were mainly associated with phage 80/81 strains, whereas phages with elongated heads were predominantly found in SWP (ΦSa2958 and ΦTCH60) and Con (ΦSa2USA) strains. Nine single-nucleotide polymorphisms were identified in the lukSF-PV gene, with six isolates harboring the R variant that has been previously associated with CA-MRSA strains. Interestingly, all six R variant strains belonged to the same Con CC30 clone and carried a ΦSa2USA-like phage. Similar chromosomal phage insertion sites were also identified in all 52 PVL-harboring CC30 strains. These analyses provide important insights into the microepidemiology of PVL-harboring CC30 strains, while the discovery of ΦSa2USA-associated R variant strains sheds further light on the evolution of PVL-positive CA-MRSA.