Skip to main content

ORIGINAL RESEARCH article

Front. Microbiol., 25 September 2018
Sec. Infectious Agents and Disease

Molecular Epidemiology of Staphylococcus aureus Bacteremia: Association of Molecular Factors With the Source of Infection

\r\nDafne Prez-Montarelo&#x;Dafne Pérez-Montarelo1Esther Viedma*&#x;Esther Viedma1*Nieves LarrosaNieves Larrosa2Carmen Gmez-GonzlezCarmen Gómez-González3Enrique Ruiz de GopeguiEnrique Ruiz de Gopegui4Irene Muoz-GallegoIrene Muñoz-Gallego1Rafael San JuanRafael San Juan1Nuria Fernndez-HidalgoNuria Fernández-Hidalgo2Benito AlmiranteBenito Almirante2Fernando ChavesFernando Chaves1
  • 1Hospital Universitario 12 de Octubre, Madrid, Spain
  • 2Hospital Universitario Vall d'Hebron, Barcelona, Spain
  • 3Hospital Galdakao-Usansolo, Galdakao, Spain
  • 4Hospital Universitario Son Espases, Palma de Mallorca, Spain

Staphylococcus aureus bacteremia (SAB) is associated with high morbidity and mortality, which varies depending on the source of infection. Nevertheless, the global molecular epidemiology of SAB and its possible association with specific virulence factors remains unclear. Using DNA microarrays, a total of 833 S. aureus strains (785 SAB and 48 colonizing strains) collected in Spain over a period of 15 years (2002–2017) were characterized to determine clonal complex (CC), agr type and repertoire of resistance and virulence genes in order to provide an epidemiological overview of CCs causing bloodstream infection, and to analyze possible associations between virulence genes and the most common sources of bacteremia. The results were also analyzed by acquisition (healthcare-associated [HA] and community-acquired [CA]), methicillin-resistant (MRSA) and methicillin-susceptible (MSSA) strains, and patient age (adults vs. children). Our results revealed high clonal diversity among SAB strains with up to 28 different CCs. The most prevalent CCs were CC5 (30.8%), CC30 (20.3%), CC45 (8.3%), CC8 (8.4%), CC15 (7.5%), and CC22 (5.9%), which together accounted for 80% of all cases. A higher proportion of CC5 was found among HA strains than CA strains (35.6 vs. 20.2%, p < 0.001). CC5 was associated with methicillin resistance (14.7 vs. 79.4%, p < 0.001), whereas CC30, CC45, and CC15 were correlated with MSSA strains (p < 0.001). Pathogen-related molecular markers significantly associated with a specific source of bacteremia included the presence of sea, undisrupted hlb and isaB genes with catheter-related bacteremia; sed, splE, and fib genes with endocarditis; undisrupted hlb with skin and soft tissue infections; and finally, CC5, msrA resistance gene and hla gene with osteoarticular source. Our study suggests an association between S. aureus genotype and place of acquisition, methicillin resistance and sources of bloodstream infection, and provides a valuable starting point for further research insights into intrinsic pathogenic mechanisms involved in the development of SAB.

Introduction

Staphylococcus aureus is an opportunistic pathogen that can potentially cause a wide range of infections. It is a leading cause of bacteremia and represents a significant global health problem (Weiner et al., 2016). S. aureus bacteremia (SAB) is often associated with severe metastatic infections, such as infective endocarditis, septic arthritis and osteomyelitis and complications, such as sepsis and septic shock, which lead to adverse outcomes that are challenging to manage (Shorr et al., 2006; Wyllie, 2006).

The incidence of SAB is difficult to determine and there are major geographical differences that reflect discrepancies in health care systems and infection control practices. In developed countries, the estimated incidence 80–190 cases per 100,000 inhabitants per year (Laupland, 2013; Le Moing et al., 2015). Despite the improvements in SAB management, including greater understanding of this infection and mandatory surveillance implemented in several countries over recent decades, SAB still causes significant morbidity and mortality, with an associated early mortality that appears to have plateaued at approximately 20–30% (van Hal et al., 2012). Certainly, little is known about global SAB epidemiology in terms of the circulating clones causing SAB in different patient subgroups, such as adults and children, or those most commonly found in the community or hospital settings. Because it is becoming progressively more difficult to differentiate between healthcare-associated and community-acquired infections due to changes in the complexity of present health care systems, it is important to identify the specific clones that are traditionally associated with the community but may be entering hospitals and replacing common nosocomial clones, and vice versa. Moreover, it would be especially interesting to study clonality taking into account that bacterial phenotype and genotype have been shown to have a possible influence on infection outcome, since different clones can adopt different strategies to overcome host responses and cause severe pathology (Recker et al., 2017). The overall mortality rate from SAB varies depending on the primary focus of infection (the highest rates occur in patients with infective endocarditis and pulmonary infections, and the lowest in patients with catheter-related infections) and on the complications deriving from SAB. This association makes it necessary to regard SAB not as a single entity, but as a heterogeneous group of infections that can evolve differently and therefore require source-specific management (van Hal et al., 2012). However, the characteristics of the most common clones causing SAB according to source of infection remain unknown. Furthermore, determining the role of particular genetic backgrounds (clonality and virulence) in bloodstream infections caused by S. aureus has become a real challenge due to the diversity, redundancy and host specificity of the virulence factors.

The aim of the present study was to explore the molecular characteristics of S. aureus strains causing bacteremia in order to provide an epidemiological overview of the circulating clones causing bloodstream infection and to analyze the possible association between virulence and the most common sources of bacteremia.

Materials and Methods

Data Collection

