Skip Navigation

Spread of Carbapenem-Resistant Enterobacteriaceae Among Illinois Healthcare Facilities: The Role of Patient Sharing

  1. William E. Trick2,3
  1. 1Division of Patient Safety and Quality, Illinois Department of Public Health
  2. 2Rush University Medical Center
  3. 3Cook County Health and Hospitals System, Chicago, Illinois
  1. Correspondence: M. J. Ray, 122 S Michigan Ave, 7th flr, Chicago, IL 60603 (michael.j.ray{at}illinois.gov).

Abstract

Background. Carbapenem-resistant Enterobacteriaceae (CRE) spread regionally throughout healthcare facilities through patient transfer and cause difficult-to-treat infections. We developed a state-wide patient-sharing matrix and applied social network analyses to determine whether greater connectedness (centrality) to other healthcare facilities and greater patient sharing with long-term acute care hospitals (LTACHs) predicted higher facility CRE rates.

Methods. We combined CRE case information from the Illinois extensively drug-resistant organism registry with measures of centrality calculated from a state-wide hospital discharge dataset to predict facility-level CRE rates, adjusting for hospital size and geographic characteristics.

Results. Higher CRE rates were observed among facilities with greater patient sharing, as measured by degree centrality. Each additional hospital connection (unit of degree) conferred a 6% increase in CRE rate in rural facilities (relative risk [RR] = 1.056; 95% confidence interval [CI], 1.030–1.082) and a 3% increase among Chicagoland and non-Chicago urban facilities (RR = 1.027; 95% CI, 1.002–1.052 and RR = 1.025; 95% CI, 1.002–1.048, respectively). Sharing 4 or more patients with LTACHs was associated with higher CRE rates, but this association may have been due to chance (RR = 2.08; 95% CI, .85–5.08; P = .11).

Conclusions. Hospitals with greater connectedness to other hospitals in a statewide patient-sharing network had higher CRE burden. Centrality had a greater effect on CRE rates in rural counties, which do not have LTACHs. Social network analysis likely identifies hospitals at higher risk of CRE exposure, enabling focused clinical and public health interventions.

Key words

(See the Editorial Commentary by Polgreen and Segre on pages 894–5.)

Extensively drug-resistant organisms (XDROs) are a public health threat because they are nearly untreatable with currently available antibiotics and can spread rapidly throughout a region by interfacility transfer of colonized patients [17]. Carbapenem-resistant Enterobacteriaceae (CRE) are increasingly prevalent XDROs associated with up to 50% mortality in infected persons [8, 9]. Klebsiella pneumoniae carbapenemase is the most common carbapenemase reported in the United States and has become a concern because of its propensity for epidemic spread [10]. Infection control interventions [11, 12] and mathematical models [13] have demonstrated that control of such antimicrobial-resistant pathogens requires regional coordination across healthcare facilities [8, 14].

CRE are endemic in Illinois healthcare facilities, especially among patients of a hospital subtype referred to as long-term acute care hospitals (LTACHs) [1]. LTACHs care for chronically critically ill patients who are at high risk of CRE carriage because of greater severity of illness, longer lengths of stay, higher antimicrobial exposure, and increased facility CRE colonization pressure compared with patients in traditional hospitals [15]. LTACHs have played an apparent role in the regional dissemination of CRE to short-term acute care hospitals [2, 3, 16]. To track CRE, the Illinois Department of Public Health (IDPH) made CRE reportable to the XDRO registry (www.xdro.org) in 2013 [14, 17].

Assessing a given hospital's risk of CRE exposure is critical to targeting surveillance and prevention efforts. Such risk may be driven by patterns of patient movement throughout the healthcare network [13]. The extent to which facilities are linked within a network can be quantified using measures of “centrality,” that is, interconnectedness, derived from social network analysis [18]. Studies have evaluated these connections through simulations; however, relatively few have been able to associate measures of centrality with clinical microbiology results. Degree is a simple count of the number of equally weighted links that a given hospital has with all other hospitals in the network, where a link is defined as 2 hospitals sharing at least 1 patient in the calendar year [18]. Eigenvector centrality is an extension of degree centrality; rather than equally weighting all connections, it gives more weight to connections with more highly connected hospitals [19]. We hypothesized that increased patient sharing with all other hospitals (as defined by degree or eigenvector centrality) and increased patient sharing with LTACHs would each predict higher risk of CRE burden for a given hospital [2022]. We used the XDRO registry and a state-wide hospital admission history database to evaluate the extent that connectedness to other facilities and, in particular, LTACHs influenced rates of CRE.

