Skip to main content

Advertisement

Download PDF

Acquisition and persistence of strain-specific methicillin-resistant Staphylococcus aureus and their determinants in community nursing homes

BMC Infectious Diseases volume 17, Article number: 752 (2017) Cite this article

Abstract

Background

Nursing home residents are frequently colonized with various strains of methicillin-resistant Staphylococcus aureus (MRSA) but the intra-facility dynamics of strain-specific MRSA remains poorly understood. We aimed at identifying and quantifying the associations between acquisition and carriage of MRSA strains and their potential risk factors in community nursing homes using mathematical modeling.

Methods

The data was collected during a longitudinal MRSA surveillance study in six nursing homes in South Central Wisconsin. MRSA cultures were obtained from subjects every 3 months for up to one year. MRSA isolates were subsequently strain-typed by pulsed-field gel electrophoresis (PFGE), and their genetic similarity was established based on the Dice coefficients. Bayesian network analysis, logistic regression and elastic net were used to quantify the associations between acquisition and carriage of MRSA strains discriminated at 80% and 95% strain similarity thresholds and potentially modifiable resident characteristics including previous antibiotic exposure, comorbidity, medical devices, chronic wounds, functional and cognitive status and recent hospitalizations.

Results

Absence of severe cognitive impairment as well as presence of a wound, device and severe comorbidity was associated with elevated probability of USA100 carriage although there was a variation based on the combination of those risk factors. Residents with severe comorbidity and cognitive status and presence of device and wound were identified as certain carriers of USA100 in our sample. Residents with a chronic wound were more likely to carry USA100 MRSA (OR = 2.77, 95% CI = 1.37–5.87). Functional status was identified as an important determinant of carriage of USA100 and USA300 strains. Comorbidity and cognitive status were the two factors associated with carriage of all clonal groups in the study (USA100, USA300 and USA1200).

Conclusions

The combination of Bayesian network analysis, logistic regression and elastic net can be used to identify associations between acquisition and carriage of MRSA strains and their potential risk factors in the face of scarce data. The revealed associations may be used to generate hypothesis for further study of determinants of acquisition and carriage of selected MRSA subtypes and to better inform infection control efforts in community nursing homes.

Peer Review reports

Background

Methicillin-resistant Staphylococcus aureus (MRSA) is a major cause of infection in health-care facilities and the community [1]. Infections due to MRSA are associated with higher morbidity, mortality and costs compared with infections caused by methicillin-susceptible strains [2,3,4]. Colonization with MRSA is known to increase the risk of a subsequent infection in hospital patients and nursing homes residents [5,6,7]. Nursing homes, where MRSA is endemic [8,9,10,11], may play an important role in the regional spread of MRSA [12,13,14]. The emergence of USA300 clone of MRSA in the community and its subsequent spread to healthcare settings in recent years [11, 15, 16] has become a serious problem due to its higher potential for virulence, persistence and transmissibility compared with healthcare-associated strains [17,18,19].

The dynamics of MRSA within nursing homes could be influenced by a number of factors. Studies, many of which were conducted in Veteran Affairs nursing homes or skilled-care facilities, demonstrated that poor functional status, comorbidity, invasive medical devices, chronic wounds, recent antibiotic exposure and hospitalization were risk factors for MRSA colonization in long-term care facilities [9, 10, 20,21,22,23]. However, the impact of these factors on the strain-specific dynamics of MRSA in community nursing homes remains poorly understood. Previously, we used Markov chain models to predict the steady-state distribution of residents colonized with USA300 and non-USA300 MRSA and to assess the impact of potential risk factors on strain-specific acquisition of MRSA in community nursing homes [24]. We found that antibiotic use in the previous 3 months significantly increased acquisition rates of strain-independent MRSA and non-USA300. The effect of antibiotic exposure on the acquisition rate of USA300 and the influence of other potentially modifiable resident characteristics on MRSA acquisition were quite pronounced, but the statistical significance at the 95% confidence level was not achieved. In addition, we employed compartmental and stochastic models in order to evaluate the epidemic potential of strain-specific MRSA and to assess the conditions for MRSA reduction and elimination from community nursing homes [25]. Our results suggested that, while MRSA elimination from nursing homes was unlikely in practice, considerable reductions in MRSA prevalence could be achieved through decolonization therapy that could sustain higher clearance rates over time. Based on our models, antibiotic use in the past 3 months elevated the prevalence of non-USA300 and USA300 MRSA in the facilities, but was unlikely to lead to an outbreak. Moreover, large-scale MRSA outbreaks were not predicted to occur in this setting. However, due to low numbers of observations in some subgroups, especially among subjects colonized with USA300 and exposed to risk factors, our capacity to ascertain whether the candidate risk factors for acquisition and durability of colonization differed by MRSA strain type was limited. Furthermore, while our compartmental and stochastic models assumed that acquisition of one strain was independent from carriage of the other strain, whether this assumption holds in reality is not known.

Bayesian networks have become an increasingly common method for modeling complex and uncertain data in many fields, including medicine [26,27,28]. A Bayesian network is a graphical representation of a joint probability distribution between entities of interest, the number of which can be too large to be modeled successfully using traditional approaches (e.g., regression analysis) [29, 30]. The elastic net (ENET), a new regularization and variable selection method, is another novel approach that gains popularity in modeling scarce data with many predictors [31]. In this study, we sought to use the combination of Bayesian networks, logistic regression and elastic net modeling to identify and quantify the associations between acquisition and carriage of MRSA strains and potentially modifiable resident characteristics in community nursing homes. The identified associations would aid in determining resident characteristics suggestive of a higher risk for acquisition or carriage of specific MRSA strains. The specific aims of this study were: (1) to identify subtypes of MRSA and candidate risk factors associated with acquisition and carriage of the pathogens using Bayesian network analysis; (2) to quantify the strength of the associations between acquisition and carriage of the MRSA subtypes and their potential risk factors using the combination of Bayesian network, elastic net and mixed effects logistic regression modeling.