A total of 785 strains causing bacteremia with different source and 48 colonization strains collected over a period of 15 years (2002–2017) were analyzed. These strains were obtained from different sources in hospitals geographically distant from each other spread across the territory of Spain (Table 1). Specifically, these strains were identified in 10 different collections: six were single-center studies developed at the Hospital 12 de Octubre in Madrid, and the remaining four corresponded to multi-center studies developed at various Spanish hospitals (Muñoz-Gallego et al., 2017; San-Juan et al., 2017; Fernández-Hidalgo et al., 2018). The main focus and objective of the studies for which these strains were collected was source of staphylococcal bacteremia, mainly endocarditis (N = 214), catheter-related bacteremia (CRB) (N = 212), skin and soft tissue infections (SSTI) (N = 66), and bone and joint infections (N = 100). Eight of these collections corresponded to SAB infections in adults, and two in children (< 15 years of age). The percentage of MRSA strains included in each collection varied. The studies were approved by the ethics committee of the University Hospital 12 de Octubre (Madrid, Spain). It was not considered necessary to obtain written informed consent because the participants were anonymized (IRH-ANT-2013-01).

TABLE 1
www.frontiersin.org

Table 1. Description of the bacteremia-producing strain collections included in the study.

Cases were classified according to acquisition: healthcare-associated (HA) or community-acquired (CA). HA included both nosocomial cases with a positive blood culture obtained from patients who had been hospitalized for 48 h or longer (Garner et al., 1988) and healthcare-associated cases following Friedman et al.'s criteria (Friedman et al., 2002). CA cases were those with a positive blood culture obtained at the time of hospital admission or within 48 h after hospital admission.

Methicillin resistance was defined on the basis of results of microdilution techniques, cefoxitin susceptibility testing and/or the presence of the mecA gene.

Molecular Studies

Blood cultures were processed with an automated blood culture system (BACTEC 9240, Becton Dickinson Microbiological System, USA). Automatic microdilution techniques were used for identification and susceptibility testing of isolates. Bacterial DNA was extracted using commercial extraction kits (Qiagen, Germany) according to the manufacturer's recommendations. DNA microarrays (Alere, Germany; Monecke et al., 2008) covering 334 target sequences and approximately 187 different genes that included species-specific markers, antimicrobial resistance genes, exotoxins, genes encoding microbial surface components recognizing adhesive matrix molecules (MCSCRAMMs), capsule genes, clonal complexes (CC) and agr group typing markers were run on the whole collection of strains. Those cases with ambiguous array results were considered as missing values for further analysis.

Only genes found with a frequency of between 5 and 95% in the whole collection were considered for statistical analysis.

Statistical Analysis

Categorical variables were compared using the chi-squared or Fisher's exact test, as appropriate. Significance levels of DNA microarray results were corrected using the Bonferroni correction for multiple tests. Pairwise comparisons of the main CCs, agr types and virulence genes were performed with source of bacteremia. Potential associations were investigated by univariate and multivariate logistic regression, in which CCs, agr types and virulence factors were considered as independent dichotomous variables, and source of bacteremia as the dependent variable. For multivariate analysis, variables with a p-value <0.1 in the univariate analysis were included in a backward stepwise algorithm. All statistical tests were two-tailed and a p-value of <0.05 was considered statistically significant. Analyses were performed using the SPSS statistical package, version 21.0 (SPSS Inc., Chicago, IL).

Results

A total of 785 S. aureus strains causing bacteremia and 48 S. aureus colonizing strains were characterized in this study with the DNA microarray. Of a total of 187 genes included in the array, 67 genes were excluded from further analysis. Twenty one genes because they were found in almost all of the strains (>95%): sarA(99.5%), saeS(98.4%), vraS(99.7%), lukF(96.7%), hl(97.5%), hld(99.3%), sspA(99.8%), sspB(100%), sspP(99.8%), icaA(99.6%), icaC(98.4%), icaD(99.5%), clfA(99.6%), clfB(99.4%), ebpS(99.4%), eno(99.5%), fnbA(98.1%), map(97.6%), sdrC(99.3%) and isdA(99.6%); and 47 genes because they were absent from almost all of the strains (< 5%): ermA(4.7%), ermB(0.3%), lnuA(0.7%), mefA(0%), vatA(0%), vatB(0%), vgaA(2.6%), vgaB(0%), aacA-aphD(3.9%), dfrS1(2.5%), fusB(0%), fusC(0.8%), tetK(2.2%), tetM(2.1%), cat(0.8%), cfr(0.1%), fexA(0.2%), qacA(0.2%), qacC(2.2%), vanA(0%), vanB(0%), vanZ(0%), seb(3.6%), see(0%), seh(4.6%), sej(3.6%), sek(3.4%), seq(3.1%), ser(3.4%), pvl(1.8%), etA(2.5%), etB(0.8%), etD(3.1%), edinA(0.6%), edinB(3.3%), edinC(0.7%), ACME cluster(0.5%), arcA-SCC(0.2%), arcB-SCC(0.1%), arc-SCC(0.1%), arcD-SCC(3.6%), cap1(0.7%), and bap(0.1%). Nine additional genes were discarded due to an unacceptable number of missing values (>30%): mecC(55%), mecR(56.2%), merA(53.3%), merB(39.5%), fosBplasmid(62.8%), lukM(52.8%), lukY(51.3%), setC-selX(52.3%), setB3(51.5%), setB2(51.7%), setB1(52.6%), fnbB(51.3%), and sdrD(51.6%).

