ORIGINAL RESEARCH article

Front. Microbiol., 24 November 2022

Sec. Evolutionary and Genomic Microbiology

Volume 13 - 2022 | https://doi.org/10.3389/fmicb.2022.1016829

Genomic analysis of methicillin-resistant Staphylococcus aureus clonal complex 239 isolated from Danish patients with and without an international travel history

  • 1. Laboratory of Medical Microbiology, Vaccine and Infectious Disease Institute, University of Antwerp, Antwerp, Belgium

  • 2. Department of Bacteria, Parasites and Fungi, Statens Serum Institut, Copenhagen, Denmark

  • 3. Cambridge Institute of Therapeutic Immunology and Infectious Disease, University of Cambridge School of Clinical Medicine, Cambridge Biomedical Campus, Cambridge, United Kingdom

Article metrics

View details

4

Citations

2,4k

Views

745

Downloads

Abstract

Introduction:

International travel has been a major determinant for the introduction of pathogens such as methicillin-resistant Staphylococcus aureus (MRSA) into naïve geographic areas. MRSA clonal complex 239 (CC239) is a highly virulent clone that is predominant in Asia. The objective of this study was to determine the geographic origin of MRSA CC239 isolates recovered from Danish cases with or without a history of international travel during 2004–2016.

Materials and methods:

Human MRSA isolates with spa types t030 and t037 (n = 60) were obtained from the National Reference Laboratory for Antimicrobial Resistance. For each case, the following data were collected from notification forms: sex, age, isolation year, specimen source (screening swab or clinical sample), infection type, and international travel history. All isolates were whole-genome sequenced, and a comparative genome and phylogenetic analysis was performed.

Results:

The majority of isolates originated from skin and soft tissue (SST) infections and screening swabs. In 31 out of 60 cases reported international travel to different parts of the world. Fifty-four isolates belonged to CC239, including sequence type 239 (ST239) (n = 43), ST241 (n = 5), ST4377 (n = 2), ST4378 (n = 1), ST1465 (n = 1), ST343 (n = 1), and ST592 (n = 1). The majority of the CC239 MRSA isolates (40/54) belonged to well-known geographic clades, including the Asian (n = 12), Serbian (n = 11), South American (n = 2), and Turkish (n = 15). Most MRSA ST239 isolates belonging to the highly virulent Asian clade carried sasX and were recovered from individuals who had travelled to Asia, Africa and the Middle East.

Conclusion:

Our data reveal multiple introductions of MRSA CC239 into Denmark through international travel, which highlights the importance of continued genomic surveillance of MRSA in persons returning from international travel to areas where MRSA is endemic.

Introduction

Methicillin-resistant Staphylococcus aureus (MRSA) is a major cause of difficult-to-treat infections around the world. Some MRSA clones have evolved into multiple clades which can be distinguished by core genome single nucleotide polymorphisms (SNPs) and are correlated with geographical origin (Lee et al., 2018). In Northern Europe, especially in Scandinavian countries, infections due to MRSA are relatively rare, accounting for ~2% of all S. aureus bloodstream infections (One Health Trust, 2022), but in recent times a steady increase in MRSA infections has been noted Danish Programme for surveillance of antimicrobial consumption and resistance in bacteria from food animals, food and humans DANMAP 2021.1 In 2020, ~17% (215/1281) of all Danish patients with MRSA infection had likely acquired MRSA outside Denmark, supporting that international travel is a risk factor for the introduction of MRSA into Denmark (Larsen et al., 2008a; Zhou et al., 2014). A longitudinal study from Sweden showed a ~ 10-fold rise in MRSA cases between 2000 and 2010, of which ~76% of the infecting MRSA isolates originated from outside Sweden (Larsson et al., 2014). The most common route of acquisition is through direct skin contact with colonized individuals, fomites or animals during an overseas trip (Zhou et al., 2014).

