Hormone Replacement: July-August 2024

Hormone replacement therapy and cancer mortality in women with 17 site-specific cancers: a cohort study using linked medical records

Chris R. CardwellTom A. RangerAlexander M. LabeitCarol A. C. CouplandBlánaid HicksCarmel Hughes

Úna McMenaminXue W. MeiPeter Murchie & Julia Hippisley-Cox

British Journal of Cancer  volume 131, pages 737–746 (2024), June 2024

British Journal of Cancer (Br J Cancer) ISSN 1532-1827 (online) ISSN 0007-0920 (print)

Abstract

Background

There is limited evidence on the safety of Hormone Replacement Therapy (HRT) in women with cancer. Therefore, we systematically examined HRT use and cancer-specific mortality in women with 17 site-specific cancers.

Methods

Women newly diagnosed with 17 site-specific cancers from 1998 to 2019, were identified from general practitioner (GP) records, hospital diagnoses or cancer registries in Scotland, Wales and England. Breast cancer patients were excluded because HRT is contraindicated in breast cancer patients. The primary outcome was time to cancer-specific mortality. Time-dependent Cox regression models were used to calculate adjusted hazard ratios (HR) and 95% confidence intervals (95% CIs) for cancer-specific mortality by systemic HRT use.

Results

The combined cancer cohorts contained 182,589 women across 17 cancer sites. Overall 7% of patients used systemic HRT after their cancer diagnosis. There was no evidence that HRT users, compared with non-users, had higher cancer-specific mortality at any cancer site. In particular, no increase was observed in common cancers including lung (adjusted HR = 0.98 95% CI 0.90, 1.07), colorectal (adjusted HR = 0.79 95% CI 0.70, 0.90), and melanoma (adjusted HR = 0.77 95% CI 0.58, 1.02).

Conclusions

We observed no evidence of increased cancer-specific mortality in women with a range of cancers (excluding breast) receiving HRT.

 

Background

Hormone replacement therapy (HRT) is widely used to reduce menopausal related vasomotor symptoms (hot flushes, and night sweats) [1], urogenital atrophy [1] and postmenopausal osteoporosis [2]. It has also been shown to reduce joint pain, mood swings, sleep disturbances and improve quality of life [1]. The earlier detection and improved survival of patients with cancer has led to increasing numbers of women with cancer experiencing menopausal symptoms. In the United Kingdom (UK) HRT is contraindicated in patients with breast cancer and oestrogen-dependent cancers [2,3,4] but is not contraindicated in patients with other cancers. Clinicians may be reluctant to treat menopausal symptoms in patients with cancer using HRT given uncertainty around the impact of HRT on cancer outcomes, but denial of HRT without clear indication has been criticised by some authors as it could lead to unnecessary suffering [56].

Some researchers summarising evidence on the safety of HRT use in patients with cancer have advised caution when prescribing HRT to patients with several cancers including bladder, gastric and lung cancer [78]. These recommendations were partly based on preclinical studies suggesting that oestrogen stimulates growth in bladder [9] and gastric cancer cell lines [10] and lung cancer mouse models [11]. Further, observational studies have shown increases in the risk of glioma and meningioma with use of oestrogen alone HRT [12] and the Women’s Health Initiative randomised controlled trial (RCT) showed a marked increase in death from lung cancer in the oestrogen plus progestin group compared with placebo [13]. Also, an early cohort study showed reduced survival with prior HRT use in patients with lung cancer [14] but this was not replicated in two later studies [1516]. Few epidemiological studies have investigated HRT use after cancer diagnosis and cancer-specific mortality, and to the best of our knowledge none has investigated HRT use after a bladder, gastric or lung cancer diagnosis.

We determined the association between HRT use after cancer diagnosis and the risk of cancer-specific mortality in women with a range of 17 cancers (excluding breast cancer because HRT is contraindicated in patients with breast cancer), to help inform the decision to use HRT in women with cancer.

Methods

The main methods are described in previously published protocols [1718]. The study utilised data sources from QResearch (version 44, England) [19] the Scottish National Prescribing Information System (Scotland) [20] and the SAIL databank (Wales) [2122]. QResearch is a general practice database including over 1000 practices and over 10 million patients [19]. QResearch is linked at individual patient level to a range of data sources including hospital admissions data and national mortality records. The Scottish National Prescribing Information System dataset was utilised from 2009 as after this timepoint it is estimated to capture over 95% of prescriptions [20] and was linked to Scottish hospital admissions data which has been shown to be accurate for a wide range of conditions [23]. The SAIL Databank is a repository of health data which allows linkage of data from a number of sources including general practice and hospital admissions data with an accuracy of over 99% [2122]. Analyses in Scotland, England and Wales all utilised UK cancer registry data which have high levels of completeness and accuracy over the study period [24]. The study has been reported in accordance with the STROBE guidelines [25].

Cohorts

Population-based cohorts of women, aged 40–79, newly diagnosed with an incident cancer were identified solely from cancer registry records in Scotland (Scottish Cancer Registry) and Wales (Welsh Cancer Intelligence and Surveillance Unit) and from three data sources in England (general practice (GP) diagnosis codes, hospital diagnoses and cancer registry records from QResearch database). The dates of diagnosis included were: January 1998 to September 2019 in England; January 2009 to December 2016 in Scotland; and January 2000 to December 2016 in Wales. Seventeen of the most common female cancers (excluding breast cancer) were investigated (see Supplementary Table 1 for ICD10 codes used). Patients previously diagnosed with other invasive cancers (apart from non-melanoma skin cancer) were excluded. In a deviation from the published protocols, women with a diagnosis of thyroid cancer were excluded due to small numbers of diagnoses across all three datasets.

