# Evolution of the Population Structure of Staphylococcus pseudintermedius in France

Marine Bergot, Patricia Martins-Simoes, Hélène Kilian, Pierre Châtre, Kate A Worthing, Jacqueline M Norris, Jean-Yves Madec, Frédéric Laurent, Marisa Haenni

► To cite this version:

Marine Bergot, Patricia Martins-Simoes, Hélène Kilian, Pierre Châtre, Kate A Worthing, et al. Evolution of the Population Structure of Staphylococcus pseudintermedius in France. Frontiers in Microbiology, 2018, 9 (3055), 10.3389/fmicb.2018.03055. anses-03993271

HAL Id: anses-03993271

https://anses.hal.science/anses-03993271

Submitted on 16 Feb 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for
updates

# Evolution of the Population Structure of *Staphylococcus pseudintermedius* in France

Marine Bergot<fnref n="1" />, Patricia Martins-Simoes<fnref n="1" />,², Hélène Kilian³, Pierre Châtre³, Kate A. Worthing⁴,⁵, Jacqueline M. Norris⁴, Jean-Yves Madec², Frédéric Laurent<fnref n="1" />,² and Marisa Haenni³*

*Staphylococcus pseudintermedius* is a colonizer as well as an important pathogen of dogs where it is responsible for skin, ear and post-operative infections. The emergence of methicillin-resistant *S. pseudintermedius* (MRSP) in the early 2000s, which were additionally resistant to most veterinary-licensed antibiotics, drew specific attention to these pathogens due to the limitations created in veterinary therapeutic options. Multiple studies showed that the sequence type (ST)71 was the most frequently identified clone in Europe. A few years ago, several publications have suggested a decline of the ST71 clone and the emergence of the ST258 lineage in Northern Europe. In this study, we show that ST71 is also decreasing over time in France and that the non-ST71 population is highly heterogeneous. Globally, the non-ST71 clones are more susceptible to antibiotics, which might be good news for veterinarians. Two other lineages, ST258 and ST496, seem to be successful in France. These isolates, as well as representatives of the ST71 clone, underwent whole-genome sequence. This study shows that the ST71 and ST496 clusters are highly homogenous while the ST258 cluster is more diverse. Each ST possesses a specific pattern of resistance and virulence genes. The reasons for the apparent and simultaneous success of the ST258 and ST496 clones remain unclear. But the emergence of the ST496 clone will require monitoring given its multi-resistant genotype and threat to canine health.

**Keywords:** dog, MRSP, MLST, ST258, ST496

**OPEN ACCESS**

**Edited by:**
Jack Wong,
The Chinese University of Hong Kong,
China

**Reviewed by:**
Jason Sahl,
Northern Arizona University,
United States
Vasvi Chaudhry,
University of Tubingen, Germany

*Correspondence:
Marisa Haenni
marisa.haenni@anses.fr

**Specialty section:**
This article was submitted to *Antimicrobials, Resistance and Chemotherapy*, a section of the journal *Frontiers in Microbiology*

**Received:** 21 August 2018
**Accepted:** 27 November 2018
**Published:** 13 December 2018

**Citation:**
Bergot M, Martins-Simoes P, Kilian H, Châtre P, Worthing KA, Norris JM, Madec J-Y, Laurent F and Haenni M (2018) Evolution of the Population Structure of *Staphylococcus pseudintermedius* in France. *Front. Microbiol.* 9:3055. doi: 10.3389/fmicb.2018.03055

## INTRODUCTION

*Staphylococcus pseudintermedius* (SP), initially described in 2005, is a normal commensal bacterium in dogs. However, it is also recognized as the most frequent opportunistic pathogen responsible for skin, ear and surgical site infections and complications in dogs. While its natural hosts are dogs, human infections typically associated with exposure to dogs have been increasingly reported in the last few years (Lozano et al., 2017; Pomba et al., 2017), highlighting the zoonotic potential of *S. pseudintermedius* (Somayaji et al., 2016).

<footnote n="1">¹ Centre International de Recherche en Infectiologie (CIRI), Equipe Pathogénie des Staphylocoques, Faculté de Médecine Laennec, Lyon, France, ² Institut des Agents Infectieux – Centre de Biologie Nord, Centre National de Référence des Staphylocoques, Hospices Civils de Lyon, Lyon, France, ³ Unité Antibiorésistance et Virulence Bactériennes, Université de Lyon – Agence Nationale de Sécurité Sanitaire (Anses), Lyon, France, ⁴ Sydney School of Veterinary Science, University of Sydney, Camden, NSW, Australia, ⁵ Department of Microbiology and Immunology, Peter Doherty Institute for Infection and Immunity, The University of Melbourne, Parkville, VIC, Australia</footnote>
Similar to S. *aureus*, methicillin-resistance in S. *pseudintermedius* has been recognized as a major issue in veterinary medicine. Since the first description of methicillin-resistance S. *pseudintermedius* (MRSP) in 1999 (Gortel et al., 1999), its global distribution has expanded and the number of clinical infections associated with MRSP in companion animals (mostly dogs) has dramatically increased. With the reported prevalence as high as 60% (Kawakami et al., 2010; Haenni et al., 2014; Worthing et al., 2018), the situation is worsened by the high frequency of multi and extensively drug resistant strains markedly restricting the therapeutic arsenal for veterinarians (Perreten et al., 2010).

Methicillin-resistance S. *pseudintermedius* has been spreading worldwide through the dissemination of specific lineages from different genetic backgrounds. The most prevalent clones are ST71 in Europe, ST68 in the United States and ST45/ST112 in Asia (Pires Dos Santos et al., 2016). However, several studies in Northern countries of Europe showed that ST71 MRSP may be challenged by other less frequent clones such as ST258 MRSP (Kjellman et al., 2015; Damborg et al., 2016; Duim et al., 2016; Gronthal et al., 2017). In 2010, we detected non-ST71 MRSP lineages in French veterinarian isolates, including two ST258 isolates (Haenni et al., 2014). This ST258 MRSP lineage was reported to be more frequently susceptible than the ST71 MRSP to enrofloxacin, gentamicin, and chloramphenicol, while also being marginally more resistant to trimethoprim-sulfonamides and tetracyclines (Pires Dos Santos et al., 2016). The same trend toward susceptibility to veterinary-licensed antibiotics has also been observed for other lineages, so that the emergence of a polyclonal non-ST71 MRSP population may be good news for the treatment of such infections (Haenni et al., 2014; Kjellman et al., 2015). Moreover, non-ST71 MRSP lineages also should be less virulent as they form less biofilm and show lower adherence to canine corneocytes (Osland et al., 2012; Latronico et al., 2014).

