- Research article
- Open Access
Identification and molecular epidemiology of methicillin resistant Staphylococcus pseudintermedius strains isolated from canine clinical samples in Argentina
BMC Veterinary Researchvolume 15, Article number: 264 (2019)
Staphylococcus pseudintermedius is the leading cause of pyoderma in dogs and the frequent use of antimicrobial treatment is associated to the development of resistance to nearly all classes of antibiotics. Despite S. pseudintermedius significance, our understanding of the molecular mechanism of β-lactam resistance and its genetic diversity remains limited. We aimed to: i) determine the phenotypic resistance profile of methicillin resistant Staphylococcus pseudintermedius (MRSP) isolated from infected dogs in three different veterinary hospitals in Buenos Aires, Argentina; ii) identify the SCCmec elements and resistance genes; and iii) analyze the clonal relationship between isolates and in regard of dominant lineages found in the world.
In addition to the differential levels of β-lactam resistance, MRSP isolates (n = 10) showed resistance to 5–6 families of antibiotics, and were therefore categorized as multidrug-resistant. All the isolates were variant of SCCmec V homologous to S. aureus; additional SCCmecFinder analysis classified five of the genomes as SCCmec type V (5C2&5) with mecA (encodes for PBP2a), mecRI and mecI and all the genes closely related to the reference SCCmec type V S. aureus TSGH17 strain. In the remaining five strains, mecA was present, although other genes associated with SCCmec V including mecR1 and mecI were missing. PBP2a was inducible in low level resistance strains (MRSP 8151), and constitutively expressed in MRSP 8150, suggesting different mecA regulatory mechanisms. MRSP isolates showed significant genetic diversity: eight PFGE clonal types and six multilocus-sequence typing (MLST) sequence types (STs) (339, 649, 919, 920, 921 and 922), including four new STs genetically distinct from STs reported in other geographic areas. Comparative genomics and phylogenetic analyses of the MRSP showed a correlation between the genetic content and the phenotypes, and established the genetic relationship between the isolates.
MRSP could be a threat to animal health due to it concerning level of antimicrobial resistance. Our study highlights genetic and epidemiological aspects of multidrug-resistant MRSP strains from Argentina showing high degree of correlation between the resistance genes and the phenotype of the isolates and, furthermore, they appeared evolutionary closer to major worldwide reported ST68 and ST71.
S. pseudintermedius is an important pathogen in dogs and cats and is sporadically associated with human infections . Over the past decade, methicillin resistant S. pseudintermedius (MRSP) has emerged in different parts of the world and has become one of the most important bacterial pathogens in small-animal-veterinary-medicine [2, 3]. Based on data from the World Health Organization (WHO) and the United Nations in 2016, antimicrobial resistance in humans, like in companion animals, represents a problem for public health.
The β-lactam resistance of MRSP is due to penicillin-binding protein 2a (PBP2a), a protein encoded by the methicillin resistant gene mecA. This gene is known to reside in a mobile genetic element, a staphylococcal cassette chromosome designated SCCmec that contains the mec gene complex, mecA and some additional genes, and the cassette chromosome recombinase (ccr) gene complex, which is responsible for insertion of the SCCmec cassette into the core genome. So far as many as thirteen different structural types of SCCmec have been described in S. aureus based on the different combinations of class of mec complexes according to the presence/absence of regulatory genes and insertion sequences, and ccr allotypes (ccrAB and ccrC) . Eleven main types, subtypes, and variants have already been described in the database of the International Working Group on the Staphylococcal Cassette Chromosome (IWG-SCC) . Some S. aureus and coagulase-negative Staphylococcus isolates carry a mecA homolog, mecC, which has been recently reported carried by SCCmec XI [6, 7]. The classification of SCCmec elements is complex, given that there are composite cassettes and pseudo-SCCmec elements that do not harbor ccr genes . While in S. aureus the structure of SCCmec elements has been shown to be relatively stable, in MRSP the SCCmec elements showed high genetic diversity .
The cefoxitin disk is considered as the main method for methicillin resistance detection in S. aureus, nevertheless this is not an accurate method of screening for methicillin resistance in S. pseudintermedius [9, 10], that must be detected with the oxacillin disk.
Several dominant MRSP lineages have been identified in the world, including ST45, ST68 and ST71 , but the molecular epidemiology of MRSP clones circulating in Argentina has not been examined.
The aims of this study were to determine the phenotypic resistance profile of MRSP, to identify the SCCmec elements and resistance genes, to analyze the clonal relationship between isolates, and to compare these isolates with the dominant lineages found globally.
Antibiotic resistance profiles
All the S. pseudintermedius isolates were considered resistant to oxacillin based on the recently revised Clinical and Laboratory Standards Institute (CLSI) breakpoints of S. pseudintermedius, the presence of mecA gene and a SCCmec element. We identified from the total of 10 MRSP strains two distinct groups with different expression of β-lactam resistance. Six of the ten isolates displayed low-level of oxacillin resistance with minimal inhibitory concentration (MIC) values in the range of 0.5 to 2 mg/L despite being mecA positive, and the rest (4 of 10 isolates) exhibited high-level resistant (≥8 mg/L) (Table 1). All the strains were resistant to oxacillin, penicillin, streptomycin and kanamycin. In addition to the β-lactam resistance gene mecA, all the isolates contained the β-lactamase gene blaZ, the kanamycin and neomycin phosphotransferase gene aph (3′)-III and the streptomycin adenylyl-nucleotidyltransferase gene ant(6)-Ia. Resistance to macrolides, lincosamides, and streptogramins-B seen in eight isolates was due to the methylase gene ermB, and all of them displayed constitutive resistance to clindamycin. The nine trimethoprim/sulfamethoxazole-resistant isolates contained the dihydrofolate reductase gene dfrG, and the three tetracycline resistant isolates carried the tetracycline and minocycline resistance gene tetM. Only one strain was resistant to gentamicin, and it had the aac(6′)-Ie–aph(2′)-Ia gene. Table 2 summarizes the genes associated with resistance and the antimicrobial resistance phenotype displayed by the isolates.
In addition, mutations in the quinolone resistance-determining region (QRDR) of the topoisomerase genes of the seven ciprofloxacin-resistant isolates were found, including amino acid substitution S84 L in the topoisomerase GyrA and S80I in GrlA. Additional amino acid substitutions were identified outside the QRDR of the topoisomerase genes (Table 3), but their role in fluoroquinolone resistance was not determined. The isolate MRSP 8472 has only one amino acid substitution in grlA gene and remained susceptible to ciprofloxacin. No mutations were found in the gyrB and grlB genes in any of the isolates. All the isolates were susceptible to chloramphenicol, rifampicin, nitrofurantoin and linezolid. All the MRSP strains in the present study were resistant to more than three antimicrobial classes and therefore were classified as multidrug-resistant (MDR) .
Differential levels of PBP2a expression in MRSP correlates with their β-lactam resistance
Resistance to methicillin in S. pseudintermedius, as well as in S. aureus is due to the presence of the mecA gene, which encodes penicillin-binding protein 2a (PBP2a); this protein shows reduced affinity for all β-lactam antimicrobials. S. aureus strains that have the mecA gene incorporated in the genetic element SCCmec are considered to be resistant to all β-lactam antimicrobials. Analysis of PBP2a by Western blot was performed by using specific anti-PBP2a antibodies in protein lysates of MRSP strains expressing either high level of oxacillin resistance (MRSP 8150, oxacillin MIC > 256 mg/L) or low (MRSP 8151, oxacillin MIC 2 mg/L) grown without and with sub-inhibitory concentrations of oxacillin and cephalexin at 0.5μg/ml. As shown in Fig. 1, while PBP2a was inducible in strains expressing low levels of resistance (i.e MRSP 8151), it appeared to be constitutively expressed in MRSP 8150, suggesting different regulatory mechanisms in MRSP. Furthermore, the strain MRSP 8150 has both blaI/blaR1 and mecI/mecR1 genes unlike MRSP 8151 that only has blaI/blaR1. These results may indicate that β-lactam expression in MRSP strains could be linked to mecA differential regulation.
Characterization of the SCCmec element in MRSP strains
Characterization of the SCCmec cassette was performed by multiplex PCR showing that all the isolates displayed the same pattern of bands, but different from the control isolates SCCmec I to VI. The pattern shared two bands with the SCCmec type V element corresponding to mecA and ccr complex, but differed in the band corresponding to J1 region suggesting that it could be a variant of SCCmec V (Fig. 2).
Each of the strain genomes was submitted to SCCmecFinder , an in silico web-based bioinformatic tool that identifies and types SCCmec elements. Using this, the SCCmec type V (5C2 and 5) / SCCmec type Vb (5C2 and 5) was identified in five of the genomes, indicating that there was significant homology with S. aureus AB512767.1 (TSGH17), which was used as a reference. Further BLAST analysis between all the genes annotated in AB512767.1 showed that all the SCCmec V genes found in the cassette were present on a single contig in MRSP 8472, while the other isolates that had the genes on more than one contig. The fastq files from nine isolates were mapped against the MRSP 8472 genome to see if any reads mapped to these genes, or to other genes known to be present in the SSCmec cassette. This comparison showed that five of the genomes (MRSP isolates 8150, 8468, 8469, 8472 and 8473) have mecA, mecR1 and mecI, and the majority of the genes that are present in AB512767.1 (Fig. 3). Two of these genomes (MRSP isolates 8150 and 8472) had all of the genes present in AB512767.1 and three of them (MRSP 8468, 8469 and 8473) have mecA, mecR1 and mecI, but were missing the last two genes on the 5′ end. The remaining five genomes (MRSP 8148, 8151, 8470, 8471 and 8474) have mecA, but were missing many of the other genes associated with SCCmec V, including mecR1 and mecI (Fig. 3). These results suggest that MRSP strains harboring mecA are differentiated into two distinct groups.
Population structure analysis
To identify the relatedness of the MRSP strains PFGE was performed. We found that the ten MRSP strains were differentiated in eight clonal types (A to H). Moreover, identical patterns were observed between MRSP 8468 and MRSP 8469 strains (type A) while MRSP 8470 and MRSP 8471 strains displayed type B (Additional file 1: Figure S1). Additionally, MRSP type characterization was performed by using MLST as described in Methods. We found that the MRSP strains were differentiated into six different STs.
Two of them MRSP 8150 and MRSP 8472 belong to ST339 while MRSP 8468 and MRSP 8469 were related to ST649. These ST types have been described and already included in PubMLST database. The six remaining isolates had previously undescribed allelic profiles and were assigned new sequence types by the S. pseudintermedius MLST database curator. MRSP 8148, MRSP 8470 and MRSP 8471 belong to ST919, MRSP 8151 to ST920, MRSP 8473 to ST921 and MRSP 8474 to ST922.
To determine the clonal relationship between the STs detected in this study with those found in the global Pub-MLST S. pseudintermedius database, all the entries available on June 2018 were clustered using the same goeBURST procedure (Fig. 4). The clonal complex (CC) consisted of allelic profiles with five or more allele matches, while singletons were unrelated to any other within the collection. The goeBURST algorithm showed that three of them (STs 919, 920 and 922) are single locus variant with each other, two (STs 649 and 921) are singletons and ST339 is part of a branch located very far from the other isolates. As shown in Fig. 4, none of the isolates in Argentina were related to ST68 or ST71. This data clearly indicates that the MRSP strains are not identical to ST68 and ST71 but evolutionary related.
Comparative genomics and phylogenetic analyses
Analysis of the MRSP genomes in PATRIC showed an average size of 2,344.838 Mb (range 2.550.634 to 2.650.119), with an average 2779 genes annotated (range 2528 to 2647). The genome composition of MRSP was found similar between the different genomes analyzed (Table 4). A phylogenetic tree that included 152 S. pseudintermedius previously sequenced genomes and the 10 strains in this study was generated (data not shown). Nineteen strains from the representative clades in this tree were selected for a new tree that included the 10 genomes (Table 5 and Fig. 5). The phylogenetic analysis shows that the 10 genomes in this study were not monophyletic, with MRSP 8472 and MRSP 8473 separated from the other eight isolates (Fig. 5) that clustered together.
S. pseudintermedius is a canine and feline commensal and opportunistic pathogen, analogous to S. aureus in humans. MRSP has recently emerged in small animals worldwide and represents a serious threat to animal health due to its characteristic multidrug resistance phenotype .
In this study we compared the resistance phenotype presented by ten clinical strains of MRSP with the genotypic analysis including WGS data. The cefoxitin disk test, which has been shown in several studies to be the most reliable predictor of the presence of mecA in both S. aureus and coagulase-negative Staphylococcus, does not identify MRSP isolates of human and veterinary origins [9, 10, 16] and would therefore not have been appropriate for the isolates in this study. Although all the isolates included in the present study were resistant to oxacillin using the current CLSI breakpoints, they could be separated into two groups, with either high or low oxacillin MICs resistance levels. Western blot analyses of PBP2a included the strains MRSP 8150 and MRSP 8151. MRSP 8150 demonstrated a high level of resistance, which was consistent with constitutive membrane levels of PBP2. The levels of PBP2a in MRSP 8151 were only inducible after exposure to oxacillin or cephalexin. The differential expression of β-lactam resistance has been observed in MRSA . In previous studies, we have demonstrated that PBP2a can be co-regulated by both mecI and blaI regulators . The observation that blaI/blaR1 regulatory system, present in all the MRSP strains in this study (Table 3), led us to speculate that blaI/blaR1 is indispensable to mediate mecA regulation in MRSP. In accordance with this observation, constitutive levels of PBP2a seen in MRSP 8150 may be due to defective function of blaI/blaR1 despite the presence of mecR1/mecI, however this hypothesis awaits experimental confirmation and further studies are needed to demonstrate this finding.
In addition to resistance to β-lactams, MRSP isolates showed resistance to other antibiotics such as aminoglycosides, macrolides, lincosamides, tetracyclines, trimethoprim/sulfamethoxazole and fluoroquinolones. Based on these results, all the MRSP that showed resistance to five or six families of antibiotics, were categorized as MDR.
The genes responsible for the resistance to penicillin, tetracycline, erythromycin and trimethoprim/sulfamethoxazole (blaZ, tetM, ermB and dfrG genes, respectively) have also been found to be the predominating genes encoding for resistance in S. pseudintermedius from dogs in other studies [18, 19]. Resistance to aminoglycosides was associated with the adenyl nucleotidyltransferase gene ant(6)-Ia and the phosphotransferase gene aph(3′)-III. The bifunctional acetyltransferase/phosphotransferase gene aac(6′)-Ie–aph(2′)-Ia was present only in one gentamicin resistant strain. Comparable results were found in a similar study performed in a veterinary medical teaching hospital in Texas, where the most common aminoglycoside resistance gene found was aph(3′)-IIIa, followed by aac(6′)/aph(2′′) and ant(4′)-Ia genes . Most of the resistance genes detected in S. pseudintermedius have also been identified in other species of staphylococci. Only the seven isolates phenotypically resistant to fluoroquinolones had mutations in topoisomerase II (gyrA S84 L) and IV (grlA S80I) simultaneously. In a previous study Descloux et al. identified numerous base pair exchanges in the genes gyrA, gyrB, grlA and grlB of S. pseudintermedius resistant to fluoroquinolones . The same substitutions: S84 L and E88G in gyrA as well as S80I and D84N in grlA, were seen in fluorquinolone resistant MRSP isolates from Japan  and two others S84 L and S80R in gyrA identified in eight ciprofloxacin-resistant MRSP isolates from Spain . Our results show excellent correlation between the resistance phenotypes and the resistance genes detected, similarly to reported recently by other authors .
The information gleaned by full genome sequencing of staphylococcal species allows to determine the diversity of SCCmec elements, the structural organization and the genetic content including genetics variants (e.g. inserts, deletions) . SCCmec of S. pseudintermedius displayed some degree of homology to those of S. aureus, but sometimes they are untypeable using SCCmec typing schemes developed for S. aureus. The homology between SCCmec cassettes of different species is considered as an indication of horizontal gene transfer between isolates. SCCmec V is largely homologous to SCCmec type V (5C2&5), previously named VI or VII from S. aureus. Conversely, SCCmec II-III consists in a combination of SCCmec II from S. epidermidis and of SCCmec III from S. aureus and has lacked the cadmium resistance operon, and SCCmec VII-241 is a newly described element that is not related to SCCmec VII from S. aureus . The SCCmec type III variants were found in a study carried out in a veterinary hospital from Japan, where SCCmec type II-III represented 85.2% of S. pseudintermedius isolates . In our study, according to the multiplex PCR method designed by Milheirico et al for S. aureus, all the isolates appear to be a variant of SCCmec V. SCCmecFinder  could only classify five genomes as SCCmec type V(5C2&5)/SCCmec type Vb(5C2&5) due to the fragmented assemblies of the new isolates. A more detailed analysis of the genomes allowed us to observe that these five genomes not only have mecA, mecR1 and mecI, but also have a majority of the genes that are present in the isolate used as reference of SCCmec V. Moreover, two of these genomes had all of the genes present in the S. aureus TSGH17 that was used as a reference. The other five isolates have mecA, but were missing many of the other genes associated with SCCmec V, including mecR1 and mecI. Further studies are warranted to characterize the SCCmec element displayed by these isolates, which have the same ccrC recombinase but appears to be a different combination of genes than those described up to now. Although we observed differences in the SCCmec elements between the isolates, the differences in oxacillin MICs seem not to be associated to different SCCmec types, as has been recently described [8, 28].
The dissemination of MRSP isolates tended to be associated with a limited number of clones, unlike methicillin susceptible S. pseudintermedius isolates that presented great genetic diversity , similarly to the situation observed in human S. aureus. ST68 clone SCCmec V and ST71 SCCmec II-III are the dominant clones that have spread in North America since 2003–04 and in Europe since 2005–06, respectively, but now have a global distribution . A study performed in Brazil was the first in South America to have detected the European clone ST71 of MRSP colonizing companion animals . The isolates of our collection showed significant genetic variation between the population manifested by eight clonal types differentiated by PFGE and six sequence types (STs) by MLST (339, 649, 919, 920, 921 and 922), including four new STs that were genetically distinct from the previous STs in other geographic regions. The analysis by goeBURST of our isolates showed that they were not related to ST68 or ST71. However, ST68 and ST919 are double locus variant from ST677, which could indicate that our isolates are evolutionarily closer to ST68 than ST71. Notably the diversity of STs shown by our MRSP isolates indicates high clonal diversity in our country. We found that both internationally reported as well as previously unreported MRSP STs are present in Argentina. Giving that the clones ST919, ST920, ST921 and ST922 had not been previously reported, it is likely that they represent locally evolved clones.
In summary, this is the first report addressing the phenotypic and genotypic characterization of canine MRSP isolated in Argentina between 2008 and 2011. The ability of MRSP to acquire and maintain resistance genes, and its propensity for horizontal transfer of resistance determinants have shown to represent a potential threat on both the veterinary and Public Health settings.
Bacterial strains and identification
Clinical samples were collected from the infected dogs at three different veterinary hospitals in Buenos Aires, Argentina between 2008 and 2011. Three strains (MRSP 8148, MRSP 8150 and MRSP 8151) were the only methicillin resistant detected in a previously studied strain collection of 28 S. pseudintermedius . The other seven strains were recovered in two laboratories from Buenos Aires city during 2011 and fully characterized at the Antimicrobial Division, INEI-ANLIS “Dr. Carlos G. Malbrán”, Regional Reference Laboratory on Antimicrobial Resistance, Buenos Aires, Argentina The ten MRSP isolates included in the present study were isolated from infections in different body sites (Table 1). Species identification was performed by conventional biochemical tests and confirmed by mass spectrometry MALDI-TOF (Bruker Daltonics Microflex LT, Billerica MA, USA). The isolates were pheno- and genotypically characterized at the Antimicrobial Division, INEI-ANLIS “Dr. Carlos G. Malbrán”.
Susceptibility testing of MRSP isolates
S. pseudintermedius strains were tested by disk diffusion to evaluate their antimicrobial susceptibility to the following antibiotics (disk concentration in brackets): oxacillin (1 μg), cefoxitin (30 μg), penicillin (10 units), erythromycin (15 μg), clindamycin (2 μg), tetracycline (30 μg), chloramphenicol (30 μg), rifampicin (5 μg), gentamicin (10 μg), streptomycin (10 μg), kanamycin (30 μg), trimethoprim-sulfamethoxazole (1.25/23.75 μg), nitrofurantoin (300 μg), linezolid (30 μg) and ciprofloxacin (5 μg). Oxacillin MIC was determined by agar dilution in MH agar + 2% NaCl (CLSI) with a range of antibiotic concentrations from 0.03 to 8 mg/L, and by Etest (bioMérieux, France). All antimicrobial susceptibility tests were carried out according to the CLSI guidelines [32, 33]. S. pseudintermedius strains were categorized as susceptible, intermediate, or resistant, when the applicable breakpoint was available in CLSI documents VET01S-3rd ed., 2015 or M100S-27th ed., 2017. Despite the lack of CLSI-approved interpretative criteria for streptomycin and kanamycin, the isolates for which the inhibition zone was 6 mm were considered as resistant. Isolates were considered as multidrug-resistant when they exhibited resistance to three or more different classes of antimicrobial agents .
PBP2a analysis in MRSP strains
Western blot analysis was used to determine changes in PBP2a levels as previously described . Briefly, membrane proteins (15 μg) will be extracted from MRSP strains growing with/without sub-inhibitory concentrations of oxacillin and cephalexin (CFL) at 0.5 μg/ml in MHB until mid-exponential phase; cell pellets were resuspended in 600 μl of phosphate-buffered saline (PBS), disrupted by adding glass beads and using a FastPrep cell disrupter (MP Biomedicals, Santa Ana, CA, USA); the lysate was centrifuged at 8,000×g for 10 min at 4 °C. The supernatant fraction was centrifuged for an additional 5 min at 8,000×g at 4 °C to remove the beads, and the supernatant transferred to ultracentrifuge tubes and ultracentrifuged at 45,000 rpm for 1 h/4 °C. The membrane pellet was resuspended in PBS, total membrane proteins quantified and stored at − 80 °C. Lysates were separated on 4 to 12% bis-Tris gels, blot transferred onto pure nitrocellulose blotting membranes, and after blocking (5% low-fat milk in PBS), PBP2a was probed with monoclonal anti-PBP2a antibody (Slidex MRSA detection kit; bioMérieux, France).
mecA gene PCR
All the MRSP strains were tested for the presence of the mecA gene by PCR. PCRs were performed as previously described , S. aureus ATCC 43300 and S. aureus ATCC 29213 were used as positive and negative control, respectively.
Identification of SCCmec elements among MRSP strains
MRSP isolates were first screened for typical SCCmec elements by multiplex PCR as previously described . S. aureus collection strains were used as control of each SSCmec type: COL, PER34, BK2464, USA100, ANS46, HU25, USA400, a clinical strain and HDE288 were used as positive control of SCCmec types I, IA, II, III, IV, V and VI respectively . The genome sequence from each of the 10 new isolates was examined in the SSCmecFinder resource  to determine the SSCmec type. In addition, the SSCmec V (GenBank Id AY894416)  nucleotide and individual protein sequences were compared by BLAST  to the new genomes in PATRIC . A careful examination of the region containing the SCCmec V genes and its flanking regions was conducted using the Proteome Comparison and Compare Region View  tools found in PATRIC. A broad examination for the presence or absence of the protein families that contain mecA, mecR1 and mecI genes across all S. pseudintermedius genomes was conducted using PATRIC’s Protein Family Sorter .
As the MRSP 8472 genome had all the genes in the SSCmec element present on a single contig, the reads from the remaining nine genomes were mapped to that genome using PATRIC’s variation service to confirm the presence or absence of the genes in the SCCmec V region. Genes were considered present when reads were present that overlapped both the 5′ and 3′ ends of the genes in MRSP 8472, as well as covering more than 60% of the total length of the gene.
Pulsed-field gel electrophoresis (PFGE) analysis
Chromosomal DNA of the MRSP strains digested with Smal was analyzed by PFGE, as described previously . PFGE was carried out by clamped homogeneous electric field electrophoresis with a CHEF DR III System (Bio-Rad Laboratories, Richmond, CA, USA). PFGE was performed under the following conditions: switch time, 2.0 to 20.0 s and run time, 20 h; temperature 11.3 °C, angle 120° and voltage 6 V/cm. Separated DNA fragments were stained with ethidium bromide and visualized with a UV transilluminator. Banding patterns were evaluated by visual inspection and interpreted according to Tenover criteria . Isolates were considered unrelated when the PFGE patterns differed in seven or more bands, consistent with three or more independent genetic events.
Genomic DNA was extracted by using DNeasy Blood and Tissue Kit (QIAGEN, Valencia, CA, USA) as per manufacturer’s instructions; concentration was measured by QubitTM assay (Invitrogen, Carlsbad, CA, USA). Illumina library preparation was carried out by Nextera XT DNA Library Preparation Kit (Illumina, San Diego, CA, USA). Hi-seq sequencing was carried out in our affiliated Weill Cornell University (New York, NY, USA) institution at the Epigenetics and Genomic Laboratory, using an Illumina HiSeq 2000. Assembly, annotation and analysis of genomes were done through the PATRIC software (https://www.patricbrc.org). The detection of resistance genes was carried out with PATRIC using the available ResFinder (genomicepidemilogy.org) and CARD (Comprehensive Antimicrobial Resistance Database, card.macmaster.ca) databases, the gene content were compared with the phenotype presented by them.
Population structure analysis
Sequence types were determined using MLST software (https://bio.tools/mlst). Sequence types were assigned by comparison with the allele sequences present in the PubMLST database (http://pubmlst.org/spseudintermedius) and isolates with a novel combination of alleles were submitted to the MLST database curator Vincent Perreten (firstname.lastname@example.org). We determine the clonal relationships of the sequence types obtained in this study with entries in the global PubMLST S. pseudintermedius database. All entries available at the time of analysis were clustered using the same goeBURST procedure database (http://www.phyloviz.net/goeburst/).
An initial tree including 152 S. pseudintermedius genomes, including the 10 isolates sequenced in this study, was created to select appropriate genomes to represent the phylogenetic diversity found across the species. Genomes were selected based on clusters identified in this tree, and as a result, nineteen previously sequenced genomes representing these branches were selected to be compared with the 10 Argentinian genomes.
Protein families from genes that were present as a single copy per genome were selected, and 1000 of these Global protein families (PGFams)  were used. Both the protein (amino acid) and gene (nucleotide) sequences were used for each of the selected genes. Protein sequences were aligned using MUSCLE , and the nucleotide coding gene for each was aligned using the Codon_align function of BioPython . A concatenated alignment of all proteins and nucloetides were written to a phylip formatted file, and then a partitions file for RaxML  was generated, describing the alignment in terms of the proteins and then the first, second and third codon positions. Support values were generated using 100 rounds of the “Rapid” bootstrapping option  of RaxML. The resulting newick file was viewed in FigTree .
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
Clinical and Laboratory Standards Institute
International Working Group on the Staphylococcal Cassette Chromosome
Minimal inhibitory concentration
Methicillin-resistant S. aureus
Methicillin-resistant S. pseudintermedius
Pulsed-field gel electrophoresis
Quinolone resistance-determining region
- SCCmec :
Staphylococcal cassette chromosome
S. intermedius group
World Health Organization
Duim B, Verstappen KM, Broens EM, Laarhoven LM, van Duijkeren E, Hordijk J, et al. Changes in the population of methicillin-resistant Staphylococcus pseudintermedius and dissemination of antimicrobial-resistant phenotypes in the Netherlands. J Clin Microbiol. 2016;54:283–8. https://doi.org/10.1128/JCM.01288-15.
Gronthal T, Eklund M, Thomson K, Piiparinen H, Sironen T, Rantala M. Antimicrobial resistance in Staphylococcus pseudintermedius and the molecular epidemiology of methicillin-resistant S. pseudintermedius in small animals in Finland. J Antimicrob Chemother. 2017;72:1021–30. https://doi.org/10.1093/jac/dkw559.
Pomba C, Rantala M, Greko C, Baptiste KE, Catry B, van Duijkeren E, et al. Public health risk of antimicrobial resistance transfer from companion animals. J Antimicrob Chemother. 2016;72:957–68. https://doi.org/10.1093/jac/dkw481.
Miragaia M. Factors contributing to the evolution of mecA-mediated β-lactam resistance in staphylococci: update and new insights from whole genome sequencing (WGS). Front Microbiol. 2018;9:2723. https://doi.org/10.3389/fmicb.2018.02723.
International Working Group on the classification of Staphylococcal Cassette Chromosome Elements (IWG-SCC). Classification of staphylococcal cassette chromosome mec (SCCmec): guidelines for reporting novel SCCmec elements. Antimicrob Agents Chemother. 2009;53:4961–7. https://doi.org/10.1128/AAC.00579-09.
Liu J, Chen D, Peters BM, Li L, Li B, Xu Z, et al. Staphylococcal chromosomal cassettes mec (SCCmec): a mobile genetic element in methicillin-resistant Staphylococcus aureus. Microb Pathog. 2016;101:56–67. https://doi.org/10.1016/j.micpath.2016.10.028.
Shore AC, Deasy EC, Slickers P, Brennan G, O'Connell B, Monecke S, et al. Detection of staphylococcal cassette chromosome mec type XI carrying highly divergent mecA, mecI, mecR1, blaZ, and ccr genes in human clinical isolates of clonal complex 130 methicillin-resistant Staphylococcus aureus. Antimicrob Agents Chemother. 2011;55:3765–73. https://doi.org/10.1128/AAC.00187-11.
Worthing KA, Schwendener S, Perreten V, Saputra S, Coombs GW, Pang S, et al. Characterization of staphylococcal cassette chromosome mec elements from methicillin-resistant Staphylococcus pseudintermedius infections in Australian animals. mSphere. 2018;3:e00491–18. https://doi.org/10.1128/mSphere.00491-18.
Bemis DA, Jones RD, Frank LA, Kania SA. Evaluation of susceptibility test breakpoints used to predict mecA-mediated resistance in Staphylococcus isolated from dogs. J Vet Diagn Investig. 2009;21:53–8. https://doi.org/10.1177/104063870902100108.
Bemis DA, Jones RD, Videla R, Kania SA. Evaluation of cefoxitin disk diffusion breakpoint for detection of methicillin resistance in Staphylococcus pseudintermedius isolates from dogs. J Vet Diagn Investig. 2012;24:964–7. https://doi.org/10.1177/1040638712452112.
Vigo GB, Giacoboni GI, Gagetti PS, Pasteran FG, Corso AC. Antimicrobial resistance and molecular epidemiology of Staphylococcus pseudintermedius strains isolated from dog clinical samples. Rev Argent Microbiol. 2015;47:206–11. https://doi.org/10.1016/j.ram.2015.06.002.
Giacoboni GI, Vinocur F, Fauret N, Grandinetti J, Manzuc P. Detection of Staphylococcus pseudintermedius resistant to methicillin and to other antimicrobials commonly used in canine pyodermias. Analecta Vet. 2017;37:19–24. https://doi.org/10.24215/15142590e012.
Magiorakos AP, Srinivasan A, Carey RB, Carmeli Y, Falagas ME, Giske CG, et al. Multidrug-resistant, extensively drug-resistant and pandrug-resistant bacteria: an international expert proposal for interim standard definitions for acquired resistance. Clin Microbiol Infect. 2012;18:268–81. https://doi.org/10.1111/j.1469-0691.2011.03570.x.
Kaya H, Hasman H, Larsen J, Stegger M, Johannesen TB, Allesøe RL, et al. SCCmecFinder, a web-based tool for typing of staphylococcal cassette chromosome mec in Staphylococcus aureus using whole-genome sequence data. mSphere. 2018;3:e00612–7. https://doi.org/10.1128/mSphere.00612-17.
Paul NC, Moodley A, Ghibaudo G, Guardabassi L. Carriage of methicillin-resistant Staphylococcus pseudintermedius in small animal veterinarians: indirect evidence of zoonotic transmission. Zoonoses Public Health. 2011;58:533–9. https://doi.org/10.1111/j.1863-2378.2011.01398.x.
Wu MT, Burnham CA, Westblade LF, Dien Bard J, Lawhon SD, Wallace MA, et al. Evaluation of oxacillin and cefoxitin disk and MIC breakpoints for prediction of methicillin resistance in human and veterinary isolates of Staphylococcus intermedius group. J Clin Microbiol. 2016;54:535–42. https://doi.org/10.1128/JCM.02864-15.
Rosato AE, Kreiswirth BN, Craig WA, Eisner W, Climo MW, Archer GL. mecA-blaZ corepressors in clinical Staphylococcus aureus isolates. Antimicrob Agents Chemother. 2003;47:1460–3. https://doi.org/10.1128/AAC.47.4.1460-1463.2003.
Kadlec K, Schwarz S. Antimicrobial resistance of Staphylococcus pseudintermedius. Vet Dermatol. 2012;23:276–82. https://doi.org/10.1111/j.1365-3164.2012.01056.x.
Norstrom M, Sunde M, Tharaldsen H, Mørk T, Bergsjø B, Kruse H. Antimicrobial resistance in Staphylococcus pseudintermedius in the Norwegian dog population. Microb Drug Resist. 2009;15:55–9. https://doi.org/10.1089/mdr.2009.0865.
Gold RM, Cohen ND, Lawhon SD. Amikacin resistance in Staphylococcus pseudintermedius isolated from dogs. J Clin Microbiol. 2014;52:3641–6. https://doi.org/10.1128/JCM.01253-14.
Descloux S, Rossano A, Perreten V. Characterization of new staphylococcal cassette chromosome mec (SCCmec) and topoisomerase genes in fluoroquinolone- and methicillin-resistant Staphylococcus pseudintermedius. J Clin Microbiol. 2008;46:1818–23. https://doi.org/10.1128/JCM.02255-07.
Onuma K, Tanabe T, Sato H. Antimicrobial resistance of Staphylococcus pseudintermedius isolates from healthy dogs and dogs affected with pyoderma in Japan. Vet Dermatol. 2012;23:17–22. https://doi.org/10.1111/j.1365-3164.2011.00995.x.
Gomez-Sanz E, Torres C, Lozano C, Sáenz Y, Zarazaga M. Detection and characterization of methicillin-resistant Staphylococcus pseudintermedius in healthy dogs in La Rioja, Spain. Immunol Microbiol Infect Dis. 2011;34:447–53. https://doi.org/10.1016/j.cimid.2011.08.002.
Wegener A, Broens EM, Zomer A, Spaninks M, Wagenaar JA, Duim B. Comparative genomics of phenotypic antimicrobial resistances in methicillin-resistant Staphylococcus pseudintermedius of canine origin. Vet Microbiol. 2018;225:125–31. https://doi.org/10.1016/j.vetmic.2018.09.013.
Perreten V, Chanchaithong P, Prapasarakul N, Rossano A, Blum SE, Elad D, et al. Novel pseudo-staphylococcal cassette chromosome mec element (ψSCCmec57395) in methicillin-resistant Staphylococcus pseudintermedius CC45. Antimicrob Agents Chemother. 2013;57:5509–15. https://doi.org/10.1128/AAC.00738-13.
Perreten V, Kadlec K, Schwarz S, Gronlund Andersson U, Finn M, Greko C, et al. Clonal spread of methicillin-resistant Staphylococcus pseudintermedius in Europe and North America: an international multicentre study. J Antimicrob Chemother. 2010;65:1145–54. https://doi.org/10.1093/jac/dkq078.
Ishihara K, Shimokubo N, Sakagami A, Ueno H, Muramatsu Y, Kadosawa T, et al. Occurrence and molecular characteristics of methicillin-resistant Staphylococcus aureus and methicillin-resistant Staphylococcus in academic veterinary hospital. Apppl Environ Microbiol. 2010;25:5165–74. https://doi.org/10.1128/AEM.02780-09.
Kasai T, Saegusa S, Shirai M, Murakami M, Kato Y. New categories designated as healthcare-associated and community-associated methicillin-resistant Staphylococcus pseudintermedius in dogs. Microbiol Immunol. 2016;60:540–51. https://doi.org/10.1111/1348-0421.12401.
Pires dos Santos T, Damborg P, Moodley A, Guardabassi L. Systematic review on global epidemiology of methicillin-resistant Staphylococcus pseudintermedius: inference of population structure from multilocus sequence typing data. Front Microbiol. 2016;7:1599. https://doi.org/10.3389/fmicb.2016.01599.
McCarthy AJ, Harrison EM, Stanczak-Mrozek K, Leggett B, Waller A, Holmes MA, et al. Genomic insights into the rapid emergence and evolution of MDR in Staphylococcus pseudintermedius. J Antimicrob Chemother. 2015;70:997–1007. https://doi.org/10.1093/jac/dku496.
Quitoco IM, Ramundo MS, Silva-Carvalho MC, Souza RR, Beltrame CO, de Oliveira TF, et al. First report in South America of companion animal colonization by the USA1100 clone of community-acquired meticillin-resistant Staphylococcus aureus (ST30) and by the European clone of methicillin-resistant Staphylococcus pseudintermedius (ST71). BMC Res Notes. 2013;6:336. https://doi.org/10.1186/1756-0500-6-336.
Clinical and Laboratory Standards Institute (CLSI). VET01-S 3rd Edition. Performance standards for antimicrobial disk and dilution susceptibility tests for bacteria isolated from animals. Wayne, PA, USA; 2015.
Clinical Laboratory Standards Institute (CLSI). M100S 27th Edition: performance standards for antimicrobial susceptibility testing. Clinical Laboratory Standards Institute, Wayne, PA, USA. 2017.
Fernandez R, Paz LI, Rosato RR, Rosato AE. Ceftaroline is active against heteroresistant methicillin-resistant Staphylococcus aureus clinical strains despite associated mutational mechanisms and intermediate levels of resistance. Antimicrob Agents Chemother. 2014;58:5736–46. https://doi.org/10.1128/AAC.03019-14.
Vannuffel P, Gigi J, Ezzedine H, Vandercam B, Delmee M, Wauters G, et al. Specific detection of methicillin-resistant Staphylococcus species by multiplex PCR. J Clin Microbiol. 1995;33:2864–7 PMCID: PMC228596 PMID: 8576335.
Milheirico C, Oliveira DC, de Lencastre H. Update to the multiplex PCR strategy for assignment of mec element types in Staphylococcus aureus. Antimicrob Agents Chemother. 2007;51:3374–7. https://doi.org/10.1128/AAC.00275-07.
Boyle-Vavra S, Ereshefsky B, Wang CC, Daum RS. Successful multiresistant community-associated methicillin-resistant Staphylococcus aureus lineage from Taipei, Taiwan, that carries either the novel staphylococcal chromosome cassette mec (SCCmec) type VT or SCCmec type IV. J Clin Microbiol. 2005;43:4719–30. https://doi.org/10.1128/JCM.43.9.4719-4730.2005.
Altschul SF. BLAST algorithm. 2001.
Wattam AR, Davis JJ, Assaf R, Boisvert S, Brettin T, Bun C, et al. Improvements to PATRIC, the all-bacterial bioinformatics database and analysis resource center. Nucleic Acids Res. 2016;45(D1):D535–42. https://doi.org/10.1093/nar/gkw1017.
Overbeek R, Olson R, Pusch GD, Olsen GJ, Davis JJ, Disz T, et al. The SEED and the rapid annotation of microbial genomes using subsystems technology (RAST). Nucleic Acids Res. 2013;42(D1):D206–14. https://doi.org/10.1093/nar/gkt1226.
Wattam AR, Gabbard JL, Shukla M, Sobral BW. Comparative genomic analysis at the PATRIC, a bioinformatic resource center, in Host-Bacteria Interactions Methods. Mol Biol. 2014;1197:287–308. https://doi.org/10.1007/978-1-4939-1261-2_17.
Chung M, de Lencastre H, Matthews P, Tomasz A. The multilaboratory project collaborators. Molecular typing of methicillin resistant Staphylococcus aureus (MRSA) by pulsed field gel electrophoresis: comparison of results obtained in a multilaboratory effort using identical protocols and MRSA strains. Microb Drug Resis. 2000;6:189–98. https://doi.org/10.1089/mdr.2000.6.189.
Tenover FC, Arbeit RD, Goering RV, Mickelsen PA, Murray BE, Persing DH, et al. Interpreting chromosomal DNA restriction patterns produced by pulsed-field gel electrophoresis: criteria for bacterial strain typing. J Clin Microbiol. 1995;33:2233–9 PMCID: PMC228385 PMID: 7494007.
Davis JJ, Gerdes S, Olsen GJ, Olson R, Pusch GD, Shukla M, et al. PATtyFams: protein families for the microbial genomes in the PATRIC database. Front Microbiol. 2016;7:118. https://doi.org/10.3389/fmicb.2016.00118.
Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32:1792–7. https://doi.org/10.1093/nar/gkh340.
Cock PJ, Antao T, Chang JT, Chapman BA, Cox CJ, Dalke A, et al. Biopython: freely available Python tools for computational molecular biology and bioinformatics. Bioinformatics. 2009;25:1422–3. https://doi.org/10.1093/bioinformatics/btp163.
Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014;30:1312–3. https://doi.org/10.1093/bioinformatics/btu033.
Stamatakis AP, Hoover P, Rougemont J. A rapid bootstrap algorithm for the RAxML web servers. Syst Biol. 2008;57:758–71. https://doi.org/10.1080/10635150802429642.
Rambaut A. FigTree, a graphical viewer of phylogenetic trees. 2007.
We acknowledge the Epigenomics Core Facility at Weill Cornell Medicine (New York, NY, United States) for their resources and assistance with whole genome sequencing experiments. We thank to Ms. Anna Brock and Dr. Roberto Rosato for contributing to the critical reading and editing of the manuscript.
This work was supported by the regular federal budget of the Ministry of Health of Argentina and NIH-R56 AI118756 grant to AER. The sponsors of this study have not being involved in the study design, collection, analysis, data interpretation writing of the manuscript, or in the decision to submit the article for publication.
Ethics approval and consent to participate
This study was conducted with MRSP strains that were referred to our laboratory with the purpose of identifying their mechanisms of resistance. We did not have any involvement in the samples collection, handling and treatment. For the previous clinical purposes, the samples were collected in accordance with the guidelines for the care and use of animals of the veterinary hospitals of Argentina and informed consent was obtained from dogs owners.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. PFGE of 10 MRSP strains digested with SmaI. Lane 1: MW marker; lane 2: MRSP 8148; lane 3: MRSP 8150; lane 4: MRSP 8151; lane 5: MRSP 8468; lane 6: MRSP 8469; lane 7: MRSP 8470; lane 8: MRSP 8471; lane 9: MRSP 8472; lane 10: MRSP 8473; lane 11: MRSP 8474. (TIF 463 kb)