Introduction
Burns patients are at an increased risk of colonization and subsequent infections by nosocomial pathogens due to the disruption of skin protective barrier and reduction of immune responses, which can lead to poor clinical outcomes and increased morbidity and mortality rates [
1]. More than 70% of deaths in burns patients results from infection with nosocomial microorganisms. Based on the literature, many pathogens are responsible for infection in burns patients, but it is well established that
Staphylococcus aureus, in particular, methicillin-resistant
S. aureus (MRSA) strains, are one of the most common nosocomial pathogens of burn wounds [
2,
3]. They could originate from the patient or transmitted easily by direct (i.e. contaminated hands or droplet) and indirect contact (i.e airborne infection). Many of the published data in Iran indicates that multidrug-resistant (MDR) bacteremia and wound infections are the most important causes of mortality in burns patients [
2,
4,
5]. Our previous study showed that there is a high rate of infection with MRSA in burns patients, which is a serious threat for the individual, and also a hazard to public health. Infection with MDR MRSA in burns patients may increase the economic burden of the healthcare system, and causes limitations to the therapeutic options available for treatment [
1–
4].
A major concern with managing MRSA infections in burns patients is the lack of awareness of the molecular and resistance patterns of MRSA. Various molecular typing techniques have been employed for
S. aureus isolates include pulsed-field gel electrophoresis, staphylococcal cassette chromosome
mecA (SCC
mec) typing, the accessory gene regulator (
agr) typing, the mec-associated hypervariable region (
dru), multilocus sequence typing, coagulase (
coa) typing, and staphylococcal protein A (
spa) typing [
6]. Although the pulsed-field gel electrophoresis technique is a standard method, it is acknowledged that polymerase chain reaction based methods such as
coa and
spa typing due to its cost-benefit, rapidity and high throughput ability could be an effective method for routine typing of MRSA isolates [
7,
8].
The coagulase protein is genetically and antigenically divergent and is an important virulence factor for
S. aureus. Sequence analysis of staphylocoagulase (SC) showed 6 regions including a signal sequence, N-terminal D1 region and D2 region, central region, 27-amino acid-repeat region and C-terminal sequence. Ten main types of SC have been described as a result of sequence diversity in the SC region [
8–
10].
The staphylococcal protein A (
spa) gene contains 3 distinct regions including Fc, X and C regions; typing is based on the number of tandem repeats and the sequence variation in region X of the protein A gene [
7]. The distribution of
spa types of MRSA strains isolated from different geographical areas of the world, show different patterns [
11].
In addition,
coa typing is also simple to perform, easily interpreted, and requires minimal laboratory skills and equipment, and therefore could be a useful addition to
spa typing for genotyping of
S. aureus, including MRSA strains [
8,
9]. In Iran, there is very little information about the genetic diversity of
S. aureus isolated from burns patients, therefore, the present study aimed to determine the resistance pattern, carriage of resistance determinants, and molecular characteristics of MRSA isolates collected from burns patients based on
coa and
spa gene polymorphism analysis.
Results
1. Sampling and antibiotic susceptibility
There were 89 S. aureus strains obtained from burns patients. All the strains under study were confirmed as MRSA due to the presence of the mecA gene as well as phenotypically. The results of antimicrobial susceptibility showed that all the isolates were susceptible to linezolid, teicoplanin, and vancomycin, but resistant to ampicillin. The highest rate of resistance was observed to erythromycin (86.5%), followed by gentamicin (84.2%), tetracycline (65.2%), amikacin (55.1%), clindamycin (51.7%), rifampicin (33.7%), and quinupristin-dalfopristin (25.8%).
The results of the microdilution method showed that 28.1% were resistant to mupirocin. Of these, 12 (13.5%) and 13 (14.6%) were HLMUPR and LLMUPR, respectively. The results of mupirocin MIC showed that 10 had MIC values of 16 μg/mL (40%), 3 had MIC of 32 μg/mL (12%), 11 (44%) had MIC ≥ 512 μg/mL, and 1 (4%) had MIC ≥ 1,024 μg/mL.
The MICs of vancomycin showed 31 (34.8%) isolates had MIC of 0.5 μg/mL, 18 (20.3%) had MIC of 1 μg/mL, 23 (25.8%) had MIC of 2 μg/mL and 17 (19.1%) had MIC of 8 μg/mL. Five different resistant phenotypes were identified among our MRSA isolates. The predominant resistance pattern included resistance to 4 antibiotics (52.8%), followed by 6 antibiotics (15.7%), 8 antibiotics (13.5%), 5 antibiotics (11.3%), and 7 antibiotics (6.7%), simultaneously.
2. Antimicrobial resistance and toxin genes
All MRSA strains in the current study harbored at least 2 antibiotic resistance genes. The most commonly detected antibiotic resistance genes were ant (4′)-Ia (60.7%), aac (6′)-Ie/aph (2″) (55.1%), aph (3′)-IIIa (46.7%), tet(M) (43.8%), erm(A) (29.2%), erm(C) (23.6%), erm(B) (20.2%), msr(B) (18%), msr(A) (13.5%), and mupA (13.5%). vanA, vanB and mupB genes were not detected among tested isolates. Analyzing toxin genes revealed that the most prevalent gene was tst (14; 15.7%), followed by pvl (9; 10.1%), eta (2; 2.2%) and etb (2; 2.2%).
3. spa typing
A total of 12 different spa genotypes were identified. t790 accounting for 18% (16/89 each) was found to be the predominant spa type, followed by t064 (13.5%; 12), t030 (11.3%; 10), t044 (10.1%; 9), t852 (9%; 8), t223 (7.9%; 7), t421 (6.7%; 6), t008 (6.7%; 6), t019 (6.7%; 6), t021 (4.5%; 4), t005 (3.4%; 3), and t10795 (2.2%; 2). PVL (lukS-lukF)-encoding genes were identified in 5 strains (55.6%, 5/9) with spa types t044, 2 strains with t852 (22.2%, 2/9), and 2 strains with t008 (22.2%, 2/9), while tst was identified among spa types t790 (64.3%, 9/14) and t223 (35.7%, 5/14). Of the 13 LLMUPR strains, 11 strains belonged to t790 (84.6%), and 2 strains (15.4%) to t008. spa types t064 (10 isolates; 83.3%) and t008 (2 isolates; 16.7%) were observed among 12 HLMUPR-MRSA strains.
4. SC typing
In the present study, 10 different types were detected and the most common coa type was type III (42; 47.2%) with spa types t790 (16; 18%), t064 (10; 11.2%), t030 (10; 11.2%) and t008 (6; 6.8%). The coa type II was observed in 13 isolates (14.6%) with spa types t019 (5; 5.6%), t044 (4; 4.5%) and t021 (4; 4.5%). The coa type VIII with spa types t852 (5; 5.6%) and t044 (5; 5.6%), were detected in 10 isolates (11.2%). From the 9 isolates with coa type I (10.1%), 4 isolates belonged to t223 (4.5%), 3 isolates to t852 (3.4%), and 2 isolates to t005 (2.2%). Both coa type IVb with spa type t421 and coa type VI with spa types t223 (3.4%) and t421 (1.1%) were detected in 4 isolates (4.5%). coa type Iva were identified in 2 isolates with spa type t10795 (2.2%). Two isolates with spa type t064 (2.2%) and 1 isolate with spa type t421 (1.1%) were found in strains with coa type X (3.4%).
PVL-encoding genes were distributed among
coa types I (1 isolate; 1.1%), II (3 isolates; 3.4%), III (2 isolates; 2.2%), and VIII (3 isolates; 3.4%), while
tst was identified among
coa types I (2 isolates; 2.2%), III (9 isolates; 10.1%), VI (3 isolates; 3.4%).
eta and
etb were detected in
spa types t019 (2.2%) and t030 (2.2%), respectively. Surprisingly, all the LLMUPR and HLMUPR-MRSA strains belonged to
coa type III. The characteristics of the 89 MRSA strains isolated from burns patients have been summarized in
Table 1.
Discussion
Several studies have revealed different findings of resistance rate of
S. aureus isolated from burns patients which may be linked to various bacterial detection methods [
10]. Although the prevalence rate of mupirocin-resistant
S. aureus strains depends on various geographic areas and/or patient populations, a relatively high prevalence of mupirocin-resistant
S. aureus strains (28.1%) was found in this study, which was similar to our previously reported rate from burns patients (28.3%) [
1]. These results however, show a different prevalence of mupirocin-resistant
S. aureus strains than reported by Abbasi-Montazeri et al (6%) [
5], Shahsavan et al (40%) [
14] and Saderi et al (2.7%) [
15] in Iran. A higher prevalence rate of mupirocin resistance was detected among the isolates examined in this study, compared with that found in Greece (2%) [
16], Egypt (17.8%) [
17], and France (2.1%) [
18]. Given the many factors affecting mupirocin resistance, the high rate of mupirocin resistance among
S. aureus strains may be related to type and volume of our clinical samples, geographic or socioeconomic factors, infection control policies in healthcare and burn centers, and unrestricted policies in the use of mupirocin. Likewise, the present study showed that the frequency of HLMUPR (13.5%) was lower than LLMUPR (14.6%). The findings of a study in Egypt showed that 61.5% of mupirocin resistance isolates had the HLMUPR pattern, while low-level resistance was reported in 38.5% of mupirocin-resistant isolates [
17].
Antimicrobial susceptibility testing revealed a high prevalence rate of resistance to erythromycin (86.5%), gentamicin (84.2%), tetracycline (65.2%), amikacin (55.1%), and clindamycin (51.7%). This finding is similar to that reported by Ko et al from 12 Asian countries [
19].
S. aureus isolates with resistance to vancomycin have emerged in many parts of the world and their prevalence is increasing [
20,
21], and in accordance with other studies [
1,
5,
13], all the isolates in our study were susceptible to vancomycin.
The predominant resistance genes in this study were
ant (
4′)
-Ia (60.7%),
aac (
6′)
-Ie/aph (
2″) (55.1%), and
aph (
3′)
-IIIa (46.7%). Similarly, other researchers showed that aminoglycoside resistance genes are dominant in
S. aureus clinical isolates [
4,
13,
22]. In addition, in line with other studies, the prevalence of
ant(
4′)
-Ia gene was higher than that of the 2 other aminoglycoside resistance genes,
aac(
6′)/aph(
2′′) and
aph(
3′)
-IIIa [
13]. In a study conducted by Ardic et al [
23] to determine aminoglycoside modifying enzyme genes in MRSA strains isolated in Turkey, 50 staphylococci from in-patients were assessed. They reported the majority of the isolates carrying
aac (
6′)
-Ie/aph (
2″) (60.5%) and followed
ant (
4′)
-Ia (24%) and
aph (
3′)
-IIIa (8%). This increased rate among the isolates tested could be related to the source of the isolated strains (burn wound), inappropriate use of aminoglycoside in the treatment of serious staphylococcal infections, and horizontal gene transfer among the strains.
The
tet(
M) gene was the 4
th most detected antibiotic resistance gene among our isolates, and was detected in 43.8% of the strains. The frequency of
tet(
M) in the study of Rashidi et al on 105
S. aureus strains collected from an intensive care unit was 56.8% [
13]. In the Emaneini et al study, this frequency was reported as 32.4% [
4].
In the present study, 13.5% of isolates carried the
mupA gene and were confirmed as HLMUPR MRSA. This finding is lower than the rates reported by Abbasi-Montazeri et al from Iran (34%) [
5], González-Domínguez et al from Spain (27.2%) [
24], and is higher than India (5%) [
25] and Korea (1.8%) [
26]. In general, variations in resistance genes were not unexpected, as these genes are subjected to a high, but variable selective pressure.
The findings in the present study provide insight into the epidemiology of
coa and
spa types in Iran. There were 89 MRSA strains isolated from burns patients with 10
coa genes identified by multiplex PCR, similar to the results of Afrogh et al [
10], who showed 9 different patterns of
coa genes among staff nose, and patients’ clinical samples in Iran. The same results were reported by Janwithayanuchit et al [
27] from Thailand, and Omar et al [
28] from Egypt. Our study has reported that types III, II, I and VIII were the most common, and accounted for 47.2%, 14.6%, 10.1% and 10.1% of the isolates, respectively. Hirose et al [
8], using the same primer, showed 7 different types of
coa gene. They showed that the majority of MRSA isolates belonged to
coa type II (91.9%), followed by type VII (3.9%), type I (1.7%), type Iva (1.1%), types III and X (0.6%) and type VI (0.2%). There were 2 MRSA isolates that were untypeable. The difference in coagulase types was attributed to geographical variation.
The present study reported the existence of t790 as the predominant
spa type in 18% of isolates which were all
coa type III. The low frequency of t790 in the study was in line with Udo’s study from Kuwait, 2016 [
21], and a study conducted by Goudarzi et al [
1] who analyzed
S. aureus strains isolated from burns patients in Iran (13.2%). More than half of
spa t790 isolates in this study were carrying the
tst gene (56.3%). LLMUPR was detected in 11 strains with t790 (12.4%). Although, virulence markers and drug resistance patterns in t790 isolates may vary, resistance to mupirocin and
tst carriage in t790 isolates has been reported by several investigators [
1,
20].
In the current study, the second most frequent
spa type was t064 (13.5%), distributed in
coa types III (11.2%) and X (2.3%). All these strains were carrying the
mupA gene and confirmed as HLMUPR strains. In agreement with our study, mupirocin resistance in t064 strains has been reported previously from Nigeria [
29], Kuwait [
21] and Ireland [
20].
The third most common
spa type was t030 detected in 11.3% of isolates, that all belonged to
coa type III. This finding contrasts with our previous study in burns patients that reported t030 as one of the most common
spa types amongst the tested isolates (24.5%) [
1]. Similarly, the authors of a study of MRSA in Kuwait reported a low prevalence of
spa type t030 (2.7%) [
21].
Another
spa type in our study was t044 (10.1%; 9 isolates), that of those with 5 of those isolates found to be positive for the
pvl encoding gene (55.6%). Of 9 isolates with the t030
spa type, 4 isolates belonged to
coa type II (44.4%) and 5 isolates to
coa type VIII (55.6%). A review performed by Asadollahi and colleagues in 2018, to determine distribution of the most prevalent
spa types among clinical isolates of methicillin-resistant and -susceptible
S. aureus, indicated that t044, distributed in 11 countries, was one of the most common
spa types [
11]. This
spa type was also previously reported [
30].
In agreement with studies performed in Kuwait [
31], and Saudi Arabia [
32], PVL-positive, multi-resistant
spa type t852 isolates (9%) were found in our survey, pointing to their spread in the Gulf Persian Cooperative Council countries.
In this study, the
spa type t223 was found in 7.9% of isolates for the first time in Iran. The resistance and virulence profiles of the t223 isolates (PVL-negative,
tst-positive) was similar to t223 isolates reported in
S. aureus strains isolated from children and parents in the Gaza Strip [
33], and also t223 recovered from healthy individuals in Jordan [
34], which may suggest the origin of
spa type t223.
The t421 isolates was detected in small numbers, similar to Wang et al [
35], who revealed that out of 99 MRSA isolates classified as ST239, 92 were ST239-
spa t037 and 4 were ST239-
spa t421.
In the present study,
spa type t008 was detected at an incidence of 6.7%. The majority of the t008 MRSA isolates were resistant to mupirocin, as has been reported in t008 isolates obtained in Kuwait, the USA, and many European countries [
20,
21].
The
spa type t019 was observed among 6 isolates (6.7%) with 5 isolates belonging to SC II and 1 isolate to SC VII. Two isolates harbored the
eta encoding gene and belonged to SC II.
spa type t019 has also been reported in the United States, and certain Asian and European countries like Egypt, Japan, Poland, and Taiwan [
20]. The findings of the present study are in line with a multicenter study conducted during a 6-year period in 17 countries, that reported a low frequency of
spa type t019 (3.6%) among tested isolates [
36].
In line with our findings that showed the existence of
spa type t021 in 4.5% of isolates, several researchers in Lebanon, Ireland, Romania, and Portugal reported a low frequency of the t021 among MRSA clinical isolates in comparison to other
spa types, but not to the same extent [
20].
t10795 with SC type IVa accounted for 2.2% of all strains. This was supported by Boswihi’s report from Kuwait [
21] that
spa type t10795, was one of the
spa types associated with the ST772-V, Bengal Bay clone, and was detected in small numbers. The Bengal Bay clone was previously reported in Italy, UK, Kuwait, Malaysia, Saudi Arabia, and UAE [
21,
35,
37].