Methods

Overview

A Bayesian network approach was used to examine the associations between candidate risk factors and acquisition and carriage of MRSA strains discriminated at the 80% and 95% strain similarity thresholds. Mixed effects logistic regression modeling and elastic net were subsequently applied to quantify the strength of the associations revealed by the network. All of the analyses were performed in R version 3.1.2 or higher [32]. The study was reviewed and approved by the Health Sciences Institutional Review Board of the University of Wisconsin-Madison.

Data

The data were collected for a prospective longitudinal study of MRSA colonization in six community nursing homes (size, ≥ 60 beds) in south-central Wisconsin between February 2008 and October 2010 [33]. Four hundred and forty nine of the 851 residents approached (53%) provided written informed consent and participated in the study. Facility-level characteristics and MRSA trends in the study facilities have been previously described [24, 33]. The study subjects were screened for MRSA colonization at multiple anatomical locations at baseline and every three months for up to one year, provided the subjects remained in the facility. Surface cultures of the nares, skin of the axilla, groin and peri-rectal region were obtained using sterile Dacron-tipped swabs. Additional surface cultures were collected from open wounds and the insertion sites of non-urinary invasive devices, when applicable; urine specimens were obtained from subjects with indwelling urinary devices [24, 33]. MRSA specimens were enriched in trypticase soy broth supplemented with 6.5% NaCl and allowed to incubate for 24 h before plating onto Mannitol Salt agar (Remel, Lenexa, KS) containing cefoxitin (4 μg/mL) [33]. The methods used to construct the data set employed in this study have been previously described [24]. Subject exposures that might impact their colonization status were also obtained. These candidate risk factors were characterized as static, that is, ascertained at baseline only, and time-varying which were collected every 3 months. The static risk factors comprised comorbidity (Comorb), functional status (Func) and cognitive status (Cogn). For the purpose of this study, they were dichotomized into non-severe (coded as 0) and severe (coded as 1). The subjects with Charlson Comorbidity Index [34] score ≥ 3, Katz Activities of Daily Living (ADL) [35] score < 2 and Minimum Data Set (MDS) Cognitive Performance Scale (CPS) [36] score ≥ 5 were classified at a severe level of the corresponding factor, and at a non-severe level otherwise. Time-varying risk factors included antibiotic exposure within the previous 3 months (AB), hospitalizations in the previous 3 months (Hosp), presence of a chronic wound (Wnd) and presence of an invasive medical device such as indwelling urinary catheter, percutaneous feeding tube, central venous catheter, or tracheostomy (Dev). The time-varying risk factors were also used as dichotomous variables in our study (Non-exposed was coded as 0 and Exposed was coded as 1). The data from the six study nursing homes were aggregated for the analysis, so that the findings of our study would be representative of the hypothetical “average” nursing home in Wisconsin.

For the purpose of this study, MRSA strains recovered from the study subjects were considered genetically distinct at the 80% strain similarity threshold if their pulsed-field gel electrophoresis (PFGE) banding patterns differed by 4–6 bands [37]. Seven unique strains at the 80% similarity threshold were identified in our study. For each of them, two dichotomous outcome variables were created. These variables indicated the occurrence of strain acquisition (transition from non-colonized to colonized) and carriage (continuous colonization) events. The variables were denoted by T (for acquisition) and C (for carriage) followed by 2 digits indicative of respective MRSA clonal groups in agreement with the CDC classification (i.e., “01” represented USA100 MRSA, “03” specified USA300, and “12” stood for USA1200 clonal group). For example, outcome variables T01 and C01 were created to denote acquisition and carriage of USA100, respectively. For each strain, an acquisition event was considered to occur at time t > 0 if the strain was recovered at time t but not at time t-1. A strain-specific carriage event was assumed to occur at time t = 0 if the strain was recovered during the baseline examination. For subsequent examinations, a carriage event for a strain was considered to occur at time t if the subject was colonized with the strain at time t-1 and remained colonized with the same strain at time t. For each outcome variable, the event occurrence was coded as 1 and non-occurrence as 0. This data set was reduced by removing strain-specific acquisition and carriage variables for which the event occurred in less than 1% of the subjects (this corresponded to fewer than 4 event occurrences in the data set). The reduced data set retained acquisition variables for USA100 (T01) and USA300 (T03), and carriage variables for USA100 (C01), USA300 (C03) and USA1200 (C12). The counts of these events per facility and risk factor exposure are presented in Additional file 1: Table S1.