The main CCs detected in this study were: CC5 (30.8%), CC30 (20.3%), CC45 (8.3%), CC8 (8.4%), CC15 (7.5%) and CC22 (5.9%). In addition, up to 22 minor CCs were also detected: CC398 (2.4%), CC121 (2.4%), CC25 (2.2%), CC9 (1.7%), CC97 (1.6%), CC6 (1.3%), CC1 (1.1%), CC7 (1.1%), CC188 (1.0%), CC101, CC10, CC49, CC59, CC509, CC20, CC12, CC75, CC96, CC395, CC522, CC707 and CC1021 (<1%; Table 2). agrII type was the most common (41.7%), followed by agrI (33.7%) and agrIII (22.2%).

TABLE 2
www.frontiersin.org

Table 2. Distribution of S. aureus clonal complexes according to place of acquisition, methicillin resistance and age of population.

Distribution of S. aureus Strains by CC and agr Type According to Acquisition, Methicillin Resistance and Age of Population

Healthcare-Associated vs. Community-Acquired

In our collection, there was a higher proportion of HA strains (68.4%) compared to CA strains, which accounted for 31.6%. Healthcare-associated strains were assigned to 26 different CCs, with CC5 (35.6%) being the most common, followed by CC30 (18.6%). The CC diversity was slightly lower among CA strains, which were assigned to 23 CCs, with CC30 (23.4%) and CC5 (20.2%) being the ones most commonly found (Table 2). Although most clones circulated in the healthcare and community settings, a higher proportion of CC5 was found among HA than among CA strains (35.6 vs. 20.2%, p < 0.001; Table 2).

The distribution of agr types is presented in Table 3. Note that while agrI was the main agr type among CA strains (31.0% vs. 39.2%, p: 0.030), agrII was associated with HA acquisition, (46.6 vs. 31.8%, p < 0.001).

TABLE 3
www.frontiersin.org

Table 3. Distribution of agr types according to place of acquisition, methicillin resistance and age of population.

Methicillin Resistance

Methicillin resistance was observed in 24.7% of the strains. Eleven different CCs were detected in this group, with CC5 being the main clone (79.4%). By contrast, among the MSSA strains, all the CCs detected in our collection were represented (N = 28) and showed a higher clonal diversity than MRSA strains. The major clones detected among MSSA strains were CC30 (26.5%), CC5 (14.5%), CC45 (10.5%) and CC15 (10.0%). A comparison of the two groups revealed that only CC5 was associated with methicillin resistance (14.7 vs. 79.4%, p < 0.001), whereas CC30, CC45 and CC15 were correlated with MSSA strains (p < 0.001; Table 2).

With respect to agr type, agrII type was observed to be significantly associated with MRSA strains (29.1 vs. 80.6%, p < 0.001), whereas, agrI was associated with MSSA strains (38.6 vs. 18.4%, p < 0.001; Table 3).

Adult vs. Child Population

While 28 CCs detected in this study were represented among strains isolated from adults, lower CC diversity was detected in children, with only 19 CCs identified. A comparison of clonality in the adult and children populations revealed that while CC5 was associated with strains from adults (34.2 vs. 20.0%, p < 0.001), CC30 was significantly related to strains from the child population (18.2 vs. 27.2%, p < 0.008; Table 2).

The distribution of agr types by patient age is shown in Table 3. This analysis showed a significant association between agrII and strains from adults (44.7 vs. 32.7%, p < 0.001), while agrIII was associated with the child population (19.9 vs. 29.1%, p < 0.008).

Clonal Complex Diversity, agr Type and Virulence Genes Among S. aureus Strains From Different Sources of Bacteremia

The main objective was to explore the distribution of CCs and virulence genes according to source of bacteremia. A collection of S. aureus strains from healthy carriers was also added to the analysis in order to evaluate potential differences between colonizing and bacteremic strains.

Remarkably, CC5 and agrII predominated in SAB from osteoarticular infections (Tables 4A,B). When we focused on antibiotic resistance genes, significant differences were identified in strains from different sources of infection for the mecA, msrA, aadD, aphA3, and sat genes (Tables 4A,B). Our results seemed to indicate a higher proportion of these resistance genes among osteoarticular infections compared with other bacteremia sources. In general, significant differences for source of bacteremia and colonization were also detected in virulence genes, such as sea, sed, hla, undisrupted hlb, splE, cna, fib, and isaB among others.

TABLE 4A
www.frontiersin.org

Table 4A. Frequency of clonal complexes, agr type, resistance and virulence genes according to different sources of infection.

TABLE 4B
www.frontiersin.org

Table 4B. Pairwise comparison of clonal complexes, agr types and virulence genes according to source of infection.

A comparison of colonizing vs. SAB strains revealed that the sed gene was most frequently found in colonizing strains (20.8 vs. 6.7%, p < 0.001; Table 4B). This gene also had greater representation among endocarditis strains (11.4 vs. 6.1%, p: 0.031). In addition, the undisrupted hlb gene mostly presented in CRB strains (41.0 vs. 20.5%, p < 0.001; Table 4B). The presence of MSCRAMM genes, such as cna, fib, vwb and isaB also varied significantly according to source of SAB. While cna and vwb genes had greater representation among strains from an endocarditis source (48.5 vs. 37.8%, p: 0.012 and 99.1 vs. 94.4%, p: 0.003, respectively), there was significant detection of fib genes (77.3 vs. 62.5%, p: 0.025) in SSTI sources. Finally, the isaB gene was more commonly found among CRB strains (88.5 vs. 78.1%, p: 0.002) and in the SSTI bacteremia group (100 vs. 79.4%, p < 0.001; Table 4B).

Various regression models were performed in order to measure the role of pathogen-related molecular markers (CC, agr type and virulence genes) adjusted for different sources of bacteremia and colonization (Table 5). All adjusted models in multivariate analysis showed that these variables were the presence of agrIV type and sed gene for colonizing strains; the presence of sea, undisrupted hlb and isaB genes for CRB; sed, splE and fib genes for an endocarditis source; undisrupted hlb for the SSTI group; and finally, CC5, msrA resistance gene and hla gene with respect to bacteremia from an osteoarticular source (Table 5).

