Skip to main content
  • Research article
  • Open access
  • Published:

A comparative genomic analysis between methicillin-resistant Staphylococcus aureus strains of hospital acquired and community infections in Yunnan province of China

Abstract

Background

Currently, Staphylococcus aureus is one of the most important pathogens worldwide, especially for methicillin-resistant S. aureus (MRSA) infection. However, few reports referred to patients’ MRSA infections in Yunnan province, southwest China.

Methods

In this study, we selected representative MRSA strains from patients’ systemic surveillance in Yunnan province of China, performed the genomic sequencing and compared their features, together with some food derived strains.

Results

Among sixty selective isolates, forty strains were isolated from patients, and twenty isolated from food. Among the patients’ strains, sixteen were recognized as community-acquired (CA), compared with 24 for hospital-acquired (HA). ST6-t701, ST59-t437 and ST239-t030 were the three major genotype profiles. ST6-t701 was predominated in food strains, while ST59-t437 and ST239-t030 were the primary clones in patients. The clinical features between CA and HA-MRSA of patients were statistical different. Compared the antibiotic resistant results between patients and food indicated that higher antibiotic resistant rates were found in patients’ strains. Totally, the average genome sizes of 60 isolates were 2.79 ± 0.05 Mbp, with GC content 33% and 84.50 ± 0.20% of coding rate. The core genomes of these isolates were 1593 genes. Phylogenetic analysis based on pan-genome and SNP of strains showed that five clustering groups were generated. Clustering ST239-t030 contained all the HA-MRSA cases in this study; clustering ST6-t701 referred to food and CA-MRSA infections in community; clustering ST59-t437 showed the heterogeneity for provoking different clinical diseases in both community and hospital. Phylogenetic tree, incorporating 24 isolates from different regions, indicated ST239-t030 strains in this study were more closely related to T0131 isolate from Tianjin, China, belonged to ‘Turkish clade’ from Eastern Europe; two groups of ST59-t437 clones of MRSA in Yunnan province were generated, belonged to the ‘Asian-Pacific’ clone (AP) and ‘Taiwan’ clone (TW) respectively.

Conclusions

ST239-t030, ST59-t437 and ST6-t701 were the three major MRSA clones in Yunnan province of China. ST239-t030 clonal Yunnan isolates demonstrated the local endemic of clone establishment for a number of years, whereas ST59-t437 strains revealed the multi-origins of this clone. In general, genomic study on epidemic clones of MRSA in southwest China provided the features and evolution of this pathogen.

Peer Review reports

Background

Staphylococcus aureus is the most commonly isolated pathogen from human and often lead to endovascular infections, endocarditis, osteomyelitis, septic arthritis and pneumonia [1]. Specifically, methicillin-resistant S. aureus (MRSA) infection has exploded in both healthcare facility and community since 1990s [2]. Although the availability of antimicrobial agents to treat MRSA infections, it still to be the dominant cause of mortality and morbidity all over the world [3]. The spread of different clones from different geographic regions have been reported, MRSA has been verified as one of the most important pathogen in hospital or community setting [4, 5]. Hospital-acquired MRSA (HA-MRSA) infections often lead to ventilator associated pneumonias, intravenous catheter associated infections or surgical wound infections, while community-acquired MRSA (CA-MRSA) infections commonly cause skin and soft tissue infection, sometimes are invasive and life threatening [6, 7].

The whole genome sequencing method has become the gold standard for understanding the genetic diversity of bacteria [8]. Nowadays, large numbers of genomic reports have revealed the molecular epidemiological consequences of clone changes or different infectious syndromes in S. aureus, especially for MRSA [9, 10]. Thus, study on genomic level of MRSA will promote the understanding of the evolution and biology of this pathogen, improving the prevention and management of MRSA. The phylogenetic investigation revealed the dissemination and hospital transmission of ST239 clone through Europe, North America, South America and Asia [11]. Therefore, ST239 clone of strains was the most important HA-MRSA around the world. However, ST59 clone of isolates were the primary CA-MRSA infections in Asia [12, 13]. Furthermore, our previous study [14] on molecular characteristics of S. aureus from food surveillance indicated that ST6-t701 was the prevalent genotypes in southwest China, and some strains had the same genotype profiles with patients’ isolates. The different clonal structures of MRSA strains maybe represent the different infection sources or disease spectrums. In this study, we selected representative MRSA strains from patients’ systemic surveillance in Yunnan province of China, performed the genomic sequencing and compared their features, together with some food derived strains.

Methods

Bacterial source and case definitions

A total of 60 MRSA strains were selected from 2013 to 2016 in both patients and food isolates in Yunnan province, southwest China. All strains were collected consecutively during the same period in both hospital and food surveillance. CA-MRSA strains of patients were isolated from The First People’s Hospital of Yunnan province (hospital A, HA). In contrast, HA-MRSA were isolated from People’s Hospital of Kunming City (hospital B, HB). CA-MRSA case was defined as MRSA isolated within 48 h after admission without the following risk factors: history of hospitalization, surgery, dialysis or stay in a long-term care unit within 1 year; dependence on an indwelling catheter, intravenous line or a percutaneous device when cultured; or previous isolation of MRSA. HA-MRSA cases were defined as MRSA isolation within 48 h with at least one risk factor, as mentioned above, or beyond 48 h, regardless of risk factors [15]. Hospitals A and B are tertiary care academic medical centers in southwest China and cover the entire Kunming area. The basic clinical information of patients from two hospitals were collected. For all these patients’ isolates, we have performed pulsed field gel electrophoresis (PFGE), multilocus sequence typing (MLST) and spa typing and selected them from previous systemic surveillance results (Additional files 1 and 2). The selection criterions of strains in this study were shown as follows: a. three major clones of patients’ MRSA, including ST239-t030, ST59-t437 and ST6-t701 were involved in the study, also contained some highly similarity strains with food derived strains (ST188-t189 and ST965-t062 etc); b. ST239-t030, ST59-t437, ST6-t701 and ST188-t189 generated four cluster groups based on PFGE patterns (Additional file 1), and we selected different PFGE pattern strains for sequencing in each genotype profile, except some highly similar isolates with food strains, such as seven ST59-t437, four ST6-t701 and three ST965-t062 strains.