MRSA clonal complex 239 (CC239) comprises the sequence type 239 (ST239) clone, which is widespread globally (Oliveira et al., 1998; Ko et al., 2005; Harris et al., 2010; Li et al., 2011). Several ST239 variants have been previously described, including the Brazilian, Hungarian, Portuguese, Viennese, Indian, Asian, and Eurasian clades (Harris et al., 2010; Wang et al., 2014; Monecke et al., 2018; De Backer et al., 2019). ST239 is endemic in Southeast Asia where it causes more than 70% of hospital-acquired MRSA infections (Ko et al., 2005). Asian ST239 clade isolates are highly virulent and harbour different virulence factors, such as the sasX gene present on a φSPβ-like prophage and encoding a virulence factor involved in colonization and immune evasion of the bacterium (Wang et al., 2014; Nakaminami et al., 2017; De Backer et al., 2019), which might have facilitated its pandemic spread across Asia.

Here, we investigated the population structure and origin of MRSA CC239 isolates recovered from 60 Danish cases with or without a history of international travel during 2004–2016.

Materials and methods

Clinical information and study isolates

The study population included all Danish cases who tested positive for MRSA with the CC239-associated spa types t030 or t037 during 2004–2016. For each case, the following data were collected from notification forms: sex, age, isolation year, specimen source (screening swab or clinical sample), infection type, and international travel history 12 months prior to diagnosis. Detailed data on the characteristics of MRSA isolates were extracted from the national MRSA database. The regional clinical microbiology laboratories performed S. aureus identification and established methicillin resistance by use of standard laboratory methods and forwarded the isolates to the National Reference Laboratory for Antimicrobial Resistance at Statens Serum Institut, where MRSA identification was confirmed by PCR detection of mecA and spa typing (Larsen et al., 2008b). Collection and use of data were approved by the Danish Data Protection Agency (protocol no. 2001-14-0021).

Genotypic characterization

Multilocus sequence typing (MLST) was performed as previously described (Enright et al., 2000). STs and CCs were assigned through the Staphylococcus aureus PubMLST database. All isolates were tested for the presence of sasX by PCR as described previously (Li et al., 2012). Positive isolates were further investigated for the presence of the φSPβ-like prophage by PCR and Sanger sequencing as described elsewhere (De Backer et al., 2019).

Whole-genome sequencing and bioinformatics analyses

All isolates were subjected to DNA extraction using the Master Pure Complete DNA and RNA purification kit (Epicentre, Madison, WI, United States), according to the manufacturer’s protocol. We used the Illumina Nextera XT kit for library construction and paired-end reads of 250 bp each were generated on a Miseq platform (Illumina Inc., San Diego, United States). The raw sequences were analyzed using the bacterial whole-genome sequencing analysis pipeline Bacpipe v.1.2 (Xavier et al., 2020). The online server2 was used to identify prophages (Arndt et al., 2016). To study the genetic relatedness of isolates, we have utilized gene-by-gene approach using chewBBACA (Silva et al., 2018). we generated a study-specific ST239 scheme based on 27 publicly available MRSA ST239 genomes representing each geographic clades (Supplementary Table S1). Phylogenetic relationship was visualized using iTOL (Letunic and Bork, 2019), and geographic clades were defined based on published MRSA ST239 isolates from (Harris et al., 2010; De Backer et al., 2019). All sequence data from this study have been deposited under BioProject ID: PRJNA674016.

Results

Clinical information

The 60 isolates originated from SST (n = 22), respiratory tract (n = 4), blood (n = 2), urine (n = 1), mammary gland (n = 1), eye (n = 1), screening swabs (n = 22), and unknown source (n = 7) (Table 1). Of the 60 cases, 31 had a history of travel to Balkan (n = 5), Africa (n = 3), North Africa (n = 2), South America (n = 2), South Asia (n = 4), Southeast Asia (n = 3) and the Middle East (n = 12), whereas 26 had no history of international travel. Three cases had unknown travel histories (Figures 1A,B).

Table 1

