Assessing the role of undetected colonization and isolation precautions in reducing Methicillin-Resistant Staphylococcus aureustransmission in intensive care units
© Kypraios et al; licensee BioMed Central Ltd. 2010
Received: 19 May 2009
Accepted: 16 February 2010
Published: 16 February 2010
Screening and isolation are central components of hospital methicillin-resistant Staphylococcus aureus (MRSA) control policies. Their prevention of patient-to-patient spread depends on minimizing undetected and unisolated MRSA-positive patient days. Estimating these MRSA-positive patient days and the reduction in transmission due to isolation presents a major methodological challenge, but is essential for assessing both the value of existing control policies and the potential benefit of new rapid MRSA detection technologies. Recent methodological developments have made it possible to estimate these quantities using routine surveillance data.
Colonization data from admission and weekly nares cultures were collected from eight single-bed adult intensive care units (ICUs) over 17 months. Detected MRSA-positive patients were isolated using single rooms and barrier precautions. Data were analyzed using stochastic transmission models and model fitting was performed within a Bayesian framework using a Markov chain Monte Carlo algorithm, imputing unobserved MRSA carriage events.
Models estimated the mean percent of colonized-patient-days attributed to undetected carriers as 14.1% (95% CI (11.7, 16.5)) averaged across ICUs. The percent of colonized-patient-days attributed to patients awaiting results averaged 7.8% (6.2, 9.2). Overall, the ratio of estimated transmission rates from unisolated MRSA-positive patients and those under barrier precautions was 1.34 (0.45, 3.97), but varied widely across ICUs.
Screening consistently detected >80% of colonized-patient-days. Estimates of the effectiveness of barrier precautions showed considerable uncertainty, but in all units except burns/general surgery and one cardiac surgery ICU, the best estimates were consistent with reductions in transmission associated with barrier precautions.
Patient isolation measures (such as barrier precautions or physical isolation of patients in single rooms or cohorts) and screening to detect MRSA colonization are core components of almost all policies for preventing the nosocomial transmission of methicillin-resistant Staphylococcus aureus (MRSA). However, there are few well-designed intervention studies considering screening and isolation precautions in the absence of other containment measures, and little is known about the extent to which they reduce transmission. Obtaining estimates of any such reduction is important for both quantifying the value of existing practices and the potential benefit of new rapid screening technologies.
Screening and isolation measures have two complementary aims: to minimize the number of days MRSA-positive patients are left unisolated; and to minimize MRSA transmission from isolated patients. In this paper we consider the extent to which we can quantify the success of control programs in achieving both aims using routine surveillance data. This presents a number of challenges. First, the number of unisolated MRSA patient days usually cannot be directly observed using surveillance data, because in most settings those who are known to be MRSA-positive are immediately isolated. Moreover, screening cultures are typically taken no more than weekly. This, combined with less than perfect swab sensitivity, makes it impossible to directly observe the total number of colonized patient days, the number of transmission events and the precise times at which these events occur. Furthermore, estimating the effectiveness of isolation measures in reducing transmission requires quantifying the relative transmission rates of both known and isolated and unobserved and unisolated MRSA positive patients. This cannot be done without estimates of the number of days MRSA-positive patients are left unisolated. Together, these obstacles prevent direct estimates of the efficacy of patient isolation in reducing transmission using routine surveillance data.
Here we show that data augmentation methods allow these problems to be overcome by treating the unobserved times at which patients acquire MRSA as unknown parameters to be estimated. By analysing such surveillance data using stochastic transmission models and making use of data augmentation techniques it is possible to simultaneously estimate the number of undetected MRSA positive patient days and the efficacy of isolation measures in reducing transmission using routinely collected MRSA surveillance data.
We collected detailed census and microbiologic data from eight adult 10-bed intensive care units (ICUs) in a tertiary academic medical center in Boston, Massachusetts where routine admission and weekly bilateral nares screening for MRSA was occurring with high compliance (90%). Types of ICUs included medical, cardiac, general surgery (2), neurosurgery, thoracic surgery, and cardiac surgery (2). One of the two general surgery wards also received burn and trauma patients. Dates of MRSA-positive clinical cultures (all sources), as well as positive and negative MRSA nares screening cultures were collected during a 17-month period from September 2003 through January 2005. Negative clinical cultures were not assessed. All specimens were processed by routine bacterial culturing techniques. This study was approved by the Institutional Review Board at Brigham and Women's Hospital with a waiver of informed consent.
Newly-identified and previously known MRSA-positive patients were placed under contact precautions consisting of gown and glove use as well as use of single rooms (all ICU rooms were single occupant). Once identified as MRSA-positive, contact precautions were applied on that admission and all subsequent admissions. No nurse cohorting was utilized. Dates of each ICU admission and discharge were obtained, along with the date on which contact precautions were initially applied (for MRSA or other highly antibiotic resistant pathogens). The first institutional date of a MRSA-positive culture was also recorded, even if it preceded the study period.
Our baseline analysis used a previously described dynamic stochastic single-ward transmission model to analyse the data . At any point in time, each patient is assumed by this model to be in one of two states, colonized (defined here as the presence of MRSA at any body site, regardless of symptoms) or uncolonized. Each patient entering the ward has a probability ϕ of being already colonized, unless they are already known to be so (i.e. having already had a previous positive test).
During a stay on the ward, an initially uncolonized patient is assumed to have a risk of becoming colonized that can be described by a Poisson process: this means that the probability of a patient becoming colonized between time t and time t+τ can be approximated by λ (t) τ (where τ is a small interval of time, the approximation becoming exact as τ approaches zero). Consequently, a patient's chances of becoming colonized increase with length of stay on the ward. The rate λ (t) can vary through time; since we are interested in comparing different potential sources of colonization we assume that λ (t) = β 0 + C(t)β 1 + I(t)β 2 . Here, β0 is the rate of background transmission, β1 is the rate of transmission due to colonized but nonisolated patients, β2 is the rate of transmission due to colonized and isolated patients, and C(t) and I(t) are the numbers of nonisolated and isolated colonized patients at time t, respectively. We thus assume that the "colonization pressure" on an uncolonized individual increases linearly with the number of colonized patients who are isolated and nonisolated. Although it is common practice to use a scaling such as β1/N instead of β1, with N the typical number of patients on the ward, this did not seem necessary since there was very little variation in ward occupancy both between wards and between different times during the study period. The inclusion of the background rate β0 models the assumption of a constant risk of colonization irrespective of the presence of colonized patients. Such background transmission could, for example, result from environmental contamination or contact with staff carriers of MRSA [2–4]. Note that the model takes no explicit account of HCW compliance with barrier precautions.
Once a patient is colonized, we assumed that he or she remains so for a three month period. Thus, a patient colonized on discharge is assumed to still be colonized if readmitted within three months. This assumption was not critical: the results presented below were virtually unchanged when different periods of colonization were assumed (e.g. 6 months, 9 months). We assumed that the nares culture had a specificity of 100% and a sensitivity of p × 100%, where p (the sensitivity of a nasal swab for detecting MRSA carriage at any site) was estimated from the data.
The above model has five main parameters, namely ϕ, p, β0,β1, and β2. We wish to estimate these from the data, which consist of admission times, discharge times, and the times and outcomes of tests. Since our approach involves estimating unseen colonization times, we can also estimate quantities such as the duration of carriage prior to a positive test result. Assessment of the effectiveness of isolation is performed by comparing estimates of β1 and β2: evidence that β1 > β2would support the hypothesis that isolation is effective in reducing transmission. Specifically, we estimate P(β1 > β2) and the ratio β1/β2 . Estimation was carried out within a Bayesian statistical framework using Markov chain Monte Carlo (MCMC) methods .
The five model parameters were assigned uninformative and independent prior distributions: Uniform(0,1) distributions for ϕ and p; and exponential distributions with rate 10-6 for each of the β parameters. Since the likelihood of the observed data (dates and outcomes) given the model parameters is computationally intractable, we used a data augmentation method in which the unobserved colonization times were included as additional parameters. This yields an augmented posterior density that is known up to proportionality, which we explored using an MCMC algorithm. Within the algorithm, infection rate parameters were updated using Metropolis-Hastings steps, while both ϕ and p were updated using Gibbs steps .
Our methods enable us to estimate the percent of patient days that are attributed to patients who are colonized but not yet detected (p hidden ), and the percent of patient days attributed to patients who are colonized and have been tested, but who are awaiting test results (p wait ). Note that p hidden > p wait since all "waiting" patients are also "hidden".
Some of the parameter estimates were combined using a random effects model with inverse variance weights [, Sec. 5.2] to derive pooled estimates and corresponding standard errors. These computations were carried out using the rmeta package for R http://www.r-project.org. All other analysis, such as the implementation of the MCMC algorithms, was performed using programs we wrote in C.
We refer to the model described above as the full model. To assess the impact of our assumptions, we repeated the analyses using two alternative models. In the first (the no-background model), β0 is assumed to be zero with high probability (specifically, β0 was assigned an exponential prior distribution with mean 10-4 so P(β0 < 0.001) > 0.99). This corresponds to the assumption that virtually all of the MRSA acquisitions result from patient-to-patient transmission (much of which is presumably mediated by transiently colonized healthcare workers, though this is not explicitly modeled). Background transmission unrelated to colonization pressure, e.g. environmental contamination, is assumed to be of negligible importance. In the second (the non-linear model), the assumption of linearly increasing colonization pressure is relaxed, and it is instead assumed that colonization pressure due to both isolated and nonisolated individuals remains constant provided at least one colonized individual is present. Although biologically less plausible than the full model, this model represents an extreme case and by fitting it we obtain insights into the impact of the assumptions underlying the full model. For both the non-linear and no-background models, assessment of the effect of isolation is measured via comparison of the two rates β1 and β2 as described above.
We consider two definitions of swab sensitivity: swab sensitivity in detecting colonization of the nares, and swab sensitivity for detecting colonization at any body site. The former is calculated from a subset of the data (without making use of the model) by comparing subsequent swabs in patients who have a first positive swab. Precisely, the sensitivity is estimated by the ratio TP/(TP+FN) using serial nares cultures alone, where TP and FN denote total numbers of true positive and false negative tests on a ward excluding the first positive swab from each patient. Since patients, once colonized, are assumed to remain so throughout their ICU stay, TP is simply the number of positive swabs excluding the first for each patient episode, and FN the number of negative swabs that follow an earlier positive for the same patient. The second reported sensitivity is the parameter p estimated from the model making use of cultures from all sites (this is necessarily lower than the sensitivity for detecting nares carriage).
Summary statistics for the study data
Number of patients
Length of stay Mean(SD)
Percent in contact precautionb
Number of swab tests per person
Number of MRSA+ swab tests per person
Number of swab tests taken after first positivec
Observed prevalence and incidence for the study data
Monthly MRSA incidence density (new cases per 1,000 patient days)
The impact of undetected colonized patients on transmission
Assessment of undetected colonization
Is isolation effective? Comparison of β1 and β2
Assessment of isolation efficacy
P(β1 > β2)
P(β1 > β2)
P(β1 > β2)
Pooled estimate of the effectiveness of isolation
Two measures of nares test sensitivity based upon serial nares cultures
Sensitivity of Detecting Nares Colonization (%)a
Sensitivity of detecting colonization at any body site (%): Median (95% CI)b
64 (56, 71)
56 (47, 65)
61 (54, 69)
50 (42, 57)
62 (54, 68)
68 (60, 76)
55 (45, 64)
54 (44, 65)
59 (55, 63)
Percentage of new admissions already colonized
Each patient admitted to a ward during the study period may be known to be colonized due to a previous positive test. If not, colonization status is unknown, since the patient may be either new to the study, or else previously admitted but with no positive test result. The parameter ϕ in our model represents the proportion of such patients who are colonized on admission. Note that ϕ should not be interpreted as a measure of admission prevalence, since it excludes known positives.
Percentage of patients who are colonized on admission to a ward
Percent having a previous positive test a
Estimated percent of MRSA carriers among admissions with unknown statusb(ϕ × 100%): Median (95% CI)
9.4 (7.3, 12)
9.6 (7.5, 12.2)
14.1 (10.3, 18.6)
14.4 (10.7, 19)
10.4 (7.8, 13.4)
10.7 (8.1, 13.9)
8.3 (5.4, 12.1)
9.0 (6, 13)
8.7 (5.9, 12.2)
9.2 (6.4, 12.7)
5.7 (3.8, 8.2)
6.1 (4.2, 8.6)
3.4 (2.1, 5.3)
3.8 (2.4, 5.8)
4.6 (2.9, 6.9)
5.9 (4.0, 8.4)
We have shown that a data-augmentation approach allows us to estimate the extent of undetected MRSA carriage and the effectiveness of barrier precautions (gown and gloves) in preventing transmission using only routinely collected high-compliance admission and weekly MRSA screening cultures. Such estimates would have been difficult or impossible with more conventional approaches and would have required assuming rather than estimating the unobserved carriage data.
Our results have shown that admission and weekly nares screening with conventional culture methods are likely to have detected at least 80% of MRSA colonized patient days in the ICUs studied, despite 48-hour culture turn-around times, and sensitivities to detect MRSA carriage of only 85% for detecting nares carriage and 59% for carriage at any body site. The explanation for this high detection of MRSA patient days despite a low sensitivity for non-nares sites is that clinical cultures of other sites are routinely being performed for medical reasons. Thus, nares screening in conjunction with usual routine culturing performed for medical reasons captures a large majority of positive MRSA carriers. While generalizability to non-ICU wards remains unclear, this level of detection in ICUs is reassuring since risk of MRSA infection is highest in these settings.
We estimate that approximately 10% of potential contact precaution days are missed due to delays in obtaining test results. If contact precautions cause a 28.1% reduction in the transmission rate (i.e. 1 - (β2/β1), using the pooled estimate from all the wards), then an additional 5% reduction in total transmission rates would be gained with instantaneous swab results, suggesting a limit to the potential benefit in using rapid screening technologies. However, the potential benefit of rapid test results would increase in line with the effectiveness of isolation measures.
While our best estimate is that gown and glove barrier precautions produce a 28% reduction in MRSA transmission, we recognize that these estimates have considerable uncertainty despite including over 11 ward-years of data and there is a substantial chance that the effect is both much lower and much larger. Thus we estimate that there is 30% chance that the reduction in transmission is above 44% and a 30% chance that isolation actually increases transmission.
This uncertainty may indicate that many more years of data are required, and/or may be explained by the fact that data from different types of ICUs were pooled. However the Woolf's test for heterogeneity did not reach significance at a 5% level (p = 0.982), indicating that these differences are consistent with chance. Moreover since comparisons between ward-types were not planned they may be of value for generating, but not testing hypotheses.
Recognizing these caveats, it is plausible that barrier precautions have differential effectiveness in different types of ICUs. The estimates found in this study suggest that medical ICUs may derive the most benefit from barrier precautions. This is important since the MRSA prevalence in medical ICUs is often higher than in other ICU types. Interestingly, our results suggest that no benefit was derived in general surgery ICUs serving burn and trauma patients. It is possible that barrier precautions may be less effective when the burden of MRSA is profuse, such as with surgical wounds or burns. Such an interpretation is supported by the fact that the rate of background transmission (β0) is substantially higher on these units. This is also supported by other work in these units suggesting that contamination with MRSA is higher in surgical versus medical ICUs . Another possible explanation for the lack of benefit seen in select wards is low compliance with contact precautions, although this was not assessed. Finally, if barrier precautions provide a false sense of security and lead to poorer hand hygiene, then transmission might effectively worsen with this practice.
It is also possible that the differing effectiveness of barrier precautions among ICUs could be due to unmeasured effects. Specifically, we did not measure or account for variations in host risk factors for acquisition (comorbidities, severity of illness, wounds, devices), nor did we measure the effect of other ongoing infection control interventions such as hand hygiene campaigns. We do note that infection control interventions were not differentially applied to certain ICUs during this time, but levels of compliance are not known. In the absence of reliable data about these factors, none were considered in our analysis.
Our results may underestimate the full effect of barrier precautions since all ICU rooms were single occupancy; thus, we could not assess the impact of this component. In addition, studies have suggested that compliance with barrier precautions is around 70% [9, 10]. Our estimates represent estimates of "effectiveness" (under routine operation) rather than a theoretical efficacy that would apply only under perfect compliance. Thus, while our estimates can be considered conservative, we believe they are the most practically relevant outcome measures. More precise estimates would require more prolonged surveillance data, data with more frequent swabbing intervals, and/or extensive typing information to help identify transmission routes.
The model-based approach that we have adopted for data analysis has several advantages over standard statistical methods, yielding more powerful analyses based on more realistic assumptions. In particular, it accounts for unobserved but important events (e.g. the time at which a patient becomes colonized), and attempts to estimate these events and their uncertainty. It also is based upon biologically meaningful assumptions regarding transmission, addressing both environmental sources and patient-to-patient transfer.
In summary, we have used a stochastic model to estimate the importance of undetected MRSA carriage and the impact of gown and glove precautions on MRSA transmission from 11-ICU years of active MRSA surveillance data. Such isolation precautions are the foundation of infection control guidelines for MRSA control. Nevertheless, there are concerns that isolation techniques may reduce the quality of patient care and incur risks related to inattention [11–13]. Thus, quantifying the beneficial effects and understanding the types of hospital wards in which such measures are likely to be insufficient is a critical part of weighing the need for additional or alternative measures. This work shows that even with conventional MRSA detection technology and weekly screening, a very large proportion of the total patient MRSA reservoir can be detected. It also suggests that barrier precautions may afford substantial benefit in medical and possibly specialty surgical ICUs, but that barriers to its effectiveness may need further study. The large uncertainty in effectiveness estimates, however, illustrate the limitations of routine surveillance data and highlight the need for well-designed prospective intervention studies to evaluate such interventions with greater precision.
In a study of 8 ICUs performing admission and weekly nares surveillance, the average percentage of colonized-patient-days attributed to undetected carriers was 14%, while the percentage of colonized-patient-days attributed to patients awaiting test results was 8%. This suggests that nares surveillance identifies a large majority of carriers and that pcr testing may confer only a small benefit over routine culture. Estimates of the effectiveness of barrier precautions showed an overall benefit of 25%, but this benefit varied widely across different types of ICUs.
TK would like to acknowledge funding from Wellcome Trust (Ref: 076850). Data collection was funded by the CDC Prevention Epicenters Program and the National Institutes of Health (K23AI64161, Huang).
- Forrester ML, Pettitt AN, Gibson GJ: Bayesian inference of hospital-acquired infectious diseases and control measures given imperfect surveillance data. Biostatistics. 2007, 8 (2): 383-401. 10.1093/biostatistics/kxl017.View ArticlePubMed
- Snyder GM: Detection of methicillin-resistant Staphylococcus aureus and vancomycin-resistant enterococci on the gowns and gloves of healthcare workers. Infect Control Hosp Epidemiol. 2008, 29 (7): 583-589. 10.1086/588701.PubMed CentralView ArticlePubMed
- Huang SS, Datta R, Platt R: Risk of acquiring antibiotic-resistant bacteria from prior room occupants. Arch Intern Med. 2006, 166: 1945-1951. 10.1001/archinte.166.18.1945.View ArticlePubMed
- Boyce JM, Havill NL, Otter JA, Adams NM: Widespread environmental contamination associated with patients with diarrhea and methicillin-resistant Staphylococcus aureus colonization of the gastrointestinal tract. Infect Control Hosp Epidemiol. 2007, 28: 1142-1147. 10.1086/520737.View ArticlePubMed
- O'Neill PD: A tutorial introduction to Bayesian inference for stochastic epidemic models using Markov chain Monte Carlo methods. Mathematical Biosciences. 2002, 180: 103-114. 10.1016/S0025-5564(02)00109-8.View ArticlePubMed
- Gilks WR, Richardson S, Spiegelhalter DJ, Eds: Markov chain Monte Carlo in practice. 1996, London: Chapman and Hall
- Sutton AJ, Abrams KR, Jones DR, Sheldon TA, Song F: Methods for meta-analysis in medical research. 2000, Chichester: Wiley
- Goodman ER, Platt R, Bass R, Onderdonk AB, Yokoe DS, Huang SS: Impact of an environmental cleaning intervention on cultures for MRSA and VRE. Infect Control and Hosp Epidemiol. 2008, 29 (7): 593-599. 10.1086/588566.View Article
- Weber DJ, Sickbert-Bennett EE, Brown VM, Brooks RH, Kittrell IP, Featherstone BJ, Adams TL, Rutala WA: Compliance with isolation precautions at a university hospital. Infect Control Hosp Epidemiol. 2007, 28 (3): 358-361. 10.1086/510871.View ArticlePubMed
- Manian FA, Ponzillo JJ: Compliance with routine use of gowns by healthcare workers (HCWs) and non-HCW visitors on entry into the rooms of patients under contact precautions. Infect Control Hosp Epidemiol. 2007, 28 (3): 337-340. 10.1086/510811.View ArticlePubMed
- Kirkland KB, Weinstein JM: Adverse effects of isolation. Lancet. 1999, 354: 1177-1178. 10.1016/S0140-6736(99)04196-3.View ArticlePubMed
- Saint S, Higgins LA, Nallamothu BK, Chenoweth C: Do physicians examine patients in contact isolation less frequently? A brief report. Am J Infect Control. 2003, 31: 354-356. 10.1016/S0196-6553(02)48250-8.View ArticlePubMed
- Catalano G, Houston SH, Catalano MC, Butera AS, Jennings SM, Hakala SM, Burrows SL, Hickey MG, Duss CV, Skelton DN, Laliotis GJ: Anxiety and depression in hospitalized patients in resistant organism isolation. South Med J. 2003, 96: 141-145. 10.1097/01.SMJ.0000050683.36014.2E.View ArticlePubMed
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2334/10/29/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.