Antibiotic-resistant test

Minimum inhibitory concentrations (MICs) for the 13 antibiotics were determined through the broth microdilution method using customized microtiter plates (Sensititre, UK) according to the manufacturers’ instructions. The antibiotics tested were penicillin (PCN), oxacillin (OXA), gentamicin (GEN), ciprofloxacin (CIP), levofloxacin (LEV), moxifloxacin (MXF), erythromycin (EM), clindamycin (CM), linezolid (LZD), vancomycin (VAN), tetracycline (TET), rifampicin (RIF) and trimethoprim/sulfamethoxazole (SXT). The tests were performed and interpreted in accordance with the Clinical and Laboratory Standards Institute (CLSI) guidelines (M100, 2018) [16], and S. aureus ATCC 29213 was used as a quality control. The mec gene was detected by PCR as described by previous study [17].

DNA extraction and genome sequencing

All MRSA isolates were recovered on Brain Heart Agar (BHI) (Luqiao, Beijing) at 37 °C for 24 h. Total genomic DNA of the isolated bacteria was extracted using a bacterial total genomic DNA extraction kit (Tiangen, Beijing) following the manufacturer’s instructions. All DNA samples were stored at − 20 °C for complete genome analysis.

Bacterial genome sequencing of all isolates was performed by our laboratory on the Illumina HiSeq platform using 2 × 150 bp paired-end reads. The libraries were built using a Nextera XT DNA Library Prep Kit following the manufacturer’s reference guidelines. Generally, 1 ng genomic DNA of each strain was used. After segment and purification, index PCR was performed to add the Illumina Nextera barcodes using i5 and i7 primers, and then the purification was executed again to remove non-target fragments. Finally, the libraries were normalized, pooled and sequenced using an Illumina HiSeq sequencing system (Illumina, San Diego, USA).

Assembly and annotation

The raw data were trimmed for quality control, and low-quality (<Q40) reads were filtered by Trimmomatic (version 0.38) [18]. Draft genomes were assembled using SOAPdenovo (version 2.04), with k-mer values (25, 31, 37, 47, 59, 71, 83 and 95) optimized to the best assembly results [19, 20]. Gapcloser (version 1.12r6) was used to fill the genomic gaps. GenemarkS software (version 4.28) [21] with default parameters was used to predict the open reading frame (ORF) of each genome, and the predicted amino acid sequences were aligned and annotated by DIAMOND [22] (E-value: 1e-5, top 5) to NCBI non-redundant nucleotide database (NR), SwissProt, TrEMBL, KEGG, GO, COG, CDD, Pfam, Pathogen-Host Interaction database (PHI), Antibiotic Resistance Genes Database (ARDB), and Virulence Factor Database (VFDB).

Comparative genomic analyses

We performed core- and pan-genome analysis using PGAP [23], which employed an all-versus-all BLAST search to group all CDSs into homologous clusters with default settings. All of the 60 strains sharing CDSs were considered to be core genes, and CDSs without orthologs were considered to be unique genes. Orthologs Clusters were annotated in COG functional distributions. Twenty-four MRSA references from other global strains in publically available datasets were used for phylogenetic analysis to determine the dissemination and evolution of these clonal lineages in southwest China (Additional file 3). Pairwise colinear comparisons of all genome sequences were performed using Mummer3 [24]. The pan-based tree was generated based on the presence or absence of genes in the pan-genome using the neighbor-joining method, and the SNP-based tree was generated based on the total genomic sequence results by using neighbor-joining method as well. The phylogenetic trees were visualized and ordered with FigTree (version 1.4.3) and iTol.

Statistics

Statistical analysis was performed using the SPSS software package (version 16.0, IBM, USA). χ2 or Kruskal-Wallis H test were used if appropriate. A P value of < 0.05 was recognized as statistically significant.

The sequence data have been deposited into the National Center for Biotechnology Information (NCBI), https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/ with accession numbers VCEL00000000 to VCGT00000000 (PRJNA543691).

Results

General information of MRSA in this study