Exposure

HRT use was ascertained from electronic GP prescribing records (Wales and England) or dispensing records (Scotland) which were available from the date of cohort entry. The main HRT definition included systemic oestrogen-containing products (and tibolone) used for menopausal symptoms based upon the British National Formulary [4] classification (Section 6.8.1). Vaginal oestrogen therapy was also identified based upon the British National Formulary classification (contained in Section 7.6.2).

Outcome

The primary outcome was cancer-specific mortality from national mortality records (based upon the corresponding cancer as the underlying cause of death in England, Scotland and Wales, see Supplementary Table 1 for ICD10 codes used) for each of the 17 cancer sites. Linked national mortality records were available until March 2020 in England, December 2020 in Scotland, and June 2020 in Wales. We investigated the cancer-specific hazard because our primary interest was in the aetiological effect of HRT on cancer-specific mortality in those who were event free. A secondary analysis was conducted on all-cause death which increased power and avoided any potential misclassification of the cause of death.

Covariates

Cancer treatment (including radiotherapy, chemotherapy and surgery) was determined from cancer registry records in Scotland and Wales and from Hospital Episode Statistics (HES) in England. Cancer stage grouping was determined from cancer registry records except, to minimise missing data, Duke’s stage was used for colorectal cancer in Scotland, Breslow thickness was used for melanoma in Scotland, Figo stage was used for cervical cancer in Scotland and Wales and Figo stage was used for ovarian cancer in Scotland. Charlson comorbidities recorded before cancer diagnosis were identified from GP records and hospital admissions (available from 2000) in Wales, GP records and HES (available from 1998) in England and from hospital admissions alone (available from 1999) in Scotland (apart from diabetes which was identified from diabetes medications in Scotland). Other medication use (including aspirin, statins, metformin and oral contraceptives) was determined at any time before cancer diagnosis from GP prescribing (in England from 1989 and in Wales from 2000) or dispensing records (in Scotland from 2009). Hysterectomy/oophorectomy was determined from hospital admissions data in Scotland and hospital admissions and GP records in Wales and England. Deprivation of home address postcode was determined based upon the relevant 2011 Index of Multiple Deprivation in Scotland and Wales [2021] and the 2011 Townsend deprivation score in England [26]. Smoking and BMI was determined from the most recent GP records before cancer diagnosis (not available in Scotland).

Statistical analysis

In the primary analysis of HRT use after diagnosis, patients were followed from 6 months after cancer diagnosis to cancer-specific mortality and censored on death from other causes, end of follow-up (the latest date at which mortality records were complete) and additionally in England and Wales end of GP records. Consequently, patients who died in the first 6 months after cancer diagnosis were excluded as it seemed unlikely that HRT use after diagnosis could impact these deaths (in sensitivity analyses, this period was increased to 1 year). HRT was modelled as a time varying covariate to avoid immortal time bias [27], i.e. patients were first considered non-users and then users after a lag of 6 months following their first prescribed/dispensed HRT. A lag is recommended in studies of medication use and cancer survival [28]. A time-varying duration-response analysis was conducted with individuals considered a non-user prior to 6 months after first prescribed/dispensed HRT, a short term user from 6 months after first prescribed/dispensed HRT to 6 months after their fifth prescribed/dispensed HRT, and a longer term user after this time. The fifth prescribed/dispensed HRT was used because in a preliminary analysis this roughly corresponded to 1 year of HRT prescriptions. Time-dependent Cox regression models were used to calculate hazard ratios (HRs), and 95% confidence intervals (CIs), comparing users of HRT with non-users after cancer diagnosis adjusting for age at diagnosis, year of diagnosis, deprivation, Charlson comorbidity (before diagnosis), anaemia (before diagnosis) other medication use (including statins, aspirin, metformin and oral contraceptives before diagnosis), cancer treatment (including surgery, chemotherapy and radiotherapy), and hysterectomy/oophorectomy (in the periods up to 1 month before cancer diagnosis and 1 month before to 6 months after cancer diagnosis). The Cox proportional hazards assumption was checked by visual inspection of log(−log) plots and appeared to be largely satisfied. Estimates were calculated within each cohort and then pooled using random effects meta-analysis models [29].

