# Comprehensive genomic analysis of antibiotic resistance plasmids in animal-associated Staphylococcus aureus in France

Rachel Contarin, Séverine Murri, Antoine Drapeau, Tom Cayssials, Jean-Yves Madec, Emilie Dordet-Frisoni, Marisa Haenni

► To cite this version:

Rachel Contarin, Séverine Murri, Antoine Drapeau, Tom Cayssials, Jean-Yves Madec, et al. Comprehensive genomic analysis of antibiotic resistance plasmids in animal-associated Staphylococcus aureus in France. Microbiology Spectrum, 2025, 13 (10), 10.1128/spectrum.00772-25. anses-05316144

HAL Id: anses-05316144

https://anses.hal.science/anses-05316144v1

Submitted on 15 Oct 2025

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
# Comprehensive genomic analysis of antibiotic resistance plasmids in animal-associated *Staphylococcus aureus* in France

Rachel Contarin,¹,² Séverine Murri,² Antoine Drapeau,² Tom Cayssials,¹ Jean-Yves Madec,² Emilie Dordet-Frisoni,¹ Marisa Haenni²

**AUTHOR AFFILIATIONS** See affiliation list on p. 16.

**ABSTRACT** In *Staphylococcus aureus*, an animal pathogen and zoonotic agent, plasmids play a pivotal role in the acquisition and spread of antibiotic resistance genes (ARGs). This study investigated the plasmid content of 329 *S. aureus* isolates from livestock and companion animals collected in France between 2010 and 2021. Plasmids (n = 211) were identified from 139 isolates. The major families identified—*rep7a*, *rep20*, and *rep10*—were associated with specific resistance genes (*str*, *cat*, *blaZ*, *erm(C)*) and exhibited widespread horizontal transfer across different *S. aureus* sequence types (STs) and animal hosts. In temporal analysis, the *rep7a/str* and *rep7a/cat* plasmids circulating in horses were progressively replaced by a *rep7a* plasmid carrying both *str* and *cat* genes. The study also highlighted the presence of mosaic plasmids, which combined elements from different bacterial species/genera, confirming the broad host range of *S. aureus* plasmids and their ability to acquire ARGs from diverse sources. Moreover, the occurrence of hybrid plasmids (carrying multiple *rep* genes) underscores the plasticity of these vectors of ARGs. This study emphasizes the need to investigate the mechanisms driving the spread and persistence of antibiotic-resistant plasmids in *S. aureus*, with a view to developing strategies aimed at combating antibiotic resistance.

**IMPORTANCE** The spread of antibiotic resistance in *Staphylococcus aureus* is a growing concern, particularly in animals that can serve as reservoirs for resistant strains. This study highlights the crucial role of plasmids in transmitting resistance genes among different animal hosts and *S. aureus* lineages. The characterization of 329 isolates collected over 10 years revealed how certain plasmid families are associated with specific resistance genes and how they evolve over time. The occurrence of mosaic and hybrid plasmids further underscores the ability of *S. aureus* to acquire resistance from diverse bacterial sources. These findings provide key insights into the mechanisms shaping antibiotic resistance in this pathogen and emphasize the fact that understanding plasmid-driven resistance is essential for developing effective interventions to limit the spread of multidrug-resistant *S. aureus* in both veterinary and human medicine.

**KEYWORDS** *Staphylococcus aureus*, plasmids, antibiotic resistance, MGE, HGT, veterinary

*Staphylococcus aureus* is an opportunistic pathogen that triggers infections and diseases in numerous animal hosts (1) and poses a significant zoonotic threat. As a result, it is classified among high-priority bacteria by the World Health Organization (WHO) and the World Organization for Animal Health (WOAH). Under selective pressure, *S. aureus* can acquire numerous antibiotic resistance genes (ARGs) (2), the best known being the *mecA* gene, which confers resistance to all β-lactam antibiotics. Resistance gene dissemination is often mediated by mobile genetic elements (MGEs), which account for 15–20% of the *S. aureus* genome (3). In a recent study, plasmids, as key components of the *S. aureus* mobilome, were identified as the primary carriers of ARGs

**Editor** Nilton Lincopan, Universidade de Sao Paulo,
Sao Paulo, Brazil

Address correspondence to Marisa Haenni,
marisa.haenni@anses.fr, or Emilie Dordet-Frisoni,
emilie.dordet-frisoni@inrae.fr.

Emilie Dordet-Frisoni and Marisa Haenni contributed
equally to this article.

The authors declare no conflict of interest.

See the funding table on p. 16.

**Received** 13 March 2025

**Accepted** 22 August 2025

**Published** 18 September 2025

Copyright © 2025 Contarin et al. This is an open-
access article distributed under the terms of the
Creative Commons Attribution 4.0 International
license.
(4). These MGEs are small extrachromosomal circular DNA molecules that can repli-
cate independently of the host chromosome and efficiently transfer resistance and
virulence determinants among bacteria, even across species (3). Although extensively
studied in Enterobacterales, they have long been underestimated in Gram-positive
bacteria, including staphylococci.

Until recently, conjugation was considered a relatively rare event in *S. aureus*, as only about 5% of staphylococcal plasmids encode the machinery required for autonomous transfer (5); however, conjugative plasmids can facilitate the transfer of mobilizable plasmids present in the same cell (6, 7). Mobilizable plasmids either carry an *oriT* and encode a corresponding relaxase or only carry an *oriT* that is recognized by the relaxase of a conjugative plasmid. In order to be transferred, these relaxases must recognize the coupling protein of a conjugative plasmid, which actively promotes their transfer (2, 8, 9). Small plasmids < 45 kb, which are by far the most numerous in *S. aureus*, can also be transferred by generalized transduction, while natural genetic competence, although of lesser importance, might also play a role in plasmid transfer (10, 11).

Staphylococcal plasmids are classified into three classes according to their size and replication mechanism: (i) small (1.3–4.6 kb) multicopy plasmids, cryptic or carrying a single resistance determinant; (ii) larger (15–30 kb) low copy plasmids, which usually carry several resistance determinants; and (iii) conjugative (30–60 kb) multi-resistance plasmids (12–14). Larger plasmids undergo theta replication, whereas small plasmids usually replicate via an asymmetric rolling-circle mechanism (15, 16). Staphylococcal plasmids are also classified based on the nucleotide sequence of conserved regions of replication-initiating (*rep*) genes. Recurrent associations between *rep* genes and ARGs are described, such as *rep7* associated with tetracycline resistance, *rep10* associated with macrolide resistance, and *rep16* associated with β-lactam resistance (17). A PCR-based replicon typing (PBRT) scheme has been developed to classify plasmids into enterococci and staphylococci, encompassing 26 *rep* families and 10 unique families (18, 19). This classification is a useful tool for investigating plasmid dynamics within bacterial populations across different ecological niches (20).

The present study, in line with recent findings on the role of plasmids as key players
in the spread of ARGs, aimed to characterize the diversity of ARG-carrying plasmids
circulating among pathogenic *S. aureus* isolates from livestock and companion animals in
France. Plasmids from field *S. aureus* isolates collected through the Resapath network
(21) were fully sequenced, and genomic data were analyzed alongside associated
metadata, including host, date of sampling, and geographical origin. A subset of 81
isolates displaying multiple ARGs was further whole genome-sequenced to localize
precisely all of the identified ARGs. This approach enabled the identification of predomi-
nant and sporadic plasmid-ARG associations and provided insights into their dissemina-
tion dynamics.

MATERIALS AND METHODS

Bacterial isolates