For all the MRSA used in this study, forty strains (66.7%) were isolated from patients, and 20 (33.3%) isolated from food. Among the patients’ strains, sixteen (27%) were recognized as CA-MRSA, compared with 24 (40%) for HA-MRSA (Fig. 1a). ST6-t701 (18, 30%), ST59-t437 (16, 26.7%) and ST239-t030 (14, 23.3%) were the three major genotype profiles in this study. Other genotypes contained ST965-t062, ST188-t189 and ST5-t14723, as Fig. 1b shown. ST6-t701 was predominated in food strains, while ST59-t437 and ST239-t030 were the primary clones in patients (Fig. 1b). The comparison of clinical features between CA and HA-MRSA from patients revealed that except for sex, the factors age, clinical manifestations, and sample types were statistically different (P < 0.05), as Table 1 shown. Most of the CA-MRSA cases were children from pediatrics under the 5 years. The genotype profiles of CA-MRSA strains showed highly heterogeneity (Simpson’s index = 0.750) compared with HA-MRSA (Simpson’s index = 0.487). All the genotypes could be found in CA-MRSA isolates, except ST239-t030. However, the HA-MRSA cases were mostly from old patients, and the patients had different primary diseases, such as hypertension, cardiovascular disease, trauma, tumor or other infections. 58.3% of HA-MRSA were ST239-t030, and 41.7% were ST59-t437 genotype. All the food strains were isolated from eight cities of Yunnan province, and the majority (21.7%) was isolated from grain products. The details of MRSA strains information in this study were shown in additional file 4.

Fig. 1
figure 1

The source and genotype profiles of MRSA used in this study. a The source of MRSA strains. b The distributions of genotype profiles of the isolates (* indicated the distribution of genotypes between patients and food was statistically significant)

Table 1 Clinical features between community- and hospital-acquired MRSA in this study

Antibiotic resistant results

All the MRSA in this study were resistant to OXA and PCN, but sensitive to VAN. The mec gene was positive for all the isolates. In general, 55% of all the strains were resistant to MXF, followed by SXT (41.7%) and EM (36.7%). Compared the antibiotic resistant results between patients and food indicated that higher antibiotic resistant rates were found in patients’ strains, such as LZD, CIP, LEV, RIF, TET and MXF (Table 2). The antibiotic resistant results between genotypes of strains were also different. ST59-t437 genotype strains showed higher antibiotic resistant rates in CM, CIP, EM, SXT and TET, while ST239-t030 had higher resistant rate in MXF (P < 0.05). ST6-t701 clone had the lowest resistant rate for most of antibiotics, as Table 2 shown.

Table 2 The antibiotic resistant results in this study

The molecular determinants based on ARDB annotations revealed several resistant genes were commonly found between these strains, such as liaR, mprF, pgsA and rpoC resistant to daptomycin; UhpT, murA and GlpT resistant to fosfomycin; gyrA and parC resistant to fluoroquinolones; gyrB and parE resistant to aminocoumarin; rpoB resistant to rifampicin; rpsL resistant to Streptomycin. The majorities of ErmC, ErmT, ErmY, ErmG, ErmA, ErmB resistant to erythromycin and tetM, tetO, tetS, tet32, tetW, tetQ resistant to tetracycline were found in ST59-t437 isolates.

Genome features

Totally, 682,961,826 raw reads were obtained for all the MRSA strains, 643,281,744 valid reads were retained after quality control and the effective rate was 94.19%. The average genome sizes of MRSA were 2.79 ± 0.05 Mbp, with GC content 33% (Additional file 5). The average predicted genes were 2693.35 ± 65.19, with average length 875.05 ± 7.04 (bp), and the coding rate of all predicted genes were 84.50 ± 0.20% for all MRSA in this study. In addition, 50.72 ± 8.38 average tRNA-coding genes and 7.27 ± 0.92 average rRNA operons were identified. The average sequencing depth was 548.48 ± 92.51, and the average coverage of genome was 93.79 ± 1.28%, as Additional file 5 shown. In general, 2572.03 ± 60.90, 2122.35 ± 34.99, 2037.05 ± 16.12, and 2571.48 ± 60.81 average genes were annotated in NR, PFAM, Swissprot and TrEMBL database respectively. 1896.53 ± 17.79 and 1484.57 ± 8.98 average genes were annotated in GO and KEGG pathways, while 1963.68 ± 24.11 average genes were annotated in COG categories.

Core and pan-genome analysis

All CDSs from all MRSA sequenced in this study were used for core and pan-genome analyses. Core genes were considered as the most conserved and shared by all strains, whereas dispensable genes were shared by unique genes identified in only one of the genomes. The most of the CDSs were highly conserved among MRSA strains, and the core genomes of these isolates were 1593 genes (Fig. 2a). Dispensable genes were distributed ranging from 16 to 50 genes in different strains (Fig. 2a). According to homologous whole clusters COG expansions, 1040 genes were assigned to metabolism, e.g., E (amino acid transport and metabolism), G (carbohydrate transport and metabolism) and P (inorganic ion transport and metabolism). Five hundred sixty-four genes were assigned to information storage and processing category, e.g., K (transcription), L (replication, recombination and repair) and J (translation, ribosomal structure and biogenesis). Three hundred ninety-six genes were assigned to cellular processes and signaling, e.g., M (cell wall/membrane/envelope biogenesis), O (posttranslational modification, protein turnover, chaperones) and V (defense mechanisms), as Fig. 2b shown. Generally, the size of the core and pan-genome depended on the number of strains were considered in the analyses. It was evident that the pan-genome curve converges rapidly, indicating the variation in MRSA strains. However, the core genome size was stable, decreasing initially and reached a plateau at 2000 genes (Fig. 2c). This revealed that MRSA has an open pan-genome in this study.

Fig. 2
figure 2

The core and pan-genome analyses of 60 MRSA in this study. a. The flower diagram of numbers of orthologs cluster for all the strains. The core number indicated the genes shared with all bacteria. Other numbers represented the specific genes in each MRSA strain. b COG categories distributions of whole orthologs cluster. c Core-pan gene family boxplot in this study