Notably, acquisition and carriage events of closely-related strains may get different designations when considered at different discriminatory thresholds. For example, if a subject carried strain X from USA100 clonal group at time t, but only strain Y from the same clonal group was recovered at time t + 1, this event at time t + 1 would be classified as carriage at the 80% strain similarity threshold and as acquisition at a higher discriminatory threshold. To assess the impact of this phenomenon on the associations between strain-specific acquisition and carriage events and candidate risk factors, the analysis of the strains discriminated at the 95% similarity threshold was also performed. In this study, MRSA isolates were considered genetically distinct at the 95% strain similarity thresholds if their PFGE banding patterns differed by 1–3 bands [37]. A total of 75 unique strains at the 95% similarity threshold were identified in the study. As done with the isolates at the 80% strain similarity threshold, two dichotomous outcome variables were created for each unique MRSA isolate identified at the 95% threshold. The variables were coded as T (for acquisition) or C (for carriage) followed by a 4-digit sequence. The first two digits were indicative of the respective MRSA clonal group, and the other two digits were assigned based on the genetic similarity to isolates maintained in the University of Wisconsin Infectious Disease Research Laboratory [37]. For example, outcome variables T0101 and C0101 were created to denote acquisition and carriage of the strain from USA100 clonal group. Likewise, the data set was reduced by removing the outcome variables that occurred in less than 1% of the subjects. The reduced data set contained 11 strain-specific acquisition variables with 4–12 event occurrences and 16 carriage variables with 4–38 event occurrences (Additional file 2: Table S2).

Associations at the 80% similarity threshold

To assess probabilistic dependencies between strain-specific acquisition and carriage events and potential risk factors, we used Bayesian network models. Bayesian networks are probabilistic graphical models that represent conditional dependencies between a set of random variables using directed acyclic graphs [29]. In Bayesian networks, random variables are represented by nodes (or vertices), and the conditional dependencies are shown by directed arcs that connect the nodes. For structure learning, we used model averaging to build a network containing only significant arcs [29]. The averaged network was constructed by applying bootstrap resampling to the data to learn a set of 1000 network structures. The arcs that appeared with higher frequency than the significance threshold estimated from the graphs were retained in the averaged network. The graph structure of each of the 1000 networks was learned with the score-based hill-climbing algorithm [29]. In this algorithm, the goodness of fit is assessed by the network score. The search begins with an empty network, and arcs are added, removed or reversed one at a time until the network score is no longer improving. Akaike information criterion (AIC) was used as a score function in the algorithm. The established network structure was used to estimate conditional probabilities of the strain-specific acquisition and carriage events. Candidate risk factors closely related to the strain-specific acquisition and carriage events were identified from the network with the so-called Markov blankets. The Markov blanket of a node in a Bayesian network comprises its parent nodes, children nodes and all other nodes that share a child with this node [29] (e.g., in Fig. 1, Func is the parent node of C03, Dev and Cogn are the children nodes of C03, and the nodes AB, Comorb, Func, and Hosp share a child with C03). The Markov blanket of each strain-specific acquisition and carriage event was determined from the averaged network and was used to quantify the strength of the associations between MRSA strains and potential risk factors. Focusing further analysis on the Markov blankets of the strain-specific events reduced the dimension of the model by removing non-informative factors [38]. The bnlearn package in R version 3.1.2 or higher [32] was used to fit Bayesian networks and to determine Markov blankets.

Fig. 1
figure1

Averaged Bayesian network that includes all potential risk factors and events for MRSA clonal groups. The network was built by averaging 1000 networks learned from bootstrap resampling of the data. The averaged network included only significant arcs. T and C followed by a 2-digit sequence indicate acquisition and carriage events for the associated strain, respectively (e.g., T01 and C01 denote acquisition and carriage for USA100, T03 and C03 for USA300). AB, antibiotic use in the previous 3 months (0 = Non-exposed, 1 = Exposed); Hosp, hospitalizations in the previous 3 months (0 = Non-exposed, 0 = Exposed); Dev, invasive device (0 = Non-exposed, 1 = Exposed); Wnd, wound (0 = Non-exposed, 1 = Exposed); Comorb, comorbidities (0 = Non-severe, 1 = Severe); Func, functional status (0 = Non-severe, 1 = Severe); Cogn, cognitive status (0 = Non-severe, 1 = Severe)

Full size image

To quantify the strength of the association between strain-specific acquisition and carriage events and candidate risk factors, mixed effects logistic regression (MELR) models with logit link function were fitted. The outcome variables were acquisition and carriage events for MRSA clonal groups. Specifically, those were acquisition and carriage of USA100 (T01 and C01, respectively), acquisition and carriage of USA300 (T03 and C03, respectively), and carriage of USA1200 (C12). The independent variables in each original model were the Markov blankets of the corresponding outcome variable (the network is shown in Fig. 1). The models accounted for variability between the facilities by including the facilities variable, Facility, as a random effects term for the intercept. Different observations on the same resident were assumed independent in our study, where 128 subjects (29%) had only baseline observation [24]. The likelihood ratio test was used in the model building, and statistical significance was declared at 95% confidence level. Multilevel bootstrapping with 5000 replications was performed to derive 95% percentile confidence intervals (CI) for the estimated model coefficients. The coefficients of the significant predictors and their bootstrapped 95% CI’s were exponentiated to obtain odds ratios (OR) and their 95% CI’s. Exponentiated model intercepts denoted the odds of the outcome occurrence in the reference group (that is, when all predictors in the model were fixed at 0). The R (version 3.1.2 or higher) package lmer was used for model fitting [32].

In addition, the associations between strain-specific acquisition and carriage events at the 80% similarity and their candidate risk factors were quantified using the elastic net modeling approach. The elastic net is a regularization and variable selection method that produces sparse models with good prediction accuracy when the number of predictors is much higher than the number of observations [31]. In our ENET models, the strain-specific acquisition and carriage events served as outcome variables, while the covariates included the Markov blankets of the respective outcome variables and Facility. Cross-validation with 485-folds and α = 0.5 (the parameter controlling the penalty) was used to find the optimal tuning parameter λ. The number of folds in cross-validation was chosen to be about a third of the number of observations. The model coefficients were estimated based on the optimal tuning parameter at which the minimal mean-squared error (MSE) was achieved. The ENET models were fitted by using the glmnet package in R version 3.1.2 or higher [32].