TABLE 5
www.frontiersin.org

Table 5. Multivariate analysis according to different sources of SAB and colonization.

Discussion

The present study describes and gives a global epidemiological overview of the molecular epidemiology of a large collection of S. aureus strains focused on bloodstream infections in Spain over a 15-years period. In this scenario, our study provides important findings regarding the distribution of clonality and virulence genes and their association with specific sources of SAB.

Our results revealed high clonal diversity among SAB strains, although the most prevalent CCs were CC5, CC30, CC45, CC8, CC15, and CC22, which together represented 80% of all cases. Additionally, substantial differences were found between strains causing MRSA and MSSA bacteremia, which indicated that MSSA strains were much more genetically diverse than their MRSA counterparts, which is consistent with studies developed in Europe (Aamot et al., 2012; Grundmann et al., 2014) and the USA (Miko et al., 2013; Park et al., 2017). The most common clone among MSSA strains was CC30, followed by CC45 and CC15, whereas among MRSA strains, there was a significant representation of CC5 in more than 75% of strains. Similar results have been reported in Latin America (Arias et al., 2017) and Germany (Schaumburg et al., 2012), where the CC5-MRSA clone was the most prevalent in the setting of bloodstream infections. Furthermore, in our collection, CC5 was found to be significantly associated with HA acquisition and the adult population, a finding which lends support to the interest aimed at investigating the pathogenic and molecular characteristics of CC5 and those factors that enhance its spread.

Several studies have suggested that while the agrI type is the most common one among clinical isolates (van Leeuwen et al., 2000; Moore and Lindsay, 2001), others (Sakoulas et al., 2003) have determined that more than half of clinical MRSA bloodstream isolates belong to agr group II. In our collection, agrII was also associated with MRSA, which may explain the higher percentage of agrII found in the nosocomial setting and among adults, in whom the prevalence of MRSA was higher. By contrast, agrI was related to MSSA strains and CA acquisition. Interestingly, a statistically significant association was also found between agrII and agrIII and adult and child populations, respectively. These associations are probably due to the correlation between agr type and CC, since CC5 (agrII) was the majority clone among adults and CC30 (agrIII) among children.

To date, different studies have explored the association between bacterial genotype, especially S. aureus virulence genes, and various clinical syndromes (Gillet et al., 2002; Jarraud et al., 2002; Peacock et al., 2002). This study focuses specifically on bacteremia. Our collection included S. aureus strains from the most common primary clinical sources of infection: CRB, SSTI, osteoarticular infection and endocarditis, as well as nasal carriage strains. We found no major differences between colonizing and bacteremia-producing strains of S. aureus, which supports the fact that most strains of S. aureus are capable of causing bacteremia. Nevertheless, and in accordance with other studies (Fowler et al., 2007; Giulieri et al., 2016), we identified specific clonal backgrounds and various molecular markers that have been associated with bloodstream infections and certain sources of bacteremia in particular. In this regard, our findings showed that CC5 in addition to hla and msrA genes were more frequently present in strains causing osteoarticular bacteremia. The association of the hla gene, present in most S. aureus strains, with different types of infection has already been reported (Stulik et al., 2014; Sharma-Kuinkel et al., 2015). Further studies are needed to elucidate the role of this important virulence factor in the pathogenesis of bacteremia. With respect to the adhesin genes (MSCRAMMs), which play an essential role in the pathogenesis of intravascular, osteoarticular and device-associated S. aureus infections (Foster et al., 2014), our study revealed an association between the fib and isaB genes and endocarditis and CRB sources, respectively. Other adhesins like clfA/B, fnbA/B, and cna and their linkage with bacteremia, endocarditis and CRB, have also been reported (Giulieri et al., 2016; San-Juan et al., 2017).

Another finding of note in this study was the presence of the undisrupted β-hemolysin (undisrupted hlb) which was significantly related to sources, such as CRB and SSTI. Different studies have demonstrated its contribution to SSTI (Hedström and Malmqvist, 1982; Lebughe et al., 2017) and biofilm-related infections (Salgado-Pabón et al., 2014). Although β-toxin is encoded in S. aureus, most strains are reported not to secrete β-toxin because the bacteriophage (ϕSa3) inserts into the hlb gene (Winkler et al., 1965; Coleman et al., 1991), inactivating it in the majority of S. aureus strains recovered from humans. Moreover, the ϕSa3 bacteriophage encodes the immune evasion cluster (IEC) sak-chip-scn (Coleman et al., 1989; de Haas et al., 2004). Coinciding with other studies (Pantucek et al., 2004; Van Wamel et al., 2006), these genes were relatively abundant in our collection, ranging between 73% (sak) and 87% (scn). Interestingly, the absence of the intact hlb gene (or which amounts to the same thing, the presence of hlb truncated by the IEC-carrying ϕSa3 phage) was significantly associated with an osteoarticular source. This intriguing association should be investigated further since other studies have reported the association between these phage-integrated genes and less severe staphylococcal infections (Jin et al., 2003).