Phylogenetic analysis

Phylogenetic analysis based on pan-genome (core and dispensable genes) of MRSA showed that five clustering groups were generated with brown, yellow, green, blue and gray color (Fig. 3a). Interestingly, the different genotype or clones of strains clustered into different groups, thus, clustering group brown was ST6-t701 strains, green was ST59-t437 isolates, blue was ST239-t030, yellow was ST188-t189, gray was ST965-t062 and ST5-t14723 clone. The SNP based phylogenetic tree also generated identical five clustering groups showing with different color in this study (Fig. 3b). Pan-based phylogenetic tree of 84 MRSA strains (60 isolates in this study and 24 references) revealed that different clusters were generated according to the genotype profiles as mentioned above (Fig. 3c). All the ST239-t030 strains in this study were more closely related to T0131 isolate from Tianjin, China, belonged to ‘Turkish clade’ from Eastern Europe, but far away from TW20 of ‘Asian clade’ or other global reference isolates. Two groups of ST59-t437 clones of MRSA in Yunnan province were generated with some reference strains. YNSA7 to YNSA167 group was more closely related to SA40TW isolate, belonged to the ‘Asian-Pacific’ clone (AP), whereas YNSA154 to YNSA210, YNSA53 and YNSA11 group were closely related to SA957, M013 and SA268 strains, belonged to the ‘Taiwan’ clone (TW). Most of the CC5 clonal MRSA (ST5 and ST965 clone; YNSA34 to YNSA31) were separated from reference strains except YNSA453, as Fig. 3c shown.

Fig. 3
figure 3

Phylogenetic trees of MRSA used in this study. a Pan-genome based neighbor-joining tree of all the MRSA. ST6-t701, clustering brown; ST59-t437, clustering green; ST239-t030, clustering blue; ST188-t189, clustering yellow; ST965-t062 and ST5-t14723, clustering gray. b SNP based neighbor-joining tree of all the MRSA. ST6-t701, clustering brown; ST59-t437, clustering green; ST239-t030, clustering blue; ST188-t189, clustering yellow; ST965-t062 and ST5-t14723, clustering gray. c Pan-genome based neighbor-joining tree of 84 MRSA, containing 60 strains in this study and 24 references. ST6-t701, clustering brown; ST59-t437, clustering green; ST239-t030, clustering blue; ST188-t189, clustering yellow; ST965-t062 and ST5-t14723, clustering gray. The reference strains were shown with red font

Clustering ST239-t030 clone comprised all the HA-MRSA cases in this study; it was the primary clone for hospital acquired infection in southwest China. Colinear comparison of strains YNSA73 and YNSA94 in this genotype showed that the chromosome of YNSA73 was collinear with YNSA94 with 99.73% identity (Fig. 4a). ST6-t701 clone (clustering brown) mainly referred to food related diseases in community, since the patients and food strains clustered into identical group, and close similarities could be found among these strains. In this group, YNSA21 was a patient isolates, and YNSA365 was a food strain, the chromosome of YNSA21 was collinear with YNSA365 showing 99.71% identity (Fig. 4b). In addition, YNSA311 and YNSA480 were both isolated from food strains of different sources, their genomes were also collinear with each other showing 99.76% identity (Fig. 4c). ST59-t437 clone represented the most heterogeneous cluster in this study. The CA and HA-MRSA cases both contained ST59-t437 genotype, and some isolates from patients were related with food. Colinear comparison between YNSA115 (food) and YNSA279 (patient) revealed the 99.79% identity of two isolates (Fig. 4d). Therefore, ST59-t437 clone had a high degree of diversity for provoking different clinical spectrum of diseases in both community and hospital.

Fig. 4
figure 4

Colinear comparison of MRSA strains in this study. a Colinear comparison of YNSA73 and YNSA94. The red line and blue line showed matches that were indicated by the upward and downward slopes, respectively. b Colinear comparison of YNSA21 and YNSA365. The red line and blue line showed matches that were indicated by the upward and downward slopes, respectively. c Colinear comparison of YNSA311 and YNSA480. The red line and blue line showed matches that were indicated by the upward and downward slopes, respectively. d Colinear comparison of YNSA115 and YNSA279. The red line and blue line showed matches that were indicated by the upward and downward slopes, respectively

Discussion

Studies on MRSA epidemic or outbreaks have originated from developed countries in Europe, North America and Japan [3, 25, 26]. However, MRSA infections became serious problem in developing countries at present, such as China and India. In China, 60% of S. aureus infections were caused by MRSA [27], among them, ST239 was the most abundant clone for hospital infections. In addition, several MRSA researches reported ST239 isolates recovered from hospital patients in Asia, South America and Eastern Europe [28,29,30,31,32]. Therefore, ST239 strain may represent the predominant MRSA clone lineage for hospital infections worldwide. In this study, we found ST239-t030 was the major clone for hospital infection in southwest China. This strain was prevalent in the local health-care facilities, and it had also been isolated from different departments within the hospitals. Consequently, the strain ST239-t030 should be treated as the most important MRSA clone by the pathogen surveillance and hospital infection and prevention control teams.

