Abstract
Scandinavian electronic health-care registers provide a unique setting to investigate potential unidentified side effects of drugs. We analysed the association between prescription drugs dispensed in Norway and Sweden and the short-term risk of developing pulmonary embolism. A total of 12,104 pulmonary embolism cases were identified from patient- and cause-of-death registries in Norway (2004–2014) and 36,088 in Sweden (2005–2014). A case-crossover design was used to compare individual drugs dispensed 1–30 days before the date of pulmonary embolism diagnosis with dispensation in a 61–90 day time-window, while controlling for the receipt of other drugs. A BOLASSO approach was used to select drugs that were associated with short-term risk of pulmonary embolism. Thirty-eight drugs were associated with pulmonary embolism in the combined analysis of the Norwegian and Swedish data. Drugs associated with increased risk of pulmonary embolism included certain proton-pump inhibitors, antibiotics, antithrombotics, vasodilators, furosemide, anti-varicose medications, corticosteroids, immunostimulants (pegfilgrastim), opioids, analgesics, anxiolytics, antidepressants, antiprotozoals, and drugs for cough and colds. Mineral supplements, hydrochlorothiazide and potassium-sparing agents, beta-blockers, angiotensin 2 receptor blockers, statins, and methotrexate were associated with lower risk. Most associations persisted, and several additional drugs were associated, with pulmonary embolism when using a longer time window of 90 days instead of 30 days. These results provide exploratory, pharmacopeia-wide evidence of medications that may increase or decrease the risk of pulmonary embolism. Some of these findings were expected based on the drugs' indications, while others are novel and require further study as potentially modifiable precipitants of pulmonary embolism.
Similar content being viewed by others
Introduction
Pulmonary embolism (PE) arises when blood clots, often originating in deep veins, obstruct pulmonary arteries1. The development of PE encompasses an intricate interplay of pathophysiological, inflammatory, and chemical processes. These processes involve factors such as stasis, endothelial injury, hypercoagulability, inflammatory responses, damage to the endothelium, insults to the coagulation cascade (involving fibrin, thrombin, and tissue factor) as well as additional factors such as immobility and genetic predispositions1. The incidence of pulmonary embolism has increased globally over the last decades, particularly in high-income countries, while pulmonary embolism-related in-hospital death rate and age-standardized mortality from pulmonary embolism have decreased or plateaued2,3,4,5. Improved survival in patients with conditions predisposing to pulmonary embolism, such as cancer, chronic obstructive pulmonary disease, and autoimmune diseases, may contribute to the increased incidence rates, while a greater proportion of low-risk cases being diagnosed and improvements in the management of pulmonary embolism may have contributed to the reduced mortality rates 2. Several risk factors have been identified for pulmonary embolism, including age6, sex6, smoking7, obesity8,9, hypertension7, and physical inactivity10. In addition, surgery11, hospitalization, or trauma that leads to extended immobilization12, cancer12,13, pregnancy14, kidney disease15, heart failure12, or previous history of deep venous thromboembolism or pulmonary embolism are established risk factors. Some medications including oral contraceptives16, hormone replacement therapy17, tamoxifen18,19,20, bisphosphonates21, coagulants22, antidepressants23, and antipsychotics24 have been associated with increased risk, while statins have been associated with reduced risk25.
In the US, approximately 90% of adults aged 65 years and older use at least one prescription medication26. Side effects from use of drugs are a major public health concern and have been estimated to cost around $30 billion in the US27. Randomized trials used to test drugs are costly and therefore generally just large enough to detect an effect on the primary outcome of interest, which may be a physiological effect (e.g., blood pressure) rather than a hard clinical outcome. Less common, but important side effects may therefore be missed in many trials. In addition, pre-approval trials have tended to include fewer women, especially in their reproductive age (although some improvement has been observed in later years28), patients with comorbidities, children, and older adults, which may limit the generalizability of the findings from such studies to the general population29. Thus, systematic monitoring of the pharmaceutical effects of all approved drugs in clinical practice could advance clinical care and public health.
Our research group has previously conducted systematic examinations of all potential associations between prescribed drugs and short-term risk of acute myocardial infarction30 and ischemic stroke31, and others have used this approach to identify drugs that may improve or impair COVID-19 prognosis32,33. In the current study, we extended this approach to examine the association between prescribed medications and short-term risk of pulmonary embolism using comprehensive nationwide data in two countries.
Methods
Study design
A case-crossover design was used and methods have been described in detail previously30,31. This design included all cases of pulmonary embolism and applied self-matching by comparing exposure before disease onset with disease-free time in the past as a comparison control. A major advantage of the case-crossover design is that person-specific characteristics that are stable over the typically short-time periods such studies are conducted over (e.g. age, sex, lifestyle, chronic conditions) do not confound the observed associations.
Dispensed medications
We assessed the risk of pulmonary embolism associated with every drug dispensed to patients that had a first-time pulmonary embolism within the study period. Data on dispensed medications prior to the diagnosis or death were extracted from the nation-wide registries of dispensed drugs in both Norway and Sweden. The Norwegian Prescription Database was established on January 1st 2004, and all Norwegian pharmacies are required to supply information on drug prescriptions including type, dose and date of dispensation34. Sweden established a similar registry on July 1st 2005, the Swedish National Prescribed Drug Register35. National personal identifiers attached to these data were used to link the information on drug prescriptions to other health-related registers existing in these countries. The databases do not include information on drugs purchased over-the-counter or drugs administered to institutionalized patients in nursing homes or hospitals. In Norway, it was possible to exclude participants who, at the time of their pulmonary embolism, were institutionalized and for whom registration of dispensed medications was not available. In Sweden, this information was not available, and therefore we included only those patients to whom at least one drug was dispensed during the year preceding the occurrence of pulmonary embolism. In Sweden, prescriptions are in general valid for one year after date of being issued36, and for chronic medication are refilled every third month. A similar system is in place in Norway.
Outcome assessment
The Norwegian Patient Register (established in 2008)37, the Swedish National Patient Register (established in 1964)38, and the cause of death registries in Norway (established in 1951)39 and Sweden (established in 1952)40 were used to identify cases of pulmonary embolism. A total of 12,104 pulmonary embolism cases were identified from patient- and cause-of-death registries in Norway (2004–2014) and 36,088 in Sweden (2005–2014). The quality of the information in the Norwegian and Swedish patient registers has been shown to be very high for other cardiovascular outcomes38,41,42, and acceptable for acute pulmonary embolism in Sweden43. In Norway, all patients registered with either a primary hospital discharge diagnosis or underlying cause of death of ICD-10 I26 from 1 January 2008 to 31 December 2014 were included, and in Sweden, the corresponding dates for hospital diagnosis and cause of death were between 1 November 2005 and 31 December 2014. Only the first registered episode of pulmonary embolism was included in the analysis for each participant.
Statistical analysis
For each patient, the occurrence of drug dispensing within 1–30 days before the date of pulmonary embolism occurrence (case period) was compared to a time window of 61–90 days before pulmonary embolism diagnosis (control period) for each drug individually. A 30 day wash-out period between the case- and control-periods was used to minimize the carryover effects of drugs. We calculated odds ratios together with 95% confidence intervals comparing the odds of drug dispensed in the case period to that in the control period using conditional logistic regression. The statistical models were adjusted for every other drug. Adjustments for other drugs were made at the fifth ATC level—so each drug investigated was adjusted for every other drug investigated. Additional analyses were conducted to examine the robustness of the findings where we extended the case-, control- and wash-out periods from 30 to 90 days (case period = 1–90 days before event occurrence, control period 181–270 days before event occurrence) and repeated all analyses. The 30 days44 and 90 days44,45,46 exposure periods were based on prior publications44,45,46.
We assessed the association between all dispensed medications and pulmonary embolism risk. Because we aimed to estimate the most likely effect size for drugs with true associations while accounting for simultaneous prescriptions, we did not use Bonferroni correction or similar conventional methods to address the issue of multiple comparisons, as they may fail to estimate the size of these associations correctly47. Instead, we applied a version of the least absolute shrinkage and selection operator (LASSO) regression analysis47,48,49,50,51 called bootstrap-enhanced least absolute shrinkage operator (BOLASSO)52. Several bootstrap samples are drawn from the dataset, where each bootstrap sample is generated by sampling N pairs (N is the total number of drugs in the dataset) with replacement. We drew 1000 bootstrap samples in this analysis. Of note, confidence intervals generated via the BOLASSO approach are not optimal, because each bootstrap sample is estimated on different penalty parameters. We include confidence intervals for ease of interpretation. As a result, drugs selected by this approach may include one (i.e., the null) within the confidence intervals. In BOLASSO, we obtain multivariable-adjusted estimates as the effect of each selected drug is controlled for the effects of all other selected drugs. In Supplementary Material, Online Appendix A, we present in detail the background of the method and how we implemented BOLASSO in conditional logistic regression models for case-crossover data.
Separate analyses were conducted for the Norwegian and Swedish data and we present both country-specific and combined estimates, using meta-analysis to combine the country-specific estimates. Drugs selected by BOLASSO with risk estimates in the same direction for both countries were considered a common hit and were included in the final analysis. A fixed-effects model was used to calculate the combined estimate53. All statistical analyses were conducted using R (version 3.2.3; R foundation for Statistical Computing, Vienna, Austria) and Stata/IC 16 (Stata Corp, College Station, Texas, USA).
Ethical approval
The studies were approved by the Regional Committee for Medical and Health Research Ethics (REC) in Central Norway and Regional Ethical Review Board in Sweden. In addition, Norwegian data was also approved by the Norwegian Data Protection Authority. All methods were performed in accordance with the relevant guidelines and regulations by the respective ethical committees from both Norway and Sweden. Exemption from the requirement of obtaining informed consent from the registered individuals was given by REC. In all data files the personal identifier was replaced with a study allocation number as part of the data preparation by the data providers.
Results
A total of 48,192 pulmonary embolism cases were included across the two countries, 36,088 from Sweden (33,678 identified via the patient register and 2410 via the cause of death register) and 12,104 from Norway (11,947 and 157 identified via the corresponding Norwegian registers). Characteristics of these individuals are presented in Table 1.
A total of 1100 and 1260 distinct pharmaceutical drugs were dispensed among patients who experienced a subsequent pulmonary embolism in Norway and Sweden, respectively. Out of these, 773 unique drugs were dispensed in either the case- or control-period, in Norway and 1091 in Sweden. Of these, BOLASSO selected 117 unique drugs in Norway (Supplementary Table 1) and 191 in Sweden (Supplementary Table 2). Finally, a total of 59 drugs were selected in common from both countries (Fig. 1). Table 2 presents the country-specific and combined estimates of these dually-selected drugs.
Case-crossover analysis of dispensed prescription medication use and risk of pulmonary embolism. The above plot illustrates (A) unique drug types which were selected in Norway, (B) unique drug types which were selected in Sweden, and (C) 59 drugs which were common hits from both the countries. Y-axis displays relative risk on the log scale, X-axis displays all the prescribed drugs studied grouped by the anatomical therapeutic chemical (ATC) classification.
Cardiovascular drugs
Several cardiovascular drugs were associated with higher risk of pulmonary embolism, some of which could be due to their indication of use. The antithrombotic agents dalteparin, enoxaparin, and clopidogrel were associated with increased risk, while the angiotensin 2 receptor blocker ‘candesartan and diuretics’ was associated with reduced risk. Diuretics yielded mixed results, with lower risk for hydrochlorothiazide and potassium-sparing agents, while furosemide was associated with higher risk. The vasodilator glyceryl trinitrate was associated with increased risk. An inverse association was observed for the beta-blocking agent carvedilol.
Antibiotics
Several antibiotics examined were associated with increased risk. This includes nystatin, doxycycline, amoxicillin, pivmecillinam, phenoxymethylpenicillin, trimethoprim, sulfamethoxazole and trimethoprim combined, erythromycin, ciprofloxacin, and nitrofurantoin.
Opioids
Several opioids were associated with higher pulmonary embolism risk, including morphine, oxycodone, oxycodone and naloxone combined, codeine combinations excluding psycholeptics, fentanyl, buprenorphine, and tramadol.
Antineoplastic and immunomodulating agents
A positive association was observed for pegfilgrastim, while lower risk was observed for methotrexate.
Other medications
Analgesics (paracetamol), antiprotozoals (metronidazole), anti-varicose therapy (organo-heparinoid), anxiolytic (oxazepam, but not diazepam), proton pump inhibitors (pantaprazole, but not omeprazole) and corticosteroids for systemic use (prednisolone), but not corticosteroids for dermatological use (clobetasol) were associated with increased risk of pulmonary embolism. Positive associations were observed for cough and cold preparations, such as acetylcysteine and opium derivatives and expectorants. An inverse association was observed for calcium supplements combined with vitamin D or other drugs/supplements.
Additional analyses with a longer exposure window
In additional analyses, when we extended the case-, control and wash-out periods from 30 to 90 days, many of the results remained similar, while others were either stronger or weaker than in the primary analysis (Supplementary Table 3). Additional positive associations were observed for ACE inhibitors (enalapril, ramipril), adrenergics/drugs for obstructive airway diseases (salbutamol, terbutaline, indacaterol, formoterol and budenoside), angiotensin II receptor blockers (losartan), antianemics (ferrous sulfate), antibiotics (azithromycin), antidepressants (sertraline, mianserin), antiemetics and antinauseants (aprepitant), antifibrinolytics (tranexamic acid), anti-inflammatory and anti-rheumatic agents (etoricoxib), antimyotics (fluconazole), antineoplastic and immunomodulating agents (capecitabine, erlotinib, tamoxifen), antipropulsives (loperamide), antipsychotics (haloperidol, risperidone), anxiolytics (diazepam), cardiac glycosides (digoxin), corticosteroids for systemic use (dexamethasone), drugs for treating constipation (sodium picosulfate), estrogens (estriol), hormonal contraceptives (levonorgestrel and ethinylestradiol, drospirenone and ethinylestradiol), hypnotics (zopiclone, zolpidem), intravaginal ring (with progestogen and estrogen), opioids (tramadol), and proton pump inhibitors (esomeprazole) (Supplementary Table 3). In addition, the positive associations between prednisolone and pulmonary embolism was substantially strengthened (Supplementary Table 3). Additional inverse associations were observed between anti-inflammatory and anti-rheumatic agents (glucosamine), antithrombotics (warfarin), calcium channel blockers (felodipine), and drugs for urinary frequency and incontinence (tolterodine) and pulmonary embolism (Supplementary Table 3). Supplementary Tables 4, 5 shows the unique drugs selected from each country using the 90 day time window, but also includes hits that were not common.
Discussion
This study systematically investigated the associations between all drugs requiring a prescription and short-term risk of developing pulmonary embolism. Using nationwide registry data from Norway and Sweden, and employing an exposure-wide approach, we identified 38 drugs that were either associated with higher (31 drugs) or lower (seven drugs) short-term risk of pulmonary embolism. In additional analyses, using a longer exposure period, which may be more appropriate for some drugs, 66 drugs were associated with higher risk and seven with lower risk.
Antibiotics tended to be most consistently associated with higher risk of pulmonary embolism, and we observed higher risk across several types of antibiotics examined. However, these associations are likely to be confounded by underlying infection, which is a strong risk factor for pulmonary embolism11,54,55, as it might lead to increased inflammation and triggering of platelets, resulting in fibrin deposition and thrombus formation56. We also cannot exclude the possibility of reverse causation, where premonitory symptoms of pulmonary embolism like dyspnea are empirically presumed to reflect the more common occurrence of respiratory infection.
Prescription of various opioids was also quite consistently associated with higher risk of pulmonary embolism. This may reflect an increased risk of pulmonary embolism during and after surgery13,57. Almost 25% of all cases of venous thromboembolism occur during or shortly after surgery12. Similarly, the positive associations observed for certain antithrombotics also likely reflects confounding by indication, as these drugs are used to reduce the risk of pulmonary embolism during or after surgery.
The positive association between the antidepressant mirtazapine and pulmonary embolism is consistent with the Million Women’s Study, which reported a positive association between antidepressant use and risk of venous thromboembolism with pulmonary embolism23. However, a case-control study reported no association58.
We observed a positive association between drugs used to treat cough and cold (e.g., acetylcysteine, and opium derivatives and expectorants) and obstructive airway disease (e.g., salbutamol, formoterol and budesonide) and risk of pulmonary embolism. An acute infection or exacerbation could potentially explain these associations.
Blood pressure-lowering drugs had differing associations with pulmonary embolism risk. While inverse associations were observed for carvedilol, hydrochlorothiazide, and for candesartan in combination with diuretics, furosemide was associated with increased risk. Some large cohort studies have suggested an inverse association between blood pressure and pulmonary embolism risk59,60. It is possible that the increased pulmonary embolism risk with use of furosemide could be due to activation of the renin–angiotensin–aldosterone system61, which has a pro-thrombotic effect62. This, on the other hand, could explain why angiotensin-2 receptor blockers reduce the risk of pulmonary embolism, especially in a possible subpopulation that requires a combination of angiotensin-2 receptor blockers and diuretics and may have higher risk of atherosclerotic disease63. Beta blockers, and mainly unselected beta blockers, such as carvedilol, can reduce increased sympathetic activity with their prothrombotic effect and thereby reduce pulmonary embolism risk64. However, other potential explanations may be reverse causation, whereby premonitory dyspnea is presumed to reflect congestive heart failure, confounding by indication, perhaps by kidney disease65, or acute dehydration by furosemide and increased risk of pulmonary embolism as part of the classical Virchow’s triad66.
The inverse association observed between dispensed statins (simvastatin and pravastatin) and pulmonary embolism is consistent with results from a cohort study25 showing reduced incidence from pulmonary embolism with use of statins. The results are also partly consistent with meta-analyses of randomized controlled trials, showing reduced incidence from venous thromboembolism with use of statins67,68. However, in these meta-analyses, only one trial specifically focused on pulmonary embolism and reported an imprecise effect estimate (HR = 0.77, 0.41–1.46). The association between serum cholesterol level with risk of pulmonary embolism is not clear6,69,70,71, and thus it is possible that the observed benefit with statins is due to other non-specific effects unrelated to lipid-lowering, such as anti-inflammatory effects72.
There was a weak inverse association observed between calcium supplementation combined with vitamin D or other drugs and risk of pulmonary embolism (HR = 0.89, 0.81–0.97). The Women’s Health Initiative randomized controlled trial reported a similar, but less precise, association between calcium and vitamin D supplementation in relation to pulmonary embolism risk (HR = 0.92, 95% CI 0.73–1.16)73.
The strength of this study includes the large-scale registry-based databases from Norway and Sweden and the comprehensive assessment of a wide range of drugs in relation to short-term risk of pulmonary embolism. These registers are of high quality and near complete. Case-crossover studies are vulnerable to misclassification74, however, in our study, we could minimize some of the sources of bias due to misclassification. Both countries have a tax-financed universal healthcare system that is accessible to all residents and reporting to the health registers is obligatory. Recall bias with regard to the exposure was eliminated by relying on recorded data. Moreover, the outcome, i.e. pulmonary embolism, can be identified with high degree of accuracy due to the modern diagnostic methods and the universal healthcare in Scandinavia43, reducing the potential for selection bias. Use of nationwide databases also supports the generalisability of the findings. The large sample size provided sufficient statistical power to detect even modest associations.
Several of the observed findings likely reflect confounding by indication, and we were not able to separate the effect of a drug and its indication in the current study. It is therefore unclear whether the observed associations were due to the drugs themselves or due to the underlying or prodromal conditions the drugs were meant to treat. However, this type of confounding may be less likely to explain the results when observing inverse associations or when different medications within the same group of medications are differentially associated with pulmonary embolism. In such cases, a direct effect of a drug on pulmonary embolism is more likely. Another limitation of the study is that we did not have detailed information on lifestyle-related risk factors or comorbidities that potentially may have acted as effect modifiers. Most previous studies on drugs and pulmonary embolism were observational cohort and case-control studies23,25,58,75. We did not identify prior case-crossover studies on medications and pulmonary embolism risk. Thus, our comparison with previous literature is also limited by differences in design and potential latency periods between drug prescription and the development of pulmonary embolism. The current case-crossover design is best suitable to identify acute or triggering effects of drugs on outcomes with a sudden onset, but may misestimate the effects of long-term drug use on risk(76). It was beyond this project to assess potential drug interactions, but this may be a point for further investigation in future studies. The databases did not include information on over-the-counter drugs or drug prescriptions in hospitals or nursing homes, however, in Norway institutionalized subjects without registration of dispensed medications were excluded, while in Sweden only subjects with at least one drug dispensed in the year preceding the occurrence of pulmonary embolism were included. Lastly, the prescription databases only contain information on the date of dispensing, and not information on the date of actual administration of drugs or whether the drug was taken, but this likely resulted in non-differential misclassification and bias toward the null.
In conclusion, we found both positive and inverse associations between a wide range of dispensed prescription drugs and short-term risk of pulmonary embolism. Further studies are needed to explore these findings in more detail. This method presents an opportunity for systematic and thorough monitoring of a wide range of drugs that could trigger or prevent the onset of serious health conditions, including pulmonary embolism.
Data availability
The data used in the current analysis are available from the Norwegian Patient Register, Norwegian Prescription Database, Norwegian Cause of Death Registry, and Swedish National Patient Register, Swedish National Prescribed Drug Register, Swedish National Cause of Death Register, but restrictions apply to the availability of these data, which were used under license for the current study, and are therefore not publicly available. Requests for data access should be addressed to Abhijit Sen. Email: abhijit.sen@ntnu.no.
References
Turetz, M., Sideris, A. T., Friedman, O. A., Triphathi, N. & Horowitz, J. M. Epidemiology, pathophysiology, and natural history of pulmonary embolism. Semin. Intervent. Radiol. 35, 92–98 (2018).
Barco, S. et al. Global reporting of pulmonary embolism-related deaths in the World Health Organization mortality database: Vital registration data from 123 countries. Res. Pract. Thromb. Haemost. 5, e12520 (2021).
Dentali, F. et al. Time trends and case fatality rate of in-hospital treated pulmonary embolism during 11 years of observation in Northwestern Italy. Thromb. Haemost. 115, 399–405 (2016).
Arshad, N., Isaksen, T., Hansen, J. B. & Brækkan, S. K. Time trends in incidence rates of venous thromboembolism in a large cohort recruited from the general population. Eur. J. Epidemiol. 32, 299–305 (2017).
Raptis, D. G., Gourgoulianis, K. I., Daniil, Z. & Malli, F. Time trends for pulmonary embolism incidence in Greece. Thromb. J. 18, 1 (2020).
Tsai, A. W. et al. Cardiovascular risk factors and venous thromboembolism incidence: The longitudinal investigation of thromboembolism etiology. Arch. Intern. Med. 162, 1182–1189 (2002).
Goldhaber, S. Z. et al. A prospective study of risk factors for pulmonary embolism in women. JAMA 277, 642–645 (1997).
Kabrhel, C., Varraso, R., Goldhaber, S. Z., Rimm, E. B. & Camargo, C. A. Prospective study of BMI and the risk of pulmonary embolism in women. Obesity 17, 2040–2046 (2009).
Rahmani, J. et al. Relationship between body mass index, risk of venous thromboembolism and pulmonary embolism: A systematic review and dose-response meta-analysis of cohort studies among four million participants. Thromb. Res. 192, 64–72 (2020).
Kabrhel, C., Varraso, R., Goldhaber, S. Z., Rimm, E. & Camargo, C. A. Jr. Physical inactivity and idiopathic pulmonary embolism in women: Prospective study. BMJ 343, d3867 (2011).
Masrouha, K. Z., Musallam, K. M., Rosendaal, F. R., Hoballah, J. J. & Jamali, F. R. Preoperative pneumonia and postoperative venous thrombosis: A cohort study of 427,656 patients undergoing major general surgery. World J. Surg. 40, 1288–1294 (2016).
Heit, J. A. et al. Relative impact of risk factors for deep vein thrombosis and pulmonary embolism: A population-based study. Arch. Intern. Med. 162, 1245–1248 (2002).
Caron, A. et al. Risk of pulmonary embolism more than 6 weeks after surgery among cancer-free middle-aged patients. JAMA Surg. 154, 1126–1132 (2019).
Dado, C. D., Levinson, A. T. & Bourjeily, G. Pregnancy and pulmonary embolism. Clin. Chest Med. 39, 525–537 (2018).
Kuo, T. H., Li, H. Y. & Lin, S. H. Acute kidney injury and risk of deep vein thrombosis and pulmonary embolism in Taiwan: A nationwide retrospective cohort study. Thromb. Res. 151, 29–35 (2017).
Pearce, H. M., Layton, D., Wilton, L. V. & Shakir, S. A. Deep vein thrombosis and pulmonary embolism reported in the prescription event monitoring study of Yasmin. Br. J. Clin. Pharmacol. 60, 98–102 (2005).
Barlow, D. H. HRT and the risk of deep vein thrombosis. Int. J. Gynaecol. Obstet. 59(1), S29–S33 (1997).
Hernandez, R. K., Sorensen, H. T., Pedersen, L., Jacobsen, J. & Lash, T. L. Tamoxifen treatment and risk of deep venous thrombosis and pulmonary embolism: A Danish population-based cohort study. Cancer 115, 4442–4449 (2009).
Lin, H. F. et al. Correlation of the tamoxifen use with the increased risk of deep vein thrombosis and pulmonary embolism in elderly women with breast cancer: A case-control study. Medicine 97, e12842 (2018).
Sun, R., Chu, Y., Gao, Y., Cheng, W. & Gao, S. Efficacy and safety of endocrine therapy for breast-cancer prevention in high-risk premenopausal or postmenopausal women: A Bayesian network meta-analysis of nine randomized controlled trials. Menopause 28, 589–600 (2021).
Vestergaard, P., Schwartz, K., Pinholt, E. M., Rejnmark, L. & Mosekilde, L. Use of bisphosphonates and raloxifene and risk of deep venous thromboembolism and pulmonary embolism. Osteoporos. Int. 21, 1591–1597 (2010).
Myers, S. P. et al. Tranexamic acid administration is associated with an increased risk of posttraumatic venous thromboembolism. J. Trauma Acute Care Surg. 86, 20–27 (2019).
Parkin, L. et al. Antidepressants, depression, and venous thromboembolism risk: Large prospective study of UK women. J. Am. Heart Assoc. 6, e005316 (2017).
Arasteh, O. et al. Antipsychotic drugs and risk of developing venous thromboembolism and pulmonary embolism: A systematic review and meta-analysis. Curr. Vasc. Pharmacol. 18, 632–643 (2020).
Lassila, R., Jula, A., Pitkaniemi, J. & Haukka, J. The association of statin use with reduced incidence of venous thromboembolism: A population-based cohort study. BMJ Open 4, e005862 (2014).
Charlesworth, C. J., Smit, E., Lee, D. S., Alramadhan, F. & Odden, M. C. Polypharmacy among adults aged 65 years and older in the United States: 1988–2010. J. Gerontol. A Biol. Sci. Med. Sci. 70, 989–995 (2015).
Sultana, J., Cutroneo, P. & Trifiro, G. Clinical and economic burden of adverse drug reactions. J. Pharmacol. Pharmacother. 4, S73–S77 (2013).
Dekker, M. J. H. J. et al. Sex proportionality in pre-clinical and clinical trials: An evaluation of 22 marketing authorization application Dossiers submitted to the European medicines agency. Front. Med. 8, 643028 (2021).
Van Spall, H. G., Toren, A., Kiss, A. & Fowler, R. A. Eligibility criteria of randomized controlled trials published in high-impact general medical journals: A systematic sampling review. JAMA 297, 1233–1240 (2007).
Sen, A. et al. Systematic assessment of prescribed medications and short-term risk of myocardial infarction—A pharmacopeia-wide association study from Norway and Sweden. Sci. Rep. 9, 8257 (2019).
Janszky, I. et al. Assessing short-term risk of ischemic stroke in relation to all prescribed medications. Sci. Rep. 11, 21673 (2021).
Lerner, I. et al. Mining electronic health records for drugs associated with 28 day mortality in COVID-19: Pharmacopoeia-wide association study (PharmWAS). JMIR Med. Inform. 10, e35190 (2022).
MacFadden, D. R. et al. Screening large population health databases for potential coronavirus disease 2019 therapeutics: A pharmacopeia-wide association study of commonly prescribed medications. Open Forum Infect. Dis. 9, ofac156 (2022).
Furu, K. Establishment of the nationwide Norwegian prescription database (NorPD)—New opportunities for research in pharmacoepidemiology in Norway. Nor. J. Epidemiol. 18, 129–136 (2008).
Wettermark, B. et al. The new Swedish prescribed drug register–opportunities for pharmacoepidemiological research and experience from the first six months. Pharmacoepidemiol. Drug Saf. 16, 726–735 (2007).
Bakken, I. J., Ariansen, A. M. S., Knudsen, G. P., Johansen, K. I. & Vollset, S. E. The Norwegian patient registry and the Norwegian registry for primary health care: Research potential of two nationwide health-care registries. Scand. J. Public Health 48, 49–55 (2020).
Ludvigsson, J. F. et al. External review and validation of the Swedish national inpatient register. BMC Public Health 11, 450 (2011).
Pedersen, A. G. & Ellingsen, C. L. Data quality in the causes of death registry. Tidsskr. Nor. Laegeforen. 135, 768–770 (2015).
Brooke, H. L. et al. The Swedish cause of death register. Eur. J. Epidemiol. 32, 765–773 (2017).
Govatsmark, R. E. S. et al. Completeness and correctness of acute myocardial infarction diagnoses in a medical quality register and an administrative health register. Scand. J. Public Health 48, 5–13 (2020).
Varmdal, T. et al. Comparison of the validity of stroke diagnoses in a medical quality register and an administrative health register. Scand. J. Public Health 44, 143–149 (2016).
Andersson, T. et al. Validation of the Swedish national inpatient register for the diagnosis of pulmonary embolism in 2005. Pulm. Circ. 12, e12037 (2022).
Walker, R. F. et al. Association of testosterone therapy with risk of venous thromboembolism among men with and without hypogonadism. JAMA Intern. Med. 180, 190–197 (2020).
Grimnes, G. et al. C-reactive protein and risk of venous thromboembolism: Results from a population-based case-crossover study. Haematologica 103, 1245–1250 (2018).
Grimnes, G., Isaksen, T., Tichelaar, Y. I. G. V., Brækkan, S. K. & Hansen, J. B. Acute infection as a trigger for incident venous thromboembolism: Results from a population-based case-crossover study. Res. Pract. Thromb. Haemost. 2, 85–92 (2018).
Greenland, S. & Robins, J. M. Empirical-Bayes adjustments for multiple comparisons are sometimes useful. Epidemiology 2, 244–251 (1991).
Zhao, P. & Yu, B. On model selection consistency of Lasso. J. Mach. Learn. Res. 7, 2541–2563 (2006).
Nee, M. et al. Prescription medicine use by pedestrians and the risk of injurious road traffic crashes: A case-crossover study. PLoS Med. 14, e1002347 (2017).
Avalos, M. et al. Prescription-drug-related risk in driving: Comparing conventional and lasso shrinkage logistic regressions. Epidemiology 23, 706–712 (2012).
Steenland, K., Bray, I., Greenland, S. & Boffetta, P. Empirical Bayes adjustments for multiple results in hypothesis-generating or surveillance studies. Cancer Epidemiol. Biomark. Prev. 9, 895–903 (2000).
Bach, F. Bolasso: Model consistent Lasso estimation through the bootstrap. In ICML '08: Proceedings of the 25th international conference on Machine learning, 33–40. https://doi.org/10.1145/1390156.1390161 (Ithaca, New York, 2008).
Borenstein, M., Hedges, L. V., Higgins, J. P. & Rothstein, H. R. A basic introduction to fixed-effect and random-effects models for meta-analysis. Res. Synth. Methods 1, 97–111 (2010).
Clayton, T. C., Gaskin, M. & Meade, T. W. Recent respiratory infection and risk of venous thromboembolism: Case-control study through a general practice database. Int. J. Epidemiol. 40, 819–827 (2011).
Chen, Y. G. et al. Association between pneumococcal pneumonia and venous thromboembolism in hospitalized patients: A nationwide population-based study. Respirology 20, 799–804 (2015).
Beristain-Covarrubias, N. et al. Understanding infection-induced thrombosis: Lessons learned from animal models. Front. Immunol. 10, 2569 (2019).
Sweetland, S. et al. Duration and magnitude of the postoperative risk of venous thromboembolism in middle aged women: Prospective cohort study. BMJ 339, b4583 (2009).
Lacut, K. et al. Association between antipsychotic drugs, antidepressant drugs and venous thromboembolism: Results from the EDITH case-control study. Fundam. Clin. Pharmacol. 21, 643–650 (2007).
Gregson, J. et al. Cardiovascular risk factors associated with venous thromboembolism. JAMA Cardiol. 4, 163–173 (2019).
Nazarzadeh, M. et al. Blood pressure and risk of venous thromboembolism: A cohort analysis of 5.5 million UK adults and Mendelian randomization studies. Cardiovasc. Res. 119, 835–842 (2022).
Potter, B. M., Ames, M. K., Hess, A. & Poglitsch, M. Comparison between the effects of torsemide and furosemide on the renin–angiotensin–aldosterone system of normal dogs. J. Vet. Cardiol. 26, 51–62 (2019).
Bekassy, Z., Lopatko Fm, I., Bader, M. & Karpman, D. Crosstalk between the renin–angiotensin, complement and kallikrein–kinin systems in inflammation. Nat. Rev. Immunol. 22, 411–428 (2022).
Chae, Y. K. et al. Inhibition of renin angiotensin axis may be associated with reduced risk of developing venous thromboembolism in patients with atherosclerotic disease. PLoS One 9, e87813 (2014).
De Peuter, O. R. et al. Non-selective vs. selective beta-blocker treatment and the risk of thrombo-embolic events in patients with heart failure. Eur. J. Heart Fail. 13, 220–226 (2011).
Singh, J. et al. Pulmonary embolism in chronic kidney disease and end-stage renal disease hospitalizations: Trends, outcomes, and predictors of mortality in the United States. SAGE Open Med. 9, 20503121211022996 (2021).
Ku, E., Lee, B. J., Wei, J. & Weir, M. R. Hypertension in CKD: Core curriculum 2019. Am. J. Kidney Dis. 74, 120–131 (2019).
Li, L., Zhang, P., Tian, J. H. & Yang, K. Statins for primary prevention of venous thromboembolism. Cochrane Database Syst. Rev. 2014, CD008203 (2014).
Kunutsor, S. K., Seidu, S. & Khunti, K. Statins and primary prevention of venous thromboembolism: A systematic review and meta-analysis. Lancet Haematol. 4, e83–e93 (2017).
Goldhaber, S. Z. et al. Risk factors for pulmonary embolism: The Framingham study. Am. J. Med. 74, 1023–1028 (1983).
Wattanakit, K. et al. Association between cardiovascular disease risk factors and occurrence of venous thromboembolism. Thromb. Haemost. 108, 508–515 (2012).
Hu, M., Li, X. & Yang, Y. Causal associations between cardiovascular risk factors and venous thromboembolism. Semin. Thromb. Hemost. 49, 679–687 (2023).
Rodriguez, A. L. et al. Statins, inflammation and deep vein thrombosis: A systematic review. J. Thromb. Thrombolysis 33, 371–382 (2012).
Blondon, M. et al. The effect of calcium plus vitamin D supplementation on the risk of venous thromboembolism. From the women’s health initiative randomized controlled trial. Thromb. Haemost. 113, 999–1009 (2015).
Mittleman, M. A. & Mostofsky, E. Exchangeability in the case-crossover design. Int. J. Epidemiol. 43, 1645–1655 (2014).
Liu, Y. et al. Current antipsychotic agent use and risk of venous thromboembolism and pulmonary embolism: A systematic review and meta-analysis of observational studies. Ther. Adv. Psychopharmacol. 11, 2045125320982720 (2021).
Hallas, J., Pottegård, A., Wang, S., Schneeweiss, S. & Gagne, J. J. Persistent user bias in case-crossover studies in pharmacoepidemiology. Am. J. Epidemiol. 184, 761–769 (2016).
Funding
This work was supported by Central Norway Regional Authority (project number 46060913).
Author information
Authors and Affiliations
Contributions
AS had full access to the data and takes responsibility for the integrity of the data and the accuracy of the data analysis, and conducted the statistical analysis. IJ, KJM, RL, and PRR conceived and designed the study. IJ, RL, IJB, PRR, and BMB acquired the data. DA, IV, BHL, BMB, LBS, JWH, IJB, PRR, KJM, RL, IJ, and AS interpreted the data. DA drafted the manuscript. DA, IV, BHL, BMB, LBS, JWH, IJB, PRR, KJM, RL, IJ, and AS critically revised the manuscript for important intellectual content. IJB, BMB, and RL provided administrative, technical or logistic support. IJ obtained funding. IJ, KLM, and RL supervised the study. AS and IJ are guarantors of the Norwegian data and RL is the guarantor of the Swedish data. Exemption from the requirement of obtaining informed consent from the registered individuals was given by REC.
Corresponding author
Ethics declarations
Competing interests
Rickard Ljung is employed at the Swedish Medical Products Agency, Uppsala, Sweden. The views expressed in this paper do not necessarily represent the views of the Government agency. The remaining authors have nothing to disclose. The interpretation and reporting of these data are the sole responsibility of the authors, and no endorsement by the Department of Health Registries is intended nor should be inferred.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Aune, D., Vardaxis, I., Lindqvist, B.H. et al. Dispensed prescription medications and short-term risk of pulmonary embolism in Norway and Sweden. Sci Rep 14, 20054 (2024). https://doi.org/10.1038/s41598-024-69637-4
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-024-69637-4