OBJECTIVE: To determine whether childhood antianaerobic antibiotic exposure is associated with the development of inflammatory bowel disease (IBD).
METHODS: This retrospective cohort study employed data from 464 UK ambulatory practices participating in The Health Improvement Network. All children with ≥2 years of follow-up from 1994 to 2009 were followed between practice enrollment and IBD development, practice deregistration, 19 years of age, or death; those with previous IBD were excluded. All antibiotic prescriptions were captured. Antianaerobic antibiotic agents were defined as penicillin, amoxicillin, ampicillin, penicillin/β-lactamase inhibitor combinations, tetracyclines, clindamycin, metronidazole, cefoxitin, carbapenems, and oral vancomycin.
RESULTS: A total of 1 072 426 subjects contributed 6.6 million person-years of follow-up; 748 developed IBD. IBD incidence rates among antianaerobic antibiotic unexposed and exposed subjects were 0.83 and 1.52/10 000 person-years, respectively, for an 84% relative risk increase. Exposure throughout childhood was associated with developing IBD, but this relationship decreased with increasing age at exposure. Exposure before 1 year of age had an adjusted hazard ratio of 5.51 (95% confidence interval [CI]: 1.66–18.28) but decreased to 2.62 (95% CI: 1.61–4.25) and 1.57 (95% CI: 1.35–1.84) by 5 and 15 years, respectively. Each antibiotic course increased the IBD hazard by 6% (4%–8%). A dose-response effect existed, with receipt of >2 antibiotic courses more highly associated with IBD development than receipt of 1 to 2 courses, with adjusted hazard ratios of 4.77 (95% CI: 2.13–10.68) versus 3.33 (95% CI: 1.69–6.58).
CONCLUSIONS: Childhood antianaerobic antibiotic exposure is associated with IBD development.
- aHR —
- adjusted hazard ratio
- CI —
- confidence interval
- HR —
- hazard ratio
- IBD —
- inflammatory bowel disease
- IQR —
- interquartile range
- THIN —
- The Health Improvement Network
What’s Known on This Subject:
Inflammatory bowel disease pathogenesis is incompletely understood. Previous pediatric studies suggested associations between antibiotic use and inflammatory bowel disease development but were limited by recall bias, lack of controls, incomplete antibiotic capture, or included exposures between symptom onset and diagnosis.
What This Study Adds:
Our population-based cohort study suggests that certain childhood antibiotic exposures are associated with an increased risk of developing inflammatory bowel disease. Our findings have implications for understanding the condition’s pathogenesis and provide additional stimulus for reducing unnecessary childhood antibiotic use.
Inflammatory bowel disease (IBD) is a chronic intestinal inflammatory condition without clear etiology. Although genome-wide association studies have identified novel genetic defects associated with IBD, these defects account for only one-half of IBD cases.1,2 Childhood IBD incidence and prevalence seem to have doubled over the last decade.3 A leading hypothesis of IBD pathogenesis is that gut bacterial community alterations, with either increases in pathogenic bacteria or decreases in protective bacteria, trigger inflammation.4
Consistent with this hypothesis, multiple studies have noted associations between intestinal infections and IBD onset, although these associations may be secondary to ascertainment bias and increased stool culture testing in those undergoing evaluation for intestinal symptoms.5–7 In addition, patients with IBD have reduced diversity of intestinal microbiota (specifically, anaerobic bacteria) relative to healthy controls, although the causal relationship is unclear.8 Moreover, murine IBD models demonstrate that the IBD phenotype can be induced or ameliorated through flora sharing alone.9,10 Lastly, antibiotics alter the human gut microbiome, with transiently decreased bacterial diversity returning to its preantibiotic state within ∼1 month.11
Previous studies have shown associations between antibiotic use and IBD development but were limited due to recall bias,12–15 lack of controls,16 incomplete antibiotic capture,17,l8 or potential selection bias due to case-control design19; examined predominantly adults20,21 or only certain antibiotic classes22; or included antibiotic exposure between IBD symptom onset and eventual IBD diagnosis.23
We sought to examine the association between childhood antianaerobic antibiotic exposure and subsequent IBD development using a large population-based cohort, hypothesizing that exposure to antibiotics with anaerobic activity would be associated with the development of IBD.
Study Design and Data Source
We performed a retrospective cohort study by using The Health Improvement Network (THIN) database. THIN is a UK collaborative effort collecting de-identified electronic medical record data from 464 participating general practices during the study period. THIN contains data on >9 million (>3 million active) patients, representing ∼5.7% of UK patient visits.24 The Vision general practice computing system (InPractice Systems, London, UK) is an electronic medical record capturing patient demographic characteristics, diagnoses, prescriptions, and some laboratory data, and was adopted by THIN practices from 1994 to 2008.
This study was approved by the institutional review board of the University of Pennsylvania and the THIN Scientific Review Committee.
All children ≤17 years of age enrolled in THIN practices using the Vision electronic medical record between 1994 and 2009 and with ≥2 years’ follow-up were included, because some adults experience symptoms for 2 years before diagnosis of IBD.20 Subjects with an IBD diagnosis occurring before their practice’s computerized medical record adoption were excluded. In addition, research in the UK’s General Practice Research Database, which extracts data from the same medical practice software and overlaps some THIN practices, demonstrated that chronic conditions are overdiagnosed in the first 6 months after registration (as new providers record previous conditions), so we excluded subjects with first IBD diagnosis within 6 months of cohort entry to ensure capture of incident IBD cases only.25 Subjects contributed person-time until first IBD diagnosis, practice deregistration, 19 years of age, or death. The last available data were from November 25, 2009.
Data on all systemic antibiotic prescriptions were captured; antibiotics are unavailable in the United Kingdom without prescription. Antianaerobic antibiotics were defined as penicillin, amoxicillin, ampicillin, penicillin/β-lactamase inhibitor combinations, tetracyclines, clindamycin, metronidazole, cefoxitin, carbapenems, and oral vancomycin. We first classified subjects as ever- or never-exposed. An antibiotic course was defined as continuous antibiotic exposure with <3 days’ interruption. Antibiotic courses were measured in weeks.
THIN does not capture inpatient prescribing. To determine whether confounding by hospital admissions might exist (when subjects might receive unobserved antibiotics and which might represent treatment of early IBD symptoms), hospital admissions were defined as any diagnoses or procedures recorded on unique days with an inpatient location assigned. We did not assess probiotic exposure because probiotics are available without prescription.
Because IBD symptoms exist before IBD diagnosis, the potential for misclassification bias and reverse causality exists, when subjects might receive antibiotics as treatment of symptoms arising from as-yet undiagnosed IBD. This time between symptom onset and diagnosis was termed the “latency period.” To address this potential for reverse causality, we measured the latency period for subjects who developed IBD, defined as the time between the earliest potential IBD symptom diagnosis (eg, abdominal pain, diarrhea, abnormal weight loss) within 5 years before the first IBD diagnosis and the first IBD diagnosis. We then excluded all antibiotic exposures occurring from the median latency period until censorship for all subjects, regardless of outcome.
The primary outcome was IBD development, identified by using diagnosis codes (Read codes). Using the general practitioner survey report as the gold standard, IBD codes have been validated as having 92% positive predictive value for identifying patients with probable IBD in the General Practice Research Database.26 THIN data elements have been previously validated against this database for diverse associations.27
Demographic Characteristics and Comorbidities
We obtained demographic information, including age at cohort entry, gender, and Townsend score (a measure of socioeconomic deprivation quintile provided by THIN). The Townsend score was dichotomized as those with most deprivation versus all others. Race and ethnicity data are unavailable in THIN. We identified subjects diagnosed with conditions known to be associated with IBD that might also be associated with antibiotic exposure (ie, chronic granulomatous disease, primary sclerosing cholangitis, chronic osteomyelitis, Behçet’s syndrome, hyperimmunoglobulin E syndrome, common variable immunodeficiency).
Demographic variables were described by using frequencies, mean, median, and interquartile range (IQR) as appropriate. The association between each covariate and IBD was examined by using χ2 tests and logistic regression as appropriate, and all covariates associated with IBD development were examined, with P < .2 used for inclusion in the multivariable model. Confounders were defined as covariates whose adjustment produced an adjusted hazard ratio (aHR) >10% different than the unadjusted hazard ratio (HR).
A Cox proportional hazards model, which jointly considers IBD status and age of IBD onset, was used to examine the association between antibiotic exposure and IBD development. Each analysis was stratified according to primary care site, to allow separate hazard functions for each site because prescribing practices and patient populations may differ across sites. Because we were interested in estimating IBD risk at different ages, age was used as the time scale for cohort entry, exit, and exposure time points in multivariate analysis. Proportional hazards assumptions were checked statistically, and graphically if violated, and covariates were treated as time varying as appropriate.
Our primary a priori analysis examined the association between exposure to antianaerobic antibiotics and IBD development. Secondary analyses examined the association between antibiotic exposure and IBD development according to any antibiotic use and according to specific individual antibiotic classes. We additionally evaluated the primary association of interest for Crohn’s disease and ulcerative colitis separately. Sensitivity analyses were performed by using a 1-year latency period (when the proportion of subjects with previous diagnoses for possible IBD symptoms tapered to an apparent baseline) and assigning missing Townsend score data first all in the highest deprivation category and then all in the lower deprivation category.
A P value < .05 was considered significant for our primary analysis, and precise P values are reported. Stata version 11 (Stata Corp, College Station, TX) was used for all analyses.
There were 1 072 668 children eligible for inclusion. A total of 242 children were excluded: IBD was diagnosed before or within 6 months of cohort entry in 201 and 38, respectively; 3 had no IBD diagnosis date. The resulting cohort contained 1 072 426 subjects for analysis, followed up for 6.6 million person-years; 748 subjects (0.07%) developed IBD, for an overall incidence rate of 1.2/10 000 person-years (Tables 1 and 2). Only 30 of the 225 100 subjects followed up from birth developed IBD, precluding meaningful analyses in this subgroup.
Among the 748 subjects who developed IBD, the median latency period between first visit with a gastrointestinal diagnosis potentially consistent with IBD and first IBD diagnosis was 3.9 months (IQR: 0.5–17.9 months), and 68.2% had latencies ≤1 year (Fig 1). Fewer than 25% of other subjects had such diagnoses within 5 years before censorship.
Antibiotic exposures were common, with 57.7% of subjects exposed to at least 1 antianaerobic antibiotic and 64.0% exposed to any antibiotic. Subjects received a median of 1 week of antianaerobic antibiotics (IQR: 0–2), with 42.3% receiving none, 31.9% receiving 1 to 2 weeks, and 25.8% receiving >2 weeks. Exposure distribution according to courses received was similar, with 42.3% receiving none, 34.9% receiving 1 to 2 courses, and 22.8% receiving >2 courses. Excluding antianaerobic antibiotic courses during the latency period excluded exposures for 12.8% of subjects who developed IBD and 6.9% of those who did not.
Demographic Characteristics and Covariates Associated With IBD
Female gender and extreme socioeconomic deprivation were negatively associated with developing IBD (Table 3). Age at cohort entry was associated with developing IBD (HR: 1.21 per year [95% confidence interval (CI): 1.19–1.24]). Family history of IBD, chronic granulomatous disease, and primary sclerosing cholangitis were all significantly associated with an increased hazard of developing IBD, whereas other immunodeficiencies were not (P > .99). Calendar year of study entry was not associated with IBD (HR: 0.99 per year [95% CI: 0.97–1.02]). Hospital admissions were associated with increased IBD risk (HR: 1.08 per admission [95% CI: 1.06–1.09]). In multivariable analysis, country and hospital admissions altered the association between antianaerobic antibiotics and IBD by <1% and were not included in final models.
Primary Outcome: Association Between Antianaerobic Antibiotics and IBD
Although 0.07% of subjects developed IBD in both unexposed and exposed groups, the IBD incidence rates (calculated by using observed person-time as denominator) in these groups were 0.83 and 1.52/10 000 person-years, respectively, an absolute risk increase of 0.69 case/10 000 person-years and an 84% relative risk increase. In univariate analysis, any antianaerobic antibiotic exposure was associated with developing IBD (P < .001, log-rank test; Fig 2A), a dose-response effect existed (P < .001, log-rank test; Fig 2B), and this relationship remained significant throughout childhood (Fig 3).
Because the survival and hazard curves suggested that the relationship between any exposure to antianaerobic antibiotics and IBD decreased with increasing age, we incorporated an interaction term of log(age) to reflect this time-varying association. In multivariate analysis, exposure before 1 year of age was associated with a 5.5-fold increased IBD risk (aHR: 5.51 [95% CI: 1.66–18.28]) compared with those unexposed at that age, with decreasing IBD risks by 5 years (aHR: 2.62 [95% CI: 1.61–4.25]) and 15 years (aHR: 1.57 [95% CI: 1.35–1.84]).
Each antianaerobic antibiotic course was associated with a 6% increase in IBD hazard (aHR: 1.06 [95% CI: 1.04–1.08]), and each week of exposure was associated with a 1% increased hazard (aHR: 1.01 [95% CI: 1.00–1.02]). Exposure to >2 antianaerobic antibiotic courses by 1 year of age was more highly associated with IBD development than exposure to 1 to 2 courses, with an aHR of 4.77 (95% CI: 2.13–10.68) versus 3.33 (95% CI: 1.69–6.58).
Secondary Associations With IBD and Sensitivity Analyses
Adjusted associations between individual antibiotic class exposure measures and IBD are summarized without the other model covariates in Tables 4 and 5. Exposure to any antibiotic, penicillins, broad-spectrum penicillins, and cephalosporins were associated with IBD development, whereas macrolide, sulfonamide, and tetracycline exposures were not. Antibiotic classes to which <0.05% of subjects were exposed were not included.
The 748 subjects who developed IBD received 1642 antibiotic prescriptions; prescription day diagnoses included 440 (26.8%) head and neck infections, 292 (17.8%) chest infections, 163 (9.9%) skin and soft tissue infections, 82 (5.0%) genitourinary infections, 27 (1.6%) gastrointestinal infections, and 638 (38.9%) without specific infectious diagnoses. Excluding from the analysis those subjects who received antianaerobic antibiotics for gastrointestinal illnesses and developed IBD did not alter the primary association of interest.
Metronidazole or fluoroquinolone exposure was associated with the development of IBD. To determine whether outcome timing misclassification caused these associations, 2 additional sensitivity analyses were performed. For those who developed IBD and were prescribed metronidazole or fluoroquinolones in the year before diagnosis, we reassigned IBD outcome at the initial date of metronidazole or fluoroquinolone prescription. Doing so decreased the magnitude, but did not eliminate, the initial associations: any fluoroquinolone exposure (aHR: 2.09 [95% CI: 1.10–3.98]) and any metronidazole exposure (aHR: 186.25 [95% CI: 10.86–3193.65]) by 1 year of age remained significant. Lastly, among all 14 506 subjects who received metronidazole, prescription day diagnoses were available for 11 936 (82.2%): 34.9% had vaginitis/vaginosis, 15.1% had oral infections (eg, dental abscess), 7.5% had gastrointestinal symptoms, 3.1% had specific infections treated with metronidazole (eg, giardiasis), and 39.5% were for unclassifiable reasons (eg, “telephone encounter”).
Results among the 2 IBD subgroups, Crohn’s disease and ulcerative colitis, were similar to the main results but did not vary with age (Table 6). Sensitivity analysis assigning all missing Townsend scores first in the highest, then the lower, deprivation category did not alter the primary outcome. Sensitivity analysis using a 1-year latency period altered the primary outcome’s precision and magnitude but not its direction: any antianaerobic antibiotic exposure by 1 year of age (aHR: 3.73 [95% CI: 1.17–11.84]); each course (aHR: 1.03 [95% CI: 0.99–1.06]); exposure to 1 to 2 courses (aHR: 2.18 [95% CI: 1.07–4.50]); and exposure to >2 courses (aHR: 4.14 [95% CI: 1.74–9.87]).
Our large, population-based cohort study demonstrates that childhood antianaerobic antibiotic exposure was associated with IBD development; increasing early and cumulative exposure was directly associated with IBD risk. Our results suggest that for every 14 300 antianaerobic prescriptions given to children annually, 1 child will develop IBD. An estimated 49 million outpatient pediatric antibiotic prescriptions, approximately one-half for penicillins, occur in the United States annually; our data suggest those prescriptions would be associated with 1700 additional IBD cases yearly.28
Other findings in our multivariable model are consistent with previous epidemiologic IBD studies. Girls were 25% less likely to develop IBD as previously observed, and our overall IBD incidence rate (1.2/10 000 person-years) was similar to previous population-based incidence rate estimates among children.29 Those with extreme social deprivation were less likely to develop IBD, consistent with the “hygiene hypothesis” (persons living in cleaner environments may be more likely to develop autoimmune disorders).30 Chronic granulomatous disease and primary sclerosing cholangitis, both known to be associated with IBD, were highly but imprecisely associated with IBD, reflecting their rarity.31,32
Our findings are consistent with previous studies demonstrating increased IBD risk with earlier and cumulative antibiotic exposure.18,19,23 The odds of 1.3 for the association between antibiotic use and IBD among adults is consistent with the hazard we noted in older adolescents.20 Unlike Margolis et al,22 however, we found no association between tetracycline exposure and IBD, although our subject population was younger and less exposed to tetracycline.
Our findings are consistent with the hypothesis that antianaerobic antibiotic exposure might alter gut flora and trigger increased inflammation in some individuals. Alternately, if specific pathogens are directly associated with IBD, antibiotic exposure might simply be a surrogate for those infectious agents. One might expect such infections to be primarily gastrointestinal, yet only 1.6% of antibiotic prescriptions to our subjects who developed IBD were for gastrointestinal infections, and excluding these subjects from analysis did not change our results. Conversely, gut colonization with certain organisms might protect against IBD (as in the “hygiene hypothesis”), and antibiotic exposure could remove these protective organisms. If, for example, giardiasis protected against IBD, treatment with metronidazole would be associated with developing IBD. This outcome seems less likely, however, because of the multiple antibiotic classes associated with IBD and the diversity of diseases they treat.
A final explanation of our findings could be reverse causality; that is, misclassification of IBD diagnosis timing caused the association, with subjects receiving antianaerobic antibiotics as treatment of as-yet undiagnosed IBD. We addressed this outcome misclassification in multiple ways, however. First, we used validated IBD diagnosis codes.26,27 Second, we excluded antibiotic exposures for all subjects during the median diagnosis latency period. Sensitivity analysis extending this period to 1 year altered the findings’ precision and magnitude but not direction. Third, the strong association between metronidazole and fluoroquinolones (2 potential IBD treatments) and IBD development increased our concern about IBD diagnosis misclassification. However, these class-specific findings remained significant after re-setting the IBD outcome at the first metronidazole or fluoroquinolone prescription in the year before IBD diagnosis, suggesting that outcome misclassification was not significantly present. Moreover, metronidazole exclusively targets anaerobes and therefore may be more associated with developing IBD if altering anaerobes alone is more inflammatory than altering anaerobes and aerobes simultaneously. Lastly, penicillins and cephalosporins are not routine IBD treatments, yet were associated with IBD development. The associations between nonanaerobic antibiotic exposure and IBD may be because those antibiotics have activity against difficult-to-culture anaerobes or because IBD pathogenesis does not require perturbations in anaerobic flora alone.
Strengths of our study include its large, population-based nature; comprehensive exposure ascertainment; validated outcome ascertainment; and the results’ robustness to multiple sensitivity analyses. Although this study would be neither feasible nor ethical in prospective or randomized fashions, administrative data have limitations, with misclassification most concerning. We addressed outcome misclassification as mentioned earlier. We also addressed exposure misclassification by accounting for potential unobserved inpatient antibiotic exposures, without altering the main findings. We cannot account for antibiotics received without prescription, but this potential exposure misclassification should be nondifferential across our exposed and unexposed groups, biasing the results to the null. Finally, we did not observe all subjects’ entire childhood. The earliest subject enrollment was 1994, so subjects enrolled at birth could become as old as 15 years. Because most childhood antibiotics are prescribed in those <5 years of age,33,34 and median childhood IBD diagnosis age is 12 years,3 those enrolled from birth were younger, more likely to be antibiotic exposed, and less likely to develop IBD. Those enrolled later in childhood would be more likely to be classified as unexposed and to develop IBD. This possible exposure misclassification would therefore bias our results toward the null and may help explain the decreasing IBD risk we noted with increasing age.
Exposure to antianaerobic antibiotics during childhood was associated with development of the lifelong autoimmune condition IBD. Our study suggests that reduction in childhood antianaerobic antibiotic use may have the potential to help curb the rising incidence of childhood IBD. Many unanswered questions remain, however, such as whether specific difficult-to-culture organisms could play roles in either IBD pathogenesis or protection against IBD, and whether alteration of flora through antibiotic exposure alters the immune system directly.
The authors thank James D. Lewis, MD, MSCE (Department of Medicine, Department of Biostatistics and Epidemiology, and Center for Clinical Epidemiology and Biostatistics, Perelman School of Medicine at the University of Pennsylvania, Philadelphia, PA) for his thoughtful help with the conception and design of this work.
- Accepted May 17, 2012.
- Address correspondence to Matthew P. Kronman, MD, MSCE, Division of Infectious Diseases, Seattle Children’s Hospital, 4800 Sandpoint Way NE, M/S R-5441, Seattle, WA 98105. E-mail:
Drs Kronman, Zaoutis, Haynes, Feng, and Coffin participated in study conception and design; Dr Haynes conducted acquisition of data; Drs Kronman, Zaoutis, Feng, and Coffin assisted with analysis and interpretation of the data; Dr Kronman drafted the manuscript; Drs Zaoutis, Haynes, Feng, and Coffin performed critical revision of the manuscript for important intellectual content; Drs Kronman and Feng conducted statistical analysis; Drs Kronman, Zaoutis, and Coffin obtained funding; Dr Haynes contributed administrative, technical, or material support; and Drs Zaoutis and Coffin provided supervision of study conduct and analysis.
The content is solely the responsibility of the authors and does not necessarily represent the official views of the Agency for Healthcare Research and Quality or the National Institutes of Health. All authors had full access to all of the data in the study and the corresponding author (Dr Kronman) takes responsibility for the integrity of the data and the accuracy of the data analysis.
FINANCIAL DISCLOSURE: Dr Zaoutis has received research support from Merck and Company, Inc for other projects and has served as a consultant for Pfizer, Inc, HemoCue, Inc, Merck and Company, Inc, Cubist Pharmaceuticals, and Astellas Pharma US, Inc. He has also received a lecture honorarium from Merck and Company, Inc. Dr Coffin previously received research support from Sage Pharmaceuticals, Inc for other projects. Drs Kronman, Haynes, and Feng have indicated they have no financial relationships relevant to this article to disclose.
FUNDING: Dr Kronman received support from a National Research Service Award Institutional Research Training grant T32-HD043021 at The Children's Hospital of Philadelphia as well as a University of Pennsylvania Center for Education and Research on Therapeutics grant U18HS016946 from the Agency for Healthcare Research and Quality. This work was also supported by a National Institutes of Health Clinical and Translational Science Award (grant UL1-RR024134 to Dr Haynes). No funding agency participated directly in the study design or conduct; data collection, management, analysis, or interpretation; or manuscript preparation, review, or approval. Funded by the National Institutes of Health (NIH).
- Barrett JC,
- Hansoul S,
- Nicolae DL,
- et al.,
- NIDDK IBD Genetics Consortium,
- Belgian-French IBD Consortium,
- Wellcome Trust Case Control Consortium
- Tysk C,
- Lindberg E,
- Järnerot G,
- Flodérus-Myrhed B
- Jess T,
- Simonsen J,
- Nielsen NM,
- et al
- Ott SJ,
- Musfeldt M,
- Wenderoth DF,
- et al
- Sokol H,
- Pigneur B,
- Watterlot L,
- et al
- Virta L,
- Auvinen A,
- Helenius H,
- Huovinen P,
- Kolho KL
- Card T,
- Logan RF,
- Rodrigues LC,
- Wheeler JG
- Hviid A,
- Svanström H,
- Frisch M
- ↵CSD Medical Research UK. THIN data statistics. Available at: http://csdmruk.cegedim.com/. Accessed December 15, 2011
- Hersh AL,
- Shapiro DJ,
- Pavia AT,
- Shah SS
- ↵Lichtenstein DR. Hepatobiliary Complications of inflammatory bowel disease. Curr Gastroenterol Rep. 2011;13(5):495–505
- ↵Marks DJ, Miyagi K, Rahman FZ, Novelli M, Bloom SL, Segal AW. Inflammatory bowel disease in CGD reproduces the clinicopathological features of Crohn's disease. Am J Gastroenterol. 2009;104(1):117–124
- Copyright © 2012 by the American Academy of Pediatrics