In Asian countries, specifically in China, ST239 clone has been identified as the predominant HA-MRSA at present [33]. Previous study showed that ST239 MRSA lineage at least contained five clades worldwide, such as Europe and Australia, North and South America, and Asian clade [34]. In this study, several published ST239 clonal MRSA complete genomes were used for comparison, namely T0131, TW20 and JKD6008. Specifically, T0131 reference strain was isolated from old patient in Tianjin, China in 2006, and closely related to the strains of the ‘Turkish clade’ and the ‘Russia variant’ (16 K) [35]. Wang et al. [36] compared the genomics of MRSA from Beijing and Hong Kong; they revealed that strains from HK clustered into the ‘Asian clade’, while BJ isolates were clustered with strains of the ‘Turkish clade’. They concluded that origins of ST239 lineage in southern and northern China were different. In our study, Yunnan isolates from southwest China were closely related to T0131 and clustered with strains of the ‘Turkish clade’ from Eastern Europe. We considered the similar ST239-t030 clonal Yunnan isolates in this study were more likely to demonstrate the local endemic of primary clone establishment for a number of years. The result also indicated the possibility that ST239-t030 MRSA from Yunnan originated from a common recent ancestor and spread during this period.

Our previous study [14] compared the S. aureus isolates from patients and food, the results showed some strains from patients had identical PFGE patterns, ST and spa types with food isolates. These strains were dominated with ST6-t701, ST5-t14723, ST59-t437 and ST965-t062, specifically for ST6-t701. Currently, ST6-t701 strain has become the predominant clone for Staphylococcal food poisoning (SFP) in China. Yan et al. [37] investigated seven outbreaks of 52 S. aureus in Shen-zhen, China. They found 63.5% of all the strains belonged to ST6-t701 genotype. Li et al. [38] performed the molecular epidemiology of seven SFP outbreaks caused by S. aureus in northwest China; the results also indicated ST6-t701was the primary clone, followed by ST5-t002 and ST59-t172. All these evidences suggested that ST6-t701 of S. aureus has an important role in SFP in China. Similar results could be found in our study, ST6-t701 clone was the most important MRSA for food related diseases in Yunnan province. The close genomic relations between food isolates with patients’, and the similarity among different sources of strains in this study both indicated that ST6-t701 was related to foodborne disease in southwest China.

CA-MRSA were characterized by a diversity of clones at the early stage, and then numbers of predominant clones have established around the world. In the United States, the USA300 epidemic clone emerged as major CA-MRSA strain over the past decades [39]. Molecular typing further identified this isolate as ST8, and most commonly spa type t008. Other studies reported that highly virulent ST80 strain dominated in Europe [40], and ST93 MRSA clone in Australia [41]. All these studies indicated the diversity of CA-MRSA strains from different regions of the world. Recently, the genome sequence of predominant CA-MRSA ST59 isolate was reported in Taiwan, these strains were frequently multi-drug resistant isolates [42]. Li et al. [43] performed the molecular typing of clonal complex 59 MRSA in seven major cities across mainland China. Their results revealed ST59-t437 was the predominant genotype, 81.8% of the strains were CA-MRSA, and 18.2% were HA-MRSA. They concluded CC59 MRSA could lead to CA and HA infections and the majority of infections were in children. The genomic analysis of CA-MRSA strains in our study also showed highly genetic diversity of isolates. ST6-t701, ST59-t437, ST188-t189, ST5-t14723 and ST965-t062 were all referred to CA-MRSA cases. Specifically, ST59-t437 clone exhibited the multiple roles in MRSA infections, such as community acquired, hospital acquired and even related with food strain.

ST59 clonal MRSA has been prevalent in the Asia-Pacific region, including mainland China, Taiwan, Vietnam, Japan and Singapore [13]. Previous studies revealed two ST59 clones of MRSA were circulated in Taiwan: the ‘Taiwan’ clone (TW) and the ‘Asian-Pacific’ clone (AP) [44, 45]. TW clone caused sepsis and infection, especially in children, whereas AP clone usually colonized as commensal in healthy individuals. Particularly, reference SA957 and SA40 isolates were the typical strains for TW and AP clone respectively, according to several studies [44, 46]. On the mainland China, ST59 was also the major MRSA clone; ST59-t437 was the predominant clone among CA-MRSA on the Chinese mainland [43]. Several studies have been reported the complete genomic analysis of representative strains isolated from patients. For example, the reference SA268 [47], a ST59 CA-MRSA strain, was isolated from a case of pneumonia in a young patient in Zhejiang Province of China; the genome analysis revealed it was almost identical to that of the Taiwanese ST59 CA-MRSA strains SA957 and M013. Cheng et al. [48] investigated an outbreak of CA-MRSA in Hong Kong by using whole genome sequencing. P8, P11 and HCW strains were analyzed in their study, and the results revealed complicated transmission between patients, healthcare worker and environment. In our study, two groups of ST59-t437 clones of MRSA in Yunnan province were generated. One group was more closely related to AP clone, whereas another was closely related to TW clone, and the result probably supported the multi-origin theory of CA-MRSA.

This study was the systemic surveillance for MRSA by using whole genome sequencing in Yunnan province, southwest China, and it gave us an eye on MRSA infections both in hospital and community in local areas. However, several limitations were found in this study. Firstly, strains used in this study were selected from our previous molecular typing database, not all the surveillance isolates were sequenced and analyzed; perhaps there was selection bias of the strains. Secondly, the short surveillance period of time for MRSA in Yunnan province only reflected the prevalence during these years. Thirdly, the scale of MRSA monitoring was relatively limited; only two hospitals were involved in this study. Therefore, further studies referred to strains from older collections and across other parts of southwest China, transmission between communities and hospitals, and even the longitudinal surveillance were required to evaluate the widespread or origin of the MRSA in Yunnan province, southwest China.