This study presents several limitations that should be mentioned. First, the heterogeneity and non-continuity of the SAB collection (geographical origin, time points and hosts) precluded us from adjusting for these variables in multivariate analysis. Moreover, the proportion of colonization strains was small in comparison with the number of SAB strains. The results therefore should be interpreted with caution. At the same time, our study includes a large number of S. aureus strains causing bacteremia, with relevant information on place of acquisition, methicillin resistance and source of infection. Second, the lack of clinical data regarding the outcome of the bacteremic episodes makes it impossible to make inferences about the prognostic importance of the molecular factors. Other studies evaluating associations between bacterial genotype and virulence have led to conflicting results (Day et al., 2001, 2002; Feil et al., 2003; Melles et al., 2004), due in part to the heterogeneous nature of the S. aureus infections included, as well as the absence of a large, well-characterized collection of isolates. Third, our study methodology was based on DNA microarrays, which should be noted in the case of the hla gene. Despite the fact that hla is present in virtually all S. aureus strains, some studies, such as Sharma-kuinkel et al. (JCM) have reported up to 12 different variants of hla. In our study, the hla gene was detected in 92.2% of strains. We think that the low frequency of this gene observed in our collection may have been due to the DNA microarray technology, which may underestimate the presence of certain minority hla variants due to lack of sensitivity. Whole genome sequencing may be a more effective genotypic characterization approach for detecting different genetic variants that may not be detected by hybridization procedures, although previous studies have shown good agreement between the genotypic results obtained using a DNA array-based methodology and those using high-throughput sequencing (Strauß et al., 2016). Finally, we did not perform gene expression studies, which would be key to determining whether a particular gene or set of genes was responsible for the specific pathogenic behavior observed in SAB from particular clinical sources. Nevertheless, our findings offer a valuable starting point for further research insights into intrinsic pathogenic mechanisms involved in the development of SAB.

In conclusion, the current study suggests a potential association between S. aureus genotype and acquisition, methicillin resistance and bloodstream infection sources. The results of this study reinforce the view that SAB continues to represent a major clinical challenge. Thus, a better understanding of S. aureus epidemiology and pathogenesis is crucial to the detection of prognostic biomarkers as well as to the development of potential therapeutic targets aimed at improving patient outcomes.

Author Contributions

DP-M, EV, and FC conceived and designed the experiments. EV, CG-G, ERG, NL, NF-H, RS, and IM-G collected the isolates. Funding was obtained by FC and BA. Experiments were performed by EV, CG-G, and IM-G. The data were analyzed by EV and DP-M. EV and DP-M prepared the manuscript draft. All authors agreed to be accountable for all aspects of the work. EV, DP-M, and FC contributed in giving final approval of the version to be published. All authors reviewed and approved the final manuscript.

Funding

This work was supported by the Health Research Fund, Department of Health, Spain; Agency for Health Technology Assessment and Research (PI15/02013 and PI15/02125) and Instituto de Salud Carlos III, Subdirección General de Redes y Centros de Investigación Cooperativa, Ministerio de Economía y Competitividad, Spanish Network for Research in Infectious Diseases (REIPI RD16/0016-0002; 0003; 0004) and cofunded by the European Regional Development Fund (FEDER).

Conflict of Interest Statement

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

Acknowledgments

We thank Mercedes Murcia for excellent technical assistance, Jaime Lora-Tamayo for statistical analysis support and Janet Dawson for language support. Moreover, we thank all clinical and microbiological researchers that participated in studies that led to the collection of the Staphylococcus aureus strains presented in this study.

References

Aamot, H. V., Blomfeldt, A., and Eskesen, A. N. (2012). Genotyping of 353 Staphylococcus aureus bloodstream isolates collected between 2004 and 2009 at a Norwegian university hospital and potential associations with clinical parameters. J. Clin. Microbiol. 50, 3111–3114. doi: 10.1128/JCM.01352-12

PubMed Abstract | CrossRef Full Text | Google Scholar

Arias, C. A., Reyes, J., Carvajal, L. P., Rincon, S., Diaz, L., Panesso, D., et al. (2017). A prospective cohort multicenter study of molecular epidemiology and phylogenomics of Staphylococcus aureus bacteremia in nine latin american countries. Antimicrob. Agents Chemother. 61, e00816–e00817. doi: 10.1128/AAC.00816-17

PubMed Abstract | CrossRef Full Text | Google Scholar

Coleman, D., Knights, J., Russell, R., Shanley, D., Birkbeck, T. H., Dougan, G., et al. (1991). Insertional inactivation of the Staphylococcus aureus beta-toxin by bacteriophage phi 13 occurs by site- and orientation-specific integration of the phi 13 genome. Mol. Microbiol. 5, 933–939. doi: 10.1111/j.1365-2958.1991.tb00768.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Coleman, D. C., Sullivan, D. J., Russell, R. J., Arbuthnott, J. P., Carey, B. F., and Pomeroy, H. M. (1989). Staphylococcus aureus bacteriophages mediating the simultaneous lysogenic conversion of beta-lysin, staphylokinase and enterotoxin A: molecular mechanism of triple conversion. J. Gen. Microbiol. 135, 1679–1697.

PubMed Abstract | Google Scholar

Day, N. P., Moore, C. E., Enright, M. C., Berendt, A. R., Smith, J. M., Murphy, M. F., et al. (2001). A link between virulence and ecological abundance in natural populations of Staphylococcus aureus. Science 292, 114–116. doi: 10.1126/science.1056495

PubMed Abstract | CrossRef Full Text | Google Scholar

Day, N. P. J., Moore, C. E., Enright, M. C., Berendt, A. P., Smith, J. M., Murphy, M. F., et al. (2002). Retraction. Science 295:971. doi: 10.1126/science.295.5557.971b

PubMed Abstract | CrossRef Full Text | Google Scholar

de Haas, C. J., Veldkamp, K. E., Peschel, A., Weerkamp, F., Van Wamel, W. J., Heezius, E. C., et al. (2004). Chemotaxis inhibitory protein of Staphylococcus aureus, a bacterial antiinflammatory agent. J. Exp. Med. 199, 687–695. doi: 10.1084/jem.20031636