In France, the occurrence of MRSP in clinical samples from dogs is continuously monitored through Resapath, the French network for surveillance of antimicrobial resistance in clinical pathogens. We utilized this network to compare two collections (2012–2013 and 2015–2016) of MRSP isolates collected from canine skin and soft tissue infections. The aims of our study were firstly to assess if a shift in the population structure of MRSP, from ST71 toward non-ST71 lineages, has occurred in France, and secondly to characterize the resistome and virulome of the most prevalent clones by using whole genome sequencing.

# MATERIALS AND METHODS

## Bacterial Sampling and Identification

Methicillin-resistance S. *pseudintermedius* isolated from skin and soft tissue infections in dogs were collected from the veterinary laboratories of Resapath, the network for surveillance of antimicrobial resistance in clinical pathogens in France<fnref n="1" />. Two collections were constituted: the first one presented 95 isolates sampled in 2012–2013, and the second one

87 isolates collected in 2015–2016. Each clinical isolate corresponded to a unique animal. All isolates were sent to the Anses laboratory in Lyon, France, for further analysis. Species identification was performed by PCR-RFLP using TaqI restriction analysis on the *kat* gene (Blaiotta et al., 2010). The presence of the *mecA* and 16S rRNA genes were detected by PCR (Maes et al., 2002) and used to confirm methicillin resistance and as an internal quality control, respectively.

## Antimicrobial Susceptibility Tests

All S. *pseudintermedius* isolates were tested for their antimicrobial susceptibility by disk diffusion method according to the guidelines of the French Society for Microbiology (CA-SFM<fnref n="2" />). Antibiotics tested were cefovecin 30 µg (which was used as the phenotypic marker of methicillin resistance), penicillin 6 µg, fusidic acid 10 µg, kanamycin 30 µg, gentamicin 15 µg, tobramycin 10 µg, erythromycin 15 µg, spiramycin 100 µg, lincomycin 15 µg, tetracycline 30 µg, chloramphenicol 30 µg, florfenicol 30 µg, enrofloxacin 5 µg, vancomycin 30 µg, and teicoplanin 30 µg (Mast Diagnostics, Amiens, France). Bacteria were classified as susceptible, intermediate or resistant according to the clinical breakpoints approved by the veterinary part of the CA-SFM. *Staphylococcus aureus* ATCC 25923 was used as the quality control strain.

## Multi-Locus Sequence Typing (MLST) and Clonal Complexes (CC)

Multi-locus sequence typing was performed for all MRSP isolates according to the 7-loci scheme (Solyman et al., 2013). STs were determined using the MLST website<fnref n="3" /> and new STs were assigned by the curator Vincent Perreten (vincent.perreten@vbi.unibe.ch). In addition, for each collection (2012–2013 and 2015–2016), STs were assigned to CCs by PhyloViz v 2.0 (Nascimento et al., 2017), using goeBURST (Francisco et al., 2009) algorithm with default parameters. A clonal complex is composed by all strains sharing 5 or more identical profile alleles (double-loci variants – DVL parameter).

## Whole Genome Sequencing and Analysis

Whole genome sequencing (WGS) was performed on a subset of 23 isolates belonging to the three most prevalent clones in both surveillance studies: all ST258 (n = 6), ST496 (n = 9), and a subset of ST71 (n = 8, among the 110 ST71 identified) were whole-genome sequenced. WGS was performed on genomic DNA extracted using the Microbial DNA extraction kit (Macherey Nagel, Hoerdt, France). The genomic library was prepared using the Nextera XT DNA Library Preparation kit (Illumina Inc., Cambridge, United Kingdom), according to the manufacturer's protocol, and 2x 250-bp paired-end reads were generated with Illumina sequencing technology (San Diego, CA, United States).

<footnote n="1">¹https://www.resapath.anses.fr/</footnote>

<footnote n="2">²www.sfm-microbiologie.org</footnote>

<footnote n="3">³http://pubmlst.org/spseudintermedius</footnote>
## Assembly and Annotation

We used the Nullarbor pipeline (T. Seemann<fnref n="4" />) to perform paired-end reads trimming using Trimmomatic v 0.33 (Bolger et al., 2014), and the cleaned reads were then assembled using SPAdes v 3.5 (Bankevich et al., 2012). Quality of genome assemblies was assessed with QUAST v 4.6.3 (Mikheenko et al., 2016). Always in the Nullarbor pipeline, the genomes annotation was then performed using Prokka v 1.10 (Seemann, 2014).

## In silico MLST

All publicly available genomes (NCBI genomes database, accession date 14.02.2018) *S. pseudintermedius* (*n* = 19) were *in silico* MLST typed using the MLST pipeline<fnref n="5" />.

## Resistome and Virulome

The resistome and virulome of all sequenced ST71, ST258, and ST496 genomes were characterized using the ABRicate v 0.8.1 software (Seeman T<fnref n="6" />) against the ResFinder database as updated on February 20, 2018, and VFDB database (Chen et al., 2005; Zankari et al., 2012). Two curated in-house databases (DB) were subsequently built and the genomes were screened against these databases, using ABRicate v 0.8.1, to search for homologous genes. As resistance genes of *S. pseudintermedius* are poorly described, and because *S. aureus* share common resistance genes with *S. pseudintermedius*, we built the first database (DB_SA) using 336 markers used in the *S. aureus* DNA microarray Genotyping kit v2 (Monecke et al., 2008). Furthermore, we performed a state of the art of reported virulence-associated genes in *S. pseudintermedius* and built the second database (DB_SP) using a total of 49 genes : spsA-R, clpP, siet, speta, se-int, lukF-I, lukS-I, sec-canine, exi, expB, agrA-D, icaA-D, nanB, coa, clpX, saeR, saeS, htrA, nucC, hlb, sigB, srrA, sarA, rot, traP, expA, and epbS (Ben Zakour et al., 2011; Bannoehr et al., 2012; Garbacz et al., 2013; Pietrocola et al., 2015; Couto et al., 2016; Ruzauskas et al., 2016). All gene sequences used in DB_SP, except the *sec-canine* gene, were obtained from the complete genome of *S. pseudintermedius* strain ED99 (GenBank accession number: CP002478.1).