Strain IDClinical informationIsolate information
Isolation yearSpecimenbTravelRegionspa typeSTCCCC239 clade csasX
560182007Respiratory tractNo travelNo travelt037ST239CC239Asian−
718552010BloodMalaysiaSoutheast Asiat037ST239CC239Asian+
569672007EyeNo travelNo travelt037ST239CC239Asian+
611332008SSTNo travelNo travelt037ST239CC239Asian+
611342008SSTNo travelNo travelt037ST239CC239Asian+
716792010SSTNo travelNo travelt037ST239CC239Asian+
719712010SSTNo travelNo travelt037ST239CC239Asian+
1030522015SSTIranMiddle Eastt037ST239CC239Asian+
719272010UnknownNo travelNo travelt037ST239CC239Asian+
563892007ScreeningGambiaAfricat037ST239CC239Asian+
818942011ScreeningIranMiddle Eastt037ST239CC239Asian+
7196872010ScreeningNo travelNo travelt037ST239CC239Asian+
891022013Respiratory tractNo travelNo travelt037ST239CC239Serbian−
891102013Respiratory tractNo travelNo travelt037ST239CC239Serbian−
789382011SSTLebanonMiddle Eastt037ST239CC239Serbian−
913542014SSTNo travelNo travelt037ST239CC239Serbian−
822002012SSTSerbia/MontenegroEuropet037ST239CC239Serbian−
1098802016SSTSyriaMiddle Eastt037ST239CC239Serbian−
1094782016ScreeningBosnia and HerzegovinaEuropet037ST239CC239Serbian−
680782009ScreeningEgyptNorth Africat037ST239CC239Serbian−
787702011ScreeningLebanonMiddle Eastt037ST239CC239Serbian−
1076732016ScreeningSerbia/MontenegroEuropet037ST239CC239Serbian−
1062272016ScreeningThailandSoutheast Asiat037ST239CC239Serbian−
564292007Mammary glandPeruSouth Americat037ST239CC239South America−
756332011SSTBrazilSouth Americat037ST239CC239South America−
610012008BloodNo travelNo travelt030ST239CC239Turkish−
892692013Respiratory tractTurkeyMiddle Eastt030ST239CC239Turkish−
855932013SSTNo travelNo travelt030ST239CC239Turkish−
1082852016SSTIndiaSouth Asiat030ST239CC239Turkish−
993592015SSTIranMiddle Eastt030ST239CC239Turkish−
576642007SSTPakistanSouth Asiat030ST239CC239Turkish−
1006652016SSTRomaniaEuropet030ST239CC239Turkish−
922142014UnknownIndiaSouth Asiat030ST239CC239Turkish−
564272007UnknownRomaniaEuropet030ST239CC239Turkish−
451372005UnknownUnknownUnknownt030ST239CC239Turkish−
1106142016UnknownUnknownUnknownt030ST239CC239Turkish−
1104802016UrineNo travelNo travelt030ST239CC239Turkish−
894482013ScreeningNo travelNo travelt030ST239CC239Turkish−
277902009ScreeningPakistanSouth Asiat030ST239CC239Turkish−
669102009ScreeningTurkeyMiddle Eastt030ST239CC239Turkish−
836272012SSTNo travelNo travelt037ST239CC239NA−
731982010ScreeningIraqMiddle Eastt037ST239CC239NA−
1009402015ScreeningLebanonMiddle Eastt037ST239CC239NA+
550172007SSTThailandSoutheast Asiat037ST343CC239ND+
555262007UnknownNo travelNo travelt037ST4377CC239ND−
408362004UnknownUnknownUnknownt037ST592CC239ND+
700602010ScreeningBahrainMiddle Eastt037ST241CC239ND−
571932007ScreeningIranMiddle Eastt030ST1465CC239ND−
554612007ScreeningEgyptNorth Africat037ST241CC239ND+
864962013ScreeningKenyaAfricat037ST241CC239ND+
1030462015ScreeningKenyaAfricat037ST241CC239ND+
561932007ScreeningNo travelNo travelt037ST241CC239ND+
568352007ScreeningNo travelNo travelt037ST4377CC239ND+
579632007ScreeningNo travelNo travelt037ST4378CC239ND+
579592007ScreeningNo travelNo travelt037ST80CC80ND−
574012007SSTNo travelNo travelt030ST1456CC30ND−
578912007SSTNo travelNo travelt030ST8CC80ND−
618752008SSTNo travelNo travelt037ST30CC30ND−
619082008SSTNo travelNo travelt037ST30CC30ND−
822592012SSTNo travelNo travelt037ST30CC30ND−

