The impact of benzodiazepines on occurrence of pneumonia and mortality from pneumonia: a nested case-control and survival analysis in a population-based cohort
- 1Division of Epidemiology and Public Health, University of Nottingham, Nottingham, UK
- 2Department of Leucocyte Biology and Anaesthetics, Intensive Care and Pain Medicine, Imperial College London, London, UK
- 3Wellcome Department of Imaging Neuroscience, University College London, London, UK
- Correspondence to Dr Robert D Sanders, Wellcome Department of Imaging Neuroscience, Institute of Cognitive Neuroscience, London, WC1N 3AR, UK;
- Received 3 July 2012
- Accepted 5 October 2012
- Published Online First 5 December 2012
Objectives Benzodiazepines have been associated with an increased incidence of infections, and mortality from sepsis, in the critically ill. Here, we determined the effect of community use of benzodiazepines on the occurrence of, and mortality following, pneumonia.
Methods A nested case-control study using 29 697 controls and 4964 cases of community-acquired pneumonia (CAP) from The Health Improvement Network, a UK primary care patient database (2001–2002), investigated the association between benzodiazepines and pneumonia occurrence using conditional logistic regression. Cox regression was then used to determine the impact of benzodiazepines on mortality in the 4964 cases of CAP. Results are presented as adjusted OR, adjusted HR and 95% CI.
Results Exposure to benzodiazepines was associated with an increased risk of pneumonia (OR 1.54, 95% CI 1.42 to 1.67). Individually diazepam, lorazepam and temazepam, but not chlordiazepoxide, were associated with an increased incidence of CAP. As a class, benzodiazepines were associated with increased 30-day (HR 1.22 (95% CI 1.06 to 1.39)) and long-term mortality (HR 1.32 (95% CI 1.19 to 1.47)) in patients with a prior diagnosis of CAP. Individually diazepam, chlordiazepoxide, lorazepam and temazepam affected long-term mortality in these patients.
Conclusions Benzodiazepines were associated with an increased risk of, and mortality from, CAP. These hypothesis generating data suggest further research is required into the immune safety profile of benzodiazepines.
What is the key question?
Do benzodiazepines increase the occurrence of, or mortality from, pneumonia?
What is the bottom line?
Benzodiazepine exposure was associated with an increased risk of developing pneumonia and dying from pneumonia.
Why read on?
Our data indicate a significant risk of benzodiazepine exposure on infectious lung disease which may be considered a modifiable risk factor in future studies.
In the UK and USA, approximately 2% of the population have taken benzodiazepines for 12 months or more1 and their use is even more prevalent in the elderly patients (up to 10%).2 Benzodiazepines are widely used for anxiety, epilepsy, muscle spasm, alcohol withdrawal, palliation, insomnia and to provide sedation. These clinical effects are transduced by activation of inhibitory γ-amino butyric acid type A receptors (GABAA) in the brain.
Little is known about the immune effects of benzodiazepines. Animal studies suggest that diazepam increases susceptibility to infection, including pneumonia,3–7 perhaps via activation of GABAA receptors on immune cells.7 Limited clinical data support a hypothesis that benzodiazepines increase susceptibility to infection. It has been shown that avoidance of benzodiazepine sedation (lorazepam) by use of a non-GABAergic sedative improves survival in critically ill septic patients by 70%.8 ,9 A similar study showed that sedation with the benzodiazepine midazolam doubled the risk of secondary infections in critically ill patients compared with non-GABAergic sedation.10 However, a more recent study did not observe the difference in infection rates.11 Furthermore, these studies used high doses of benzodiazepines to achieve sedation and so it is unclear whether these data have relevance to community use of these drugs.
The influence of subsedative doses of benzodiazepines on infectious outcomes is also incompletely evaluated in humans. One group reported that benzodiazepines and/or antidepressants increase mortality from pneumonia in over 60-year-olds12 but not over 80-year-olds.13 Additionally, a randomised controlled trial of diazepam therapy for acute stroke reported an increased incidence of pneumonia with diazepam treatment.14 In contrast, another group has suggested that benzodiazepines may be associated with a decreased incidence of pneumonia; however, this study used questionnaires to determine drug exposure and thus is prone to recall and reporting bias.15 A case-control study of patients over 65-years-old also did not report an increased risk of pneumonia with exposure to benzodiazepines though they report wide CI due at least in part to the low numbers of patients taking benzodiazepines in this study.16 Furthermore, all these studies were small and the outcomes varied and limited; for example, currently it is unclear whether there is any long-term effect on mortality.
We hypothesised that benzodiazepines would increase the odds of pneumonia and mortality from pneumonia. We also performed analyses of the non-benzodiazepine, zopiclone, as it also acts by activation of GABAA receptors but lacks many of the other molecular effects of benzodiazepines.17 Testing zopiclone was intended to provide preliminary information about the safety of GABAergic drugs as a group. We focus here on the risk of pneumonia, and 30-day and long-term mortality from pneumonia.
The Health Improvement Network (THIN) database is a large, longitudinal collection of records of patients registered to various primary care facilities all over the UK. This contains about 9.1 million registered patients, of which 3.4 million are categorised as ‘active patients’.18 Data are entered using Read codes that map onto the International Classification of Disease-9 codes. Ethical approval for this study was obtained from the Scientific Review Committee (SRC) (11-010R). We conducted two analytical studies within the THIN population-based cohort: (i) a matched nested case-control analysis to address whether the occurrence of pneumonia was associated with benzodiazepine exposure using conditional logistic regression and (ii) a survival analysis using Cox regression to investigate the association of benzodiazepine exposure with mortality in pneumonia patients.
Cases were patients of all ages in the database with a diagnosis of pneumonia as identified by the patient's medical records occurring between 1 July 2001 and 1 July 2002. We used only the first recorded pneumonia diagnosis within this period for each case. This ensured that only ‘incident’ cases were picked. Identification of cases was done using specific medical Read codes corresponding to a pneumonia diagnosis (available on request) with a designated index date. For each case, six controls were matched by practice, sex and age at index date (within 3 years). The generally accepted rule is that more than four to five controls rarely add to the efficiency of a study. We included six controls because the cost of including additional controls is negligible (an advantage provided by THIN), and there is concern for sufficient numbers in stratified analyses.19 Both cases and controls required at least a year's worth of prior data to be included in the study to allow for adequate capture of confounding variables. We did not exclude controls or cases if they had a recorded pneumonia diagnosis before our selected study period and included previous pneumonia episodes as a covariate in our analyses.
All recorded prescriptions of benzodiazepines and zopiclone were used. Exposure to drug was classified as ‘current’ when the most recent prescription was within 30 days before the pneumonia index date. Prescriptions within 31–90 days before the index were treated as ‘recent’ exposures, while prescriptions of 90 days or more before the index date were treated as ‘past’ exposure. An additional category was created for ‘no use’ where subjects had never been prescribed any benzodiazepine. These categories were determined by the British National Formulary prescription guidelines and confirmed by plotting histograms to check prescription lengths and patterns. A supplementary analysis was conducted to understand the odds of pneumonia in patients who were chronically on benzodiazepines. These patients were defined as having prescriptions both in the 30- and 90-day periods before the pneumonia event date.
We evaluated the following comorbidities as potential confounders: ischaemic heart disease, pulmonary disease (including chronic obstructive pulmonary disease and asthma) and previous pneumonia episodes. In addition, we used a combined weighted comorbidity index, the Charlson index adapted for use with International Classification of Disease-9 codes.20 ,21 Other potential confounders considered include alcohol consumption (categorised as above or at/below recommended weekly units for both male and female subjects), diagnosis of depression and psychosis, current smoking (the most recent record of smoking status was used) and socioeconomic status measured using the Townsend deprivation score quintiles (the first quintile being the least deprived, while the fifth quintile being the most deprived).22
Conditional logistic regression was used to investigate the association between exposure to the different drugs and risk of pneumonia. Cox regression was used to assess the association between drug exposure and mortality following confirmation that the proportional hazard assumptions were met. Two mortality endpoints were investigated, 30-day and long-term mortality. Our primary interest was the effect of benzodiazepines as a group on the occurrence of, and mortality following, pneumonia. However, we also looked at the effects of individual benzodiazepines and zopiclone. Results have been expressed as OR or HR with 95% CI. The multivariable model included all variables that were either a significant risk factor for pneumonia in the univariate analysis (p<0.05) or were found to modify the OR for drug association by at least 10% when included in a model with the main exposure variable. Given that age, comorbidity12 ,13 and gender23 may influence the effect of benzodiazepines, we tested for interactions among benzodiazepines and age, gender or Charlson comorbidity status. All analyses were carried out in Stata11.24 The data are reported in accordance with STROBE guidance.
The sample used in this research had a total of 34 661 patients, of whom there were 29 697 controls and 4964 cases. The characteristics of cases and controls are summarised in table 1. Cases were more likely to have had previous episodes of pneumonia, myocardial infarction, a history of depression, a history of psychosis related diagnosis, be current smokers and have a higher Charlson's comorbidity index score as compared with controls.
Benzodiazepines and risk of pneumonia
Table 2 shows the adjusted and unadjusted OR for benzodiazepine as a class, and as individual benzodiazepines. After adjusting for current smoking, presence of lung disease, Townsend deprivation index, diagnosis of depression or psychosis, Charlson's comorbidity index, myocardial infarction and previous episode of pneumonia a significant association was seen between benzodiazepine use and increase in pneumonia risk (adjusted OR 1.54, 95% CI 1.42 to 1.67). After adjustment for confounders, prescriptions of diazepam, lorazepam and temazepam were associated with an increased risk of pneumonia (table 2). However, we did not find a statistically significant association between current use of chlordiazepoxide and pneumonia risk (table 2). Furthermore, use of zopiclone, a non-benzodiazepine drug acting on GABAA receptors, similarly showed a higher risk of pneumonia (table 2). Refinement of the analysis to look at current prescriptions showed that diazepam, temazepam and zopiclone were all associated with higher odds of pneumonia. A supplementary analysis was conducted to understand the odds of pneumonia in patients who were chronically on benzodiazepines. These patients were defined as having prescriptions both in the 30- and 90-day periods before the pneumonia event date. They were compared with non-benzodiazepine users. Chronic users represented the majority of the cohort and therefore it is unsurprising that the results were similar in this analysis to the primary analysis (see online supplementary table 1a).
A test for interaction was carried out with age, gender and Charlson's comorbidity index score for the main drug exposure, benzodiazepine use. No interactions were found with age or gender for drug exposure; however, a statistically significant interaction was found by Charlson's comorbidity index score (p<0.001). Table 3 shows the stratified analysis results. The association between benzodiazepine use and pneumonia incidence was stronger in people with lower comorbidity scores.
Benzodiazepines and mortality following pneumonia
We next analysed whether benzodiazepine exposure increased long-term mortality (median follow-up 2.8 years) following pneumonia and found that current and recent prescriptions exerted the greatest effect (table 4). Individually all drugs tested, except zopiclone, were associated with higher long-term mortality. However, 30-day mortality was only affected by diazepam and lorazepam (table 4). When individual benzodiazepine use was divided by the timing of prescription, only current diazepam prescription appeared to influence mortality (table 4). Again, we conducted a supplementary analysis to understand whether the odds of death following pneumonia in patients who had prescriptions in ‘chronic users’ were different compared with those who did not take benzodiazepines. The results were similar and largely unchanged by this analysis (see online supplementary table 1b). A supplementary propensity score analysis based on the medical history was also conducted to help limit concerns over residual confounding. Using this technique, the results were largely the same; however, zopiclone also increased long-term mortality (see online supplementary table 1c). Using the whole dataset we also investigated whether the mortality effect of the class of benzodiazepines was affected by age, gender or comorbidity; similar to the effect on the incidence of pneumonia, we found that increasing comorbidity reduced the impact of benzodiazepines on mortality (table 5).
Herein we have demonstrated that benzodiazepine and zopiclone exposure was associated with an approximately 50% increase in risk of pneumonia with a similar associated increase in mortality from pneumonia. Interestingly, the non-benzodiazepine zopiclone was associated with a higher risk of pneumonia. We adjusted for presence of lung disease in cases and controls, previous pneumonia, depression, myocardial infarction, socio-economic factors, comorbidity, alcohol use and smoking status. However we cannot exclude the effect of unmeasured confounding, though we have adjusted for many factors in this analysis. For example, while we have adjusted for psychiatric and medical comorbidities, we cannot exclude further confounding by indication. In order to better understand the benzodiazepine effect we investigated whether age, gender or comorbidity influenced the results, finding that increasing comorbidity reduced the impact of the drug on risk of pneumonia and mortality from pneumonia. Finally, we must stress that our data cannot be used to definitively link cause and effect, though the findings are in line with our hypothesis and preclinical data. Prospective cohort studies are now required to investigate our hypothesis further.
This is a large population-based case-control and cohort study with 29 697 controls and 4946 cases. We have used THIN database that is representative of the UK population; previous results from research using THIN have been in keeping with research from other large databases.25 This makes our findings applicable to the general population and ensures no recall bias as exposures were recorded prospectively before the diagnosis of pneumonia. The mode of collection of data (ie, routinely, non-interventional and in a prospective manner prior to the diagnosis of pneumonia) also ensures that reporting bias is eliminated at the point of data collection. The prospective data collection ensures the avoidance of temporal bias as relating to the timing of the exposure to benzodiazepines. For the case-control study, cases and controls were individually matched by age, sex and location of practice, therefore, reducing the confounding effects of these variables. In particular, due to the prescribing policies and variations in different regions, matching by practice harmonises this factor.
A few limitations to our study should be highlighted. One is the possibility of selection bias; while we have adjusted for psychiatric and medical comorbidities, other patient factors associated with benzodiazepine use may account for the difference. Likewise, we cannot be totally confident that we have removed all residual confounders. Another possibility is an error of misclassification of pneumonia that may arise from differences in coding practices. While these data reflect infectious lung disease,26–28 and thus are adequate to investigate the effect of these drugs on a clinical infection, the absence of a confirmatory chest x-ray questions the actual diagnosis of pneumonia. However, if this happens, the results would be biased towards unity because the misclassification should be independent of exposure status. Codes for ‘acute lower respiratory tract infection’ (ALRTI) have been included in the case definition as some cases of pneumonia may be coded as ALRTI in the absence of a confirmatory radiograph (see supplementary file). Nonetheless, our findings suggest that the clinical diagnosis of pneumonia and the associated mortality are both affected by the drugs. A related concern may be whether severity of pneumonia should be adjusted for; however, these data are not available within the dataset. Nonetheless, severity of pneumonia may well be dependent on exposure status to benzodiazepines (given that risk of, and mortality following, pneumonia appear to be dependent) and therefore even if the data on the severity of pneumonia were available it would not be sensible to adjust for this. In order to avoid counting a single case multiple times (eg, if the patient reconsults for the same episode of illness), only the first ever record of pneumonia diagnosis occurring within the study period was considered.
Differences in coding practices of some of the confounding variables might also bring in some form of bias in the study. However, this misclassification should be independent of exposure and outcome status; thus, bias will be non-differential. The Charlson's comorbidity index score is calculated by weighting disease burdens with mortality being the outcome;20 ,21 ,29 however, the index has previously been validated in the literature for use in morbidity studies.20 ,29 Furthermore, we augmented this index with further variables at our disposal to limit residual confounding. Another possible limitation it that we included covariates in the final model based on our initial bivariate analysis; this may be inferior to including all variables identified a priori.
Prescription data were used as a proxy for exposure to the drugs; however, we cannot be certain that the drugs were taken as recommended by the general practitioner. This might lead to an overestimation of exposure to these drugs in the analysis. This limitation, however, would be non-differential and would again bias the results towards unity (ie, reduce the benzodiazepine effect). However, analysis of ‘chronic users’ showed similar effects to the overall cohort (see online supplementary table 1a,b).
We suspect that chlordiazepoxide did not affect the incidence of pneumonia as it is used in the treatment of alcohol dependence and we adjusted for this factor; hence the impact of chlordiazepoxide was dwarfed by the underlying disease. When current prescriptions (those issued within 30 days of the pneumonia) were analysed, diazepam, temazepam and zopiclone were associated with a higher risk of pneumonia; however, lorazepam was not. Nonetheless, the point estimate for the effect of a current prescription of lorazepam on the incidence of pneumonia remained above one and lorazepam was associated with increased mortality when we focused on chronic users of the drug (see online supplementary table 1b). Analysis of the long-term mortality data showed a similar pattern of increased mortality for each benzodiazepine. However, zopiclone did not affect mortality in the primary analysis, and we suspect that this relates to the low usage of zopiclone in the population of ‘cases’ though confirmation in future studies is required. Further investigation is required because zopiclone was associated with an increase in long-term mortality following pneumonia in the propensity score matched analysis (see online supplementary table 1c).
The effect of comorbidity on the benzodiazepine-associated change in risk of, and mortality from, pneumonia may be explained by the influence of cumulative disease processes outweighing the effects of the drug. Another explanation is that inflammation from comorbid diseases influences the expression of the GABAA receptor, the putative immune target for benzodiazepines. We have recently observed that certain inflammatory cytokines and stimuli reduce GABAA receptor expression on alveolar macrophage (Sanders et al, unpublished observations). A reduction in expression of the receptor target, by comorbid inflammation, would be expected to reduce the immune effects of benzodiazepines as we have observed.
Our findings are complemented by some previous studies. The increase in risk of pneumonia is similar to a previous study by Knol et al30 which recorded a 60% increased risk during the first week of exposure to ‘anti-psychotic’ medication. Vergis et al31 showed that the use of ‘tranquilisers’ increased mortality from pneumonia in patients in long-term care. Using a retrospective cohort study design, Hak et al12 followed up 229 people over 60 years of age with pneumonia over a period of 2 years and found that the use of benzodiazepines and/or antidepressants increased the risk of pneumonia by 89%. The same group performed a similar study in over 80-year-olds which showed a conflicting result; benzodiazepine and/or antidepressant use was not associated with hospital admission or death in this cohort with pneumonia after multivariable analysis.13 These two contrasting findings from the same group12 ,13 may be consistent with our observations, as elderly patients likely have greater comorbidity, reducing the impact of benzodiazepines on mortality (table 5). In contrast, our finding that patients with less comorbid disease are vulnerable to benzodiazepine-induced susceptibility to pneumonia are in line with effects on a relatively younger population.12
Further studies are needed to confirm our findings. Prospective cohort studies are warranted to understand the impact of benzodiazepines on outcomes from infection. Further studies should also confirm the safety of other drugs, such as the antiepileptic medication topiramate that targets GABAA receptors.32 Randomised controlled trials are potentially difficult given the lack of therapeutic alternatives at present. However, benzodiazepines that are selective for GABAA receptor subtypes are under development33 and may have a superior immune profile to the current non-selective drugs (Sanders et al, unpublished observations). In intensive care, where patients are at high risk of infection and large doses of benzodiazepines are used, it is possible to randomise patients to sedation with a GABAergic or non-GABAergic sedative8–10 and randomised controlled trials should be considered in this context.
Benzodiazepines and zopiclone are commonly prescribed medications that have significant immune effects. Our data herein suggest that they may increase both the risk of and mortality from pneumonia. This is consistent with data from clinical trials8–10 ,34 and concerns expressed over the intensive care unit effects of these drugs leading to movement away from benzodiazepine sedation.35 Nonetheless, given the widespread use of benzodiazepine drugs, further studies are required to evaluate their safety in the context of infection.
RDS devised the hypothesis. All authors took part in study design and interpretation. EO and PM performed the analyses. PM had full access to all of the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis. RDS and EO wrote the paper with input from RH and PM. All authors approved the manuscript. We thank Christopher JP Smith, Division of Epidemiology and Public Health, University of Nottingham, for the initial data extraction. RDS is supported by the Medical Research Council, Swindon, UK (RDS: G0802353). RH is supported as the GSK/British Lung Foundation Chair of Respiratory Epidemiology. RDS has received an honorarium from Hospira, Lake Forest, Illinois, USA, for speaking on the immune effects of sedative agents at the Canadian Society of Anesthesiologists meeting 2011. There is no ongoing relationship. There was no industrial involvement in this study.
RDS and PRM contributed equally to this work.
Funding Funding came from departmental resources. RDS is supported by the Medical Research Council, Swindon, UK (RDS: G0802353). RH is supported as the GSK/British Lung Foundation Chair of Respiratory Epidemiology.
Competing interest All authors have completed the Unified Competing Interest form at www.icmje.org/coi_disclosure.pdf (available on request from the corresponding author) and declare: no support from any organisation for the submitted work; RDS received an honorarium for a talk at the Canadian Society of Anaesthesiologists from Hospira in the previous 3 years (there is no ongoing relationship and there was no industrial involvement in this work); no other relationships or activities that could appear to have influenced the submitted work.
Ethics approval SRC (11-010R).
Provenance and peer review Not commissioned; externally peer reviewed.