Analyses were repeated for all-cause mortality. An additional analysis was conducted at all sites, additionally adjusting for BMI (at diagnosis, based upon complete case), from GP records before diagnosis, restricted to England and Wales. Additional analyses were conducted for the following more common cancers: colorectal, melanoma, ovarian, lung and cervical cancer. Analyses were conducted with type of HRT coded as a single time-varying variable, with a lag of 6 months, into the following hierarchical categories: combined HRT (with or without other HRT), oestrogen-only HRT (with or without tibolone) and tibolone only. HRT type was compared with no HRT use. A number of sensitivity analyses were conducted. First, the lag was increased to 1 year (with follow-up starting at one year after cancer diagnosis). Second, analyses were repeated with vaginal oestrogen therapy (mainly oestradiol pessaries and oestriol creams) included within the HRT definition. Third, systemic users of HRT were compared with users of vaginal oestrogen therapy, as users of vaginal oestrogen therapy are likely to share indications but receive much lower amounts of oestrogen than systemic users. In this analysis HRT was included using a single time-varying covariate (lagged by 6 months) coded as systemic HRT use (with or without vaginal oestrogen therapy use), vaginal oestrogen therapy use and no HRT use. Fourth, analyses was conducted restricting to stage 1 to 3 disease and restricting to stage 1 and 2 disease. Fifth, analyses were conducted varying the age range: restricting to women aged over 55 (who are more likely to be post-menopausal and widening to women aged 18–79 years. Sixth, a new user analysis was conducted restricted to women who had not used HRT in the period 18 to 6 months before diagnosis. Seventh, an analysis was conducted adjusting for cancer stage using multiple imputation. Stage was imputed in 10 imputed datasets using ordinal logistic regression models with cancer-specific death status, cumulative hazard along with all confounders from the adjusted model included in imputation models [30], and results were combined using Rubin’s rules [31]. Eighth, analyses were conducted additionally adjusting for cancer stage (based upon complete case) and additionally adjusting for cancer stage and smoking (based upon complete case, restricted to England and Wales). Finally, an analysis was conducted for death from cardiovascular disease (based upon ICD10 codes I20 to I99 or G45 as the underlying cause of death).

A separate analysis was conducted to investigate HRT use before diagnosis. In this analysis patients were followed from the date of cancer diagnosis to cancer-specific mortality (censored as previously) and HRT use was defined as one or more prescribed/dispensed systemic HRT in the period 18 months to 6 months before diagnosis. HRT use in the 6 months immediately before diagnosis was ignored because medication use can increase in this period due to cancer symptoms [32] which may be more marked in patients with advanced cancer. HRs (and 95% CIs) were calculated using Cox regression adjusting for age, year, deprivation, Charlson comorbidity (before diagnosis), other medication use (before diagnosis) and hysterectomy/oophorectomy (up to 1 month before cancer diagnosis). STATA 16/17 was used for all analyses. Analysis code is available from the authors upon request.

Ethical approval

Ethical approval for English data was obtained from the QResearch scientific committee (Ref: OX24, project title ‘Use of hormone replacement therapy and survival from cancer’). Ethical approval for the QResearch database is obtained annually from East Midlands – Derby Research Ethics Committee (Ref:18/EM/0400). Approval for analysis of the Welsh data has been obtained from the SAIL Databank Information Governance Review Panel (Reference: 0965) and approval for the analysis of the Scottish data has been obtained from the Public Benefit and Privacy Panel for Health and Social Care (Reference: 2021-0014).

Results

The final cohort contained 182,589 patients with cancer, who survived more than 6 months after their cancer diagnosis, across 17 cancer sites followed for 840,133 person years. There were 54,861 cancer-specific deaths during follow-up. Overall, 7% (11,972) of patients with cancer used systemic HRT after cancer diagnosis. For instance, 5% of patients with colorectal cancer, 4% of patients with lung cancer and 11% of patients with malignant melanoma used systemic HRT after cancer diagnosis.

Characteristic of users and non-users of HRT

The characteristics of users and non-users of HRT are shown in Table 1 and Supplementary Table 2. In general, users of HRT, compared with non-users, were younger at diagnosis, had higher rates of hysterectomy/oopherectomy, and smoking. Also, a lower proportion of users of HRT had diabetes or chronic kidney disease and had been prescribed statins, aspirin or metformin and a higher proportion had been prescribed oral contraceptives. A lower proportion of users of HRT received chemotherapy than HRT non-users. The distribution of stage was fairly similar for women with colorectal and lung cancer in users of HRT compared with HRT non-users, but a greater proportion of HRT users had stage 1 and 2 disease for women with cervical and ovarian cancer. Other characteristics of users of HRT and non-users were largely similar.

Table 1 Patient characteristics by hormone replacement therapy use after cancer diagnosis in women with at least 6 months of follow-up.

Full size table

HRT use after cancer diagnosis and cancer-specific mortality

The pooled associations between HRT use after diagnosis and cancer-specific mortality are shown in Table 2 and Fig. 1. There was no evidence of higher cancer-specific mortality in users of HRT after diagnosis, compared with non-users, at any of the 17 cancer sites studied. In contrast, use of HRT compared with non-use was associated with a lower rate of cancer-specific mortality for colorectal cancer (adjusted HR = 0.79 95% CI 0.70, 0.90), ovarian cancer (adjusted HR = 0.60 95% CI 0.39, 0.93), uterus cancer (adjusted HR = 0.43 95% CI 0.27, 0.67), kidney cancer (adjusted HR = 0.55 95% CI 0.40, 0.76), oral cancer (adjusted HR = 0.58 95% CI 0.42, 0.80) and non-Hodgkin lymphoma (adjusted HR = 0.77 95% CI 0.60, 0.99). However, the analysis of uterus and kidney was based upon relatively small numbers of cancer-specific deaths in users of HRT (less than 40) and in further analysis of oral cancer and non-Hodgkin lymphoma the association did not follow a dose-response as there was no association in patients with 5 or more prescriptions compared with no prescriptions (adjusted HR = 0.87 95% CI 0.57, 1.32 and adjusted HR = 0.89 95% CI 0.65, 1.21, respectively). Associations were generally similar in analyses additionally adjusting for BMI (in England and Wales) shown in Supplementary Table 3.

Table 2 Pooled analyses of hormone replacement therapy use after diagnosis and cancer-specific mortality in England, Scotland and Wales.

Full size table

Fig. 1

pastedGraphic_3.png

Adjusted hazard ratios for the association between hormone replacement therapy use after diagnosis and cancer-specific mortality in England (pastedGraphic_4.png), Scotland (pastedGraphic_5.png), Wales (pastedGraphic_6.png) and pooled (pastedGraphic_7.png), by site.

Full size image

Additional analyses conducted for colorectal, lung, melanoma, cervical and ovarian cancer are shown in Table 3 and Supplementary Table 4. The inverse association between HRT and cancer-specific mortality in patients with colorectal cancer was similar in most sensitivity analyses. However, the association was attenuated when systemic users of HRT were compared with users of vaginal oestrogen therapy (adjusted HR = 0.93 95% CI 0.78, 1.10) and in the analysis restricted to stage 1–3 patients with colorectal cancer (adjusted HR = 0.87 95% CI 0.68, 1.11). The null association between HRT and cancer-specific mortality in patients with lung cancer (adjusted HR = 0.98 95% CI 0.90, 1.07) was generally consistent across sensitivity analyses, except for the analysis comparing systemic to vaginal oestrogen therapy users in which a slight increase was observed (adjusted HR = 1.18 95% CI 1.01, 1.39). The null association between HRT and cancer-specific mortality in patients with melanoma and cervical cancer was also fairly similar across sensitivity analyses. The inverse association between HRT use and cancer-specific mortality in patients with ovarian cancer was similar in most sensitivity analyses. Further, analyses by HRT type are shown in Supplementary Table 5. There were no marked differences in associations by HRT type for patients with colorectal cancer, lung cancer or melanoma. In patients with cervical cancer and ovarian cancer an inverse association was observed solely in patients on oestrogen-only HRT compared with no HRT use (adjusted HR = 0.65 95% CI 0.47, 0.91 and adjusted HR = 0.55 95% CI 0.38, 0.79, respectively).

Table 3 Pooled sensitivity analyses of hormone replacement therapy after diagnosis and cancer-specific mortality in England, Scotland and Wales.

Full size table

HRT use after cancer diagnosis and all-cause mortality

Analysis of the association between HRT use after diagnosis and all-cause mortality, shown in Supplementary Table 6, did not reveal any evidence of an increase in all-cause mortality in patients with cancer using HRT at any of the 17 sites studied.

HRT use before diagnosis and cancer-specific mortality

The pooled association between HRT use before diagnosis and cancer-specific mortality is shown in Table 4 and Supplementary Figure 1. There was no evidence that HRT before diagnosis was associated with higher cancer-specific mortality in patients with cancer at any of the 17 sites studied. There were inverse associations between HRT use, compared with HRT non-use, before diagnosis and cancer-specific mortality for patients with colorectal cancer (adjusted HR = 0.86 95% CI 0.77, 0.97), cervical cancer (adjusted HR = 0.69 95% CI 0.49, 0.98), oral cancer (adjusted HR = 0.72 95% CI 0.55, 0.95) and non-Hodgkin’s lymphoma (adjusted HR = 0.79, 0.66, 0.95).

Table 4 Pooled analysis of hormone replacement therapy use before diagnosis and cancer-specific mortality in England, Scotland and Wales.

Full size table

Discussion

Overall, there was no evidence that patients with any of the 17 cancers studied who took HRT following their cancer diagnosis had higher rates of cancer-specific or all-cause mortality. Use of HRT was associated with reductions in cancer-specific mortality in women with colorectal, ovarian, uterus, kidney, oral and non-Hodgkin lymphoma, but these associations were based upon relatively small numbers or were generally not consistent across sensitivity analyses.

There has been limited previous research on the safety of HRT use after diagnosis in cancer patients. Small randomised controlled trials have been conducted investigating HRT use and survival in patients with ovarian cancer (showing reduced mortality but not progression free survival) [33], endometrial cancer (showing no association) [34] and breast cancer [35], but not at other cancer sites. Observational studies have investigated HRT and survival for patients with ovarian, endometrial, colorectal, melanoma and lung cancer. These observational studies showed a reduced risk of mortality in users of HRT with ovarian cancer [36], a reduced cancer recurrence in users of HRT with endometrial cancer [37] and a reduced cancer-specific mortality with current use of HRT in patients with colorectal cancer [38]. Mixed associations were observed between HRT and survival in previous smaller studies of patients with melanoma [3940] and lung cancer [14,15,16]. To our knowledge, observational studies have not been conducted investigating HRT use after diagnosis and survival for the other cancers sites studied.

Previous reviews of the oncologic safety of HRT have recommended, based upon preclinical and other evidence, that patients with bladder [7], lung [78], brain [7] and gastric cancer [7] avoid HRT. Our study does not provide evidence of increased cancer-specific mortality in users of HRT with bladder, brain or gastric cancer. In most analyses of lung cancer there was no evidence of association but in a sensitivity analysis comparing users of systemic HRT with users of vaginal oestrogen therapy, there was a slight association with increased cancer-specific mortality. Consequently, further research on HRT use in patients with lung cancer is merited.

Our study has several strengths and limitations. The study utilised data from three independent population-based data sources containing over 180,000 patients with cancer with follow-up of up to 21 years. The safety of HRT after diagnosis has not been previously investigated at many of the cancer sites studied. The use of prescribing/dispensing records will have eliminated recall bias and should capture all HRT use because, at the time of the study, HRT was only available by prescription in the UK. However, these data sources do not contain information on actual adherence to HRT.

The main weakness of our study is that HRT was not randomly allocated and hence HRT users may differ from HRT non-users in ways which influence cancer-specific mortality resulting in confounding. We accounted for potential confounding by adjusting for a wide range of confounders, but we cannot rule out residual confounding from unavailable variables (such as parity, age at menopause or alcohol consumption) or incomplete variables (such as smoking). Further, we did not have comprehensive information on contraindications for HRT (e.g. abnormal liver functions tests) and could not reliably adjust for these. Also, across the three countries there will be differential capture of some confounders due to the sources used (for instance, rates of cancer treatments were different between countries). It is possible that users of HRT have healthier lifestyles in general [41] and may have lower body mass index and be more physically active both of which have been shown to be associated with better cancer survival [4243]. It is also possible that users of HRT will have more contact with their GP and be more likely to attend screening and other diagnostic examinations also leading to improved outcomes. Family history of cancer may also impact on the decision to use HRT and cancer-specific mortality. Patients with cancer who have a better prognosis may be more likely to receive HRT, because such patients may be more concerned about their quality of life and clinicians may be more inclined to prescribe HRT to them, leading to artificially better survival in patients on HRT. There was no association with increased cancer-specific mortality when investigating HRT before diagnosis (in which this latter bias would not occur). In cancer sites where stage was available, there was no evidence of an association of increased cancer specific mortality with HRT use after diagnosis following adjustment for stage, and a range of potential confounders. We cannot rule out residual confounding by stage for cancer sites where stage was not complete or confounding by stage for cancer sites where stage was not available. The analyses based upon the multiple imputation of stage relies upon the assumption that stage data are missing at random, which we cannot test, and which could be violated if patients with missing data were more likely to have worse stage even after adjusting for variables in the imputation models. We conducted an active comparator analysis comparing systemic with vaginal oestrogen therapy users (who are likely to share many indications and risk factors but will have markedly lower exposure to oestrogen) to attempt to reduce confounding by indication [44]. In our analysis some women using combined HRT consisting of a prescription for oestrogen-only HRT and a separately prescribed progestogen (such as the levonorgestrel-releasing intrauterine system) may have been misclassified as oestrogen-only users of HRT. There remains the possibility of Type 2 error particularly at rarer cancer sites and sites with limited HRT use.

Importantly, there is a particular risk of confounding by indication in women with oestrogen-sensitive cancers such as uterus, ovarian and cervical cancer. This bias could incorrectly lead to null or even inverse associations if women perceived to be at lower risk of recurrence are more likely to receive HRT. Consequently our findings for these cancer sites should be interpreted particularly cautiously. Also, we could not investigate rarer or specific subtypes of certain cancers which may be particularly oestrogen-sensitive, because of small numbers and/or lack of data, and our results cannot be extrapolated to these groups; for instance oestrogen receptor positive gastric [7], oestrogen receptor positive bladder cancer [7], endometrial stromal sarcoma [45], granulosa cell ovarian cancer [45], low-grade serous ovarian [45] and cervical adenocarcinoma [45].

Many limitations of our study reflect the use of routinely collected data which does not contain sufficient detail on tumour type or covariates. There is a need for prospective cohort studies which can capture detailed information on specific potential confounders such as family history, physical activity and alcohol intake as well as ensuring comprehensive data on tumour characteristics including histological classification and stage.

Our study may provide some reassurance to clinicians and patients of the safety of systemic HRT in women with one of the 17 cancers studied, but as stated above should be interpreted cautiously in women with oestrogen-sensitive cancers. Along with other known risks associated with HRT use [4647], our findings may contribute to the decision of cancer patients, and their prescribers, to use HRT.

In conclusion, in this large observational study we observed little consistent evidence of an association between HRT use and increased cancer-specific mortality in women with any of the 17 cancers studied.

Data availability

The datasets from Scotland (eDRIS; Public Health Scotland, https://www.isdscotland.org/products-and-services/edris/), Wales (SAIL databank; Swansea University, https://saildatabank.com/) and England (QResearch; University of Oxford, https://www.qresearch.org/) were obtained under strict data access conditions that allowed the study to be conducted but do not allow direct data sharing. However, the data analysed in this study would in principle be available to a researcher who applied to the data custodians and obtained the same approvals.

References
  1. Baber RJ, Panay N, Fenton A. 2016 IMS Recommendations on womens midlife health and menopause hormone therapy. Climacteric. 2016;19:109–50.

Article CAS PubMed Google Scholar 

  1. Goodman NF, Cobin RH, Ginzburg SB, Katz IA, Woode DE, Endocrinologists AAoC. American Association of Clinical Endocrinologists Medical Guidelines for Clinical Practice for the diagnosis and treatment of menopause. Endocr Pr. 2011;17:1–25.

Article Google Scholar 

  1. Cobin RH, Goodman NF. American Association of Clinical Endocrinologists and American College of Endocrinology Position Statement on Menopause–2017 Update. Endocr Pr. 2017;23:869–80.

Article Google Scholar 

  1. Joint Formulary Committee. British National Formulary 80: September 2020-March 2021. (BMJ Group and Pharmaceutical Press: London, UK, 2020).
  2. Manson JE, Kaunitz AM. Menopause management — getting clinical care back on track. N. Engl J Med. 2016;374:803–6.

Article PubMed Google Scholar 

  1. Newson LR. Best practice for HRT: Unpicking the evidence. Br J Gen Pr. 2016;66:597–8.

Article Google Scholar 

  1. Deli T, Orosz M, Jakab A. Hormone replacement therapy in cancer survivors – review of the literature. Pathol Oncol Res. 2020;26:63–78.

Article CAS PubMed Google Scholar 

  1. Kuhle CL, Kapoor E, Sood R, Thielen JM, Jatoi A, Faubion SS. Menopausal hormone therapy in cancer survivors: a narrative review of the literature. Maturitas. 2016;92:86–96.

Article CAS PubMed Google Scholar 

  1. Shen SS, Smith CL, Hsieh JT, Yu J, Kim I, Jian W, et al. Expression of estrogen receptors-α and -β in bladder cancer cell lines and human bladder tumor tissue. Cancer. 2006;106:2610–6.

Article CAS PubMed Google Scholar 

  1. 10.Wang X, Chen Q, Huang X, Zou F, Fu Z, Chen Y, et al. Effects of 17β-estradiol and tamoxifen on gastric cancer cell proliferation and apoptosis and ER-α36 expression. Oncol Lett. 2017;13:57–62.

Article PubMed Google Scholar 

  1. 11.Hammoud Z, Tan B, Badve S, Bigsby RM. Estrogen promotes tumor progression in a genetically defined mouse model of lung adenocarcinoma. Endocr Relat Cancer. 2008;15:475–83.

Article CAS PubMed Google Scholar 

  1. 12.Benson VS, Kirichek O, Beral V, Green J. Menopausal hormone therapy and central nervous system tumor risk: Large UK prospective study and meta-analysis. Int J Cancer. 2015;136:2369–77.

Article CAS PubMed Google Scholar 

  1. 13.Chlebowski RT, Schwartz AG, Wakelee H, Anderson GL, Stefanick ML, Manson JE. Oestrogen plus progestin and lung cancer in postmenopausal women (Women’s Health Initiative trial): a post-hoc analysis of a randomised controlled trial. Lancet. 2009;374:1243–51.

Article CAS PubMed PubMed Central Google Scholar 

  1. 14.Ganti AK, Sahmoun AE, Panwalkar AW, Tendulkar KK, Potti A. Hormone replacement therapy is associated with decreased survival in women with lung cancer. J Clin Oncol. 2006;24:59–63.

Article CAS PubMed Google Scholar 

  1. 15.Clague J, Reynolds P, Henderson KD, Sullivan-Halley J, Ma H, Lacey JV, et al. Menopausal hormone therapy and lung cancer-specific mortality following diagnosis: The California teachers study. PLoS One. 2014;9:e103735.

Article PubMed PubMed Central Google Scholar 

  1. 16.Huang B, Carloss H, Wyatt SW, Riley E. Hormone replacement therapy and survival in lung cancer in postmenopausal women in a rural population. Cancer. 2009;115:4167–75.

Article PubMed Google Scholar 

  1. 17.McMenamin U, Hicks B, Hughes C, Murchie P, Hippisley-Cox J, Ranger T, et al. Hormone replacement therapy in women with cancer and risk of cancer-specific mortality and cardiovascular disease: a protocol for a cohort study from Scotland and Wales. BMC Cancer. 2021;21:313.

Article CAS PubMed PubMed Central Google Scholar 

  1. 18.Ranger TA, Burchardt J, Clift AK, Mei W, Coupland C, Tan PS, et al. Hormone replacement therapy and cancer survival: a longitudinal cohort study: protocol paper. BMJ Open. 2021;11:e046701.

Article PubMed PubMed Central Google Scholar 

  1. 19.Hippisley-Cox J, Stables D, Pringle M. QRESEARCH: a new general practice database for research. Inf Prim Care. 2004;12:49–50.

Google Scholar 

  1. 20.Alvarez-Madrazo S, McTaggart S, Nangle C, Nicholson E, Bennie M. Data resource profile: the Scottish national prescribing information system (PIS). Int J Epidemiol. 2016;45:714–715.

Article PubMed PubMed Central Google Scholar 

  1. 21.Ford DV, Jones KH, Verplancke JP, Lyons RA, John G, Brown G, et al. The SAIL Databank: building a national architecture for e-health research and evaluation. BMC Health Serv Res. 2009;9:157.

Article PubMed PubMed Central Google Scholar 

  1. 22.Lyons RA, Jones KH, John G, Brooks CJ, Verplancke JP, Ford DV, et al. The SAIL databank: linking multiple health and social care datasets. BMC Med Inf Decis Mak. 2009;9:3.

Article Google Scholar 

  1. 23.Information Services Division (ISD) Scotland: Assessment of SMR01 Data 2010-2011. Available from: https://webarchive.nrscotland.gov.uk/20180620031846/https://www.isdscotland.org/Health-Topics/Hospital-Care/Publications/2012-05-08/2012-05-08-SMR01-Summary.pdf?27016848326 (Accessed Feb 2024).
  2. 24.The United Kingdom and Ireland Association of Cancer Registries: Performance Indicators. 2019. Available from: https://www.ukiacr.org/kpis (Accessed Feb 2024).
  3. 25.von Elm E, Altman DG, Egger M, Pocock SJ, Gotzsche PC, Vandenbroucke JP, et al. The strengthening the reporting of observational studies in epidemiology (STROBE) statement: guidelines for reporting observational studies. PLoS Med. 2007;4:e296.

Article Google Scholar 

  1. 26.Hippisley-Cox J, Coupland C. Development and validation of QMortality risk prediction algorithm to estimate short term risk of death and assess frailty: cohort study. BMJ. 2017;358:j4208.

Article PubMed PubMed Central Google Scholar 

  1. 27.Levesque LE, Hanley JA, Kezouh A, Suissa S. Problem of immortal time bias in cohort studies: example using statins for preventing progression of diabetes. BMJ. 2010;340:b5087.

Article PubMed Google Scholar 

  1. 28.Chubak J, Boudreau DM, Wirtz HS, McKnight B, Weiss NS. Threats to validity of nonrandomized studies of postdiagnosis exposures on cancer recurrence and survival. J Natl Cancer Inst. 2013;105:1456–62.

Article PubMed PubMed Central Google Scholar 

  1. 29.Stukel TA, Demidenko E, Dykes J, Karagas MR. Two-stage methods for the analysis of pooled data. Stat Med. 2001;20:2115–30.

Article CAS PubMed Google Scholar 

  1. 30.White IR, Royston P. Imputing missing covariate values for the Cox model. Stat Med. 2009;28:1982–98.

Article PubMed PubMed Central Google Scholar 

  1. 31.Sterne JA, White IR, Carlin JB, Spratt M, Royston P, Kenward MG, et al. Multiple imputation for missing data in epidemiological and clinical research: potential and pitfalls. BMJ. 2009;338:b2393.

Article PubMed PubMed Central Google Scholar 

  1. 32.Pottegard A, Hallas J. New use of prescription drugs prior to a cancer diagnosis. Pharmacoepidemiol Drug Saf. 2017;26:223–7.

Article PubMed Google Scholar 

  1. 33.Saeaib N, Peeyananjarassri K, Liabsuetrakul T, Buhachat R, Myriokefalitaki E. Hormone replacement therapy after surgery for epithelial ovarian cancer. Cochrane Database Syst Rev. 2020;1:CD012559.

PubMed Google Scholar 

  1. 34.Barakat RR, Bundy BN, Spirtos NM, Bell J, Mannel RS, Gynecologic Oncology Group S. Randomized double-blind trial of estrogen replacement therapy versus placebo in stage I or II endometrial cancer: a Gynecologic Oncology Group Study. J Clin Oncol. 2006;24:587–92.

Article CAS PubMed Google Scholar 

  1. 35.Poggio F, Del Mastro L, Bruzzone M, Ceppi M, Razeti MG, Fregatti P, et al. Safety of systemic hormone replacement therapy in breast cancer survivors: a systematic review and meta-analysis. Breast Cancer Res Treat. 2022;191:269–75.

Article CAS PubMed Google Scholar 

  1. 36.Li D, Ding CY, Qiu LH. Postoperative hormone replacement therapy for epithelial ovarian cancer patients: a systematic review and meta-analysis. Gynecol Oncol. 2015;139:355–62.

Article CAS PubMed Google Scholar 

  1. 37.Londero AP, Parisi N, Tassi A, Bertozzi S, Cagnacci A. Hormone replacement therapy in endometrial cancer survivors: A meta-analysis. J Clin Med. 2021;10:3165.

Article CAS PubMed PubMed Central Google Scholar 

  1. 38.Jang YC, Huang HL, Leung CY. Association of hormone replacement therapy with mortality in colorectal cancer survivor: a systematic review and meta-analysis. BMC Cancer. 2019;19:1199.

Article CAS PubMed PubMed Central Google Scholar 

  1. 39.Hicks BM, Kristensen KB, Pedersen SA, Holmich LR, Pottegard A. Hormone replacement therapy and the risk of melanoma in post-menopausal women. Hum Reprod. 2019;34:2418–29.

Article CAS PubMed Google Scholar 

  1. 40.MacKie RM, Bray CA. Hormone replacement therapy after surgery for stage 1 or 2 cutaneous melanoma. Br J Cancer. 2004;90:770–2.

Article CAS PubMed PubMed Central Google Scholar 

  1. 41.Gleason CE, Dowling NM, Friedman E, Wharton W, Asthana S. Using predictors of hormone therapy use to model the healthy user bias: how does healthy user status influence cognitive effects of hormone therapy? Menopause. 2012;19:524–33.

Article PubMed PubMed Central Google Scholar 

  1. 42.Friedenreich CM, Stone CR, Cheung WY, Hayes SC. Physical activity and mortality in cancer survivors: a systematic review and meta-Analysis. JNCI Cancer Spectr. 2020;4:pkz080.

Article PubMed Google Scholar 

  1. 43.Petrelli F, Cortellini A, Indini A, Tomasello G, Ghidini M, Nigro O, et al. Association of obesity with survival outcomes in patients with cancer: A systematic review and meta-analysis. JAMA Netw Open. 2021;4:e213520.

Article PubMed PubMed Central Google Scholar 

  1. 44.Lund JL, Richardson DB, Sturmer T. The active comparator, new user study design in pharmacoepidemiology: historical foundations and contemporary application. Curr Epidemiol Rep. 2015;2:221–8.

Article PubMed PubMed Central Google Scholar 

  1. 45.Harris BS, Bishop KC, Kuller JA, Ford AC, Muasher LC, Cantrell SE, et al. Hormonal management of menopausal symptoms in women with a history of gynecologic malignancy. Menopause. 2020;27:243–8.

Article PubMed Google Scholar 

  1. 46.Vinogradova Y, Coupland C, Hippisley-Cox J. Use of hormone replacement therapy and risk of venous thromboembolism: Nested case-control studies using the QResearch and CPRD databases. BMJ. 2019;364:k4810.

Article PubMed PubMed Central Google Scholar 

  1. 47.Vinogradova Y, Coupland C, Hippisley-Cox J. Use of hormone replacement therapy and risk of breast cancer: nested case-control studies using the QResearch and CPRD databases. BMJ. 2020;371:m3873.

Article PubMed PubMed Central Google Scholar 

Download references

Acknowledgements

We would like to acknowledge the support of the eDRIS team (Public Health Scotland) for their involvement in obtaining approvals, provisioning and linking data and the secure analytical platform within the National Safe Haven. We would also like to acknowledge support of SAIL Databank for facilitating access to the dataset from Wales. We acknowledge the contribution of EMIS practices who contribute to the QResearch database and the Chancellor, Masters and Scholars of the University of Oxford for continuing to develop and support the QResearch database. The Hospital Episode Statistics data used in the English portion of this analysis are re-used by permission from NHS Digital who retain the copyright. We thank the Office for National Statistics (ONS) for providing the mortality data for the English analyses. The ONS bears no responsibility for the analysis or interpretation of the data. The authors would also like to thank the PPI representatives for providing a patient and public perspective on the study design, findings, interpretation of the study and lay summary materials.

Funding

This work was supported by Cancer Research UK (reference C37316/A29656). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Author information

Authors and Affiliations

  1. Centre for Public Health, Queen’s University Belfast, Belfast, Northern Ireland, UK

Chris R. Cardwell, Alexander M. Labeit, Blánaid Hicks & Úna McMenamin

  1. Nuffield Department of Primary Care Health Sciences, University of Oxford, Oxford, UK

Tom A. Ranger, Xue W. Mei & Julia Hippisley-Cox

  1. Centre for Academic Primary Care, University of Nottingham, Nottingham, UK

Carol A. C. Coupland

  1. School of Pharmacy, Queen’s University Belfast, Belfast, Northern Ireland, UK

Carmel Hughes

  1. Division of Applied Health Sciences Section, Academic Primary Care, University of Aberdeen, Foresterhill, Aberdeen, UK

Peter Murchie

Contributions

Funding acquisition: CRC, CACC, BH, CH, UMcM, PM, JHC; Data acquisition: CRC, CACC, BH, CH, UMcM, PM, JHC; Study design: CRC, CACC, BH, CH, UMcM, PM, JHC; Data analysis: CRC, TAR, AML, CACC, XWM; Data interpretation: All authors; Writing original draft: CRC; Writing, review and editing: All authors; All authors approved the final version of the manuscript.

Corresponding author

Correspondence to Chris R. Cardwell.

Ethics declarations

Competing interests

JHC is an unpaid director of QResearch, a not-for-profit organisation which is a partnership between the University of Oxford and EMIS Health who supply the QResearch database used for this work. JHC has a 50% shareholding in ClinRisk Ltd, co-owning it with her husband, who is a director. As a shareholder and spouse of a director she has a financial and family interest in the ongoing and future success of the company. The company licences software both to the private sector and to NHS bodies or bodies that provide services to the NHS (through GP electronic health record providers, pharmacies, hospital providers and other NHS providers). This software implements algorithms (outside the scope of this research) developed from access to the QResearch database during her time at the University of Nottingham. The other authors have declared no competing interests.

Ethics approval and consent to participate

Approval for analysis of the Welsh data has been obtained from the SAIL Databank Information Governance Review Panel (Reference: 0965) and approval for the analysis of the Scottish data has been obtained from the Public Benefit and Privacy Panel for Health and Social Care (Reference: 2021–0014). In both Scotland and Wales only de-identified data was used and therefore written informed consent was not obtained. Approval for the analysis of the English data was obtained from the QResearch scientific committee (Project ref: OX24, title ‘Use of hormone replacement therapy and survival from cancer’). Ethical approval for the QResearch database is obtained annually from East Midlands—Derby Research Ethics Committee (Ref: 18/EM/0400).

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

 

Copyright © 2024 Springer Nature Limited. All rights reserved.