METHODS

Data Sources

We obtained facility-level CRE case information reported from 1 January 2014 through 31 December 2014 from the Illinois XDRO registry. To minimize surveillance bias, we removed cases only reported from screening cultures (ie, rectal or stool source for the isolate) as not every facility actively screens patients for CRE. Then, we removed duplicate case reports by retaining the earliest CRE report for each patient. Facility-level characteristics such as bed size and number of admissions were obtained from the 2013 IDPH annual hospital questionnaire [23]. The facility list was restricted to short-term acute care hospitals and LTACHs. We focused our evaluation of CRE rates on short-term acute care hospitals.

Social Network Analysis

The social network analysis included all patient admissions in IDPH's hospital discharge dataset during the first quarter of 2013. An adjacency matrix was built in R (www.r-project.org), with each cell representing the number of unique patients who visited both corresponding facilities during the 3-month period. Because individual patients can be colonized for prolonged periods of time and because most interfacility patient sharing is not through direct transfers, we did not require same-day interfacility transfers. We calculated degree and eigenvector centrality using UCINET (Harvard Analytic Technologies, Inc., Harvard, Massachusetts).

Model Specification

The response variable was the number of facility-level CRE cases per 10 000 patient-days. Predictor variables included the 2 centrality measures (degree and eigenvector), number of beds for each facility, number of patients shared with LTACHs, county type (Chicagoland [Cook and adjacent counties], urban non-Chicago, and rural) based on the Illinois Behavioral Risk Factor Surveillance System designations (Figure 1) [24], with rural counties serving as the referent group. Bivariate associations were examined before constructing multivariable models. We used the Kruskal–Wallis nonparametric test on continuous covariates, and calculated Spearman rank correlations to test for ordinal trends. Both centrality measures were categorized into quintiles, and degree centrality was also evaluated as a discrete variable. The number of patients shared with LTACHs was split into the following 5 interval categories: zero, 1–3, 4–6, 7–9, or 10+ patients. These categories were eventually dichotomized as (<4 or ≥4) based on identifying a readily apparent departure from linearity through locally weighted scatterplot smoothing. Number of beds was evaluated as a discrete variable, and indicator variables were coded for county type. We roughly compared degree and eigenvector centrality using the Akaike information criterion (AIC) to see which measure of centrality better predicted CRE rates.

Figure 1.

Illinois county types. Note: County types are based on Illinois Behavioral Risk Factor Surveillance System (BRFSS) designations [24]. We defined Chicagoland as Chicago + suburban Cook County + Collar County BRFSS locations.

To account for overdispersion of values, we used negative binomial regression and contrast statements to generate rate ratios across the covariates. We used an offset of log patient-days, which allowed us to interpret our response variable as a rate. It was decided a priori to adjust for number of beds in the final models to better estimate the impact of patient sharing on CRE rates. We constructed separate models for degree and eigenvector centrality and examined all potential 2-way interactions. All analyses, with the exception of the social network analysis, were performed using SAS software, version 9.4 (SAS Institute, Cary, North Carolina).

RESULTS

Of 185 Illinois hospitals, 99 (54%) reported at least 1 CRE case during the study period. There was a mean of 3.5 CRE cases per hospital, or a crude rate of 1.0 CRE case per 10 000 patient-days (Table 1). There was an average of 64 patient-sharing connections (degree) per facility, with a minimum of 1 connection and a maximum of 145 connections (Table 1). The majority of facilities (63%) reporting CRE cases were within Chicagoland counties; 80 (45%) had shared at least 1 patient with an LTACH, and 39 (21%) shared 4 or more patients with LTACHs. Among CRE cases, the mean age was 64 years. Cases were equally distributed among males (50.1%) and females (49.9%). Among case patients for whom race information was reported to the registry (61% of reports), most were African American (44%) or non-Hispanic white (43%).