Clinical information and characteristics of isolates used in this study.

Figure 1

Figure 1

Geographic origin and sequence type (ST) distribution (A) Travel destinations of individuals (n = 31) entering Denmark used in this study. Colors represent geographic regions of travel destinations. (B) ST type distribution of isolates (n = 60) by specimen type. SST, skin and soft tissue.

Population structure

Of the 60 MRSA isolates, 54 belonged to STs within CC239, including ST239 (n = 43), ST241 (n = 5), ST343 (n = 1), ST592 (n = 1), ST1465 (n = 1), ST4377 (n = 2), and ST4378 (n = 1) (Table 1). The remaining isolates belonged to ST30 (n = 3) and ST1456 (n = 1) within CC30, and to ST8 within CC8 (n = 1) and ST80 (n = 1) within CC80.

Population structure of MRSA ST239

Allelic-loci differences were calculated among the 43 MRSA ST239 genomes from this study and the collection of 27 publicly available MRSA ST239 genomes, and a core-genome phylogeny was constructed (Figure 2). In total, 3,123 loci were determined, from which 921 accessory-genome loci were removed and 2,161 core-genome loci were used to calculate the differences between the strains. The MRSA ST239 isolates were assigned to the Turkish (n = 15), Asian (n = 12), Serbian (n = 11) and South American (n = 2) clades (Figure 2). Three isolates could not be assigned to a clade.

Figure 2

Figure 2