Associations at 95% similarity threshold

To obtain a preliminary understanding of the overall dependence structure between strain-specific transition and carriage events and candidate risk factors, we constructed a minimal BIC forest. The minimal BIC forest is a graphical model that represents the optimal forest for the data in which the penalized likelihood criterion, Bayesian Information Criterion (BIC), is optimized [39]. The nodes of the forest denote random variables. The edges between the nodes identify variables that are conditionally dependent given the other variables in the model. Likewise, the absence of edges between nodes identifies conditionally independent variables given the other variables. The minimal BIC forest was found and plotted by using the gRapHD package in R version 3.1.2 or higher [32]. Each vertex in the plot was determined by using the algorithm proposed by Fruchterman & Reingold [40] with 5000 iterations. To further our understanding of the associations between the potential risk factors and strain-specific acquisition and carriage events for MRSA strains at the 95% similarity threshold, we built an averaged Bayesian network, similarly to the network that involved clonal groups.

In order to explore the associations between each individual potential risk factor and strain-specific acquisition and carriage events while attenuating the impact of other risk factors, we considered 7 Bayesian networks, each of which included a single potential risk factor and all acquisition and carriage events. Similarly to constructing the Bayesian network that included all potential risk factors, each of the 7 networks was built by applying bootstrap resampling to data to learn a set of 1000 network structures and building an averaged network that included only significant arcs. In each averaged network, the Markov blanket of the potential risk factor was determined. To quantify the strength of the association between each potential risk factor and the strain-specific events that formed its Markov blanket, we attempted to employ logistic regression. However, logistic regression models could not be fitted successfully due to small number of observations and relatively large number of predictors. Elastic net models were fitted instead. In those models, the potential risk factor was considered a dependent variable, while the independent variables were strain-specific acquisition and carriage events that formed the Markov blanket of the risk factor, and Facility. Similarly, to derive the optimal tuning parameter λ, cross-validation with 485-folds and α = 0.5 was performed. The optimal tuning parameter at which the minimal mean-squared error (MSE) was achieved was employed to estimate model coefficients. The elastic net models were fit by using the glmnet package in R version 3.1.2 or higher [32].

Results

Associations at the 80% similarity threshold

Our Bayesian network revealed common patterns of conditional dependencies between acquisition and carriage events of MRSA clonal groups and their candidate risk factors (Fig. 1). The conditional probabilities of carriage of USA100, C01, estimated from the network are presented in Fig. 2 (the conditions associated with a zero probability are not shown). None of the conditional probabilities of occurrence of other events (T01, T03, C03, C12) exceeded 0.05 (data not shown). Thus, a resident with severe comorbidity and cognitive status who had an invasive device and chronic wound was certain to carry USA100 in our sample. Absence of severe cognitive impairment as well as presence of a wound, device and severe comorbidity influenced the probability of USA100 carriage although there was a variation based on the combination of those risk factors. For example, a resident with non-severe comorbidity and cognitive status who had a chronic wound but not a device was estimated to be a carrier of USA100 with probability of 0.42 in our sample.

Fig. 2
figure2

Conditional probabilities of carriage of USA100 derived from the Bayesian network depicted in Fig. 1. Conditions that correspond to a probability of 0 are not displayed. The secondary axis shows the total number of observations with respective combinations of risk factors. Pr, estimates of conditional probabilities; # Obs., the total number of observations with respective combinations of risk factors; Comorb, comorbidity (0 = Non-severe, 1 = Severe); Cogn, cognitive status (0 = Non-severe, 1 = Severe); Dev, presence of invasive device (0 = Non-exposed, 1 = Exposed); Wnd, presence of wound (0 = Non-exposed, 1 = Exposed)

Full size image

This Bayesian network (Fig. 1) was also used to identify Markov blankets of the strain-specific acquisition and carriage events at the 80% strain similarity threshold. The Markov blankets for T01 and T03 consisted of a single factor, C01 and Comorb, respectively. The Markov blanket for C01 comprised T01, Comorb, Cogn, Dev and Wnd, while AB, Comorb, Cogn, Dev, Func, and Hosp formed the Markov blanket for C03. Lastly, the Markov blanket for C12 included Comorb and Cogn.

Estimates for the significant predictors of acquisition and carriage events for MRSA clonal groups identified by MLER and ENET models from the respective Markov blankets are presented in Table 1. In this table, MLER estimates identify predictors that were statistically significant at the 95% confidence level. While logistic regression suggested statistically significant associations between the strain-specific acquisition and carriage events and their candidate risk factors, most of the 95% confidence intervals were extremely wide. One finding stood out: carriage of USA100 was more likely to occur in residents with a chronic wound than in those without it (OR = 2.77, 95% CI = 1.37–5.87).

Table 1 Coefficients derived from MELR and ENET models for strains discriminated at the 80% similarity threshold
Full size table

In elastic net models, presence of chronic wound or invasive device and severe comorbidity were positively associated with carriage of USA100, while severe cognitive status appeared to be a protective factor. Notably, presence of wound showed a stronger association with carriage of USA100 compared to other factors. Among the factors positively associated with carriage of USA300, functional status appeared to be more pronounced. Carriage of USA300 was negatively associated with hospitalizations in the past 3 months. Acquisition of neither strain was associated with carriage of another strain.

Associations at the 95% similarity threshold

