Abstract
BACKGROUND: The risk of skin cancer associated with antihypertensive medication use is unclear, although thiazides have been implicated in regulatory safety warnings. We aimed to assess whether use of thiazides and other antihypertensives is associated with increased rates of keratinocyte carcinoma and melanoma.
METHODS: We conducted a population-based inception cohort study using linked administrative health data from Ontario, 1998–2017. We matched adults aged ≥ 66 years with a first prescription for an antihypertensive medication (thiazides, angiotensin-converting enzyme inhibitors, angiotensin II receptor blockers, calcium channel blockers, β-blockers) by age and sex to 2 unexposed adults who were prescribed a non-antihypertensive medication within 30 days of the index date. We evaluated each antihypertensive class in a separate cohort study. Our primary exposure was the cumulative dose within each class, standardized according to the World Health Organization’s Defined Daily Dose. Outcomes were time to first keratinocyte carcinoma, advanced keratinocyte carcinoma and melanoma.
RESULTS: The inception cohorts included a total of 302 634 adults prescribed an antihypertensive medication and 605 268 unexposed adults. Increasing thiazide exposure was associated with an increased rate of incident keratinocyte carcinoma (adjusted hazard ratios [HRs] per 1 Defined Annual Dose unit 1.08, 95% confidence interval [CI] 1.03–1.14), advanced keratinocyte carcinoma (adjusted HR 1.07, 95% CI 0.93–1.23) and melanoma (adjusted HR 1.34, 95% CI 1.01–1.78). We found no consistent evidence of association between other antihypertensive classes and keratinocyte carcinoma or melanoma.
INTERPRETATION: Higher cumulative exposure to thiazides was associated with increased rates of incident skin cancer in people aged 66 years and older. Consideration of other antihypertensive treatments in patients at high risk of skin cancer may be warranted.
Keratinocyte carcinoma (also called nonmelanoma skin cancer) includes basal and squamous cell carcinoma, the 2 most common malignancies in North America.1,2 More than 76 000 and 5.4 million cases are diagnosed annually in Canada and the United States, respectively.1,2 Melanoma is the fifth most common malignancy in the US, with 83 362 annual cases, 10 000 deaths and incidence rates rising by 10% between 2005 and 2015.3–5 It is estimated that there were 8000 new melanoma cases and 1300 related deaths in Canada in 2020.6 Given the high and increasing population burden of these skin cancers, identifying modifiable risk factors is important.7,8
Ultraviolet (UV) radiation exposure is the most important environmental risk factor for skin cancer.9,10 Medication-induced phototoxicity, in which medications interact with UV radiation to cause cellular damage in the skin, may increase the carcinogenic potential of sun exposure.11 Antihypertensive medications are used by about 1 in 5 adults,12 and thiazide diuretics, calcium channel blockers, β-blockers, angiotensin II receptor blockers (ARBs) and angiotensin-converting enzyme (ACE) inhibitors have all been reported to be phototoxic.11,13 Hydrochlorothiazide is first-line pharmacotherapy for hypertension14 and is considered the most phototoxic antihypertensive medication.11,13 It has been implicated in increasing the risk of keratinocyte carcinoma, with 2 recent case–control studies15,16 prompting Health Canada, the European Medicines Agency and the US Food and Drug Administration to issue warnings regarding prolonged use.17–19 Weaker associations have been reported between thiazides and melanoma, 20 and there is conflicting evidence on the association between other antihypertensive classes and skin cancer.21 To determine whether exposure to thiazides or other antihypertensive medications is associated with skin cancer, we conducted a population-based inception cohort study.
Methods
We followed the Reporting of Studies Conducted Using Observational Routinely Collected Health Data Statement for Pharmacoepidemiology (RECORD-PE) guideline to ensure adequate reporting.22–24
Study population
We used linked population-based administrative health data for Ontario, Canada. These data sets were linked using unique encoded identifiers and analyzed at ICES (Appendix 1, Supplemental Methods, available at www.cmaj.ca/lookup/doi/10.1503/cmaj.201971/tab-related-content).
Fees for physician visits and hospital admissions in Ontario are covered for all residents under the Ontario Health Insurance Plan (OHIP),25 and medication costs for residents aged 65 years and older are covered by the Ontario Drug Benefit (ODB) plan. The ODB database has an accuracy rate > 99%.26 We limited the study to patients aged 66 years and older to allow ascertainment of incident medication exposures.
We included incident users of antihypertensive medications (exposed group) if they had a first-ever prescription for a thiazide, ACE inhibitor, ARB, calcium channel blocker or β-blocker between Jan. 1, 1998, and Dec. 31, 2016 (Appendix 1, Figure S1). The index date was defined as the date of their first filled antihypertensive prescription in ODB. We used a lookback window of at least 1 year, measured back to age 65 years or 1991 (whichever was most recent) to exclude people who had previous prescriptions for an antihypertensive, phototoxic or immunosuppressive medication for at least 90 days within any 180-day period. Because solid organ transplant recipients are at substantially higher risk for skin cancer,27 we looked back at least 17 years to 1981 to exclude any patient who received a transplant before the index date based on the Canadian Organ Replacement Register. We matched each exposed antihypertensive user to 2 nonusers (unexposed group) of the same age (± 30 d) and sex who filled a prescription for any non-antihypertensive medication (excluding known phototoxic or immunosuppressive medications; listed in Appendix 1, Supplemental Methods) within 30 days of the exposed patient’s index date; the date of the filled non-antihypertensive medication was the index date for the unexposed group. Nonusers could become classified as exposed if they subsequently filled an antihypertensive prescription after cohort entry.
We excluded people with a history of keratinocyte carcinoma or melanoma before the index date (lookback window of 6–25 years for keratinocyte carcinoma based on OHIP data availability since 1992; lookback window of at least 34 years to 1964 for melanoma based on Ontario Cancer Registry data availability). To account for latency and induction periods between medication exposure and its potential effect on skin cancer risk, we also excluded people with keratinocyte carcinoma or melanoma diagnosed within 365 days after the index date or with less than 1 year of follow-up owing to death or loss of OHIP coverage. The outcome observation window began after this 1-year lag.
From the overall cohort of incident users of antihypertensive medications, we created separate inception cohorts for each antihypertensive class based on the specific medication prescribed to each exposed user on their index date (i.e., their first prescribed antihypertensive medication). If multiple antihypertensive medications were prescribed on the index date, then we included the patient in the inception cohorts for each corresponding antihypertensive class. Patients contributed person-time only to the antihypertensive inception cohorts they entered on the index date; that is, although we examined each antihypertensive class as a time-varying covariate within each cohort (see Exposures and Statistical analysis sections), participants did not cross over between the separate inception cohorts after the index date.
Exposures
For each inception cohort, our exposure was the time-varying cumulative dose of the index antihypertensive class (thiazide diuretics, ACE inhibitors, ARB, calcium channel blockers, β-blockers; medications listed in Appendix 1, Supplemental Methods). To standardize the cumulative dose, we calculated a Defined Annual Dose for each medication using the World Health Organization’s Defined Daily Dose (www.whocc.no/atc_ddd_index/). We divided each patient’s cumulative dose by the total Defined Annual Dose (365 d × Defined Daily Dose) for each medication.
Outcomes
The primary outcome was time to first diagnosis of keratinocyte carcinoma (including basal and squamous cell carcinoma) identified using a validated OHIP claims–based algorithm (sensitivity 83%–85%, specificity 93%–99%).28 Our secondary outcomes were time to first advanced keratinocyte carcinoma and first melanoma. We defined advanced keratinocyte carcinoma as keratinocyte carcinoma causing death (International Classification of Diseases, 9th Revision 173 or International Statistical Classification of Diseases and Related Health Problems, 10th Revision C44 listed as a cause of death in the Ontario Registrar General Death database)29 or requiring complex surgery (Appendix 1, Supplemental Methods); this definition of advanced keratinocyte carcinoma has been used in previous studies but has not been validated.30 Cases of invasive melanoma were identified from the Ontario Cancer Registry, which captures 94% of melanomas diagnosed in Ontario.31
Covariates
We obtained baseline age, sex, location of residence (urban or rural) and socioeconomic status based on residential postal code (income quintile) from the Ontario Registered Persons Database. 32,33 We assessed history of hypertension before or on the index date using a validated administrative data algorithm.34 We derived the Charlson–Deyo comorbidity index using a 2-year look-back window to identify hospital discharge records in the Canadian Institute for Health Information.35,36 To account for detection bias associated with health care utilization,37 we assessed the number of outpatient and inpatient physician visits in the 12 months before the index date. To account for secular trends in prescription patterns and skin cancer diagnosis, we included the time period of the index date (1998–2000; 2001–2005; 2006–2010; 2011–2016) as a covariate. We used ODB to identify time-varying exposure (updated monthly) to immunosuppressant and other phototoxic medications for a clinically meaningful duration of at least 180 days in any 365-day period during follow-up.
Statistical analysis
We conducted analyses separately for each of the 5 antihypertensive class inception cohorts. We used multivariable cause-specific hazards regression and restricted cubic spline functions with 5 knots (5th, 27.5th, 50th, 72.5th and 95th percentiles) to model the associations between the index drug exposure and each skin cancer outcome (keratinocyte carcinoma, advanced keratinocyte carcinoma, melanoma), while accounting for the competing risk of death. The primary exposure in each analysis was the cumulative dose of the index antihypertensive class of interest, modelled as a continuous time-varying covariate. We used 0 to 0.5 Defined Annual Dose units as the lowest dose category (rather than 0) to avoid collinearity with an ever or never use covariate and because short-term exposure is unlikely to meaningfully increase skin cancer risk. Patients who were unexposed at baseline contributed person-time in the unexposed group; if they were later prescribed the exposure antihypertensive class (not as their first antihypertensive, by definition), subsequent person-time was attributed to the exposed group. Patients were followed from their index date until their first skin cancer, death, loss of OHIP eligibility for ≥ 90 days (outmigration from Ontario) or Dec. 31, 2017, whichever was earliest.
Each model evaluating cumulative medication exposure included ever or never use of each antihypertensive class as a binary, time-varying covariate. The inclusion of this covariate controls for treatment selection bias (confounding by indication) due to unmeasured factors associated with the propensity to be prescribed the medications.38 In every model, we adjusted for the time-varying cumulative dose of each of the other 4 nonindex antihypertensive classes.38 We adjusted each model for the fixed baseline and time-varying covariates listed above. We estimated associations as adjusted hazard ratios (adjusted HRs) and 95% confidence intervals (CIs). We considered 2-sided p < 0.05 to be statistically significant. We tested the proportional hazards assumption for time-independent covariates in each model by checking their interaction with the logarithm of time.
In additional analyses to promote clinical interpretability, we fit the above models with cumulative Defined Annual Dose units for each medication class, as follows: 0 to 0.5, > 0.5 to 1.5, increasing in single unit increments up to > 7.5. We used contrast statements with orthogonal polynomial coefficients to calculate p for linear trend for increasing cumulative dose categories. We also evaluated these ordinal exposures as cumulative duration of use, irrespective of the dose. In sensitivity analyses, we extended the lag time to exclude patients with their first keratinocyte carcinoma diagnosis or end of follow-up occurring within 2 and 3 years after the index date.
Ethics approval
Owing to use of routinely collected data, this project was exempt from ethics review, as authorized under section 45 of Ontario’s Personal Health Information Protection Act.
Results
We identified 2 630 956 patients who were prescribed an antihypertensive medication between 1998 and 2016. After applying our exclusion criteria, we included a total of 302 634 patients exposed to antihypertensive medications at baseline and 605 268 matched unexposed patients in our inception cohorts (Figure 1). Patients exposed to antihypertensive medications but excluded because they could not be matched (n = 218 614) tended to be older, have more preceding physician visits and have a more recent index date, but were otherwise similar to those included (Appendix 1, Table S1). Angiotensin-converting enzyme inhibitors were the most commonly first prescribed antihypertensive medication class (n = 112 723, 37%), followed by thiazide diuretics (n = 87 725, 29%), β-blockers (n = 68 218, 23%), calcium channel blockers (n = 50 533, 17%) and ARB (26 317, 9%).
Baseline characteristics for the matched thiazide inception cohort are given in Table 1. The number of physician visits and Charlson–Deyo comorbidity index at baseline were similar between exposed and matched unexposed patients, but the former were more likely to have been diagnosed with hypertension (70% v. 23%). Baseline characteristics for the other antihypertensive inception cohorts are in Appendix 1, Tables S2–S5.
In the thiazide inception cohort, 8768 of 87 525 (10%) originally exposed and 18 706 of 175 050 (11%) matched, originally unexposed patients were diagnosed with a keratinocyte carcinoma (Appendix 1, Table S6), 1302 (1.5%) exposed and 2559 (1.5%) unexposed were diagnosed with an advanced keratinocyte carcinoma and 296 (0.3%) exposed and 564 (0.3%) unexposed were diagnosed with a melanoma. There were 2 426 593 total person-years of follow-up (median 8.9 years, interquartile range 4.6–13.7), with 40 deaths attributed to keratinocyte carcinoma. In each baseline group, 34% of patients died of causes unrelated to keratinocyte carcinoma during follow-up.
Increasing thiazide exposure was associated with increased rates of incident keratinocyte carcinoma (adjusted HR per Defined Annual Dose unit increase 1.08, 95% CI 1.03–1.14) and melanoma (adjusted HR 1.34, 95% CI 1.01–1.78) (Appendix 1, Table S7). The association between thiazide exposure and advanced keratinocyte carcinoma was of similar magnitude to the association with keratinocyte carcinoma overall, but with wider confidence intervals (adjusted HR 1.07, 95% CI 0.93–1.23). Restricted cubic spline curves show strengthening of the association between thiazides and keratinocyte carcinoma (Figure 2) and advanced keratinocyte carcinoma (Figure 3) with increasing cumulative exposure, whereas effect estimates became increasingly imprecise for melanoma beyond 10 Defined Annual Dose units (Figure 4). In the analysis with thiazides modelled as an ordinal exposure, patients exposed to > 7.5 Defined Annual Dose units had a 44% higher rate of keratinocyte carcinoma (adjusted HR 1.44, 95% CI 1.35–1.54, Appendix 1, Table S8), 52% higher rate of advanced keratinocyte carcinoma (adjusted HR 1.52 [95% CI 1.29 to 1.79], Appendix 1, Table S9) and 60% higher rate of melanoma (1.60 [95% CI 1.10–2.32], Appendix 1, Table S10) compared with 0–0.5 Defined Annual Dose units.
Along with increasing dose of thiazides, other variables independently associated with keratinocyte carcinoma in the inception cohort for thiazides included age (per additional year, adjusted HR 1.02, 95% CI 1.02 adjusted 1.02), male sex (adjusted HR 1.42, 95% CI 1.38–1.45), rural residence (adjusted HR 1.12, 95% CI 1.08–1.16), being in the highest 2 income quintiles (adjusted HR 1.32, 95% 1.29–1.36), immunosuppressant drug use (adjusted HR 1.38, 95% CI 1.29–1.48) and phototoxic drug use (adjusted HR 1.09, 95% CI 1.06–1.12) (Appendix 1, Table S11).
In the inception cohorts for the other antihypertensive classes, associations between other antihypertensive classes and keratinocyte carcinoma, advanced keratinocyte carcinoma and melanoma were small, inconsistent or statistically nonsignificant (Figures 2–⇑4; Appendix 1, Table S7). Results were similar in secondary analyses where exposures were modelled as ordinal cumulative duration variables (Appendix 1, Tables S8–10), and in sensitivity analyses that excluded patients with keratinocyte carcinoma diagnoses or end of follow-up within 2 and 3 years after starting their first antihypertensive medication (Appendix 1, Table S12).
Interpretation
In this population-based inception cohort study among older adults in Ontario, we found a significant increase in the rate of keratinocyte carcinoma and melanoma associated with increasing cumulative thiazide diuretic use. Exposure to other antihypertensive classes, which are not thought to be as phototoxic,11,13 either did not affect the rate of keratinocyte carcinoma and melanoma or had inconsistent or weak associations.
Our finding that thiazide diuretics increased the rate of keratinocyte carcinoma is consistent in direction with recent case–control studies from Denmark.15,16 The authors found that cumulative hydrochlorothiazide dosages over 50 000 mg (equivalent to 5.5 Defined Annual Dose units) were associated with 29% and 400% increased odds of basal cell carcinoma and squamous cell carcinoma, respectively.15 In a meta-analysis evaluating the association between thiazide diuretics and keratinocyte carcinoma, thiazide diuretics were associated with increased risk of both basal and squamous cell carcinoma.20 Our finding of an association between thiazide use and melanoma is also consistent in direction with that meta-analysis, which found an odds ratio of 1.14.20 Many of the included studies were limited by deficiency of data on cumulative thiazide dosage, immortal time bias, a lack of accounting for polypharmacy and treatment selection bias.15,16,20,39–54 In a systematic review of keratinocyte carcinoma risk associated with all antihypertensive classes, use of β-blockers and calcium channel blockers was associated with a slightly increased risk for basal (odds ratio 1.09 for β-blockers and 1.15 for calcium channel blockers) but not squamous cell carcinoma. 21 Findings for ACE inhibitors and ARBs were heterogeneous.21
Our study addresses limitations of previous studies, including having relevant comparison groups, time-varying cumulative exposure variables, large sample size and standardized measurement of cumulative dose. We also adjusted for concomitant use of other antihypertensive medications55 and adjusted for treatment selection bias by including a binary ever or never use variable for each antihypertensive class.38
Limitations
We assessed 3 different outcomes for 5 separate antihypertensive drug cohorts, and the multiple comparisons could have inflated the chance of a false-positive result. Inherent in the observational design is the potential for unmeasured confounding, although we adjusted for treatment selection bias, as noted above. For example, our administrative data do not include UV radiation exposure, race or skin type, yet in the ethnically diverse population of Ontario 29% of the population identify as visible minorities.56 We used a new-user design, but it is possible that patients were exposed to antihypertensive medications before our minimum 1-year lookback window, although these are chronic disease medications that are unlikely to be stopped for such a prolonged time. Our study was restricted to older adults; associations between antihypertensive medications and skin cancer in younger populations may be different owing to their lower lifetime UV radiation exposure and keratinocyte carcinoma incidence. We were unable to identify squamous and basal cell carcinomas separately. Detection bias may partly explain the association between thiazides and skin cancer;37 we adjusted for physician visits and are reassured by the lack of association between the other antihypertensive classes and skin cancer.
Conclusion
We found consistent dose-dependent increases in skin cancer risk associated with thiazides but not with other antihypertensive classes. Although not evaluated in our study, interventions that may decrease the burden of skin cancer for patients taking thiazide diuretics include sun protection and monitoring for early detection of skin cancer. Clinicians may consider alternatives to thiazide diuretics to treat hypertension in patients at high risk for skin cancer. Future studies should replicate our analyses in other populations and geographies, as skin cancer risk differs by latitude even within Canada.57,58 Studies on the impact of current regulatory safety warnings and whether they or other interventions decrease the burden of skin cancer associated with thiazide use are needed.
Acknowledgement
The authors thank Dr. Allan Detsky for his input on a draft of this manuscript.
Footnotes
Competing interests: Aaron Drucker reports receiving consulting fees from Sanofi, RTI Health Solutions, Eczema Society of Canada and Canadian Agency for Drugs and Technology in Health and honoraria from CME Outfitters. Dr. Drucker’s institutions have received educational grants from Sanofi (Women’s College Hospital ) and research grants from Sanofi and Regeneron (Brown University). Martin Weinstock has received consulting fees from AbbVie and Almirall. Husam Abdel-Qadir reports receiving personal fees from Amgen Canada and from the Canadian Vigour Centre, for time spent on an end-point adjudication committee for the THEMIS clinical trial. Dr. Abdel-Qadir also reports receiving grants from the Heart and Stroke Foundation of Canada, Canadian Institutes of Health Research, and the Canadian Cardiovascular Society. No other competing interests were declared.
This article has been peer reviewed.
Contributors: All of the authors contributed to the conception and design of the work. All of the authors contributed to the acquisition, analysis and interpretation of data. Aaron Drucker drafted the manuscript. All of the authors revised the manuscript critically for important intellectual content, gave final approval of the version to be published and agreed to be accountable for all aspects of the work.
Funding: This work was funded by an operating grant from the Canadian Dermatology Foundation; the funder had no role in the design, conduct, interpretation of or decision to publish the study.
Data sharing: The data set from this study is held securely in coded form at ICES. While data sharing agreements prohibit ICES from making the data set publicly available, access may be granted to those who meet prespecified criteria for confidential access, available at www.ices.on.ca/DAS. The full dataset creation plan and underlying analytic code are available from the authors upon request, understanding that the computer programs may rely upon coding templates or macros that are unique to ICES and are therefore either inaccessible or may require modification.
Disclaimer: This study was supported by ICES, which is funded by an annual grant from the Ontario Ministry of Health and Long-Term Care (MOHLTC). This study also received funding from an operating grant from the Canadian Dermatology Foundation. The analyses, conclusions, opinions and statements expressed herein are solely those of the authors and do not reflect those of the funding or data sources; no endorsement is intended or should be inferred. Parts of this material are based on data and/or information compiled and provided by CIHI. However, the analyses, conclusions, opinions and statements expressed in the material are those of the author(s), and not necessarily those of CIHI. We thank IMS Brogan Inc. for use of their Drug Information Database. Parts of this material are based on data and information provided by Cancer Care Ontario (CCO). The opinions, results, view, and conclusions reported in this paper are those of the authors and do not necessarily reflect those of CCO. No endorsement by CCO is intended or should be inferred. Parts of this report are based on Ontario Registrar General information on deaths, the original source of which is ServiceOntario. The views expressed therein are those of the author and do not necessarily reflect those of Ontario Registrar General or Ministry of Government Services. As a prescribed entity under Ontario’s privacy legislation, ICES is authorized to collect and use health care data for the purposes of health system analysis, evaluation and decision support. Secure access to these data is governed by policies and procedures that are approved by the Information and Privacy Commissioner of Ontario.
- Accepted February 22, 2021.
This is an Open Access article distributed in accordance with the terms of the Creative Commons Attribution (CC BY-NC-ND 4.0) licence, which permits use, distribution and reproduction in any medium, provided that the original publication is properly cited, the use is noncommercial (i.e., research or educational use), and no modifications or adaptations are made. See: https://creativecommons.org/licenses/by-nc-nd/4.0/