key: cord-0732108-z464aibq
authors: Gedeborg, Rolf; Styrke, Johan; Loeb, Stacy; Garmo, Hans; Stattin, Pär
title: Androgen deprivation therapy and excess mortality in men with prostate cancer during the initial phase of the COVID-19 pandemic
date: 2021-10-07
journal: PLoS One
DOI: 10.1371/journal.pone.0255966
sha: ca97648a36a556287b925c29726ee4b7e018f7db
doc_id: 732108
cord_uid: z464aibq

BACKGROUND: Men have a higher risk of death from COVID-19 than women and androgens facilitate entrance of the SARS-CoV-2 virus into respiratory epithelial cells. Thus, androgen deprivation therapy may reduce infection rates and improve outcomes for COVID-19. In the spring of 2020, Sweden was highly affected by COVID-19. The aim was to estimate the impact of androgen deprivation therapy on mortality from COVID-19 in men with prevalent prostate cancer by comparing all-cause mortality in the spring of 2020 to that in previous years. PATIENTS AND METHODS: Using the Prostate Cancer data Base Sweden all men with prostate cancer on March 1 each year in 2015–2020 were followed until June 30 the same year. Exposure to androgen deprivation therapy was ascertained from filled prescriptions for bicalutamide monotherapy, gonadotropin-releasing hormone agonists (GnRH), or bilateral orchidectomy. RESULTS: A total of 9,822 men died in March-June in the years 2015–2020, of whom 5,034 men were on androgen deprivation therapy. There was an excess mortality in 2020 vs previous years in all men. The crude relative mortality rate ratio for 2020 vs 2015–2019 was 0.93 (95% confidence interval (CI) 0.83 to 1.04) in men on GnRH, and 0.90 (95% CI 0.78 to 1.05) in men on bicalutamide monotherapy. After multivariable adjustment these ratios were attenuated to 1.00 (95% CI 0.89 to 1.12) and 0.97 (95% CI 0.84 to 1.12), respectively. When restricting the analysis to the regions with the highest incidence of COVID-19 or to the time period between 2 April to 10 June when mortality in 2020 was increased >30% compared to previous years, the results were similar to the main analysis. CONCLUSIONS: In this large national population-based cohort of men with prevalent prostate cancer, there was no clear evidence in support for an effect of androgen deprivation therapy on COVID-19 mortality.

a1111111111 a1111111111 a1111111111 a1111111111 a1111111111

Men have consistently been shown to have approximately 50% higher risk of COVID-19 compared to women, and the increase in risk remained after adjustment for risk factors [1, 2] . A similarly increased risk in men has also been observed for critical care admission and death in COVID-19 [2] [3] [4] [5] . These differences associated with sex have also been observed in children [6, 7] .

One proposed mechanism for the susceptibility of men to COVID-19 is that androgen signaling facilitates entry of the virus into host cells [8] . Internalization of the SARS-CoV-2 virus relies on proteases such as the Transmembrane Protease Serine 2 (TMPRSS2) surface protein expressed by epithelial cells (e.g in the lung). In vitro, inhibition of the TMPRSS2 protease activity resulted in decreased entry of SARS-CoV-2 into pulmonary epithelial cells [8] . The regulation of TMPRSS2 expression is androgen-dependent, suggesting that androgen deprivation therapy (ADT) could be protective against COVID-19 [9] . The response to gonadotropinreleasing hormone (GnRH) agonists would then from a mechanistic perspective be of primary interest. While this is an interesting hypothesis, further support from patient-relevant clinical outcomes is needed. An early study reported that men with prostate cancer not on ADT had a four-fold increased risk of COVID-19 and a more severe disease compared with men on ADT [10] . However, the results in this study were based on 118 men with COVID-19 of whom merely four were on ADT and 114 were not. Two subsequent small studies failed to replicate this association [11, 12] . Based on these findings, several clinical trials have been initiated to evaluate the benefit of ADT in men with COVID-19 [13] . However, in men with prostate cancer, use of ADT has been associated with adverse effects including an increased risk of cardiovascular disease [14] even within six months after initiation, and this increase has been particularly strong in men with prior cardiovascular disease [15] . Men hospitalized with COVID-19 often have cardiovascular risk factors [16] , but it is unknown if very short-term ADT for COVID-19 affects cardiovascular risk. There are also numerous other potential mechanisms that could lead to the observed male susceptibility to COVID-19 [17] .