A preliminary dependence structure between strain-specific acquisition and carriage events and their candidate risk factors revealed by the minimal BIC forest is shown in Additional file 3: Figure S1. The graph suggests that the candidate risk factors play an important role in the spread and persistence of most MRSA strains discriminated at the 95% similarity threshold. A Bayesian network that shows conditional dependencies between strain-specific acquisition and carriage events of MRSA at the 95% similarity threshold and their potential risk factors is presented in Fig. 3. The Bayesian network and minimal BIC forest revealed similarities in the network structures. Moreover, the Bayesian networks for both 80% and 95% similarly thresholds appear qualitatively similar (Fig. 1 and Fig. 3, respectively).

Fig. 3
figure3

Averaged Bayesian network of potential risk factors and strain-specific acquisition and carriage events. The network was built by averaging 1000 networks learned from bootstrap resampling of the data. The averaged network included only significant arcs. T and C followed by a 4-digit sequence indicate acquisition and carriage events for the associated strains discriminated at the 95% similarity threshold, respectively. AB, antibiotic use in the previous 3 months (0 = Non-exposed, 1 = Exposed); Hosp, hospitalizations in the previous 3 months (0 = Non-exposed, 1 = Exposed); Dev, invasive device (0 = Non-exposed, 1 = Exposed); Wnd, wound (0 = Non-exposed, 1 = Exposed); Comorb, comorbidities (0 = Non-severe, 1 = Severe); Func, functional status (0 = Non-severe, 1 = Severe); Cogn, cognitive status (0 = Non-severe, 1 = Severe)

Full size image

An example of the Bayesian network that involved a single candidate risk factor (AB) and strains discriminated at the 95% similarity threshold is presented in Additional file 4: Figure S2. Elastic net analysis that was used to identify associations between strain-specific acquisition and carriage events and a single candidate risk factor resulted in the selection of covariates and their coefficients shown in Table 2. The total of 5 covariates was not selected by the ENET models from the Markov blankets of the candidate risk factors (C0104 and T0113 in relation to AB, C0114 in relation to Hosp, T0157 in relation to Dev, and T0100 in relation to Wnd). Hosp was associated with acquisition and carriage events of strains from USA100 clonal group only, while other risk factors were associated with the events for strains from both USA100 and USA300 clonal groups. Carriage of the strain from USA1200 clonal group, C1201, was negatively associated with AB and positively associated with Comorb, Func, and Cogn.

Table 2 Coefficients of ENET models for strains discriminated at the 95% strain similarity threshold
Full size table

Discussion

Our study used the combination of Bayesian network analysis, logistic regression and elastic net modeling approaches to determine and quantify the associations between the acquisition and carriage of MRSA strains and their potential risk factors in community nursing homes in Wisconsin. The Bayesian networks that considered strains at the 80% and 95% similarity thresholds revealed a qualitatively similar structure of conditional dependencies between strain-specific acquisition and carriage events and potentially modifiable resident characteristics. To our knowledge, it is the first study that used predictive models to investigate the relation between acquisition and carriage of strain-specific MRSA and their determinants in community nursing homes.

Our results indicated that residents with a chronic wound were more likely to carry USA100 MRSA compared with residents free of wounds. This finding may suggest the presence of a chronic wound as a significant predictor of carriage of USA100 MRSA. The presence of wound was found to be a risk factor for MRSA carriage in other studies [41]. It may also be probable, however, that carriage of USA100 increases the risk for developing a chronic wound. More research is needed to study the direction of the association between carriage of USA100 MRSA and the presence of chronic wounds. Functional status appeared to be an important determinant for carriage of USA300 and USA100 strains in our study. This is in agreement with other studies conducted in nursing homes that identified functional status as a risk factor for MRSA colonization as well [22, 42]. The conditional dependencies revealed by the Bayesian networks highlight the associations between potentially modifiable resident characteristics and strain-specific acquisition and carriage of MRSA. Thus, the knowledge about resident cognitive status, comorbidity, presence of wound and device informs the likelihood of carriage of USA100 in our network. One of the important utilities of this approach can be informing targeted screening for MRSA in nursing homes (i.e., screening residents who are at a higher risk for carrying MRSA). Many hospitals are now pursuing targeted screening [43] in place of universal screening to reduce the costs. More research is needed to establish risk factors for MRSA colonization in nursing homes.

The limitations of our study stem largely from the scarcity of data and from modeling assumptions. The small number of observations available for strain-specific acquisition and carriage events, especially at the higher discriminatory threshold of 95%, coupled with scarcity of events for most risk factors, made some of the association estimates highly uncertain. To improve the predictive power of the models, the data from the six facilities were combined for the analysis. Our results therefore represent the hypothetical “average” nursing home. Combining the data did not allow us to describe the dynamics of MRSA within each facility. We did not account for variation between the residents within each facility due to the limited amount of observations at a resident level. Moreover, the proportion of the participants (53% of the approached residents) in the study nursing homes and subsequent attrition over time may limit the generalizability of the conclusions. Furthermore, the risk factor variables were dichotomized for the purpose of this study. Variable dichotomization was likely to reduce the discriminatory power of the models aimed at detecting the associations between the candidate risk factors and strain-specific events. For example, residents who had a single dose of antibiotics and those with routine antibiotic exposure over the previous 3 months were classified as exposed to antibiotics. Similarly, residents were classified as having severe or non-severe comorbidity, functional status and cognitive status based on the chosen cutoffs for the associated scores. That is, residents with a similar level of exposure to the risk factor of interest could be classified into different levels of this factor. However, we used the best available data to quantify the associations between MRSA strains and candidate risk factors with regards to their acquisition and persistence. More studies are needed to investigate such associations further. While we learned the network structure from the data, incorporating expert knowledge in the network may shed additional light into the determinants of acquisition and carriage of MRSA strains.