Core-genome phylogeny of all MRSA ST239 isolates isolated from Denmark (n = 40) and published context sequences (n = 27). Clustering is done in line with Harris et al. (2010) and De Backer et al. (2019) clades are indicated by colored overlays. Metadata is visualized in the lanes in compliance to the legend, from left to right: (i) Origin of the isolate; (ii) Year of isolation; (iii) Region of travel; (iv) sasX carriage; and (v) Carriage of antimicrobial resistance genes spc, blaZ, tet(M), dfrG, and erm(A) with black signifying carriage and white non-carriage. The context sequences were downloaded from Genbank and Sanger Institute. Branch lengths represent the substitution rate per position in the core-genome alignment as indicated in the scale bar. The tree is publicly available on MicroReact (https://microreact.org/project/oxRfLg4vYw4KqeXRdUXb33).

MRSA ST239 isolates belonging to the Asian clade clustered closely with published Asian clade isolates (S2, Z172, DEN907, TW20, V605, AMRC01, XN108, and V521). Isolates from the same country and different years were closely related, for instance isolates originating from Kenya and Iran, respectively (Figure 2). Interestingly, several isolates obtained from patients with no international travel history (n = 8) belonged to the Asian clade, and these clustered according to the isolation year (Figure 2). None of the study isolates were closely related to DEN907 (Figure 2), which is a known Asian clade reference strain isolated in Denmark from a Thai patient (Harris et al., 2010).

Additionally, 15 isolates from individuals who reported travelling to Turkey (n = 2), Iran (n = 1), India (n = 2), Pakistan (n = 2), and Romania (n = 2), no international travel history (n = 4) and unknown (n = 2) clustered in the Turkish clade (Figure 2). Interestingly, Turkish clade isolates obtained from cases who reported travel to Pakistan and India formed a distinct substructure, showing higher genetic relatedness to the Turkish clade reference strain, T0131 from China, than the TUR-1 isolated from Turkey.

Isolates from cases who reported travel to Serbia (n = 2), Lebanon (n = 2), Thailand (n = 1), Bosnia Herzegovina (n = 1), Syria (n = 1), Egypt (n = 1) and no travel (n = 3) clustered within the Serbian clade (Figure 2). These were closely related to Serbian isolates (SE1153NT1, SE1895NT2, SE1890NT3, SE1582NT2). Three monophyletic substructures were observed within the Serbian clade grouping non-European isolates separately from European and no travel isolates. The study isolates from cases who reported travel to Brazil and Peru clustered according to their geographical origin, namely with the South American clade reference strains BMB9393 and AGT9 (Figure 2).

We found that the sasX gene was present in 11 out of 12 MRSA ST239 isolates from Asia, but not any other ST239 clade isolates (n = 28) (Figure 2 and Table 1). In addition, sasX was present in eight of the 14 MRSA CC239 belonging to STs other than ST239 (Table 1).

Accessory genomic content of the MRSA ST239 isolates

Within each clade, isolates showed similarities in their core (n = 2,161) and accessory (n = 962) allelic-loci. All isolates from the Asian cluster (n = 12) showed presence of antimicrobial resistance genes spc, blaZ, tet(M), dfrG and erm(A), encoding resistance to, respectively, aminoglycoside, beta-lactam, tetracycline, trimethoprim and macrolide antibiotics (Figure 2; Supplementary Table S2). Isolates in this clade harboured 1–3 plasmid replicons and multiple virulence genes, including lukD, lukE, hlgB, sek and seq (Supplementary Table S2). The rep21 plasmid replicon was predominantly present, except for one isolate (81894) (Supplementary Table S2).

In contrast, the Serbian clade showed only presence of the blaZ (11/11) and tet(M) (7/11) antimicrobial resistance genes (Figure 2; Supplementary Table S2), and contained fewer virulence genes than isolates from the Asian clade, namely lukD, lukE, hlgB, sek and seq (Supplementary Table S2). Within the Serbian clade, plasmid presence was similar to the Asian clade with individual isolates carrying up to three different plasmid rep genes. rep10 was the most prevalent rep type (8/11), followed by rep7 (6/11) and rep21 (4/11) (Supplementary Table S2).

Finally, isolates from the Turkish (n = 15) and South American (n = 2) clades harboured, on average, fewer plasmids replicons (0–2 rep types). All isolates carried spc, blaZ and tet(M) antimicrobial resistance genes and virulence genes lukD, lukE, sek and seq (Figure 2; Supplementary Table S2).

Discussion

This study reveals multiple introductions of MRSA CC239 into Denmark, and a predominance of the highly virulent Asian clade having been introduced by travel to Asia but also to Africa and the Middle East. The finding of genetically similar isolates from cases without any recent travel history indicates either an underreporting of travel or the spread of imported MRSA in the community. This study highlights the distinct clustering identified among CC239 MRSA isolated from Danish patients with an international travel history during 2004–2016. In total, isolates in four geographically defined clades were identified, i.e., the Serbian, Turkish, Asian and South American clades defined by previous studies (Harris et al., 2010; De Backer et al., 2019). Plasmid replicon abundance was higher in the Asian and Serbian clades as compared to the Turkish and South American ones, but some virulence factors were found across all clades. The latter two clades shared the same antimicrobial resistance genes while the other two had distinct antimicrobial resistance gene patterns.

ST239 S. aureus is an atypical member of canonical CC8, as it has unique features that distinguish it from other members of the clonal complex. These features include the presence of the arcC-allele, as well as other genes that are not typically found in CC8 strains. The difference in allele frequencies between ST239 and other CC8 strains suggests that this strain may have experienced positive selection, which could be due to the adaptive advantages conferred by the acquired genes. Further research is needed to determine the precise role that these genes play in the ecology and pathogenesis of ST239. The analysis of clonal complex sub-groups revealed distinct genotypic signatures that were used to classify individual CCs. ST239 strains were classified as CC239 rather than CC8 (Al Laham et al., 2015; Monecke et al., 2018).

In this study, the Turkish (Eurasian) clade was revealed to be most dominant with isolates from cases who had travelled across Asia, Africa and the Middle East. Importation of CC239 from these regions into Scandinavia has been described before (Stenhem et al., 2010). Interestingly, CC239 was also isolated from non-traveling cases (n = 19) with the majority (n = 15) being part of either the Asian, Turkish and Serbian clades. Our Asian clade isolates did not cluster closely with DEN907, the Asian clade reference from Denmark, supporting that isolates might have been imported from the individual’s country of travel. These findings also imply a need for careful interpretation of clades named after the region where they were originally identified.

We identified an Asian clade isolate from a case who reported travel to Africa (56389). In recent years, there have been reports of Asian clade ST239 isolates in Africa (Abdulgader et al., 2015; Founou et al., 2018). This is not surprising, given the increasing global interconnectedness of the world today. This finding also suggests that there could be more ST239 isolates in Africa than have been found to date. This might be because whole-genome sequencing data of African MRSA is scarce and these isolates were not yet analyzed up to the clade level (Abdulgader et al., 2015). A recent study found ST239 to be prevalent in Egypt although whole genome sequencing was not performed making it difficult to assess the presence of the Asian clade in the African continent (Abd El-Hamid et al., 2022).

Furthermore, we identified a substructure within the Turkish clade that appears to be associated with travel to Asia. These isolates, originating from cases who had travelled to India and Pakistan, clustered with T0131, originating from China, indicating that there might be a Turkish subclade of the clade prevalent in the Asian region (Wang et al., 2014; Liao et al., 2020). However, the presence of isolates from cases without an international travel history within the Asian clade might reflect local transmission events in Denmark.

Our study identified sasX positive isolates outside the Asian clade. sasX is a known marker of the Asian clade and highly prevalent in South and East Asia, while also present in countries with low ST239 prevalence such as Japan (Nakaminami et al., 2017). Carriage of sasX has been shown to be vital to MRSA colonization and pathogenesis and was identified as a driver of MRSA epidemic waves (Li et al., 2012; Al Laham et al., 2015; Monecke et al., 2018) MRSA ST239 is currently not one of the major MRSA clones causing infections in Europe (Murray et al., 2022), but based on its infection-causing potential and its pandemic spread in Asia, our findings highlight the importance of continued genomic surveillance of MRSA in persons travelling to countries where MRSA is prevalent.

Funding

JC was supported by the University of Antwerp doctoral assistant funds. BX was supported by University of Antwerp Research funds (BOF-DOCPRO 2012–27450).

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Statements

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found at: https://www.ncbi.nlm.nih.gov/genbank/, BioProject ID: PRJNA674016.

Author contributions

JC, JL, AL, and SM-K: conceptualization. JC, BX, JV, and CL: experimentation. JC, BX, JV, and SP: bioinformatics and visualization. JC, JV, BX, JL, AL, and SM-K: manuscript writing and data interpretation. SM-K: supervision. All authors contributed to the article and approved the manuscript.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary material

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

References

  • 1

    Abd El-HamidM. I.SewidA. H.SamirM.HegazyW. A. H.BahnassM. M.MosbahR. A.et al. (2022). Clonal diversity and epidemiological characteristics of ST239-MRSA strains. Front. Cell. Infect. Microbiol.12:782045. doi: 10.3389/fcimb.2022.782045

  • 2

    AbdulgaderS. M.ShittuA. O.NicolM. P.KabaM. (2015). Molecular epidemiology of methicillin-resistant Staphylococcus aureus in Africa: a systematic review. Front. Microbiol.6:348. doi: 10.3389/fmicb.2015.00348

  • 3

    Al LahamN.MediavillaJ. R.ChenL.AbdelateefN.ElamreenF. A.GinocchioC. C.et al. (2015). MRSA clonal complex 22 strains harboring toxic shock syndrome toxin (TSST-1) are endemic in the primary hospital in Gaza, Palestine. PLoS One10:e0120008. doi: 10.1371/journal.pone.0120008

  • 4

    ArndtD.GrantJ. R.MarcuA.SajedT.PonA.LiangY.et al. (2016). PHASTER: a better, faster version of the PHAST phage search tool. Nucleic Acids Res.44, W16–W21. doi: 10.1093/nar/gkw387

  • 5

    De BackerS.XavierB. B.VanjariL.CoppensJ.LammensC.VemuL.et al. (2019). Remarkable geographical variations between India and Europe in carriage of the staphylococcal surface protein-encoding sasX/sesI and in the population structure of methicillin-resistant Staphylococcus aureus belonging to clonal complex 8. Clin. Microbiol. Infect.25, 628.e1–628.e7. doi: 10.1016/j.cmi.2018.07.024

  • 6

    EnrightM. C.DayN. P.DaviesC. E.PeacockS. J.SprattB. G. (2000). Multilocus sequence typing for characterization of methicillin-resistant and methicillin-susceptible clones of Staphylococcus aureus. J. Clin. Microbiol.38, 1008–1015. doi: 10.1128/JCM.38.3.1008-1015.2000

  • 7

    FounouL. L.AmoakoD. G.FounouR. C.EssackS. Y. (2018). Antibiotic resistance in food animals in Africa: a systematic review and meta-analysis. Microb. Drug Resist.24, 648–665. doi: 10.1089/mdr.2017.0383

  • 8

    HarrisS. R.FeilE. J.HoldenM. T.QuailM. A.NickersonE. K.ChantratitaN.et al. (2010). Evolution of MRSA during hospital transmission and intercontinental spread. Science327, 469–474. doi: 10.1126/science.1182395

  • 9

    KoK. S.LeeJ. Y.SuhJ. Y.OhW. S.PeckK. R.LeeN. Y.et al. (2005). Distribution of major genotypes among methicillin-resistant Staphylococcus aureus clones in Asian countries. J. Clin. Microbiol.43, 421–426. doi: 10.1128/JCM.43.1.421-426.2005 PMID:

  • 10

    LarsenA. R.BöcherS.SteggerM.GoeringR.PallesenL. V.SkovR. (2008a). Epidemiology of European community-associated methicillin-resistant Staphylococcus aureus clonal complex 80 type IV strains isolated in Denmark from 1993 to 2004. J. Clin. Microbiol.46, 62–68. doi: 10.1128/JCM.01381-07

  • 11

    LarsenA. R.SteggerM.SørumM. (2008b). Spa typing directly from a mecA, spa and pvl multiplex PCR assay-a cost-effective improvement for methicillin-resistant Staphylococcus aureus surveillance. Clin. Microbiol. Infect.14, 611–614. doi: 10.1111/j.1469-0691.2008.01995.x

  • 12

    LarssonA. K.GustafssonE.JohanssonP. J.OdenholtI.PeterssonA. C.MelanderE. (2014). Epidemiology of MRSA in southern Sweden: strong relation to foreign country of origin, health care abroad and foreign travel. Eur. J. Clin. Microbiol. Infect. Dis.33, 61–68. doi: 10.1007/s10096-013-1929-2

  • 13

    LeeA. S.De LencastreH.GarauJ.KluytmansJ.Malhotra-KumarS.PeschelA.et al. (2018). Methicillin-resistant Staphylococcus aureus. Nat. Rev. Dis. Primers.4:18033. doi: 10.1038/nrdp.2018.33

  • 14

    LetunicI.BorkP. (2019). Interactive tree of life (iTOL) v4: recent updates and new developments. Nucleic Acids Res.47, W256–W259. doi: 10.1093/nar/gkz239

  • 15

    LiY.CaoB.ZhangY.ZhouJ.YangB.WangL. (2011). Complete genome sequence of Staphylococcus aureus T0131, an ST239-MRSA-SCCmec type III clone isolated in China. J. Bacteriol.193, 3411–3412. doi: 10.1128/JB.05135-11

  • 16

    LiM.DuX.VillaruzA. E.DiepB. A.WangD.SongY.et al. (2012). MRSA epidemic linked to a quickly spreading colonization and virulence determinant. Nat. Med.18, 816–819. doi: 10.1038/nm.2692

  • 17

    LiaoF.MoZ.GuW.XuW.FuX.ZhangY. (2020). A comparative genomic analysis between methicillin-resistant Staphylococcus aureus strains of hospital-acquired and community infections in Yunnan province of China. BMC Infect Dis.20:137. doi: 10.1186/s12879-020-4866-6

  • 18

    MoneckeS.SlickersP.GawlikD.MüllerE.ReissigA.Ruppelt-LorzA.et al. (2018). Molecular typing of ST239-MRSA-III from diverse geographic locations and the evolution of the SCCmec III element during its intercontinental spread. Front. Microbiol.9:1436. doi: 10.3389/fmicb.2018.01436

  • 19

    MurrayC. J. L.IkutaK. S.ShararaF.SwetschinskiL.Robles AguilarG.GrayA.et al. (2022). Global burden of bacterial antimicrobial resistance in 2019: a systematic analysis. Lancet399, 629–655. doi: 10.1016/S0140-6736(21)02724-0

  • 20

    NakaminamiH.ItoT.HanX.ItoA.MatsuoM.UeharaY.et al. (2017). First report of sasX-positive methicillin-resistant Staphylococcus aureus in Japan. FEMS Microbiol. Lett.364:fnx171. doi: 10.1093/femsle/fnx171

  • 21

    OliveiraD.Santos-SanchesI.MatoR.TamayoM.RibeiroG.CostaD.et al. (1998). Virtually all methicillin-resistant Staphylococcus aureus (MRSA) infections in the largest Portuguese teaching hospital are caused by two internationally spread multiresistant strains: the 'Iberian' and the 'Brazilian' clones of MRSA. Clin. Microbiol. Infect.4, 373–384. doi: 10.1111/j.1469-0691.1998.tb00081.x

  • 22

    One Health Trust (2022). ResistanceMap: Antibiotic Resistance. Available at: https://resistancemap.onehealthtrust.org/AntibioticResistance.php (Accessed May 31, 2022).

  • 23

    SilvaM.MachadoM. P.SilvaD. N.RossiM.Moran-GiladJ.SantosS.et al. (2018). chewBBACA: a complete suite for gene-by-gene schema creation and strain identification. Microb. Genome.4:e000166. doi: 10.1099/mgen.0.000166

  • 24

    StenhemM.OrtqvistA.RingbergH.LarssonL.Olsson LiljequistB.HaeggmanS.et al. (2010). Imported methicillin-resistant Staphylococcus aureus, Sweden. Emerg. Infect. Dis.16, 189–196. doi: 10.3201/eid1602.081655

  • 25

    WangZ.ZhouH.WangH.ChenH.LeungK. K.TsuiS.et al. (2014). Comparative genomics of methicillin-resistant Staphylococcus aureus ST239: distinct geographical variants in Beijing and Hong Kong. BMC Genomics15:529. doi: 10.1186/1471-2164-15-529

  • 26

    XavierB. B.MysaraM.BolzanM.Ribeiro-GonçalvesB.AlakoB. T. F.HarrisonP.et al. (2020). BacPipe: a rapid, user-friendly whole-genome sequencing pipeline for clinical diagnostic bacteriology. iScience23:100769. doi: 10.1016/j.isci.2019.100769

  • 27

    ZhouY. P.Wilder-SmithA.HsuL. Y. (2014). The role of international travel in the spread of methicillin-resistant Staphylococcus aureus. J. Travel Med.21, 272–281. doi: 10.1111/jtm.12133

Summary

Keywords

MRSA, CC239, ST239, sasX, Denmark, WGS

Citation

Coppens J, Xavier BB, Vlaeminck J, Larsen J, Lammens C, Van Puyvelde S, Goossens H, Larsen AR and Malhotra-Kumar S (2022) Genomic analysis of methicillin-resistant Staphylococcus aureus clonal complex 239 isolated from Danish patients with and without an international travel history. Front. Microbiol. 13:1016829. doi: 10.3389/fmicb.2022.1016829

Received

11 August 2022

Accepted

07 November 2022

Published

24 November 2022

Volume

13 - 2022

Edited by

Ludmila Chistoserdova, University of Washington, United States

Reviewed by

Shakeel Mowlaboccus, Murdoch University, Australia; Edet E. Udo, Kuwait University, Kuwait

Updates

Copyright

*Correspondence: Surbhi Malhotra-Kumar,

†These authors have contributed equally to this work

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.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics