Genotyping and characterization of prophage patterns in clinical isolates of Staphylococcus aureus

Objective Staphylococcus aureus is considered an important pathogen with a variety of virulence factors in communities and hospitals all around the world. Prophage typing is a practical technique for categorizing this bacterium. In this study, we focused on the detection of prophage patterns in methicillin-resistant S. aureus (MRSA) and methicillin-sensitive S. aureus (MSSA) strains based on their virulence factors, antimicrobial resistance patterns, and molecular typing by rep-PCR. Results Out of 126 S. aureus isolates, 45 (35.7%) were identified as MRSA. In total, 17 different prophage types were detected and 112 strains out of 126 strains contained at least one prophage. There was a statistically significant relationship between hld, hlg, eta and SGA, SGA, and SGFb, respectively. The results of the rep-PCR analysis revealed 14 different patterns among the MRSA and MSSA isolates. In conclusion, the presence of different prophage-encoded virulence factors and antibiotic-resistant genes among MRSA strains enables them to produce a broad range of diseases. Thus, diverse MRSA strains which have these prophages can be considered as a potential threat to the patient’s health in either the hospital or the community.


Introduction
Staphylococcus aureus is a gram-positive coccus and is considered a normal flora or an opportunistic pathogen which can cause community-acquired (CA) and hospitalacquired (HA) infections. S. aureus can cause a different range of diseases from soft tissue and skin infections to severe life-threatening infections such as toxic shock syndrome (types 1 and 2) [1,2]. The capability of producing different virulence factors is the main reason for the high pathogenicity of S. aureus. In addition to its high pathogenicity, this bacterium also has high adaptive powers against environmental changes [3,4]. Bacteriophages can convert non-pathogenic strains to pathogenic strains through the horizontal transduction of the virulence gene [5].
The production of virulence factors is the result of the phenotypic changes caused by lysogenic conversion which plays an important role in the biology of these species [6][7][8][9][10][11]. S. aureus phages that are involved in human diseases are classified into six categories: 3 SGA, 11 SGB, 77 SGF (with the two subgroups SGFa and SGFb) and 187 SGL which belong to the siphoviridae family, while SGD or the Twort-like phage belongs to the myoviridae family [2].
In recent years, in epidemiological studies and in assessing the genetic linkages of bacteria, common methods of bacterial typing have been replaced by molecular methods [12]. The rep-PCR technique is an alternative technique for producing fingerprint directly and without using endonuclease enzymes. Staphylococcal bacteriophages are widely used in the typing of Staphylococcus strains associated with human diseases through the rep-PCR technique [13][14][15]. Given that there are few studies about investigating the relationship between different prophage patterns, virulence factors, Dini et al. BMC Res Notes (2019) 12:669 antibiotic resistance, SCCmec, and rep types in Iran and some other countries, the current study focuses on this field.

Sampling and isolation of S. aureus
A total of 126 samples isolated in a previous study by Tahmasebi et al. were used [16]. This study was approved by the ethics committee of Hamadan University of Medical Sciences (No: IRUMSHA. REC. 1397.193).

DNA extraction
The total genomic DNAs of the isolates were extracted by the boiling method. The quality of the extracted DNAs was assessed using a Nanodrop device (A & E Lab Nano-200, UK).

Detecting the prophages of S. aureus
Multiplex PCR assay was used for identifying SGA, SGB, SGD, SGL, and SGF (and its two subgroups SGFa and SGFb) genes by their particular primers. The primers used for S. aureus bacteriophages are given in Additional file 1: Table S1. The PCR protocol used was based on the study of Pantuck et al. [1].

Detection of SCCmec types
SCCmec types I to V were detected using a multiplex-PCR assay reported previously [9].

Sequencing
One sample of each prophage PCR product (amplicon) was sequenced by Bioneer Co. (South Korea) and the data were analyzed using the Chromas software.

Rep-PCR
Rep-PCR was performed on 55 isolates chosen from the MRSA and MSSA strains based on their prophage profiles using the primers indicated in (Additional file 1: Table S1).

Analysis of the rep-PCR results
The rep-PCR band patterns were compared and clustered by Dice and unweighted paired group (UPGMA) methods using an online service (the inslico.ehu.es database), respectively.

Statistical analysis
The data were analyzed using SPSS software version 16 (Chicago, SPSS Inc., IBM, USA). Descriptive statistical methods were used to determine the frequency, percentage, and mean, while Fisher's test was used to compare the qualitative results. P ≤ 0.05 was considered as statistically significant.

Prevalence of MRSA strains
Among 126 tested S. aureus isolates, 45 strains (35.7%) showed resistance to cefoxitin and were identified as MRSA isolates. All the MRSA isolates carried the mecA gene.

Prevalence of virulence factors
The prevalence of the virulence factors is given in Additional file 2: Figure S1. The hla gene was the prevalent virulence factor (12.7%) and pvl1, pvl2, and sea showed the least (7.9%) distribution among the S. aureus isolates.

Prevalence of prophages
The results of this study showed that out of 126 strains, 112 contained at least one prophage (Additional file 3: Figure S2). SGB (88%) was the most prevalent prophage, while SGL and SGA were the least prevalent prophages (3.2%). All of the phage types were detected in the MRSA and MSSA strains, except for SGL which was detected only in the MRSA strains. In total, 17 different prophage types were detected among the S. aureus isolates in this study ( Table 2).

The relationship between prophages and virulence factors
By using Fisher's test, it was noticed that the P. values of hld/SGA, hlg/SGA, and eta/SGFb were 0.05, 0.04, and 0.04, respectively. This shows that there is a statistically significant relationship between these virulence factors and their related prophages.

The relationship between prophages and SCCmec types
The comparison of the antibiotic resistance genes and the presence of prophages showed that there was a statistically significant relationship between: SGB and mecA; SGFa and mecA, mecR, and mecC; SGFb and mecA and SCCmec IV; SGL and mecR and SCCmec IVa; and SGF and mecA ( Table 3).

The results of rep-PCR
The analysis of genetic linkage among the isolates using rep-PCR showed a similarity of 50-100% among the S. aureus isolates (Additional file 4: Figure S3). Genetic diversity was established among the S. aureus isolates by detecting 14 different rep-PCR fingerprints with the similarity cutoff of ≥ 95%. Fourteen different rep-PCR profiles, including seven common types and seven unique types, were identified.

Discussion
In the present study, the frequency of resistance to methicillin was 35.7%. Most studies in Iran have reported the prevalence of MRSA from 43 to 56% [19,20]. However, the rate of MRSA in our study is lower than those of some studies in Iran and other countries.
Our findings revealed that all the prophage types (SGA, SGB, SGD, SGF, SGFa, SGFb, and SGL) were detected among the MRSA and MSSA strains and the SGB prophage type was the predominant one. However, SGF prophage has been the predominant prophage type in previous reports from 2012 to 2018 in Iran [2,[19][20][21][22].   [20,23]. However, in our study, seventeen different prophage patterns and seven different prophage types were detected. In contrast to a previous study in Iran in which the prophage pattern SGB/SGF/GFa/SGFb (81%) was reported as the dominant pattern, in the current study, SGF/SGFa was identified as the dominant pattern in 38% of the isolates [22]. Another difference between the findings of the present research and those of other researches in Iran is that SGL was found among the studied strains of this study while it was not reported in previous studies.
According to the results of the current study, there was not a significant relationship between resistance to methicillin and prophage types. All prophage types were detected among both MRSA and MSSA strains except for SGL which was detected only among the MRSA isolates. In the present study, the hla gene detected in 12% of the strains was the predominant hemolysin gene. The hlb gene detected in 11% of the S. aureus strains was the second dominant hemolysin gene, though it was the predominant hemolysin gene and was detected in 100% of the MRSA strains in a study by Rahimi et al. in Tehran hospitals in 2018 [23]. According to the results of the current study, both MRSA and MSSA strains had the hla gene, however, the prevalence of this gene was significantly more in the MRSA strains which proves the relationship between this virulence factor and methicillin resistance in S. aureus. On the other hand, there was a statistically significant relationship between SGA and hld which is another hemolysin gene.
Another statistically significant relationship between prophages and virulence factors was the one between eta and SGFb. In contrast to Rahimi et al. [23] who showed that there was a relationship between eta and SGB, the authors of the current study did not find the same result. In the current study, the frequency of sea, hemolysin, and tst genes was lower than those of other studies in Iran and other countries [19,[24][25][26]. One of the reasons for the low abundance of these genes may be the low number of samples tested in this study compared to those of similar studies. However, the frequency of pvl genes was a little higher than those of other studies in Iran. PVL is one of the most important virulence factors and a marker for SCCmec type II which is associated with SGA phages [2,27]. In the current study, pvl genes were found among strains with SCCmec (types I, II, III, IV), and SGF, SGD, and SGB phages. This shows the community origins of these strains because the presence of the pvl gene is only limited to CA-MRSA strains. According to the PCR results, 28 out of 126 (22.22%) S. aureus isolates were pvl carriers. In 2013, Abimanyu et al. in India reported that 10 out of 25 (40%) MRSA isolates were positive for pvl genes [28] which was higher than the results of our study. In another study in Iran, Ohadian Moghadam et al.
In the current study, SCCmec I, II, III were detected among the HA-MRSA strains and SCCmec IVa, IVb, IVc, IVd were found among the CA-MRSA strains. Consistent with previous studies in Iran and in most Asian countries, SCCmec type III was found as the dominant SCCmec type among the HA-MRSA isolates [19,30,31]. In the current study, SCCmec IVa showed a higher frequency among the SCCmec types. This type was reported by Edselve et al. in Denmark. In the current study, SCCmec type V was detected in none of the isolates.
In the current study, S. aureus isolates were categorized by SCCmec typing, phage typing, and rep-PCR. In the study of Sommerhalter et al. it has been shown that the performance of rep-PCR typing is comparable to those of other typing methods. The analysis of the linkage among the isolates by rep-PCR showed a genetic diversity among the S. aureus isolates (14 different rep-PCR profiles including seven common types and seven unique types). The common types including A and G harbored 10 and 14 isolates, respectively. Isolates in the same rep profiles showed different phage patterns, virulence genes, SCCmec types, and pvl genes. Therefore, it can be concluded that there is no significant relationship between the rep patterns and the phage profiles, virulence genes, and SCCmec types. The results of studies from Iran and other countries indicated diversity among MRSA isolates from hospitals and the community. In a study by Manafi et al. in Iran [32], six main clusters of S. aureus were detected by rep-PCR and no significant relationship was observed between PVL-positive S. aureus and rep-PCR patterns. This is in line with our findings.

Conclusion
Our findings illustrated the high diversity of different prophage types, rep types, and SCCmec types in both MRSA and MSSA strains. According to the results of the rep-PCR technique, we also face diverse S. aureus isolates in hospitals in Hamedan and their circulation from the community to the hospitals and vice versa.

Limitation
The results of this study had some similarities and discrepancies with those of other studies. It is necessary to conduct studies with more samples in different areas and to use robust typing methods to obtain better results for better interpretations.