Conclusions

We employed Bayesian networks, logistic regression and elastic net modeling approaches to study the associations between acquisition and persistence of MRSA strains and potentially modifiable resident characteristics in community nursing homes. The discovered associations may be used to generate hypotheses regarding risk factors for acquisition and persistence of MRSA strains in this setting for further research. Bigger data sets are needed to test these hypotheses in order to identify risk factors that contribute to acquisition and persistence of strain-specific MRSA in nursing homes.

Abbreviations

CI:

Confidence interval

ENET:

Elastic net

MELR:

Mixed effects logistic regression

MRSA:

Methicillin-resistant Staphylococcus aureus

OR:

Odds ratio

PFGE:

Pulsed-field gel electrophoresis

References

  1. 1.

    Dantes R, Mu Y, Belflower R, Aragon D, Dumyati G, Harrison LH, Lessa FC, Lynfield R, Nadle J, Petit S, et al. National burden of invasive methicillin-resistant Staphylococcus Aureus infections, United States, 2011. JAMA Intern Med. 2013;173(21):1970–8.

  2. 2.

    Capitano B, Leshem OA, Nightingale CH, Nicolau DP. Cost effect of managing methicillin-resistant Staphylococcus Aureus in a long-term care facility. J Am Geriatr Soc. 2003;51(1):10–6.

  3. 3.

    Cosgrove SE, Qi Y, Kaye KS, Harbarth S, Karchmer AW, Carmeli Y. The impact of methicillin resistance in Staphylococcus Aureus bacteremia on patient outcomes: mortality, length of stay, and hospital charges. Infect Control Hosp Epidemiol. 2005;26(2):166–74.

  4. 4.

    Anderson DJ, Kaye KS, Chen LF, Schmader KE, Choi Y, Sloane R, Sexton DJ. Clinical and financial outcomes due to methicillin resistant Staphylococcus Aureus surgical site infection: a multi-center matched outcomes study. PLoS One. 2009;4(12):e8305.

  5. 5.

    Muder RR, Brennen C, Wagener MM, Vickers RM, Rihs JD, Hancock GA, Yee YC, Miller JM, VL Y. Methicillin-resistant staphylococcal colonization and infection in a long-term care facility. Ann Intern Med. 1991;114(2):107–12.

  6. 6.

    Shurland SM, Stine OC, Venezia RA, Zhan M, Furuno JP, Miller RR, Roghmann MC. USA300 methicillin-resistant S. Aureus (USA300 MRSA) colonization and the risk of MRSA infection in residents of extended-care facilities. Epidemiol Infect. 2012;140(3):390–9.

  7. 7.

    Fukuta Y, Cunningham CA, Harris PL, Wagener MM, Muder RR. Identifying the risk factors for hospital-acquired methicillin-resistant Staphylococcus Aureus (MRSA) infection among patients colonized with MRSA on admission. Infect Control Hosp Epidemiol. 2012;33(12):1219–25.

  8. 8.

    Reynolds C, Quan V, Kim D, Peterson E, Dunn J, Whealon M, Terpstra L, Meyers H, Cheung M, Lee B, et al. Methicillin-resistant Staphylococcus Aureus (MRSA) carriage in 10 nursing homes in Orange County, California. Infect Control Hosp Epidemiol. 2011;32(1):91–3.

  9. 9.

    Mody L, Kauffman CA, Donabedian S, Zervos M, Bradley SF. Epidemiology of Staphylococcus Aureus colonization in nursing home residents. Clin Infect Dis. 2008;46(9):1368–73.

  10. 10.

    Garazi M, Edwards B, Caccavale D, Auerbach C, Wolf-Klein G. Nursing homes as reservoirs of MRSA: myth or reality? J Am Med Dir Assoc. 2009;10(6):414–8.

  11. 11.

    Murphy CR, Hudson LO, Spratt BG, Quan V, Kim D, Peterson E, Tan G, Evans K, Meyers H, Cheung M, et al. Predicting high prevalence of community methicillin-resistant Staphylococcus Aureus strains in nursing homes. Infect Control Hosp Epidemiol. 2013;34(3):325–6.

  12. 12.

    Lee BY, Bartsch SM, Wong KF, Singh A, Avery TR, Kim DS, Brown ST, Murphy CR, Yilmaz SL, Potter MA, 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.

  13. 13.

    Macal CM, North MJ, Collier N, Dukic VM, Wegener DT, David MZ, Daum RS, Schumm P, Evans JA, Wilder JR, et al. Modeling the transmission of community-associated methicillin-resistant Staphylococcus Aureus: a dynamic agent-based simulation. J Transl Med. 2014;12:124.

  14. 14.

    Barnes SL, Harris AD, Golden BL, Wasil EA, Furuno JP. Contribution of interfacility patient movement to overall methicillin-resistant Staphylococcus Aureus prevalence levels. Infect Control Hosp Epidemiol. 2011;32(11):1073–8.

  15. 15.

    Tattevin P, Diep BA, Jula M, Perdreau-Remington F. Methicillin-resistant Staphylococcus Aureus USA300 clone in long-term care facility. Emerg Infect Dis. 2009;15(6):953–5.

  16. 16.

    Diekema DJ, Richter SS, Heilmann KP, Dohrn CL, Riahi F, Tendolkar S, McDanel JS, Doern GV. Continued emergence of USA300 methicillin-resistant Staphylococcus Aureus in the United States: results from a nationwide surveillance study. Infect Control Hosp Epidemiol. 2014;35(3):285–92.

  17. 17.

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

  18. 18.

    Desai R, Pannaraj PS, Agopian J, Sugar CA, Liu GY, Miller LG. Survival and transmission of community-associated methicillin-resistant Staphylococcus Aureus from fomites. Am J Infect Control. 2011;39(3):219–25.

  19. 19.

    Crnich CJ. Impact and management of MRSA in the long-term care setting. Curr Transl Geriatr and Exp GErontol Rep. 2013;2:125–35.

  20. 20.

    Stone ND, Lewis DR, Johnson TM, 2nd, Hartney T, Chandler D, Byrd-Sellers J, McGowan JE, Jr., Tenover FC, Jernigan JA, Gaynes RP: Methicillin-resistant Staphylococcus Aureus (MRSA) nasal carriage in residents of veterans affairs long-term care facilities: role of antimicrobial exposure and MRSA acquisition. Infect Control Hosp Epidemiol 2012, 33(6):551–557.

  21. 21.

    Murphy CR, Quan V, Kim D, Peterson E, Whealon M, Tan G, Evans K, Meyers H, Cheung M, Lee BY, et al. Nursing home characteristics associated with methicillin-resistant Staphylococcus Aureus (MRSA) burden and transmission. BMC Infect Dis. 2012;12:269.

  22. 22.

    Bradley SF, Terpenning MS, Ramsey MA, Zarins LT, Jorgensen KA, Sottile WS, Schaberg DR, Kauffman CA. Methicillin-resistant Staphylococcus Aureus: colonization and infection in a long-term care facility. Ann Intern Med. 1991;115(6):417–22.

  23. 23.

    Trick WE, Weinstein RA, DeMarais PL, Kuehnert MJ, Tomaska W, Nathan C, Rice TW, McAllister SK, Carson LA, Jarvis WR. Colonization of skilled-care facility residents with antimicrobial-resistant pathogens. J Am Geriatr Soc. 2001;49(3):270–6.

  24. 24.

    Batina NG, Crnich CJ, Anderson DF, Döpfer D. Models to predict prevalence and transition dynamics of methicillin-resistant Staphylococcus Aureus in community nursing homes. Am J Infect Control. 2016;44(5):507–14.

  25. 25.

    Batina NG, Crnich CJ, Anderson DF, Dopfer D. Identifying conditions for elimination and epidemic potential of methicillin-resistant Staphylococcus Aureus in nursing homes. Antimicrob Resist Infect Control. 2016;5:32.

  26. 26.

    Gomez SA, Kulow M, Anklam KS, Park D, Kaspar CW, Ivanek R, Dopfer D. Gene markers of generic Escherichia Coli associated with colonization and persistence of Escherichia Coli O157 in cattle. Prev Vet Med. 2014;117(1):140–8.

  27. 27.

    Kahn CE Jr, Roberts LM, Shaffer KA, Haddawy P. Construction of a Bayesian network for mammographic diagnosis of breast cancer. Comput Biol Med. 1997;27(1):19–29.

  28. 28.

    O'Brien KM, Cole SR, Engel LS, Bensen JT, Poole C, Herring AH, Millikan RC. Breast cancer subtypes and previously established genetic risk factors: a bayesian approach. Cancer Epidemiol Biomarkers Prev. 2014;23(1):84–97.

  29. 29.

    Nagarajan R, Scutari M, Lèbre S. Bayesian networks in R with applications in systems biology. New York, NY: Springer; 2013.

  30. 30.

    Darwiche A: Modeling and reasoning with Bayesian networks. In. Cambridge; New York: Cambridge University Press; 2010: 1 online resource.

  31. 31.

    Zou H, Hastie T: Regularization and variable selection via the elastic net. Journal of the Royal Statistical Society 2005, 67, Part 2:301–320.

  32. 32.

    R Core Team: R: A language and environment for statistical computing. In.: R Foundation for Statistical Computing, Vienna, Austria. URL https://www.r-project.org/. 2014.

  33. 33.

    Crnich CJ, Duster M, Hess T, Zimmerman DR, Drinka P. Antibiotic resistance in non-major metropolitan skilled nursing facilities: prevalence and interfacility variation. Infect Control Hosp Epidemiol. 2012;33(11):1172–4.

  34. 34.

    Charlson ME, Pompei P, Ales KL, MacKenzie CR. A new method of classifying prognostic comorbidity in longitudinal studies: development and validation. J Chronic Dis. 1987;40(5):373–83.

  35. 35.

    Katz S, Ford AB, Moskowitz RW, Jackson BA, Jaffe MW. Studies of illness in the aged. The index of Adl: a standardized measure of biological and psychosocial function. JAMA. 1963;185:914–9.

  36. 36.

    Morris JN, Fries BE, Mehr DR, Hawes C, Phillips C, Mor V, Lipsitz LA, Cognitive Performance MDS. Scale. J Gerontol. 1994;49(4):M174–82.

  37. 37.

    Crnich CJ, Duster M, Warrack S, Maki D, Safdar N. Comparison of pulsed-gel electrophoresis and a commercial repetitive-element PCR method for assessment of methicillin-resistant Staphylococcus Aureus clustering in different health care facilities. J Clin Microbiol. 2014;52(6):2027–32.

  38. 38.

    Scutari M, Mackay I, Balding D. Improving the efficiency of genomic selection. Stat Appl Genet Mol Biol. 2013;12(4):517–27.

  39. 39.

    Edwards D, de Abreu GC, Labouriau R. Selecting high-dimensional mixed graphical models using minimal AIC or BIC forests. BMC bioinformatics. 2010;11:18.

  40. 40.

    Fruchterman TMJ, Reingold EM: Graph Drawing by Force-directed Placement. Software: Practice and Experience 1991, 21(11)(11):1129–1164.

  41. 41.

    von Baum H, Schmidt C, Svoboda D, Bock-Hensley O, Wendt C. Risk factors for methicillin-resistant Staphylococcus Aureus carriage in residents of German nursing homes. Infect Control Hosp Epidemiol. 2002;23(9):511–5.

  42. 42.

    Furuno JP, Shurland SM, Zhan M, Johnson JK, Venezia RA, Harris AD, Roghmann MC. Comparison of the methicillin-resistant Staphylococcus Aureus acquisition among rehabilitation and nursing home residents. Infect Control Hosp Epidemiol. 2011;32(3):244–9.

  43. 43.

    Furuno JP, McGregor JC, Harris AD, Johnson JA, Johnson JK, Langenberg P, Venezia RA, Finkelstein J, Smith DL, Strauss SM, et al. Identifying groups at high risk for carriage of antibiotic-resistant bacteria. Arch Intern Med. 2006;166(5):580–5.