PubMed Abstract | CrossRef Full Text | Google Scholar

Feil, E. J., Cooper, J. E., Grundmann, H., Robinson, D. A., Enright, M. C., Berendt, T., et al. (2003). How clonal is Staphylococcus aureus? J. Bacteriol. 185, 3307–3316. doi: 10.1128/JB.185.11.3307-3316.2003

PubMed Abstract | CrossRef Full Text | Google Scholar

Fernández-Hidalgo, N., Ribera, A., Larrosa, M. N., Viedma, E., Origüen, J., de Alarcón, A., et al. (2018). Impact of Staphylococcus aureus phenotype and genotype on the clinical characteristics and outcome of infective endocarditis. A multicenter, longitudinal, prospective, observational study. Clin. Microbiol. Infect. 24, 985–991. doi: 10.1016/j.cmi.2017.12.002

CrossRef Full Text | Google Scholar

Foster, T. J., Geoghegan, J. A., Ganesh, V. K., and Höök M. (2014). Adhesion, invasion and evasion: the many functions of the surface proteins of Staphylococcus aureus. Nat. Rev. Microbiol. 12, 49–62. doi: 10.1038/nrmicro3161

PubMed Abstract | CrossRef Full Text | Google Scholar

Fowler, V. G., Nelson, C. L., McIntyre, L. M., Kreiswirth, B. N., Monk, A., Archer, G. L., et al. (2007). Potential associations between hematogenous complications and bacterial genotype in Staphylococcus aureus infection. J. Infect. Dis. 196, 738–747. doi: 10.1086/520088

PubMed Abstract | CrossRef Full Text | Google Scholar

Friedman, N. D., Kaye, K. S., Stout, J. E., McGarry, S. A., Trivette, S. L., Briggs, J. P., et al. (2002). Health care–associated bloodstream infections in adults: a reason to change the accepted definition of community-acquired infections. Ann. Intern. Med. 137, 791–797. doi: 10.7326/0003-4819-137-10-200211190-00007

PubMed Abstract | CrossRef Full Text | Google Scholar

Garner, J. S., Jarvis, W. R., Emori, T. G., Horan, T. C., and Hughes, J. M. (1988). CDC definitions for nosocomial infections, 1988. Am. J. Infect. Control 16, 128–140. doi: 10.1016/0196-6553(88)90053-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Gillet, Y., Issartel, B., Vanhems, P., Fournet, J. C., Lina, G., Bes, M., et al. (2002). Association between Staphylococcus aureus strains carrying gene for Panton-Valentine leukocidin and highly lethal necrotising pneumonia in young immunocompetent patients. Lancet 359, 753–759. doi: 10.1016/S0140-6736(02)07877-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Giulieri, S. G., Holmes, N. E., Stinear, T. P., and Howden, B. P. (2016). Use of bacterial whole-genome sequencing to understand and improve the management of invasive Staphylococcus aureus infections. Expert Rev. Anti Infect. Ther. 14, 1023–1036. doi: 10.1080/14787210.2016.1233815

PubMed Abstract | CrossRef Full Text | Google Scholar

Grundmann, H., Schouls, L. M., Aanensen, D. M., Pluister, G. N., Tami, A., Chlebowicz, M., et al. (2014). The dynamic changes of dominant clones of Staphylococcus aureus causing bloodstream infections in the European region: results of a second structured survey. Euro Surveill. 19:20987.

PubMed Abstract | Google Scholar

Hedström, S. A., and Malmqvist, T. (1982). Sphingomyelinase activity of Staphylococcus aureus strains from recurrent furunculosis and other infections. Acta Pathol. Microbiol. Immunol. Scand. B 90, 217–220. doi: 10.1111/j.1699-0463.1982.tb00108.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Jarraud, S., Mougel, C., Thioulouse, J., Lina, G., Meugnier, H., Forey, F., et al. (2002). Relationships between Staphylococcus aureus genetic background, virulence factors, agr groups (alleles), and human disease. Infect. Immun. 70, 631–641. doi: 10.1128/IAI.70.2.631-641.2002

PubMed Abstract | CrossRef Full Text | Google Scholar

Jin, T., Bokarewa, M., McIntyre, L., Tarkowski, A., Corey, G. R., Reller, L. B., et al. (2003). Fatal outcome of bacteraemic patients caused by infection with staphylokinase-deficient Staphylococcus aureus strains. J. Med. Microbiol. 52, 919–923. doi: 10.1099/jmm.0.05145-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Laupland, K. B. (2013). Incidence of bloodstream infection: a review of population-based studies. Clin. Microbiol. Infect. 19, 492–500. doi: 10.1111/1469-0691.12144

PubMed Abstract | CrossRef Full Text | Google Scholar

Le Moing, V., Alla, F., Doco-Lecompte, T., Delahaye, F., Piroth, L., Chirouze, C., et al. (2015). Staphylococcus aureus bloodstream infection and endocarditis–a prospective cohort study. PLoS ONE 10:e0127385. doi: 10.1371/journal.pone.0127385

CrossRef Full Text | Google Scholar

Lebughe, M., Phaku, P., Niemann, S., Mumba, D., Peters, G., Muyembe-Tamfum, J. J., et al. (2017). The impact of the Staphylococcus aureus virulome on infection in a developing country: a cohort study. Front. Microbiol. 8:1662. doi: 10.3389/fmicb.2017.01662

PubMed Abstract | CrossRef Full Text | Google Scholar

