We’re sorry, something doesn't seem to be working properly.
Please try refreshing the page. If that doesn't work, please contact support so we can address the problem.
Not just a matter of size: a hospital-level risk factor analysis of MRSA bacteraemia in Scotland
BMC Infectious Diseasesvolume 16, Article number: 222 (2016)
Worldwide, there is a wealth of literature examining patient-level risk factors for methicillin-resistant Staphylococcus aureus (MRSA) bacteraemia. At the hospital-level it is generally accepted that MRSA bacteraemia is more common in larger hospitals. In Scotland, size does not fully explain all the observed variation among hospitals. The aim of this study was to identify risk factors for the presence and rate of MRSA bacteraemia cases in Scottish mainland hospitals. Specific hypotheses regarding hospital size, type and connectivity were examined.
Data from 198 mainland Scottish hospitals (defined as having at least one inpatient per year) were analysed for financial year 2007-08 using logistic regression (Model 1: presence/absence of MRSA bacteraemia) and Poisson regression (Model 2: rate of MRSA bacteraemia). The significance of risk factors representing various measures of hospital size, type and connectivity were investigated.
In Scotland, size was not the only significant risk factor identified for the presence and rate of MRSA bacteraemia. The probability of a hospital having at least one case of MRSA bacteraemia increased with hospital size only if the hospital exceeded a certain level of connectivity. Higher levels of MRSA bacteraemia were associated with the large, highly connected teaching hospitals with high ratios of patients to domestic staff.
A hospital’s level of connectedness within a network may be a better measure of a hospital’s risk of MRSA bacteraemia than size. This result could be used to identify high risk hospitals which would benefit from intensified infection control measures.
Healthcare-associated infections (HAIs) have been an unwanted aspect of healthcare systems for several decades. Staphylococcus aureus is a common pathogen associated with HAIs which results in a spectrum of clinical conditions ranging from mild and often self-limiting skin and soft tissue infections to more serious illnesses including bacteraemia, pneumonia and toxic shock syndrome. Despite representing only a small proportion of all symptomatic S. aureus infections, S. aureus bacteraemias (SAB) have a disproportionate impact on the burden of disease owing to their high associated mortality rates from life-threatening complications such as infective endocarditis and metastatic infections . Therefore, SAB is considered a manifestation of serious disease and of high clinical importance. Furthermore, when methicillin-resistant S.aureus (MRSA) is the causative agent of these bacteraemias, the risk of treatment failure and mortality is greater, with increased costs of managing these patients [2–4].
Implementing effective infection control measures to reduce the number of cases and the nosocomial spread of the bacteria has been a priority in recent years, with measures including enhanced surveillance and infection control policies, as well as government targets. In order to design successful and targeted intervention programs, risk factor analyses are often carried out to identify behaviours and characteristics associated with transmission, morbidity and mortality. Much of the literature on risk factors for MRSA bacteraemia infections has been performed at the patient-level, identifying characteristics such as previous hospitalisation, recent antimicrobial therapy, indwelling vascular devices, prior invasive or surgical procedures, often within risk groups including the elderly and those with latent health conditions [5–7]. Patient-level characteristics, however, do not fully explain the observed heterogeneity in the numbers of MRSA bacteraemia cases among hospitals. It has been suggested that hospital size is associated with the hospital-level incidences [8, 9], however, recent research suggests that size may not be not the only factor [8–13]. Hospital connectedness (or connectivity) through shared transfer patients, was found to influence MRSA bacteraemia incidences [8, 10, 12, 13]. Hospital type and hence the procedures that would be carried out in such facilities was also previously found to be an important hospital-level risk-factor [8, 9].
The aim of this study was to determine whether hospital connectivity and hospital type significantly influence hospital-level MRSA bacteraemia morbidity rates in mainland Scottish hospitals once hospital size has been accounted for. Our approach was twofold; firstly to identify risk factors for the presence of MRSA bacteraemia in all Scottish mainland hospitals and then for General hospitals only, to identify risk factors for high rates of MRSA bacteraemia. Controlling the spread of MRSA bacteraemia depends on knowledge of where the bulk of disease burden lies and understanding why infections occur where they do.
Initially all 264 facilities listed in the financial year 2007-08 Information Services Division (ISD) cost book reports (R020 and R020LS, http://www.isdscotlandarchive.scot.nhs.uk/isd/6058.html) were considered. For the purposes of analysis, a hospital was defined as a secondary healthcare facility with at least one inpatient per year and with at least one hospital specialty (previously known as departments). Hospitals from island Health Boards (NHS Shetland, NHS Orkney and NHS Western Isles) were excluded since hospitals there were considered atypical and not comparable to other ‘mainland’ hospitals included in the study as connectivity is difficult to assess. This resulted in 198 National Health Service (NHS) hospitals being included in the study for the financial year 2007-08 (Fig. 1) and 66 facilities being excluded.
(1) Case data. Methicillin-resistant S. aureus (MRSA) bacteraemia case data by hospital for the financial year 2007-08 (6 April 2007 – 5 April 2008) were received from Health Protection Scotland (HPS). This financial year was chosen because at the time it was the only year for which data were available for all risk factors (hospital size, type and connectivity). (2) Risk factors. Data of potential risk factors were obtained from Information Services Division (ISD) for each hospital for the financial year 2007-08 (a full list of variables in Additional file 1). As the aim was to test specific hypotheses regarding the effect of hospital size, type and connectivity; the majority of the data that were obtained represent different measures of those characteristics. Although various other measures including average occupancy rate, total number of whole time equivalent (medical and dentistry, nursing and midwifery, domestic, and support services) staff, patient-staff ratios (number of patients to medical and dentistry, nursing and midwifery, domestic, and support services staff), and cleaning by hospital size were considered (Additional file 1). (2a) Hospital size. Measures of size included occupied bed days (OBD), surface area (m2), average staffed beds, total inpatients discharged, total staff members and others (Additional file 1). (2b) Hospital type. ISD group Scottish into the following categories: Category A (General (mainly acute) including Teaching (A1), large General (A2), General (A3), and Sick children’s (A4)); Category B (Long stay), Category C (Mental); Category D (psychiatry of learning difficulties); Category E (Maternity); and Category J (Community). Information on the presence and absence of 46 specialties was also acquired. For a complete list of these specialties see Additional file 2. (2c) Connectivity. To quantify movements of patients between hospitals, patient admission data were obtained from ISD. The patient admission data covered all admissions to healthcare facilities in Scotland for the calendar year 2007 (1 January – 31 December 2007). From this dataset the movements of all patients between hospitals were extracted either as direct transfers, i.e. from one hospital directly to another hospital, or as indirect transfers, i.e. when a patient was discharged from one hospital into the community and subsequently (within the period covered by the data) admitted to another hospital. From these data, a movement matrix was derived for all connected hospitals in Scotland  and then various summary measures of hospital connectivity were generated (Table 1).
Descriptive analyses were carried out to describe hospital type and the number of cases per hospital. Several variables were log10(x + 1) transformed to correct for non-normal distribution of the data while others were categorised. A full list of all 39 explanatory variables is shown in Additional file 1. The hospital connectivity variables were examined as both continuous and categorical variables. Cut-off levels of connectivity variables were determined using Receiver Operating Characteristic (ROC) curve analysis , above which hospitals were considered positive for MRSA bacteraemia. The cut-off chosen was one that maximised the sum of the sensitivity and specificity. Confidence intervals (25th and 75th percentiles) for the cut-off values were generated from 10,000 bootstrap simulations. Data regarding hospital specialties (n = 46) were summarised using nonmetric multi-dimensional scaling (NMS) in order to reduce the number of variables. The NMS was performed in PC-ORD version 6.08 (MJM software Design, Gleneden Beach, OR) using the 36 specialties that were sufficiently represented across all hospitals (i.e. present in at least 5 % of the hospitals, therefore 10 specialties were excluded from this analysis). Multi-response Permutation Procedures (MRPP) analysis  was used a posteriori to test the hypothesis of no difference between hospital status (presence/absence of MRSA bacteraemia) and hospital type (Category A, B, C, D, E, J).
For both the single variable and multivariable analyses, two models were considered. Model 1 was a logistic regression model with binomial distribution fitted to presence-absence data for all hospitals (n = 198) to identify risk factors for having cases versus not. Model 2 was a generalised linear model with a Poisson distribution fitted to count data, offset by the logarithm of OBD, to identify risk factors associated with higher rates of MRSA bacteraemia in General (Category A) hospitals only (n = 38). A multiplicative over dispersion parameter was added to address over dispersion in the Poisson model.
A hypothesis driven approach to model selection was performed. Before undertaking the multivariable analysis, groups of variables representing hospital size, connectivity and hospital type were examined individually in a single variable analysis. As size was viewed a priori as important, a size-related variable was entered into the model first. With size in the model, variables associated with connectedness and hospital type were tested, including interactions. Likelihood-ratio tests and the Akaike Information Criterion (AIC) were used to select the most parsimonious model. Statistical significance was set at p < 0.05. Unless otherwise stated all analyses were carried out using Proc Glimmix in SAS version 9.3.1 (SAS Institute Inc., Cary, NC).
Figure 1a-c and Additional file 3 shows the 198 hospitals with respect to the number of MRSA bacteraemia cases (Fig. 1a), hospital size (Fig. 1b) and connectivity (Indegree) (Fig. 1c). Hospitals are colour coded by “type” as defined by the NMS analysis. The majority of hospitals (n = 152, 77 %) reported no MRSA bacteraemia cases (Fig. 1a). The remaining 46 hospitals reported 662 MRSA bacteraemia cases in total (range 0-72). From these 662 cases, 47 cases died from all-cause mortality within 30 day, giving a 30-day case-fatality ratio of 7.1 %.
Single variable analysis of risk factors
Figure 1b shows the distribution of MRSA bacteraemia cases with respect to hospital size. Hospitals were ranked by size, as measured by occupied bed days. Hospital size was a statistically significant predictor of the presence of cases (p < 0.001) and number of cases in General hospitals (p < 0.001).
Table 1 shows the best cut-off levels for the connectivity measures identified using ROC curve analysis. For all cut-offs found, hospitals above the threshold were at significantly increased risk for the presence of MRSA bacteraemia (Table 1).
Figure 2 shows the results of NMS of the specialties for the 198 hospitals included in this study. Two significant axes explained a total of 83.9 % of the variation in the data. NMS axis 1 which explains 78.7 % of the variation, significantly separated hospitals with and without MRSA (MRPP, p < 0.001). Progression along NMS Axis 1 was associated with an increasing number of MRSA bacteraemia cases, increasing hospital size (OBD) and an increasing number of specialties. Three significantly (MRPP, p < 0.001) different hospital groups were identified, distributed along NMS Axis 1 (Fig. 2). The “Other” group corresponds to the ISD hospital designated categories B, C, D, E, J. Hospitals designated “Other” were at lower risk than General Hospitals with no Teaching (Categories A2, A3, A4) which in turn had a lower risk than Teaching hospitals (Category A1). In general, Teaching hospitals were the largest, most connected and had the highest number of MRSA cases (Fig. 1a-c).
Multivariable analysis of risk factors
Model 1: Risk factors for presence of MRSA bacteraemia cases (n = 198 hospitals). Results showed (Fig. 3 ) that the probability of a hospital having at least one MRSA bacteraemia case increased with hospital size (log10OBD), particularly for hospitals that exceeded the threshold for the connectivity variable Indegree (Log10OBD*Indegree threshold p = 0.05; odds ratio (95%CI): 136 (7- > 999)).
Model 2: Risk factors for higher rates of MRSA bacteraemia (n = 38 hospitals). After adjusting for size, results showed (Table 2) that the Indegree (number of hospitals that transferred patients to a given hospital), being a teaching hospital and the ratio of patients-to-domestic staff were all significant predictors of the level of MRSA bacteraemia in General hospitals. For each hospital sending at least one patient (Indegree), there is a 3 % increase in the expected number of MRSA cases per OBD. The expected number of MRSA bacteraemia cases per OBD is 64 % higher in teaching than non-teaching hospitals. Higher ratios of patients to domestic staff is associated with a 28 % increase in the expected number of MRSA cases per OBD for a 10 fold increase in the ratio, corresponding to a 2.8 % increase for a 26 % increase in the ratio (Table 2).
In this study, data on MRSA bacteraemia were obtained for all Scottish mainland hospitals, defined as a secondary healthcare facility with at least one inpatient per year and with at least one hospital specialty. Only 23 % of the hospitals examined in this study reported at least one case of MRSA bacteraemia. This proportion was small however; the hospital definition included long stay, mental and community hospitals which were typically small facilities, providing very few specialties. In general, hospitals that reported at least one case of MRSA bacteraemia over the study period were large hospitals. However, even large hospitals rarely had large numbers of reported MRSA bacteraemia cases unless they were also well connected to the wider hospital network. For all connectivity measures generated in this study, there seemed to be a cut-off below which a hospital had no or lower MRSA bacteraemia cases than expected from a model that included size alone. Although further testing would be required, such a measure is highly useful in identifying hospitals at high risk for MRSA bacteraemia (and potentially other HAIs) and could direct where to concentrate surveillance efforts at a national level.
Hospital connectivity was also a significant risk factor in the level of MRSA bacteraemia observed in General hospitals thus providing further evidence that size alone cannot fully predict a hospital’s MRSA bacteraemia status. Connectivity has also been found to be an important risk factor for the spread of nosocomial or antibiotic-resistant pathogens in other studies [8, 10, 11, 13] and described in a recent review . Well-connected hospitals received volumes of transfer patients from several other facilities. Previous studies have shown that transfer patients have a comparatively higher risk of nosocomial infection than non-transferred patients owing to different demographics, health status at time of transfer (often acutely unwell and critical) and other patient-level risk characteristics .
When examining hospital type, Teaching hospitals were highlighted as having a large number of MRSA bacteraemia cases. Teaching hospitals are large, well connected and have a high number of specialties. High patient-to-staff ratios were also identified as a risk factor for higher levels of MRSA bacteraemia. Patient-to-staff ratio of domestic workers was used as a proxy measure for the ratio of patients to cleaning staff. Hospital cleaning remains a vital strategy for effective infection control which was highlighted by this small but significant risk from our model. This result, however, supports the previously reported effect of heightened hygiene awareness and enhanced cleaning on reducing the number of new infections and microbial contamination of the hospital [13, 17–19].
This study has some important limitations that must be considered. Our definition of a hospital was intended to capture facilities with overnight patients receiving medical care. Nursing homes were excluded since not all patients would be receiving medical care, the demographics and risk factors of patients would likely differ from true hospitals, and nursing homes operate in a different way to hospitals. Nursing homes are thought to harbour a substantial number of S. aureus colonisations and infections are therefore likely to be an important reservoir for hospital infections [20–23]. In addition, this study only used a single year of data. However, almost 90 % of the hospitals had the same status (either MRSA bacteraemia positive or negative) over the 5 years (financial years 2006-07 to 2010-11) for which data were available. There were three hospitals that were only positive for the year of this study. All models were run with and without these hospitals with no change in the results. The same applies to the hospital network; there were no changes in the network over the studied period. Although there may be day-to-day variation in connections between hospitals, the network is stable when measured over longer periods (van Bunnik, personal communication). Lastly, since we did not have information on individual level risk factors we could not take these into account in the models presented in this paper. However, future models should ideally take into account both individual level effects and hospital level effects to identify a complete set of risk factors for the presence of MRSA bacteraemia cases in Scottish hospitals.
This is the first Scotland-wide study utilising hospital-level characteristics to examine differences in MRSA bacteraemia morbidity among hospitals. This study showed that 23 % of Scottish mainland hospitals reported cases of MRSA bacteraemia. Effort to control MRSA bacteraemia should therefore be focused on large, well-connected hospitals. Further research should be conducted with respect to reducing the ratio of patients to domestic staff as a means of reducing the levels of MRSA bacteraemia infections, especially in teaching hospitals where the levels are very high.
Health Protection Scotland (HPS) and Scottish Government prioritise the implementation of effective infection control measures to reduce the spread of pathogens in nosocomial settings. Mandatory enhanced surveillance for all S. aureus was introduced in Scotland in 2014. Programmes that prevent MRSA bacteraemia transmission will likely be relevant to other epidemiologically important healthcare pathogens that spread by patient-to-patient transmission. Here, we presented a study that adds to existing literature since our multivariate models could readily be applied as public health tools that estimate hospital-level risk of MRSA bacteraemia and that identify high risk hospitals in Scotland which would benefit from intensified infection control measures. With further investigation, the application could potentially be extended to other HAIs and other countries.
Ethics approval and consent to participate
Approval for this project and release of these data from HPS was received from the Caldicott Guardian for NHS Greater Glasgow and Clyde (PAC 54/11).
Availability of data and materials
The data within this manuscript is confidential and cannot be shared. Data of interest has been provided in aggregated form in Additional file 3.
Staphylococcus aureus bacteraemia
methicillin-resistant Staphylococcus aureus
National Health Service
Information Services Division
Health Protection Scotland
nonmetric multidimensional scaling
receiver operating characteristic
Akaike Information Criterion
occupied bed days
Johnson AP. Methicillin-resistant Staphylococcus aureus: the European landscape. J Antimicrob Chemother. 2011;66 suppl 4:iv43–8.
Hunt AC, Edwards B, Girvan EK, Cosgrove B, Edwards GFS, Gould IM. Methicillin-resistant Staphylococcus aureus in Northeastern Scotland in 2003 to 2007: evolving strain distribution and resistance patterns. J Clin Microbiol. 2011;49(5):1975–8.
Reilly JS, Stewart S, Christie P, Allardice G, Smith A, Masterton R, Gould IM, Williams C. Universal screening for meticillin-resistant Staphylococcus aureus: interim results from the NHS Scotland pathfinder project. J Hosp Infect. 2010;74(1):35–41.
Wyllie DH, Peto TEA, Crook D. MRSA bacteraemia in patients on arrival in hospital: a cohort study in Oxfordshire 1997-2003. BMJ. 2005;331(7523):992.
Elston JWT, Barlow GD. Community-associated MRSA in the United Kingdom. J Infect. 2009;59(3):149–55.
Kock R, Becker K, Cookson B, van Gemert-Pijnen JE, Harbarth S, Kluytmans J, et al. Methicillin-resistant Staphylococcus aureus (MRSA): burden of disease and control challenges in Europe. Euro Surveill. 2010;15(41):19688.
Cookson B, Bonten MJM, MacKenzie FM, Skov RL, Verbrugh HA, Tacconelli E. Meticillin-resistant Staphylococcus aureus (MRSA): screening and decolonisation. Int J Antimicrob AG. 2010;37(3):195–201.
Donker T, Wallinga J, Grundmann H. Patient referral patterns and the spread of hospital-acquired infections through National Health Care Networks. PLoS Comput Biol. 2010;6(3):e1000715.
Wakefield DS, Pfaller M, Massanari RM, Hammons GT. Variation in methicillin-resistant Staphylococcus aureus occurrence by geographic location and hospital characteristics. Infect Control. 1987;8(4):151–7.
Donker T, Wallinga J, Slack R, Grundmann H. Hospital networks and the dispersal of hospital-acquired pathogens by patient transfer. PLoS ONE. 2012;7(4):25.
Ciccolini M, Donker T, Grundmann H, Bonten MJM, Woolhouse MEJ. Efficient surveillance for healthcare-associated infections spreading between hospitals. Proc Nat Acad Sci USA. 2014;111(6):2271–6.
Eveillard M, Quenon JL, Rufat PMD, Mangeol APD, Fauvelle FPD. Association between hospital-acquired infections and patients’ transfers. Inf Cont Hosp EP. 2001;22(11):693–6.
Donker T, Wallinga J, Grundmann H. Dispersal of antibiotic-resistant high-risk clones by hospital networks: changing the patient direction can make all the difference. J Hosp Infect. 2014;86(1):34–41.
Sing T, Sander O, Beerenwinkel N, Lengauer T. ROCR: visualizing classifier performance in R. Bioinformatics. 2005;21:3940–1.
McCune B, Grace JB. Analysis of Ecological Communities. Gleneden Beach, Oregon: MjM Software Design; 2002.
Ciccolini M, Donker T, Köck R, Mielke M, Hendrix R, Jurke A, Rahamat-Langendoen J, Becker K, Niesters HGM, Grundmann H, Friedrich AW. Infection prevention in a connected world: The case for a regional approach. Int J Med Microbiol. 2013;303:380–7.
Dancer SJ, White L, Robertson C. Monitoring environmental cleanliness on two surgical wards. Int J Environ Heal R. 2008;18(5):357–64.
Dancer SJ, White LF, Lamb J, Girvan EK, Robertson C. Measuring the effect of enhanced cleaning in a UK hospital: a prospective cross-over study. BMC Med. 2009;7(28):1741–7015.
Petti S, Polimeni A, Dancer SJ. Effect of disposable barriers, disinfection, and cleaning on controlling methicillin-resistant Staphylococcus aureus environmental contamination. Am J Infect Control. 2013;41(9):836–40.
Hudson LO, Reynolds C, Spratt BG, Enright MC, Quan V, Kim D, et al. Diversity of methicillin-resistant Staphylococcus aureus strains isolated from residents of 26 nursing homes in Orange County, California. J Clin Microbiol. 2013;51(11):3788–95.
McKinnell JA, Miller LG, Eells SJ, Cui E, Huang SS. A systematic literature review and meta-analysis of factors associated with methicillin-resistant Staphylococcus aureus colonization at time of hospital or intensive care unit admission. Infect Control Hosp EP. 2013;34(10):1077–86.
Jans B, Schoevaerdts D, Huang TD, Berhin C, Latour K, Bogaerts P et al. Epidemiology of multidrug-resistant microorganisms among nursing home residents in Belgium. PLoS ONE. 2013; 8(5):1–8.
Lee BY, Bartsch SM, Wong KF, Singh A, Avery TR, Kim DS, et al. The importance of nursing homes in the spread of methicillin-resistant Staphylococcus aureus (MRSA) among hospitals. Med Care. 2013;51(3):205–15.
Wassernam S, Faust K. Social Network Analysis: Methods and Applications. Cambridge: Cambridge University Press; 1994.
Zou G. A modified poisson regression approach to prospective studies with binary data. Am J Epidemiol. 2004;159:702–6.
We acknowledge staff at Health Protection Scotland (HPS) and Information Services Division (ISD Scotland) (both part of Public Health and Intelligence (PHI), under NHS National Services Scotland), and the Scottish MRSA reference laboratory (SMRSARL) for providing data and background information that facilitated this study.
This work was supported by the Scottish Infection Research Network (SIRN), Scottish University Life Sciences Alliance (SULSA), MRC, Zoetis and European Commission’s Seventh Framework the European Union Seventh Framework Programme (FP7-HEALTH-2011-single-stage) “Evolution and Transfer of Antibiotic Resistance” (EvoTAR). T.P. is supported by the Centre of Expertise on Animal Disease Outbreaks (EPIC) funded by the Scottish Government (Programmes of Research 2011-2016 no. UEH/851/11). J.R.F. is supported by institute strategic grant funding from the Biotechnology and Biological Sciences Research Council. M.C.T. was supported by a Wellcome Trust strategic award to the Centre for Immunity, Infection and Evolution (WT095831).
The authors declare that they have no competing interests.
CG, MECT, JRF and MEJW conceived the study. CG acquired and analysed the data and wrote the first draft of the manuscript. BADvB generated the connectivity measures used in this study. OB, JW, LI and CR provided the case data for hospitals in Scotland. CR and TP provided additional statistical support. MECT prepared the final draft for submission. All authors read and approved the final manuscript.
List of all variables considered in the risk factor analysis (Model 1: Logistic regression + Model 2: Poisson regression. (DOCX 17 kb)
List of all 46 hospital specialities that were obtained from Information Services Division (ISD). Results of univariate screening using Chi-square or Fisher’s exact where appropriate and Odds Ratios (for specialities that with statistically significant chisquare tests and no cells where n = 0). Specialities in bold print (n = 36) were used in the Nonmetric Multidimensional Scaling (NMS) carried out (Fig. 3). (DOCX 18 kb)
Summary table of hospital characteristics based on one measure of hospital size, Occupied Bed Days (OBD). OBD is divided into 4 categories based on 1st quartile (5550), median (12168) and 3rd quartile (34135). Unless specified the values represent the range within a given OBD category (min – max). (DOCX 14 kb)