Table 1.

Characteristics of Short-Term Acute Care Hospitals That Reported at Least 1 Carbapenem-Resistant Enterobacteriaceae Case to the Illinois Extensively Drug-Resistant Organism Registry During, 2014

The association between degree centrality and CRE rate appeared to be linear (Loess smoothing parameter = 1.0), which indicated that specifying degree centrality as a discrete variable was appropriate. In bivariate analysis, each additional hospital connection, defined by degree, corresponded to an absolute rate increase of approximately 0.09 cases per 10 000 patient-days for each additional connection (Table 2). Sharing patients with LTACHs (4 or more patients shared in 3 months) was associated with an increase in CRE rate. CRE rates were significantly different for each county type, with the highest crude rates in Chicagoland counties (Table 2).

Table 2.

Bivariate Associations Between Facility Characteristics and Carbapenem-Resistant Enterobacteriaceae Rates (per 10 000 Patient-days) Among Illinois Acute Care Hospitals During, 2014

In multivariable modeling adjusted for hospital bed size and LTACH sharing, we found that hospital degree centrality was an independent predictor CRE rate and that the strength of association varied significantly by county type (Table 3). In Chicagoland, each additional hospital connection (defined by degree centrality) was associated with approximately a 3% increase in CRE rate (Table 3); the association was similar in non-Chicago urban counties. The effect of additional connections was twice as strong in rural counties, with about a 6% increase in CRE rate with each additional connection (Table 3). The rate ratio for sharing 4 or more patients with LTACHs strongly suggested an increase rate for a facility's exposure to LTACH patients. However, LTACH sharing and degree are strongly correlated variables, and the association with LTACH sharing may have been due to chance (P = .11). We did find that degree centrality was strongly associated with high LTACH sharing (pooled t-test, P < .0001), which suggests collinearity as a reason for the nonsignificant LTACH variable in our model. LTACH sharing was also significant in our multivariate model when we omitted the centrality variable (P = .0004). Degree and eigenvector centrality predicted CRE equally well according to AIC.

Table 3.

Adjusted Associations Between Hospital Characteristics (Centrality and Long-Term Acute Care Hospital Sharing) and Carbapenem-Resistant Enterobacteriaceae Rates (per 10 000 Patient-days) Among Short-Term Acute Care Hospitals in Illinois

DISCUSSION

Using social network analysis, we evaluated the influence of patient sharing on hospital-specific rates of CRE. We found that greater hospital centrality was independently associated with higher CRE rates in the state-wide patient-sharing network. High LTACH connectedness trended toward higher hospital CRE rates, although the association was not statistically significant after adjusting for centrality; this may have been related to the strong correlation between centrality and LTACH sharing. Based on our results and those reported by other investigations, we believe that increased patient sharing with LTACHs confers an increased risk of CRE exposure. The centrality association was stronger for rural counties than for urban counties. Our results provide an important way of identifying hospitals at highest risk of CRE exposure. These findings have immediate public health implications in that early interventions should be focused on the most connected facilities, as well as those with strong connections to LTACHs. Also, centrality could be used to measure a facility's risk for encountering CRE patients and potentially generate an expected rate of CRE using prediction models as we continue with this work. As such, even single facility outbreaks should prompt early interventions among those facilities that are strongly coupled in terms of patient sharing to try to preempt regional spread, emphasizing the key value of prospectively determining the networks.

Measures of centrality are the key quantitative expression of social network analysis. There are several centrality measures that can be used. Since we could not ascertain directionality with a fully de-identified dataset, we chose to evaluate degree and eigenvector centrality. Both measures reliably predicted an increased CRE rate. Because degree centrality is relatively easy to comprehend and performed as well as eigenvector centrality, we present the results for the degree measure.