Download references

Acknowledgements

We would like to express our great appreciation to Drs. Vicki Bier, Molly Carnes and David Zimmerman for their constructive suggestions, invaluable feedback, and support during different stages of the study. We would like to thank Ms. Helena Tsotsis for her help with data retrieval.

Funding

Christopher J. Crnich activity on this project was supported by the Clinical and Translational Science Award (CTSA) program, through the NIH National Center for Advancing Translational Sciences (NCATS) (Grant UL1TR000427; the content is solely the responsibility of the authors and does not necessarily represent the official views of the NIH); the University of Wisconsin School of Medicine and Public Health, through the Wisconsin Partnership Program; and a research scholarship from the ASP/IDSA T. Franklin Williams Program.

Availability of data and materials

The dataset analyzed during this study is not publicly available to protect data confidentiality and subjects’ privacy. A limited dataset may be available from the corresponding author on reasonable request.

Author’s contributions

NGB contributed to the conception and design of the study, performed mathematical modeling, interpreted the results, and wrote the manuscript. CJC provided the data for this study, contributed to the conception and design of the study, interpretation of the results and their implications to the field, and oversaw the revision of the manuscript. DD conceived and designed the study, and supervised modeling effort, interpretation of the results and writing the manuscript. All authors reviewed the manuscript and approved the final version of the manuscript for publication.