López-aguilera, S., Go, M., Barrado, L., González-rodríguez-salinas, M. C., Otero, J. R., and Chaves, F. (2013). Colonización nasal por Staphylococcus aureus en estudiantes de medicina: importancia en la transmisión hospitalaria. Enfermedades Infecciosas Microbiol. Clín. 31, 500–505. doi: 10.1016/j.eimc.2012.12.005

CrossRef Full Text | Google Scholar

Melles, D. C., Gorkink, R. F. J., Boelens, H. A. M., Snijders, S. V., Peeters, J. K., Moorhouse, M. J., et al. (2004). Natural population dynamics and expansion of pathogenic clones of Staphylococcus aureus. J. Clin. Invest. 114, 1732–1740. doi: 10.1172/JCI200423083

PubMed Abstract | CrossRef Full Text | Google Scholar

Miko, B. A., Hafer, C. A., Lee, C. J., Sullivan, S. B., Hackel, M. A., Johnson, B. M., et al. (2013). Molecular characterization of methicillin-susceptible Staphylococcus aureus clinical isolates in the united states, 2004 to 2010. J. Clin. Microbiol. 51, 874–879. doi: 10.1128/JCM.00923-12

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Moore, P. C., and Lindsay, J. A. (2001). Genetic variation among hospital isolates of methicillin-sensitive Staphylococcus aureus: evidence for horizontal transfer of virulence genes. J. Clin. Microbiol. 39, 2760–2767. doi: 10.1128/JCM.39.8.2760-2767.2001

PubMed Abstract | CrossRef Full Text | Google Scholar

Muñoz-Gallego, I., Lora-Tamayo, J., Pérez-Montarelo, D., Bra-as, P., Viedma, E., and Chaves, F. (2017). Influence of molecular characteristics in the prognosis of methicillin-resistant Staphylococcus aureus prosthetic joint infections: beyond the species and the antibiogram. Infection 45, 533–537. doi: 10.1007/s15010-017-1011-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Pantucek, R., Doskar, J., Ruzicková, V., Kaspárek, P., Orácová, E., Kvardová, V., et al. (2004). Identification of bacteriophage types and their carriage in Staphylococcus aureus. Arch. Virol. 149, 1689–1703. doi: 10.1007/s00705-004-0335-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Park, K. H., Greenwood-Quaintance, K. E., Uhl, J. R., Cunningham, S. A., Chia, N., Jeraldo, P. R., et al. (2017). Molecular epidemiology of Staphylococcus aureus bacteremia in a single large Minnesota medical center in 2015 as assessed using MLST, core genome MLST and spa typing. PLoS ONE 12:e0179003. doi: 10.1371/journal.pone.0179003

PubMed Abstract | CrossRef Full Text | Google Scholar

Peacock, S. J., Moore, C. E., Justice, A., Kantzanou, M., Story, L., Mackie, K., et al. (2002). Virulent combinations of adhesin and toxin genes in natural populations of Staphylococcus aureus. Infect. Immun. 70, 4987–4996. doi: 10.1128/IAI.70.9.4987-4996.2002

PubMed Abstract | CrossRef Full Text | Google Scholar

Recker, M., Laabei, M., Toleman, M. S., Reuter, S., Saunderson, R. B., Blane, B., et al. (2017). Clonal differences in Staphylococcus aureus bacteraemia-associated mortality. Nat. Microbiol. 2, 1381–1388. doi: 10.1038/s41564-017-0001-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Sakoulas, G., Eliopoulos, G. M., Moellering, R. C., Novick, R. P., Venkataraman, L., Wennersten, C., et al. (2003). Staphylococcus aureus accessory gene regulator (agr) group II: is there a relationship to the development of intermediate-level glycopeptide resistance? J. Infect. Dis. 187, 929–938. doi: 10.1086/368128

PubMed Abstract | CrossRef Full Text | Google Scholar

Salgado-Pabón, W., Herrera, A., Vu, B. G., Stach, C. S., Merriman, J. A., Spaulding, A. R., et al. (2014). Staphylococcus aureusβ-toxin production is common in strains with theβ-toxin gene inactivated by bacteriophage. J. Infect. Dis. 210, 784–792. doi: 10.1093/infdis/jiu146

CrossRef Full Text | Google Scholar

San-Juan, R., Pérez-Montarelo, D., Viedma, E., Lalueza, A., Fortún, J., Loza, E., et al. (2017). Pathogen-related factors affecting outcome of catheter-related bacteremia due to methicillin-susceptible Staphylococcus aureus in a Spanish multicenter study. Eur. J. Clin. Microbiol. Infect. Dis. 36, 1757–1765. doi: 10.1007/s10096-017-2989-5

PubMed Abstract | CrossRef Full Text | Google Scholar

San-Juan, R., Viedma, E., Chaves, F., Lalueza, A., Fortún, J., Loza, E., et al. (2016). High MICs for vancomycin and daptomycin and complicated catheter-related bloodstream infections with methicillin-sensitive Staphylococcus aureus. Emerg. Infect. Dis. 22, 1057–1066. doi: 10.3201/eid2206.151709

PubMed Abstract | CrossRef Full Text | Google Scholar

Schaumburg, F., Köck, R., Mellmann, A., Richter, L., Hasenberg, F., Kriegeskorte, A., et al. (2012). Population dynamics among methicillin-resistant Staphylococcus aureus isolates in Germany during a 6-year period. J. Clin. Microbiol. 50, 3186–3192. doi: 10.1128/JCM.01174-12

PubMed Abstract | CrossRef Full Text | Google Scholar

Sharma-Kuinkel, B. K., Wu, Y., Tabor, D. E., Mok, H., Sellman, B. R., Jenkins, A., et al. (2015). Characterization of alpha-toxin hla gene variants, alpha-toxin expression levels, and levels of antibody to alpha-toxin in hemodialysis and postsurgical patients with Staphylococcus aureus. J. Clin. Microbiol. 53, 227–236. doi: 10.1128/JCM.02023-14

