- 1Department of Veterinary Sciences, Universidade Federal de Uberlândia, Uberlândia, Brazil
- 2Graduate Program in Genetics and Biochemistry, Universidade Federal de Uberlândia, Uberlândia, Brazil
- 3Department of Preventive Veterinary Medicine and Animal Health, Universidade de São Paulo, São Paulo, Brazil
- 4Department of Veterinary Sciences, Universidade Federal de Minas Gerais, Belo Horizonte, Brazil
- 5Institute of Veterinary Medicine, Georg-August-Universitat, Göttingen, Germany
- 6Department of Veterinary Sciences, Universidade Federal do Vale do São Francisco, Petrolina, Brazil
The study’s objective was to isolate Staphylococcus pseudintermedius from dogs with superficial pyoderma and/or recurrent otitis, and their guardians, to determine oxacillin, cefovecin, and gentamicin resistance, associated risk factors for infection, and genetic similarity between isolates from dogs and guardians. Prevalence of S. pseudintermedius in dogs was 76% and in humans was 56%, with concurrent identification in 44%. Oxacillin resistance occurred in 6.38% of dog isolates and 11.36% of isolates from guardians, with relatively strong (disk diffusion test) or moderate (Minimum Inhibitory Concentration) association between phenotypic testing and mecA gene presence. For cefovecin, dog isolates presented resistance in 8.15% disk diffusion and 23.40% broth microdilution. In humans, 6.81 and 36.36% showed cefovecin resistance in disk diffusion and Minimum Inhibitory Concentration tests. Gentamicin resistance in dogs was identified in broth microdilution testing in 2.12%. fnbB was identified in 4.39% of isolates with relatively strong association between results of dogs and humans. Dogs sleeping with S. pseudintermedius-positive humans were 6 times more likely to test positive, and dogs attending grooming sessions were 4 times more likely positive for S. pseudintermedius. The whole genome sequencing revealed transmission between dogs and humans in 3 cases. Resistance genes for 8 antibiotic classes were found in chromosomal and plasmid genomes.
1 Introduction
Staphylococcus pseudintermedius (SP), along with two other species, S. intermedius and S. delphini, both coagulase-positive, form the group known as the Staphylococcus intermedius group (SIG) (1). Like S. aureus, SP can produce a variety of virulence factors, including coagulase, thermonuclease, proteases, surface proteins such as clumping factors and protein A, as well as cytotoxins, exfoliative toxins, and enterotoxins (2).
SP is the most frequently isolated species in cases of pyoderma in dogs, often identified in ears and infected wounds, and can complicate immunomodulatory responses (3, 4). Another critical aspect regarding SP is its ability to develop resistance to multiple drugs. In this context, Methicillin-resistant Staphylococcus pseudintermedius (MRSP), characterized by the presence of the mecA gene, has emerged as a significant challenge for small animal veterinarians, especially veterinary dermatologists, due to its extensive drug resistance and capacity to cause nosocomial infections (5). MRSP is increasingly isolated from dogs and has also been identified as a causative agent of infections in humans (6–9).
Other genes are associated with antimicrobial resistance in Staphylococcus species and may complicate treatment in cases of infection, such as the msr(A) gene, which encodes a macrolide efflux protein (10). Additionally, other mechanisms may be involved in β-lactam resistance, such as the production of β-lactamases. In this context, the blaZ gene, whose expression is regulated by the blaI and blaR1 genes, is frequently reported (11–13).
Regarding its zoonotic potential, reports of MRSP infections in humans are increasing worldwide, with isolates found in wounds from bites or other skin infections, as well as external otitis, sinusitis, catheter-related bacteremia, pneumonia, and brain abscess (14–16). Furthermore, it has been demonstrated that guardians of dogs with deep pyoderma are often carriers of SP (17), and studies suggest transmission of MRSP from dogs to their guardians (18, 19) and to veterinary hospital employees (20).
The whole genome sequencing (WGS) has been employed in the investigation of SP in animals and humans, contributing to the characterization of virulent strains and aiding in the tracking of multidrug-resistant strains (21, 22). However, this technology is not universally available, as countries with lower and middle development have limitations in accessing such tools. Few studies in Brazil have utilized these technologies in the investigation of SP in dogs (23, 24). Therefore, the objective was to isolate, characterize and correlate risk factors for infection of SP and associate the isolates from dogs with superficial pyoderma and/or recurrent otitis and their guardians.
2 Methods
2.1 Samples from the dogs and their guardians
A cross-sectional study was conducted in which clinical samples from 50 animals (Supplementary Table S0) with otitis and/or pyoderma treated at public and private veterinary hospitals, as well as from their guardians, were evaluated to determine the simultaneous presence of positivity, hypothetical risk factors, and prevalence. The criteria for the sample size were limited to the time available for data collection within the study.
The inclusion criteria for the study were dogs that had recurrent superficial otitis and/or pyoderma on at least three separate occasions, regardless of previous use of topical and/or systemic antimicrobial medications. As an exclusion criterion, only animals that had not received antibiotics in the previous 15 days were selected. Animals and/or humans who were using antimicrobials at the time of clinical evaluation were excluded from the study.
For the identification of superficial pyoderma lesions, we considered the presence of circular alopecia, papules, pustules, erythema, and epidermal collarettes. To characterize otitis, we considered the presence of otological itching or pain, as well as erythema, edema, and any type of ear discharge.
The guardians also responded to a questionnaire regarding the dogs’ lifestyle habits, interactions with humans, progression of the animals’ disease, and treatments already administered (Supplementary Table S1).
Material from the external ears was collected using a veterinary otoscope and sterile swabs, which were then stored in Stuart medium (Absorve®). Skin samples were collected similarly, preferably from intact pustules (previously punctured with a hypodermic needle) or epidermal collarettes at the lower portion of vesicobullous or pustular lesions. Nasal samples from guardians were self-collected using swabs under research supervision. Following collection, samples were immediately transported to the laboratory for culture and identification.
This study was approved by the Animal Ethics Committee (CEUA) under protocol number 055/18, and by the Human Research Ethics Committee (CEP) under protocol number 2,980,894.
2.2 Culture and identification of S. pseudintermedius isolates
The samples were enriched in Brain Heart Infusion broth (Oxoid®) for 12 h at 37 °C and subsequently streaked onto Baird-Parker agar (Oxoid®) with incubation at 37 °C for 24 h. Three colonies showing typical black appearance, with or without a halo indicative of lecithinase activity, were subcultured onto Baird-Parker agar (Oxoid®) and further incubated at 37 °C for 24 h. Pure colonies were phenotypically tested for coagulase and catalase enzyme production and morphologically identified by Gram staining (25). Colonies displaying phenotypic characteristics of coagulase-positive Staphylococcus were genotypically analyzed using PCR for SP identification, and simultaneously preserved in UHT milk (Ninho®, Nestlé) at −80 °C in an ultrafreezer.
DNA extraction was performed using a heat-based method. A loopful of colonies grown on Baird-Parker agar (Oxoid®) was transferred into 100 μL of TE buffer (10 mM Tris HCl, 5 mM EDTA, pH 8.0) in a 1.5 mL microcentrifuge tube. The mixture was heated at 90 °C for 15 min to release DNA into the buffer. The tube was then frozen completely for approximately 30 min and thawed at room temperature. After thawing, the mixture was vortexed to homogenize and centrifuged at 14,000 rpm for 5 min. The supernatant containing the DNA was carefully collected using a 20 μL pipette and stored at −20 °C.
For the identification of SP, primers previously described by Sasaki et al. (26) targeting the nuc gene were used (Forward primer: F2 TRGGCAGTAGGATTCGTTAA; Reverse primer: R5 CTTTTGTGCTYCMTTTTGG). The PCR reactions were prepared using 12.5 μL of GoTaq® Green Master Mix, 2.5 μL of each primer (10 pmol/μL), 7 μL of GoTaq® Green water, and 3 μL of total DNA. The amplification protocol included an initial denaturation cycle at 95 °C for 2 min, followed by 35 cycles of denaturation at 95 °C for 30 s, annealing at 56 °C for 35 s, and extension at 72 °C for 1 min. A final extension step was performed at 72 °C for 2 min.
Subsequently, agarose gel electrophoresis was conducted using a 1.5% agarose gel with a 100 bp molecular weight marker at 80 V, 80 mA, and 100 W for 120 min (Thermo Scientific®). After gel electrophoresis, gel imaging was performed using a UV transilluminator (L-PIX EX®). Presence of a band at 926 base pairs (926 bp) indicated the sample as SP (26).
The MRSP 3279 strain of SP served as positive control, while an amplification reaction without added DNA served as a negative control. We selected only samples where SP was identified in both dogs and humans for analysis.
2.3 Virulence and antimicrobial specific genes
The same previously extracted DNA was used, and amplifications were carried out using 12.5 μL of GoTaq® Green Master Mix, 2.5 μL of each primer for the respective genes (10 pmol/μL) (Table 1), 7 μL of GoTaq® Green water, and 3 μL of total DNA.
Table 1. Sequence and base pairs of the genes investigated in the isolates of SP from dogs and their guardians.
The amplification was conducted under the following conditions: an initial denaturation cycle at 95 °C for 5 min; followed by 35 amplification cycles (denaturation at 94 °C for 2 min, annealing at 53 °C for 2 min, and extension at 72 °C for 1 min); and a final extension at 72 °C for 7 min. After gel electrophoresis, gel imaging was performed using a UV transilluminator (L-PIX EX®). The presence of each gene was determined by the size of the bands in base pairs as per Table 1. The MRSP 3279 strain of SP was used as positive control, and a reaction without added DNA served as a negative control.
2.4 Antimicrobial susceptibility (disk diffusion test)
SP (from dogs and humans) were thawed and reactivated in Brain Heart Infusion broth (Oxoid®) for 12 h at 37 °C. Subsequently, they were streaked onto Tryptic Soy Agar (TSA, BD®) plates and incubated for 24 h at 37 °C. Three colonies were selected using a sterile loop, suspended in sterile 0.9% NaCl solution until reaching turbidity equivalent to 0.5 McFarland standard (27), and then streaked onto Muller-Hinton agar (Oxoid®).
After the agar drying, antimicrobial disks containing oxacillin (1 μg – Cefar®), gentamicin (10 μg – Cefar®), and cefovecin (30 μg – Oxoid®) were applied. After 18 h of incubation at 35 °C, the inhibition zones were measured and compared against reference values provided by CLSI (27). The resistance zone diameter breakpoints for SP were defined as: gentamicin ≤12, oxacillin ≤17, and cefovecin ≤20 (Table 2).
2.5 Minimum inhibitory concentration
The same reactivated and adjusted cultures used for the disk diffusion test were employed for the determination of Minimum Inhibitory Concentration (MIC), using the antimicrobials gentamicin (Novafarma®), oxacillin (Novafarma®), and cefovecin (Convenia®, Zoetis). Microdilution method was used in Mueller-Hinton broth (Oxoid®) adjusted with 0.37 g/100 mL of Ca2+ and 0.835 g/100 mL of Mg2+ (27). Concentrations tested were 0.25, 0.5, 1, 2, 4, 8, 16, 32 μg/mL for the different antibiotics. The resistance breakpoints for SP were defined as: gentamicin ≥16, oxacillin ≥0.5, and cefovecin ≥2 (27, 28).
For analysis, 180 μL of adjusted MH broth (Oxoid®) was added into 96-well microdilution plates containing the target antimicrobials at previously established concentrations. Subsequently, 20 μL of adjusted bacterial suspension was added (28). The plates were then incubated at 35 °C for 18 h for gentamicin and cefovecin samples, and for 24 h for oxacillin (27, 28). After incubation, visual reading was conducted, defining MIC as the lowest antibiotic concentration at which no visible bacterial growth was observed, indicated by the absence of turbidity in the medium (Table 2).
2.6 DNA extraction and sequencing
The DNA was extracted using CTAB (cetyltrimethylammonium bromide) buffer (Sigma-Aldrich). Colonies of SP grown in one petri plate were resuspended in 10 mL of PBS. Then centrifuged at 12.000 g, at 25 °C for 20 min. A total of 500uL of CTAB buffer (100 mM Tris–HCl pH 8.0. 20 mM sodium EDTA pH 8.0; 1.4 M NaCl; 2% CTAB) and 10uL of proteinase K was added in the pellet which was incubated at 65 °C per 60 min and homogenized every 10 min. After, the sample was centrifuged at 12.000 g, at 25 °C for 10 min. The pellet was discarded, and 500uL of a solution chloroform/isoamyl alcohol (24v:1v) was added, homogenized and incubated for 5 min at 25 °C. The supernatant was centrifuged at 6000 g for 10 min at 25 °C and cold isopropanol added (0.65x the volume of supernatant). Then it was incubated at −20 °C for 1 h and centrifuged at 12.000 g for 30 min at 25 °C. The supernatant was discarded, the pellets washed with cold ethanol 70% and centrifuged at 12000 g for 10 min at 25 °C. Finally, the supernatant was dry, eluted in water and lyophilized until the use.
A 450 bp read library was constructed using NEBNext Fast DNA Fragmentation and Library Preparation Kit (New England Biolabs, Ipswich, NE, United States) and quality checked by the Agilent 2,100 Bioanalyzer. The gDNA library was then sequenced by Illumina HiSeq 2,500 sequencing platform (Illumina, San Diego, CA, United States), 2 × 150 bp paired end, insert size 500 bp. Read quality was assessed using FastQC (29) and read trimming was performed using Fastp v. 0.23.4 (30).
2.7 Genome assembly, annotation, pan-genome and phylogenomic analysis
The genomes were assembled using SPAdes v. 3.15.5 (31). Contigs with sizes smaller than 200 bp were removed using Seqkit v. 2.8.0 (32). Completeness and contamination were evaluated using CheckM2 v. 1.0.2 (33) while fragmentation was evaluated using QUAST v. 5.0.2 (34). Contaminant contigs were removed using BV-BRC’s MetagenomicBinning (35). Taxonomy was checked using GTDB-Tk v. 2.3.2 with database r220 (36) and comparison against the reference genome S. pseudintermedius SP_11304-3A (GCA_016126715.1), using a threshold of >95% Average Nucleotide Identity (ANI) calculated with FastANI v. 1.32 (37). Additionally, genome-wide ANI were computed for all pairwise comparisons of genomes using FastANI.
Genome annotation was performed using Prokka v. 1.14.6 (38). To characterize the entirety of genes in the dataset (pan-genome), Panaroo v. 1.2.7 (39) was employed with the strict option to ensure inclusion of only high-quality genes. Nucleotide sequences of individual gene families were aligned using Prank (40). Core genes were defined as those present in ≥ 99% of the genomes, while accessory genes were those present in ≥ 0% and < 99% of the genomes. Single Nucleotide Polymorphisms (SNPs) were extracted from the alignments generated by Panaroo using snp-sites v. 2.5.1 (41). Subsequently, a maximum likelihood phylogenetic tree using RAxML v. 8.2.12 (42) was constructed with a general time reversible (GTR) model (43) and a Gamma model of nucleotide substitution. The phylogenetic tree was visualized using the Interactive Tree of Life (iTOL) online platform (44). For pairwise comparison of genomes, genetic distances were calculated based on SNPs from the concatenated alignment of core genes using snp-dists v. 0.8.2.1
2.8 In silico identification of sequence type, antimicrobial resistance and virulence genes, and plasmid reconstruction
The sequence types (ST) of each isolate were determined using MLST v.2.19.0,2 which extracts seven single copy housekeeping genes (ack, cpn60, fdh, pta, pura, sar and tuf) and compares their sequence identity to previously deposited allele combinations in the S. pseudintermedius PubMLST database.3 All newly identified Sequence Types (STs) were submitted to pubMLST (Supplementary Table S3).
ABRicate v.1.0.04 was used to determine the presence of virulence genes using threshold values of >80% sequence identity and >80% sequence coverage. Virulence genes from the isolates were compared to those found on the Virulence Factor Database (VFDB) (45). AMRFinderPlus (46) was also used to determine the presence of antibiotic and biocide resistance genes, using threshold values of >80% sequence identity and >80% sequence coverage.
MOB-suite v3.1.4 (47) was used to predict plasmid sequences from the assembly genome and identify their replicon types and mobility. AMRfinderPlus (46) was then employed to identify antibiotic and biocide resistance genes in the plasmid sequences, with threshold values of >80% sequence identity and >80% sequence coverage.
2.9 Statistical analysis
The data was tabulated and subjected to descriptive statistics. Calculations were performed using GraphPad Prism 8.3. To assess the association between positivity in dogs and their guardians and the influence of risk factors “bathing location” and “sleeping location of the animal, the Fisher’s exact (p < 0.05) was used, followed by calculation of the Odds Ratio (OR). Comparisons related to gender, age, and race were made using Fisher’s test (p < 0.05) (48).
To determine if there was an association between the presence of resistance and virulence genes in dogs and humans, the Phi (ϕ) coefficient test (p < 0.05) was employed. The strength of association was classified according the ϕ coefficient value (49, 50): very strong (0.8–1), strong (0.6–0.8), relatively strong (0.4–0.6), moderate (0.2–0.4), weak (0.1–0.2), insignificant or nonexistent (0–0.1) (49).
The association between susceptibility results for oxacillin and cefovecin by disk diffusion method and Minimum Inhibitory Concentration (gold standard) was evaluated using the Fisher’s exact test (p < 0.05). Susceptibility, specificity, positive predictive value, negative predictive value, and likelihood odds were also calculated (48).
All 21 sequenced strains had negative results for antimicrobial resistance using the disk diffusion method. The concordance between genotypic and phenotypic resistance, as determined by MIC, was assessed for oxacillin, cefovecin, and gentamicin. For the evaluation of diagnostic performance, phenotypic resistance was considered the gold standard, with genotypic results being compared to this reference.
Associations were evaluated using Fisher’s exact test (two-tailed, p < 0.05) for 21 sequenced strains. Diagnostic performance metrics included susceptibility, specificity, PPV, NPV, Odds Ratio, true positives and true negatives. For gentamicin, all three aminoglycoside resistance genes (aac(6′)-Ie/aph(2″)-Ia, ANT(6)-Ia/aad(6), and aph(3′)-IIIa) were analyzed independently, while β-lactam resistance was evaluated through the complete blaIZR1 operon. Statistical analyses were performed using Rstudio 2024.0.9.1 (51) with the tidyverse (52) epiR (53) and caret (54) packages.
2.10 Ata availability
The genomes of the 21 SPs are isolated from dogs and their guardians that have been sequenced are available on the National Center for Biotechnology Information (NCBI) database, bioproject PRJNA1135159 (Supplementary Table S3).
3 Results
3.1 Distribution of sample profiles according to descriptive and microbiological characteristics
The animals were grouped according to characteristics related to isolated strains, antimicrobial sensitivity, and the presence of resistance genes. Considering the 50 animals used in the study, females accounted for 26/50 (52%) and males for 24/50 (48%) of the sample. The most predominant breed was Shih Tzu 17/50 (p = 0.0338 - Fisher’s exact test), followed by mixed-breed 7/50 (14%), Bulldog 5/50 (10%), Pug 4/50 (8%), Golden 4/50 (8%), Maltese 4/50 (8%), Yorkshire 3/50 (6%), Cocker Spaniel 2/50 (4%), and others 4/50 (8%). The age range was 1 to 5 years 32/50 (p = 0.0012 - Fisher’s exact test), followed by 6 to 10 years 15/50 (30%) and 11 to 15 years 3/50 (6%) (Supplementary Table S0).
3.2 The occurrence of SP in dogs and their guardians was high, and the resistance levels varied according to the phenotypic and genotypic methods employed
The 100 clinical samples from dogs and their guardians yielded a total of 289 colonies, which were tested for the presence of SP, with one sample from each owner’s nostril, one skin and/or ear sample from each dog, and two samples from animals with simultaneous otitis and pyoderma, with three colonies selected per animal sample. Among colonies from dogs and humans, 96/153 (62.94%) and 96/136 (70.59%), respectively, tested positive for SP. The prevalence of SP isolates in dogs was 38/50 (76.00%) and in humans it was 28/50 (56%). Simultaneous identification in samples from dogs and their guardians occurred in 22/50 (44.00%) of cases. There was a positive association between the presence of SP in the dog and in its owner, with humans having five times greater odds of acquiring this bacterium when in contact with a positive dog (p < 0.025, OR = 5).
From the 22 dogs positive for SP on the skin or ear, 47 strains were obtained, among which 11 (23.40%) were resistant to oxacillin by disk diffusion testing (Supplementary Table S2). Among these, 3/11 (27.27%) harbored the mecA gene. Using MIC, 14/47 (29.78%) were resistant to oxacillin, and of these, 4/14 (28.57%) strains were positive for the mecA gene.
Eight out of the 11 strains (72.72%) showing oxacillin resistance in the disk diffusion test and 10/14 (71.42%) in MIC did not carry the mecA gene. Conversely, 2/36 (5.56%) samples susceptibility to oxacillin in the disk test and 1/33 (3.03%) in MIC were positive for the mecA gene.
Among the 44 strains isolated from 22 tutors, 7 (15.90%) were resistant to oxacillin in the disk diffusion test, with 3/7 (42.85%) positive for the mecA gene. In MIC, 23/44 (52.27%) were resistant, and 5/23 (21.74%) carried the mecA gene. Among the resistant strains, four in the disk diffusion test and 17 in MIC were negative for the mecA gene. Meanwhile, 2/44 (4.5%) strains susceptible to oxacillin in both in vitro tests were positive for the mecA gene.
Only 5/47 (10.64%) SP strains isolated from dogs were positive for any mecA gene, with two strains positive for mecA and mecA1 genes, two for mecA, mecA1, and mecA2 genes, and one positive for the mecA2 gene. From the tutor-isolated strains, 5/44 (11.36%) were positive for any mecA gene, with three positives for mecA and mecA1 genes, one for mecA, mecA1, and mecA2 genes, and one for the mecA1 gene. In two strains susceptibility in disk diffusion and resistant in MIC, the mecA gene was identified.
Among SP strains isolated from dogs, 3/47 (6.38%) showed resistance to cefovecin in the disk diffusion test and 9/47 (19.14%) in MIC. In the 44 tutor-isolated strains, 3 (6.81%) and 16 (36.36%) strains were resistant in the disk diffusion and MIC tests, respectively.
For oxacillin-resistant strains in phenotypic tests (set of isolates from dogs and guardians), there was a relatively strong association between oxacillin resistance in the disk diffusion test and the presence of the mecA gene (ϕ coefficient = 0.430), and this association was moderate for MIC (ϕ coefficient = 0.294). In other words, in vitro oxacillin resistance was not an important or sole parameter for identifying the mecA gene.
In 42/47 (89.36%) dog isolates, previous use of some antimicrobial was reported, with 29/42 (69.04%) belonging to the β-lactam class. There was a strong association between previous use of β-lactam antimicrobials and resistance in both disk susceptibility testing (ϕ coefficient 0.682) and MIC (ϕ coefficient = 0.694) for oxacillin.
None of the SP isolates from dogs were resistant to gentamicin in the disk diffusion test; however, in the MIC test, 1/47 (2.12%) demonstrated resistance. Among the isolates from the humans, 2/44 (4.55%) showed resistance to gentamicin in the disk diffusion test, and 4/44 (9.09%) in the MIC (Supplementary Table S2). There was no association between the MIC values of SP isolated from dogs and their tutors for gentamicin.
The msrA gene was identified in 2/47 (4.25%) of the dog isolates and in 2/44 (4.54%) of the tutor isolates. There was no association between the presence of the msrA genes in dogs and their guardians.
The indices obtained from the disk diffusion test when compared to the gold standard (MIC) are described in Table 2.
3.3 The occurrence of the fnbB gene was low, and no samples exhibited gene related to Panton-valentine leucocidin (PVL) production
The fnbB gene was identified in only 4/91 (4.39%) isolates, with 3 (75%) from dogs and 1 (25%) from a tutor, showing relatively strong association between the presence of this gene in isolates from dogs and tutors (ϕ coefficient = 0.563). The PVL-producing gene was not identified in any SP isolates from samples of dogs or their tutors.
3.4 There was a low association between lifestyle habits and the presence of SP in dogs and their tutors
The association between where the dog sleeps (whether with the owner or in its own bed) and the isolation of SP in dogs and tutors did not represent a risk factor for colonization in humans (p = 0.3869, OR = 1.857). However, we observed that dogs that sleep with SP-positive humans are six times more likely to become SP-positive themselves compared to those sleeping in their own bed (p = 0.0264, OR = 6).
Regarding the location where the dog is bathed, we observed that this is not a factor influencing SP positivity in tutors (p = 0.7734, OR = 0.7618). However, despite the association not being statistically significant (p = 0.0519), the p value showed a tendence, and dogs bathed at grooming salons have a four-fold higher chance (OR = 4.091) of being colonized by SP compared to those bathed at home.
3.5 WGS showed a diverse population, with similarity in some SP genomes isolated from dogs and their guardians
Strains were kept frozen until DNA extraction for WGS. When we reactivated the strains, only 21 were viable for DNA extraction. Then, we sequenced the genomes of 21 SP isolates collected from humans (11/21–52.38%) and their dogs (10/21–47.62%). The genomes contained between 23 and 186 contigs, N50 values range between 45,986 and 590,378 bp, and the sequences sizes range from 2.47 to 2.96 Mb (mean 2.60 Mb) (Supplementary Table S3). The ANI values of each genome compared to the reference genome ranged from 98.91 to 99.27%. We identified 3,351 genes families in the pangenome of all SP genomes. A total of 2,045/3,351(61.03%) genes comprised the core genome (present in ≥ 99% genomes), and the number of accessory genes per genome ranged from 233 to 826. To examine the phylogenetic relationships of the isolates, we constructed a maximum likelihood tree based on SNPs extracted from nucleotide alignments of 2,045 core genes (Figure 1a).
Figure 1. Genomic relationships, ANI, and SNP distances among Staphylococcus pseudintermedius isolates from dogs and their guardians. (a) Core-genome phylogeny of S. pseudintermedius isolates showing sequence types (STs) and the presence/absence of antimicrobial resistance (AMR) and virulence genes. Colors indicate host source (dogs, guardians) and STs. (b) Heatmap of average nucleotide identity (ANI, %) among isolates, with dendrogram clustering to highlight genomic similarity within and between hosts. (c) Core single nucleotide polymorphism (SNP) distance matrix with hierarchical clustering, illustrating genetic diversity and potential host-sharing relationships.
Multi-locus sequence typing (MLST) analysis allowed us to identify only 1 STs previously described (C019B), ST 1005. We identified 10 new STs, which were submitted to pubMLST database (Supplementary Table S3; Figure 1a).
We identified resistance genes to 8 classes of antibiotics and 1 heavy metal class (Supplementary Table S3; Figure 1a). One genome isolated from a dog showed at least one resistance gene related to all classes (C21), and no genome was negative for all genes. The most frequent genes were those characterizing resistance to beta lactams, blaI, blaR1 and blaZ, and tetracyclines, tetM, which were present in 18/21 (85.71%) of the samples. Resistance to trimethoprim, defined by the dfrG gene, was the fourth most frequent, present in 13/21 (61.90%) of the samples. About quinolones, two points of mutations were found in the grlA (S80I) and gyrA (S84L) genes, both in 9/21 (42.86%) of the genomes. Resistance to aminoglycosides was defined by the genes aac(6′)-Ie-aph(2″)-Ia (4/21–19.05%), aph(3′)-IIIa (4/21–19.05%) and anti(6)-Ia/aad(6) (4/21–19.05%). Regarding macrolide resistance, the mefE (4/21–19.05%) and ermB (4/21–19.05%) genes were found. The sat4 gene, related to streptothricin resistance, was identified in 4/21 (19.05%) of the samples. One SP genome (4.76%) isolated from a dog had the catA gene, which is related to resistance to phenocols. Finally, qacG gene, related to resistance to quaternary ammonium compounds, was present in 15/21 (71.43%) genomes. The mecA gene was not identified in the genomes by AMRfinderplus, and it was also not detected by PCR in the sequenced samples. The msrA gene was also not found in the genomes, but it was identified by PCR in one guardian SP genome (T019).
The analysis of the relationship between resistance genes and phenotypic profiles (Supplementary Table S4) demonstrated that the aac(6′)-Ie/aph(2″)-Ia gene showed modest performance in predicting gentamicin resistance, with a susceptibility of 33.3%, specificity of 83.3%, and an odds ratio of 2.37, although without statistical significance. In contrast, the ANT(6)-Ia/aad(6) and aph(3′)-IIIa genes showed no association with the resistant phenotype (susceptibility = 0%). For oxacillin and cefovecin, the blaIZR1 operon exhibited 100% susceptibility, but very low specificity (<20%), resulting in a high rate of false positives, which significantly limits its positive predictive value (PPV < 30%). No statistical significance was observed.
We utilized the mob-suite (47) to identify plasmids, followed by verification of the presence of antibiotic and biocide resistance genes in the plasmid sequences using AMRfinderPlus (46). A total of 54 plasmids were identified, with 28/54 (51.85%) found in 10/10 (100.00%) genomes of dog-associated SP and 26/54 (48.15%) in 10/11 (90.90%) genomes of guardian-associated SP (Supplementary Table S5). We identified 4/54 (7.41%) mobilizable plasmids, of which 3/4 (75.00%) were from dogs and 1/4 (25.00%) from humans, and 50/54 (92.59%) non-mobilizable plasmids, evenly distributed between dogs and humans SP (Figures 2a,b). The plasmids were present across all STs, with between 0 and 5 resistance genes identified per plasmid. One dog SP plasmid (C21B, Supplementary Table S5) carried five resistance genes (ant(6)-Ia, aph(3′)-IIIa, catA, erm(B), sat4). Two dog SP plasmids and one human SP plasmid carried 3 resistance genes (ant(6)-Ia, aph(3′)-IIIa, sat4). One dog plasmid harbored the erm(B) gene, which was absent in human plasmids. The qacG gene was found in 9 human and 6 dog SP plasmids.
Figure 2. Mobility and antimicrobial resistance genes of S. pseudintermedius plasmids from dogs and their guardians. (a) Sankey diagram showing distribution of plasmid mobility and AMR genes among plasmids detected in S. pseudintermedius isolates from dogs. (b) Sankey diagram showing plasmid mobility and AMR gene profiles among plasmids identified in isolates from dog guardians. The width of the flows is proportional to the number of plasmids in each category.
As for virulence, 4/21 (19.05%) genomes had the sea and seh genes, all of them isolated from humans. Two (9.52%) had sec and sell genes, isolated from a dog and his owner. The pvl and fnbB genes were not found by VFDB, which coincided with the PCR results.
Using the ANI values and core SNP distance, we constructed matrices to show the similarity between the samples (Figures 1b,c). Genome-wide ANI values for every possible pair of SP genomes ranged from 99.04 to 100.00% (Figure 1b; Supplementary Table S4). Core SNP distance for every possible pair of genomes ranged from 0 to 8,276 (Figure 1c). We observed greater similarity between samples of the dog and its owner in three cases (T05 and C05, T027 and C027, and T028 and C028). For the other genomes, we observed high similarity between samples isolated from the same owner (T022A and T022B, T020A, T020B and T020C) or from the same dog (C020A and C020C, C022A and C022B). We also observed a high similarity between the guardian genomes T019B T020, which were clustered together and showed the same resistance and virulence profile.
4 Discussion
The prevalence of SP in dog samples varies; however, most studies show percentages between 80 and 91% positivity for this pathogen (26, 55, 56). Considering that this study evaluated animals with skin or ear infections, the high prevalence found (76.00%) was expected. Interestingly, a high prevalence of SP in humans was found as well. Although cases of infections in humans have been reported (8, 16, 20, 57), the carrier condition attributed to humans is considered occasional. Like this study, in most reported cases, there was previous contact with dogs (8, 56, 58). The WGS analysis showed that in three cases, isolates from humans and their dogs were highly similar, suggesting transmission between hosts. Regarding the other sequenced samples, there was no evidence of transmission, and it was hypothesized that this could have occurred due to the acquisition of other strains shed by the animals, as dogs can eliminate SP through different routes such as the nose, pharynx, urogenital tract, rectum, and even via ocular discharge (58–61), or due the acquisition of a strain present in the environment, as the pathogen can persist in the environment for a long time (62).
The new MLST allelic profiles identified in the study show the occurrence of strains different from those deposited in pubMLST. The only known ST (1005) was associated with an SP isolated from a dog with otitis, originating from the state of Mato Grosso.5 We did not detect globally recognized MRSP STs such as ST71, which has been previously detected in Brazil (23, 60). Overall, the study indicates a broad genetic diversity of SP in the region.
The methicillin resistance (oxacillin) in vitro is a strong indicator of the presence of the mecA gene in SP, located on the mobile genetic element SCCmec (63). MRSP is recognized as a complicating bacterium in veterinary hospitals (64). The prevalence of MRSP in dogs with clinical infections varies widely, with reports ranging from 2 to 94.28% (14, 26, 65–67), and in the study, MRSP was identified in less than 50% of samples from dogs and their tutors. In vitro resistance to cefovecin, a third generation long-acting β-lactam antibiotic developed for use in dogs and cats (68), has also been correlated with the presence of mecA (69). Literature reports indicate good susceptibility of SP strains to cefovecin (69–71), as observed in our study. We observed a moderate or insignificant association between phenotypic resistance to oxacillin and cefovecin and the presence of the mecA gene, respectively. The resistance to these drugs may have occurred due to alternative mechanisms, such as β-lactamase production, as evidenced by the high prevalence of β-lactamase genes in the sequenced samples, similar to findings in previous studies (21, 24).
Although human infection by MRSP has been correlated with the presence of dogs in the household (9), the occurrence of MRSP in both dogs and their tutors simultaneously is underexplored. We found that dogs and their guardians exhibited similar percentages of MRSP, as described by Findik et al. (56), suggesting that contact with dogs is indeed a risk factor for humans becoming carriers of this bacterium. However, regarding the presence of the mecA, mecA1, and mecA2 genes, the association between the simultaneous presence of these genes in a strain isolated from the dog and in its guardian was moderate or insignificant. We believe several factors may explain these results: (i) genetic exchange between strains did not actively occur; (ii) loss of resistance genes by bacteria during transmission from dog to tutor or vice versa; (iii) genotypic differences in persistently colonized SP isolates from dogs and tutors reflected in antimicrobial susceptibility (71). Adiguzel et al. (72), when analyzing 3,482 clinical isolates of methicillin-resistant Staphylococcus (MRS) in dogs and cats between 2012 and 2019, revealed a high prevalence of S. pseudintermedius and coagulase-negative Staphylococcus with frequent cross-resistance to antibiotics.
The use of antimicrobials is considered a risk factor for increasing SP resistance (73–75), a fact also observed by us for β-lactam antimicrobials. However, the association between antimicrobial use and resistance to oxacillin or cefovecin was strong or moderate, suggesting once again that other factors may be influencing this relationship.
Although SP has shown good susceptibility to gentamicin (66, 76–78), resistant strains have been detected (79, 80). We found a low percentage of resistance in isolation from dogs and tutors, and there was no association between MIC values among SP isolated from dogs and their tutors. Of the sequenced samples, 6/21 (%) exhibited at least one aminoglycoside resistance gene, but only one sample from a tutor showed phenotypic resistance. This may be related to different selection pressures, as even though the pathogen carries resistance genes, their expression is regulated by challenges imposed on the pathogen in the environment (81). Two other human samples exhibited phenotypic resistance without the presence of genes related to this class, suggesting that other mechanism(s) may have been responsible for the resistance. Although no gene showed a significant association with phenotypic resistance, the aac(6′)-Ie/aph(2″)-Ia gene is commonly associated with gentamicin resistance (82), unlike ANT(6)-Ia/aad(6) and aph(3′)-IIIa, which are typically linked to resistance to streptomycin and kanamycin/neomycin, respectively (83, 84).
Regarding resistance to macrolides, we found a low number of strains with the msr(A) gene. This gene encodes an ABC transporter and rRNA methylase that can confer resistance to streptogramin B and macrolide antibiotics (85, 86). However, of the sequenced samples, 7/21 (33.33%) exhibited at least one gene related to resistance to this class, indicating that these strains possess mechanisms that qualify them to resist drugs in this class. Of the samples sequenced, one guardian sample was positive for the msrA gene in the PCR, and negative after sequencing. This gene has already been found exclusively in plasmids (87), and we believe that this isolate may have lost a plasmid containing this gene, since the sample was reactivated at two different times: one for PCR and the other for DNA extraction to send for sequencing. Another plausible hypothesis is that the gene was fragmented across multiple contigs during sequencing, which may have impeded its detection by AMRfinderplus due to coverage-based filtering thresholds.
In addition to genes related to resistance to aminoglycosides, β-lactams, and macrolides, we identified genes conferring resistance to phenicols, streptothricin, tetracyclines, trimethoprim, and point mutations in genes capable of conferring resistance to quinolones. This finding demonstrates a diverse genetic repertoire in SP in Brazil, consistent with previous observations (24). The near-universal presence of the blaIZR1 operon across the sequenced genomes has been reported in previous studies (12, 13) and suggests a potential for β-lactam resistance. However, phenotypic resistance was observed in only six isolates, and the lack of a strong genotype–phenotype association may be attributed to (1) differences in gene expression (88) and (2) polymorphisms within the blaZ gene, which may have deleterious effects on its function (89).
Additionally, a gene related to resistance to heavy metals (qacG) was found in genomes of SP isolated from dogs and tutors, suggesting that these strains have the capacity to resist disinfectants such as chlorhexidine digluconate (90), potentially becoming a persistent problem in veterinary environments or even within the household environment.
The process of conjugation between plasmids is a crucial mechanism of horizontal gene transfer (HGT), particularly significant for the spread of resistance genes (91). We have identified resistance genes associated with aminoglycosides, macrolides, phenicols, streptothricin and quarternary ammonium, especially in non-mobilizable plasmids, i.e., plasmids incapable of initiating conjugation themselves (47) that relies on transformation or transduction for propagation (92, 93). Overall, SP from both humans and dogs demonstrated the ability to exchange resistance genes to important antibiotics and biocides via HGT.
The absence of the PVL gene and the low prevalence of fnbB gene were unexpected, as these genes are important for establishing cytotoxic effects (2) and promoting adhesion to host cells (94, 95). Conversely, in the sequenced genomes, we detected the presence of virulent genes sea, sec, seh, and sell, related to five individuals and one animal sample. These toxins exhibit emetic and superantigenic activities in primate models, categorizing them as SEs (superantigens that cause food poisoning), with sea, sec, and seh genes indicating potential to cause foodborne illness (96). The sea and sec genes have been previously described in human samples, suggesting higher carriage rates in human isolates compared to animal and environmental isolates (97). These genes were also found in samples from both healthy and diseased animals, implying that these enterotoxins may not significantly contribute to disease pathogenesis in animals (98, 99).
Human contact with SP-positive dogs increases the likelihood of humans becoming a carrier of the bacteria by fivefold. However, sleeping with the dog increased the likelihood of the dog becoming SP-positive by sixfold compared to those that slept in their own bed. Our findings suggest that, depending on the tutor’s immunity and the colonization abilities of specific bacterial strains, the tutor may serve as a reservoir for the bacterium, potentially facilitating its transmission to the dog. Hanselman et al. (55) found that 50% of veterinarians were nasal carriers of SP, indicating that repeated contact with dogs may lead to carrier status. Additionally, dogs bathed in grooming facilities were four times more likely to be colonized by SP than those bathed at home, and we hypothesize that exposure to other dogs, the environment, and grooming professionals may promote SP colonization in these animals.
The diagnostic test provides values for susceptibility, specificity, positive predictive value, and negative predictive value, which are indicators of the viability of an analytical method. Comparing the results of disk diffusion tests, commonly used in most commercial laboratories, with those of MIC, considered the gold standard, can assist veterinarians, especially those dealing with small animals, in clinical practice. The susceptibility and specificity results showed that the disk diffusion method is not a good option for determining strains resistant to these antimicrobials. This method, while having good specificity, may result in a higher number of false negatives and classify resistant strains as susceptibility (48). In assessing antimicrobial resistance, it is ideal for tests to have high susceptibility so that the number of strains falsely classified as susceptibility is also low. Using an antimicrobial to which the bacterium is resistant will increase selective pressure for these resistant strains, leading to treatment failure or recurrences (74, 75).
5 Conclusion
There was a high prevalence of SP isolates in samples from dogs, tutors, and both concurrently, with contact with positive dogs being a risk factor for human colonization. The main risk factors for SP infection in dogs were sleeping with a positive tutor and frequenting grooming facilities. The isolates showed low rates of resistance to oxacillin and cefovecin, and the presence of the mecA gene was not the most important mechanism. Previous use of β-lactams moderately influenced strain resistance. The disk diffusion diagnostic method was not effective for identifying strains resistant to oxacillin and cefovecin.
Given the genetic variability, diversity of resistance genes and plasmids, and evidence of SP transmission between dogs and humans reported in this study, the results demonstrated that SP could pose a risk to humans, like other Staphylococcus species. Continuous surveillance is necessary to understand the bacterial population in animals and to minimize transmission and disease in both animals and humans.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary material.
Ethics statement
The studies involving humans were approved by Human Research Ethics Committee (CEP) under protocol number 2,980,894. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study. The animal studies were approved by Local Animal Ethics Committee (CEUA) under protocol number 055/18. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent was obtained from the owners for the participation of their animals in this study.
Author contributions
CM: Conceptualization, Investigation, Writing – original draft, Writing – review & editing. RM: Investigation, Writing – review & editing. CB: Investigation, Writing – review & editing. RR: Investigation, Writing – review & editing. AP: Investigation, Writing – review & editing. RB: Investigation, Writing – review & editing. GL: Investigation, Writing – review & editing. BF: Formal analysis, Methodology, Writing – review & editing. AB: Formal analysis, Visualization, Writing – review & editing. MH: Investigation, Writing – review & editing. MV: Formal analysis, Writing – review & editing. FA: Investigation, Writing – review & editing. VA: Investigation, Writing – review & editing. BB: Investigation, Resources, Writing – review & editing. MC: Investigation, Writing – review & editing. DR: Project administration, Resources, Writing – review & editing.
Funding
The author(s) declare that financial support was received for the research and/or publication of this article. This work was supported by the Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) and Fundação Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES 001).
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.
Generative AI statement
The authors declare that no Gen AI was used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
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.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2025.1652593/full#supplementary-material
Footnotes
1. ^https://github.com/tseemann/snp-dists
2. ^https://github.com/tseemann/mlst
3. ^https://pubmlst.org/bigsdb?db=pubmlst_spseudintermedius_seqdef
4. ^https://github.com/tseemann/abricate
5. ^https://pubmlst.org/bigsdb?page=info%26db=pubmlst_spseudintermedius_isolates%26id=1125
References
1. Devriese, LA, Vancanneyt, M, Baele, M, Vaneechoutte, M, De Graef, E, Snauwaert, C, et al. Staphylococcus pseudintermedius sp. nov., a coagulase-positive species from animals. Int J Syst Evol Microbiol Microbiol Soc. (2005) 55:1569–73. doi: 10.1099/ijs.0.63413-0
2. Futagawa-Saito, K, Makino, S, Sunaga, F, Kato, Y, Sakurai-Komada, N, Ba-Thein, W, et al. Identification of first exfoliative toxin in Staphylococcus pseudintermedius. FEMS Microbiol Lett. (2009) 301:176–80. doi: 10.1111/j.1574-6968.2009.01823.x
3. Wladyka, B, Bista, M, Sabat, AJ, Bonar, E, Grzeszczuk, S, Hryniewicz, W, et al. A novel member of the thermolysin family, cloning and biochemical characterization of metalloprotease from Staphylococcus pseudintermedius. Acta Biochim Pol. (2008) 55:525–36. doi: 10.18388/abp.2008_3059
4. Breathnach, RM, Baker, KP, Quinn, PJ, Mcgeady, TA, Aherne, CM, and Jones, BR. Clinical, immunological and histopathological findings in a subpopulation of dogs with pododermatitis. Vet Dermatol. (2005) 16:364–72. doi: 10.1111/j.1365-3164.2005.00471.x
5. Ruscher, C, Lübke-Becker, A, Wleklinski, CG, Şoba, A, Wieler, LH, and Walther, B. Prevalence of methicillin-resistant Staphylococcus pseudintermedius isolated from clinical samples of companion animals and equidaes. Vet Microbiol. (2009) 136:197–201. doi: 10.1016/j.vetmic.2008.10.023
6. Frank, LA, and Loeffler, A. Meticillin-resistant Staphylococcus pseudintermedius: clinical challenge and treatment options. Vet Dermatol. (2012) 23:283–e56. doi: 10.1111/j.1365-3164.2012.01047.x
7. Ross Fitzgerald, J. The Staphylococcus intermedius group of bacterial pathogens: species re-classification, pathogenesis and the emergence of meticillin resistance. Vet Dermatol. (2009) 20:490–5. doi: 10.1111/j.1365-3164.2009.00828.x
8. Kadlec, K, Schwarz, S, Perreten, V, Grönlund Andersson, U, Finn, M, Greko, C, et al. Molecular analysis of methicillin-resistant Staphylococcus pseudintermedius of feline origin from different European countries and North America. J Antimicrob Chemother. (2010) 65:1826–8. doi: 10.1093/jac/dkq203
9. van Duijkeren, E, Kamphuis, M, van der Mije, IC, Laarhoven, LM, Duim, B, Wagenaar, JA, et al. Transmission of methicillin-resistant Staphylococcus pseudintermedius between infected dogs and cats and contact pets, humans and the environment in households and veterinary clinics. Vet Microbiol. (2011) 150:338–43. doi: 10.1016/j.vetmic.2011.02.012
10. Ojo, KK, Striplin, MJ, Ulep, CC, Close, NS, Zittle, J, Luis, H, et al. Staphylococcus efflux msr(a) gene characterized in Streptococcus, Enterococcus, Corynebacterium, and Pseudomonas isolates. Antimicrob Agents Chemother. (2006) 50:1089–91. doi: 10.1128/AAC.50.3.1089-1091.2006
11. Kadlec, K, Feßler, AT, Hauschild, T, and Schwarz, S. Novel and uncommon antimicrobial resistance genes in livestock-associated methicillin-resistant Staphylococcus aureus. Clin Microbiol Infect. (2012) 18:745–55. doi: 10.1111/j.1469-0691.2012.03842.x
12. Saengsawang, P, Tanonkaew, R, Kimseng, R, Nissapatorn, V, Wintachai, P, Rodríguez-Ortega, MJ, et al. Whole genome sequence analysis of multidrug-resistant Staphylococcus aureus and Staphylococcus pseudintermedius isolated from superficial pyoderma in dogs and cats. Antibiotics. (2025) 14:643. doi: 10.3390/antibiotics14070643
13. Miszczak, M, Korzeniowska-Kowal, A, Wzorek, A, Prorok, P, Szenborn, L, Rypuła, K, et al. Staphylococcus aureus and Staphylococcus pseudintermedius isolated from humans and pets - a comparison of drug resistance and risk factors associated with colonisation. J Vet Res. (2025) 69:199–211. doi: 10.2478/jvetres-2025-0036
14. Somayaji, R, Priyantha, MAR, Rubin, JE, and Church, D. Human infections due to Staphylococcus pseudintermedius, an emerging zoonosis of canine origin: report of 24 cases. Diagn Microbiol Infect Dis. (2016) 85:471–6. doi: 10.1016/j.diagmicrobio.2016.05.008
15. Talan, DA, Goldstein, E, Staatz, D, and Overturf, GD. Staphylococcus intermedius: clinical presentation of a new human dog bite pathogen. Ann Emerg Med. (1989) 18:410–3. doi: 10.1016/S0196-0644(89)80582-7
16. Barnham, M, and Holmes, B. Isolation of CDC group M-5 and Staphylococcus intermedius from infected dog bites. J Infect. (1992) 25:332–4. doi: 10.1016/0163-4453(92)91759-5
17. Kempker, R, Eaton, M, Mangalat, D, and Kongphet-Tran, T. Beware of the pet dog: a case of Staphylococcus intermedius infection. Am J Med Sci. (2009) 338:425–7. doi: 10.1097/MAJ.0b013e3181b0baa9
18. Guardabassi, L, Loeber, ME, and Jacobson, A. Transmission of multiple antimicrobial-resistant Staphylococcus intermedius between dogs affected by deep pyoderma and their owners. Vet Microbiol. (2004) 98:23–7. doi: 10.1016/j.vetmic.2003.09.021
19. Vincze, S, Paasch, A, Walther, B, Ruscher, C, and Barbara, K. Multidrug- and methicillin resistant Staphylococcus pseudintermedius as a cause of canine pyoderma: a case report. Berl Munch Tierarztl Wochenschr. (2010) 123:353–8.
20. Soedarmanto, I, Kanbar, T, Ülbegi-Mohyla, H, Hijazin, M, Alber, J, Lämmler, C, et al. Genetic relatedness of methicillin-resistant Staphylococcus pseudintermedius (MRSP) isolated from a dog and the dog owner. Res Vet Sci. (2011) 91:e25–7. doi: 10.1016/j.rvsc.2011.01.027
21. Paul, NC, Moodley, A, Ghibaudo, G, and Guardabassi, L. Carriage of methicillin-resistant Staphylococcus pseudintermedius in small animal veterinarians: indirect evidence of zoonotic transmission. Zoonoses Public Health. (2011) 58:533–9. doi: 10.1111/j.1863-2378.2011.01398.x
22. Bruce, SA, Smith, JT, Mydosh, JL, Ball, J, Needle, DB, Gibson, R, et al. Accessory genome dynamics of local and global Staphylococcus pseudintermedius populations. Front Microbiol. (2022) 13:798175. doi: 10.3389/fmicb.2022.798175
23. Srednik, ME, Perea, CA, Giacoboni, GI, Hicks, JA, and Schlater, LK. First report of Staphylococcus pseudintermedius ST71-SCCmec III and ST45-ΨSCCmec57395 from canine pyoderma in Argentina. BMC Res Notes. (2023) 16:19. doi: 10.1186/s13104-023-06285-3
24. Penna, B, Silva, MB, Botelho, AMN, Ferreira, FA, Ramundo, MS, Silva-Carvalho, MC, et al. Detection of the international lineage ST71 of methicillin-resistant Staphylococcus pseudintermedius in two cities in Rio de Janeiro state. Braz J Microbiol. (2022) 53:2335–41. doi: 10.1007/s42770-022-00852-9
25. De Souza, TGV, Santana, JA, MMS, C, Pereira, ST, RGC, X, do Amarante, VS, et al. Occurrence, genetic diversity, and antimicrobial resistance of methicillin-resistant Staphylococcus spp. in hospitalized and non-hospitalized cats in Brazil. PLoS One. (2024) 19:e0309711. doi: 10.1371/journal.pone.0309711
26. Takashi, S, Sae, T, Yoshikazu, T, Arihito, S, Masayuki, O, Shintaro, H, et al. Multiplex-PCR method for species identification of coagulase-positive staphylococci. J Clin Microbiol. (2010) 48:765–9. doi: 10.1128/jcm.01232-09
27. Clinical and Laboratory Standards Institute (CLSI). Performance standards for antimicrobial disk and dilution susceptibility tests for Bacteria isolated from animals. 7th ed. Wayne, PA: (2024).
28. Clinical and Laboratory Standards Institute (CLSI). Performance standards for antimicrobial susceptibility testing. 34th ed. Wayne, PA: (2024).
29. Andrews, S. FastQC: a quality control tool for high throughput sequence data. Interne; (2010). Available online at: https://www.bioinformatics.babraham.ac.uk/projects/fastqc/ (Accessed March 15, 2025).
30. Chen, S, Zhou, Y, Chen, Y, and Gu, J. Fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics. (2018) 34:i884–90. doi: 10.1093/bioinformatics/bty560
31. Bankevich, A, Nurk, S, Antipov, D, Gurevich, AA, Dvorkin, M, Kulikov, AS, et al. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol. (2012) 19:455–77. doi: 10.1089/cmb.2012.0021
32. Shen, W, Le, S, Li, Y, and Hu, F. SeqKit: a cross-platform and ultrafast toolkit for FASTA/Q file manipulation. PLoS One. (2016) 11:e0163962. doi: 10.1371/journal.pone.0163962
33. Parks, DH, Imelfort, M, Skennerton, CT, Hugenholtz, P, and Tyson, GW. CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res. (2015) 25:1043–55. doi: 10.1101/gr.186072.114
34. Gurevich, A, Saveliev, V, Vyahhi, N, and Tesler, G. QUAST: quality assessment tool for genome assemblies. Bioinformatics. (2013) 29:1072–5. doi: 10.1093/bioinformatics/btt086
35. Olson, RD, Assaf, R, Brettin, T, Conrad, N, Cucinell, C, Davis, JJ, et al. Introducing the bacterial and viral bioinformatics resource center (BV-BRC): a resource combining PATRIC, IRD and ViPR. Nucleic Acids Res. (2023) 51:D678–89. doi: 10.1093/nar/gkac1003
36. Chaumeil, PA, Mussig, AJ, Hugenholtz, P, and Parks, DH. GTDB-Tk v2: memory friendly classification with the genome taxonomy database. Bioinformatics. (2022) 38:5315–6. doi: 10.1093/bioinformatics/btac672
37. Jain, C, Rodriguez-R, LM, Phillippy, AM, Konstantinidis, KT, and Aluru, S. High throughput ANI analysis of 90K prokaryotic genomes reveals clear species boundaries. Nat Commun. (2018) 9:5114. doi: 10.1038/s41467-018-07641-9
38. Seemann, T. Prokka: rapid prokaryotic genome annotation. Bioinformatics. (2014) 30:2068–9. doi: 10.1093/bioinformatics/btu153
39. Tonkin-Hill, G, MacAlasdair, N, Ruis, C, Weimann, A, Horesh, G, Lees, JA, et al. Producing polished prokaryotic pangenomes with the Panaroo pipeline. Genome Biol. (2020) 21:180. doi: 10.1186/s13059-020-02090-4
40. Löytynoja, A. Phylogeny-aware alignment with PRANK In: DJ Russell, editor. Multiple sequence alignment methods. Totowa, NJ: Humana Press (2014). 155–70.
41. Page, AJ, Taylor, B, Delaney, AJ, Soares, J, Seemann, T, Keane, JA, et al. Snp-sites: rapid efficient extraction of SNPs from multi-FASTA alignments. Microb Genom Microbiol Soc. (2016) 2:e000056. doi: 10.1099/mgen.0.000056
42. Stamatakis, A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. (2014) 30:1312–3. doi: 10.1093/bioinformatics/btu033
43. Tavaré, S, and Miura, RM. Some probabilistic and statistical problems in the analysis of DNA sequences. Providence, RI: American Mathematical Society (1986).
44. Letunic, I, and Bork, P. Interactive tree of life (iTOL) v4: recent updates and new developments. Nucleic Acids Res. (2019) 47:W256–9. doi: 10.1093/nar/gkz239
45. Liu, B, Zheng, D, Jin, Q, Chen, L, and Yang, J. VFDB 2019: a comparative pathogenomic platform with an interactive web interface. Nucleic Acids Res. (2019) 47:D687–92. doi: 10.1093/nar/gky1080
46. Feldgarden, M, Brover, V, Gonzalez-Escalona, N, Frye, JG, Haendiges, J, Haft, DH, et al. AMRFinderPlus and the reference gene catalog facilitate examination of the genomic links among antimicrobial resistance, stress response, and virulence. Sci Rep. (2021) 11:12728. doi: 10.1038/s41598-021-91456-0
47. Robertson, J, and Nash, JHE. MOB-suite: software tools for clustering, reconstruction and typing of plasmids from draft assemblies. Microb Genom. Microbiol. Soc. (2018) 4:e000206. doi: 10.1099/mgen.0.000206
49. Sapra, RL, and Saluja, S. Understanding statistical association and correlation. Curr Med Res Pract. (2021) 11:31–8. doi: 10.4103/cmrp.cmrp_62_20
50. Lian Duan, W, Street, N, Liu, Y, Songhua, X, and Brook, W. Selecting the right correlation measure for binary data. ACM Trans Knowl Discov. (2014) 9:28. doi: 10.1145/2637484
52. Wickham, H, Averick, M, Bryan, J, Chang, W, LD, MG, François, R, et al. Welcome to the Tidyverse. J Open Source Soft. (2019) 43:1686. doi: 10.21105/joss.01686
53. Stevenson, M, and Sergeant, E. epiR: tools for the analysis of epidemiological data. R package version 2.0.85 [internet]. (2025). Available online at: https://CRAN.R-project.org/package=epiR (Accessed March 15, 2025).
54. Kuhn, M. Building predictive models in R using the caret package. J Stat Softw. (2008) 28:1–26. doi: 10.18637/jss.v028.i05
55. Hanselman, BA, Kruth, SA, Rousseau, J, and Weese, JS. Coagulase positive staphylococcal colonization of humans and their household pets. Can Vet J. (2009) 50:954–8.
56. Findik, A, Çiftci, A, Önyay, T, Sezener, MG, Koçak, Y, and Gülhan, T. Determination of methicillin resistance and some genotypic characteristics of staphylococci isolated from dogs and their owners. Turk J Vet Anim Sci. (2018) 42:549–55. doi: 10.3906/vet-1611-50
57. Stegmann, R, Burnens, A, Maranta, CA, and Perreten, V. Human infection associated with methicillin-resistant Staphylococcus pseudintermedius ST71. J Antimicrob Chemother. (2010) 65:2047–8. doi: 10.1093/jac/dkq241
58. Frank, LA, Kania, SA, Kirzeder, EM, Eberlein, LC, and Bemis, DA. Risk of colonization or gene transfer to owners of dogs with meticillin-resistant Staphylococcus pseudintermedius. Vet Dermatol. (2009) 20:496–501. doi: 10.1111/j.1365-3164.2009.00826.x
59. Rubin, JE, and Chirino-Trejo, M. Prevalence, sites of colonization, and antimicrobial resistance among Staphylococcus pseudintermedius isolated from healthy dogs in Saskatoon, Canada. J Vet Diagn Invest. (2011) 23:351–4. doi: 10.1177/104063871102300227
60. Viegas, FM, Santana, JA, Silva, BA, Xavier, RGC, Bonisson, CT, Câmara, JLS, et al. Occurrence and characterization of methicillin-resistant Staphylococcus spp. in diseased dogs in Brazil. PLoS One. (2022) 17:e0269422. doi: 10.1371/journal.pone.0269422
61. Cengiz, S, Okur, S, Oz, C, Turgut, F, Gumurcinler, B, Sevuk, NS, et al. Prevalence and clonal diversity of methicillin-resistant Staphylococcus aureus and methicillin-resistant Staphylococcus pseudintermedius isolated from dogs and cats with eye discharge. Acta Microbiol Immunol Hung. (2023) 70:134–41. doi: 10.1556/030.2023.01899
62. Røken, M, Iakhno, S, Haaland, AH, Wasteson, Y, and Bjelland, AM. Transmission of methicillin-resistant Staphylococcus spp. from infected dogs to the home environment and owners. Antibiotics. (2022) 11:637. doi: 10.3390/antibiotics11050637
63. Kunyan, Z, Jo-Ann, MC, Sameer, E, Thomas, L, and Conly, JM. 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. (2005) 43:5026–33. doi: 10.1128/jcm.43.10.5026-5033.2005
64. Takashi, S, Ken, K, Yoshikazu, T, Namiko, T, Shinichi, K, and Keiichi, H. Methicillin-resistant Staphylococcus pseudintermedius in a veterinary teaching hospital. J Clin Microbiol. (2007) 45:1118–25. doi: 10.1128/jcm.02193-06
65. Perreten, V, Kadlec, K, Schwarz, S, Grönlund 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. doi: 10.1093/jac/dkq078
66. Scherer, CB, Botoni, LS, Coura, FM, Silva, RO, dos Santos, RD, Heinemann, MB, et al. Frequency and antimicrobial susceptibility of Staphylococcus pseudintermedius in dogs with otitis externa. Cienc Rural. (2018) 48:1–7. doi: 10.1590/0103-8478cr20170738
67. Bourguignon, E, Viçosa, GN, Corsini, CMM, Moreira, MAS, and Nero, LA. Description of methicillin-resistant Staphylococcus pseudintermedius from canine pyoderma in Minas Gerais state, Brazil. Arq Bras Med Vet Zootec. (2016) 68:299–306. doi: 10.1590/1678-4162-8150
68. EMEA (European Medicine Agency). Convenia: Scientific discussion. (2006). Available online at: https://www.ema.europa.eu (Accessed February 9, 2025).
69. Iyori, K, Toyoda, Y, Ide, K, Iwasaki, T, and Nishifuji, K. Usefulness of cefovecin disk-diffusion test for predicting mecA gene-containing strains of Staphylococcus pseudintermedius and clinical efficacy of cefovecin in dogs with superficial pyoderma. Vet Dermatol. (2013) 24:162–e36. doi: 10.1111/j.1365-3164.2012.01100.x
70. Hariharan, H, Gibson, K, Peterson, R, Frankie, M, Matthew, V, Daniels, J, et al. Staphylococcus pseudintermedius and Staphylococcus schleiferi subspecies coagulans from canine pyoderma cases in Grenada, West Indies, and their susceptibility to Beta-lactam drugs. Vet Med Int. (2014) 2014:850126. doi: 10.1155/2014/850126
71. Hartmann, FA, White, DG, West, SEH, Walker, RD, and DeBoer, DJ. Molecular characterization of Staphylococcus intermedius carriage by healthy dogs and comparison of antimicrobial susceptibility patterns to isolates from dogs with pyoderma. Vet Microbiol. (2005) 108:119–31. doi: 10.1016/j.vetmic.2005.03.006
72. Adiguzel, MC, Schaefer, K, Rodriguez, T, Ortiz, J, and Sahin, O. Prevalence, mechanism, genetic diversity, and cross-resistance patterns of methicillin-resistant Staphylococcus isolated from companion animal clinical samples submitted to a veterinary diagnostic laboratory in the Midwestern United States. Antibiotics. (2022) 11:609. doi: 10.3390/antibiotics11050609
73. Schwarz, S, and Noble, WC. Aspects of bacterial resistance to antimicrobials used in veterinary dermatological practice. Vet Dermatol. (1999) 10:163–76. doi: 10.1046/j.1365-3164.1999.00170.x
74. Schwarz, S, and Chaslus-Dancla,. Use of antimicrobials in veterinary medicine and mechanisms of resistance. Vet Res. (2001) 32:201–25. doi: 10.1051/vetres:2001120
75. Holzbauer, S, and Chiller, T. Antimicrobial resistance in bacteria of animal origin. Emerg Infect Dis. (2006) 12:1180–12. doi: 10.3201/eid1207.060503
76. Osland, AM, Vestby, LK, Fanuelsen, H, Slettemeås, JS, and Sunde, M. Clonal diversity and biofilm-forming ability of methicillin-resistant Staphylococcus pseudintermedius. J Antimicrob Chemother. (2012) 67:841–8. doi: 10.1093/jac/dkr576
77. Metiner, K, Bagcigil, AF, and Ilgaz, A. Determination of the diversity and antibiotic resistance profiles of Staphylococcus species from dogs with otitis externa and examination of mecA gene occurrence. Vet Med (Praha). (2015) 60:261–7. doi: 10.17221/8178-VETMED
78. Rubin, JE, and Chirino-Trejo, M. Antimicrobial susceptibility of canine and human Staphylococcus aureus collected in Saskatoon, Canada. Zoonoses Public Health. (2011) 58:454–62. doi: 10.1111/j.1863-2378.2011.01392.x
79. Nienhoff, U, Kadlec, K, Chaberny, IF, Verspohl, J, Gerlach, GF, Kreienbrock, L, et al. Methicillin-resistant Staphylococcus pseudintermedius among dogs admitted to a small animal hospital. Vet Microbiol. (2011) 150:191–7. doi: 10.1016/j.vetmic.2010.12.018
80. Casagrande Proietti, P, Bietta, A, Coletti, M, Marenzoni, ML, Scorza, AV, and Passamonti, F. Insertion sequence IS256 in canine pyoderma isolates of Staphylococcus pseudintermedius associated with antibiotic resistance. Vet Microbiol. (2012) 157:376–82. doi: 10.1016/j.vetmic.2011.12.028
81. Hipólito, A, García-Pastor, L, Blanco, P, Trigo da Roza, F, Kieffer, N, Vergara, E, et al. The expression of aminoglycoside resistance genes in integron cassettes is not controlled by riboswitches. Nucleic Acids Res. (2022) 50:8566–79. doi: 10.1093/nar/gkac662
82. Caldwell, SJ, Huang, Y, and Berghuis, AM. Antibiotic binding drives catalytic activation of aminoglycoside kinase APH(2″)-Ia. Structure. (2016) 24:935–45. doi: 10.1016/j.str.2016.04.002
83. Ramirez, MS, and Tolmasky, ME. Aminoglycoside modifying enzymes. Drug Resist Updat. (2010) 13:151–71. doi: 10.1016/j.drup.2010.08.003
84. Boehr, DD, Thompson, PR, and Wright, GD. Molecular mechanism of aminoglycoside antibiotic kinase APH(3′)-IIIa: roles of conserved active site residues*. J Biol Chem. (2001) 276:23929–36. doi: 10.1074/jbc.M100540200
85. Boerlin, P, Burnens, AP, Frey, J, Kuhnert, P, and Nicolet, J. Molecular epidemiology and genetic linkage of macrolide and aminoglycoside resistance in Staphylococcus intermedius of canine origin. Vet Microbiol. (2001) 79:155–69. doi: 10.1016/S0378-1135(00)00347-3
86. Lüthje, P, and Schwarz, S. Molecular basis of resistance to macrolides and lincosamides among staphylococci and streptococci from various animal sources collected in the resistance monitoring program BfT-GermVet. Int J Antimicrob Agents. (2007) 29:528–35. doi: 10.1016/j.ijantimicag.2006.12.016
87. Sánchez-Osuna, M, Pedrosa, M, Bierge, P, Gómez-Sánchez, I, Alguacil-Guillén, M, Espasa, M, et al. Genomic analysis of Staphylococcus aureus isolates from bacteremia reveals genetic features associated with the COVID-19 pandemic. iScience. (2024) 27:402. doi: 10.1016/j.isci.2024.110402
88. Florence, D, Isabelle, P, Roland, L, Ekkehard, C, and Patrice, C. Modes and modulations of antibiotic resistance gene expression. Clin Microbiol Rev. (2007) 20:79–114. doi: 10.1128/cmr.00015-06
89. Rocha, GD, Nogueira, JF, Gomes dos Santos, MV, Boaventura, JA, Nunes Soares, RA, de Simoni, J, et al. Impact of polymorphisms in blaZ, blaR1 and blaI genes and their relationship with β-lactam resistance in S. aureus strains isolated from bovine mastitis. Microb Pathog. (2022) 165:105453. doi: 10.1016/j.micpath.2022.105453
90. Walker, MA, Singh, A, Gibson, TW, Rousseau, J, and Weese, JS. Presence of Qac genes in clinical isolates of methicillin-resistant and methicillin-susceptible Staphylococcus pseudintermedius and their impact on chlorhexidine digluconate susceptibility. Vet Surg. (2020) 49:971–6. doi: 10.1111/vsu.13413
91. Ramsay, JP, and Firth, N. Diverse mobilization strategies facilitate transfer of non-conjugative mobile genetic elements. Mob Genet Elem HGT Prokaryotes Microbiota. (2017) 38:1–9. doi: 10.1016/j.mib.2017.03.003
92. Ankita, K, Yu-Wei, W, John-Marc, C, Marimikel, C, Lara, R, Rocha Andrea, M, et al. Large circular plasmids from groundwater Plasmidomes span multiple incompatibility groups and are enriched in multimetal resistance genes. mBio. (2019) 10:18. doi: 10.1128/mbio.02899-18
93. Bardiau, M, Yamazaki, K, Ote, I, Misawa, N, and Mainil, JG. Characterization of methicillin-resistant Staphylococcus pseudintermedius isolated from dogs and cats. Microbiol Immunol. (2013) 57:496–501. doi: 10.1111/1348-0421.12059
94. Heilmann, C, Schweitzer, O, Gerke, C, Vanittanakom, N, Mack, D, and Götz, F. Molecular basis of intercellular adhesion in the biofilm-forming Staphylococcus epidermidis. Mol Microbiol. (1996) 20:1083–91. doi: 10.1111/j.1365-2958.1996.tb02548.x
95. Arciola, CR, Campoccia, D, Ravaioli, S, and Montanaro, L. Polysaccharide intercellular adhesin in biofilm: structural and regulatory aspects. Front Cell Infect Microbiol. (2015) 5:7. doi: 10.3389/fcimb.2015.00007
96. Lina, G, Bohach, GA, Nair, SP, Hiramatsu, K, Jouvin-Marche, E, and Mariuzza, R. Standard nomenclature for the superantigens expressed by Staphylococcus. J Infect Dis. (2004) 189:2334–6. doi: 10.1086/420852
97. Phumthanakorn, N, Fungwithaya, P, Chanchaithong, P, and Prapasarakul, N. Enterotoxin gene profile of methicillin-resistant Staphylococcus pseudintermedius isolates from dogs, humans and the environment. J Med Microbiol Microbiol Soc. (2018) 67:866–73. doi: 10.1099/jmm.0.000748
98. Tanabe, T, Toyoguchi, M, Hirano, F, Chiba, M, Onuma, K, and Sato, H. Prevalence of staphylococcal enterotoxins in Staphylococcus pseudintermedius isolates from dogs with pyoderma and healthy dogs. Microbiol Immunol. (2013) 57:651–4. doi: 10.1111/1348-0421.12069
Keywords: mecA, msrA, fnbB, Panton-valentine leukocidin, resistance
Citation: Martins CR, de Melo RT, Bastos CM, Rosa RO, Pereira AG, Figueiredo Braz R, Loures GdP, Fonseca BB, Buiatte ABG, Heinemann MB, Viana MVC, Aburjaile FF, Azevedo VAdC, Brenig B, da Costa MM and Rossi DA (2025) Genetic connections and antimicrobial resistance in dogs and owners Staphylococcus pseudintermedius isolates. Front. Vet. Sci. 12:1652593. doi: 10.3389/fvets.2025.1652593
Edited by:
Jose L. Gonzales, Wageningen University and Research, NetherlandsCopyright © 2025 Martins, de Melo, Bastos, Rosa, Pereira, Figueiredo Braz, Loures, Fonseca, Buiatte, Heinemann, Viana, Aburjaile, Azevedo, Brenig, da Costa and Rossi. 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) and the copyright owner(s) 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: Roberta Torres de Melo, cm9iZXJ0YS1tZWxvQGhvdG1haWwuY29t
Christina Resende Martins1