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

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 oestrogendependent 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 [5,6].
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 [7,8].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 [15,16].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 [17,18].The study utilised data sources from QResearch (version 44, England) [19] the Scottish National Prescribing Information System (Scotland) [20] and the SAIL databank (Wales) [21,22].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% [21,22].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 1 for ICD10 codes used).Patients previously diagnosed with other invasive cancers (apart from nonmelanoma 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 oestrogencontaining 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 [20,21] 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 timevarying 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 timevarying 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.

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 nonuse 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.
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 cancerspecific 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).

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).

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 [39,40] 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 [7,8], 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  Stage imputed using multiple imputation as described in methods.Adjusted model contains all terms in a along with stage.f Restricted to patients with available stage.Adjusting for stage using complete case, model contains all terms in a along with stage. 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 [42,43].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 levonorgestrelreleasing 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], lowgrade 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 [46,47], 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.

d
Restricted to patients stage 1-3.Adjusted model contains all terms in a along with stage.e -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

Table 1 .
Patient characteristics by hormone replacement therapy use after cancer diagnosis in women with at least 6 months of follow-up.
a Hysterectomy/oophorectomy in the following time periods: before cancer (up to 1 month before cancer diagnosis), at cancer diagnosis (from 1 month before cancer diagnosis to 6 months after cancer diagnosis) and after cancer diagnosis (more than 6 months after cancer diagnosis).b Range shown to maintain statistical disclosure control.c BMI available for 84,285 HRT non-users and 6685 HRT users in England and 19,161 HRT non-users and 1097 HRT users in Wales.C.R. Cardwell et al.

Table 2 .
Pooled analyses of hormone replacement therapy use after diagnosis and cancer-specific mortality in England, Scotland and Wales.
a Adjusted model contains age, year of diagnosis, deprivation, cancer treatment (surgery, radiotherapy, chemotherapy), Charlson comorbidities (before diagnosis), anaemia (before diagnosis), medication use (before diagnosis: statin, aspirin, metformin and oral contraceptive) and hysterectomy/oophorectomy (before or at diagnosis).b P value from adjusted Cox regression model.c Non-Hodgkin lymphoma.d Range shown to maintain statistical disclosure control.

Table 3 .
Pooled sensitivity analyses of hormone replacement therapy after diagnosis and cancer-specific mortality in England, Scotland and Wales.
bIndividuals not using HRT or vaginal oestrogen therapy excluded.c Restricted to patients not using systemic HRT in the period 18 to 6 months before cancer diagnosis.

Table 4 .
Pooled analysis of hormone replacement therapy use before diagnosis and cancer-specific mortality in England, Scotland and Wales.
a Adjusted model contains age, year, deprivation, Charlson comorbidities (before diagnosis), medication use (before diagnosis: statin, aspirin, metformin, oral contraceptive) and hysterectomy/oophorectomy (before diagnosis).b P value from unadjusted Cox regression model.c P value from adjusted Cox regression model.C.R. Cardwell et al.