PubMed Abstract | CrossRef Full Text | Google Scholar

Shorr, A. F., Tabak, Y. P., Killian, A. D., Gupta, V., Liu, L. Z., and Kollef, M. H. (2006). Healthcare-associated bloodstream infection: a distinct entity? Insights from a large U.S. database. Crit. Care Med. 34, 2588–2595. doi: 10.1097/01.CCM.0000239121.09533.09

PubMed Abstract | CrossRef Full Text | Google Scholar

Strauß, L., Ruffing, U., Abdulla, S., Alabi, A., Akulenko, R., Garrine, M., et al. (2016). Detecting Staphylococcus aureus virulence and resistance genes: a comparison of whole-genome sequencing and DNA microarray technology. J. Clin. Microbiol. 54, 1008–1016. doi: 10.1128/JCM.03022-15

PubMed Abstract | CrossRef Full Text | Google Scholar

Stulik, L., Malafa, S., Hudcova, J., Rouha, H., Henics, B. Z., Craven, D. E., et al. (2014). α-hemolysin activity of methicillin-susceptible Staphylococcus aureus predicts ventilator-associated pneumonia. Am. J. Respir. Crit. Care Med. 190, 1139–1148. doi: 10.1164/rccm.201406-1012OC

PubMed Abstract | CrossRef Full Text | Google Scholar

van Hal, S. J., Jensen, S. O., Vaska, V. L., Espedido, B. A., Paterson, D. L., and Gosbell, I. B. (2012). Predictors of mortality in Staphylococcus aureus bacteremia. Clin. Microbiol. Rev. 25, 362–386. doi: 10.1128/CMR.05022-11

PubMed Abstract | CrossRef Full Text | Google Scholar

van Leeuwen, W., van Nieuwenhuizen, W., Gijzen, C., Verbrugh, H., and van Belkum, A. (2000). Population studies of methicillin-resistant and -sensitive Staphylococcus aureus strains reveal a lack of variability in the agrD gene, encoding a staphylococcal autoinducer peptide. J. Bacteriol. 182, 5721–9. doi: 10.1128/JB.182.20.5721-5729.2000

PubMed Abstract | CrossRef Full Text | Google Scholar

Van Wamel, W. J. B., Rooijakkers, S. H. M., Ruyken, M., Van Kessel, K. P. M., and Van Strijp, J. A. G. (2006). The innate immune modulators staphylococcal complement inhibitor and chemotaxis inhibitory protein of Staphylococcus aureus are located on β-hemolysin-converting bacteriophages. J. Bacteriol. 188, 1310–1315. doi: 10.1128/JB.188.4.1310-1315.2006

PubMed Abstract | CrossRef Full Text | Google Scholar

Viedma, E., Sanz, F., Orellana, M. A., San Juan, R., Aguado, J. M., Otero, J. R., et al. (2014). Relationship between agr dysfunction and reduced vancomycin susceptibility in methicillin-susceptible Staphylococcus aureus causing bacteraemia. J. Antimicrob. Chemother. 69, 51–58. doi: 10.1093/jac/dkt337

PubMed Abstract | CrossRef Full Text | Google Scholar

Weiner, L. M., Webb, A. K., Limbago, B., Dudeck, M. A., Patel, J., Kallen, A. J., et al. (2016). Antimicrobial-resistant pathogens associated with healthcare-associated infections: summary of data reported to the national healthcare safety network at the centers for disease control and prevention, 2011–2014. Infect. Control Hosp. Epidemiol. 37, 1288–1301. doi: 10.1017/ice.2016.174

PubMed Abstract | CrossRef Full Text | Google Scholar

Winkler, K. C., de Waart, J., Grootsen, C., Zegers, B. J. M., Tellier, N. F., and Vertregt, C. D. (1965). Lysogenic conversion of staphylococci to loss of beta-toxin. J. Gen. Microbiol. 39, 321–333. doi: 10.1099/00221287-39-3-321

PubMed Abstract | CrossRef Full Text | Google Scholar

Wyllie, D. H. (2006). Mortality after Staphylococcus aureus bacteraemia in two hospitals in Oxfordshire, 1997–2003: cohort study. BMJ 333, 281–280. doi: 10.1136/bmj.38834.421713.2F

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: Staphylococcus aureus bacteremia, bacteremia source, molecular epidemiology, clonal complex, virulence factors

Citation: Pérez-Montarelo D, Viedma E, Larrosa N, Gómez-González C, Ruiz de Gopegui E, Muñoz-Gallego I, San Juan R, Fernández-Hidalgo N, Almirante B and Chaves F (2018) Molecular Epidemiology of Staphylococcus aureus Bacteremia: Association of Molecular Factors With the Source of Infection. Front. Microbiol. 9:2210. doi: 10.3389/fmicb.2018.02210

Received: 22 June 2018; Accepted: 29 August 2018;
Published: 25 September 2018.

Edited by:

Jorge Blanco, Universidade de Santiago de Compostela, Spain

Reviewed by:

Frieder Schaumburg, Universitätsklinikum Münster, Germany
María Guembe, Hospital General Universitario Gregorio Marañón, Spain
QingZhong Liu, Shanghai Jiao Tong University, China

Copyright © 2018 Pérez-Montarelo, Viedma, Larrosa, Gómez-González, Ruiz de Gopegui, Muñoz-Gallego, San Juan, Fernández-Hidalgo, Almirante and Chaves. 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: Esther Viedma, ester.viedma@salud.madrid.org

These authors have contributed equally to this work

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.