A spatial approach for the epidemiology of antibiotic use and resistance in community-based studies: the emergence of urban clusters of Escherichia coli quinolone resistance in Sao Paulo, Brasil
© Kiffer et al; licensee BioMed Central Ltd. 2011
Received: 22 November 2010
Accepted: 28 February 2011
Published: 28 February 2011
Population antimicrobial use may influence resistance emergence. Resistance is an ecological phenomenon due to potential transmissibility. We investigated spatial and temporal patterns of ciprofloxacin (CIP) population consumption related to E. coli resistance emergence and dissemination in a major Brazilian city. A total of 4,372 urinary tract infection E. coli cases, with 723 CIP resistant, were identified in 2002 from two outpatient centres. Cases were address geocoded in a digital map. Raw CIP consumption data was transformed into usage density in DDDs by CIP selling points influence zones determination. A stochastic model coupled with a Geographical Information System was applied for relating resistance and usage density and for detecting city areas of high/low resistance risk.
E. coli CIP resistant cluster emergence was detected and significantly related to usage density at a level of 5 to 9 CIP DDDs. There were clustered hot-spots and a significant global spatial variation in the residual resistance risk after allowing for usage density.
There were clustered hot-spots and a significant global spatial variation in the residual resistance risk after allowing for usage density. The usage density of 5-9 CIP DDDs per 1,000 inhabitants within the same influence zone was the resistance triggering level. This level led to E. coli resistance clustering, proving that individual resistance emergence and dissemination was affected by antimicrobial population consumption.
Mortality due to infections represents approximately 85% of all deaths worldwide , and community acquired ones are highly prevalent [1, 2]. Furthermore, most of the community infections present a risk of resistance acquisition to a first line antimicrobial drug . The increasing burden of antimicrobial resistance associated to the relative drought on the antimicrobial developmental pipeline has urged for strategies to reduce antibiotic consumption and to maximize treatment approaches [3, 4]. These strategies, as stated by Lipsitch and Samore , "rest on the well-supported idea that the use of antimicrobial agents is a powerful selective force that promotes the emergence of resistant strains". Although infections are ecological phenomena [3, 5–7] and require a new analytical paradigm, researches have generally taken a case-by-case approach to understanding the emergence of new infections or of its consequences . Nevertheless, for any infectious disease, the infection or colonization status of a subject affects the acquisition risk of this very same condition by others and this is due to the transmissibility risk , or the so-called ecological phenomenon. Based on this, antibiotic usage as a driver of resistance shall be equally affected. Should a person or a population be submitted to an antibiotic usage, other people around who have not received antibiotics may have increased colonization or infection risk of acquiring a resistant organism . Thus, it is reasonable to assume that population antimicrobial usage may be a major driver in determining the emergence of resistance in community. However, a potential population pressure affecting antimicrobial resistance emergence has not been quantitatively demonstrated so far.
The EUREQA Project (Epidemiologia do Uso e da Resistência de Antibióticos e Quimioterápicos na população)  is a multi-phase spatial and temporal project aimed at correlating population risk factors, mainly antimicrobial usage density, with community bacterial resistance. The present study aimed at establishing quantitative spatial correlations among antimicrobial usage mean density as measured for the entire population and resistance to ciprofloxacin in community urinary tract infections (UTI) caused by E. coli isolated from adult women in a major urban environment, São Paulo city, which occupies an area of only 0.02% of the Brazilian Territory with a population of approximately 10.5 million inhabitants.
The present study has been submitted and approved by the Ethics Committee of the Federal University of Sao Paulo (UNIFESP, process CEP 0123/05) and it was based on data observation without patient identification.
The EUREQA database stores two essential informations: i) UTI E. coli events by individual address; ii) antimicrobial selling points. E. coli events had the following case definition: routine urine culture results positive for E. coli (single isolate per culture, with ≥ 100.000 UFC/ml), collected on two large healthcare outpatient facilities, covering a public and a private sector unit. Since E. coli occurs most commonly in female patients, age and sex-specific restrictions were adopted and only data from female patients with age higher than 16 years-old were included in this first analysis. In order to detect resistance clusters, cases were categorized in susceptible (antibiogram = sensitive) or resistant (antibiogram = resistant or intermediate) to ciprofloxacin. Ciprofloxacin was the antimicrobial chosen based on its rapid resistance development by E. coli and other Gram-negative bacteria [9–12]. The antimicrobials selling points were obtained from a pharmaceutical market auditing company, IMS Health Brazil, covering 98% of the studied area. They were categorized as drugstore, delivery pharmacy, pharmacy chain, hospital, or major distributor. Those categorized as hospitals and major distributors were excluded from the final database, which focused on the community market exclusively. Each selling point contained monthly Defined Daily Doses (DDD) information of the all antimicrobials included in the EUREQA Project (all with ATC code ). The DDD by WHO definition  was the basic unit chosen for building a population antimicrobial consumption measurement. For the present study, DDDs for ciprofloxacin and other quinolones (nalidixic acid, norfloxacin, levolfoxacin, gatifloxacin, moxifloxacin) were kept in the database. For each selling point, a mean cumulative DDD of a selected antimicrobial (i.e. ciprofloxacin) for the previous three months (DDDA) was calculated. Exceptions included the months of the January and February, which used replicated values, according to the following: January, DDDA Jan = 3*(DDD Jan )/3, and February, DDDA Fev =(1.5*(DDD Jan ) +1.5*(DDD Fev )/3. The selection of three months window accumulated data was based on time-series studies in hospitals where resistance associated to increased antimicrobial consumption emerged usually no later than this timeframe in gram-negatives .
Selling points influence zone determination
Influence zones were used as the best explanation for consumption pattern, since pharmacy consumption, either by direct buy or delivery, originates from a zone closer to the subject`s residence. This way, a DDDA value initially belonging to a point, was distributed over an area, basically considering that all individuals over that area were submitted to that consumption influence. This procedure was irrespective of individuals actually using or not the antimicrobials measured. It only considered that individuals were under the same consumption influence and, thus, under the same risk for resistance acquisition due to transmissibility or indirect effect .
Antimicrobial consumption density determination
where DDDA k,t is the average of the DDD cumulative for three months for k-th influence zone and for t-th month (t = 1,2,..., 12), Pop k / 1,000 is the estimated population by a thousand inhabitants restricted to k-th influence zone and 30 is the number of days in the month.
where n is the number of the intersections areas that occurs between the k-th influence zone and the census tract map, i is the i-th census tract, CTPop i is a population of the i-th census tract, CTArea i is the area of the i-th census tract and INT_Area k,i is the intersection area between of the k-th influence zone and the i-th census tract.
Finally, each resistant or susceptible case contained within an influence zone of the respective usage density, received a value of the respective D_DDDA k,t according to date of occurrence (closest possible data of the D_DDDA k,t measurement).
Generalized additive model: a semi-parametric binary regression approach
To detect resistance clusters in the study area, a spatial Generalized Additive Model (GAM) [18, 19] was adopted relating E. coli resistance and relevant covariates, spatial and temporal effects. Both the global test for the spatial variation of resistance risk and the identification of areas of high/low risk were performed using the Monte Carlo simulation method proposed by Kelsall and Diggle .
The modeling structure used is based upon a spatial point process  under which a risk measure can be defined and estimated continuously over the study region. This procedure firstly defines a risk measure to detect ciprofloxacin resistant clusters in the study area A. For this, the susceptible cases are treated as a sample of the population at risk: the ciprofloxacin non resistant (susceptible) residents of São Paulo city.
Consider the locations of the n1 resistant cases and the locations of the (n - n1)susceptible cases in the study region A as observations of two Poisson processes I and II, with intensities λ1 (s) and λ2 (s), respectively. A log risk measure at location s can be defined as ρ (s) = log(λ1 (s)/λ2 (s)) and the aim is to investigate the spatial variation of ρ(s) on A.
A semi-parametric approach to estimating ρ(s) is used, called the Generalized Additive Model (GAM) . The GAM method allows estimating the spatial risk whilst controlling for potential individual or local environmental factors (such as exposure to D_DDDA, for instance).
Where c = log(q1/q2), so c is simply an additive constant, therefore it does not modify the overall characteristics of the spatial distribution of the risk over the region. Hence, one can obtain estimates of ρ(s), apart from an additive constant, with a model fitted to a binary outcome. For easy of interpretation the logarithmic scale is used with a base 2, as in this scale a unit increase in the log-risk surface from one location to another implies the doubling of risk.
where x is the vector of individual/local prognostic factors, β are the effects and g(s) is a smooth (unknown) function of s.
If the risk is constant then g(s) = 0 and model (Equation I) becomes the usual logistic regression model . Hence, model (Equation I) is the usual logistic regression model extended by an additive component g(s), so the interpretation of factor effects and the spatial effect are made accordingly.
where β0: regression model intercept; β1, β2, and β3: covariates effects; D_DDDA: antimicrobial usage mean density; I1 = 1 if T1< D_DDDA <T2 and I1 = 0 otherwise; I2 = 1 if D_DDDA > T2; and I2 = 0 otherwise. E. coli resistance change points (T1 and T2)_were identified after performing an initial analysis on the data by fitting an ordinary logistic regression model without the spatial term but allowing for an additive (smooth) effect of D_DDDA upon the risk of E. coli ciprofloxacin resistance.
Additionally a sensitivity analysis was performed and although correlation between spatial effects and regression coefficients for spatially varying covariates may be a general concern, it did not prove relevant for the present data analysis. A generalized additive model allowing for a smooth function of D_DDDA in the linear predictor replacing the change point model together with a smooth spatial effect was also applied (solid line shown in Figure 4) and results maintained their significance at 5% level. Furthermore, the presented model has a comparable goodness of fit with the GAM model, thus being chosen for the present study due to its clearer interpretation.
Discussion and Conclusions
Different epidemiological determinants may favor the emergence or establishment of specific resistances in given environments [5, 23–25], in particular after shown that resistance lies within the human commensal bacterial species, which were once considered relatively harmless residents of the human microbiome . It also seems reasonable to speculate that a human-driven increase in antibiotic concentrations of a given ecosystems, such as a city, may influence both antibiotic resistance and the microbial population dynamics . As previously noted , a population could experience indirect effects of antimicrobial use, such as an enhanced risk for resistance acquisition, because of antimicrobial usage by others in the population. Nevertheless, it is still unknown how these different environmental and individual determinants are distributed over space and time and their possible influences on a resistance emergence or clonal spread. It has been shown by time-series analysis that antimicrobial usage in a restricted and contained environment, such as a hospital, is temporally linked to the emergence of bacterial resistance [9, 10]. However, no correlation of antimicrobial population usage and bacterial resistance in the community had been established before the present study.
The present model used selling points as the basis for influence zone determination and this concept was supported by the premise that pharmacy consumption patterns, either by direct buy or delivery, originates from a zone closer to the subject`s residence. Nevertheless, authors do acknowledge that potential biases may exist in this concept, since consumption patterns may also be related to a subject`s habit of using pharmacies elsewhere in the city and that an overestimation of CIP usage may occur in commercial areas of a city. However, based on this model, it was found that the level of D_DDDA which triggered the occurrence of resistance in this study was established between 5 and 9 DDDs of CIP/1,000 inhabitants-day and that CIP D_DDDA was found to be related to the clusters. It means that in a population of 1,000 inhabitants in a given day, if a mean of 5 to 9 DDDs of ciprofloxacin had been consumed over the past three months, E. coli resistance would occur in subjects of the same geographical area (not necessarily the same individuals taking the medication). However, the model could not rule out other population covariates which might be associated to the elevated risk areas of clustered E. coli CIP resistance. Nevertheless, the model allows for the inclusion of many different covariates, or risk factors, to be investigated, as long as they can be geo-referenced, once other population covariates could be associated to the CIP resistant E. coli clusters.
We do envisage a combined approach to explore associated risk factors, since molecular epidemiology approaches have given insights into the area. Recently, the occurrence of an ecological phenomenon has been demonstrated  by a single clonal group accounting for nearly half of community-acquired urinary tract infections in women caused by E. coli with resistance to trimethoprim-sulfamethoxazole in three geographically diverse communities. In a different study, a geographical information system and a Bernoulli regression model were applied to detect clusters of higher risk of acquisition of S. aureus soft tissue abscesses . Also, there have been strategies used with spatial scan statistics to identifying clusters of samples and to detecting areas with significantly high or low sampling rates of a national antimicrobial resistance monitoring program . Jones et al.  applied a spatial analysis methodology and found more than half of emerging infectious diseases events to be caused by bacteria or rickettsia, with a large number of drug-resistance and significant correlation with socio-economic, environmental and ecological factors. There have been other significant efforts to understanding resistance dynamics and transmission with different approaches. It should be noticed that at least two other strategies [30, 31] were based on bacterial surveillance data: McCormick et al  have shown patterns of geographic variation explained by the differences between the proportions of resistance in specific S. pneumoniae serotypes. Stelling et al  used an electronic laboratory data system (WHONET) and a space-time permutation scan statistic for generating a semi-automated detection of disease outbreak with Shigella spp . Both latter approaches seem to be complimentary to the present one. Furthermore, all these approaches, including the present one, adopt a population-level perspective to elucidate the infection emergence problem and risk factors associated, in particular the resistance issue.
On the one hand individual-level studies do not provide the opportunity to observe the ecological phenomena, or indirect effects of treatment on resistance [3, 32]. On the other, although avoiding the statistical non-independence issue between samples, population-level studies have not been able to provide a replicable framework on which one may test the hypothesis of the indirect effects of treatment on resistance.
The present model has been applied on a longitudinal observational study, and as such, conclusions are limited to the environment used, even considering the elevated sample size and the goodness of fit of model. As a cluster-randomized study of antimicrobial resistance, we applied the ratio of resistant to susceptible as an endpoint. The spatial clustering of antimicrobial resistance in this population points to an ecological effect. Nevertheless, the adopted framework allows for the inclusion of any number of temporal and spatial variables. In particular, for the study of resistance emergence, it is possible to fit in the model virtually any antimicrobial-microorganism combination, in pairs or in groups, if the objective is to observe class effects. We anticipate our spatial-temporal modelling strategy embedded into a geographic information system  analytical framework to be the start-up point for validating, replicating, and understanding antimicrobial resistance in community and for public health interventions.
We thank the EUREQA group - Caio MF Mendes, Jussara O Ortiz, Virginia R Correia, Corina C Freitas, Soraya Andrade, Paula Koga, Amilton Mouro, Miriam Santos - for the participation and the support in developing the study. We thank also INPE, UNIFESP and the Fleury Institute for Research & Development for institutional support and cooperation.
- World Health Organization: Department of Communicable Disease Surveillance and Response. WHO Global Strategy for Containment of Antimicrobial Resistance (WHO/CDS/CSR/DRS/2001.2). 2001, WHO, Switzerland, 1-96.Google Scholar
- Black RE, Morris SS, Bryce J: Where and why are 10 million children dying every year?. Lancet. 2003, 361: 2226-2234. 10.1016/S0140-6736(03)13779-8.View ArticlePubMedGoogle Scholar
- Lipsitch M, Samore MH: Antimicrobial use and antimicrobial resistance: A population perspective. Emerging Infectious Diseases. 2002, 8: 347-354. 10.3201/eid0804.010312.PubMed CentralView ArticlePubMedGoogle Scholar
- Infectious Diseases Society of America: Bad bugs, no drugs, as antibiotic discovery stagnates, a public health crisis brews.http://www.idsociety.org/BBNDWhitePaper04.htmhttp://www.idsociety.org/BBNDWhitePaper04.htm
- Harbarth S, Samore MH: Antimicrobial resistance determinants and future control. Emerg Infect Dis. 2005, 11: 794-801.PubMed CentralView ArticlePubMedGoogle Scholar
- Woolhouse ME: Epidemiology: emerging diseases go global. Nature. 2008, 451: 898-899. 10.1038/451898a.View ArticlePubMedGoogle Scholar
- Jones KE, Patel NG, Levy MA, Storeygard A, Balk D, Gittleman JL, Daszak P: Global trends in emerging infectious diseases. Nature. 2008, 451: 990-U994. 10.1038/nature06536.View ArticlePubMedGoogle Scholar
- EUREQA: Epidemiologia do Uso e da Resistência de Antibióticos e Quimioterápicos na população. EUREQA Project, FAPESP Research Grant 2006/53922-9.http://www.dpi.inpe.br/eureqahttp://www.dpi.inpe.br/eureqa
- Aldeyab MA, Monnet DL, Lopez-Lozano JM, Hughes CM, Scott MG, Kearney MP, Magee FA, McElnay JC: Modelling the impact of antibiotic use and infection control practices on the incidence of hospital-acquired methicillin-resistant Staphylococcus aureus: a time-series analysis. Journal of Antimicrobial Chemotherapy. 2008, 62: 593-600. 10.1093/jac/dkn198.View ArticlePubMedGoogle Scholar
- Lopez-Lozano JM, Monnet DL, Yague A, Burgos A, Gonzalo N, Campillos P, Saez M: Modelling and forecasting antimicrobial resistance and its dynamic relationship to antimicrobial use: a time series analysis. International Journal of Antimicrobial Agents. 2000, 14: 21-31. 10.1016/S0924-8579(99)00135-1.View ArticlePubMedGoogle Scholar
- Cizman M: The use and resistance to antibiotics in the community. International Journal of Antimicrobial Agents. 2003, 21: 297-307. 10.1016/S0924-8579(02)00394-1.View ArticlePubMedGoogle Scholar
- Livermore DM: Minimising antibiotic resistance. Lancet Infectious Diseases. 2005, 5: 450-459. 10.1016/S1473-3099(05)70166-3.View ArticlePubMedGoogle Scholar
- WHO: Collaborating Centre for Drug Statistics Methodology Norwegian Institute of Public Health.http://www.whocc.no/atcddd/http://www.whocc.no/atcddd/
- Okabe A: Spatial tessellations: concepts and applications of Voronoi diagrams. 2000, Chichester; New York: Wiley, 2View ArticleGoogle Scholar
- Aurenhammer F: Voronoi Diagrams - a Survey of a Fundamental Geometric Data Structure. Computing Surveys. 1991, 23: 345-405. 10.1145/116873.116880.View ArticleGoogle Scholar
- Bowyer A: Computing Dirichlet Tessellations. Computer Journal. 1981, 24: 162-166. 10.1093/comjnl/24.2.162.View ArticleGoogle Scholar
- Hastie T, Tibshirani R: Generalized Additive-Models - Some Applications. Journal of the American Statistical Association. 1987, 82: 371-386. 10.2307/2289439.View ArticleGoogle Scholar
- Hastie T, Tibshirani R: Generalized additive models. 1990, London; New York: Chapman and Hall, 1Google Scholar
- Kelsall JE, Diggle PJ: Spatial variation in risk of disease: a nonparametric binary regression approach. Journal of the Royal Statistical Society Series C-Applied Statistics. 1998, 47: 559-573. 10.1111/1467-9876.00128.View ArticleGoogle Scholar
- Diggle PJ: Statistical Analysis of Spatial Point Patterns. 1983, Academic PressGoogle Scholar
- Hosmer DW, Lemeshow S: Applied logistic regression. 1989, New York: WileyGoogle Scholar
- Harbarth S, Albrich W, Goldmann DA, Huebner J: Control of multiply resistant cocci: do international comparisons help?. Lancet Infect Dis. 2001, 1: 251-261. 10.1016/S1473-3099(01)00120-7.View ArticlePubMedGoogle Scholar
- Manges AR, Johnson JR, Foxman B, O'Bryan TT, Fullerton KE, Riley LW: Widespread distribution of urinary tract infections caused by a multidrug-resistant Escherichia coli clonal group. N Engl J Med. 2001, 345: 1007-1013. 10.1056/NEJMoa011265.View ArticlePubMedGoogle Scholar
- Morten OAS, Dantas G, Church GM: Functional Characterization of the Antibiotic Resistance Reservoir in the Human Microflora. Science. 2009, 325: 1128-10.1126/science.1176950.View ArticleGoogle Scholar
- Hanage WP, Fraser C, Tang J, Connor TR, Corander J: Hyper-recombination, diversity, and antibiotic resistance in pneumococcus. Science. 2009, 324: 1454-1457. 10.1126/science.1171908.View ArticlePubMedGoogle Scholar
- Martinez JL: Antibiotics and antibiotic resistance genes in natural environments. Science. 2008, 321: 365-367. 10.1126/science.1159483.View ArticlePubMedGoogle Scholar
- Tirabassi MV, Wadie G, Moriarty KP, Garb J, Konefal SH, Courtney RA, Sachs BF, Wait R: Geographic information system localization of community-acquired MRSA soft tissue abscesses. Journal of Pediatric Surgery. 2005, 40: 962-966. 10.1016/j.jpedsurg.2005.03.010.View ArticlePubMedGoogle Scholar
- Vieira AR, Houe H, Wegener HC, Wong DMALF, Bodker R, Emborg HD: Spatial Scan Statistics to Assess Sampling Strategy of Antimicrobial Resistance Monitoring Program. Foodborne Pathogens and Disease. 2009, 6: 15-21. 10.1089/fpd.2008.0132.View ArticlePubMedGoogle Scholar
- McCormick AW, Whitney CG, Farley MM, Lynfield R, Harrison LH, Bennett NM, Schaffner W, Reingold A, Hadler J, Cieslak P: Geographic Diversity and Temporal Trends of Antimicrobial Resistance in Streptococcus pneumoniae in the United States. Nature Medicine. 2003, 9: 424-430. 10.1038/nm839.View ArticlePubMedGoogle Scholar
- Stelling J, Yih WK, Galas M, Kulldorff M, Pichel M, Terragno R, Tuduri E, Espetxe S, Binsztein N, O'Brien TF: Automated use of WHONET and SaTScan to detect outbreaks of Shigella spp. using antimicrobial resistance phenotypes. Epidemiology and Infection. 2010, 138: 873-883. 10.1017/S0950268809990884.PubMed CentralView ArticlePubMedGoogle Scholar
- Tacconelli E: Antimicrobial use: risk driver of multidrug resistant microorganisms in healthcare settings. Current Opinion in Infectious Diseases. 2009, 22: 352-358. 10.1097/QCO.0b013e32832d52e0.View ArticlePubMedGoogle Scholar
- Camara G, Vinhas L, Ferreira KR, Queiroz GR, Souza RCM, Monteiro AMV, Carvalho TM, Casanova MA, Freitas UM: TerraLib: An Open Source GIS Library for Large-scale Environmental and Socio-economic Applications. 2008, Berlin: SpringerGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.