Because *sec-canine* gene was not available for *S. pseudintermedius*, we used the *sec-canine* gene of the closest phylogenetic species *S. intermedius* (strain 95-011195, GenBank accession number: U91526.1) and we included it in DB_SP.

Presence of both virulence and resistance-associated genes was initially detected: (i) if ≥50% of the gene length and ≥50% nucleotide identity were detected when using the DB_SA; (ii) if ≥90% of the gene length and ≥70% nucleotide identity were detected when using the DB_SP. For both databases (DB_SP and DB_SA), we used ABRicate to identify the presence of virulence and resistance-associated genes. As DB_SA was built based on *S. aureus* genes, we used a lower identity threshold (i.e., 50%) than the one used for DB_SP (i.e., 90%) to be able to detect homologous genes. Then, for the DB_SA, a second validation

of positive hits was performed manually using a threshold of at least 80% of the query sequence length. We also screened for the presence of putative CRISPR regions, using CRISPRcasfinder (Grissa et al., 2007).

## Fluoroquinolone Resistance

Mutations in topoisomerase II (*gyrA*) and IV (*grlA*) are known to confer fluoroquinolone resistance in *S. pseudintermedius* (Descloux et al., 2008). Detection of mutations Ser84Leu in *gyrA* and Ser80Ile in *grlA* sequences was performed using Snippy v 3.1⁶. The fluoroquinolone-susceptible *S. pseudintermedius* strain CCUG49543T (also known as LMG 22219T) *gyrA* and *grlA* sequences as reference (accession numbers AM262968 and AM262971, respectively) (Devriese et al., 2005). Visual confirmation of non-synonymous mutations was performed with Seaview (Gouy et al., 2010).

## SNP Based Phylogeny

SNP detection was performed with Snippy v.3⁷ to construct a SNP based phylogeny for ST71, ST258, and ST496. Strain ED99 (ST25) was used as a reference outgroup/external ST to root the tree. A raw phylogenetic tree using the approximately-maximum-likelihood (ML) model and generalized time-reversible (GTR) models of nucleotide evolution were built using FastTree2 (Price et al., 2010). Strain ED99 was used to root the phylogenetic tree. To visualize and annotate the trees we used the bioconductor package *ggtree* (Yu et al., 2017). In addition, a pairwise SNP distance matrix was built with snp-dists v 0.6⁸, as implemented in Snippy v.3.

## Statistical Analysis

Proportions of antimicrobial resistance in ST71 isolates vs. non-ST71 isolates were compared using the chi-squared test. Differences were considered significant if *P* < 0.05.

## Ethics Statement

No ethical approval was needed since this study did not involve any experimentation on animals. Only bacteria collected from clinical samples were used in the frame of this study.

## Data Availability

The data sets supporting the results of this article are available from the Sequence Read Archive (SRA) under the BioProject no. PRJNA498009.

# RESULTS

## Population Structure of French MRSP

ST71 MRSP remains the dominant methicillin-resistant clone in France (Table 1) but its prevalence over the years revealed a decreasing trend, even not statistically significant (*p* = 0.1643), from 65.3% in the 2012–2013 collection to 55.2% in the 2015–2016 collection. The non-ST71 MRSP group is largely polyclonal,

<footnote n="4">⁴https://github.com/tseemann/nullarbor</footnote>

<footnote n="5">⁵https://github.com/tseemann/mlst</footnote>

<footnote n="6">⁶https://github.com/tseemann</footnote>

<footnote n="7">⁷https://github.com/tseemann/snippy</footnote>

<footnote n="8">⁸https://github.com/tseemann/snp-dists</footnote>
<table><caption>**TABLE 1** | Distribution of the STs per collection.</caption><thead><tr><th rowspan="2">Sequence type (ST)</th><th colspan="2">No. of isolates (%)</th></tr><tr><th>2012–2013 (n = 95)</th><th>2015–2016 (n = 87)</th></tr></thead><tbody><tr><td>71</td><td>62 (65.3)</td><td>48 (55.2)</td></tr><tr><td>45</td><td>2 (2.1)</td><td>0</td></tr><tr><td>68</td><td>1 (1.1)</td><td>0</td></tr><tr><td>85</td><td>1 (1.1)</td><td>3 (3.4)</td></tr><tr><td>136</td><td>0</td><td>1 (1.1)</td></tr><tr><td>181</td><td>1 (1.1)</td><td>1 (1.1)</td></tr><tr><td>196</td><td>1 (1.1)</td><td>0</td></tr><tr><td>258</td><td>1 (1.1)</td><td>5 (5.7)</td></tr><tr><td>261</td><td>8 (8.4)</td><td>3 (3.4)</td></tr><tr><td>294</td><td>2 (2.1)</td><td>2 (2.3)</td></tr><tr><td>298</td><td>0</td><td>1 (1.1)</td></tr><tr><td>339</td><td>0</td><td>1 (1.1)</td></tr><tr><td>342</td><td>0</td><td>2 (2.3)</td></tr><tr><td>406</td><td>1 (1.1)</td><td>0</td></tr><tr><td>407</td><td>1 (1.1)</td><td>0</td></tr><tr><td>408</td><td>1 (1.1)</td><td>0</td></tr><tr><td>410</td><td>1 (1.1)</td><td>0</td></tr><tr><td>413</td><td>1 (1.1)</td><td>0</td></tr><tr><td>414</td><td>4 (4.2)</td><td>1 (1.1)</td></tr><tr><td>415</td><td>1 (1.1)</td><td>0</td></tr><tr><td>416</td><td>1 (1.1)</td><td>0</td></tr><tr><td>419</td><td>1 (1.1)</td><td>0</td></tr><tr><td>420</td><td>1 (1.1)</td><td>0</td></tr><tr><td>421</td><td>1 (1.1)</td><td>0</td></tr><tr><td>424</td><td>1 (1.1)</td><td>0</td></tr><tr><td>425</td><td>1 (1.1)</td><td>0</td></tr><tr><td>450</td><td>0</td><td>1 (1.1)</td></tr><tr><td>496</td><td>0</td><td>9 (10.3)</td></tr><tr><td>551</td><td>0</td><td>1 (1.1)</td></tr><tr><td>1026</td><td>0</td><td>1 (1.1)</td></tr><tr><td>1027</td><td>0</td><td>1 (1.1)</td></tr><tr><td>1028</td><td>0</td><td>1 (1.1)</td></tr><tr><td>1029</td><td>0</td><td>1 (1.1)</td></tr></tbody></table>