We took advantage of a unique setting in which we had access to diagnosis, treatment and vital status up to July 2020 for virtually all men with prostate cancer in Sweden, with a source population of more than 10 million people that was strongly affected by the COVID-19 pandemic in the spring of 2020.

The aim of this study was to investigate if hormonal treatment for prostate cancer is associated with a reduction of mortality risk related to COVID-19.

Men with prostate cancer were identified in The National Prostate Cancer Register (NPCR) of Sweden, which since 1998 captures 98% of all cases of prostate cancer in Sweden as compared to the Swedish Cancer Register to which registration is mandatory. NPCR contains comprehensive data on cancer characteristics and primary treatment [18] .

Welfare (https://www.socialstyrelsen.se/en/) and Statistics Sweden (https://www.scb.se/en/) which are Government Agencies providing access to the linked healthcare registers. The data can, however, be retrieved after application made to any of the steering groups of NPCR and PCBaSe. For detailed information, please see www.npcr.se/in-english, where registration forms, manuals, and annual reports from NPCR are available alongside a full list of publications from PCBaSe. The statistical program code used for the present study analyses can be provided on request (contact hans. garmo@rccmellan.se).

By use of the Personal Identity Number unique for every Swedish citizen, data for men in NPCR diagnosed up to 31 Dec, 2019 was linked in July 2020 (i.e. with follow-up on vital status until 1 July 2020) to the Cancer Register [19] , statistics on deaths submitted by the Swedish Tax Agency to Statistics Sweden, the Prescribed Drug Register [20] , and the National Patient Register [21] held at The National Board of Health and Welfare [22] , to create the Prostate Cancer data Base Sweden (PCBaSe) RAPID 2019. The study was approved by the Ethical Review Authority. Since this was a study exclusively based on existing large-scale registry data the requirement for individual informed consent was waived by the Authority.

In PCBaSe an open cohort was defined with men with prevalent prostate cancer alive on 1 March 2020 and similar subsets for corresponding months in 2015-2019. Men from three regions in which GnRH agonists were provided directly from the hospital without prescription were excluded. The remaining 18 Swedish regions forming the source population for the study covered 91% of the Swedish population as of 31 December 2019. For complementary analyses, a more restricted time interval was defined from the first to the last date in 2020 with a mortality rate more than 30% above the average mortality rate in 2015-2019.

Separate analyses were also performed within the ten regions with the highest cumulative number of cases of COVID-19 as determined on 28 June 2020 from data published by the Public Health Agency of Sweden (S1 Fig).

Exposure to ADT was ascertained from fillings in The Prescribed Drug Register for bicalutamide monotherapy (Anatomical Therapeutic Chemical (ATC) code L02BB03 and also the rare use (~2%) of flutamide ATC code L02BB01), gonadotropin-releasing hormone (GnRH) agonists (L02AE). We have previously demonstrated a very high adherence to ADT and that it is rare that a man prescribed GnRH agonists does not continue his medication [23] [24] [25] [26] [27] . Men who in addition to GnRH agonist received 30 days of flare prophylaxis with bicalutamide were classified as exposed to GnRH agonist only. Data on bilateral orchidectomy (procedure code KFC10 or KFC15) were extracted from the National Patient Register.

The number of deaths per day from all causes was counted by use of the date of death. These counts were used to compare mortality in 1 March to 30 June in 2020 to corresponding months in 2015-2019. Mortality for the entire Swedish population during these time periods is also provided to contextualize the results.

Age at death (stratified), county of residence, prostate cancer risk category (absence/presence of distant metastases at the date of prostate cancer diagnosis), history of curative treatment (radical prostatectomy or radiotherapy), the Charlson Comorbidity Index and a Drug Comorbidity Index [28] were used as covariates in regression models. The Charlson Comorbidity Index was calculated from hospital discharge diagnoses registered in the National Patient Register during the 10-year period prior to start of follow-up [29, 30] . The Drug Comorbidity Index was calculated from filled prescriptions during a 365-day interval preceding the start of follow-up. Prescribed medications were categorized by chemical subgroup, i.e. the first five positions of the ATC code [28, 31] . Only ATC codes prevalent in at least 1% of observed deaths were used to calculate the Drug Comorbidity Index. Calendar year and day in the study period were also used as covariates to generate more comparable time periods.

The mortality rates by calendar day for each exposure group and calendar time period were calculated and plotted with locally weighted smoothing (LOESS). Using the log of numbers at risk as offset in a Poisson regression the rate of death was compared between 2020 and the average from previous years. Complementary analyses were performed in which the analyses were restricted to the 10 regions with the highest incidence of COVID-19, and the calendar period between 2 April to 10 June when excess mortality was above 30%.

Baseline characteristics and exposure for study men are described in Table 1 .

There was a higher number of deaths from all causes in the general Swedish population during the period 1 March to 30 June 2020 compared to the mean for corresponding periods in the previous five years (Fig 1) . In total, 9,822 men with prostate cancer died during March-June the years 2015-2020, of whom 5,034 men were on ADT and 4,788 men were not on ADT (Fig 2 and Table 2 ). The increase in mortality in March-June 2020 was slightly higher in men on GnRH vs men on bicalutamide. There were large differences in mortality rate according to the type of ADT, but the difference in relative mortality rate comparing the year 2020 to the previous five years within exposure groups was small.

Crude comparisons indicated a 7-14% lower excess mortality in men on any ADT compared to men not on ADT and adjustment for potential confounders attenuated these estimates. The crude relative mortality rate ratio comparing the increase in March-June 2020 vs the same period in 2015-2019 for men on any ADT compared to men not on ADT was 0.86 (95% confidence interval (CI) 0.78 to 0.95). This difference was attenuated after adjustment for covariates to 0.96 (95% CI 0.87 to 1.06) ( Table 2) .

When restricting the analysis to the regions with the highest incidence of COVID-19, the results were essentially identical as in the main analysis ( Table 3 ). In the analyses restricted to the time period between 2 April to 10 June when mortality in 2020 was increased >30% compared to previous years, men on GnRH had an adjusted relative mortality rate ratio of 0.89 (95% CI 0.77 to 1.03) and men on bicalutamide monotherapy 0.94 (95% CI 0.78 to 1.14).

In this population-based cohort of men with prostate cancer, we found no difference in relative excess mortality in men with prostate cancer on GnRH agonists or bicalutamide monotherapy vs those not on ADT during the peak of the first wave of the COVID-19 pandemic in March-June 2020, compared to the average from corresponding months in 2015-2019. There was only a weak reduction in one of the restricted subgroup analyses, seen in both exposure groups. Considering the mechanistic hypothesis an effect would mainly be expected from exposure to GnRH agonists but no notable difference was seen in the outcome compared to bicalutamide monotherapy. Thus, there was no clear support for the hypothesis that ADT mitigates the course of COVID-19.

The hypothesis that ADT mitigates the course of COVID-19 has been corroborated by in vitro studies showing that androgens facilitate cellular uptake of the SARS-CoV-2 virus. In humans, men have a 50% higher risk of infection and hospitalization, critical care admission, and death from COVID-19 compared with women [32] [33] [34] [35] . This increase in risk remained after adjustment for age, ethnicity, socio-economic factors, smoking, and concomitant and previous diseases [1] .

It is therefore of interest to investigate if ADT mitigates COVID-19. An early study suggested a lower risk of COVID-19 in men with prostate cancer on ADT compared with men not on ADT [10] . Similarly, in a small observational study of 22 men on ADT these men had lower rates of hospitalization, supplemental oxygen requirements, endotracheal intubation, and mortality than 36 control men not on ADT [36] . In contrast, other pilot studies have reported no protective effect with ADT [11, 12, 37, 38] . However, all of these studies were severely hampered by their small number of exposed men and restricted adjustment for putative confounding factors. In our main analysis of almost 10,000 men with prostate cancer who died during the study period, and of whom half were on ADT, there was no difference in excess mortality in men on any type of ADT vs men not on ADT. In a restricted analysis of ten regions with the highest COVID-19 burden the results were similar to that in the full study group. In a analysis restricted to the period with the highest mortality, we observed a 12% decrease in excess mortality in men on ADT.

Despite our efforts to adjust for putative confounders, there was likely residual confounding. One way to assess confounding is to use the E-value which is defined as the weakest association that an unmeasured confounder would need to have with both exposure and outcome to fully eliminate an association [39] . The E-value was 1.53 for a relative mortality rate ratio of 0.88, which was the most extreme point estimate in all our analyses, including restricted subgroup analyses. An imbalance of exposure to ADT related to prostate cancer risk category [40, 41] as well as an association of risk category with mortality [42] of that magnitude is possible. The finding of no association in the main analysis and a weak association in only one out of While there is a hypothetical mechanism for ADT to be potentially beneficial in the early phase of COVID-19, there are also suggestions and preliminary data that ADT could be detrimental in later stages of severe COVID-19, since it could exacerbate or activate an excessive and harmful inflammatory response [43] . It has even been hypothesized that testosterone replacement therapy could be beneficial [44] .

There are some other mechanisms, unrelated to suppression of androgen signaling, for an apparent impact of ADT on incidence and severity of COVID-19. For instance, men with advanced prostate cancer on ADT are, compared to men with less advanced cancer who are not on ADT, more likely to impose stricter protective measures such as physical distancing, self-isolation, strict hand hygiene, and wearing a face mask.

A limitation of our study is that the excess in all-cause mortality we observed may not be exclusively caused by COVID-19. The increased mortality during this pandemic has also to some part been attributed to reluctance to seek medical attention for other serious conditions [45] . In a recent study based on data from 29 high-income countries the difference between reported COVID-19 deaths and excess mortality in Sweden was among low, suggesting a high rate of correct adjudication of cause of death in Sweden [46] . Another limitation of our study is the lack of measurements to verify castrate testosterone levels. We found, however, high adherence to GnRH in a recent study in PCBaSe [23] and previous studies have shown that a large proportion of men on GnRH have castrate levels [47] . It would also have been of interest to study men with castration-resistant prostate cancer but we could not identify these men with the data available. Finally, it is unknown if the results from this observational study of old men with prostate cancer on ADT is applicable to the populations in current clinical trials employing short-term ADT in mostly middle-aged men diagnosed with COVID-19.

The strengths of this study are that it was based on a national, population-based cohort, with almost 10,000 deaths and that 5,000 of the deaths were among men on ADT, making it substantially larger than previous human studies combined. The study was performed in a national population of over 10 million people, with high mortality from COVID-19 in the spring of 2020. Furthermore, the study was based on a clinical cancer register and other health care registers with virtually complete capture of comprehensive data with known high quality [14, 48, 49] .

In conclusion, there was no difference in relative excess mortality in men with prostate cancer on ADT vs those not on ADT during the peak of the first wave of the COVID-19 pandemic in March-June 2020 compared to the average from corresponding time periods in previous years. Thus, we found no clear support for the hypothesis that androgen deprivation therapy mitigates the disease course of COVID-19. Analyses restricted to regions with highest incidence of COVID-19 by end of June 2020, and to the time period 2 April to 10 June 2020 when there was a peak in excess mortality.

https://doi.org/10.1371/journal.pone.0255966.t003

Disclaimer

Risk factors for SARS-CoV-2 among patients in the Oxford Royal College of General Practitioners Research and Surveillance Centre primary care network: a cross-sectional study

Male sex identified by global COVID-19 meta-analysis as a risk factor for death and ITU admission

A clinical risk score to identify patients with COVID-19 at high risk of critical care admission or death: An observational cohort study

Risk Factors Associated With Mortality Among Patients With COVID-19 in Intensive Care Units in Lombardy, Italy

Factors Associated With Death in Critically Ill Patients With Coronavirus Disease 2019 in the US

Epidemiology and Clinical Presentation of Children Hospitalized with SARS-CoV-2 Infection in Suburbs of Paris

COVID-19 in children and adolescents in Europe: a multinational, multicentre cohort study. Lancet Child Adolesc Health

SARS-CoV-2 Cell Entry Depends on ACE2 and TMPRSS2 and Is Blocked by a Clinically Proven Protease Inhibitor

Targeting Metastatic Hormone Sensitive Prostate Cancer: Chemohormonal Therapy and New Combinatorial Approaches

Androgen-deprivation therapies for prostate cancer and risk of infection by SARS-CoV-2: a population-based study (N = 4532)

On the relationship between androgen-deprivation therapy for prostate cancer and risk of infection by SARS-CoV-2

Androgen deprivation and SARS-CoV-2 in men with prostate cancer

The National Library of Medicine (NLM) at the National Institutes of Health (NIH)

Quantifying observational evidence for risk of fatal and nonfatal cardiovascular disease following androgen deprivation therapy for prostate cancer: a meta-analysis

Risk and timing of cardiovascular disease after androgen-deprivation therapy in men with prostate cancer

Preexisting Comorbidities Predicting COVID-19 and Mortality in the UK Biobank Community Cohort

Sex differences in SARS-CoV-2 infection rates and the potential link to prostate cancer

Cohort Profile Update: The National Prostate Cancer Register of Sweden and Prostate Cancer data Base-a refined prostate cancer trajectory

The completeness of the Swedish Cancer Register: a sample survey for year 1998

The First Decade with the Swedish Prescribed Drug Register-A Systematic Review of the Output in the Scientific Literature

External review and validation of the Swedish national inpatient register

PCBaSe Sweden: A register-based resource for prostate cancer research

Long-term adherence to GnRH agonists in men with prostate cancer. A nation-wide population-based study in prostate cancer data base Sweden

Anti-androgen prescribing patterns, patient treatment adherence and influencing factors; results from the nationwide PCBaSe Sweden

Persistence and adherence to androgen deprivation therapy in men with prostate cancer: an administrative database study

Adherence to guidelines for androgen deprivation therapy after radical prostatectomy: Swedish population-based study

Adherence Measures for Patients with Metastatic Castration-Resistant Prostate Cancer Treated with Abiraterone Acetate plus Prednisone: Results of a Prospective, Cluster-Randomized Trial. Cancers (Basel)

An aggregated comorbidity measure based on the history of filled drug prescriptions-development and application in two separate cohorts

Coding algorithms for defining comorbidities in ICD-9-CM and ICD-10 administrative data

A new method of classifying prognostic comorbidity in longitudinal studies: development and validation

Prescription-based prediction of baseline mortality risk among older men

Mortality rate and gender differences in COVID-19 patients dying in Italy: A comparison with other countries

Presenting Characteristics, Comorbidities, and Outcomes Among 5700 Patients Hospitalized With COVID-19 in the New York City Area

Gender Differences in Patients With COVID-19: Focus on Severity and Mortality. Front Public Health

Analysis on 54 Mortality Cases of Coronavirus Disease 2019 in the Republic of Korea from

Accounting for missing data in statistical analyses: multiple imputation is not always the answer

Incidence and outcomes of severe acute respiratory syndrome coronavirus 2 infection in patients with metastatic castration-resistant prostate cancer

Androgen Deprivation Therapy in Men With Prostate Cancer Does Not Affect Risk of Infection With SARS-CoV-2

Sensitivity Analysis in Observational Research: Introducing the E-Value

Guidelines on Prostate Cancer. Part II-2020 Update: Treatment of Relapsing and Metastatic Prostate Cancer

Regional cancer centers in collaboration

Long-term outcomes among noncuratively treated men according to prostate cancer risk category in a nationwide, population-based study

Androgen-deprivation therapy and SARS-Cov-2 infection: the potential double-face role of testosterone

Role of testosterone in COVID-19 patients-A double-edged sword? Med Hypotheses

Excess Deaths From COVID-19 and Other Causes

Excess deaths associated with covid-19 pandemic in 2020: age and sex disaggregated time series analysis in 29 high income countries

Impact of surgical and medical castration on serum testosterone level in prostate cancer patients

Capture rate and representativity of The National Prostate Cancer Register of Sweden

Evaluation of data quality in the National Prostate Cancer Register of Sweden