Antipsychotic drug exposure and risk of pulmonary embolism: a population-based, nested case–control study
© Conti et al.; licensee BioMed Central. 2015
Received: 7 January 2015
Accepted: 23 April 2015
Published: 29 April 2015
Only three observational studies investigated whether exposure to antipsychotics is associated with an increased risk of pulmonary embolism, with conflicting results. This study was therefore carried out to establish the risk of pulmonary embolism associated with antipsychotic drugs, and to ascertain the risk associated with first- and second-generation antipsychotic drugs, and with exposure to individual drugs.
We identified 84,253 adult individuals who began antipsychotic treatment in a large Italian health care system. Cases were all cohort members who were hospitalized for non-fatal or fatal pulmonary embolism during follow-up. Up to 20 controls for each case were extracted from the study cohort using incidence density sampling and matched by age at cohort entry and gender. Each individual was classified as current, recent or past antipsychotic user. The occurrence non-fatal or fatal pulmonary embolism was the outcome of interest.
Compared to past use, current antipsychotic use more than double the risk of pulmonary embolism (odds ratio 2.31, 95% confidence interval 1.16 to 4.59), while recent use did not increase the risk. Both conventional and atypical antipsychotic exposure was associated with an increase in risk, and the concomitant use of both classes increased the risk of four times (odds ratio 4.21, 95% confidence interval 1.53 to 11.59).
Adding the results of this case–control study to a recent meta-analysis of three observational studies substantially changed the overall estimate, which now indicates that antipsychotic exposure significantly increases the risk of pulmonary embolism.
Pulmonary embolism (PE) is a significant cause of morbidity and mortality, occurring at an estimated 95 cases per 100,000 patient-years [1,2]. The severity of PE ranges from asymptomatic, incidentally discovered subsegmental thrombi to hemodynamically unstable, life-threatening episodes, and sudden death [3,4].
In 1997 Walker and colleagues conducted an epidemiological study suggesting that exposure to clozapine, an agent belonging to the group of the so-called second-generation antipsychotic (AP) drugs, significantly increased the risk of PE mortality . Since then, other epidemiological cohort and case control studies provided additional data on this association . A recent systematic review and meta-analysis identified 13 studies, of which only three investigated PE outcomes, while ten estimated the association between AP exposure and risk of a composite outcome which included deep venous thrombosis, femoral vein thrombi, popliteal vein thrombi, iliac vein thrombi, deep vessels of lower extremity thrombophlebitis, and pulmonary embolism . Analysis of the three PE studies failed to detect a significant association between exposure to APs and risk of PE outcomes (odds ratio 4.90, 95% confidence interval 0.77 to 30.98), but the confidence interval was very wide and included the possibility of substantial harm .
Therefore, we investigated whether exposure to AP drugs is associated with an increased risk of PE, and ascertained the risk associated with first- and second-generation AP drugs, and with exposure to individual drugs. We carried out a nested case–control study using a large administrative database in the health system of Lombardy, a region located in northern Italy. As secondary aim, we updated the published meta-analysis with the results of this study, using the same methodology.
Setting and data source
Lombardy is the largest and the most affluent region in Italy, with a population of around 10 million inhabitants. Lombardy is located in the northernmost part of the country and includes the metropolitan area of Milan, Italy’s second largest city.
The data used for this study were retrieved from the Regional Health Service (RHS) databases of Lombardy [8,9]. These databases include: i) demographic and administrative information on all residents in the Lombardy region; ii) all community (outside the hospital) drug prescriptions reimbursed by the RHS; iii) all public and private hospital discharge forms, with diagnoses according to the ICD-9-CM. According to local regulations on administrative database analyses, no formal approval of the study protocol was required; however, as previously reported for similar analyses , to preserve patient privacy, the identification codes from all the databases were converted to anonymous codes, and the conversion table was then destroyed.
Study population and design
From the regional prescription database we identified all individuals aged ≥ 18 starting a new treatment with AP drugs from 1 January 2012 to 31 December 2013. AP drugs were defined as all drugs belonging to the N05A Anatomical Therapeutic Chemical Classification (ATC) group (with the exception of N05AN, lithium). New AP users were individuals with no AP prescriptions in the 12 months before the first prescription issued in the study period. From the group of new AP users, patients with a recorded diagnosis from hospital discharge forms of neoplasm (ICD9-CM 140–239), hip fracture (81.xx), pulmonary embolism or deep vein thrombosis (415.xx, 453.xx, 451.1x) in the year before the first AP prescription were excluded.
Each member of the cohort was followed from the first AP prescription until the earliest of the following event: outcome of interest (PE), death for any cause, emigration and end of follow-up.
Selection of cases and controls
Cases were all cohort members who were hospitalized for non-fatal or fatal PE during follow-up (ICD9-CM codes 415.xx). For cases with more than one hospital admission for PE during follow-up, we selected the first record of PE. Up to 20 controls for each case were extracted from the study cohort using incidence density sampling and matched by age at cohort entry and gender. All controls were alive and at risk of the outcome of interest at the date of the first PE in their matched case (index date). Cases and controls with pregnancy (ICD9-CM codes 630–677), leg/hip fracture or a diagnosis of neoplasm in the three months before the index date were excluded. In agreement with Parker and colleagues , controls with any prescriptions of warfarin or other antithrombotic drug use (ATC code B01AA03/B01AA07) before the index date were excluded. For cases, prescriptions for warfarin or other antithrombotic drugs in the six weeks before the index date did not trigger exclusion as they could plausibly be treatment for the index event itself ; however we excluded cases with any use of warfarin or other antithrombotic drugs earlier than this six week period.
Assessment of AP exposure
We assessed exposure to AP drugs on the basis of prescriptions on or before the index date. Each individual was classified for exposure to AP as current user (one or more prescriptions within three months before index date); recent user (one or more prescriptions in the period between 4 and 12 months before index date); past user (one or more prescriptions 13 months or more before index date). Exposure during the past three months was classified according to AP type (conventional only, atypical only, both), formulation (oral, acute injection, long-acting injection), number of AP scripts, number of different AP prescribed, and overall amount of AP exposure, calculated by dividing the total amount of prescribed drug by the recommended daily dose, according to each agent’s defined daily dose (DDD). The DDD is a theoretical unit of measurement defined as the assumed average maintenance daily dose for a drug, used for its main indication in adults .
Using the whole study cohort, we estimated the rate ratios with 95% confidence intervals (CIs) for the association between PE and AP exposure, taking past use as the reference category.
For the case–control analysis, cases and controls were compared using chi-square test, for categorical data, and t test for continuous data. We undertook multiple conditional logistic regression to estimate the odds ratio (OR), with 95% CIs, for the association between PE and AP use. Subgroup analyses were performed for age (≥65 years, < 65 years) and gender. ORs were adjusted for the following confounding variables: psychiatric hospitalization (ICD9-CM code 290–319); hospitalization for coronary heart disease (410.xx-414.xx), heart failure (428.xx), stroke (431.xx, 433.xx-436.xx), chronic liver disease (570–573), chronic renal disease (584–586), Parkinson’s disease (ICD9-CM code 332 or drug prescription with ATC code N04), diabetes (ICD9-CM code 249–250 or drug prescription with ATC code A10); prescription of lithium (ATC code N05AN01), statins (C10AA), NSAIDs (M01A), antiplatelet treatment (A01AD05, B01AC06, N02BA01, B01AB, B01AC04, B01AC05), and antihypertensive agents (C02, C03,C07,C08,C09). All covariates were considered in the 12 months before index date. We did not adjust for varicose veins, inflammatory bowel disease, and use of hormone replacement therapy, as suggested by Parker , as only very few cases with these disorders were recorded.
Numbers needed to harm and extra cases per 10 000 treated were estimated, using the incidence of PE in the past users and the adjusted ORs from current users.
Finally, in order to provide a summary of the whole evidence base, we added the results of the present case–control study to a recently published meta-analysis of observational studies on the risk of PE outcomes associated with AP exposure (Additional file 1) . A meta-analysis influence test that eliminated each of the included studies one at a time was carried out to test for possibly excessive influence of individual studies.
Data management and analyses were carried out using SAS statistical package (SAS Institute, Cary, NC), version 9.2 and Comprehensive Meta-Analysis version 2: the PHREG procedure with STRATA statement was used to fit the conditional logistic regression models. All tests were two-tailed, and results were considered significant at p < 0.05.
Cases and controls
Characteristics of cases (patients with first record of pulmonary embolism) and matched controls at index date
Cases (n = 232)
Controls (n = 4353)
Mean age ± SD
77.1 ± 13.5
76.9 ± 13.7
Age ≥ 65
Age < 65
Coronary heart disease
Chronic liver disease
Chronic renal disease
Risk associated with AP drugs
Risk of pulmonary embolism associated with antipsychotic drug use
Cases (n = 232)
Controls (n = 4353)
Unadjusted odds ratio (95% CI)
Adjusted* odds ratio (95% CI)
Category of antipsychotic exposure:
2.91 (1.48 - 5.73)
2.31 (1.16 - 4.59)
1.24 (0.58 - 2.64)
0.96 (0.44 - 2.07)
Class of antipsychotic received: (within past 3 months)
4.67 (2.27 - 9.61)
3.52 (1.69 - 7.35)
2.49 (1.26 - 4.93)
2.01 (1.01 - 4.03)
5.66 (2.11 - 15.17)
4.21 (1.53 - 11.59)
No of antipsychotic scripts received: (within past 3 months)
2.82 (1.41 - 5.65)
2.19 (1.08 - 4.42)
3.01 (1.5 - 6)
2.44 (1.21 - 4.93)
No of different antipsychotic drugs received: (within past 3 months)
2.74 (1.39 - 5.4)
2.19 (1.1 - 4.36)
7.68 (3.32 - 17.75)
5.87 (2.48 - 13.89)
DDDs: (within past 3 months)
2.48 (1.23 - 5)
1.96 (0.96 - 3.99)
3.31 (1.66 - 6.6)
2.62 (1.3 - 5.25)
Formulation: (within past 3 months)
2.95 (1.5 - 5.81)
2.34 (1.17 - 4.64)
2.28 (0.75 - 6.95)
2.05 (0.66 - 6.37)
As the vast majority of cases and controls received oral formulations, we were not able to calculate precise estimates for acute and long-acting injections.
Risk of pulmonary embolism associated with antipsychotic use: subgroup analysis for age and gender
Adjusted* odds ratios (95% confidence interval) for current use vs past use for each subgroup
age < 65
age ≥ 65
0.77 (0.08 - 7.17)
4.96 (1.55 - 15.9)
3.94 (1.43 - 10.88)
0.5 (0.13 - 1.93)
1.7 (0.39 - 7.31)
1.03 (0.4 - 2.68)
0.6 (0.19 - 1.89)
3.65 (1.48 - 9.04)
2.31 (1.16 - 4.59)
Numbers needed to harm and extra cases per 10 000 current users of antipsychotics over one year
Number needed to harm (95% confidence interval)
Extra cases per 10 000 treated (95% confidence interval)
624 (404 to 1,365)
16 (7 to 25)
age ≥ 65
244 (178 to 388)
41 (26 to 56)
326 (237 to 520)
31 (19 to 42)
age ≥ 65 and women
160 (121 to 236)
62 (42 to 83)
Meta-analysis of observational studies
This nested case–control study, based on administrative and clinical records of individuals resident in the region of Lombardy in northern Italy, found a more than double increased risk of PE for individuals prescribed AP drugs in the past three months. The increase in risk was more than three-fold for individuals aged 65 or above, almost four times for women, and almost five time for women aged 65 or above. The overall absolute risk, however, was relatively low, with an excess of 16 extra cases of PE per 10 000 patients treated over one year; however, in women aged 65 or above the risk rises up to 62 extra cases of PE per 10 000.
Of 13 studies included in a recent systematic review and meta-analysis, only three analysed PE outcomes. The present study, therefore, significantly expands previous knowledge on this compelling issue. We showed that adding these results to those of the three observational studies  substantially changed the overall estimate, which now indicates that AP exposure increases of more than three times the risk of PE, with a narrower and statistically significant confidence interval.
Strengths of the study include the following. First, it was based on a large and representative population cohort, extracted from an administrative database covering a population of almost ten million inhabitants, avoiding bias from selection, non-response, or poor recall. This administrative database has been shown to have good levels of accuracy and completeness in recording prescriptions and clinical diagnoses [8,9,12,13]. In this system underreporting is unlikely, as data are used for reimbursement reasons.
Second, unlike previous studies that used composite outcome measures that include diagnostic codes related to venous thromboembolism and PE, with the advantage of increasing statistical power, we selected PE outcomes only, based on clinical grounds, as we reasoned that this is a major healthcare problem, causing over 300,000 deaths annually in Europe alone .
A third strength is that not only we adjusted for many potential confounding factors, but also we a priori excluded individuals with characteristics that might be associated with the outcome of interest, namely those with neoplasm, hip fractures and PE or venous thromboembolism in the year preceding the index date. We also excluded controls with any prescriptions for warfarin before the index date, as this could be treatment for a previous venous thromboembolism. According to Parker and colleagues, prescriptions for warfarin in the six weeks before diagnosis of PE did not trigger exclusion for cases, as they could plausibly be treatment for the index event itself . Cases in which there was any use of warfarin earlier than this six week period, however, were excluded.
Finally, we had access to full prescription details for all AP drugs, including drug name and formulation, dose instructions, and purchase date. This enabled us to look in detail at characteristics of the drug exposure in relation to risk.
Unlike other studies, we were not able to adjust our estimates for some confounding factors such as smoking history, alcohol consumption, and body mass index, which could play a role in the development of PE outcomes [1-4]. However, in order to mitigate this limitation, we employed a nested design by incorporating a case–control approach within an established cohort of individuals who were exposed to AP drugs. This design overcomes some of the disadvantages associated with case–control studies while incorporating some of the advantages of a cohort study. In particular, confounding by indication is less likely, as all subjects were exposed to the variable of interest .
Similarly to other studies, we failed to adjust for some other potential confounding variables, including for example physical activity and inactivity which are well-known risk factor of PE . Additionally, no adjustment was possible for medicines not reimbursed by the National health System, such as for example oral contraceptives, as these are not stored into the administrative databases. Another study limitation is that we had no possibility to access psychiatric diagnoses, which might also be associated with PE outcomes, as recently reported by Parker and colleagues . We were able to include, as a confounding factor, psychiatric hospitalization, but we acknowledge that this may only be considered a proxy measure of psychiatric diagnoses.
As with all analyses of prescription databases, the lack of data on whether patients eventually took the prescribed agents should be highlighted, since a relevant proportion of the medicines prescribed for people with chronic conditions are not taken .
The biological mechanisms explaining the relationship between AP drugs and PE are unknown . It has been suggested that all conditions associated with inactivity or immobilization may be involved. Patients might be immobilized due to the sedative properties of most AP drugs, and AP drugs may also cause weight gain which can increase inactivity and immobilization. Other biological mechanisms include raised levels of antiphospholipid antibodies and hyperhomocysteinemia, which are factors associated with PE outcomes, and might also be associated with AP exposure .
Our findings add to the available body of evidence on the risk of PE in individuals exposed to AP drugs. The three studies available on this compelling issue are highly heterogeneous, as two found an increased risk of more than ten times, while one found a minimal increase in risk, yielding a very imprecise and statistically non significant overall estimate. Adding our study to the forest plot increased the precision of the estimate, and demonstrates that exposure to AP drugs significantly increases the risk of developing PE.
Clinical interpretation of this increased risk is difficult. The main finding of around 16 extra cases of PE per 10 000 patients prescribed AP drugs, or, in other words, the finding that a doctor would need to prescribe AP drugs to around 600 individuals to see one extra case of PE, would suggest that this risk may not be a major consideration when deciding whether AP drugs should be used. However, in specific populations, such as those aged 65 or above and women, the risk is less irrelevant, also taking into consideration that PE may be a life-threatening condition.
Additionally, very often a compelling decision is not whether AP drugs should be used but, rather, which one is less harmful. We found that both first- and second-generation AP drugs carry an increased risk of PE, but there might be differences among individual agents. We recognize that this analysis is statistically less powerful, so that for AP drugs that were not associated with an increased risk we cannot exclude a type-II error, that is failure to detect a difference, if a difference existed.
In conclusion, these data suggest that the expanding use of AP drugs to treat conditions other than approved indications should require careful consideration of the potential harmful consequences of AP drugs, including a small but significant increased risk of PE events. Patients need information on the balance of risks and benefits of these drugs before they start treatment. When treatment is needed for approved indications, clinicians should be aware that AP drugs should cautiously be used, especially among patients at high risk of PE, for example by using minimal effective doses and avoiding polypharmacy, as we observed that dose and use of two or more AP drugs substantially increased the risk. Clinicians should additionally be aware that, as with all other adverse events associated with AP drugs , not all drugs are equally harmful.
- Kyrle PA, Eichinger S. Deep vein thrombosis. Lancet. 2005;365:1163–74.View ArticlePubMedGoogle Scholar
- Cohen AT, Dobromirski M, Gurwith MM. Managing pulmonary embolism from presentation to extended treatment. Thromb Res. 2014;133:139–48.View ArticlePubMedGoogle Scholar
- Zamagni E, Brioli A, Tacchetti P, Zannetti B, Pantani L, Cavo M. Multiple myeloma, venous thromboembolism, and treatment-related risk of thrombosis. Semin Thromb Hemost. 2011;37:209–19.View ArticlePubMedGoogle Scholar
- Walter RJ, Moores LK, Jimenez D. Pulmonary embolism: current and new treatment options. Curr Med Res Opin. 2014;30:1975–89.View ArticlePubMedGoogle Scholar
- Walker AM, Lanza LL, Arellano F, Rothman KJ. Mortality in current and former users of clozapine. Epidemiology. 1997;8:671–7.View ArticlePubMedGoogle Scholar
- Jonsson AK, Spigset O, Hagg S. Venous thromboembolism in recipients of antipsychotics: incidence, mechanisms and management. CNS Drugs. 2012;26:649–62.View ArticlePubMedGoogle Scholar
- Barbui C, Conti V, Cipriani A. Antipsychotic drug exposure and risk of venous thromboembolism: a systematic review and meta-analysis of observational studies. Drug Saf. 2014;37:79–90.View ArticlePubMedGoogle Scholar
- Conti V, Lora A, Cipriani A, Fortino I, Merlino L, Barbui C. Persistence with pharmacological treatment in the specialist mental healthcare of patients with severe mental disorders. Eur J Clin Pharmacol. 2012;68:1647–55.View ArticlePubMedGoogle Scholar
- Lora A, Barbato A, Cerati G, Erlicher A, Percudani M. The mental health system in Lombardy, Italy: access to services and patterns of care. Soc Psychiatry Psychiatr Epidemiol. 2011;47:447–54.View ArticlePubMedGoogle Scholar
- Parker C, Coupland C, Hippisley-Cox J. Antipsychotic drugs and risk of venous thromboembolism: nested case–control study. BMJ. 2010;341:c4245.View ArticlePubMedGoogle Scholar
- WHO Collaborating Centre for Drug Statistic Methodology. Guidelines for ATC Classification and DDD Assignment. WHO Collaborating Centre for Drug Statistic Methodology edn. Oslo: WHO; 2003.Google Scholar
- Berta P, Seghieri C, Vittadini G. Comparing health outcomes among hospitals: the experience of the Lombardy Region. Health Care Manag Sci. 2013;16:245–57.View ArticlePubMedGoogle Scholar
- Barbarito F, Pinciroli F, Barone A, Pizzo F, Ranza R, Mason J, et al. Implementing the lifelong personal health record in a regionalised health information system: The case of Lombardy. Comput Biol Med: Italy; 2013.Google Scholar
- Sedgwick P. Nested case–control studies. BMJ. 2010;340:c2582.View ArticlePubMedGoogle Scholar
- Curtis D. Antipsychotics and venous thrombosis. Confounding factors may account for the association. BMJ. 2010;341:c5628.View ArticlePubMedGoogle Scholar
- Jones G. Prescribing and taking medicines. BMJ. 2003;327:819.View ArticlePubMedPubMed CentralGoogle Scholar
- Hagg S, Jonsson AK, Spigset O. Risk of venous thromboembolism due to antipsychotic drug therapy. Expert Opin Drug Saf. 2009;8:537–47.View ArticlePubMedGoogle Scholar
- Leucht S, Cipriani A, Spineli L, Mavridis D, Orey D, Richter F, et al. Comparative efficacy and tolerability of 15 antipsychotic drugs in schizophrenia: a multiple-treatments meta-analysis. Lancet. 2013;382:951–62.View ArticlePubMedGoogle 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/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.