with most STs occurring only once. Moreover, seventeen STs were newly identified and assigned to the following numbers: 406, 407, 408, 410, 413, 414, 415, 416, 419, 420, 421, 424, 425, 1026, 1027, 1028, and 1029. The ST258 MRSP lineage was detected in one isolate in 2012–2013 and in five isolates in 2015–2016, which were all isolated in different counties. Conversely, ST261 was identified in eight isolates from eight different counties in 2012–2013 but only in one isolate in 2015–2016. Finally, we observed the emergence of the ST496 MRSP clone, which was absent in 2012–2013 and was found in nine isolates from four different counties in 2015–2016.

## Antibiotic Susceptibility Profiles

ST71 MRSP was highly resistant to veterinary-licensed antibiotics, with chloramphenicol, tetracyclines, and fusidic acid being the last antibiotics with resistance rates lower than 90% in 2015–2016 (Tables 2, 3). Between 2012–2013 and

2015–2016, we observed a significant increase in the resistance to tobramycin, kanamycin, and chloramphenicol. In non-ST71 MRSP, we observed a higher rate of resistance to gentamicin (97.4% in 2015–2016), macrolides-lincosamides (84.8 and 63.6%, respectively, due to the presence of the inducible MLS<sub>P</sub> phenotype). Compared to ST71 MRSP isolates, the rate of resistance to kanamycin, tobramycin, and enrofloxacin (30.3% vs. >90%) was lower in non-ST71 MRSP but higher for tetracyclines (78.8% vs. 37.1%).

## Phylogeny of French Isolates

Whole genome sequencing was performed on all ST258 and ST496, as well as on a subset of eight ST71 selected from the two collections (four isolates from each collection) and originating from eight different districts. The phylogenetic tree of the French strains and the ones available on public databases shows that the ST71 and ST496 clades are highly homogenous (Figure 1). On the contrary, the ST258 cluster is more diverse. To support these observations, we used the pairwise SNP distance matrix (Supplementary Table S1) which showed that ST71 is a highly homogenous cluster with an average of 149 SNPs, while ST258 is much more heterogeneous with an average of 1919 SNPs, about thirteen time more than ST71, and twice more than ST496.

## Distribution of Resistance Genes According to MLST Types

Multiple resistance genes were screened (reported in Figure 2), which perfectly correlated with the phenotypic characterization. Each ST presented a specific pattern. For ST496, all isolates displayed an identical pattern, with a higher number of resistance genes when compared to ST71 and ST258. Concerning the ST258 MRSP group, isolates presented a conserved pattern of resistance genes (blaZ, tet(M), aphA3, truncated sat genes) with sporadic variations, while ST71 was more heterogeneous. ST258 and ST496 displayed the tet(M) gene, whereas ST71 (when present) displayed the tet(K) gene. ST258 only rarely displayed the *aacA-aphD* gene conferring resistance to gentamicin, kanamycin, and tobramycin. All ST496 and ST71 isolates were resistant to fluoroquinolones, and they all presented the prototypical mutations Ser84Leu and Ser80Ile in the *gyrA* and *grlA* genes, respectively. On the contrary, all ST258 isolates were susceptible and all isolates presented the intact *gyrA* and *grlA* genes, except one isolate presenting the mutated *grlA* gene.

## Detection of the CRISPR Sequences

CRIPSR sequences were detected in all isolates belonging to the ST496, but were absent in ST71 and ST258 isolates. Typing of the CRISPR elements found in ST496 isolates revealed that they all harbored the same subtype of CRISPR: the subtype IIIA, which is located within the SCC*mec* element found in this lineage.

## Distribution of Virulence Genes According to MLST Types