During regional outbreaks of multidrug-resistant organism, hospitals routinely perform self-assessments of exposure risk but lack appropriate data to adequately measure such risk. State-wide patient-sharing data, which are now increasingly available through sources like the Healthcare Cost and Utilization Project [25], provide an important way to assess hospital risk of CRE exposure based on its position in regional patient-sharing networks. Studies have shown the utility of identifying the most central facilities in a healthcare network as “sentinel” hospitals, which can facilitate surveillance efforts [26]. A hospital's awareness of its centrality within a patient-sharing network, as well as awareness of its patient sharing with LTACHs, can also provide important information about risk. Knowledge of the case rates at other facilities, coupled with information about the patient-sharing network, could motivate high-risk hospitals to actively test newly admitted patients for CRE colonization and could be used to prioritize the public health response to the highest-risk hospitals. Public health can play a critical role in identifying tightly connected hospitals and educating personnel at such facilities about their risk and need for enhanced infection control interventions. Enhanced interfacility communication can be facilitated by public health–led health information exchanges such as the Illinois XDRO registry [17].

Most connections were between LTACHs and hospitals that have higher CRE case rates and also relatively high centrality. The association with LTACHs is consistent with prior analyses that have shown the central role LTACHs have in the spread of CRE [2]. We found a threshold effect in that a hospital had to share at least 4 or more patients during our analytic time frame (3 months) to have an elevated CRE rate. This threshold makes sense because we found that many, but not all, LTACH patients are colonized with CRE [1]. Therefore, on average, an acute care hospital would need to share multiple patients with LTACHs before encountering a CRE-colonized individual. Once the threshold of 4 incidents of patient sharing with LTACHs was reached, the crude CRE rate doubled.

A major strength of our study is the novelty of using a social network variable as a predictor of CRE case rates. This has been done in some healthcare settings, for example, examination of sexually transmitted infections through networks [2730], but social network measures have been used relatively sparingly in other epidemiological areas [20]. In contrast to prior studies of centrality, we were able to access a full year of clinical CRE data from the Illinois XDRO registry, which provides facility-level CRE rates across the entire state. Access to healthcare use data for acute- and LTACHs allowed us to quantify the magnitude of a hospital's affinity to the network (ie, other facilities in the state). Another strength is that we had access to clinical microbiological data, which allowed us to evaluate the impact of centrality in a real-world setting rather than through simulations.

Our study has several limitations. First, although CRE reports were obtained from a mandatory surveillance system (XDRO registry), it is possible that case reporting was incomplete. Incomplete reporting could bias our findings either toward or away from the null, depending on whether underreporting was random or associated with centrality. Second, we did not have directional social network measures (ie, we know that the same patient visited 2 facilities during the time period but did not know which facility was visited first). Despite this potential limitation, we found that strong connections to the network were associated with a high rate of CRE, and we were able to evaluate patient sharing beyond simple direct transfers. Third, we did not have admissions data for nursing homes, which would have allowed us to quantify associations to the network across the broader continuum of care. Finally, variability in case rates may be related to differences in clinical culture-acquisition practices or seasonal differences in patient transfer patterns.

Combining surveillance data with social network measures such as centrality can guide large-scale interventions by identifying at-risk facilities in the network. Our work shows that by using widely available patient-sharing data, hospitals and public health officials can identify which facilities are at highest risk of CRE exposure. There is urgency for such analyses in order to intervene early and prevent the spread of XDROs.

Notes

Acknowledgments. M. J. R. had full access to all study data and takes responsibility for data integrity and accuracy of the data analysis. We acknowledge Strengthening Health Systems through Interprofessional Education (Project SHINE) and the Illinois Department of Public Health, Division of Patient Safety and Quality. We acknowledge Kruti Doshi, Collaborative Research Unit, Cook County Health and Hospitals System, for her help with social network analysis.

Financial support. This study was funded by the Centers for Disease Control and Prevention (CDC) cooperative agreement U54-CK000161 for the CDC Prevention Epicenter Program.

Potential conflicts of interests. All authors: No reported conflicts. All authors have submitted the ICMJE Form for Disclosure of Potential Conflicts of Interest. Conflicts that the editors consider relevant to the content of the manuscript have been disclosed.

  • Received December 14, 2015.
  • Accepted June 2, 2016.

References

| Table of Contents

Published on behalf of

cid cover HIV Medicine Association Logo

Society Members: For your free access to this journal, log in via the IDSA members area.

Impact Factor: 8.886

5-Yr impact factor: 9.206

For Reviewers

Looking for your next opportunity?

Looking for jobs...

For the Media