This study included 329 *S. aureus* isolates of veterinary origin collected between 2010 and 2021 via the French Resapath network (21) across 70 different French departments. Isolates were collected from horses (n = 90), dogs (n = 87), cats (n = 83), cattle (n = 44), sheep (n = 11), goats (n = 3), primates (n = 4), poultry and birds (n = 3), and rodents (n = 4) (Table S1). The isolates included in this study were selected from a larger collection (n > 2,000 isolates). Among these >2,000 isolates, 525 were collected between 2010 and 2021. From these 525 isolates, our final collection of 329 isolates was chosen to be representative of the French territory, different years, and different resistance profiles found in this collection. All isolates were cultured on a COS medium (Biomérieux, Marcy l'Etoile, France) at 37°C.
## Antibiotic susceptibility testing

Susceptibility testing was performed using the disk diffusion method on Mueller-Hinton agar in accordance with the guidelines of the Antibiogram Committee of the French Society for Microbiology (22). *S. aureus* ATCC 25923 was used as the quality control strain. The antibiotics tested included penicillin G (PEN), kanamycin (KAN), gentamicin (GEN), tobramycin (TOB), chloramphenicol (CHL), florfenicol (FFL), tetracycline (TET), erythromycin (ERY), fusidic acid (FD), cefoxitin (FOX), spiramycin (SPI), lincomycin (LIN), and enrofloxacin (ENR) (Mast Diagnostics, Amiens, France).

## Detection of antimicrobial resistance genes

The 30 most common ARGs in *S. aureus* of animal origin were identified using the multiplex PCR technique (Table S2) (4). Genes not detected by PCR but identified through sequencing were included in the analysis. DNA extractions were performed using the boiling method, heating colonies at 100°C for one hour in one mL of NaCl 0.9%. PCRs were performed using the Qiagen Multiplex PCR Kit (Qiagen, Courtaboeuf, France) in a final volume of 25 µL. The following program was used for all multiplex PCRs: initial lysis and denaturation step at 95°C for 15 min, followed by 35 cycles at 94°C for 30 s, 55°C for 30 s (annealing), and 72°C for one minute. A final extension step was performed at 72°C for ten minutes. The amplicon sizes ranged from 158 to 1272 bp.

## Plasmid detection

The presence of plasmids was assessed using pulsed-field gel electrophoresis (PFGE) on S1-digested DNA (Fig. S1), as detailed in the Supplementary material.

For isolates in which no plasmid was detected by PFGE (n = 114), the presence of plasmids was verified by dot blot hybridization using a probe panel targeting conserved fragments of the *rep* genes labeled with digoxigenin by PCR using the primers listed in Table S3. Drops (4 µL) of boiled genomic DNA from each isolate were spotted in duplicate onto positively charged nylon membranes (Hybond N1; Amersham Biosciences, Orsay, France). Rep probes were hybridized and detected using a DIG-HIGH Prime DNA Labeling and Detection Starter Kit II (Roche Diagnostics, Meylan, France) following the manufacturer's instructions.

## Short- and long-read sequencing

A subset of 81 isolates chosen among the isolates carrying the largest number of ARGs and with diverse PFGE and resistance profiles (identified by PCR) was short-read whole genome-sequenced. DNA was extracted using the NucleoSpin Microbial DNA Extraction Kit (Macherey-Nagel, Hoerdt, France). For all other isolates presenting at least one plasmid via S1-PFGE (n = 142) or dot blot (n = 4), plasmid DNA was extracted using a Plasmid DNA Purification NucleoBond PC20 Kit (Macherey-Nagel) and long read-sequenced. Each identified plasmid was designated with a 'p' followed by the corresponding strain number. In cases where multiple plasmids were detected within the same strain, they were numbered in a decreasing order of size, with the largest plasmid assigned to number 1 and the smallest assigned to the highest number. Library preparation, short-read whole-genome sequencing (Illumina NovaSeq6000 technology), and long-read plasmid sequencing (GridION Oxford Nanopore Technology, Oxford Nanopore Technology, United Kingdom) were outsourced to Eurofins Genomics (Ebersberg, Germany). Short reads were quality trimmed and *de novo* assembled using Shovill v1.0.4, and the quality of assemblies was assessed using QUAST v5.0.2 (23) (Table S4). With regard to plasmids, *de novo* assembly was performed using either Unicycler v0.5.0 (24) or Flye v2.9.2 (25). Plasmids were annotated using Bakta 1.9.1 software (DB: 5.0.0).
## Characterization of plasmids and genomes and phylogenetic analysis

Plasmid replicons (*rep* coding gene) were detected using PlasmidFinder v.2.1.6 software (databases: v.2021-11-29) with default parameters (26, 27). Conjugation systems and proteins involved in mobilization were identified using CONJscan (28, 29), MOBscan (30), and Bakta 1.9.1 software (DB: 5.0.0). The presence of *oriT* was detected using the latter (DB: 5.0.0). ARGs were identified using ResFinder v.4.1.7 (database: v.2021-09-22) with default parameters (26, 31) and RGI-Comprehensive Antibiotic Resistance Database (CARD) (32). MinHash Approximation of SHared k-mers (MASH) v2.2.2 was used to calculate the distance between nucleotide sequences of plasmids using default parameters (33). Plasmid clustering was performed with thresholds of 0.18 for MASH clusters (A to N) and 0.06 for MASH groups defined by a numeric identifier (34). Results were visualized using iTOL v7 (https://itol.embl.de).

Genomes and MGEs, including integrative and conjugative elements (ICEs), transposons, composite transposons, SCC*mec*, and phages, were characterized as previously described (4). The *spa* type was determined using SpaTyper v0.3.3 (35). A neighbor-joining tree based on the *S. aureus* core genome (cgMLST), with core gene data available at https://www.cgmlst.org/ncs (36), was constructed using pyMLST v.2.1.3 with default parameters (37) and visualized using iTOL v7 (https://itol.embl.de).

## RESULTS

### Diversity of resistance genes in *S. aureus* isolated from animals

Analysis of the resistome of the 329 *S. aureus* veterinary clinical isolates collected between 2010 and 2021 via the Resapath network (21) showed that 92% (n = 303) carried at least one of the tested ARGs, with a maximum of 11 ARGs in a single isolate and an average of 4.3 genes per isolate. Among the 30 genes tested, 26 were detected (Table S2). The high prevalence of *mecA* (269/329; 82%) and *blaZ* (231/329; 70%), both conferring resistance to β-lactams, reflected a sampling bias due to the laboratory collection, which mostly counted MRSA isolates (Table S1; Fig. S2). In addition to these two genes, the *tet*(M) (tetracycline) (179/329; 54%), *aac(6′-aph(2′')* (aminoglycosides) (131/329; 40%), and *dfrK* (trimethoprim) (111/329; 34%) resistance genes were the most frequently identified (Table S1; Fig. S2). Overall, no significant host specificity was observed (Fig. S2), except for the *aac(6′-aph(2′')* gene, which was overrepresented in isolates from equines (91%) and primates (75%), and the *dfrK* gene, which was significantly more abundant in equine isolates (77%) (Fig. S2). The observed genotypes were consistent with the resistance phenotypes determined by disk diffusion (Table S1; Fig. S3).

### Plasmid diversity in *S. aureus* isolates of animal origin

PFGE on S1-digested DNA was performed on 256 isolates carrying at least one ARG without taking into account the chromosomally encoded *mecA* gene. Among them, 139 isolates (139/256; 54%) harbored at least one plasmid as determined by PFGE (Table S1). Long-read sequencing (Oxford Nanopore) of these 139 isolates revealed the presence of 211 plasmids carrying replicon proteins, since 54 isolates (39%) contained multiple plasmids (up to five, including six isolates presenting two plasmids of the *rep7a* family) (Table S1). Twenty-one hybrid plasmids also referred to as "multi-replicon plasmids" were identified, which carried two (n = 19) or three (n = 2) *rep* genes (Table S5). The 211 plasmids ranged in size from 1.4 to 50.6 kb, with 156 (74%) varying from 1.4 to 8.9 kb, and 55 (26%) from 19.7 to 50.6 kb (Fig. S4). Twenty-five isolates harboring multiple plasmids carried both large (>20 kb) and small (<10 kb) plasmids.

Resistome analysis of all plasmids revealed that 91% (n = 192/211) encoded at least one ARG (Table S5). MASH-based clustering grouped all plasmids into 13 heterogeneous clusters (named A to M) (Fig. 1), with clusters A to C being predominant (see sections below for detailed description). Most plasmids from clusters E to N either lacked ARGs or represented rare and atypical plasmids (Table S5).
<figure><img src="image_2.png" /><figcaption>**FIG 1** MASH analysis of the 211 plasmids from *S. aureus* isolates. This tree based on MASH distances was constructed using the neighbor-joining method and visualized using iTOL v7 (https://itol.embl.de). From the outer to the inner circle, the isolation year and the host of the isolate carrying a plasmid, size of the plasmid, and clusters are defined with a threshold of 0.18 MASH distance. The correspondence between colors and each host or cluster type is indicated in the top left corner.</figcaption></figure>

## Plasmids of the rep20 family

Cluster A comprised 48 plasmids (23%), most of which harbored the *rep20* gene (Table S5). Among these, four plasmids carried the *rep20* gene in association with *rep21*. Additionally, four plasmids carried multiple *rep* genes: one with a *rep16/rep19/rep7a* combination, another with *rep16/rep22/rep5a*, and two with *rep16/rep5a* (Table S5). Based on genomic distances calculated by k-mer analysis using MinHash Approximation of SHared k-mers (MASH) (33), cluster A was further divided into four groups (A1, A32, A4, and A33) according to the plasmid size, associated ARGs, and accessory genes present on the plasmid backbone (Fig. 2A).
<figure><img src="image_3.png" /></figure>

<figure><img src="image_4.png" /><figcaption>FIG 2 Plasmids of the *rep20* family. (A) DNA sequence comparison of plasmids belonging to the MASH cluster A (using the Clinker web-based tool https://cagecat.bioinformatics.nl/tools/clinker). The group to which each plasmid belongs and its name are indicated in front of each sequence. Only selected representatives from each group of MASH cluster A are shown. The nucleotide sequence identity between each coding region is represented with a color gradient ranging from white (0% homology) to black (100% homology). The scales for the gradient and size are located at the bottom. Genes of interest are</figcaption></figure>
Fig 2 (Continued)

highlighted in color to specify their function, which is indicated in the legend provided at the bottom. A few genes of interest are also specified. rel: relaxase; mob: mobilization protein. (B) Histogram showing the frequency of isolates carrying each plasmid from the *rep20* family, along with their associated ARGs, according to the year the plasmid-carrying strain was isolated. Each plasmid/ARG association is represented by a different color, as detailed in the legend at the bottom. In addition to the ARGs and *rep* families carried by the plasmids, the group and the number of isolates carrying each plasmid are indicated alongside each color.

The A1 group encompassed 27 large plasmids (25–27 kb), all of which carried *rep20*, *oriT*, and *blaZ* genes, except for the p42191 plasmid, which additionally carried the *tet*(K) gene, a relaxase MOBV, and belonged to *rep16/rep19/rep7a* families (Table S5). This *tet*(K) gene was located on a composite transposon formed by two *ISSau6* elements (Table S5). All plasmids in this A1 group carried the *cadD* gene (low-level cadmium resistance). Apart from p42191, they also carried the enterotoxin genes *sed*, *ser*, and *sej* (Fig. 2A), which are known to be plasmid-borne and are frequently reported in clinical isolates (38, 39). A1 group plasmids were identified between 2010 and 2019 (Fig. 2B) in both companion (cats, dogs, and a rabbit) and livestock (cattle) animals (Table S5; Fig. 1).

The A32 group comprised eight large (32 kb) plasmids carrying the *blaZ*, *sul1*, *dfrS1*, and *aac(6')-aph(2'*)* resistance genes (Table S5; Fig. 2A), as well as an *oriT* and a gene encoding a relaxase (Fig. 2A). This relaxase gene was flanked by two *IS256* insertion sequences, suggesting the formation of a composite transposon and the potential mobility of this region. These plasmids were almost exclusively detected in the *S. aureus* isolates of equine origin (n = 6) and observed in 2017, 2019, and 2020. In 2020, A32 group plasmids were the only ones detected from cluster A, suggesting a recent and potentially expanding emergence of this plasmid lineage in horse-associated isolates across France (Fig. 2B; Table S5).

The A4 group (n = 10) comprised plasmids carrying the *mer* operon (*merR*, *merA*, *merB*, *merT*), which was consistently flanked by two *IS431* insertion sequences, suggesting the potential mobility of this heavy metal resistance region (Fig. 2A). Based on plasmid content, this group could be further split into two subgroups. The first (n = 6) comprised 24 kb plasmids (e.g., p31933_1 and p29676_1 in Fig. 2A) carrying conjugal transfer protein, the *blaZ* gene, and harbored by strains from cattle, cats, and dogs from 2010 to 2017 (Table S5). The second group (n = 4) included larger *rep20/rep21* hybrid plasmids (31 kb) (e.g., p33658 and p32810_1 in Fig. 2A) displaying an *oriT*, the *blaZ* gene, and the *qacA* gene conferring resistance to biocides (Table S5). These plasmids were exclusively detected in companion animals between 2012 and 2018 and showed no evident temporal or spatial link with other *rep20* plasmids (Table S5; Fig. 2B).

The last group of cluster A, A33, included three hybrid plasmids, namely, p48613-1 (20 kb), p49684 (21 kb), and p48267 (28 kb), which all displayed an *oriT* and proteins involved in mobilization (Fig. 2A; Table S5). The two first plasmids carried the *rep16/rep5a* and the *blaZ* genes, while p48267 additionally harbored a genetic region carrying three other antibiotic resistance genes (*aph(3')-III*, *tet(L)*, and *sat-4*), along with a replication gene from the *rep22* family (Fig. 2A).

Globally, the proportion of plasmids carrying only the *blaZ* gene (groups A1 and A4) tended to decline from 25% in 2011 to less than 5% in 2019 and was no longer detected in 2020 in this data set (Fig. 2B). In parallel, since 2017, a new *rep20* plasmid (A32 group) carrying both *blaZ* and *aac(6')-aph(2'*)* genes has emerged across France, predominantly in *S. aureus* isolated from horses (Fig. 2B; Table S5).

## Plasmids of the *rep7a* family

Cluster B encompassed 90 plasmids (43%), making it the most prevalent cluster identified in this study. All of these plasmids carried only the *rep7a* gene and were small in size, ranging from 3.5 to 5 kb, with one outlier of 7 kb. This cluster was divided into six
groups based on MASH analysis (B7, B2, B3, B6, B13, and B21), taking into account both plasmid backbone structures and the ARGs they carried (Fig. 3; Table S5).

The B7 group was the largest (n = 60) and comprised plasmids sharing a high sequence similarity (Fig. 1 and 3A). These plasmids either carried only the *str* gene (n = 47) (e.g., p32084-1 and p42836-1) or both the *str* and *cat(pC221)* genes (n = 13) (e.g., p50407-1), conferring resistance to streptomycin and chloramphenicol, respectively (Fig. 3A; Table S5). All plasmids also encoded proteins involved in mobilization (MOBP), and most of them harbored an *oriT* (40/60) (Table S5; Fig. 3A). They were identified in different animal hosts, including cats, dogs, horses, cattle, and goats, over the entire study period from 2010 to 2020 (Table S5; Fig. 3B).

The B2 group encompassed 17 plasmids mostly carrying the *tet(K)* gene. Apart from p53983-2, all of them also encoded the mobilization protein MOBV (Table S5). One plasmid, p32762-1 (7 kb), also harbored the *spd* gene (*ant* (9) gene family) conferring resistance to spectinomycin, along with the insertion sequence *ISSau6* (Table S5; Fig. 3A). The B3, B6, B13, and B21 groups comprised two to five plasmids, which carried various ARGs, including *cat(pC221)* (e.g., p24614-2 from B3 and p32970-1 from B13), *str* (e.g., p33868 and p25658 from B6), *fosB* (e.g., p54235-2 from B3), *spd* (e.g., p36827-2 and p47948 from B21), and *qacC* (biocide resistance) (e.g., p41453-2 from B3) (Fig. 3A; Table S5).

Temporal analysis of this sampling of *rep7a* plasmids showed that, until 2012, the *str* and *cat(pC221)* genes were carried on distinct plasmids, primarily in equine isolates (Fig. 3B), with the *rep7a/cat(pC221)* plasmid mostly found in the Parisian region (Table S5). From 2013 onwards, a plasmid carrying both genes emerged in the western Paris area and progressively replaced the plasmids carrying individual resistance genes (Fig. 3B; Table S5).

### Plasmids associated with erm(C) genes: rep10 family and others

Cluster C encompassed 28 plasmids (13%) of very small size (2.4–2.8 kb), all of which belonged to the *rep10* family. These plasmids exhibited highly conserved sequences with no mobilization protein and no *oriT* and consistently carried only the *erm(C)* gene, conferring resistance to macrolides (Table S5). They were identified in companion animals, horses, and cattle from 2010 to 2020.

Apart from *rep10* plasmids, the *erm(C)* gene was also detected on hybrid plasmids belonging to cluster D (n = 15). Plasmids from the D8 group as defined by the MASH analysis (n = 11, 8.7–8.9 kb) harbored a *repUS18/repUS12* replication module and carried multiple resistance determinants, including the *erm(C)*, *aadD*, and *tet(L)* genes (Table S5). Notably, a majority of these plasmids from the D8 group harbored the protein mobilization MOBV (n = 7/11); they were isolated from cattle, dogs, and goats between 2010 and 2017 (Table S5).

### Atypical plasmids

The 30 plasmids belonging to clusters E to M (14.2%) exhibited significant sequence variability. Cluster E exclusively comprised plasmids harboring *rep13* or *rep21* genes, ranging in size from 2.2 to 3.3 kb (Table S5), and encoding two to six coding sequences. Only a few of these plasmids carried resistance genes: *cat(pC194)* was found in four *rep13* plasmids (e.g., p32151-2, p40695), *Inu(A)* in one (p33623-2), and *qacC*, *qacJ*, or *qacG* in four *rep21* plasmids (e.g., p35747, p32825-4, p42835-4, and p45322-2) (Table S5). The functions of the remaining coding sequences were mostly unknown.

Two large atypical ARG-carrying plasmids each constituted a unique cluster: p42837-1 (cluster L) and p33622-1 (cluster H) (Table S5; Fig. S3). The p42837-1 plasmid resembled those carrying the *cfr* gene (97.69% identity (85% query cover) with the pSA737 plasmid [40] and 97.12% (query cover 92%) with the p12-02300 plasmid of *Staphylococcus epidermidis* [41]) and thus likely belonged to the pWBG4-type conjugative plasmid family (5). This plasmid harbored the *Tn558* transposon carrying the *fexA* gene but lacked the
A

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

B

<figure><img src="image_6.png" /><figcaption>FIG 3 *rep7a* plasmid family. (A) DNA sequence comparison of plasmids belonging to the MASH cluster B (using the Clinker web-based tool https://cagecat.bioinformatics.nl/tools/clinker). The group to which each plasmid belongs and its name are indicated in front of each sequence. Only selected representatives from each group of MASH cluster B are shown. The nucleotide sequence identity between each coding region is represented with a color gradient ranging from white (Continued on next page)</figcaption></figure>
Fig 3 (Continued)

(0% homology) to black (100% homology). The scales for the gradient and size are located at the bottom. Genes of interest are highlighted in color to specify their function, which is indicated in the legend provided at the bottom. A few genes of interest are also specified. mob: mobilization protein. (B) Histogram showing the frequency of isolates carrying each plasmid from the *rep7a* family, along with their associated ARGs according to the year the plasmid-carrying strain was isolated. Each plasmid/ARG association is represented by a different color, as detailed in the legend at the bottom. In addition to the ARGs and *rep* families carried by the plasmids, the group and the number of isolates carrying each plasmid are indicated alongside each color.

typical 3 kb region flanking the *cfr* gene, with truncated *tnrB* and *bin3* genes (Fig. S5A). The p33622-1 plasmid appeared to be a hybrid between *S. aureus* and *Enterococcus faecium* plasmids (Fig. S5B). The successful circularization of the plasmid confirmed its mosaic structure and circular form, but the PlasmidFinder software used in this study was unable to detect any replication gene families. The *S. aureus* portion spanning over 30 kb with a 31% GC content was very similar to plasmid pM084526_1 (99.7% identity, 83% coverage for this region), originating from a human ST398 isolate (RIVM_M084526) (42) and carrying the *fexA* and *cfr* genes. However, in p33622-1, only the region containing the putative T4SS system was conserved, with no resistance gene retained (Fig. S5B). Additionally, in the *S. aureus* region, the plasmid included a PemIKSa toxin-antitoxin system previously described in plasmid pCH91 (43). In p33622-1, 38% of the sequence appeared to derive from *E. faecium* (CP091208.1) or *E. durans* (CP043327.1) plasmids, with >99.9% nucleotide identity. This region harbored the *erm(B)* gene, a zeta toxin-antitoxin system, and the IS1216 insertion sequence commonly associated with the *cfr* gene in *Enterococcus* plasmids (44).

Additional atypical plasmids were detected as well, including p42206 (*rep24b*), p54000-1 (*rep5c*), p36827-1 (ND, undetermined family), and p32664 (*rep5d*), none of which were associated with ARGs (Table S5). Notably, a conjugation system similar to that of plasmid pGO1 (45) (T4SS_typeFATA) was detected in plasmid p44206 (Table S5). These plasmids generally originated from human-associated *S. aureus* or, in the case of the latter two, from coagulase-negative staphylococci, which are known reservoirs of MGEs and potential sources of ARGs.

Plasmid dissemination between sequence types

A subset of 81 isolates that presented the highest number of plasmids carrying ARGs was whole genome-sequenced. A large proportion of the isolates (52/81; 64%) belonged to ST398, followed by ST8 (12/81; 15%) and ST5 (5/81; 6%) (Table S1; Fig. 4). The ST398 strains comprised isolates of eight different *spa* types, including *t011*, which is typically associated with horses, and carried two SCC*mec* types (IV and V) (Table S1). Within each ST, several of the ARG/plasmid associations described above (e.g., *str/rep7a*, *tet(K)/rep7a*, *erm(C)/rep10*, or *aadD/tet(L)/erm(C)/repUS18/repUS12*) were identified (Fig. 4; Table S5). For instance, ST8 isolates displayed 24 plasmids, and ST398 isolates carried 94, including primarily *erm(C)/rep10* (*n* = 8) and *str/rep7a* (*n* = 36), respectively (Table S5; Fig. 4). In addition to being clonally expanded within each ST, several ARG/plasmid associations were identified across multiple STs, including *erm(C)/rep10* in six STs, *str/rep7a* and *tet(K)/rep7a* in four, *blaZ/rep20* in three, and *aadD/erm(C)/tet(L)/repUS18/repUS12* in two (Table S5; Fig. 4).

Plasmid versus chromosomal location of antibiotic resistance genes

Whole-genome sequencing revealed distinct localization patterns for ARGs. Several genes, namely, *spd*, *str*, *Inu(A)*, *erm(C)*, *erm(T)*, *cat(pC221)*, *cat(pC194)*, and *qac*, were exclusively plasmid-borne (Table 1; Fig. S6A), while others like *aac(6')-aph(2'')*, *ant(6)-la*, *ant(9)-la*, *aph(3')-III*, *fusA*, *fusC*, *SA fosB*, *Inu(B)*, *lsa(E)*, *vga(A)V*, *dfrC*, and *dfrK*, were predominantly chromosomal (Table 1). The *aadD* gene was identified either on the chromosome or on plasmids: it was mostly plasmid-borne in whole-sequenced genomes (15/26; 58%) (Table 1) but was less (15/59, 25%) associated with plasmids among the 329 analyzed isolates (Fig. S6B). When plasmid-associated, it co-occurred with *erm(C)* and
<figure><img src="image_7.png" /><figcaption>FIG 4 Phylogenetic tree of 81 S. aureus isolates sequenced using Illumina technology and occurrence of the main plasmids carrying antibiotic resistance genes (ARGs). This tree based on the cgMLST distance was constructed using the neighbor-joining method and visualized using iTOL v7 (https://itol.embl.de). From the outer to the inner circle, the annotations represent the host, sequence type (ST), and diverse ARG/plasmid associations. The correspondence between colors and each host type, ST, or ARG/plasmid association is indicated on the right. The tree was rooted using the *Staphylococcus epidermidis* strain ATCC12228.</figcaption></figure>
<table><caption>TABLE 1 Association of antibiotic resistance genes (ARGs) with plasmids or chromosome and with other mobile genetic elements (MGEs)<sup>a</sup></caption><thead><tr><th>Antibiotic family</th><th>ARGs</th><th>Chromosome</th><th>Plasmid</th><th>Other MGEs</th></tr></thead><tbody><tr><td rowspan="10">Aminoglycoside</td><td>aac(6')-aph(2'')</td><td>33</td><td>1</td><td></td></tr><tr><td></td><td></td><td>2</td><td>cn_IS256</td></tr><tr><td>aadD</td><td>5</td><td>11</td><td></td></tr><tr><td></td><td></td><td>4</td><td>cn ISSau6</td></tr><tr><td></td><td>6</td><td></td><td>cn ISSau6</td></tr><tr><td>ant (6)-la</td><td>4</td><td></td><td></td></tr><tr><td>ant (9)-la</td><td>1</td><td></td><td></td></tr><tr><td></td><td>3</td><td></td><td>Tn554</td></tr><tr><td>aph(3')-III</td><td>2</td><td></td><td></td></tr><tr><td>sat-4</td><td>1</td><td></td><td></td></tr><tr><td>spd</td><td></td><td>3</td><td></td></tr><tr><td>str</td><td></td><td>49</td><td></td></tr><tr><td rowspan="3">Fosfomycin</td><td>fusA</td><td>2</td><td></td><td></td></tr><tr><td>fusC</td><td>1</td><td></td><td></td></tr><tr><td>SA fosB</td><td>22</td><td></td><td></td></tr><tr><td rowspan="4">Lincosamide</td><td>lnu(A)</td><td></td><td>1</td><td></td></tr><tr><td>lnu(B)</td><td></td><td></td><td></td></tr><tr><td>lsa(E)</td><td></td><td></td><td></td></tr><tr><td>vga(A)V</td><td>4</td><td></td><td>Tn5406</td></tr><tr><td rowspan="5">Macrolide</td><td>erm(A)</td><td>1</td><td></td><td></td></tr><tr><td></td><td>3</td><td></td><td>Tn554</td></tr><tr><td>erm(B)</td><td>1</td><td></td><td></td></tr><tr><td></td><td>1</td><td></td><td>Tn917</td></tr><tr><td>erm(C)</td><td></td><td>35</td><td></td></tr><tr><td></td><td>erm(T)</td><td></td><td>2</td><td>cn ISSau6</td></tr><tr><td></td><td></td><td></td><td>2</td><td></td></tr><tr><td rowspan="3">Phenicol</td><td>cat(pC194)</td><td></td><td>3</td><td></td></tr><tr><td>cat(pC221)</td><td></td><td>13</td><td></td></tr><tr><td>cat(pC233)</td><td>1</td><td></td><td></td></tr><tr><td></td><td>fexA</td><td>13</td><td>1</td><td>Tn558</td></tr><tr><td rowspan="3">Tetracycline</td><td>tet(K)</td><td>8</td><td>15</td><td></td></tr><tr><td></td><td></td><td>1</td><td>cn ISSau6</td></tr><tr><td>tet(L)</td><td>1</td><td>11</td><td></td></tr><tr><td></td><td></td><td>4</td><td>cn ISSau6</td></tr><tr><td rowspan="3">Trimethoprim</td><td>tet(M)</td><td>60</td><td></td><td>Tn916</td></tr><tr><td>dfrC</td><td>11</td><td></td><td>cn ISSau6</td></tr><tr><td>dfrK</td><td>15</td><td>2</td><td></td></tr><tr><td rowspan="4">β-Lactam</td><td></td><td>27</td><td></td><td>Tn559</td></tr><tr><td>blaZ</td><td>58</td><td>17</td><td></td></tr><tr><td></td><td>1</td><td></td><td>Tn552</td></tr><tr><td>mecA</td><td>69</td><td></td><td></td></tr><tr><td></td><td></td><td>7</td><td>cn ISSau6</td></tr><tr><td rowspan="4">Antiseptic</td><td>qacA</td><td></td><td>1</td><td></td></tr><tr><td>qacC</td><td></td><td>1</td><td></td></tr><tr><td>qacG</td><td></td><td>1</td><td></td></tr><tr><td>qacJ</td><td></td><td>2</td><td></td></tr></tbody></table>

<footnote><sup>a</sup>For each ARG, the number of genes associated with either the chromosome or plasmid is indicated. Other MGEs include composite transposons (elements beginning with "cn"), transposons (TnX), or an ICE (Tn916).</footnote>
tet(L). The *blaZ* gene was also largely chromosomal (58/75; 77%), with the plasmid variant (mostly on *rep20* plasmids) declining over time. Notably, chromosomal and plasmid *blaZ* variants differed in sequence (Table 1; Fig. S6C).

Some chromosomal ARGs were embedded in MGEs facilitating their spread. For example, *fexA* and *dfrK* were consistently linked to transposons *Tn558* and *Tn559*, respectively, regardless of the genomic location. ICEs and prophages were also abundant, but none carried ARGs: 83 ICEs (*Tn916/tet(M)* and ICE6013) were identified in 64 genomes (79%) (Table 1; Table S1), and 155 prophages from 18 families were found in 77 genomes (95%).

DISCUSSION

The role of plasmids in the dissemination of ARGs among Gram-positive bacteria remains less explored than that in Gram-negative bacteria. A recent study highlighted plasmids as the main vectors of ARGs in *S. aureus* (4). Since this study relied entirely on National Center for Biotechnology Information (NCBI) data sets, conclusions about the prevalence and evolution of these elements might suffer from the inherent limitations of this type of large public database (i.e., strong sampling biases, incomplete or imprecise metadata, and the frequent embedding of plasmid sequences within draft genome assemblies), thus making accurate plasmid identification and characterization more complicated. To better describe the diversity and dynamics of circulating plasmids, our study focused on a well-defined collection of 329 *S. aureus* isolates originating from animals in France and sampled over a 10-year period. Although geographically restricted to a single country and lacking the Bayesian inference required to formally estimate the timing of plasmid emergence or disappearance, this data set nevertheless provided valuable temporal insights into the local dynamics of plasmid dissemination, including their appearance, spread, and loss across animal-associated *S. aureus* isolates. Within this framework, a large diversity of ARGs was observed, the majority of them conferring resistances to β-lactams, aminoglycosides, and tetracyclines. This finding aligns with what has been observed in the large-scale study based on NCBI data (4) and is consistent with the fact that these three antibiotic classes have been the ones most commonly used for treating animals in France over the past decade (46). We further demonstrated that, in *S. aureus*, certain ARGs were consistently plasmid-borne, such as *str*, *erm(C)*, *tet(K)*, *tet(L)*, and *cat(pC221)*, while others were only occasionally carried on plasmids (e.g., *blaZ* and *aac(6')-aph(2'')*) or were predominantly chromosomal (e.g., *tet(M)*, *dfrC*, and *mecA*).

Although rare, transfer events from the chromosome to plasmids may be hypothe-
sized, as observed with *blaZ*-carrying *rep20* plasmids from ST612 equine isolates (e.g.,
48219, 54205, and 54213), which have recently acquired the *aac(6')-aph(2'')* gene.
This aminoglycoside-resistance gene predominantly identified as chromosomal in our
whole-sequenced genome and in NCBI data (4) may have been mobilized onto plasmids,
thanks to the presence of two flanking IS256 elements. This hypothesis was further
supported by the fact that the *aac(6')-aph(2'')* gene sequences found on either the
chromosome or plasmids were identical, potentially indicating recent transfer events.
Conversely, while a decrease in plasmid-borne *blaZ* and a concomitant increase in its
chromosomal localization were observed over time, our data do not support recent
transfers from plasmids to the chromosome. Indeed, the chromosomal *blaZ* variant
differed in sequence coverage and identity (94% for both) from the plasmid-borne
version identified on the 27 kb *rep20* plasmid (e.g., 24506, 27262, and 32722). Sequence
divergences reflect different evolutionary paths depending on the plasmid or chromo-
somal location (47), and the decrease in *blaZ* observed in animals in France is most
probably due to modifications in the staphylococcal lineages circulating in the country.

The analysis of 211 plasmids from animal-associated *S. aureus* isolates from France
revealed substantial diversity, yet three replicon families clearly dominated: *rep7a* (43%),
*rep20* (19%), and *rep10* (13%). Strong associations between specific plasmid backbones
and ARGs were observed, even though some genes were identified on two distinct
plasmid families, namely, *erm(C)* on *rep10* and *repUS18/repUS12*, without any association
with other mobile genetic elements, highlighting the remarkable ability of certain genes to disseminate efficiently on different genetic determinants (48). Each of the three major plasmid clusters was found across multiple STs and animal hosts, thereby supporting the hypothesis of horizontal plasmid transfer between diverse *S. aureus* lineages. Indeed, our analysis revealed that certain ARG/plasmid associations were not restricted to specific lineages but were found across multiple STs. While the spread of some plasmids reflected the clonal success of the ST (such as ST398 or ST8) by which they are harbored, others like *erm(C)/rep10* (found in six STs) or *str/rep7a* and *tet(K)/rep7a* (found in four STs) clearly indicated inter-ST dissemination. These plasmids, especially *rep7a* and *rep10* families, appeared to be efficient vectors for ARGs, contributing to their widespread presence in both livestock- and companion animal-associated isolates. The detection of identical plasmid profiles in unrelated STs and hosts further highlights the mobility and ecological versatility of plasmids, underscoring the role of plasmids as key drivers of resistance spread across genetic and host boundaries.

Ubiquitous small plasmids (<10 kb) comprised the majority of the *S. aureus* plasmidome—in line with previous reports (49–53)—whereas large plasmids (>20 kb) were relatively rare, representing less than a third of our collection. This distribution highlights the significant role of small plasmids in the ecology and evolution of *S. aureus* likely driven by rapid replication and unexpected transfer efficiency despite the absence of dedicated conjugation machinery (19). Among the small plasmids commonly observed in *S. aureus*, the *rep10* family (cluster C) stands out as the most conserved. This plasmid family is consistently associated with the *erm(C)* gene, which is invariably linked to a replicase encoded by the *rep10* gene, as already suggested in a large-scale study of over 10,000 genomes of *S. aureus* (4). These plasmids have been repeatedly reported not only in *S. aureus* (19, 51, 52, 54) but also, albeit less frequently, in other Firmicutes, such as *S. epidermidis* (18) and *Enterococcus faecalis* (55). Despite lacking canonical mobilization genes, *rep10* plasmids have been shown to transfer via transformation or transduction between *Staphylococcus* species and even across genera (56–58). These alternative mechanisms may account for their broad distribution across diverse STs and host origins.

This study also revealed the prevalence of *rep7a* plasmids (cluster B), which, despite their small size, were associated with multiple ARGs and harbored conserved backbones encoding mobilization functions, suggesting that they may act as mobilizable elements. These plasmids were frequently detected in strains also carrying larger plasmids—most commonly *rep20* family members. They have been previously described in *S. aureus* (19) and are commonly referred to by three prototypic names in the literature: (i) pT181, which carries the *tet(K)* gene, and a MOBV relaxase (59, 60), which corresponds to plasmids of the B2 group in our study; (ii) pC221, which displays the *cat(pC221)* gene, a replication gene, and several genes of unknown function (54) and is similar to p24614-2 (B3 group); and (iii) pS194, which harbors the *str* gene and a relaxase (61) and resembles p32084-1 (B7 group). In our study, temporal analysis of *rep7a* plasmids further revealed a significant evolution in both the *S. aureus* plasmidome and its associated resistome. Interestingly, the *rep7a* plasmid carrying both *cat(pC221)* and *str* (e.g., p50407 from B7 group) emerged in equine isolates, particularly in the Paris region and surrounding areas. In contrast, other *rep7a* plasmid variants, as well as *rep10* and *rep20* families, showed no apparent geographical clustering and were distributed more broadly across France. The high sequence similarity between *rep7a/cat(pC221)* (e.g., p32970-1, B13 group) and *rep7a/str* (e.g., 32084-1, B7 group) plasmids coupled with their rolling-circle replication may have facilitated recombination events leading to the formation of multi-resistance plasmids. Recombination has already been described as a key driver of plasmid evolution and adaptation, particularly in Gram-negative bacteria (62–64).

Apart from small and highly recombinant plasmids, we also observed significant genetic plasticity among larger plasmids (20–32 kb in size), particularly within the *rep20* plasmid family (cluster A), which exhibited marked heterogeneity despite being primarily associated with the *blaZ* gene. These plasmids were identified not only in
livestock-associated *S. aureus* like ST398, but also in human-associated STs, such as ST8 and ST5, that were collected from cats and dogs. In the latter ST, *rep20* plasmids often appeared as hybrid plasmids frequently carrying an additional replicon, such as *rep5a* (49). The most abundant group within cluster A, A1, carried only *blaZ* and shared 99.9% identity with the 27 kb plasmid p19321-P03 from strain 19321, a community-associated MRSA isolate (65). It also closely resembled pSK67, an *S. aureus* plasmid historically reported in Australian hospitals between 1946 and 1981 (66). In addition to ARGs, these plasmids also harbored enterotoxin genes (*ser*, *sej*, *sep*) commonly associated with food-borne *S. aureus* outbreaks (67, 68). In this study, group A1 was exclusively detected in companion animal isolates (dogs, cats, and one rabbit) primarily belonging to ST8, suggesting possible human-to-animal transmission through close contact.

Our results showed that a single isolate often carried more than one plasmid, usually one large (mainly *rep20*, cluster A) and one small (often from the *rep10* or *rep7a* plasmids of clusters C and B, respectively). This pattern has also been observed in human-derived MRSA isolates in Malaysia, where large plasmids were found to coexist with *rep10* plasmids within the same genome (49). These findings are further supported by a large-scale analysis of *S. aureus* genomes from the NCBI database, which highlighted the frequent co-occurrence of multiple plasmid types (4). A few studies explicitly reported the coexistence of a large plasmid alongside a smaller one. However, large conjugative plasmids were capable of mobilizing smaller nonconjugative ones in *S. aureus* (5, 69). Although only one conjugative plasmid was formally identified in our data set, this absence likely reflects limitations in current bioinformatic tools rather than a true absence of such elements. Moreover, the coexistence of multiple plasmids within a single isolate facilitates the formation of hybrid plasmids, which have been reported in previous studies (70). For example, in plasmids from the D8 group, the association of *repUS12* and *repUS18* was linked with *aadD*, *erm(C)*, and *tet(L)* genes. The *rep20* replicon was also found in association with *rep21* in hybrid plasmids belonging to the A4 group, as previously described (71). The *rep21* gene encodes a Rep1-type replication protein typically found in small cryptic plasmids (<3 kb), which have been identified in clinical MRSA human isolates (49) and are prone to integration into larger plasmids (~30 kb).

Here, *rep20/rep21* plasmids closely resembled 24 kb A4 group plasmids carrying only *rep20* and *blaZ*, with the only differences being in the insertion of *rep21* near an origin of transfer and the *ISSau6* element, along with the insertion of the *qac* gene (e.g., p33658). This configuration suggests the integration of a small Rep1-type plasmid into a larger backbone. The Rep1 domain of replication proteins has been shown to act as dual-function replication-relaxase proteins supporting the mobilization of otherwise nonconjugative plasmids (72, 73). Lee et al. (72) demonstrated that in *Bacillus subtilis*, ICE elements can mediate the transfer of Rep1-type plasmids through a relaxase-independent mechanism (72). In our study, relaxase activity was also predicted in certain *repUS18/repUS12* plasmids that encode replication proteins with Inc18 and Rep1 domains, respectively, and were found to carry *aadD*, *tet(L)*, and *erm(C)* genes. The presence of such mosaic plasmids comprising multiple replicons and ARGs underscores the remarkable plasticity and adaptive potential of *S. aureus* plasmidome.

Finally, the remarkable plasticity of *S. aureus* plasmids was further evidenced by the identification of an atypical mosaic plasmid, p33622-1, composed of backbone elements originating from both *S. aureus* and *Enterococcus faecium*. Although the *cfr* gene commonly associated with resistance to linezolid was not present in the plasmids identified in our isolates, this hybrid structure suggests that plasmid-mediated exchange of the *cfr* gene between species is possible. Supporting this, a previous study demonstrated the successful transfer of a *cfr*-carrying plasmid from *Enterococcus* to *S. aureus* (74). Recombination between these two plasmids may have occurred through the IS shared by both plasmids, facilitating inter-species genetic exchange. These findings further confirm the broad host range and dynamic nature of *S. aureus* plasmids, which can expand their reservoir of ARGs through horizontal gene transfer across species. Such mosaic plasmids enhance the genetic diversity and evolutionary potential of *S. aureus*,
emphasizing the importance of understanding plasmid mobility and recombination in the context of antimicrobial resistance dissemination.

## Conclusion

This study demonstrated the significant diversity of plasmids in *S. aureus* isolates of animal origin, with the predominance of the *rep7a*, *rep20*, and *rep10* families. Despite their relatively small size, these plasmids were widespread and driven by their genetic plasticity, broad host range, and capacity for mobilization and homologous recombination.

While plasmid dissemination in *S. aureus* has traditionally been linked to epidemic clones, our results suggest that plasmids are also spreading horizontally. This was highlighted in this longitudinal analysis by the detection of specific plasmids, including small multi-resistance ones, in different STs. Overall, our analyzes underscored the notion that plasmids most likely act as dynamic vehicles for ARGs in *S. aureus*. These findings highlight the critical need to unravel the mechanisms driving the epidemic success and widespread dissemination of ARG-carrying plasmids in *S. aureus*.

## ACKNOWLEDGMENTS

This work was supported by a PhD grant from the French National Research Institute for Agriculture, Food and the Environment (INRAE) and French National Agency for Food, Environmental and Occupational Health and Safety (ANSES).

We thank Pierre Châtre (ANSES Lyon) for his valuable assistance in conducting the PFGE. We also thank Pauline François (ANSES Lyon) for her guidance with bioinformatics and sequencing and Véronique Métayer (ANSES Lyon) for her help in performing the antibiograms and PCRs. This manuscript has has been professionally proofread (Proof-reading-Service, Hertfordshire, UK).

## AUTHOR AFFILIATIONS

## AUTHOR ORCIDs

Rachel Contarin http://orcid.org/0009-0006-2215-6622

Emilie Dordet-Frisoni http://orcid.org/0000-0002-6315-5670

Marisa Haenni http://orcid.org/0000-0003-1333-5603

## FUNDING

<table><thead><tr><td>Funder</td><td>Grant(s)</td><td>Author(s)</td></tr></thead><tbody><tr><td>ANSES/INRAE</td><td>PhD grant</td><td>Rachel Contarin</td></tr></tbody></table>

## AUTHOR CONTRIBUTIONS

Rachel Contarin, Formal analysis, Investigation, Methodology, Software, Writing – original draft | Séverine Murri, Investigation | Antoine Drapeau, Software | Tom Cayssials, Investigation | Jean-Yves Madec, Funding acquisition, Validation, Writing – review and editing | Emilie Dordet-Frisoni, Conceptualization, Formal analysis, Project administration, Supervision, Validation, Writing – review and editing | Marisa Haenni, Conceptualization, Funding acquisition, Project administration, Supervision, Validation, Writing – review and editing

<footnote n="1">¹INTHERES, Université de Toulouse, INRAE, ENVT, Toulouse, France</footnote>

<footnote n="2">²Anses-Université de Lyon, Unité Antibiorésistance et Virulence Bactériennes, Lyon, France</footnote>
DATA AVAILABILITY

All plasmid and chromosomal raw data are available on the NCBI's Sequence Read Archive (SRA) under BioProject PRJNA1193534.

ADDITIONAL FILES

The following material is available online.

Supplemental Material

Figure S1 (Spectrum00772-25-s0001.pdf). PFGE profiles.

Supplemental material (Spectrum00772-25-s0002.docx). Supplemental methods, Tables S2 and S3, and Fig. S2 to S6.

Table S1 (Spectrum00772-25-s0003.xlsx). List of strains and their metadata.

Table S4 (Spectrum00772-25-s0004.xlsx). Quality control of the 81 sequenced *S. aureus* genomes.

Table S5 (Spectrum00772-25-s0005.xlsx). List of plasmids and their metadata.

REFERENCES

1. HaagAF, FitzgeraldJR. 2019. *Staphylococcus aureus* in animals. Microbiol Spectr 7. https://doi.org/10.1128/microbiolspec.GPP3-0060-2019

2. Haaber J, Penadés JR, Ingmer H. 2017. Transfer of antibiotic resistance in *Staphylococcus aureus*. Trends Microbiol 25:893-905. https://doi.org/10.1016/j.tim.2017.05.011

3. Alibayov B, Baba-Moussa L, Sina H, Zdeňková K, Demnerová K. 2014. *Staphylococcus aureus* mobile genetic elements. Mol Biol Rep 41:5005-5018. https://doi.org/10.1007/s11033-014-3367-3

4. Contarin R, Drapeau A, François P, Madec J-Y, Haenni M, Dordet-Frisoni E. 2024. The interplay between mobilome and resistome in *Staphylococcus aureus*. mBio 15:e0242824. https://doi.org/10.1128/mbio.02428-24

5. Ramsay JP, Kwong SM, Murphy RJT, Yui Eto K, Price KJ, Nguyen QT, O'Brien FG, Grubb WB, Coombs GW, Firth N. 2016. An updated view of plasmid conjugation and mobilization in *Staphylococcus*. Mob Genet Elements 6:e1208317. https://doi.org/10.1080/2159256X.2016.1208317

6. Coluzzi C, Garcillán-Barcia MP, de la Cruz F, Rocha EPC. 2022. Evolution of plasmid mobility: origin and fate of conjugative and nonconjugative plasmids. Mol Biol Evol 39:msac115. https://doi.org/10.1093/molbev/msac115

7. Smillie C, Garcillán-Barcia MP, Francia MV, Rocha EPC, de la Cruz F. 2010. Mobility of plasmids. Microbiol Mol Biol Rev 74:434-452. https://doi.org/10.1128/MMBR.00020-10

8. Projan SJ, Archer GL. 1989. Mobilization of the relaxable *Staphylococcus aureus* plasmid pC221 by the conjugative plasmid pGO1 involves three pC221 loci. J Bacteriol 171:1841-1845. https://doi.org/10.1128/jb.171.4.1841-1845.1989

9. Smith MCA, Thomas CD. 2004. An accessory protein is required for relaxosome formation by small staphylococcal plasmids. J Bacteriol 186:3363-3373. https://doi.org/10.1128/JB.186.11.3363-3373.2004

10. McCarthy AJ, Lindsay JA. 2012. The distribution of plasmids that carry virulence and resistance genes in *Staphylococcus aureus* is lineage associated. BMC Microbiol 12:104. https://doi.org/10.1186/1471-2180-12-104

11. Morikawa K, Takemura AJ, Inose Y, Tsai M, Nguyen Thi LT, Ohta T, Msadek T. 2012. Expression of a cryptic secondary sigma factor gene unveils natural competence for dna transformation in *Staphylococcus aureus*. PLoS Pathog 8:e1003003. https://doi.org/10.1371/journal.ppat.1003003

12. Malachowa N, DeLeo FR. 2010. Mobile genetic elements of *Staphylococcus aureus*. Cell Mol Life Sci 67:3057-3071. https://doi.org/10.1007/s00018-010-0389-4

13. Novick RP. 1989. Staphylococcal plasmids and their replication. Annu Rev Microbiol 43:537-565. https://doi.org/10.1146/annurev.mi.43.100189.9.02541

14. Partridge SR, Kwong SM, Firth N, Jensen SO. 2018. Mobile genetic elements associated with antimicrobial resistance. Clin Microbiol Rev 31:00088-17. https://doi.org/10.1128/CMR.00088-17

15. Khan SA. 2005. Plasmid rolling-circle replication: highlights of two decades of research. Plasmid 53:126-136. https://doi.org/10.1016/j.plasmid.2004.12.008

16. Lindsay JA. 2010. Genomic variation and evolution of *Staphylococcus aureus*. Int J Med Microbiol 300:98-103. https://doi.org/10.1016/j.ijmm.2009.08.013

17. Mores CR, Montelongo C, Putonti C, Wolfe AJ, Abouelfetouh A. 2021. Investigation of plasmids among clinical *Staphylococcus aureus* and *Staphylococcus haemolyticus* isolates from Egypt. Front Microbiol 12:659116. https://doi.org/10.3389/fmicb.2021.659116

18. Jensen LB, Garcia-Migura L, Valenzuela AJS, Løhr M, Hasman H, Aarestrup FM. 2010. A classification system for plasmids from enterococci and other Gram-positive bacteria. J Microbiol Methods 80:25-43. https://doi.org/10.1016/j.mimet.2009.10.012

19. Lozano C, García-Migura L, Aspiroz C, Zarazaga M, Torres C, Aarestrup FM. 2012. Expansion of a plasmid classification system for Gram-positive bacteria and determination of the diversity of plasmids in *Staphylococcus aureus* strains of human, animal, and food origins. Appl Environ Microbiol 78:5948-5955. https://doi.org/10.1128/AEM.00870-12

20. Orlek A, Stoesser N, Anjum MF, Doumith M, Ellington MJ, Peto T, Crook D, Woodford N, Walker AS, Phan H, Sheppard AE. 2017. Plasmid classification in an era of whole-genome sequencing: application in studies of antibiotic resistance epidemiology. Front Microbiol 8:182. http://doi.org/10.3389/fmicb.2017.00182

21. ANSES. 2024. Resapath. Available from: https://shiny-public.anses.fr/resa path2. Retrieved 07 Oct 2024.

22. Comité de l'Antibiograme de la Société Française de Microbiologie. 2024. Société Française de Microbiologie. Available from: https://www.sfm-microbiologie.org/boutique/_comite-de-lantibiogramme-de-la-sfm-c-a-sfm-v1-0-juin-2024. Retrieved 25 Nov 2024.

23. Gurevich A, Saveliev V, Vyahhi N, Tesler G. 2013. QUAST: quality assessment tool for genome assemblies. Bioinformatics 29:1072-1075. https://doi.org/10.1093/bioinformatics/btt086

24. Wick RR, Judd LM, Gorrie CL, Holt KE. 2017. Unicycer: resolving bacterial genome assemblies from short and long sequencing reads. PLoS Comput Biol 13:e1005595. https://doi.org/10.1371/journal.pcbi.1005595

25. Kolmogorov M, Yuan J, Lin Y, Pevzner PA. 2019. Assembly of long, error-prone reads using repeat graphs. Nat Biotechnol 37:540-546. https://doi.org/10.1038/s41587-019-0072-8

26. Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, Madden TL. 2009. BLAST+: architecture and applications. BMC Bioinformatics 10:421. https://doi.org/10.1186/1471-2105-10-421

27. Carattoli A, Zankari E, García-Fernández A, Voldby Larsen M, Lund O, Villa L, Møller Aarestrup F, Hasman H. 2014. *In silico* detection and typing of plasmids using PlasmidFinder and plasmid multilocus sequence typing. Antimicrob Agents Chemother 58:3895-3903. https://doi.org/10.1128/AAC.02412-14

28. Néron B, Denise R, Coluzzi C, Touchon M, Rocha EPC, Abby SS. 2023. MacSyFinder v2: improved modelling and search engine to identify molecular systems in genomes. Peer Community J 3. https://doi.org/10.24072/pcjournal.250
29. Cury J, Abby SS, Doppelt-Azeroual O, Néron B, Rocha EPC. 2020. Identifying conjugative plasmids and integrative conjugative elements with CONJscan, p 265–283. In de la Cruz F (ed), Horizontal gene transfer: methods and protocols. Springer US, New York, NY.

30. Garcillán-Barcia MP, Redondo-Salvo S, Vielva L, de la Cruz F. 2020. MOBscan: automated annotation of MOB relaxases. Methods Mol Biol 2075:295–308. https://doi.org/10.1007/978-1-4939-9877-7_21

31. Bortolaia V, Kaas RS, Ruppe E, Roberts MC, Schwarz S, Cattoir V, Philippon A, Allesoe RL, Rebelo AR, Florenza AF, et al. 2020. ResFinder 4.0 for predictions of phenotypes from genotypes. J Antimicrob Chemother 75:3491–3500. https://doi.org/10.1093/jac/dkaa345

32. Alcock BP, Huynh W, Chalil R, Smith KW, Raphenya AR, Wlodarski MA, Edalatmand A, Petkau A, Syed SA, Tsang KK, et al. 2023. CARD 2023: expanded curation, support for machine learning, and resistome prediction at the comprehensive antibiotic resistance database. Nucleic Acids Res 51:D690–D699. https://doi.org/10.1093/nar/gkac920

33. Ondov BD, Treangen TJ, Melsted P, Mallonee AB, Bergman NH, Koren S, Phillippy AM. 2016. Mash: fast genome and metagenome distance estimation using MinHash. Genome Biol 17:132. https://doi.org/10.1186/s13059-016-0997-x

34. Robertson J, Bessonov K, Schonfeld J, Nash JHE. 2020. Universal whole-sequence-based plasmid typing and its utility to prediction of host range and epidemiological surveillance. Microb Genom 6:e000435. https://doi.org/10.1099/mgen.0.000435

35. Bartels MD, Petersen A, Worning P, Nielsen JB, Larner-Svensson H, Johansen HK, Andersen LP, Jarløv JO, Boye K, Larsen AR, Westh H. 2014. Comparing whole-genome sequencing with Sanger sequencing for spa typing of methicillin-resistant *Staphylococcus aureus*. J Clin Microbiol 52:4305–4308. https://doi.org/10.1128/JCM.01979-14

36. Leopold SR, Goering RV, Witten A, Harmsen D, Mellmann A. 2014. Bacterial whole-genome sequencing revisited: portable, scalable, and standardized analysis for typing and detection of virulence and antibiotic resistance genes. J Clin Microbiol 52:2365–2370. https://doi.org/10.1128/JCM.00262-14

37. Biguenet A, Bordy A, Atchon A, Hocquet D, Valot B. 2023. Introduction and benchmarking of pyMLST: open-source software for assessing bacterial clonality using core genome MLST. Microb Genom 9:001126. https://doi.org/10.1099/mgen.0.001126

38. Varshney AK, Mediavilla JR, Robiou N, Guh A, Wang X, Gialanella P, Levi MH, Kreiswirth BN, Fries BC. 2009. Diverse enterotoxin gene profiles among clonal complexes of *Staphylococcus aureus* isolates from the Bronx, New York. Appl Environ Microbiol 75:6839–6849. https://doi.org/10.1128/AEM.00272-09

39. Shearer JES, Wireman J, Hostetler J, Forberger H, Borman J, Gill J, Sanchez S, Mankin A, Lamarre J, Lindsay JA, Bayles K, Nicholson A, O'Brien F, Jensen SO, Firth N, Skurray RA, Summers AO. 2011. Major families of multiresistant plasmids from geographically and epidemiologically diverse staphylococci. G3 (Bethesda) 1:581–591. https://doi.org/10.1534/g3.111.000760

40. Mendes RE, Deshpande LM, Castanheira M, DiPersio J, Saubolle MA, Jones RN. 2008. First report of *cfr*-mediated resistance to linololid in human staphylococcal clinical isolates recovered in the United States. Antimicrob Agents Chemother 52:2244–2246. https://doi.org/10.1128/AAC.00231-08

41. Bender J, Strommenger B, Steglich M, Zimmermann O, Fenner I, Lensing C, Dagwadordsch U, Kekulé AS, Werner G, Layer F. 2015. Linololid resistance in clinical isolates of *Staphylococcus epidermidis* from German hospitals and characterization of two *cfr*-carrying plasmids. J Antimicrob Chemother 70:1630–1638. https://doi.org/10.1093/jac/dkv025

42. Schouls LM, Veldman K, Brouwer MSM, Dierikx C, Witteveen S, van Santen-Verheuvel M, Hendrickx APA, Landman F, Hengeveld P, Wullings B, et al. 2022. *cfr* and *fexA* genes in methicillin-resistant *Staphylococcus aureus* from humans and livestock in the Netherlands. Commun Med 2:1–13. https://doi.org/10.1038/s43856-022-00200-w

43. Bukowski M, Lyzen R, Helbin WM, Bonar E, Szalewska-Palasz A, Wegrzyn G, Dubin G, Dubin A, Wladyka B. 2013. A regulatory role for *Staphylococcus aureus* toxin-antitoxin system PemlKsa. Nat Commun 4:2012. https://doi.org/10.1038/ncomms3012

44. Guerin F, Sassi M, Dejoies L, Zouari A, Schutz S, Potrel S, Auzou M, Collet A, Lecointe D, Auger G, Cattoir V. 2020. Molecular and functional analysis of the novel *cfr*(D) linololid resistance gene identified in *Enterococcus faecium*. J Antimicrob Chemother 75:1699–1703. https://doi.org/10.1093/jac/dkaa125

45. Caryl JA, O'Neill AJ. 2009. Complete nucleotide sequence of pG01, the prototype conjugative plasmid from the staphylococci. Plasmid 62:35–38. https://doi.org/10.1016/j.plasmid.2009.03.001

46. Urban D, Chevance A, Pasquereau B, Fourés F. 2024. Médicaments antimicrobiens chez l'animal. Surveillance des ventes et des utilisations. report. ANMV

47. Olsen JE, Christensen H, Aarestrup FM. 2006. Diversity and evolution of *bla*Z from *Staphylococcus aureus* and coagulase-negative staphylococci. J Antimicrob Chemother 57:450–460. https://doi.org/10.1093/jac/dki492

48. Wang H, Zhuang H, Ji S, Sun L, Zhao F, Wu D, Shen P, Jiang Y, Yu Y, Chen Y. 2021. Distribution of erm genes among MRSA isolates with resistance to clindamycin in a Chinese teaching hospital. Infect Genet Evol 96:105127. https://doi.org/10.1016/j.meegid.2021.105127

49. Al-Trad EI, Chew CH, Che Hamzah AM, Suhaili Z, Rahman NIA, Ismail S, Puah SM, Chua KH, Kwong SM, Yeo CC. 2023. The plasmidomic landscape of clinical methicillin-resistant *Staphylococcus aureus* isolates from Malaysia. Antibiotics (Basel) 12:733. https://doi.org/10.3390/antibiotics12040733

50. Geofrey MA, Sauli E, Kanje LE, Beti M, Shayo MJ, Kuchaka D, van Zwetselaar M, Wadugu B, Mmbaga B, Mkumbaye SI, Kumburu H, Sonda T. 2024. Genomic characterization of methicillin-resistant *Staphylococcus aureus* isolated from patients attending regional referral hospitals in Tanzania. BMC Med Genomics 17:211. https://doi.org/10.1186/s12920-024-01979-4

51. Gerken TJ, Roberts MC, Dykema P, Melly G, Lucas D, De Los Santos V, Gonzalez J, Butaye P, Wiegner TN. 2021. Environmental surveillance and characterization of antibiotic resistant *Staphylococcus aureus* at coastal beaches and rivers on the island of Hawai'i. Antibiotics (Basel) 10:980. https://doi.org/10.3390/antibiotics10080980

52. Neyaz L, Rajagopal N, Wells H, Fakhr MK. 2020. Molecular characterization of *Staphylococcus aureus* plasmids associated with strains isolated from various retail meats. Front Microbiol 11:223. https://doi.org/10.3389/fmicb.2020.000223

53. Tsergouli K, Karampatakis T, Kontopoulou K, Pappa S, Kampouridou P, Kallasidou G, Tsioka K, Zotou S, Farmaki E-E, Kotzamanidis C, Papa A. 2022. Spa diversity and genetic characterization of t127 methicillin-resistant *Staphylococcus aureus* in a tertiary Greek hospital. Acta Microbiol Immunol Hung 69:185–192. https://doi.org/10.1556/030.2022.01825

54. Projan SJ, Kornblum J, Moghazeh SL, Edelman I, Gennaro ML, Novick RP. 1985. Comparative sequence and functional analysis of pT181 and pC221, cognate plasmid replicons from *Staphylococcus aureus*. Mol Gen Genet 199:452–464. https://doi.org/10.1007/BF00330758

55. Wardal E, Gawryszewska I, Hryniewicz W, Sadowy E. 2013. Abundance and diversity of plasmid-associated genes among clinical isolates of *Enterococcus faecalis*. Plasmid 70:329–342. https://doi.org/10.1016/j.plasmid.2013.07.003

56. Fišarová L, Botka T, Du X, Mašlaňová I, Bárdy P, Pantůček R, Benešík M, Roudnický P, Winstel V, Larsen J, Rosenstein R, Peschel A, Doškář J. 2021. *Staphylococcus epidermidis* phages transduce antimicrobial resistance plasmids and mobilize chromosomal islands. mSphere 6:e00223-21. https://doi.org/10.1128/mSphere.00223-21

57. Humphrey S, San Millán Á, Toll-Riera M, Connolly J, Flor-Duro A, Chen J, Ubeda C, MacLean RC, Penadés JR. 2021. Staphylococcal phages and pathogenicity islands drive plasmid evolution. Nat Commun 12:5845. https://doi.org/10.1038/s41467-021-26101-5

58. Iordăneșcu S. 1977. Relationships between cotransducible plasmids in *Staphylococcus aureus*. J Bacteriol 129:71–75. https://doi.org/10.1128/jb.129.1.71-75.1977

59. Khan SA, Novick RP. 1983. Complete nucleotide sequence of pT181, a tetracycline-resistance plasmid from *Staphylococcus aureus*. Plasmid 10:251–259. https://doi.org/10.1016/0147-619x(83)90039-2

60. Wang PZ, Projan SJ, Henriquez V, Novick RP. 1993. Origin recognition specificity in pT181 plasmids is determined by a functionally asymmetric palindromic DNA element. EMBO J 12:45–52. https://doi.org/10.1002/j.1460-2075.1993.tb05630.x

61. Projan SJ, Moghazeh S, Novick RP. 1988. Nucleotide sequence of pS194, a streptomycin-resistance plasmid from *Staphylococcus aureus*. Nucleic Acids Res 16:2179–2187. https://doi.org/10.1093/nar/16.5.2179

62. Jain C, Rodriguez-R LM, Phillippy AM, Konstantinidis KT, Aluru S. 2018. High throughput ANI analysis of 90K prokaryotic genomes reveals clear species boundaries. Nat Commun 9:5114. https://doi.org/10.1038/s4146-017-0807-9
63. Partridge SR. 2011. Analysis of antibiotic resistance regions in Gram-negative bacteria. FEMS Microbiol Rev 35:820-855. https://doi.org/10.111/j.1574-6976.2011.00277.x

64. Redondo-Salvo S, Fernández-López R, Ruiz R, Vielva L, de Toro M, Rocha EPC, Garcillán-Barcia MP, de la Cruz F. 2020. Pathways for horizontal gene transfer in bacteria revealed by a global map of their plasmids. Nat Commun 11:3602. https://doi.org/10.1038/s41467-020-17278-2

65. Kennedy AD, Porcella SF, Martens C, Whitney AR, Braughton KR, Chen L, Craig CT, Tenover FC, Kreiswirth BN, Musser JM, DeLeo FR. 2010. Complete nucleotide sequence analysis of plasmids in strains of *Staphylococcus aureus* clone USA300 reveals a high level of identity among isolates with closely related core genome sequences. J Clin Microbiol 48:4504-4511. https://doi.org/10.1128/JCM.01050-10

66. Gillespie MT, May JW, Skurray RA. 1985. Antibiotic resistance in *Staphylococcus aureus* isolated at an Australian hospital between 1946 and 1981. J Med Microbiol 19:137-147. https://doi.org/10.1099/00222615-19-2-137

67. Bayles KW, Iandolo JJ. 1989. Genetic and molecular analyses of the gene encoding staphylococcal enterotoxin D. J Bacteriol 171:4799-4806. https://doi.org/10.1128/jb.171.9.4799-4806.1989

68. Zhang S, Iandolo JJ, Stewart GC. 1998. The enterotoxin D plasmid of *Staphylococcus aureus* encodes a second enterotoxin determinant (sej). FEMS Microbiol Lett 168:227-233. https://doi.org/10.1111/j.1574-6968.1998.tb13278.x

69. O'Brien FG, Ramsay JP, Monecke S, Coombs GW, Robinson OJ, Htet Z, Alshaikh FAM, Grubb WB. 2015. *Staphylococcus aureus* plasmids without mobilization genes are mobilized by a novel conjugative plasmid from

community isolates. J Antimicrob Chemother 70:649-652. https://doi.org/10.1093/jac/dku454

70. Monk IR, Foster TJ. 2012. Genetic manipulation of staphylococci-breaking through the barrier. Front Cell Infect Microbiol 2:49. https://doi.org/10.3389/fcimb.2012.00049

71. Haim MS, Zaheer R, Bharat A, Di Gregorio S, Di Conza J, Galanternik L, Lubovich S, Golding GR, Graham MR, Van Domselaar G, Cardona ST, Mollerach M. 2021. Comparative genomics of ST5 and ST30 methicillin-resistant *Staphylococcus aureus* sequential isolates recovered from paediatric patients with cystic fibrosis. Microb Genom 7:mgen000510. https://doi.org/10.1099/mgen.0.000510

72. Lee CA, Thomas J, Grossman AD. 2012. The *Bacillus subtilis* conjugative transposon ICEB51 mobilizes plasmids lacking dedicated mobilization functions. J Bacteriol 194:3165-3172. https://doi.org/10.1128/JB.00301-12

73. Soler N, Robert E, Chauvot de Beauchêne I, Monteiro P, Libante V, Maigret B, Staub J, Ritchie DW, Guédon G, Payot S, Devignes M-D, Leblond-Bourget N. 2019. Characterization of a relaxase belonging to the MOB~T~ family, a widespread family in firmicutes mediating the transfer of ICEs. Mob DNA 10:18. https://doi.org/10.1186/s13100-019-0160-9

74. Cafini F, Nguyen LTT, Higashide M, Román F, Prieto J, Morikawa K. 2016. Horizontal gene transmission of the *cfr* gene to MRSA and *Enterococcus*: role of *Staphylococcus epidermidis* as a reservoir and alternative pathway for the spread of linezolid resistance. J Antimicrob Chemother 71:587-592. https://doi.org/10.1093/jac/dkv391