All MLST types shared a core virulence profile, with a group of 24 genes (including the biofilm-associated *ica* operon – except for two ST496 isolates in which *icaB* was absent – and *sps* surface
<table><caption>TABLE 2 | Antimicrobial resistance in isolates collected in 2012–2013 (n = 95).</caption><thead><tr><th rowspan="2">Antibiotic</th><th rowspan="2">Breakpoints (mm)<sup>1</sup> R &lt; /S ≥</th><th colspan="2">ST71 (n = 62)</th><th colspan="2">non-ST71 (n = 33)</th><th rowspan="2">p</th></tr><tr><th>No. of strains</th><th>Percentage</th><th>No. of strains</th><th>Percentage</th></tr></thead><tbody><tr><td>Kanamycin</td><td>15/17</td><td>62</td><td>100.0</td><td>29</td><td>87.9</td><td>0.023</td></tr><tr><td>Gentamicin</td><td>20/20</td><td>57</td><td>91.9</td><td>10</td><td>30.3</td><td>&lt;2.10-6</td></tr><tr><td>Tobramycin</td><td>20/20</td><td>59</td><td>95.2</td><td>10</td><td>30.3</td><td>&lt;2.10-6</td></tr><tr><td>Chloramphenicol</td><td>19/22</td><td>11</td><td>17.7</td><td>7</td><td>21.2</td><td>0.681</td></tr><tr><td>Florfenicol</td><td>19/22</td><td>0</td><td>0.0</td><td>0</td><td>0.0</td><td>-</td></tr><tr><td>Tetracyclines</td><td>17/19</td><td>23</td><td>37.1</td><td>26</td><td>78.8</td><td>0.0001</td></tr><tr><td>Erythromycin</td><td>17/22</td><td>59</td><td>95.2</td><td>28</td><td>84.8</td><td>0.182</td></tr><tr><td>Spiramycin</td><td>20/20</td><td>59</td><td>95.2</td><td>28</td><td>84.8</td><td>0.182</td></tr><tr><td>Lincomycin</td><td>17/21</td><td>58</td><td>93.5</td><td>21</td><td>63.6</td><td>0.0002</td></tr><tr><td>Fusidic acid</td><td>24/24</td><td>2</td><td>3.2</td><td>1</td><td>3.0</td><td>0.573</td></tr><tr><td>Enrofloxacin</td><td>17/22</td><td>61</td><td>98.4</td><td>10</td><td>30.3</td><td>&lt;2.10-6</td></tr><tr><td>Teicoplanin</td><td>17/17</td><td>0</td><td>0.0</td><td>0</td><td>0.0</td><td>-</td></tr><tr><td>Vancomycin</td><td>17/17</td><td>0</td><td>0.0</td><td>0</td><td>0.0</td><td>-</td></tr></tbody></table>

<table><caption>TABLE 3 | Antimicrobial resistance in isolates collected in 2015–2016 (n = 87).</caption><thead><tr><th rowspan="2">Antibiotic</th><th rowspan="2">Breakpoints (mm)<sup>1</sup> R &lt; /S ≥</th><th colspan="2">ST71 (n = 48)</th><th colspan="2">non-ST71 (n = 39)</th><th rowspan="2">p</th></tr><tr><th>No. of strains</th><th>Percentage</th><th>No. of strains</th><th>Percentage</th></tr></thead><tbody><tr><td>Kanamycin</td><td>15/17</td><td>48</td><td>100.0</td><td>38</td><td>97.4</td><td>0.44</td></tr><tr><td>Gentamicin</td><td>20/20</td><td>39</td><td>81.3</td><td>16</td><td>41.0</td><td>0.0001</td></tr><tr><td>Tobramycin</td><td>20/20</td><td>42</td><td>87.5</td><td>17</td><td>43.6</td><td>0.00001</td></tr><tr><td>Chloramphenicol</td><td>19/22</td><td>2</td><td>4.2</td><td>17</td><td>43.6</td><td>0.00001</td></tr><tr><td>Florfenicol</td><td>19/22</td><td>0</td><td>0.0</td><td>0</td><td>0.0</td><td>-</td></tr><tr><td>Tetracyclines</td><td>17/19</td><td>16</td><td>33.3</td><td>35</td><td>89.7</td><td>&lt;0.00001</td></tr><tr><td>Erythromycin</td><td>17/22</td><td>46</td><td>95.8</td><td>33</td><td>84.6</td><td>0.15</td></tr><tr><td>Spiramycin</td><td>20/20</td><td>46</td><td>95.8</td><td>33</td><td>84.6</td><td>0.15</td></tr><tr><td>Lincomycin</td><td>17/21</td><td>46</td><td>95.8</td><td>28</td><td>71.8</td><td>0.0018</td></tr><tr><td>Fusidic acid</td><td>24/24</td><td>1</td><td>2.1</td><td>3</td><td>7.7</td><td>0.467</td></tr><tr><td>Enrofloxacin</td><td>17/22</td><td>47</td><td>97.9</td><td>19</td><td>48.7</td><td>&lt;0.00001</td></tr><tr><td>Teicoplanin</td><td>17/17</td><td>0</td><td>0.0</td><td>0</td><td>0.0</td><td>-</td></tr><tr><td>Vancomycin</td><td>17/17</td><td>0</td><td>0.0</td><td>0</td><td>0.0</td><td>-</td></tr></tbody></table>

proteins) being systematically present in isolates of all three STs (Figure 3). ST496 isolates were characterized by the presence of the *spsI* gene. ST258 isolates were devoid of the *spsR* gene, but were the only ones to present the *nanB* gene coding for a potential sialidase toxin.

# DISCUSSION

## Decreasing Trend of the ST71 MRSP Lineage

In this study we showed a decreasing trend of ST71 MRSP, with a proportion falling from 65.3% in 2012–2013 to 55.2% in 2015–2016, among the French MRSP associated with canine SSTIs. This clear trend is nevertheless not statistically significant (p = 0.1643), due to the number of isolates included, but also mainly to the temporal proximity of the observed populations.

A previous French study reported that the proportion of ST71 isolates within MRSP was as high as 82.9% in 2010 (Haenni et al., 2014). The decrease in the ST71 clone between the 2010 study and the 2015–2016 collection is statistically significant (p = 0.0023), even though the number of isolates in the 2010 study was low (n = 41). Of note, the sampling scheme was not identical since all pathologies were included in the first surveillance study; however, SSTIs represented 63.4% of the reported infections. Taken together, our results corroborate the downward trend of ST71 isolates within the French population of MRSP over the last decade. The decrease of the ST71 MRSP lineage, which has already been reported in several countries from Northern Europe (Duim et al., 2016; Gronthal et al., 2017), thus most likely reflects of a global evolution of the population structure of MRSP in Europe. The number of isolates included in the present study may be considered a limitation. However, all MRSP isolates collected from SSTIs over the two periods

<footnote n="1">¹Antibiograms were performed by disk diffusion according to the recommendations of the French CA-SFM.</footnote>

<footnote n="1">¹Antibiograms were performed by disk diffusion according to the recommendations of the French CA-SFM.</footnote>
<figure><img src="image_3.png" /></figure>

were included and the same sampling effort was made in both periods. Considering the robustness of the collection and the number of other studies reporting the same trend in Northern countries, the tendencies detected here can be considered reliable in France.

## Population Structure of Non-ST71 Isolates

In both the 2012–2013 and 2015–2016 collections, as many as 21 distinct non-ST71 STs were identified, the vast majority being found only once. Analysis of the MRSP population structure based on the clonal complexes (CC) represented in our sampling revealed 2 different CCs and 12 singletons in the 2012–2013 collection (Supplementary Figure S1), and 1 CC and 14 singletons in the 2015–2016 collection (Supplementary Figure S2). The presence of many singletons and only a few CCs points out the coexistence of distinct and diverse MRSP population lineages. Identical non-ST71 lineages were rarely identified in both collections. The exceptions are: (i) ST265, which was identified eight times in 2012–2013 and three times