Author information

Affiliations

  1. Department of Industrial and Systems Engineering, University of Wisconsin-Madison, 3270 Mechanical Engineering Building, 1513 University Avenue, Madison, WI, 53706, USA
    • Nataliya G. Batina
  2. Department of Medicine, University of Wisconsin-Madison, 2500 Overlook Terrace, B5112E, Madison, WI, 53705, USA
    • Christopher J. Crnich
  3. William S. Middleton Veterans Affairs Hospital, 2500 Overlook Terrace, B5112E, Madison, WI, 53705, USA
    • Christopher J. Crnich
  4. Department of Medical Sciences, School of Veterinary Medicine, University of Wisconsin-Madison, 2027 Veterinary Medicine Building, 2015 Linden Dr, Madison, WI, 53706, USA
    • Dörte Döpfer
Authors
  1. Nataliya G. Batina
    View author publications
    You can also search for this author in
  2. Christopher J. Crnich
    View author publications
    You can also search for this author in
  3. Dörte Döpfer
    View author publications
    You can also search for this author in

Corresponding author

Correspondence to Nataliya G. Batina.

Ethics declarations

Ethics approval and consent to participate

The study was reviewed and approved by the Health Sciences Institutional Review Board of the University of Wisconsin-Madison (study ID 2012–0382). All study participants provided written informed consent to participate in the study.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

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

Additional files

Additional file 1: Table S1.

Counts of event occurrences per facility and exposure to potential risk factors. (DOCX 35 kb)

Additional file 2: Table S2.

Counts of acquisition and carriage for strains identified at the 95% strain similarity threshold. (DOCX 31 kb)

Additional file 3: Figure S1.

Minimal BIC Forest for MRSA strains distinguished at 95% strain similarity threshold. The network shows the conditional dependence structure between strain-specific acquisition and carriage events and the candidate risk factors. (DOCX 15 kb)

Additional file 4: Figure S2.

Averaged Bayesian network that includes recent antibiotic use and strain-specific acquisition and carriage events. The network was built by averaging 500 networks learned from bootstrap resampling of the data. The averaged network included only significant arcs. T and C followed by a 4-digit sequence indicate acquisition and carriage events for the associated strains discriminated at the 95% similarity threshold, respectively. AB, antibiotic use in the past 3 months (0 = Non-exposed, 1 = Exposed). (DOCX 17 kb)

Rights and permissions

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

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Batina, N.G., Crnich, C.J. & Döpfer, D. Acquisition and persistence of strain-specific methicillin-resistant Staphylococcus aureus and their determinants in community nursing homes. BMC Infect Dis 17, 752 (2017). https://doi.org/10.1186/s12879-017-2837-3

Download citation

Keywords

  • USA300 MRSA
  • Non-USA300 MRSA
  • Colonization
  • Risk factors
  • Acquisition
  • Persistence
  • Nursing home

BMC Infectious Diseases

ISSN: 1471-2334

Contact us