Conclusions

In this study, a comparative genomic analysis was performed between CA and HA-MRSA infections in southwest China. ST239-t030, ST59-t437 and ST6-t701 were the three major MRSA clones in Yunnan province of China. ST239-t030 clonal Yunnan isolates demonstrated the local endemic of clone establishment for a number of years, whereas ST59-t437 strains revealed the multi-origins of this clone. In general, genomic study on epidemic clones of MRSA in southwest China provided the features and evolution of this pathogen.

Availability of data and materials

All data generated or analysed during this study are included in this published article. The sequence data have been deposited into the National Center for Biotechnology Information (NCBI), https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/ with accession numbers VCEL00000000 to VCGT00000000 (PRJNA543691).

Abbreviations

CA:

community-acquired

HA:

hospital-acquired

MLST:

multilocus sequence typing

MRSA:

methicillin-resistant Staphylococcus aureus

PFGE:

Pulsed-field gel electrophoresis

References

  1. Lowy FD. Staphylococcus aureus infections. N Engl J Med. 1998;339(8):520–32.

    Article  CAS  PubMed  Google Scholar 

  2. David MZ, Daum RS. Community-associated methicillin-resistant Staphylococcus aureus: epidemiology and clinical consequences of an emerging epidemic. Clin Microbiol Rev. 2010;23(3):616–87.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Chambers HF, Deleo FR. Waves of resistance: Staphylococcus aureus in the antibiotic era. Nat Rev Microbiol. 2009;7(9):629–41.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Fluit AC, Wielders CL, Verhoef J, Schmitz FJ. Epidemiology and susceptibility of 3,051 Staphylococcus aureus isolates from 25 university hospitals participating in the European SENTRY study. J Clin Microbiol. 2001;39(10):3727–32.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Otto M. Basis of virulence in community-associated methicillin-resistant Staphylococcus aureus. Annu Rev Microbiol. 2010;64:143–62.

    Article  CAS  PubMed  Google Scholar 

  6. DeLeo FR, Otto M, Kreiswirth BN, Chambers HF. Community-associated meticillin-resistant Staphylococcus aureus. Lancet. 2010;375(9725):1557–68.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Klevens RM, Edwards JR, Tenover FC, McDonald LC, Horan T, Gaynes R. Changes in the epidemiology of methicillin-resistant Staphylococcus aureus in intensive care units in US hospitals, 1992-2003. Clin Infect Dis. 2006;42(3):389–91.

    Article  PubMed  Google Scholar 

  8. Copin R, Shopsin B, Torres VJ. After the deluge: mining Staphylococcus aureus genomic data for clinical associations and host-pathogen interactions. Curr Opin Microbiol. 2018;41:43–50.

    Article  CAS  PubMed  Google Scholar 

  9. Li M, Diep BA, Villaruz AE, Braughton KR, Jiang X, DeLeo FR, Chambers HF, Lu Y, Otto M. Evolution of virulence in epidemic community-associated methicillin-resistant Staphylococcus aureus. Proc Natl Acad Sci U S A. 2009;106(14):5883–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Recker M, Laabei M, Toleman MS, Reuter S, Saunderson RB, Blane B, Torok ME, Ouadi K, Stevens E, Yokoyama M, et al. Clonal differences in Staphylococcus aureus bacteraemia-associated mortality. Nat Microbiol. 2017;2(10):1381–8.

    Article  CAS  PubMed  Google Scholar 

  11. Lakhundi S, Zhang K. Methicillin-Resistant Staphylococcus aureus: Molecular Characterization, Evolution, and Epidemiology. Clin Microbiol Rev. 2018;31(4):e00020-18.

  12. Chen FJ, Lauderdale TL, Huang IW, Lo HJ, Lai JF, Wang HY, Shiau YR, Chen PC, Ito T, Hiramatsu K. Methicillin-resistant Staphylococcus aureus in Taiwan. Emerg Infect Dis. 2005;11(11):1760–3.

    Article  PubMed  Google Scholar 

  13. Chuang YY, Huang YC. Molecular epidemiology of community-associated meticillin-resistant Staphylococcus aureus in Asia. Lancet Infect Dis. 2013;13(8):698–708.

    Article  PubMed  Google Scholar 

  14. Liao F, Gu W, Yang Z, Mo Z, Fan L, Guo Y, Fu X, Xu W, Li C, Dai J. Molecular characteristics of Staphylococcus aureus isolates from food surveillance in Southwest China. BMC Microbiol. 2018;18(1):91.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  15. Fridkin SK, Hageman JC, Morrison M, Sanza LT, Como-Sabetti K, Jernigan JA, Harriman K, Harrison LH, Lynfield R, Farley MM. Methicillin-resistant Staphylococcus aureus disease in three communities. N Engl J Med. 2005;352(14):1436–44.

    Article  CAS  PubMed  Google Scholar 

  16. Humphries RM, Kircher S, Ferrell A, Krause KM, Malherbe R, Hsiung A, Burnham CA. The Continued Value of Disk Diffusion for Assessing Antimicrobial Susceptibility in Clinical Laboratories: Report from the Clinical and Laboratory Standards Institute Methods Development and Standardization Working Group. J Clin Microbiol. 2018;56(8):e00437-18.

  17. Velasco V, Sherwood JS, Rojas-Garcia PP, Logue CM. Multiplex real-time PCR for detection of Staphylococcus aureus, mecA and Panton-valentine Leukocidin (PVL) genes from selective enrichments from animals and retail meat. PLoS One. 2014;9(5):e97617.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  18. Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Li R, Yu C, Li Y, Lam TW, Yiu SM, Kristiansen K, Wang J. SOAP2: an improved ultrafast tool for short read alignment. Bioinformatics. 2009;25(15):1966–7.

    Article  CAS  PubMed  Google Scholar 

  20. Luo R, Liu B, Xie Y, Li Z, Huang W, Yuan J, He G, Chen Y, Pan Q, Liu Y, et al. SOAPdenovo2: an empirically improved memory-efficient short-read de novo assembler. Gigascience. 2012;1(1):18.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Besemer J, Lomsadze A, Borodovsky M. GeneMarkS: a self-training method for prediction of gene starts in microbial genomes. Implications for finding sequence motifs in regulatory regions. Nucleic Acids Res. 2001;29(12):2607–18.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Buchfink B, Xie C, Huson DH. Fast and sensitive protein alignment using DIAMOND. Nat Methods. 2015;12(1):59–60.

    Article  CAS  PubMed  Google Scholar 

  23. Zhao Y, Wu J, Yang J, Sun S, Xiao J, Yu J. PGAP: pan-genomes analysis pipeline. Bioinformatics. 2012;28(3):416–8.

    Article  CAS  PubMed  Google Scholar 

  24. Delcher AL, Salzberg SL, Phillippy AM. Using MUMmer to identify similar regions in large sequence sets. Curr Protoc Bioinformatics. 2003;10(3):1-18.

  25. Klevens RM, Morrison MA, Nadle J, Petit S, Gershman K, Ray S, Harrison LH, Lynfield R, Dumyati G, Townes JM, et al. Invasive methicillin-resistant Staphylococcus aureus infections in the United States. JAMA. 2007;298(15):1763–71.

    Article  CAS  PubMed  Google Scholar 

  26. Udo EE, Pearman JW, Grubb WB. Genetic analysis of community isolates of methicillin-resistant Staphylococcus aureus in Western Australia. J Hosp Infect. 1993;25(2):97–108.

    Article  CAS  PubMed  Google Scholar 

  27. Xiao YH, Giske CG, Wei ZQ, Shen P, Heddini A, Li LJ. Epidemiology and characteristics of antimicrobial resistance in China. Drug Resist Updat. 2011;14(4–5):236–50.

    Article  PubMed  Google Scholar 

  28. Aires de Sousa M, Crisostomo MI, Sanches IS, Wu JS, Fuzhong J, Tomasz A, de Lencastre H. Frequent recovery of a single clonal type of multidrug-resistant Staphylococcus aureus from patients in two hospitals in Taiwan and China. J Clin Microbiol. 2003;41(1):159–63.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Kim ES, Lee HJ, Chung GT, Lee YS, Shin DH, Jung SI, Song KH, Park WB, Kim NJ, Park KU, et al. Molecular characterization of methicillin-resistant Staphylococcus aureus isolates in Korea. J Clin Microbiol. 2011;49(5):1979–82.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Ko KS, Lee JY, Suh JY, Oh WS, Peck KR, Lee NY, Song JH. Distribution of major genotypes among methicillin-resistant Staphylococcus aureus clones in Asian countries. J Clin Microbiol. 2005;43(1):421–6.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Liu Y, Wang H, Du N, Shen E, Chen H, Niu J, Ye H, Chen M. Molecular evidence for spread of two major methicillin-resistant Staphylococcus aureus clones with a unique geographic distribution in Chinese hospitals. Antimicrob Agents Chemother. 2009;53(2):512–8.

    Article  CAS  PubMed  Google Scholar 

  32. Xu BL, Zhang G, Ye HF, Feil EJ, Chen GR, Zhou XM, Zhan XM, Chen SM, Pan WB. Predominance of the Hungarian clone (ST 239-III) among hospital-acquired meticillin-resistant Staphylococcus aureus isolates recovered throughout mainland China. J Hosp Infect. 2009;71(3):245–55.

    Article  CAS  PubMed  Google Scholar 

  33. Xiao M, Wang H, Zhao Y, Mao LL, Brown M, Yu YS, O'Sullivan MV, Kong F, Xu YC. National surveillance of methicillin-resistant Staphylococcus aureus in China highlights a still-evolving epidemiology with 15 novel emerging multilocus sequence types. J Clin Microbiol. 2013;51(11):3638–44.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Harris SR, Feil EJ, Holden MT, Quail MA, Nickerson EK, Chantratita N, Gardete S, Tavares A, Day N, Lindsay JA, et al. Evolution of MRSA during hospital transmission and intercontinental spread. Science. 2010;327(5964):469–74.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Li Y, Cao B, Zhang Y, Zhou J, Yang B, Wang L. Complete genome sequence of Staphylococcus aureus T0131, an ST239-MRSA-SCCmec type III clone isolated in China. J Bacteriol. 2011;193(13):3411–2.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Wang Z, Zhou H, Wang H, Chen H, Leung KK, Tsui S, Ip M. Comparative genomics of methicillin-resistant Staphylococcus aureus ST239: distinct geographical variants in Beijing and Hong Kong. BMC Genomics. 2014;15:529.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  37. Yan X, Wang B, Tao X, Hu Q, Cui Z, Zhang J, Lin Y, You Y, Shi X, Grundmann H. Characterization of Staphylococcus aureus strains associated with food poisoning in Shenzhen, China. Appl Environ Microbiol. 2012;78(18):6637–42.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Li G, Wu S, Luo W, Su Y, Luan Y, Wang X. Staphylococcus aureus ST6-t701 isolates from food-poisoning outbreaks (2006-2013) in Xi’an, China. Foodborne Pathog Dis. 2015;12(3):203–6.

    Article  CAS  PubMed  Google Scholar 

  39. Diep BA, Gill SR, Chang RF, Phan TH, Chen JH, Davidson MG, Lin F, Lin J, Carleton HA, Mongodin EF, et al. Complete genome sequence of USA300, an epidemic clone of community-acquired meticillin-resistant Staphylococcus aureus. Lancet. 2006;367(9512):731–9.

    Article  CAS  PubMed  Google Scholar 

  40. Stegger M, Price LB, Larsen AR, Gillece JD, Waters AE, Skov R, Andersen PS. Genome sequence of Staphylococcus aureus strain 11819-97, an ST80-IV European community-acquired methicillin-resistant isolate. J Bacteriol. 2012;194(6):1625–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Chua KY, Seemann T, Harrison PF, Monagle S, Korman TM, Johnson PD, Coombs GW, Howden BO, Davies JK, Howden BP, et al. The dominant Australian community-acquired methicillin-resistant Staphylococcus aureus clone ST93-IV [2B] is highly virulent and genetically distinct. PLoS One. 2011;6(10):e25887.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Huang TW, Chen FJ, Miu WC, Liao TL, Lin AC, Huang IW, Wu KM, Tsai SF, Chen YT, Lauderdale TL. Complete genome sequence of Staphylococcus aureus M013, a pvl-positive, ST59-SCCmec type V strain isolated in Taiwan. J Bacteriol. 2012;194(5):1256–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Li J, Wang L, Ip M, Sun M, Sun J, Huang G, Wang C, Deng L, Zheng Y, Fu Z, et al. Molecular and clinical characteristics of clonal complex 59 methicillin-resistant Staphylococcus aureus infections in mainland China. PLoS One. 2013;8(8):e70602.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Chen CJ, Unger C, Hoffmann W, Lindsay JA, Huang YC, Gotz F. Characterization and comparison of 2 distinct epidemic community-associated methicillin-resistant Staphylococcus aureus clones of ST59 lineage. PLoS One. 2013;8(9):e63210.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Huh K, Chung DR. Changing epidemiology of community-associated methicillin-resistant Staphylococcus aureus in the Asia-Pacific region. Expert Rev Anti-Infect Ther. 2016;14(11):1007–22.

    Article  CAS  PubMed  Google Scholar 

  46. Feng Y, Chen HL, Chen CJ, Chen CL, Chiu CH. Genome comparisons of two Taiwanese community-associated methicillin-resistant Staphylococcus aureus ST59 clones support the multi-origin theory of CA-MRSA. Infect Genet Evol. 2017;54:60–5.

    Article  CAS  PubMed  Google Scholar 

  47. Qu T, Feng Y, Jiang Y, Zhu P, Wei Z, Chen Y, Otto M, Yu Y. Whole genome analysis of a community-associated methicillin-resistant Staphylococcus aureus ST59 isolate from a case of human sepsis and severe pneumonia in China. PLoS One. 2014;9(2):e89235.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  48. Cheng VCC, Wong SC, Cao H, Chen JHK, So SYC, Wong SCY, Sridhar S, Yuen KY, Ho PL. Whole-genome sequencing data-based modeling for the investigation of an outbreak of community-associated methicillin-resistant Staphylococcus aureus in a neonatal intensive care unit in Hong Kong. Eur J Clin Microbiol Infect Dis. 2019;38(3):563–73.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgments

None.

Funding

This work was supported by National Sci-Tech key project (2017ZX10103010) and open subject of the first people’s hospital of Yunnan province (2019LCZXKF-HX01). The funding body has no role in design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Author information

Authors and Affiliations

Authors

Contributions

X. F., and Y. Z. design of the work; F. L., Z. M., and W. G. do the experiments; X. F., and W. X. analysis the data; F. L. drafting the work. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Yunhui Zhang.

Ethics declarations

Ethics approval and consent to participate

The human sample collection and detection protocols were approved by the Ethics Review Committee [Institutional Review Board (IRB)] of The People’s Hospital of Kunming City and The First People’s Hospital of Yunnan Province. All adult subjects provided informed consent, and a parent or guardian of any child participant provided informed consent on their behalf. The informed consents were verbal for all the participants, because the samples were too large; we couldn’t get all the written ones. In addition, the Ethics Review Committee (IRB) has approved this procedure.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1.

The molecular typing database of patients’ MRSA in Yunnnan province of China.

Additional file 2.

The genotype profiles of patients’ MRSA in Yunnan province of China.

Additional file 3.

The details of reference strains used for phylogenetic analysis in this study.

Additional file 4.

The details of strain information and antibiotic resistant results in this study.

Additional file 5.

The genome features of 60 MRSA in this study.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Liao, F., Mo, Z., Gu, W. et al. 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 (2020). https://0-doi-org.brum.beds.ac.uk/10.1186/s12879-020-4866-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12879-020-4866-6

Keywords