in 2015–2016 and thus seems to be on a decreasing trend; (ii) ST258, which was identified once in 2012–2013 and five times in 2015–2016 and, thus, seems to be expanding. Interestingly, a concomitant decrease of ST71 MRSP with the emergence of ST258 MRSP has been previously observed in Finland, Denmark, and Netherlands (Damborg et al., 2016; Duim et al., 2016; Gronthal et al., 2017).

## Clones and Antimicrobial Resistance

Our data show that, despite a higher tetracycline resistance, French ST258 MRSP isolates present a higher susceptibility to gentamicin and sulfonamides (absence of the *aacA-aphD* and *cat* genes), but also to fluoroquinolones. These results are consistent with the potential dissemination of a ST258 MRSP in Europe, which may be good news for veterinarians since this lineage is more susceptible than ST71 MRSP to veterinary-licensed antibiotics (Damborg et al., 2016; Duim et al., 2016; Gronthal et al., 2017). The SNP based phylogeny of the French ST71, ST258, ST496 and publicly available ST71, ST258, and ST496 genomes shows a very low diversity within the ST71 and ST496
<figure><img src="image_4.png" /><figcaption>**FIGURE 2** | SNP based phylogeny of ST71, ST258, and ST496 with a matrix of resistance associated genes and CRISPR. The SNP based phylogeny was constructed using 19008 SNP. The tip label represents the strain name followed by the ST number. The presence or absence of resistance genes and the CRISPR region are shown as an array of colored boxes on the right: ST71, red boxes; ST258, yellow boxes; and ST496, blue boxes.</figcaption></figure>

populations and a more heterogeneous ST258 population. The success of this newly emerging ST258 MRSP lineage remains unknown and cannot be explained by a difference in virulence factors since only *nanB*, a putative sialidase toxin, was identified to be exclusively present in the ST258 isolates. However, as only a few virulence factors have been reported and validated in *S. pseudintermedius*, one cannot exclude the presence or absence of other, yet undescribed, virulence associated genes within this lineage.

## Emergence of the ST496 Clone

In parallel with the emergence of a ST258 MRSP lineage, we also observed the emergence of a novel ST496 MRSP lineage in France in the collection of 2015–2016. All isolates belonging to this new clone were systematically resistant to all

veterinary-licensed antibiotics. Of note, over the 87 isolates from the 2015–2016 collection, only the nine ST496 and the three ST85 isolates presented such a large multi-resistance phenotype, with susceptibility present only for florfenicol and fusidic acid. Interestingly, this clone has very recently been reported in Australia (Worthing et al., 2018), where it is described as one of the major MRSP clones, especially in Sydney. Whole-genome comparison of French and Australian isolates shows a very low level of population diversity within this clade, the same resistome and a very close virulome. Compared to isolates of the other two most prominent MRSP lineages in France – ST71 and ST258 – ST496 isolates were the only ones that systematically carried cell wall associated (CWA) genes *spsI* and, except for one isolate, *spsF* genes, which are hypothesized to play a role in pathogenesis as they code for
<figure><img src="image_5.png" /><figcaption>**FIGURE 3** | SNP based phylogeny of ST71, ST258 and ST496 with a matrix of virulence associated genes. The SNP based phylogeny was constructed using 19008 SNP. The tip label represents the strain name followed by the ST number. The presence or absence of virulence genes are shown as an array of colored boxes on the right: ST71, red boxes; ST258, yellow boxes; and ST496, blue boxes.</figcaption></figure>

proteins capable of adhesion to extracellular matrix (Bannoehr et al., 2011).

## CONCLUSION

Put together, our data shows a shift, over the last 6 years, within the population structure of MRSP in France. An apparent decline of the ST71 MRSP lineage is concomitant with the emergence of two novel MRSP lineages in France: the ST258 MRSP lineage from Northern Europe, and a ST496 MRSP lineage which has become the main MRSP clone in Sydney, Australia. Based on the structure of the phylogenetic tree, the ST496 isolates from Australia appear to be ancestral to the French isolates. This is also indicated by the fact that ST496 emerged in Australia in

2013, but was not found in France at that time. The emergence of the ST258 and ST496 lineages in France might be associated with the distribution and importation of animals from other countries. The reasons for their apparent and simultaneous success remain unclear. For ST496, this could be associated with an increased colonization capacity, via the presence of *spsi* and *spsf* genes, but the expression and role of these putative microbial surface components recognizing adhesive matrix molecules (MSCRAMMs) remains to be validated *in vivo*. For ST258, no known virulence factors were exclusively found in this lineage which might help elucidate its success over the ST71 MRSP lineage.

While the decrease in the prevalence of ST71 MRSP is obvious in Europe, it is unclear which other MRSP clone(s) – and for which reasons – will successfully arise. This will have to be closely
monitored since the emergence and dissemination of the ST496 MRSP would be a threat to animal health due to its concerning level of antimicrobial resistance.

## AUTHOR CONTRIBUTIONS

MH, MB, and PM-S designed the experiments and drafted the manuscript. HK and PC performed the laboratory experiments. MB and PM-S performed the genomic analysis. J-YM, FL, KW, and JN actively contributed to the manuscript's writing. All authors approved the final version of this manuscript.

## REFERENCES

Bankevich, A., Nurk, S., Antipov, D., Gurevich, A. A., Dvorkin, M., Kulikov, A. S., et al. (2012). SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. *J. Comput. Biol.* 19, 455–477. doi: 10.1089/cmb.2012.0021

Bannoehr, J., Ben Zakour, N. L., Reglinski, M., Inglis, N. F., Prabhakaran, S., Fossum, E., et al. (2011). Genomic and surface proteomic analysis of the canine pathogen *Staphylococcus pseudintermedius* reveals proteins that mediate adherence to the extracellular matrix. *Infect. Immun.* 79, 3074–3086. doi: 10.1128/IAI.00137-11

Bannoehr, J., Brown, J. K., Shaw, D. J., Fitzgerald, R. J., Van Den Broek, A. H., and Thoday, K. L. (2012). *Staphylocococcus pseudintermedius* surface proteins SpsD and SpsO mediate adherence to ex vivo canine corneocytes. *Vet. Dermatol.* 23:e126. doi: 10.1111/j.1365-3164.2011.01021.x

