Abstract
Background:
Tamoxifen (TAM) and aromatase inhibitor (AI) therapies have been associated with increased risk of thromboembolic and cardiovascular events, respectively, in addition to other side effects. This study analysed the risk of these events and the overall survival (OS) benefit in breast cancer patients treated with AI, compared with TAM-treated patients, in a large population-based cohort.
Methods:
This observational cohort study included women diagnosed with breast cancer and treated with TAM or AI. Data were extracted from primary care records in a population database (SIDIAP, System for the Development of Research in Primary Care). Incidence rates of study outcomes are reported. Survival analyses included Kaplan–Meier estimation and Cox proportional hazards models. Sensitivity analysis was carried out, using Fine and Gray models to account for competing risk of death. Confounding was minimized using propensity score adjustment and inverse probability weighting (IPW) adjustment.
Results:
Data from 3082 postmenopausal women treated with TAM, and 18,455 treated with AI, were available. Adjusted hazard ratios (HRs) [95% confidence interval (CI)] for AI users, compared with TAM group, were 0.93 (95%CI 0.69–1.26) for thromboembolic events (TEEs); 1.13 (95%CI 0.79–1.63) for cardiovascular events, and 0.76 (95%CI 0.70–0.82) for mortality. Additional analyses using competing risk analysis had similar results, while IPW adjustment showed a potential risk of pulmonary embolism (PE) [2.26 (95%CI 1.02–4.97)] in AI-treated patients.
Conclusions:
AI users had >20% lower all-cause mortality compared with TAM users, without increasing risk to experience cardiovascular and TEEs. This would locate AI therapy on the first line in clinical practice. Thus, AI might be the most preferable option in adjuvant hormonal therapy choice.
Keywords
Introduction
Aromatase inhibitors (AIs) and tamoxifen (TAM) are known to be effective adjuvant endocrine therapies for patients with hormone-receptor-positive breast cancer. Generally, these patients have good prognosis, with an overall survival (OS) rate exceeding 80%.1,2 However, these therapies have been associated with side effects that can affect quality of life and could impact on mortality, among them, cardiovascular events (CVEs) and thromboembolic events (TEEs) are emerging as competing causes of death. 3 A number of randomized controlled trials (RCT) have explored the cardiovascular effect, comparing AIs versus TAM, with heterogeneous results.4,5 These studies have provide evidence of increased CVEs associated with AI therapies, compared with TAM, likely due to high depletion of estradiol levels and alteration of lipid metabolism related to AIs, 6 or to the cardioprotective role of tamoxifen per se. 7
Although RCT and meta-analysis are gold standard experimental approaches for the study of efficacy and safety in ‘ideal’ conditions, they are sometimes not representative of clinical practice conditions or of the actual profile of the treated community, 8 and, thus, they cannot address definitively safety issues, particularly for side effects with low incidence. We therefore aimed to analyze the risk of CVE and TEE, and OS benefit during AI therapy, compared with TAM, in real-world conditions.
Thus, the present study used the SIDIAP (System for the Development of Research in Primary Care) database, which provided anonymized clinical information as coded by primary care practitioners in Catalonia, Spain, covering more than 7 million patients. 9 SIDIAP contains information on socio-demographics and extended clinical data. Moreover, SIDIAP is linked to pharmacy invoice data, which provides detailed information on drugs dispensed in community pharmacies under the universal health care system. Using this database, we performed a population-based study, including almost 28,000 women treated with AI or TAM for up to 10 years of follow up, to assess thromboembolic and cardiovascular events, and resulting OS in general clinical practice.
Methods
Data sources
SIDIAP (http://www.sidiap.org) is an anonymized clinical database of more than 7 million patient records collected from more than 370 primary care teams covering >80% of the total population of Catalonia. Among the available variables are socio-demographic data, lifestyle risk factors, prescriptions dispensed and comorbidities. Health professionals gather this information using ICD-10 codes and structured forms designed for the collection of clinical factors (alcohol use, smoking, body mass index, etc.). Migration out of the catchment area is also recorded, allowing for longitudinal follow up of patients. 10 Death is also registered in the SIDIAP database, as provided by the universal health insurance database for Catalonia (in Catalan, ‘registre central de persones assegurades’).
Study design and participants
Retrospective observational cohort study of women diagnosed with early breast cancer, defined as nonmetastatic breast cancer, stage I–III, and treated with monotherapy of TAM or AIs as registered in the SIDIAP database from January 2006 to December 2015. Therapeutic regimen in patients was identified by its anatomical therapeutic chemical classification in pharmacy dispensing records, coded as L02BG for AIs (L02BG03 for anastrozole, L02BG04 for letrozole, L02BG06 for exemestane), and L02BA01 for TAM.
Exclusion criteria were previous history of cancer (except nonmelanoma skin cancers) and patients who had received a switching therapy (TAM followed by AI or vice versa). No concomitant anticancer drugs other than TAM or AI were used.
Ethics statement
This study used only data collected routinely from the SIDIAP database. The Idiap Jordi Gol Research Ethics Committee and the SIDIAP Database Scientific Committee have approved the study protocol (P16/031). No human subjects or tissues were used in this study. Data provided by SIDIAP was anonymized and risk of identification was almost null according to Spanish law LO 15/ 1999 13 December. Thus, informed consent did not need to be obtained from participants.
Follow up
Participants were followed up from therapy initiation (first day of TAM or AI dispensing) until the earliest of three endpoints: treatment cessation (defined by a refill gap of 6 months or more with no dispensation of the index therapy), plus 1 month wash-out (for carry-over effects); evaluated outcomes date (as recorded in electronic medical records); or death, migration out of catchment area, or end of SIDIAP data availability (31 December 2015).
In the overall mortality assessment, patients were followed-up during all the study period (2006–2015).
Variables
Outcomes of the study
Analyzed outcomes were the first TEE [pulmonary embolism (PE) and deep vein thrombosis (DVT), including phlebitis and thrombophlebitis] and the first CVE [coronary artery disease (CAD), and cerebrovascular diseases (CVD), including stroke and intracerebral haemorrhage, among others] occurring during adjuvant therapy. In addition, PE, DVT, CAD and CVD were analysed separately as secondary outcomes. OS, expressed as mortality status during follow up, was also reported. ICD-10 codes used to identify the outcomes of the study are documented in Supplementary Table S1.
Confounders
A prespecified list of confounders was extracted from SIDIAP, informed by previous clinical knowledge and scientific literature. These confounding factors fell into five clusters:
Sociodemographics: age (at treatment initiation), body mass index (BMI), and socioeconomic status (assessed by MEDEA, a validated deprivation index). 11
Menopausal status: defined as women >55 years old at diagnosis in the TAM group or all patients treated with AI. Menopausal status of patients <55 years old in TAM group is unknown.
Lifestyle factors: smoking, alcohol use (defined according to The Catalan Health Care System: none/low, as a mean of 0 g of alcohol per week; moderate, not exceeding 170 g of alcohol per week; high/alcoholic, 170 g of alcohol or more per week).
Past medical history: Charlson comorbidity index, and previous history of CVE and TEE.
Concomitant use of antiplatelets or anticoagulants or statins at cohort entry (i.e. TAM/AI initiation).
Statistical analysis
Data from SIDIAP were managed using MySQL. Differences in baseline characteristics between TAM and AI participants were described and imbalances analysed using t test and Chi-square test.
Incidence rates of study outcomes (during treatment for TEE/CVE and at any time for OS) were estimated. 12
For each outcome, survival analysis was done by Kaplan–Meier estimation and Cox proportional hazards model to estimate cumulative probability plots and HRs according to treatment/exposure, respectively.
Additionally, a subanalysis using Fine and Gray regression models were fitted to estimate subdistribution hazard ratios (SHR) for TEE and CVE (separately) according to treatment arm, accounting for a competing risk of death. 13
HR and SHR are reported with 95% CI, and using TAM as a reference group (and AI as the ‘exposed’ group). Moreover, the assumption of proportionality was verified though proportional hazards assumption for a Cox regression model test.
Adjustment in survival analysis was conducted using the propensity score (PS). PS was estimated using logistic regression models, where treatment group was the outcome and the previously listed confounders were adjusted for. The final list of variables used in PS adjustments are listed in Table 1, including statins, anticoagulants and antiplatelet drugs. Missing data were imputed before the PS estimates, using multiple imputation by chained equations, obtaining 10 imputed datasets that were combined using Rubin’s rules. 14 Previous TEE and CVE history were included in their respective analyses, and in OS evaluation. Additional analysis censoring patients with previous TEE and CVE was performed to account for potential baseline higher risk.
Baseline characteristics of candidates in postmenopausal women.
Participants included in TAM group were older than 55 years old to ensure postmenopausal status.
AI, aromatase inhibitors; BMI, body mass index; CVE, cerebrovascular event; Q, quartile; QMEDEA, quintile MEDEA deprivation index; SDE, standard deviation; TAM, tamoxifen; TEE, thromboembolic event.
An additional analysis adjusting survival analysis by stabilized inverse probability weighting (IPW) and using a robust sandwich-type variance estimator was performed. 15
In order to compare our results with an analysis not accounting for menopausal status, we repeated the same models using total TAM users, including those women younger than 55 years old.
All statistical analyses were conducted using R for Windows version 3.3.3 and the following R packages: foregin, Hmisc, compareGroups, survival and mice.
Results
Of the 36,472 eligible participants, 21,537 (18,455 AI and 3082 TAM) monotherapy users were included in the analysis (see Figure 1), with a median (interquartile range) of treatment of 29 (10–53) months and a maximum of 119 months. Baseline characteristics of participants are presented in Table 1. Postmenopausal women treated with TAM were older, less likely to be current smokers, concomitant users of anticoagulants or statin therapy, and had lower prevalence of previous TEE than AI users; but had a similar BMI, Charlson comorbidity index, current alcohol drinker status, users of platelet inhibitors and prevalence of previous CVE.

Flow chart of SIDIAP cohort study.
Thromboembolic adverse events
A total of 49 patients in the TAM group experienced TEEs (1.59%), whereas these were 345 (1.87%) patients in the AI group. This is equivalent to incidence rates of 8.16/1000 person-years (95%CI 6.10–10.69) in TAM users, and 6.93/1000 person-years (95%CI 6.23–7.69) in AI patients (Figure 2a). No significant differences in thromboembolic risk were observed between both therapies [adjusted HR 0.93 (95%CI 0.69–1.26)] (Table 2). Survival analysis, adjusted for competing risk, showed similar findings (Table 2). Sensitivity analyses excluding patients with previous TEE did not change the estimates (data not shown).

Cumulative hazard plot of (a) TEEs, (b) CVEs and (c) mortality in SIDIAP cohort according to drug therapy. Graphs show Kaplan–Meier curves representing the outcomes of the study in terms of cumulative hazards.
Thromboembolic, cardiovascular and mortality risk of AI treatment compared with TAM treatment in postmenopausal women.
Adjusted results were obtained using continuous PS estimates.
AI, aromatase inhibitors; CI, confidence interval; CVEs, cardiovascular events; HR, hazard ratio; PS, propensity score; SHR, subdistribution hazard ratio; TAM, tamoxifen; TEEs, thromboembolic events.
Cardiovascular adverse events
A total of 33 (1.07%) TAM users had at least one CVE event during follow up, compared with 271 events (1.47%) in the AI user group. Incidence rates were therefore 5.50/1000 person-years (95%CI 3.85–7.63) in TAM, and 5.43/1000 person-years (95%CI 4.81–6.10) in AI users. Cumulative hazard plots of CVEs are shown in Figure 2b. No significant increase in cardiovascular risk was detected in AI-treated patients [adjusted HR 1.13 (95%CI 0.79–1.63)] (Table 2). Survival analysis, adjusted for competing risk, show similar findings (Table 2). Sensitivity analyses excluding patients with previous CVE did not change the estimates (data not shown).
Mortality
Overall mortality was 22.58% (696 participants) in the TAM group, and 19.75% (3644 subjects) in the AI group. Crude mortality rates were 40.68/1000 person-years (95%CI 37.74–43.78) in TAM, and 40.25/1000 (95%CI 38.95–41.57) in AI users. Cumulative hazard plots of mortality are shown in Figure 2c. Adjusted Cox models showed a better prognostic for AI users [HR of 0.76 (95%CI 0.70–0.82)] compared with TAM users (Table 2). Similar findings were observed after competing risk adjustment (Table 2).
Secondary outcomes
TEEs: PE and DVT
In our cohort, 100 PE events [7 in TAM group, incidence rate 1.17 (95%CI: 0.51–2.31); and 93 in AI group, incidence rate 1.87 (95%CI: 1.52–2.28)] and 294 DVTs [42 in TAM group, incidence rate 6.99 (95%CI: 5.10–9.36); and 252 in AI group, incidence rate 5.06 (95%CI: 4.47–5.72)] were reported. No differences in DVT risk were found between both groups. A nonsignificant increased risk of PE (Adjusted SHR of 2.15 [95%CI 0.99–4.64]) in AI group was observed (Table 3).
Risk of PE, DVT, CAD and CVD of AI treatment compared with TAM treatment in postmenopausal woman.
Adjusted results were obtained using continuous PS estimates.
AI, aromatase inhibitors; CAD, coronary artery disease; CI, confidence interval; CVD, cerebrovascular diseases, including stroke and intracerebral hemorrhage; CVEs, cardiovascular events; DVT, deep vein thrombosis, phlebitis and thrombophlebitis; HR, hazard ratio; PE, pulmonary embolism; PS, propensity score; SHR, subdistribution hazard ratio; TAM, tamoxifen; TEEs, thromboembolic events.
CVEs: CAD and CVD
Of 304 CVEs, 292 were CAD [32 in TAM, incidence rate 5.33 (95%CI: 3.71–7.44); and 260 in AI users, incidence rate 5.21 (95%CI: 4.60–5.87)] and 12 were CVD [1, incidence rate 0.17 (95%CI:0.01–0.82) in TAM; and 11, 0.22 (95%CI: 0.12–0.38) in AI users]. No increased risk of CAD nor CVD was detected between both groups (Table 3).
A secondary analysis where patients with previous CVE and TEE were censored showed similar results (data not shown).
Inverse probability weighting analysis
Survival analyses adjusted for IPW are reported in Table 4. CVEs were similar to those obtained in Cox regression models, with no significant difference in risk between AI and TAM users. Likewise, lower mortality risk was detected in AI users. On the contrary, an increased risk of PE was detected in AI users [stabilized IPW HR 2.26 (95%CI 1.02–4.97)].
Risk of TEEs and CVEs in AI treatment compared with TAM treatment in postmenopausal women using stabilized IPW adjustment.
AI, aromatase inhibitors; CAD, coronary artery disease; CI, confidence interval; CVD, cerebrovascular diseases, including stroke and intracerebral hemorrhage; CVEs, cardiovascular events; DVT, deep vein thrombosis, phlebitis and thrombophlebitis; IPW HR, Inverse probability weighting hazard ratio; PE, pulmonary embolism; TAM, tamoxifen; TEEs, thromboembolic events.
Survival analysis not accounting for menopausal status (including TAM women younger than 55 years old) showed an increased risk of CVE in AI patients [adjusted SHR 1.96 (95%CI 1.37–2.81)]. Nevertheless, the IPW analysis did not found differences between groups [HR 0.87 (95%CI 0.56–1.38)] (Supplementary Tables S2–5).
Discussion
In our cohort study, including 21,537 women with a breast cancer diagnosis treated in actual practice conditions, similar CVE and TEE risk was observed in both AI and TAM treatment groups. In addition, AI use seems to be a beneficial choice in terms of overall mortality reduction, with a >20% lower rate in AI users after adjusting for potential confounders.
A further analysis adjusting for IPW corroborated these findings but suggested a potential increase of PE associated with AI use.
In general, previous data on cardiovascular risk provided by RCT and meta-analyses, despite being heterogeneous, have identified a potential excess CVE risk associated with AI therapy, with slightly increasing odds of developing CVD, compared with patients receiving TAM therapy (OR 1.26, 95%CI = 1.10–1.43). 16 In this regard, Abdel-Qadir and colleagues published a similar population-based study using routinely collected data from Canada that observed an increased risk of myocardial infarction in AI users (HR 2.02; 95%CI = 1.16–3.53), but, exploring a lower-risk subgroup of patients aged <74 years, with stage I–II breast cancer and no prior ischaemic heart disease, this excess risk was not detected (HR 1.53; 95%CI = 0.41–5.71, P = 0.53). 17
In order to minimize imbalances when comparing treatment groups, we adjusted for a long list of potential confounders using propensity score equations, which is recommended in this type of study. 18 An additional analysis adjusting for IPW was performed to correct for potential attrition bias, confirming that AI-treated women did not experience an increased risk of CVE compared with postmenopausal women in the TAM group. Consistent with this, other population studies selecting older women also reported similar incidences in stroke and several heart diseases between treatment groups.19,20 All together, these findings suggest that the increased risk detected in AI users compared with total TAM users is driven mainly by menopausal status.
In contrast with previous RCT findings, 3 an increasing risk of PE was observed after IPW adjustment of our data: patients treated with AIs had twice the risk of PE compared with TAM users. However, this was a post hoc analysis based on a limited number of events, and needs further confirmation in external cohorts. Further research is needed to explore this potential association between AI treatment and PE risk.
As a side note, and similar to the results of Abdel-Qadir and colleagues, 17 our supplemental survival analysis adjusted by computing risk of death including all TAM users (lower and older than 55 years old) suggested the AI group has nearly twice the risk of CVEs compared with the total TAM group, but later selection of postmenopausal women has shown similar hazard risks for all events in both treatment groups, proposing that menopause status is not a confounder but a potential interaction.
In addition to the lack of association between AI and CVE observed in our population, the significant OS in AI treated patients places these drugs in front to TAM in terms of cardiovascular safety and efficacy on recurrence incidence. 21 It is noteworthy that selective oestrogen receptor modulators (SERMs) have been associated with higher proportion of adverse drug reaction (ADR) reports related to QT prolongation, Torsade de Pointes, and ventricular arrhythmias compared with AIs, in the European database of suspected ADR reports. Nonetheless, the overall number of these events was very small. 22
One limitation of the study is that data of previous exposure to chemotherapy or radiotherapy were not available. In any case, these treatments are given independently of endocrine therapy election. 23 Hence, potential toxic effects in heart would be allocated randomly among patients. In addition, we performed a sensitivity analysis using IPW to minimize the presence of a potential bias by indication due to other factors. Likewise, data of severity, grade and the clinical stage of breast cancer were not accessible. However, TAM and AI monotherapies are recommended and used mainly for hormone receptor-positive early breast cancer. 24 In our study, we excluded patients with sequential therapies (TAM/AI) and further studies analysing outcomes in these patients might provide additional safety data.
Additionally, the SIDIAP data were collected during routine clinical practice (not by an expert researcher), potentially limiting the validity of coded outcomes. However, high accuracy in coding for all of the study outcomes was previously validated in the SIDIAP database. 10
The main strength of the study is the sample size, almost 22,000 participants, reflecting real population conditions. Likewise, the SIDIAP dataset includes all treatment centres and has the potential to include all patients in the source population, increasing the external validity of our findings.
In summary, no difference in CVD was observed between postmenopausal AI and TAM users. Furthermore, AI users had >20% lower all-cause mortality, yielding a positive risk–benefit for long-term use of these therapies.
Supplemental Material
21-11-2019_Supplementary_data – Supplemental material for Thromboembolic, cardiovascular and overall mortality risks of aromatase inhibitors, compared with tamoxifen treatment: an outpatient-register-based retrospective cohort study
Supplemental material, 21-11-2019_Supplementary_data for Thromboembolic, cardiovascular and overall mortality risks of aromatase inhibitors, compared with tamoxifen treatment: an outpatient-register-based retrospective cohort study by Marta Pineda-Moncusí, Natalia Garcia-Giralt, Adolfo Diez-Perez, Ignasi Tusquets, Sonia Servitja, Joan Albanell, Daniel Prieto-Alhambra and Xavier Nogués in Therapeutic Advances in Medical Oncology
Footnotes
Acknowledgements
Authors Marta Pineda-Moncusí and Natalia Garcia-Giralt contributed equally to this work.
Funding
The authors disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: The research was supported by the National Institute for Health Research (NIHR) Oxford Biomedical Research Centre (BRC). The views expressed are those of the author(s) and not necessarily those of the National Health Service (NHS), the NIHR or the Department of Health. The research was supported by Centro de Investigación Biomédica en Red de Fragilidad y Envejecimiento Saludable (grant number CB16/10/00245); Formación en Investigación en Salud grants (grant numbers PI13/00444, PI16/00818) from Instituto de Salud Carlos III; and European Regional Development Funds.
Conflict of interest statement
ADP reports grants from Amgen, personal fees from UCB, Mereo, Sandoz, Gilead, and other from Active Life Sci, outside the submitted work. DPA reports unrestricted research funding from UCB, Amgen, and Servier Laboratoires outside of the submitted work; DPA’s department has received fees for speaker services from Amgen outside of the submitted work. MPM, NGG, IT, SS, JA and XN have nothing to disclose.
Supplemental material
Supplemental material for this article is available online.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
