The evolution of drug resistant bacteria is a severe public health problem, both in hospitals and in the community. Currently, some countries aim at concentrating highly specialized services in large hospitals in order to improve patient outcomes. Emergent resistant strains often originate in health care facilities, but it is unknown to what extent hospital size affects resistance evolution and the resulting spillover of hospital-associated pathogens to the community. We used two published datasets from the US and Ireland to investigate the effects of hospital size and controlled for several confounders such as antimicrobial usage, sampling frequency, mortality, disinfection and length of stay. The proportion of patients acquiring both sensitive and resistant infections in a hospital strongly correlated with hospital size. Moreover, we observe the same pattern for both the percentage of resistant infections and the increase of hospital-acquired infections over time. One interpretation of this pattern is that chance effects in small hospitals impede the spread of drug-resistance. To investigate to what extent the size distribution of hospitals can directly affect the prevalence of antibiotic resistance, we use a stochastic epidemiological model describing the spread of drug resistance in a hospital setting as well as the interaction between one or several hospitals and the community. We show that the level of drug resistance typically increases with population size: In small hospitals chance effects cause large fluctuations in pathogen population size or even extinctions, both of which impede the acquisition and spread of drug resistance. Finally, we show that indirect transmission via environmental reservoirs can reduce the effect of hospital size because the slow turnover in the environment can prevent extinction of resistant strains. This implies that reducing environmental transmission is especially important in small hospitals, because such a reduction not only reduces overall transmission but might also facilitate the extinction of resistant strains. Overall, our study shows that the distribution of hospital sizes is a crucial factor for the spread of drug resistance.
The last decades have shown that the introduction of an antibiotic agent is almost inevitably followed by the spread of resistance mutations that jeopardize the beneficial effect of this agent [1], [2]. Because of this process of bacterial adaptation to antibiotics, maintaining the benefits of antibiotic therapy requires a steady development of new drugs or drug classes. Population biological models may contribute to slowing down the required pace of this “drug treadmill” by identifying the factors that determine the adaptability of bacterial populations to antimicrobial treatment [3]. The epidemic spread of antibiotic resistance can be strongly affected by the structure of the human host population[4]. One of the most important instances of such population structure is the interaction between the hospital and the community[4]. These two settings differ with respect to several parameters that are crucial for the spread of antibiotic resistance. While the hospital environment is characterized by small population sizes, high transmission rates, fast turn over and frequent use of antibiotics, the community exhibits comparatively large population sizes, small transmission rates, slow turn-over rates and infrequent use of antibiotics. Hospitals are often the source of emergent resistant strains [5], but this spread is not unidirectional, as illustrated by outbreaks of community-acquired MRSA in health care facilities [6]. Since the spread of resistance mutations [7] increases with antibiotic usage [8], the difference in treatment frequencies may explain why hospitals mostly act as source for resistance mutations and the community acts as a sink. The size distribution of hospitals is an important determinant for the population structure generated through the hospital-community interaction. In small hospitals, bacterial population sizes and frequencies are subject to strong stochastic effects and populations may frequently become extinct. It has been empirically shown that resistance levels tend to be lower in small hospitals [9], [10], [11], [12], [13]. In principle, this can be due to two reasons: On the one hand, small hospitals might be associated with different types of patients and treatments (i. e. lower antibiotic usage [14]), which select for less resistance. On the other hand, small hospital size by itself might hinder bacterial adaptation and thereby reduce resistance levels. The fact that the analysis in [10] controlled for patient characteristics, suggests that, at least in that case, the impact of hospital size was due to the second “intrinsic” mechanism. In this study we use a simple population biological model to analyze the intrinsic effects of hospital size, i. e. we assess to what extent the stochastic effects resulting from small hospital populations may help in alleviating the burden of antibiotic resistance.
Environmental transmission in the hospital can be added to the above model via two additional, deterministic compartments corresponding to the sensitive and resistant strain: the density of bacteria of strain i in the environment, denoted Ei. Bacteria of strain i colonize the environment with a rate NH, i cE and are cleared from the environment at a constant rate rTurn, E. Bacteria of strain i from the environment can in turn infect susceptible patients with a force of infection βE Ei. . Because the total number of bacteria in the environment is presumably very large and the dynamics of the environmental compartment are not directly affected by the fluctuations in the patient population, we assume that this compartment can be adequately described deterministically.
We used an Irish surveillance dataset published by the Health Service Executive Ireland [19] to investigate to what extent the incidence of infections, especially by resistant strains, correlates with hospital size. This dataset contains information from 53 hospitals about both the total number of new infections with S. aureus as well as infections with MRSA (all positive blood-cultures were recorded). Additionally reported quantities were: length of stay, total inpatient antibiotic usage, injectable inpatient antibiotic usage, the usage of hospital-specific antibiotics, consumption of alcohol hand-gel and the frequency of blood cultures per admission. We found a significant correlation between the number of patient days/year (which is a proxy for hospital size), and the rate of both total (figure 3A) and resistant (figure 3C) S. aureus acquisitions, as well as the percentage of methicillin-resistant isolates among all S. aureus positive blood cultures (see figure 3E). This correlation remained significant even when controlling for all the above-mentioned variables (figure 3B, D, F). In the minimum adequate model chosen on the basis of the Akaike information criterion (as implemented in the function step (lm () ) in R), hospital size was the parameter which overall had the most significant impact. Unsurprisingly, the amount of overall antibiotic usage also strongly correlated with both absolute and relative resistance levels. Nevertheless, the model fits (Figure 3) suggest that the impact of hospital size on resistance level is at least of similar magnitude than the impact of antibiotics consumption. If smaller hospitals had lower resistance levels due to extinction events, the same should generally be true for the total incidence of hospital-acquired infections (although to a lesser degree, since the number of total infections is always higher than the one of resistant infections). To test our findings obtained with the Irish dataset, we used surveillance data from the Pennsylvania Health Care Cost Containment Council (PHC4) [20] to investigate in how far infection rates are proportional to hospital sizes. In total, 149 Hospitals reported nosocomial infection rates (i. e. infections that became symptomatic >48h after admission) quarterly from 2005 to 2007. Furthermore, mortality rates, mean length of stay and infection rates by disease type (e. g. bloodstream infection or pneumonia) were reported once per year. The median infection rate throughout these 12 quarters strongly correlated with hospital size (see figure 4A). This correlation remained very strongly significant even when accounting for potential confounders such as length of stay, mortality or the number of quarters in which electronic surveillance was used (see figure 4B). In the minimum adequate model, the usage of electronic surveillance also had a significant influence on infection rates (see figure 3B). This is presumably because a higher fraction of infections is reported with electronic surveillance. However, even in this analysis the most significant effect (i. e. smallest p-value) on resistance levels was due to hospital size. We observed that infection rates and hospital size were significantly correlated in any given year for the total infection rate as well as for the infection rates of most types of infections (see table 2). Apart from the average level of infection rates, one would also expect that the temporal increase of infection rates would be smaller with frequent extinction events. During the years 2005–2007, there was a slight overall increase in infection rates (see figure 4C). For each hospital, we fitted the total rate of acquiring infection with a linear model with time and presence of electronic surveillance as explanatory variables. If the change over the three years was not significantly correlated to time, we set this change to zero. Also the change in infection rates was significantly correlated to the logarithm of the hospital size (see figure 4D). We mainly consider two simple settings (see Figure 2), in which hospitals of equal size are linked to a community, which is either panmictic or strongly subdivided. In the case of a subdivided community, we allow for migration between the sub-communities (migration rate between 1%/year and 20%/year). These two settings can be considered as models of the population structure in urban and rural environments respectively. In either case we assume that the number of hospital beds per inhabitant (i. e. the fraction between the total hospital size and the community size) is constant. Specifically, we consider a community of 3*105 individuals with 1000 hospital beds. We find that, if drug use is high in the hospital and low in the community, the level of resistance increases with hospital size (Figure 5). This effect is even more pronounced if the community is subdivided as well. However, substructure in the community only seems to have a minor impact compared to the population structure of hospitals (Figure 5). This makes sense intuitively as population sizes in the community are much larger and hence stochastic effects are comparatively weak. The pattern in figure 5 represents the typical situation for a nosocomial pathogen where treatment rates are high in the hospital but low in the community. In the following we generalize this pattern to a broad range of treatment frequencies. We first describe the case of a panmictic community, which as Figure 5 suggests represents the conservative scenario concerning the effect of population size, and then consider the additional effects conferred by substructures in the community.
It has been argued that concentrating highly specialized services in large hospitals will both improve patient outcomes as well as reduce health care costs [24], [25], [26], [27]. Yet, larger hospitals may have their own disadvantages by facilitating the spread of infectious diseases and resistance genes [10]. We used two surveillance datasets from the US [20] and Ireland [19] to illustrate how the incidence of infections and resistance depends on hospital size. In line with the results of Bhavnani et al. [10], we found that the rates of acquiring both resistant infections and infections in general strongly and significantly correlated with patient-days/year (hospital size). Furthermore, the increase of total infections over time and the ratio of resistant and drug-sensitive infections also significantly correlated with hospital size. Of course, different patient populations in small and large hospitals could also explain such a pattern. Antibiotic consumption and morbidity (and thereby susceptibility to infections) might be much higher in larger hospitals caring for more severely ill patients. However, this would not explain the steeper rise of infections in large hospitals. Additionally, the correlations between hospital size and incidence levels persisted even when taking the following potentially confounding factors into account: mortality, length of stay, antibiotic usage, usage of injectable or hospital-specific antibiotics, sampling frequency, or the usage of disinfectants. Both datasets analyzed have both problems and advantages. In the dataset from Ireland, for example, only S. aureus was considered, all positive blood cultures were recorded (even those which might have been acquired outside the given hospital) and there are no data on mortality. Conversely, the dataset from the US does not include the drug-sensitivity of the infections, drug usage or disinfection. However, these drawbacks do not overlap, and our results are coherent for both datasets. This makes it somewhat less likely that such confounders are the only reason for the correlation of resistance with hospital size. Furthermore, hospital size had a more significant impact than antibiotic usage on both relative prevalence and absolute incidence rates of methicillin-resistant S. aureus. Reducing antibiotic consumption is a common recommendation for curbing the spread of resistance. The analyzed data suggest that a reduction in hospital size might therefore be a similarly successful intervention. We used a stochastic epidemiological model in order to describe both the spread of drug resistance in a hospital setting as well as the interaction between one or several hospitals and the community. Like all theoretical models, this is a very simplified description of complex interactions. For example, the epidemiological fitness costs of resistance are inherently difficult to determine. We used estimates of the epidemiological costs of resistance in tuberculosis [28] and of glycopeptide resistance in farm animals [29] as well as in vitro data [30] as guidelines, because no estimates on the transmissibility of resistant nosocomial infections exist, to our knowledge. Additionally, the model neglects co-infection, whereas in reality patients may harbor both resistant and sensitive strains, such that sufficiently long after discontinuation of antibiotic therapy the sensitive strains could dominate again. This effect can, however, be captured in our model by increasing the back-mutation rate (i. e. the rate with which a patient that is predominantly infected with the resistant strain becomes dominated by the sensitive strain in the absence of therapy). We find that increasing this rate substantially increases the magnitude of the effect of hospital size (results not shown). Intuitively, this is the case because with a higher reversion rate the level of resistance in the community (the “sink”) declines. Therefore it becomes less likely that the resistant strain is reintroduced from the community once it becomes extinct in a hospital. Thus, the assumption made here leads to rather conservative estimates for the effect of hospital size. In summary, the results from our theoretical model should be regarded as qualitative descriptions, not as quantitative predictions on how much resistance could be reduced if hospitals were smaller. Theoretically, the impact of hospital size on the evolution of antibiotic resistance can be explained by the meta-population dynamics characterized by local extinctions and re-colonizations. According to this interpretation, the beneficial effect of small hospital size is the result of a simple evolutionary mechanism: Selection typically acts in different directions in the hospital and in the community: resistance is selected for in the hospital and selected against in the community. Thus any mechanism, which weakens the effectiveness of selection in the hospital relative to that in community, will lead to an increase in the level of resistance. According to Fisher' s fundamental theorem, the effectiveness of selection acting on a trait (here trait = resistance) is proportional to the variance in that trait. Small hospital size decreases the variance in resistance through frequent extinctions of the resistance-conferring allele. In other words, if the resistance-conferring allele is extinct in a given hospital, selection (which would potentially favor this allele) is completely ineffective because the population consists only of the sensitive strain and selection can only favor one strain over the other if they coexist. In addition to the hospital-community setting described here, many more instances of population structure might cause similar meta-population effects that hinder the epidemic spread of antibiotic resistance mutations. Indeed, we made the same observations when considering inter-ward transfer. Further examples include: household structure, caring facilities, schools etc. Thus the effects described here are likely to extend beyond the hospital-community setting. In this context it should be noted that for a real hospital of a given size, the magnitude of stochastic effects will typically be considerably stronger than what might be expected from a simple population-biological model like the one presented here. The reason for this is that several processes such as population structure, fluctuating population sizes, and variability in transmission rates (e. g. super-spreaders) strongly enhance stochasticity. Thus, a real hospital with a given number of patients will typically behave like an idealized hospital, which is much smaller. (This type of problem is well known in population genetics [31], and as a consequence populations are often characterized by a so-called effective population size rather than by their census population size). Nevertheless, all other things being equal, a small hospital will be subject to stronger stochastic effects than a large hospital and will hence suffer less from antibiotic resistance. In contrast to the typical beneficial effect of small hospital size, we also found that under some circumstances small hospitals may increase the prevalence of antibiotic resistance. This effect occurs if antibiotic usage in the community reaches similar magnitude than in the hospital. Such a setting might be uncommon for nosocomial infections, which are mainly treated in the hospital. However, it might occur for other infections, which occur frequently in both the hospital and the community (e. g. E. Coli causes many opportunistic infections in the hospital and urinary tract infections in the community). Overall, the results of this study suggest the general pattern that strong population subdivision in those compartments, where antibiotic usage is high (typically: the hospital), can substantially reduce the spread of antibiotic resistance. Furthermore, we find that this beneficial effect of small hospital-size is substantially reduced if a substantial fraction of infection events is not acquired directly from other patients (which exhibit a fast turnover in the hospital), but, indirectly, via a slowly decaying environmental compartment. This latter point indicates that reducing such environmental transmission might be especially important in small hospitals as this might not only reduce the force of infection but also the burden of antibiotic resistance by promoting the stochastic extinction of resistant pathogen strains.
The increasing spread of bacteria, which are resistant to antibiotics, is a serious threat to clinical care. Currently, several countries aim at concentrating highly specialized services in large hospitals in order to improve patient outcomes. However, empirical studies have shown that resistance levels correlate with hospital size. To illustrate this correlation, we analyze two published datasets from the US and Ireland and controlled for antimicrobial usage, disinfection and length of stay. The proportion of patients acquiring both sensitive and resistant infections in hospitals strongly correlated with hospital size. Moreover, we observe the same pattern for both the percentage of resistant infections and the temporal increase of hospital-acquired infections. To investigate to what extent hospital size can directly affect the prevalence of antibiotic resistance, we use mathematical models describing the epidemic spread of resistance in hospitals and the community. We find that small hospitals typically lead to considerably lower resistance levels than large hospitals. However, this beneficial effect of small hospital size may be reduced if bacteria are transmitted indirectly via the environment. Therefore, reducing environmental transmission might be particularly important in small hospitals. Overall, our findings suggest that the short-term benefits of larger hospitals may come at the price of increasing resistance in the long term.
infectious diseases/epidemiology and control of infectious diseases
infectious diseases/antimicrobials and drug resistance
On Being the Right Size: The Impact of Population Size and Stochastic Effects on the Evolution of Drug Resistance in Hospitals and the Community