Ben Zakour, N. L., Bannoehr, J., Van Den Broek, A. H., Thoday, K. L., and Fitzgerald, J. R. (2011). Complete genome sequence of the canine pathogen *Staphylococcus pseudintermedius*. *J. Bacteriol.* 193, 2363–2364. doi: 10.1128/JB.00137-11

Blaiotta, G., Fusco, V., Ercolini, D., Pepe, O., and Coppola, S. (2010). Diversity of Staphylococcus species strains based on partial kat (catalase) gene sequences and design of a PCR-restriction fragment length polymorphism assay for identification and differentiation of coagulase-positive species (*S. aureus*, *S. delphini*, *S. hyicus*, *S. intermedius*, *S. pseudintermedius*, and *S. schleiferi* subsp. *coagulans*). *J. Clin. Microbiol.* 48, 192–201. doi: 10.1128/JCM.00542-09

Bolger, A. M., Lohse, M., and Usadel, B. (2014). Trimmomatic: a flexible trimmer for Illumina sequence data. *Bioinformatics* 30, 2114–2120. doi: 10.1093/bioinformatics/btu170

Chen, L., Yang, J., Yu, J., Yao, Z., Sun, L., Shen, Y., et al. (2005). VFDB: a reference database for bacterial virulence factors. *Nucleic Acids Res.* 33, D325–D328.

Couto, N., Belas, A., Oliveira, M., Almeida, P., Clemente, C., and Pomba, C. (2016). Comparative RNA-seq-based transcriptome analysis of the virulence characteristics of methicillin-resistant and -susceptible *Staphylococcus pseudintermedius* strains isolated from small animals. *Antimicrob. Agents Chemother.* 60, 962–967. doi: 10.1128/AAC.01907-15

Damborg, P., Moodley, A., Aalbaek, B., Ventrella, G., Dos Santos, T. P., and Guardabassi, L. (2016). High genotypic diversity among methicillin-resistant *Staphylococcus pseudintermedius* isolated from canine infections in Denmark. *BMC Vet. Res.* 12:131. doi: 10.1186/s12917-016-0756-y

Descloux, S., Rossano, A., and Perreten, V. (2008). Characterization of new staphylococcal cassette chromosome mec (SCCmec) and topoisomerase genes in fluoroquinolone- and methicillin-resistant *Staphylococcus pseudintermedius*. *J. Clin. Microbiol.* 46, 1818–1823. doi: 10.1128/JCM.02255-07

Devriese, L. A., Vancanneyt, M., Baele, M., Vaneechoutte, M., De Graef, E., Snauwaert, C., et al. (2005). *Staphylococcus pseudintermedius* sp. nov., a coagulase-positive species from animals. *Int. J. Syst. Evol. Microbiol.* 55, 1569–1573. doi: 10.1099/ijs.0.63413-0

Duim, B., Verstappen, K. M., Broens, E. M., Laarhoven, L. M., Van Duijkeren, E., Hordijk, J., et al. (2016). Changes in the population of methicillin-resistant *Staphylococcus pseudintermedius* and dissemination of antimicrobial-resistant phenotypes in the Netherlands. *J. Clin. Microbiol.* 54, 283–288. doi: 10.1128/JCM.01288-15

## FUNDING

This work was supported by the French Agency for Food, Environment and Occupational Health & Safety (Anses) and Santé Publique France (SPF).

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2018.03055/full#supplementary-material

Francisco, A. P., Bugalho, M., Ramirez, M., and Carrico, J. A. (2009). Global optimal eBURST analysis of multilocus typing data using a graphic matroid approach. *BMC Bioinformatics* 10:152. doi: 10.1186/1471-2105-10-152

Garbacz, K., Zarnowska, S., Piechowicz, L., and Haras, K. (2013). Pathogenicity potential of *Staphylococcus pseudintermedius* strains isolated from canine carriers and from dogs with infection signs. *Virulence* 4, 255–259. doi: 10.4161/ viru.23526

Gortel, K., Campbell, K. L., Kakoma, I., Whittem, T., Schaeffer, D. J., and Weisiger, R. M. (1999). Methicillin resistance among staphylococci isolated from dogs. *Am. J. Vet. Res.* 60, 1526–1530.

Goy, M., Guindon, S., and Gascuel, O. (2010). SeaView version 4: a multiplatform graphical user interface for sequence alignment and phylogenetic tree building. *Mol. Biol. Evol.* 27, 221–224. doi: 10.1093/molbev/msp259

Grissa, I., Vergnaud, G., and Pourcel, C. (2007). CRISPRFinder: a web tool to identify clustered regularly interspaced short palindromic repeats. *Nucleic Acids Res.* 35, W52–W57. doi: 10.1093/nar/gkm360

Gronthal, T., Eklund, M., Thomson, K., Piiparinen, H., Sironen, T., and Rantala, M. (2017). Antimicrobial resistance in *Staphylococcus pseudintermedius* and the molecular epidemiology of methicillin-resistant *S. pseudintermedius* in small animals in Finland. *J. Antimicrob. Chemother.* 72, 1021–1030. doi: 10.1093/jac/ ddx086

Haenni, M., Alves De Moraes, N., Chatre, P., Médaille, C., Moodley, A., and Madec, J. Y. (2014). Characterization of clinical canine methicillin-resistant and methicillin-susceptible *Staphylococcus pseudintermedius* in France. *J. Glob. Antimicrob. Resistance* 2, 119–123. doi: 10.1016/j.jgar.2014.02.002

Kawakami, T., Shibata, S., Murayama, N., Nagata, M., Nishifuji, K., Iwasaki, T., et al. (2010). Antimicrobial susceptibility and methicillin resistance in *Staphylococcus pseudintermedius* and *Staphylococcus schleiferi* subsp. coagulans isolated from dogs with pyoderma in Japan. *J. Vet. Med. Sci.* 72, 1615–1619. doi: 10.1292/jvms.10-0172

Kjellman, E. E., Slettemeas, J. S., Small, H., and Sunde, M. (2015). Methicillin-resistant *Staphylococcus pseudintermedius* (MRSP) from healthy dogs in Norway - occurrence, genotypes and comparison to clinical MRSP. *Microbiologyopen* 4, 857–866. doi: 10.1002/mbo.258

