INTRODUCTION
Urinary tract infections (UTIs) are among the mostly commonly encountered infections in community and hospital settings (
1,
2). By far the most frequently isolated pathogen in these infections is the Gram-negative, facultatively anaerobic bacteria
Escherichia coli (
2,
3). UTIs are generally treatable with a short course of antibiotics, although the prevalence of uropathogens resistant to one or more antibiotics, such as extended spectrum β-lactamase-producing (ESBL)
E. coli, is rising (
1,
4–7). Empirical therapy for UTI is often prescribed without a urine culture or before the results of susceptibility testing are available, with 20% being cited as a threshold past which empirical treatment is compromised (
8). Significant variation exists between regions (
4,
9–12). Consequently, Infectious Diseases Society of America (IDSA) guidelines recommend that treatment selection be based on local or regional resistance rates (
8); there is also a need for early identification of patients at risk for treatment-resistant infections.
Much of the previous research on risk factors for antimicrobial resistance in urinary
E. coli infections relies on crosstabular methods that present resistance proportions for various subgroups (e.g., men versus women, children versus adults), which makes it difficult to appreciate the relative contributions of each variable to observed variation in the prevalence of resistance between patients, places, times, and settings of acquisition (community-acquired or nosocomial). Additionally, resistance in an isolate to each antimicrobial is typically treated as if it is independent of others, even though mechanisms of antimicrobial resistance are rarely confined to a single drug (
13,
14). For example, the genes in ESBL-producing
E. coli conferring penicillin and cephalosporin resistance often cooccur on plasmids with genes granting resistance to other classes of antimicrobials, such as aminoglycosides, trimethoprim-sulfamethoxazole, and quinolones (
15). Appropriate statistical methods can leverage these correlations to more precisely estimate risk factors for each type of resistance. Advances in data interoperability have made it easier to perform studies integrating a greater breadth of data, rather than relying on relatively small samples collected over a short period of time, often from a single hospital.
In the province of Quebec, Canada, the great majority of microbiological testing, including antimicrobial susceptibility testing, is done in hospital rather than in private laboratories. Adult urinary E. coli isolates are regularly tested for resistance to six antimicrobials representing five classes, ampicillin (penicillin), gentamicin (aminoglycoside), ciprofloxacin (fluoroquinolone), nitrofurantoin (nitrofuran), trimethoprim-sulfamethoxazole (TMP-SMX; combination dihydrofolate reductase inhibitor/sulfonamide), and tobramycin (aminoglycoside). Several major health care networks have adopted a common infection control software (Nosokos; Nosotech, Rimouski, Canada), facilitating the aggregation of many years of susceptibility testing results from across the province using a common data dictionary. Using hierarchical logistic regression methods, it is possible to overcome the limitations of previous research by simultaneously analyzing many potential risk factors and allowing risk factors to be correlated across multiple types of resistance to different antimicrobials. The latter also permits the estimation of an overall association for a risk factor across all types of resistance. In this study, we developed models to investigate the temporal, geographic, and patient-level predictors of resistance to six antimicrobials in a large sample of community-acquired and nosocomial urinary E. coli isolates from three communities in the province of Quebec, Canada.
DISCUSSION
Empirical treatment of urinary tract infections should be based on susceptibility profiles specific to the locale and patient population. Susceptibility testing data capable of informing treatment decisions is present in hospital information systems, but has to date been underutilized for this purpose, due to a lack of appropriate statistical methods and interoperability across systems. In this study, we jointly modelled the temporal, geographic, and patient-level associations with resistance to six antimicrobials in urinary E. coli isolates from inpatients and outpatients from three cities across the province of Quebec using a shared data structure.
Our results showed commonalities in risk factors across drugs and between community-acquired and nosocomial isolates. The marked divergence in the prevalence of resistance between the three cities in our study underscores the importance making local resistance data available to prescribers, rather than relying on provincial or national estimates. Men and recently hospitalized individuals were generally at a higher risk for antimicrobial resistance, whereas age predicted higher resistance to only some antimicrobials. Our hierarchical approach accounted for correlations between types of resistance and enabled us to more precisely estimate risk factors for individual types of resistance in specific patient populations. As an illustration, consider a frontline clinician in Quebec City who uses a local antibiogram to determine that the average rate of susceptibility to TMP-SMX in urinary
E. coli is 82%. The clinician also knows that sex and hospitalization in the past 30 days influence susceptibility rates, but does not know the precise magnitude of these influences. With the type of model we fit in this study, the clinician would know to expect 84.2% (95% CI, 83.1% to 85.2%) susceptibility if the patient were a 70-year-old woman but only 75.2% (95% CI, 72.9% to 77.6%) susceptibility if the patient were a 30-year-old, recently hospitalized woman. These two predicted susceptibilities are on opposite sides of the 80% susceptibility threshold cited for empirical treatment with TMP-SMX (
8).
This study’s results support some of the findings from previous research. Both male sex and recent hospitalization have been previously identified as positive predictors of resistance in urinary
E. coli (
9,
16–20), although the association with sex is not uniform across drugs (e.g., see references
9 and
19). While male sex typically accounted for only a few percentage points of increased risk in our study, when combined with stronger risk factors, such as recent hospitalization, community, and potentially age, dramatically different resistance profiles can result, particularly for the mainline treatment TMP-SMX (e.g., see
Fig. 1 and
3). Differences in etiology between men and women may be important, as community-acquired UTIs are common in otherwise healthy women (
21) but in men are often associated with anatomical abnormalities, such as an enlarged prostate (
22). Nonetheless, the tendency for men to have higher resistance to some antimicrobials persisted in the nosocomial setting, despite the similarity in the primary route of acquisition (urinary catheterization). Rising rates of resistance over time are consistent with trends in urinary
E. coli (
23) and pathogens in general (
24).
The evidence regarding the relationship between age and resistance in urinary
E. coli suggests that the association is highly dependent on the drug in question (
9,
19,
20). Age has been consistently linked to an elevated risk of ciprofloxacin resistance in adults (
4,
19,
20,
25,
26), as observed in our fitted model. This may be linked to cumulative changes in the gut microbiome resulting from the use of broad-spectrum antimicrobials like ciprofloxacin (
27), altering an individual’s susceptibility to invasion by antimicrobial-resistant organisms (
28). Seasonal trends in community-acquired infections likely relate to patterns of antibiotic prescribing (
29,
30) and may also be influenced by other factors, like international travel (
31). Resistance rates were remarkably higher in densely urban Montreal compared to the other two sites, with the more remote Rimouski generally showing the lowest rates. Our results are consistent with the hypothesized positive association between population density and resistance (
32), although we cannot draw general conclusions from only three sites. It should be noted, however, that the McGill University Health Centre serves as a major reference center for urology in Montreal and may receive complex cases from other hospitals, which could be overrepresenting the prevalence of resistance in the city.
Surprisingly, our model does not generally predict large differences in the risk of antimicrobial resistance among isolates from elderly (65+ years of age) individuals based on setting of acquisition (community or nosocomial; see
Fig. 3). A metaanalysis of 54 observational studies by Fasugba et al. (
23) concluded that ciprofloxacin resistance was higher in hospital-acquired (38%; 95% confidence interval, 36% to 41%) compared with community-acquired (27%; 95% confidence interval, 24% to 31%) infections, although they did not directly control for confounders such as sex, previous hospitalization, or age among adult patients. Fleming et al. (
33) reported that in 156 urinary
E. coli taken from a Georgia hospital, prevalence of resistance to TMP-SMX (in addition to that to several other antimicrobials) was higher in the hospital-acquired (34.6%) than the community-acquired (25.2%) isolates, but did not adjust for demographic factors. Other studies (
4,
34,
35) also report differences in rates of resistance between these two subpopulations (interestingly, Lob et al. [
4] report these differences in the United States but not Canada), but it is also necessary to note the wide variation in estimated rates of resistance between studies. Geographic differences, as well as the demographic associations detected in our large database of isolates, help to explain this variation. It is possible that, had we compared nosocomial isolates from adults under 65 (which are much less common) to community-acquired isolates in the same age range, the results might have been different.
Our study had several limitations. The lack of supplementary clinical data (e.g., symptoms and clinical outcomes) meant that we were unable to distinguish asymptomatic bacteriuria from symptomatic UTI, or the type and severity of infections; we also lacked data on many factors potentially relevant to treatment, such as catheterization, prior antibiotic use, comorbid illnesses, and drug interactions. Since our sampling frame included only patients for which clinical specimens had been taken, and specimens may not be systematically taken for uncomplicated UTIs, our data set may not be reflective of all treated UTIs. The purpose of a standard antibiogram is to provide a baseline expected level of resistance, and our model improves these baseline predictions by making them more precisely targeted to particular patient groups. These models are intended to supplement, not replace, good clinical judgment, especially regarding the usually unnecessary treatment of asymptomatic bacteriuria and the judicious use of broad-spectrum antimicrobials.
Changes in clinician behavior over time or between locations could explain some of the variability in rates of resistance. Since we cannot confirm the length of catheter use in hospitalized patients or identify when a patient was hospitalized in a hospital network outside where their sample was tested (more likely in Montreal and Quebec City, which have many hospitals), there is a possibility of misclassifying nosocomial isolates as community-acquired, as well as of incorrectly identifying previous hospitalization status. This could obscure the differences we observed between nosocomial and community-acquired isolates and blunt the increased risk of resistance seen in recently hospitalized patients.
The value of electronic health records for surveillance, quality improvement, and epidemiological research is increasing as databases become more interlinked. This study demonstrates the utility of standardized antimicrobial resistance data from multiple institutions to produce locally relevant profiles of antimicrobial resistance. Our modeling approach allowed us to make inferences about temporal, geographic, and demographic variation in the probability of resistance to 6 types of antimicrobials used in the treatment of urinary tract infections. Clinically relevant differences in resistance between communities and patient populations in the province of Quebec could inform empirical treatment decisions. In the future, a model-based approach for antimicrobial resistance informed by local, provincial, and national trends could be incorporated into decision support systems for clinicians.