Abstract
Staphylococci can cause a wide array of infections that can be life threatening. These infections become more deadly when the isolates are antibiotic resistant and thus harder to treat. Many resistance determinants are plasmid-mediated; however, staphylococcal plasmids have not yet been fully characterized. In particular, plasmids and their contributions to antibiotic resistance have not been investigated within the Arab states, where antibiotic use is not universally regulated. Here, we characterized the putative plasmid content among 56 Staphylococcus aureus and 10 Staphylococcus haemolyticus clinical isolates from Alexandria, Egypt. Putative plasmid sequences were detected in over half of our collection. In total, we identified 72 putative plasmid sequences in 27 S. aureus and 1 S. haemolyticus isolates. While these isolates typically carried one or two plasmids, we identified one isolate—S. aureus AA53—with 11 putative plasmids. The plasmid sequences most frequently encoded a Rep_1, RepL, or PriCT_1 type replication protein. As expected, antibiotic resistance genes were widespread among the identified plasmid sequences. Related plasmids were identified amongst our clinical isolates; homologous plasmids present in multiple isolates clustered into 11 groups based upon sequence similarity. Plasmids from the same cluster often shared antibiotic resistance genes, including blaZ, which is associated with β-lactam resistance. Our analyses suggest that plasmids are a key factor in the pathology and epidemiology of S. aureus in Egypt. A better characterization of plasmids and the role they contribute to the success of Staphylococci as pathogens will guide the design of effective control strategies to limit their spread.
Introduction
Staphylococci are commensal bacteria in humans and animals with the ability to cause serious infections. Prime among the Staphylococci is Staphylococcus aureus with pathologies that range from superficial skin lesions to bacteremia, pneumonia, and endocarditis, among others (). S. aureus is a major challenge in both the hospital and community settings, with treatment made more difficult given its ability to rapidly acquire antibiotic resistance (; ; ). The archetypal and most common antibiotic resistant examples are the Methicillin-resistant S. aureus (MRSA) strains (), which are usually resistant to most β-lactam antibiotics (). Strains of Coagulase-Negative Staphylococci (CoNS), including Staphylococcus epidermidis and Staphylococcus haemolyticus, also have acquired antibiotic resistance (). For example, over the past 40 years, Methicillin-resistant S. haemolyticus (MRSH) strains have become more prevalent (; ; ; ; ). To treat MRSA and MRSH, clinicians must resort to different classes of antibiotics, but further acquisition of antibiotic resistance by Methicillin-resistant strains is a major problem (). In Egypt, antibiotic use is not regulated, and most antimicrobial agents are available without the need for a prescription. This unregulated availability increases selective pressure that leads to development of antibiotic resistance. Thus, antibiotic therapy is often inadequate, a problem compounded by the use of the wrong antibiotic and both inappropriate dosage and duration. Consequently, Egypt has seen alarming increases in antibiotic resistance and MRSA/MRSH prevalence (; ; ; ).
Horizontal gene transfer is a major source for acquisition of antibiotic resistance genes in Staphylococci, facilitated by phage transduction and plasmid conjugation (). Plasmids can be reservoirs and vectors for antibiotic resistance, virulence, and fitness genes, even in the absence of selective pressure (). Horizontal transfer of plasmids can occur not only between Staphylococcal strains of the same species but also between species, with CoNS believed to act as reservoirs for antibiotic resistance genes (). Moreover, recent evidence suggests that plasmids also can be exchanged between Staphylococci and the distantly related Enterococci (). Data on the diversity of plasmid content of Staphylococci is limited (; ; ; ), but it is known that Staphylococcal plasmids can range in size from 1 to 60 kbp and that they can carry resistance genes (or gene clusters) (). Clinical isolates often contain one or more multi-resistance plasmids (). Thus, plasmids are key to the epidemiology of Staphylococci, contributing to the rapid spread of antibiotic resistance, and play a major role in determining clinical outcomes (; ). While over 1,000 S. aureus plasmid sequences and 15 S. haemolyticus plasmid sequences are publicly available in resources such as PLSDB (), just 6 S. aureus and no S. haemolyticus plasmids are annotated as collected from the Arab states. This knowledge gap motivated our investigation of plasmids from this region, and we hypothesize that plasmids in Egyptian Staphylococcal strains encode for antibiotic resistance.
Large plasmids in bacteria can be classified by their incompatibility groups, as plasmids with similar replication initiation (Rep) proteins cannot be stably maintained in the absence of selective pressure (; ). To date, 7 distinct types of Rep proteins and 15 incompatibility groups have been identified for Staphylococcal plasmids (). Accordingly, Staphylococcal plasmids can be classified using the sequence of their Rep genes (). Antibiotic resistance genes may be linked to Rep families; for example, tetracycline resistance has been associated with the rep7 gene and chloramphenicol resistance with the rep7, rep13, and rep27 genes (). Currently, there are limited data on the Rep families present in circulating Staphylococcal isolates, including those from Egypt. Knowledge of the plasmid Rep families present in these clinical isolates could give us insight into Staphylococcal plasmid biology and thus help to better understand Methicillin-resistant strains.
In this study, we utilized genomics to (i) catalog the Rep family classification in the plasmids of Egyptian clinical isolates of S. aureus and S. haemolyticus, (ii) profile the antibiotic resistance genes present in these plasmids, and (iii) assess the prevalence of plasmid-borne antibiotic resistance in Egyptian MRSA isolates. We report on 72 plasmids present in 27 S. aureus and 1 S. haemolyticus isolates from Egypt. We also report on the antibiotic resistance genes in these plasmids, including those from MRSA and MRSH isolates. We furthermore provide bioinformatic and genetic context for these plasmid sequences. We found that the plasmids in Egyptian isolates of S. aureus and S. haemolyticus are diverse and include antibiotic resistance genes for antibiotics used clinically, indicating the potential for plasmids to act as gene reservoirs and thus facilitate the rapid horizontal transfer of antibiotic resistance in Egypt.
Materials and Methods
Identifying Plasmid Sequences
Raw sequencing reads for the 56 S. aureus and 10 S. haemolyticus genomes were previously sequenced, annotated, and typed by our group (). Here, the raw reads were assembled using plasmidSPAdes (v3.9.0) (). Plasmid sequences with a length greater than 1,000 bp were selected for further analysis. Each putative plasmid sequence was queried via megablast against the nr/nt database (April 2020). The best hit to a complete plasmid sequence record for each sequence was identified and recorded. PlasmidFinder (v2.1) hosted by the Center for Genomic Epidemiology was queried to identify replicon sequences (). BLAST and PlasmidFinder results were considered together to distinguish putative plasmids from other sequences identified by plasmidSPAdes (e.g., bacteriophage). To verify even coverage, reads were mapped to plasmid sequences using the Bowtie2 (v2.3.2) () plug-in through Geneious Prime v2019.1.1 (Biomatters Ltd., Auckland, NZ). To estimate plasmid copy numbers, genome and plasmid coverages were calculated using BBMap v38.471.
Annotating Plasmid Sequences
Putative plasmid sequences were annotated using PATRIC v3.6.3 (; ). Putative plasmid sequences also were examined for antibiotic resistance using ResFinder (v3.2) (), using default parameters (%ID threshold 90%, selected minimum length 60%).
Clustering Similar Plasmids
Putative plasmid sequences were clustered using the USEARCH algorithm (v11.0.667) () with the “cluster_fast” option for id = 0.8 to both strands. Each cluster was aligned using MAFFT v7.388 () and Mauve () through Geneious Prime and manually inspected. BLAST hits were referenced to verify that plasmids within the same cluster exhibited similarity to related (if not the same) plasmid sequence in the nr/nt database. Plasmid alignments were generated using BLAST Ring Image Generator v0.95 ().
Phylogenomic Tree of S. aureus Strains
Anvi’o () was used to identify the single copy number genes amongst the S. aureus genome sequences. A phylogenomic tree was created using concatenated amino acid sequences for these single copy number genes. These concatenated sequences were aligned using MAFFT v7.388 () through Geneious Prime (Biomatters Ltd., Auckland, New Zealand). The tree was derived using FastTree v2 () through Geneious Prime and visualized by iTOL v5.6.1 ().
Results
Upon the examination of 56 S. aureus and 10 S. haemolyticus genomes of Egyptian clinical isolates, we identified 72 putative plasmid sequences in 28 strains: 27 strains of S. aureus and 1 strain of S. haemolyticus. From our previous work (), we know that 24 of these isolates are mecA-positive and represent 11 different MLSTs and 2 unknown MLSTs. While 10 of the strains harbored just a single putative plasmid sequence, one strain—AA53—contained 11 (Table 1). The copy number of the plasmid sequences ranged from ∼1 to >300, determined by comparing the coverage of each individual plasmid sequence relative to the average coverage for the assembly (Supplementary Table 1). 50 of the 72 putative plasmid sequences contained replicon proteins identified by PlasmidFinder. Table 2 summarizes these results. The majority of the strains (85%) contained plasmid sequences associated with rep16 (n = 10), rep20 (n = 9), or rep21 (n = 9). The replicons rep22, rep24a, and rep39 were each unique to a single isolate. As expected because of incompatibility, no strain was found to have more than one plasmid for a given replicon protein. Putative plasmid sequences containing a replicon protein frequently showed sequence identity to a Staphylococcus plasmid sequence (Supplementary Table 2).
TABLE 1
| Strain | MLST | Methicillin resistance | GenBank WGS master accession | # Plasmids |
| S. aureus AA1 | ST-1 | mecA +ve | JAEOUR000000000 | 2 |
| S. aureus AA2 | ST-80 | mecA +ve | JAEOUZ000000000 | 5 |
| S. aureus AA3 | ST-80 | mecA +ve | JAEOVE000000000 | 5 |
| S. aureus AA4 | ST-80 | mecA +ve | JAEOVM000000000 | 5 |
| S. aureus AA5 | ST-22 | mecA +ve | JAEOVQ000000000 | 1 |
| S. aureus AA6 | ST-97 | mecA +ve | JAEOVX000000000 | 1 |
| S. aureus AA8 | ST-97 | mecA +ve | JAEOWM000000000 | 1 |
| S. aureus AA13 | ST-239 | mecA +ve | JAEOUV000000000 | 1 |
| S. aureus AA17 | ST-152 | mecA +ve | JAEOUX000000000 | 3 |
| S. aureus AA29 | ST-239 | mecA +ve | JAEOVD000000000 | 2 |
| S. aureus AA30 | ST-6 | mecA +ve | JAEOVF000000000 | 1 |
| S. aureus AA32 | Unknowna | mecA +ve | JAEOVH000000000 | 1 |
| S. aureus AA35 | Unknown | mecA +ve | JAEOVJ000000000 | 4 |
| S. aureus AA36 | ST-97 | mecA +ve | JAEOVK000000000 | 1 |
| S. aureus AA41 | ST-1482 | mecA -ve | JAEOVN000000000 | 2 |
| S. haemolyticus AA42 | ST-30 | mecA +ve | JAEOWY000000000 | 2 |
| S. aureus AA45 | ST-80 | mecA +ve | JAEOVO000000000 | 4 |
| S. aureus AA51 | ST-1 | mecA +ve | JAEOVR000000000 | 2 |
| S. aureus AA53 | ST-88 | mecA +ve | JAEOVT000000000 | 11 |
| S. aureus AA57 | ST-239 | mecA +ve | JAEOVV000000000 | 2 |
| S. aureus AA59 | ST-1 | mecA -ve | JAEOVW000000000 | 1 |
| S. aureus AA63 | ST-239 | mecA +ve | JAEOWB000000000 | 3 |
| S. aureus AA64 | ST-239 | mecA +ve | JAEOWC000000000 | 3 |
| S. aureus AA65 | ST-1 | mecA -ve | JAEOWD000000000 | 3 |
| S. aureus AA68 | ST-1 | mecA -ve | JAEOWF000000000 | 1 |
| S. aureus AA69 | ST-1 | mecA +ve | JAEOWG000000000 | 2 |
| S. aureus AA70 | ST-5 | mecA +ve | JAEOWH000000000 | 1 |
| S. aureus AA78 | ST-1 | mecA +ve | JAEOWK000000000 | 2 |
List of Egyptian S. aureus and S. haemolyticus strains harboring plasmids.
aIncomplete alignment.
TABLE 2
| Strain | # Replicons | Replicon protein | ||||||||
| 5aa | 7ab | 10c | 16d | 20e | 21f | 22f | 24ae | 39e | ||
| S. aureus AA1 | 2 | + | + | |||||||
| S. aureus AA2 | 1 | + | ||||||||
| S. aureus AA3 | 1 | + | ||||||||
| S. aureus AA4 | 2 | + | + | |||||||
| S. aureus AA5 | 1 | + | ||||||||
| S. aureus AA6 | 1 | + | ||||||||
| S. aureus AA8 | 1 | + | ||||||||
| S. aureus AA13 | 1 | + | ||||||||
| S. aureus AA17 | 3 | + | + | + | ||||||
| S. aureus AA29 | 2 | + | + | |||||||
| S. aureus AA30 | 1 | + | ||||||||
| S. aureus AA32 | 1 | + | ||||||||
| S. aureus AA35 | 4 | + | + | + | + | |||||
| S. aureus AA36 | 1 | + | ||||||||
| S. aureus AA41 | 1 | + | ||||||||
| S. aureus AA45 | 1 | + | ||||||||
| S. aureus AA51 | 1 | + | ||||||||
| S. aureus AA53 | 6 | + | + | + | + | + | + | |||
| S. aureus AA57 | 2 | + | + | |||||||
| S. aureus AA59 | 1 | + | ||||||||
| S. aureus AA63 | 3 | + | + | + | ||||||
| S. aureus AA64 | 3 | + | + | + | ||||||
| S. aureus AA65 | 3 | + | + | + | ||||||
| S. aureus AA68 | 1 | + | ||||||||
| S. aureus AA69 | 1 | + | ||||||||
| S. aureus AA70 | 1 | + | ||||||||
| S. aureus AA78 | 2 | + | + | |||||||
| S. haemolyticus AA42 | 2 | + | + | |||||||
Number and type of replicons identified by PlasmidFinder for plasmid sequences in the Egyptian S. aureus strains and S. haemolyticus strain.
“+” indicates the presence of a plasmid. Replication Protein Types: aRep_3, bRep_trans, cRepL, dPriCT_1, eRepA_N, fRep_1.
To ascertain if the 22 putative plasmid sequences that did not contain a replicon protein likely represented a plasmid, we referred to the BLAST hits (Supplementary Table 2). These all showed greatest sequence similarity (in some cases, even 100% sequence identity) to a previously reported Staphylococcal plasmid sequence, often one that was significantly longer than our putative plasmid sequence. Of note are the putative plasmid sequences within S. aureus strains AA2, AA3, AA4, and AA45. Each of these strains contained three or four sequences that showed sequence similarity to the same S. aureus plasmid sequence, pGR2A (Accession no. CP010403). For each of these strains, one sequence had a recognizable rep20 replicon. We downloaded the pGR2A record and mapped the raw reads for strains AA2, AA3, AA4, and AA45 to the pGR2A sequence. While the reads mapped to the plasmid sequence, there were two or three regions of low coverage (Figure 1), suggesting that either the plasmid within these strains: (1) has undergone reassortment or expansion, differentiating it from pGR2A; (2) has split into separate extrachromosomal DNAs; (3) contains a region that was not sequenced by the Illumina short-read technology; and/or integrated into the genome. The pGR2A-like plasmid sequences accounted for 14 of the 22 putative plasmid sequences without a recognized replicon protein. Thus, there are just eight putative plasmid sequences in our data set for which a replicon protein was not identified.
FIGURE 1
To identify homologous plasmids found in more than one isolate, all 72 putative plasmid sequences were clustered. In total, we identified 11 clusters while 16 plasmids did not show any sequence similarity to another plasmid in our data set. These 16 unclustered plasmids, however, did exhibit sequence similarity to plasmid sequences deposited in GenBank (Supplementary Table 2). This included the two plasmids within the S. haemolyticus A42 strain. Table 3 lists details about these clusters. The largest clusters, plasmid Cluster C and plasmid Cluster J, contained a plasmid sequence found in 10 and 5 different strains, respectively. In some cases, all members of the cluster showed greatest homology to the same GenBank plasmid record, for instance Cluster A (n = 4), Cluster B (n = 4), Cluster D (n = 2), Cluster E (n = 2), and Cluster K (n = 2). In other cases, members of the cluster exhibited greatest sequence similarity to different GenBank plasmid records, reflecting the differences between the homologous plasmids, as well as hinting at diversity in the plasmid group amongst strains of S. aureus. Figure 2 shows that instances of gene acquisition/loss can be detected between homologous plasmids in a cluster. For some clusters, the plasmids are carried by closely related S. aureus strains, while other clusters include plasmids found in strains that span the genomic diversity of the Egyptian isolates (Figure 3).
TABLE 3
| Cluster ID (replicon) | Cluster details | Closest relative (by BLAST homology) | |||||
| Strains | Avg. length (bp) | %GC | % Pairwise identity | Description | Accession # | Length (bp) | |
| A (rep7a) | AA29, AA35, AA57, AA64 | 3,909 | 29.3 | 99.5 | Unnamed3 | CP033985 | 3,785 |
| B (rep20) | AA2, AA3, AA4, AA45 | 27,938 | 29.1 | 100 | pGR2A | CP010403 | 28,895 |
| C (rep16) | AA17 | 20,970 | 28.7 | 80.7 | pPS00087.1A.1 | CP022721 | 20,730 |
| AA30 | pl1_M2024 | CP047022 | 19,760 | ||||
| AA41 | pWBG757 | GQ900397 | 20,730 | ||||
| AA1, AA59, AA65, AA68, AA69, AA78 | pWBG750 | GQ900392 | 20,653 | ||||
| AA51 | p18809-P04 | CP002146 | 28,404 | ||||
| D (rep5a) | AA5, AA32 | 34,114 | 30.3 | 100 | pIT4-R | CP028471 | 34,104 |
| E (rep20) | AA6, AA8, AA53 | 16,610 | 28.1 | 99.2 | pM121 | CP007671 | 20,409 |
| F (rep20) | AA35 | 21,196 | 28.4 | 80 | Unnamed1 | CP030394 | 25,071 |
| AA36 | Unnamed1 | CP030413 | 20,403 | ||||
| G (rep21) | AA35 | 3,112 | 29 | 99.5 | pV605 | CP013960 | 43,115 |
| AA4 | SAP104A | GQ900450 | 3,011 | ||||
| AA57, AA64 | pV605 | CP013960 | 43,115 | ||||
| H (rep10) | AA1 | 2,560 | 30.8 | 96.8 | Unnamed1 | CP029648 | 2,473 |
| AA17 | pSau-2716Lar | MH423311 | 3,267 | ||||
| AA35 | Unnamed | CP047816 | 2,473 | ||||
| AA63 | Unnamed2 | CP047819 | 2,477 | ||||
| I (rep21) | AA53 | 3,123 | 29.1 | 98.8 | 3 | LR130517 | 3,011 |
| AA63 | Unnamed1 | CP034006 | 3,011 | ||||
| J (rep10) | AA13 | 2,556 | 31 | 98.9 | pUSA05-1-SUR24 | CP014446 | 2,415 |
| AA29 | pSR265 | MN251859 | 2,413 | ||||
| AA53 | p19321-P01 | CP002148 | 2,473 | ||||
| AA64 | SAP078B | GQ900431 | 2,415 | ||||
| AA65 | pSR241 | MN251858 | 2,402 | ||||
| K (rep7a) | AA53, AA63 | 3,901 | 29.6 | 98.8 | pOC160-2 | LC012933 | 3,788 |
Clusters of plasmids based on nucleotide sequence similarity within the Egyptian S. aureus strains.
FIGURE 2
FIGURE 3
We annotated the plasmid sequences using the web service PATRIC (Supplementary Table 3). 38.6% of the coding regions within the 72 plasmid sequences are annotated as hypothetical proteins. PATRIC did identify resistances to metals, e.g., cadmium, copper, cobalt/zinc/cadmium, lead, mercury, and zinc, and antibiotics. Furthermore, the plasmid from S. aureus AA70 also encodes for exotoxin proteins. To further investigate antibiotic resistance genes, each plasmid sequence was examined with the tool ResFinder. Table 4 summarizes the antibiotic resistance determinants predicted by ResFinder, and full results can be found in Supplementary Table 4. 20 of the 27 (74%) S. aureus strains included a putative plasmid encoding blaZ, which is associated with β-lactam resistance. These included all four of the mecA-negative S. aureus strains (AA41, AA59, AA65, and AA68). S. aureus AA69 and AA41 were the only strains carrying a plasmid encoding aminoglycoside resistance. The one plasmid-carrying S. haemolyticus strain (AA42) did not encode blaZ; however, it carried plasmids that encode for lincosamide resistance, macrolide resistance, and streptogramin B resistance. This strain is the only one identified in this study that carries a plasmid that encodes for the lnu(A) gene, which is associated with lincosamide resistance.
TABLE 4
| Predicted phenotype | Resistance gene | Strains |
| Aminoglycoside resistance | aac(6′)-aph(2”) | AA69 |
| ant(6)-Ia | AA41 | |
| aph(3′)-III | AA41 | |
| Beta-lactam resistance | blaZ | AA1, AA2, AA3, AA4, AA6, AA8, AA17, AA30, AA35, AA36, AA41, AA45, AA51, AA53, AA59, AA65, AA68, AA69, AA70, AA78 |
| Fusidic acid resistance | fusB | AA2, AA3, AA4, AA45 |
| Lincosamide resistance | lnu(A) | AA42 |
| Macrolide resistance | erm(C) | AA1, AA1, AA13, AA17, AA29, AA35, AA53, AA63, AA63, AA64, AA65 |
| mph(C) | AA41, AA42 | |
| Macrolide, Lincosamide, and Streptogramin B resistance | msr(A) | AA41, AA42 |
| Phenicol resistance | cat(pC221) | AA29, AA35, AA35, AA35, AA53, AA57, AA63, AA64, AA64, AA64 |
| Tetracycline resistance | tet(K) | AA2, AA3, AA4, AA17, AA45 |
| tet(L) | AA1, AA51, AA59, AA65, AA68, AA69, AA78 |
Resistance genes carried on Egyptian S. aureus and S. haemolyticus plasmids.
Whereas 48.6% of the plasmids carried a single antibiotic resistance gene, others carried multiple antibiotic resistance genes. We therefore examined the plasmid sequences for resistance co-occurrence. Strains harboring plasmids encoding for fusidic acid resistance (fusB) (n = 4) also harbored distinct plasmids that encoded for tetracycline resistance [tet(K)] (n = 5) (r = 0.9). Furthermore, all four strains with plasmids carrying both fusB and tet(K) were MLST-80. Plasmids encoding for tet(L) were only harbored by strains of MLST-1. However, this association between tet(L) and MLST-1 was not observed more broadly; publicly available MLST-1 genomes (Supplementary Table 5) were queried for the tet(L) sequence. Only 8 of the 112 strains queried also coded for tet(L). Interestingly, six of these were Egyptian MRSA isolates (Supplementary Table 5).
Discussion
Plasmids are key reservoirs for genetic content in Staphylococci and allow the rapid propagation of antibiotic resistance (; ). Here, we screened the genomes of 56 clinical isolates of S. aureus and 10 clinical isolates of S. haemolyticus previously sequenced by our group and found 72 putative plasmid sequences (Table 1). We identified replicon proteins in 50 of the 72 plasmids. Within the S. aureus strains, eight different replicon proteins were detected, representative of six of the seven Staphylococcal replicon protein types, and the S. haemolyticus isolate included a plasmid with an additional replicon protein not found in the S. aureus strains, rep39 (Table 2). This particular plasmid (rep39) was homologous to pLNU4, a plasmid found in an S. chromogenes isolate that was associated with bovine mastitis, and can be acquired by S. aureus via transformation (). For strains containing more than 1 plasmid, 11 different combinations of the replicon proteins were found. Strains containing plasmids with rep7a and rep10 were most common (n = 6); whereas one strain, S. aureus AA13, contained just a single rep10 plasmid, no strains contained a single rep7a plasmid. Rep sites were not identified in eight of the putative plasmids, although their sequence reads had homology to plasmid sequences. These eight could represent partial plasmid sequences or may rely on vertical transfer of plasmids, as is the case of smaller plasmids with high copy number (). However, based on both plasmid copy number (Supplementary Table 1) and BLAST hits (Supplementary Table 2), we believe that these eight sequences are in fact partial plasmid sequences.
S. aureus strains AA2, AA3, AA4, and AA45 had putative plasmid sequences that showed sequence similarity to the S. aureus plasmid pGR2A (Accession no. CP010403). The pGR2A plasmid codes a unique bla system that can suppress oxacillin resistance in MRSA strains (). The non-functional BlaR1 from pGR2A does not cleave the mecA repressor BlaI; thus, mecA is expressed inefficiently, resulting in a phenotype of penicillin resistance, but oxacillin susceptibility (). It must be noted that all four strains with plasmid sequences homologous to pGR2A were mecA-positive strains (Tables 1, 3). In the clinical setting, oxacillin-susceptible mecA-positive S. aureus may pose a challenge due to possible misdiagnosis as a typical MRSA (; ).
In addition to classifying plasmids by their Rep genes, we analyzed plasmid relatedness by clustering putative plasmid-coding regions based on sequence homology (Table 3). We identified 11 plasmid clusters among the S. aureus strains; the two S. haemolyticus plasmid sequences did not resemble any of the S. aureus plasmids. Members of a cluster could have homology to different GenBank plasmid records (C, J), or all could share the same GenBank plasmid record (A, B, D, E, K). The former likely indicates homology of non-identical plasmids and the acquisition/loss of content as indicated by cluster C (Figure 2), while the latter points to multiple S. aureus isolates possessing the same plasmid, potentially due to a shared selective pressure (; ). Our results found that Rep genes are not exclusive to a single cluster, concurring with prior findings of frequent gene exchange amongst Staphylococcal plasmids (see review; ). Strains harboring plasmids from the same cluster often shared antibiotic resistance genes, e.g., macrolide resistance in Cluster A plasmids, fusidic acid resistance in Cluster B, and beta-lactam resistance in Cluster C. Thus, while replicon proteins can provide insight into plasmid incompatibilities, they have limited epidemiological information compared to the cluster to which a plasmid belongs.
Of the 28 isolates in which we identified putative plasmid sequences, 24 (86%) were mecA-positive in addition to widespread presence of resistance genes in their plasmid sequences (Table 4). In the clinical setting, the main classes of drugs used against Staphylococcus are those that target the cell wall (e.g., β-lactam antibiotics, glycopeptides), ribosome (e.g., tetracyclines, aminoglycosides, macrolides), and nucleic acid biosynthesis (fluoroquinolones and sulfamethoxazole/trimethoprim) (). Resistance to these antibiotic classes is a critical concern, especially in the acquisition of multiple resistances (; ). There have been reports of over-the-counter antibiotic use in Egypt that could select for resistance to key antibiotics, e.g., penicillins and cephalosporins (; ). The most common antibiotic resistances in the plasmid sequences were those to macrolides, tetracyclines, fusidic acid, and aminoglycosides, and 20 of our Staphylococcus putative plasmid sequences had genes linked to β-lactam antibiotics resistance (Table 4). While tetracycline resistance can be plasmid- or chromosome-mediated, macrolide resistance is almost exclusively encoded by mobile elements in Staphylococci. Nine different Egyptian S. aureus strains carried plasmids with the 23S rRNA methylase erm(C) (Table 4). Previous studies have found that erm(C) is the most abundant gene conferring resistance to macrolides, as well as to lincosamides and streptogramin B compounds (see review; ). mph(C), found in S. aureus AA41 and S. haemolyticus AA42, mediates resistance to only macrolides, whereas lnu(A), found only in S. haemolyticus AA42, encodes only lincosamide resistance. In addition, our annotation indicated resistance genes for chloramphenicol, clindamycin, and sulfamethoxazole/trimethoprim (Supplementary Table 3). The presence of these antibiotic resistance genes can greatly limit the antibiotic options to treat these clinical isolates (; ; ). Of utmost concern is that multiple resistances could be disseminated via plasmids, especially in MRSA ().
We did not identify genes involved in conjugation in our isolates (Supplementary Table 3). It has been noted that S. aureus may rely on vertical transmission of its plasmids, while conjugation genes are relatively rare (). The distribution of plasmids belonging to the same cluster across the Egyptian S. aureus isolates suggests both vertical and horizontal transfer (Figure 3). Horizontal acquisition of the plasmid in our isolates could occur via natural competence or by conjugation systems not yet characterized and therefore not present in our annotation references (). It must be highlighted that many of the plasmid sequences presented in these Egyptian isolates were homologous to GenBank entries for plasmids that have not been well-characterized (Table 3); from 37 plasmids homologous to our sequences, only seven had a publication associated with their entry. This underscores the general lack of S. aureus plasmid knowledge. Furthermore, 38.6% of the genes annotated within our plasmids are characterized as hypothetical proteins. Their function and/or the possible genetic advantage that these genes provide the bacterial strains remains an open question. The presence of genes for multiple antibiotic resistance in uncharacterized plasmids indicates not just the relevance of plasmids in the clinical treatment of Egyptian Staphylococcus, but also points to the depth of knowledge still to be discovered in this research area (; ).
In this study, we utilized genomics to report on the plasmid content of Egyptian Staphylococci, which further adds to our global knowledge of Staphylococcal plasmids. Plasmids and their contributions to antibiotic resistance have not been investigated previously within the Arab states. Putative plasmid sequences were present in over half of isolates analyzed in this study, with widespread presence of antibiotic resistance genes. Since plasmids are reservoirs of and vectors for antibiotic resistance in Staphylococci, plasmids could be a key factor in the pathology and epidemiology of Staphylococci, including MRSA and MRSH (; ). Making this knowledge available allows policy makers and infection control officers the chance to design and implement well-informed infection control and antimicrobial stewardship plans. Further characterization of plasmids is necessary to better understand the infectious potential of Staphylococci and devise optimal strategies to curb their spread and virulence.
Statements
Data availability statement
Data analyzed as part of this study include the following GenBank WGS accession numbers: JAEOUR000000000, JAEOUZ000000000, JAEOVE000000000, JAEOVM000000000, JAEOVQ000000000, JAEOVX000000000, JAEOWM000000000, JAEOUV000000000, JAEOUX000000000, JAEOVD000000000, JAEOVF000000000, JAEOVH000000000, JAEOVJ000000000, JAEOVK000000000, JAEOVN000000000, JAEOWY000000000, JAEOVO000000000, JAEOVR000000000, JAEOVT000000000, JAEOVV000000000, JAEOVW000000000, JAEOWB000000000, JAEOWC000000000, JAEOWD000000000, JAEOWF000000000, JAEOWG000000000, JAEOWH000000000, and JAEOWK000000000.
Author contributions
CM, AJW, and AA designed the study. CRM and CP conducted the analyses. CRM, CM, CP, and AA wrote the initial draft. All authors contributed to the manuscript editing and approved the final version.
Funding
We acknowledge funding from NIH (R01 DK104718 awarded to AJW), NSF (1661357 awarded to CP), and USAID (GSP-T85 awarded to AA) and DFG (ZI 665/3-1 awarded to AA). The funders did not play a part in the design or conduct of the study.
Conflict of interest
AJW discloses membership on the Advisory Boards of Pathnostics and Urobiome Therapeutics. This membership did not play a part in the design of the study, analysis or interpretation of data, writing of the manuscript, or decision to publish the results. The remaining 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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2021.659116/full#supplementary-material
Supplementary Table 1Plasmid vs. genome sequence coverage.
Supplementary Table 2BLAST-based similarity of putative plasmid sequences to records in nr/nt database.
Supplementary Table 3Annotation of putative plasmid sequences.
Supplementary Table 4Antibiotic resistance genes identified for each putative plasmid sequence.
Supplementary Table 5Publicly available MLST-1 S. aureus genomes screened for tet(L).
Footnotes
References
1
AbouelfetouhA. (2017). The status of methicillin resistance among Egyptian Staphylococcus aureus isolates: an overview.Infect. Disord. Drug Targets1767–69. 10.2174/1871526516666160802111200
2
AbouelfetouhA.KassemM.NaguibM.El-NakeebM. (2017). Investigation and treatment of fusidic acid resistance among methicillin-resistant staphylococcal isolates from Egypt.Microb. Drug Resist.238–17. 10.1089/mdr.2015.0336
3
AcmanM.van DorpL.SantiniJ. M.BallouxF. (2020). Large-scale network analysis captures biological features of bacterial plasmids.Nat. Commun.11:2452. 10.1038/s41467-020-16282-w
4
AlikhanN.-F.PettyN. K.Ben ZakourN. L.BeatsonS. A. (2011). BLAST ring image generator (BRIG): simple prokaryote genome comparisons.BMC Genomics12:402. 10.1186/1471-2164-12-402
5
Al-TamimiM.Abu-RaidehJ.HimsawiN.KhasawnehA.HawamdehH. (2020). Methicillin and vancomycin resistance in coagulase-negative Staphylococci isolated from the nostrils of hospitalized patients.J. Infect. Dev. Ctries.1428–35. 10.3855/jidc.11025
6
AntipovD.HartwickN.ShenM.RaikoM.LapidusA.PevznerP. A. (2016). plasmidSPAdes: assembling plasmids from whole genome sequencing data.Bioinformatics323380–3387. 10.1093/bioinformatics/btw493
7
BeckerK.HeilmannC.PetersG. (2014). Coagulase-negative staphylococci.Clin. Microbiol. Rev.27870–926. 10.1128/CMR.00109-13
8
BorgM. A.de KrakerM.SciclunaE.van de Sande-BruinsmaN.TiemersmaE.MonenJ.et al (2007). Prevalence of methicillin-resistant Staphylococcus aureus (MRSA) in invasive isolates from southern and eastern Mediterranean countries.J. Antimicrob. Chemother.601310–1315. 10.1093/jac/dkm365
9
BouchamiO.AchourW.MekniM. A.RoloJ.HassenA. B. (2011). Antibiotic resistance and molecular characterization of clinical isolates of methicillin-resistant coagulase-negative staphylococci isolated from bacteremic patients in oncohematology.Folia Microbiol.56122–130. 10.1007/s12223-011-0017-1
10
BrettinT.DavisJ. J.DiszT.EdwardsR. A.GerdesS.OlsenG. J.et al (2015). RASTtk: a modular and extensible implementation of the RAST algorithm for building custom annotation pipelines and annotating batches of genomes.Sci. Rep.5:8365. 10.1038/srep08365
11
CarattoliA.ZankariE.García-FernándezA.Voldby LarsenM.LundO.VillaL.et al (2014). In silico detection and typing of plasmids using PlasmidFinder and plasmid multilocus sequence typing.Antimicrob. Agents Chemother.583895–3903. 10.1128/AAC.02412-14
12
DarlingA. C. E.MauB.BlattnerF. R.PernaN. T. (2004). Mauve: multiple alignment of conserved genomic sequence with rearrangements.Genome Res.141394–1403. 10.1101/gr.2289704
13
EdgarR. C. (2010). Search and clustering orders of magnitude faster than BLAST.Bioinformatics262460–2461. 10.1093/bioinformatics/btq461
14
ErenA. M.EsenÖC.QuinceC.VineisJ. H.MorrisonH. G.SoginM. L.et al (2015). Anvi’o: an advanced analysis and visualization platform for ’omics data.PeerJ3:e1319. 10.7717/peerj.1319
15
FalagasM. E.KarageorgopoulosD. E.LeptidisJ.KorbilaI. P. (2013). MRSA in Africa: filling the global map of antimicrobial resistance.PLoS One8:e68024. 10.1371/journal.pone.0068024
16
FerreiraR. B. R.NunesA. P. F.KokisV. M.KrepskyN.FonsecaL.deS.et al (2002). Simultaneous detection of the mecA and ileS-2 genes in coagulase-negative staphylococci isolated from Brazilian hospitals by multiplex PCR.Diagn. Microbiol. Infect. Dis.42205–212. 10.1016/s0732-8893(01)00345-5
17
FeßlerA. T.WangY.WuC.SchwarzS. (2018). Mobile macrolide resistance genes in staphylococci.Plasmid992–10. 10.1016/j.plasmid.2018.05.001
18
FosterT. J. (2017). Antibiotic resistance in Staphylococcus aureus. Current status and future prospects.FEMS Microbiol. Rev.41430–449. 10.1093/femsre/fux007
19
GalataV.FehlmannT.BackesC.KellerA. (2019). PLSDB: a resource of complete bacterial plasmids.Nucleic Acids Res.47D195–D202. 10.1093/nar/gky1050
20
GrundmannH.Aires-de-SousaM.BoyceJ.TiemersmaE. (2006). Emergence and resurgence of meticillin-resistant Staphylococcus aureus as a public-health threat.Lancet368874–885. 10.1016/S0140-6736(06)68853-3
21
Hill-CawthorneG. A.HudsonL. O.El GhanyM. F. A.PiepenburgO.NairM.DodgsonA.et al (2014). Recombinations in staphylococcal cassette chromosome mec elements compromise the molecular detection of methicillin resistance in Staphylococcus aureus.PLoS One9:e101419. 10.1371/journal.pone.0101419
22
HisatsuneJ.HirakawaH.YamaguchiT.FudabaY.OshimaK.HattoriM.et al (2013). Emergence of Staphylococcus aureus carrying multiple drug resistance genes on a plasmid encoding exfoliative toxin B.Antimicrob. Agents Chemother.576131–6140. 10.1128/AAC.01062-13
23
IkonomidisA.MichailG.VasdekiA.LabrouM.KaravasilisV.StathopoulosC.et al (2008). In vitro and in vivo evaluations of oxacillin efficiency against mecA-positive oxacillin-susceptible Staphylococcus aureus.Antimicrob. Agents Chemother.523905–3908. 10.1128/AAC.00653-08
24
JensenL. B.Garcia-MiguraL.ValenzuelaA. J. S.LøhrM.HasmanH.AarestrupF. M. (2010). A classification system for plasmids from enterococci and other Gram-positive bacteria.J. Microbiol. Methods8025–43. 10.1016/j.mimet.2009.10.012
25
KatohK.StandleyD. M. (2013). MAFFT multiple sequence alignment software version 7: improvements in performance and usability.Mol. Biol. Evol.30772–780. 10.1093/molbev/mst010
26
KwongS. M.RamsayJ. P.JensenS. O.FirthN. (2017). Replication of Staphylococcal Resistance Plasmids.Front. Microbiol.8:2279. 10.3389/fmicb.2017.02279
27
LangmeadB.SalzbergS. L. (2012). Fast gapped-read alignment with Bowtie 2.Nat. Methods9357–359. 10.1038/nmeth.1923
28
LetunicI.BorkP. (2019). Interactive Tree Of Life (iTOL) v4: recent updates and new developments.Nucleic Acids Res.47W256–W259. 10.1093/nar/gkz239
29
LüthjeP.von Köckritz-BlickwedeM.SchwarzS. (2007). Identification and characterization of nine novel types of small staphylococcal plasmids carrying the lincosamide nucleotidyltransferase gene lnu(A).J. Antimicrob. Chemother.59600–606. 10.1093/jac/dkm008
30
MashalyG. E.-S.El-MahdyR. H. (2017). Vancomycin heteroresistance in coagulase negative Staphylococcus blood stream infections from patients of intensive care units in Mansoura University Hospitals, Egypt.Ann. Clin. Microbiol. Antimicrob.16:63. 10.1186/s12941-017-0238-5
31
McCarthyA. J.LindsayJ. A. (2012). The distribution of plasmids that carry virulence and resistance genes in Staphylococcus aureus is lineage associated.BMC Microbiol.12:104. 10.1186/1471-2180-12-104
32
McCarthyA. J.LoefflerA.WitneyA. A.GouldK. A.LloydD. H.LindsayJ. A. (2014). Extensive horizontal gene transfer during Staphylococcus aureus co-colonization in vivo.Genome Biol. Evol.62697–2708. 10.1093/gbe/evu214
33
MontelongoC.MoresC. R.PutontiC.WolfeA. J.AbouelfetouhA. (2021). Phylogenomic study of Staphylococcus aureus and Staphylococcus haemolyticus clinical isolates from Egypt.bioRxiv [Preprint]10.1101/2021.05.01.442118
34
MorikawaK.TakemuraA. J.InoseY.TsaiM.Nguyen ThiL. T.OhtaT.et al (2012). Expression of a cryptic secondary sigma factor gene unveils natural competence for DNA transformation in Staphylococcus aureus.PLoS Pathogens8:e1003003. 10.1371/journal.ppat.1003003
35
OttoM. (2013). Coagulase-negative staphylococci as reservoirs of genes facilitating MRSA infection: staphylococcal commensal species such as Staphylococcus epidermidis are being recognized as important sources of genes promoting MRSA colonization and virulence.BioEssays354–11. 10.1002/bies.201200112
36
PriceM. N.DehalP. S.ArkinA. P. (2010). FastTree 2 – approximately maximum-likelihood trees for large alignments.PLoS One5:e9490. 10.1371/journal.pone.0009490
37
SabatA. J.PournarasS.AkkerboomV.TsakrisA.GrundmannH.FriedrichA. W. (2015). Whole-genome analysis of an oxacillin-susceptible CC80 mecA-positive Staphylococcus aureus clinical isolate: insights into the mechanisms of cryptic methicillin resistance.J. Antimicrob. Chemother.702956–2964. 10.1093/jac/dkv210
38
ShahkaramiF.RashkiA.GhalehnooZ. R. (2014). Microbial susceptibility and plasmid profiles of methicillin-resistant Staphylococcus aureus and methicillin-susceptible S. aureus.Jundishapur J. Microbiol.7:e16984. 10.5812/jjm.16984
39
ShearerJ. E. S.WiremanJ.HostetlerJ.ForbergerH.BormanJ.GillJ.et al (2011). Major families of multiresistant plasmids from geographically and epidemiologically diverse staphylococci.G31581–591. 10.1534/g3.111.000760
40
TongS. Y. C.DavisJ. S.EichenbergerE.HollandT. L.FowlerV. G. (2015). Staphylococcus aureus infections: epidemiology, pathophysiology, clinical manifestations, and management.Clin. Microbiol. Rev.28603–661. 10.1128/CMR.00134-14
41
VestergaardM.FreesD.IngmerH. (2019). Antibiotic resistance and the MRSA problem.Microbiol. Spectr.7:GPP3-0057-2018. 10.1128/microbiolspec.gpp3-0057-2018
42
WattamA. R.DavisJ. J.AssafR.BoisvertS.BrettinT.BunC.et al (2017). Improvements to PATRIC, the all-bacterial bioinformatics database and analysis resource center.Nucleic Acids Res.45D535–D542. 10.1093/nar/gkw1017
43
XuZ.ShahH. N.MisraR.ChenJ.ZhangW.LiuY.et al (2018). The prevalence, antibiotic resistance and mecA characterization of coagulase negative staphylococci recovered from non-healthcare settings in London, UK.Antimicrob. Resist. Infect. Control7:73. 10.1186/s13756-018-0367-4
44
ZankariE.HasmanH.CosentinoS.VestergaardM.RasmussenS.LundO.et al (2012). Identification of acquired antimicrobial resistance genes.J. Antimicrob. Chemother.672640–2644. 10.1093/jac/dks261
45
ZetolaN.FrancisJ. S.NuermbergerE. L.BishaiW. R. (2005). Community-acquired meticillin-resistant Staphylococcus aureus: an emerging threat.Lancet Infect. Dis.5275–286. 10.1016/S1473-3099(05)70112-2
Summary
Keywords
antibiotic resistance, Staphylococcus, plasmids, incompatibility group, Rep, horizontal gene transfer
Citation
Mores CR, Montelongo C, Putonti C, Wolfe AJ and Abouelfetouh A (2021) Investigation of Plasmids Among Clinical Staphylococcus aureus and Staphylococcus haemolyticus Isolates From Egypt. Front. Microbiol. 12:659116. doi: 10.3389/fmicb.2021.659116
Received
27 January 2021
Accepted
07 May 2021
Published
04 June 2021
Volume
12 - 2021
Edited by
Alla Lapidus, Saint Petersburg State University, Russia
Reviewed by
Ben Pascoe, University of Bath, United Kingdom; Bartolome Moya Canellas, Instituto de Investigación Sanitaria de Palma (IdISPa), Spain
Updates
Copyright
© 2021 Mores, Montelongo, Putonti, Wolfe and Abouelfetouh.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Alaa Abouelfetouh, Alaa.abouelfetouh@pharmacy.alexu.edu.eg
This article was submitted to Evolutionary and Genomic Microbiology, a section of the journal Frontiers in Microbiology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.