Latronico, F., Moodley, A., Nielsen, S. S., and Guardabassi, L. (2014). Enhanced adherence of methicillin-resistant *Staphylococcus pseudintermedius* sequence type 71 to canine and human corneocytes. *Vet. Res.* 45:70. doi: 10.1186/1297-9716-45-70

Lozano, C., Rezusta, A., Ferrer, I., Perez-Laguna, V., Zarazaga, M., Ruiz-Ripa, L., et al. (2017). *Staphylococcus pseudintermedius* human infection cases in Spain: dog-to-human transmission. *Vector Borne Zoonotic Dis.* 17, 268–270. doi: 10.1089/vbz.2016.2048

Maes, N., Magdalena, J., Rottiers, S., De Gheldre, Y., and Struelens, M. J. (2002). Evaluation of a triplex PCR assay to discriminate *Staphylococcus aureus* from coagulase-negative staphylococci and determine methicillin resistance from blood cultures. *J. Clin. Microbiol.* 40, 1514–1517. doi: 10.1128/JCM.40.4.1514-1517.2002

Mikheenko, A., Saveliev, V., and Gurevich, A. (2016). MetaQUAST: evaluation of metagenome assemblies. *Bioinformatics* 32, 1088–1090. doi: 10.1093/bioinformatics/btv697

Monecke, S., Slickers, P., and Ehricht, R. (2008). Assignment of Staphylococcus aureus isolates to clonal complexes based on microarray analysis and pattern
recognition. *FEMS Immunol. Med. Microbiol.* 53, 237-251. doi: 10.1111/j.1574-695X.2008.00426.x

Nascimento, M., Sousa, A., Ramirez, M., Francisco, A. P., Carrico, J. A., and Vaz, C. (2017). PHYLOViZ 2.0: providing scalable data integration and visualization for multiple phylogenetic inference methods. *Bioinformatics* 33, 128-129. doi: 10.1093/bioinformatics/btw582

Osland, A. M., Vestby, L. K., Fanuelsen, H., Slettemeas, J. S., and Sunde, M. (2012). Clonal diversity and biofilm-forming ability of methicillin-resistant *Staphylococcus pseudintermedius*. *J. Antimicrob. Chemother.* 67, 841-848. doi: 10.1093/jac/dkr576

Perreten, V., Kadlec, K., Schwarz, S., Groenlund Andersson, U., Finn, M., Greko, C., et al. (2010). Clonal spread of methicillin-resistant *Staphylococcus pseudintermedius* in Europe and North America: an international multicentre study. *J. Antimicrob. Chemother.* 65, 1145-1154. doi: 10.1093/jac/dkq078

Pietrocola, G., Gianotti, V., Richards, A., Nobile, G., Geoghegan, J. A., Rindi, S., et al. (2015). Fibronectin binding proteins SpsD and SpsL both support invasion of canine epithelial cells by *Staphylococcus pseudintermedius*. *Infect. Immun.* 83, 4093-4102. doi: 10.1128/IAI.00542-15

Pires Dos Santos, T., Damborg, P., Moodley, A., and Guardabassi, L. (2016). Systematic review on global epidemiology of methicillin-resistant *Staphylococcus pseudintermedius*: inference of population structure from multilocus sequence typing Data. *Front. Microbiol.* 7:1599. doi: 10.3389/fmicb.2016.01599

Pomba, C., Rantala, M., Greko, C., Baptiste, K. E., Catry, B., Van Duijkeren, E., et al. (2017). Public health risk of antimicrobial resistance transfer from companion animals. *J. Antimicrob. Chemother.* 72, 957-968. doi: 10.1093/jac/dkw481

Price, M. N., Dehal, P. S., and Arkin, A. P. (2010). FastTree 2- approximately maximum-likelihood trees for large alignments. *PLoS One* 5:e9490. doi: 10.1371/journal.pone.0009490

Ruzauskas, M., Couto, N., Pavilonis, A., Klimiene, I., Siugzdiniene, R., Virgailis, M., et al. (2016). Characterization of *Staphylococcus pseudintermedius* isolated from diseased dogs in Lithuania. *Pol. J. Vet. Sci.* 19, 7-14. doi: 10.1515/pjvs-2016-0002

Seemann, T. (2014). Prokka: rapid prokaryotic genome annotation. *Bioinformatics* 30, 2068-2069. doi: 10.1093/bioinformatics/btu153

Solyman, S. M., Black, C. C., Duim, B., Perreten, V., Van Duijkeren, E., Wagenaar, J. A., et al. (2013). Multilocus sequence typing for characterization of *Staphylococcus pseudintermedius*. *J. Clin. Microbiol.* 51, 306-310. doi: 10.1128/JCM.02421-12

Somayaji, R., Rubin, J. E., Priyantha, M. A., and Church, D. (2016). Exploring *Staphylococcus pseudintermedius*: an emerging zoonotic pathogen? *Future Microbiol.* 11, 1371-1374. doi: 10.2217/fmb-2016-0137

Worthing, K. A., Abraham, S., Coombs, G. W., Pang, S., Saputra, S., Jordan, D., et al. (2018). Clonal diversity and geographic distribution of methicillin-resistant *Staphylococcus pseudintermedius* from Australian animals: discovery of novel sequence types. *Vet. Microbiol.* 213, 58-65. doi: 10.1016/j.vetmic.2017.11.018

Yu, G., Smith, D. K., Zhu, H., Guan, Y., and Lam, T. T. (2017). ggtree: an r package for visualization and annotation of phylogenetic trees with their covariates and other associated data. *Methods Ecol. Evol.* 8, 28-36. doi: 10.1111/2041-210X.12628

Zankari, E., Hasman, H., Cosentino, S., Vestergaard, M., Rasmussen, S., Lund, O., et al. (2012). Identification of acquired antimicrobial resistance genes. *J. Antimicrob. Chemother.* 67, 2640-2644. doi: 10.1093/jac/dks261

**Conflict of Interest Statement:** 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.

Copyright © 2018 Bergot, Martins-Simoes, Kilian, Châtre, Worthing, Norris, Madec, Laurent and Haenni. 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.