Skip to main content
  • Research article
  • Open access
  • Published:

Estrogen-modulating treatment among mid-life women and COVID-19 morbidity and mortality: a multiregister nationwide matched cohort study in Sweden

Abstract

Background

It has been repeatedly shown that men infected by SARS-CoV-2 face a twofold higher likelihood of dying, being hospitalized or admitted to the intensive care unit compared to women, despite taking into account relevant confounders. It has been hypothesized that these discrepancies are related to sex steroid hormone differences with estrogens being negatively correlated with disease severity. The objective of this study was therefore to evaluate COVID-19-related mortality and morbidity among peri- and postmenopausal women in relation to estrogen-containing menopause hormonal treatments (MHT).

Methods

This is a national register-based matched cohort study performed in Sweden between January 1 to December 31, 2020. Study participants comprised women over the age of 53 years residing in Sweden. Exposure was defined as prescriptions of local estrogens, systemic estrogens with and without progestogens, progestogens alone, or tibolone. MHT users were then compared with a matched cohort of non-users. The primary outcome consisted of COVID-19 mortality, whereas the secondary outcomes included inpatient hospitalizations/outpatient visits and confirmed SARS-CoV-2 infection. Multivariable adjusted Cox regression-derived hazard ratios (HRs) were calculated.

Results

Use of systemic estrogens alone is associated with increased COVID-19 mortality among older women (aHR 4.73, 1.22 to 18.32), but the association is no longer significant when discontinuation of estrogen use is accounted for. An increased risk for COVID-19 infection is further observed for women using combined systemic estrogens and progestogens (aHR 1.06, 1.00 to 1.13) or tibolone (aHR 1.21, 1.01 to 1.45). Use of local estrogens is associated with an increased risk for COVID-19-related death (aHR 2.02,1.45 to 2.81) as well as for all secondary outcomes.

Conclusions

Systemic or local use of estrogens does not decrease COVID-19 morbidity and mortality to premenopausal background levels. Excess risk for COVID-19 morbidity and mortality was noted among older women and those discontinuing systemic estrogens. Higher risk for death was also noted among women using local estrogens, for which non-causal mechanisms such as confounding by comorbidity or frailty seem to be the most plausible underlying explanations.

Trial registration details

Not applicable.

Peer Review reports

Background

Since the beginning of the coronavirus disease-19 (COVID-19) pandemic, several studies have aimed to comprehensively elucidate the biological pathways that underlie the vast differences in clinical course of the disease and which factors may influence these pathways. Globally, approximately 500 million people have been infected so far by the novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) with men being equally, if not slightly less, susceptible to infection than women [1,2,3,4]. However, men seem to experience almost two-fold higher mortality rates than women after controlling for potential confounding factors [5]. Similar patterns of disease severity have been observed with previous coronavirus outbreaks, including SARS-CoV (severe acute respiratory syndrome coronavirus) and MERS-CoV (Middle East respiratory syndrome coronavirus) [6,7,8].

The etiology of sex differences in COVID-19 severity and clinical course remains obscure. Several mechanisms have been hypothesized to underlie these differences including the role of sex hormones. One factor that has been suggested is the importance of estrogens in female immune responses, such as suppression of pro-inflammatory cytokines (i.e., IL-1β and IL-6), increased antibody production by B cells, and modulation of the activity of ACE2 [9, 10]. ACE2 is the receptor to which the virus’ spike protein binds when entering the cell. Although the exact mechanisms are yet to be understood, estrogens seem to play a role in preventing the cytokine storm release after a SARS-CoV-2 infection [9, 10]. In fact, accumulating data indicate that it is the cytokine storm and the subsequent respiratory distress and multiorgan failure that account for the high mortality observed among COVID-19 affected patients [11]. Biological sex differences in terms of sensitivity to severe COVID-19 may thus explain why fewer women than men have been hospitalized, have been admitted to intensive care units (ICU), and have died during the pandemic. However, only a handful of studies have yet examined the potential role of estrogen containing drugs on COVID-19 mortality [4, 12,13,14,15]. The existing studies indicate a protective effect of estrogen-containing medications against COVID-19 mortality among women already infected by SARS-CoV-2 [4, 12, 14]. However, all prior findings originate from data comprising combinations of exogenous menopausal hormone replacement therapy (MHT) with different drug preparations and administration routes, thus not allowing to delineate the potential effects of local versus systemic estrogen preparations, which often also refer to different dosages, as well as combinations with progestogens or not, against severe COVID-19. More evidence is needed to guide clinical recommendations.

The aim of this study is to evaluate COVID-19-related mortality and morbidity among peri- and postmenopausal women in relation to menopause hormonal treatments (MHT). This nationwide register-based study tests the hypothesis that MHT/estrogen-modulating treatment reduces COVID-19 mortality, disease severity, and risk for SARS-CoV-2 infection, among peri- and postmenopausal women in Sweden. Detailed register information allows for valid assessment of preparations containing estrogen, progestogen, and potential combinations, while also taking the hormone substance and administration route into account.

Methods

Study design

We performed a nationwide register-based matched cohort study in Sweden. Data were prospectively collected and retrieved pseudonymized after cross-linkages across different national socioeconomic and healthcare registries. Linkages were based on unique personal identification numbers of the study participants; all Swedish residents are assigned a unique 12-digit personal identification number at birth or upon immigration [16]. Data were assembled from the following registries: the National Patient Registry (NPR), including all in-patient care and outpatient specialists visits and related diagnoses from both private and public caregivers in Sweden [17]; the Prescribed Drug Registry (PDR), including information on the Anatomical Therapeutic Code (ATC-code), drug name, strength, date of prescription, and dispensing; the Cause of Death Registry (CDR), including information on cause of death, date, and place of death [18]; the Total Population Registry (TPR), including information on country of birth, sex, and marital status as well as immigration and emigration [19]; the Longitudinal Integration Database for Health Insurance and Labour Market Studies (LISA) [20] including information on educational attainment of the population; Sminet, a surveillance system collecting notifications for 60 communicable diseases classified either as dangerous for public health (such as HIV and hepatitis A–E) or dangerous for society (such as SARS-CoV-2 and Ebola) according to the Communicable Diseases Act [21]. Data from the National Patient Registry, the Prescribed Drug Registry, and the Cause of Death Registry were provided by the Swedish National Board of Health and Welfare, while the government agency Statistics Sweden provided data for the Total Population Registry and the Longitudinal Integration Database for Health Insurance and Labour Market Studies. The Public Health Agency of Sweden provided data for Sminet.

Patients and/or the public were not involved in this study.

Study population

The study population comprises women residing in Sweden during January 1, 2020–December 31, 2020, who were 53 years or older at the study start, in order to increase the chance that only postmenopausal women were included in the population [22]. The exposed cohort comprises all women who had filled a prescription of any type of estrogen-modulating treatment with supply covering January 1, 2020, and potentially extending beyond this date. Information on estrogen-modulating treatment was retrieved from the PDR using ATC codes (Additional file 1: Appendix 1). For each exposed woman, up to five women without supply of any of these prescriptions of the exposed individual were randomly selected from the underlying population and matched by age (± 2 years) and healthcare region. Women from the comparison group could be matched to more than one exposed woman belonging to different exposure groups. Individuals with gender dysphoria diagnosed at any time during the study period and identified by the ICD-10-SE (International Classification of Diseases, tenth revision, Swedish Edition) codes F640, F648, and F649 were excluded.

A flowchart displaying the included population is attached (Fig. 1).

Fig. 1
figure 1

A flowchart displaying the included population

Study variables

Exposure

MHT use was defined as filled prescriptions of any hormonal modulating treatment for which the estimated duration overlapped with the study start (January 1, 2020). The duration of filled prescriptions was determined based on the defined daily dose (DDD) (i.e., the assumed average maintenance dose for a drug used for its main indication in adults) of vaginal, oral, transdermal, subcutaneous, or intramuscular tablets/gel/cream/spray/suppositories/patches/injections dispensed, allowing for a gap of 30 days between dispensing dates. Systemic estrogens comprised both oral as well as transdermal preparations, while progestogens included mostly synthetic progestins, but also rarely bioidentical progesterone and dydrogesterone. The use of hormonal intrauterine device (IUD) was defined as a prescription up to 5 years before first filled prescription date (5-year look-back). For filled prescriptions with missing data on defined daily dose (DDD), the duration was imputed from recommended use. The ever exposure approach was used, i.e., each woman was allocated to the exposure group starting at January 1, 2020, with that exposure continuing until censoring, regardless of following dispensing or discontinuation. If two or more of the specified drugs were prescribed during the same time-period, the drug with higher estrogenicity was used for categorization. Non-use was defined as no modulatory hormonal treatment supply at January 1, 2020 (with the exception of hormonal IUD).

The ATC codes utilized to define exposure state comprised in short: G03C (estrogens), G03F (continuous/sequential combination estrogen-progestogen), G03D (progestogens), G03AC (progestogens), G02BA (progestogens), or G03CX01 (tibolone) (for details, see Additional file 1: Appendix 1).

Exposure groups.

  1. i)

    Women on local estrogens alone (group A)

  2. ii)

    Women on systemic estrogens without progestogens (group B)

  3. iii)

    Women on both estrogens and progestogens (group C)

  4. iv)

    Women on progestogens alone (group D)

  5. v)

    Women on tibolone alone (group E)

Follow-up and outcomes

The study participants were followed up for COVID-19 mortality or related morbidity until the first of either of the following events: death, emigration, study outcome, or the end of study (December 31, 2020).

The primary outcome was defined as death due to COVID-19 as either the main or underlying cause retrieved from the death certificates (according to the ICD-10-SE code U07.1 or U07.2). Secondary outcomes included the following: (i) inpatient hospitalizations (with or without the need for ICU admission) or outpatient visits due to COVID-19 (according to the ICD-10 codes U07.1 or U07.2) or (ii) laboratory COVID-19 confirmed events outside or inside a hospital setting (after having performed RT-PCR for SARS-CoV-2 from nasal or oropharyngeal swabs taken by healthcare personnel). Only outcomes registered after January 1, 2020, were considered valid. Data on the primary outcome of interest was retrieved from the CDR. Data on the secondary outcomes were retrieved from (i) the NPR and (ii) SmiNet, respectively. Outcomes were assessed independently; a person who was admitted in the hospital and then died from COVID-19 was included both in the death event counts (primary outcome) and in the inpatient hospitalizations/outpatient visits event counts (secondary outcome). Similar imputations were made within the secondary outcomes, i.e., all registered inpatient hospitalizations/outpatient visits events were included among the laboratory COVID-19 confirmed events.

Covariates

Data on the following sociodemographic factors were retrieved from the TPR and LISA registers: age at study inclusion (years), civil status (married/register partnership vs unmarried), income (low [below 20th percentile], middle [21st–79th percentile], or high [80th percentile and above] in relation to the income distribution in Sweden) and education level (< 9, 9–12, 13 years or more of total education attained). Regarding potential comorbidities and their effect on the outcome of interest a composite proxy, the modified Charlson Comorbidity Index was utilized. The Charlson Comorbidity Index (CCI) is a disease index often used as a proxy for comorbidity burden in statistical analysis. It is used to categorize comorbidities based on the ICD-10 codes of certain comorbid conditions recorded in the NPR (data included between January 1, 2016, through January 1, 2020). Each comorbidity category has an associated weight based on the adjusted risk of 1-year mortality, the sum of which results in a single comorbidity score for each patient. For the calculation of the modified CCI in our study, we included data based on COVID-19-related risk factors only. A complete list with the diagnoses and weights employed in the analysis is included in Additional file 1: Appendix 2. Proxies were used for alcohol abuse (ICD-10 code F10.2) and obesity (ICD-10 code E66) and were accounted for as independent risk factors for COVID-19 complications [23].

Statistical analysis

Main analysis

We first calculated unadjusted incidence rates (number of outcomes divided by accumulated person-years) of COVID-19-related mortality and morbidity among women with and without hormonal treatment. Rates of outcomes and covariates were calculated for different exposure groups. Cox regression analyses were performed in order to estimate hazard ratios with 95% confidence intervals (CIs) for the primary and secondary outcomes in relation to the exposure status. Regression models were adjusted for age at study inclusion, civil status, income, education level, obesity, alcohol abuse, and modified CCI. We have also repeated the regression analysis after adjusting for CCI categories as independent factors instead of as a composite proxy variable. Sensitivity analyses were performed in order to explore the impact of age as well as ongoing treatment on COVID-19 primary and secondary outcomes. Firstly, the analyses were repeated after stratifying the different exposure groups by age group at inclusion (i.e., 53–62 years, 63–72 years, and 73 years or over). Lastly, we studied the effect on all-cause mortality of “current ongoing treatment only” by censoring exposed women when they discontinued treatment or reached the end of supply before the end of the study (i.e., restricting the analyses among those with ongoing treatment as opposed to ever treatment). A two-sided p-value below 0.05 was considered statistically significant. All statistical analyses were performed using Statistical Analysis Software (SAS) version 9.4 (SAS Institute, Cary, NC, USA).

Further analyses

In order to evaluate the impact of different factors affecting the COVID-19 disease severity, we conducted the following explorative analyses:

  1. a)

    The impact on the primary and secondary outcomes of all preparations containing systemic estrogen with and without progestogen grouped together

  2. b)

    The competing risks effect, by assessing also all-cause mortality

  3. c)

    The effect of timing of COVID-19 morbidity and mortality in relation to the COVID-19 pandemic outbreak, by stratifying the primary and secondary outcomes in relation to the first (1 January 2020–31 August 2020) or second wave of COVID-19 pandemic (1 September 2020–31 December 2020)

  4. d)

    The impact of the route of administration of systemic estrogens, i.e., oral vs transdermal estrogen-including MHT, on the primary outcome

Results

Demographic characteristics

The exposed groups included 9981 women using local estrogens alone (group A), 3189 on systemic estrogens without progestogens (group B), 8352 on estrogens and progestogens (group C), 9323 on progestogens alone (group D), and 923 women receiving tibolone alone (group E) (Table 1). Mean age of the study population was 71.4 years (± standard deviation [SD]: 11.2) for group A, 58.2 (± 6.5) for group B, 57.2 (± 4.8) for group C, 55.8 (± 4.9) for group D, and 59.8 (± 6.0) for group E and was similar between the different exposure groups and corresponding comparison groups. A higher proportion (around 50%) of women exposed to systemic estrogen and/or progestogen (groups B-D) had attained more than 13 years of education compared with women in the matched cohort (42%), while no differences were observed between the comparison cohort and women using local estrogens alone or tibolone alone. Fewer women on systemic estrogen treatment or tibolone alone (groups B–C, E) had received a diagnosis of obesity compared with the matched cohort. Lastly, the majority (around 88%) of both exposed and unexposed women were healthy at baseline as indicated by a zero Charlson Comorbidity Index (CCI) with no numerical differences noted between the groups (Table 1).

Table 1 Demographic characteristics of the study population, exposed women, and matched unexposed women presented by drug preparation and route of administration

COVID-19 mortality

During follow-up, 114 cases of COVID-19 deaths were observed among the unexposed women and 50 among women who had used local estrogens (incidence rate, 5.3 per 1000 person-years) (Table 2), yielding a statistically significant hazard ratio of 2.02 (95% CIs: 1.45–2.81) for COVID-19 mortality in the adjusted models (Table 3). Women treated with systemic estrogens without progestogens also showed a statistically significant increased risk of COVID-19-related death compared with matched unexposed individuals (adjusted HR: 6.39, 95% CIs: 1.69–24.21). No deaths due to COVID-19 were reported for women treated with a combination of estrogens and progestogens or with tibolone alone with similar death rates estimated in the matched unexposed population. No clear difference in mortality risk was observed for women using progestogens alone compared with non-users (Table 3).

Table 2 Incidence rates (IR) for the association of sex steroid treatment with all primary and secondary study outcomes
Table 3 Crude and adjusted hazard ratios (HRs) for the association of sex steroid treatment with death due to COVID-19

COVID-19 morbidity and SARS-CoV-2 infection

Women using local estrogens alone experienced higher risk of outpatient visits/inpatient hospitalizations with or without the need for ICU admission (aHR: 1.23, 95% CIs: 1.10–1.38) and higher risk of SARS-CoV-2 infection (aHR: 1.13, 95% CIs: 1.06–1.21; Table 4). In the remaining exposed groups, no associations were found for the secondary outcomes of interest. The only exception was a significantly higher risk of COVID-19 infection for women exposed both to estrogens and progestogens (aHR: 1.06, 95% CIs: 1.00–1.13) and for women exposed to tibolone alone (aHR: 1.23, 95% CIs: 1.01–1.45; Table 4) compared with their matched non-exposed women.

Table 4 Crude and adjusted Cox regression models regarding the association of sex steroid treatment with laboratory confirmed SARS-CoV-2 infection or outpatient visits/inpatient hospitalizations with or without the need for ICU admission due to COVID-19

Sensitivity analyses

Analyses of different age groups showed that the risk of death due to COVID-19 was significantly increased for women above 63 years of age receiving local estrogens alone, with women aged 63–72 years exhibiting the highest risk (aHR 4.17, 95% CIs: 1.70–13.65). Only older women (i.e., 73 years or above) using systemic estrogens without progestogens exhibited increased risk of dying due to COVID-19 (aHR 4.73, 95% CIs: 1.22–18.32) (Tables 5 and 6). When only considering ongoing treatment, the positive association with COVID-19 mortality was no longer significant for systemic estrogens (Tables 7 and 8) and remained significant only for local estrogens alone (aHR 1.71, 95% CIs: 1.03–2.82).

Table 5 Crude and adjusted Cox regression models stratified by age regarding the association of sex steroid treatment with death due to COVID-19
Table 6 Crude and adjusted Cox regression models stratified by age regarding the association of sex steroid treatment with laboratory confirmed SARS-CoV-2 infection or outpatient visits/ inpatient hospitalizations with or without the need for ICU admission due to COVID-19
Table 7 Crude and adjusted Cox regression models evaluating the effect of treatment cessation in relation to death due to COVID-19 (i.e., restricting the analyses among those with ongoing treatment as opposed to ever treatment)
Table 8 Crude and adjusted Cox regression models evaluating the effect of treatment cessation in relation to laboratory confirmed SARS-CoV-2 infection or outpatient visits/inpatient hospitalizations with or without the need for ICU admission (i.e., restricting the analyses among those with ongoing treatment as opposed to ever treatment)

Further analyses

  1. a)

    Regardless of concomitant progestogen use or route of administration, we found that systemic estrogens conferred an increased risk of death by COVID-19 (aHR 4.22, 95% CIs: 1.37–13.05), but there was no apparent effect on the secondary outcomes (p > 0.05) (Additional file 2: Table S1).

  2. b)

    When assessing all-cause mortality, no differences were observed between exposure groups and matched comparison groups (Additional file 2: Table S2).

  3. c)

    With regard to the timing of the outbreak of the pandemic, we observed that the associations with COVID-19 morbidity and mortality were statistically significant mostly during the first wave and not consistently significant during the second wave (Additional file 2: Table S3).

  4. d)

    With regard to the route of administration of systemic estrogens without progestogens, we observed an increased risk of death from COVID-19 only among users of oral estrogens with an incidence rate of 4.6 (1.05–10.7), while no estimation could be performed among the users of transdermal preparations as there were no deaths observed in that group during the study period.

Throughout the analyses, in models using the distinct CCI contributors instead of the composite index, results remained substantially unchanged (data not shown).

Discussion

Overall, the present large Swedish cohort of peri- and postmenopausal women does not demonstrate any reduced risk for COVID-19 infection or related mortality among women using MHT compared with unexposed women, suggesting that MHT was not able to reverse theoretical increased risks among women with vasomotor symptoms or genitourinary symptoms. Initiating or continuing MHT treatment solely as a prophylactic treatment against COVID-19 should therefore be avoided. On the contrary, women on local estrogen alone experienced higher COVID-19 mortality, higher risk of outpatient visits/inpatient hospitalizations, and higher risk of SARS-CoV-2 infection. Also, there was an increased risk for death in the group on systemic estrogens without progestogens, but the excess risk was observed primarily in elderly women, and those discontinuing estrogen therapy during the follow-up period. When assessing mortality rates among women on ongoing estrogen therapy, no increased risk could be seen. Indication bias and residual confounding by co-morbidities not available to control for could thus at least partially explain these associations. In addition, the associations between estrogen use and COVID-19 outcomes were only observed during the first wave of the pandemic, further suggesting that non-causal mechanisms may explain the observed excess risks.

Interpretation

Increased risk of death from COVID-19 in women treated with systemic estrogens with or without progestogens as well as local estrogens have not been reported before. These findings, at first look, seem to disagree with what would be expected based on pathophysiological mechanisms and findings from the limited existing literature and were admittedly unanticipated. However, in a sensitivity analysis, we observed that the association between systemic estrogen therapy and COVID-19 mortality was no longer significant when exposure was defined as ongoing treatment with systemic estrogen therapy and not ever-treatment. The unfavorable effect of systemic estrogens on COVID-19 mortality may thus be explained by estrogen withdrawal rather than estrogen continuation, an observation that is in line with what other studies have so far demonstrated. Indeed, this fraction of peri- and postmenopausal women receiving MHT treatment probably represent a population of lower baseline estrogen levels who however restore their levels to those found in pre-menopause during treatment, i.e., a period in life with lower risk for severe COVID-19. However, when these women interrupted their treatment, their hormone levels fell to postmenopausal levels, and it is only then that the undesired severe effects of COVID-19 were noted. The latter observation strengthens the notion for a protective role of estrogens against COVID-19 and is in agreement with the literature, possibly implying that the relationship between COVID-19 and sex hormones is more complex than initially hypothesized. Although the exact mechanism is still unclear, it is believed that estrogen receptor a (ERα) directly interacts with the spike protein of the SARS-CoV-2 in certain tissues leading to a modified signaling pathway, which potentially affects SARS-CoV-2 infection and related pathology [24].

Estrogens were highlighted already from the beginning of the pandemic as playing a central role in COVID-19 morbidity and mortality due to the increased risk of deep vein thrombosis in COVID-19 infected individuals [25]. Thus, some exposed women at higher risk might have taken the initiative to discontinue their medications, especially systemic estrogen-containing preparations, during the study period, inducing a risk of misclassification. Likewise, women with systemic estrogen containing treatments might have been urged by their prescribing physician to seek medical attention earlier or to a greater extent if having COVID-19 like symptoms inducing a risk of over-ascertainment of some secondary outcomes (i.e., laboratory confirmed SARS-CoV-2 infection), explaining the higher risk for infection noted in some groups (surveillance bias). It is therefore plausible and important to consider that some of the associations noted between estrogen-modulating treatment and COVID-19 disease could also partly be explained by different types of bias which could not be controlled for in this study design, rather than suggesting causality and not by the treatment itself. It is of note that treatment with systemic estrogens and progestogens as well as tibolone exhibited neither a protective nor a harmful effect.

Furthermore, a higher risk for all morbidity outcomes as well as COVID-19 mortality were noted for women treated solely with local estrogens. It is generally known that the estrogen plasma levels that are reached due to local estrogen use are usually very low and do not differ greatly from the levels reported in healthy untreated postmenopausal women (i.e., 4 pg/ml increase in estradiol levels during use) [26]; one would therefore not expect a pharmacologic effect systemically. Local estrogens are usually prescribed against genitourinary symptoms of menopause. In fact, menopausal symptoms, both vasomotor and genitourinary, have been reported to reflect an underlying estrogen deficiency and endothelial dysfunction [27, 28], making these women more susceptible for cardiovascular events and complications [29, 30]. Moreover, women with menopausal symptoms, especially vaginal atrophy, are more prone to local infections (e.g., urinary tract infection (UTI), recurrent candidiasis etc.) [31, 32] with newer data indicating a link even to systemic infections, among which COVID-19 [15, 33, 34]; we could therefore hypothesize that women prescribed local estrogens were at baseline at higher risk for infections in general, and thus even for COVID-19 (potential indication bias).

It is of note that the majority of COVID-19 fatal events in Sweden among the elderly were observed in care homes and were primarily seen during the first wave of the pandemic, when doctors were not adequately trained to treat COVID-19 effectively [35]. During the same period, the use of protective equipment (facemasks/gloves) for health care workers was not mandatory or even encouraged, increasing their risk of transmitting the disease [36]. Thus, in general, women actually being prescribed local estrogens present with factors rendering them at higher risk for both infection and morbidity due to COVID-19, which could explain our findings.

Comparison with related studies

Our study findings are in line with that of Costeira et al.; in a population-based cohort from the UK including menopausal women on MHT or other hormonal therapies (e.g., combined oral contraceptive pills, COCPs), the authors demonstrated increased risks of predicted COVID-19 for MHT users alone [13]. The authors did not report on mortality, but no increased risk of hospitalization was seen among MHT users alone which was also confirmed in our study. Unfortunately, Costeira et al. [13] did not collect information on MHT type or route of administration. On the contrary, a handful of studies have suggested a protective effect of MHT use against COVID-19 severity and even the risk of contracting the disease. A recently published study from Sweden included three groups of postmenopausal women, namely women with breast cancer receiving tamoxifen or aromatase inhibitors, women receiving MHT, and control women not receiving estrogen-modulating treatment [12]. In this study, Sund et al. demonstrated a lower adjusted risk of death due to COVID-19 in women under MHT compared to the control group with a reversed relationship between estrogen levels and risk of COVID-19-related death [12]. However, it should be noted that the groups studied by Sund et al. did not arise from the general background population of perimenopausal women but originated instead from the population with laboratory confirmed SARS-CoV-2 infection. In fact, during the first wave of the pandemic in Sweden, immigrants from low- and middle-income countries and/or increased household size contracted the virus to a greater extent and were affected disproportionately by COVID-19 [37]. In line with that, the control group in the study by Sund et al. [12] included a higher proportion of individuals with worse socioeconomic status (i.e., poorest income quintile) and low education level (i.e., primary education), both risk factors that contribute to the increased mortality due to COVID-19 [38,39,40]. Thus, despite the fact that the authors adjusted for socioeconomical background factors, we believe that the intrinsic risk of the “control population” was already amplified, tipping the risk balance in favor of the study group. On the contrary, the comparison group in our study was matched according to age and healthcare region, and the regression models were also adjusted for education and socioeconomic status, limiting the risk of confounding. Likewise, to Sund et al., a retrospective cohort study from the UK by Dambha-Miller et al. showed that hormone replacement treatment was associated with lower all-cause mortality in COVID-19 positive women (aOR: 0.22, 95% CIs: 0.05–0.94), whereas no associations were found for women receiving COCPs. The latter finding indicated that the difference in COVID-19 mortality might be the result of a greater increase in hormone levels observed in perimenopausal women on MHT compared to premenopausal women on COCPs [14]. Again, the control group chosen was older, socioeconomically deprived, and with a higher rate of comorbidities, while the outcome regarded all-cause mortality and not just COVID-19-specific mortality. In addition, an analysis of electronic health records for a large (n = 68,466), international COVID-19 cohort [4] by Seeland et al., demonstrated that women older than 50 years under estradiol treatment had lower fatality risk, by more than 50%, compared to non-treated women (HR: 0.29, 95% CIs: 0.11–0.76) [4]. It should however be noted that the study population originated once more from the infected population probably affecting the intrinsic COVID-19 fatality risk of the control group [4]. Lastly, a multinational, cross-sectional retrospective study from Latin America including mid-aged women attending a routine health check-up reported that women on combined MHT (e.g., containing both estrogen and progestogen) presented a lower risk of laboratory confirmed SARS-CoV-2 infection (OR 0.62, 95% CI 0.41–0.94) [15], which comes in contrast to our findings. However, no association was found between estrogen-only containing MHT and COVID-19 infection, similarly to our study [15]. It should however be noted that participants in the study were relatively affluent (attending private clinical centers), younger (aged 40–64 years), and therefore not representative of the general population. Furthermore, due to the study design, women who experienced more severe COVID-19 disease or persistent symptoms were not included in the population and mortality could not be assessed.

Strengths and limitations

The main study strength is the fact that we used data from a large nationwide cohort in a country with high COVID-19 incidence during the first and second COVID-19 wave, lowering the risk of selection bias and making results generalizable to the whole population of peri-post menopausal women in Sweden. Moreover, data was obtained from multiple well-validated registries, thus limiting the risk for report bias. Lastly, we had access to a wide range of covariates with a potential role in COVID-19 morbidity and mortality and were therefore able to adjust for several confounders.

The study limitations include lacking information regarding compliance to pharmaceutical treatment such as MHT as well as purchases over the counter during the pandemic (e.g., local estrogens); as a consequence, misclassification cannot be entirely ruled out. Lockdowns were not implemented in Sweden and potential drug shortages during that period of the pandemic would not affect our estimates since these are based on actually dispensed medications before the pandemic outbreak.

Furthermore, we lack data on the menopausal state of the participants since menopause is not always recorded in the register through an ICD-10 diagnosis code. Thus, a proportion of the non-exposed population are expected to be premenopausal and not have the same risk for COVID-19 outcomes as those in peri- and post-menopause exposed to MHT [41]; this proportion, however, of the comparison group is expected to be lower with advancing age. The sources of exposure misclassification above should in any case not have induced false associations but would rather attenuate existing ones. In addition, despite adjusting for a variety of comorbid conditions, residual unmeasured confounding is still possible. Indication bias, especially in the case of those needing treatment with local estrogens vaginally who might have a higher infection susceptibility, cannot be ruled out. We also lacked data on the use of immunosuppressant drugs (such as corticosteroids) that could potentially alter the COVID-19 infection severity. Furthermore, free testing at healthcare facilities was established in the middle of June 2020. Due to somewhat lower testing capacity at the very beginning of the pandemic, there are certainly individuals positive for COVID-19 who were never tested and recorded in the system. This might have introduced some degree of misclassification of this outcome, which nevertheless is expected to be non-differential; this could have led to possible attenuation of some associations. In an effort to account for this, relevant imputations regarding secondary outcome definitions were performed. In addition, we lack information on the BMI of the participating women; we have instead adjusted the regression analyses for the ICD-10 diagnosis of obesity, which has a similar detrimental effect on COVID-19 morbidity and mortality. It should however be noted that registry data only capture the most severe cases, increasing the risk of missing individuals with obesity and incorrectly classifying them as not having the condition. Nevertheless, based on contraindications for MHT use, the misclassification would mostly affect the unexposed population and could have eventually attenuated our findings; it is not expected to have introduced false associations. In addition, initiation of MHT after the study start did not constitute an exclusion criterium per se; however, because of the exposure definition used in the study, women with later dated prescriptions could nevertheless not be included in our analyses, as they were censored after prescription date. Moreover, despite the population-based character of the study, the relative rarity of certain exposures (such as use of bioidentical progestogens) and outcomes restricts our capacity to detect some exposure-outcome associations and thus lowers the study’s statistical power in relation to certain research questions. Lastly, asymptomatic individuals or those with milder symptoms that never sought health care were not captured in this study, affecting the (less severe) secondary outcome events of COVID-19; however, that should not affect the events of the primary outcome (death) nor the severe secondary outcomes (inpatient hospitalization with or without ICU admission).

Conclusions

In this population-based register study in Sweden, we could not confirm that MHT succeeded in reversing the theoretically increased risks for COVID-19 complications among menopausal women, with indications for use. Initiating or continuing MHT treatment solely for this purpose is therefore advised against. On the contrary, an increased risk for COVID-19 mortality was observed among women on local estrogens compared to non-exposed women. This finding can be due to indication bias and unmeasured confounding, such as frailty and susceptibility for infections among women on local estrogen therapy. Increased risks among older women on systemic estrogens at the start of the study period were no longer apparent when addressing current use. The specific role of MHT cessation for COVID-19-related mortality in this group warrants further investigation. No increased risks for COVID-19 mortality were observed among women on MHT prescribed in ages according to existing clinical recommendations. The findings need to be interpreted with caution, as they were mostly evident during the first pandemic wave and are subject to inherent limitations in register data; nevertheless, they do need to be explored further, in order to guide clinical recommendations.

Availability of data and materials

Restrictions apply to the availability of some or all data generated or analyzed during this study to preserve patient confidentiality or because they were used under license. The corresponding author will on request detail the restrictions and any conditions under which access to some data may be provided.

Abbreviations

CCI:

Charlson Comorbidity Index

CI:

Confidence interval

COVID-19:

Coronavirus disease-19

DDD:

Defined daily dose

HR:

Hazzard ratio

ICD:

International Classification of Diseases

IUD:

Intrauterine device

MHT:

Menopause hormonal therapy

SARS-CoV-2:

Severe acute respiratory syndrome coronavirus 2

References

  1. Cagnacci A, Xholli A. Age-related difference in the rate of coronavirus disease 2019 mortality in women versus men. Am J Obstet Gynecol. 2020;223(3):453–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Williamson EJ, Walker AJ, Bhaskaran K, Bacon S, Bates C, Morton CE, et al. Factors associated with COVID-19-related death using OpenSAFELY. Nature. 2020;584(7821):430–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Meister T, Pisarev H, Kolde R, Kalda R, Suija K, Milani L, et al. Clinical characteristics and risk factors for COVID-19 infection and disease severity: a nationwide observational study in Estonia. PLoS ONE. 2022;17(6):e0270192.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Seeland U, Coluzzi F, Simmaco M, Mura C, Bourne PE, Heiland M, et al. Evidence for treatment with estradiol for women with SARS-CoV-2 infection. BMC Med. 2020;18(1):369.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Scully EP, Haverfield J, Ursin RL, Tannenbaum C, Klein SL. Considering how biological sex impacts immune responses and COVID-19 outcomes. Nat Rev Immunol. 2020;20(7):442–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Karlberg J, Chong DS, Lai WY. Do men have a higher case fatality rate of severe acute respiratory syndrome than women do? Am J Epidemiol. 2004;159(3):229–31.

    Article  CAS  PubMed  Google Scholar 

  7. Chan-Yeung M, Xu RH. SARS: epidemiology. Respirology. 2003;8 Suppl(Suppl 1):9–14.

    Google Scholar 

  8. Alghamdi IG, Hussain II, Almalki SS, Alghamdi MS, Alghamdi MM, El-Sheemy MA. The pattern of Middle East respiratory syndrome coronavirus in Saudi Arabia: a descriptive epidemiological analysis of data from the Saudi Ministry of Health. Int J Gen Med. 2014;7:417–23.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Averyanova M, Vishnyakova P, Yureneva S, Yakushevskaya O, Fatkhudinov T, Elchaninov A, et al. Sex hormones and immune system: menopausal hormone therapy in the context of COVID-19 pandemic. Front Immunol. 2022;13:928171.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Choi SW, Kim J, Lee JH, Kim SK, Lee SR, Kim SH, et al. Hormone therapy in the era of the COVID-19 pandemic: a review. J Menopausal Med. 2022;28(1):1–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Hojyo S, Uchida M, Tanaka K, Hasebe R, Tanaka Y, Murakami M, et al. How COVID-19 induces cytokine storm with high mortality. Inflamm Regen. 2020;40:37.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Sund M, Fonseca-Rodriguez O, Josefsson A, Welen K, Fors Connolly AM. Association between pharmaceutical modulation of oestrogen in postmenopausal women in Sweden and death due to COVID-19: a cohort study. BMJ Open. 2022;12(2):e053032.

    Article  PubMed  Google Scholar 

  13. Costeira R, Lee KA, Murray B, Christiansen C, Castillo-Fernandez J, Ni Lochlainn M, et al. Estrogen and COVID-19 symptoms: associations in women from the COVID Symptom Study. PLoS ONE. 2021;16(9):e0257051.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Dambha-Miller H, Hinton W, Wilcox CR, Joy M, Feher M, de Lusignan S. Mortality in COVID-19 among women on hormone replacement therapy: a retrospective cohort study. Fam Pract. 2022;39(6):1049–55.

    Article  PubMed  Google Scholar 

  15. Vallejo MS, Blumel JE, Bencosme A, Calle A, Dextre M, Diaz K, et al. Factors affecting climacteric women with SARS-CoV-2 infection: a multinational Latin America study (REDLINC XI). Maturitas. 2022;165:33–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Ludvigsson JF, Myrelid P. Swibreg–a new version of national IBD registry. Lakartidningen. 2009;106(45):3014–5.

    PubMed  Google Scholar 

  17. Ludvigsson JF, Andersson E, Ekbom A, Feychting M, Kim JL, Reuterwall C, et al. External review and validation of the Swedish national inpatient register. BMC Public Health. 2011;11:450.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Brooke HL, Talback M, Hornblad J, Johansson LA, Ludvigsson JF, Druid H, et al. The Swedish cause of death register. Eur J Epidemiol. 2017;32(9):765–73.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Ludvigsson JF, Almqvist C, Bonamy AK, Ljung R, Michaelsson K, Neovius M, et al. Registers of the Swedish total population and their use in medical research. Eur J Epidemiol. 2016;31(2):125–36.

    Article  PubMed  Google Scholar 

  20. Ludvigsson JF, Svedberg P, Olen O, Bruze G, Neovius M. The longitudinal integrated database for health insurance and labour market studies (LISA) and its use in medical research. Eur J Epidemiol. 2019;34(4):423–37.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Agency TPH. SMINET [Available from: https://www.folkhalsomyndigheten.se/sminet/om-sminet/.

  22. Sundell M, Brynhildsen J, Fredrikson M, Hoffmann M, Spetz Holm AC. Insufficient use of menopausal hormone therapy in Swedish women with early or premature menopause caused by bilateral oophorectomy: a register-based study. BJOG. 2023.

  23. Al-Lami RA, Urban RJ, Volpi E, Algburi AMA, Baillargeon J. Sex hormones and novel corona virus infectious disease (COVID-19). Mayo Clin Proc. 2020;95(8):1710–4.

    Article  CAS  PubMed  Google Scholar 

  24. Solis O, Beccari AR, Iaconis D, Talarico C, Ruiz-Bedoya CA, Nwachukwu JC, et al. The SARS-CoV-2 spike protein binds and modulates estrogen receptors. Sci Adv. 2022;8(48):eadd4150.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Ramirez I, De la Viuda E, Baquedano L, Coronado P, Llaneza P, Mendoza N, et al. Managing thromboembolic risk with menopausal hormone therapy and hormonal contraception in the COVID-19 pandemic: recommendations from the Spanish Menopause Society, Sociedad Espanola de Ginecologia y Obstetricia and Sociedad Espanola de Trombosis y Hemostasia. Maturitas. 2020;137:57–62.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Santen RJ, Mirkin S, Bernick B, Constantine GD. Systemic estradiol levels with low-dose vaginal estrogens. Menopause. 2020;27(3):361–70.

    Article  PubMed  Google Scholar 

  27. Angelou K, Grigoriadis T, Diakosavvas M, Zacharakis D, Athanasiou S. The genitourinary syndrome of menopause: an overview of the recent data. Cureus. 2020;12(4):e7586.

    PubMed  PubMed Central  Google Scholar 

  28. Waetjen LE, Crawford SL, Chang PY, Reed BD, Hess R, Avis NE, et al. Factors associated with developing vaginal dryness symptoms in women transitioning through menopause: a longitudinal study. Menopause. 2018;25(10):1094–104.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Thurston RC, Aslanidou Vlachos HE, Derby CA, Jackson EA, Brooks MM, Matthews KA, et al. Menopausal vasomotor symptoms and risk of incident cardiovascular disease events in SWAN. J Am Heart Assoc. 2021;10(3):e017416.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Cagnacci A, Gambera A, Bonaccorsi G, Xholli A, study A. Relation between blood pressure and genito-urinary symptoms in the years across the menopausal age. Climacteric. 2022;25(4):395–400.

    Article  CAS  PubMed  Google Scholar 

  31. Peters KJ. What is genitourinary syndrome of menopause and why should we care? Perm J. 2021;25.

  32. Caretto M, Giannini A, Russo E, Simoncini T. Preventing urinary tract infections after menopause without antibiotics. Maturitas. 2017;99:43–6.

    Article  PubMed  Google Scholar 

  33. Ciesielska A, Kusiak A, Ossowska A, Grzybowska ME. Changes in the oral cavity in menopausal women-a narrative review. Int J Environ Res Public Health. 2021;19(1).

  34. Gameiro CM, Romao F, Castelo-Branco C. Menopause and aging: changes in the immune system–a review. Maturitas. 2010;67(4):316–20.

    Article  CAS  PubMed  Google Scholar 

  35. Pressrum COVID-19. [cited 2022–11–08]. Available from: https://www.socialstyrelsen.se/om-socialstyrelsen/pressrum/press/ny-statistik-om-smittade-och-avlidna-i-covid-19-70-ar-och-aldre/.

  36. Ludvigsson JF. The first eight months of Sweden’s COVID-19 strategy and the key actions and actors that were involved. Acta Paediatr. 2020;109(12):2459–71.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Drefahl S, Wallace M, Mussino E, Aradhya S, Kolk M, Branden M, et al. A population-based cohort study of socio-demographic risk factors for COVID-19 deaths in Sweden. Nat Commun. 2020;11(1):5097.

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  38. Yoshikawa M, Asaba K. Educational attainment decreases the risk of COVID-19 severity in the European population: a two-sample Mendelian randomization study. Front Public Health. 2021;9:673451.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Niedzwiedz CL, O’Donnell CA, Jani BD, Demou E, Ho FK, Celis-Morales C, et al. Ethnic and socioeconomic differences in SARS-CoV-2 infection: prospective cohort study using UK Biobank. BMC Med. 2020;18(1):160.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Tang IW, Vieira VM, Shearer E. Effect of socioeconomic factors during the early COVID-19 pandemic: a spatial analysis. BMC Public Health. 2022;22(1):1212.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Wang XW, Hu H, Xu ZY, Zhang GK, Yu QH, Yang HL, et al. Association of menopausal status with COVID-19 outcomes: a propensity score matching analysis. Biol Sex Differ. 2021;12(1):16.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We are grateful for the support in retrieving registry data provided by Alva Enoksson-Wallas and Jonas Färnstrand.

Funding

Open access funding provided by Uppsala University. The authors have not received a separate grant from any funding agency for the submitted work. However, EE has a part-time research position funded by Uppsala University Hospital (grant number ALF 937815). The funder had no role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication.

Author information

Authors and Affiliations

Authors

Contributions

EE, HKK, KF, FP, and AS conceptualized the study. EE and MK drafted the manuscript. ML and DH carried out the statistical analyses. All authors interpreted the results and critically reviewed the manuscript. AS supervised the study. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Evangelia Elenis.

Ethics declarations

Ethics approval and consent to participate

The study was approved by the Swedish Ethical Review Authority (Dnr 2020/03936 with date of approval 17/08/2020). The need for written or oral informed consent was waived since all data received from the Swedish registries were pseudonymized.

Consent for publication

Not applicable.

Competing interests

All authors have completed the ICMJE uniform disclosure form at http://www.icmje.org/disclosure-of-interest/(available on request from the corresponding author) and declare the following: no support from any organization for the submitted work; EE has received lecture fee from Gedeon Richter and research grant from Ferring Pharmaceuticals and serves as the medical advisor of Tilly AB; AS has received a one-time consultant fee from Biogen Inc., none of which are in any way related to this manuscript. HKK receives honorariums for lectures on contraception, fibroid care, and endometriosis from Abbvie, Actavis, Bayer, Gedeon Richter Exeltis, Nordic Pharma, Natural Cycles, Mithra, Teva, Merck, Organon, Ferring, and Consilient Health and provides expert opinions for Bayer, Evolan, Gedeon Richter, Exeltis, Merck, Teva, TV4 och Natural Cycles, Pharmiva, Dynamic Code (past), Ellen, Estercare, Pharmiva, Gedea, Gesynta, Leia, and Essity. HKK is an investigator in trials sponsored by Bayer, MSD, Mithra, Ethicon, Azanta, Gedeon Richter, Gedea, and Pharmiva. HKK teaches in courses sponsored by Organon, Bayer, and Gedeon Richter. ML and DH report participation in research projects funded by pharmaceutical companies, all regulator-mandated phase IV-studies, all with funds paid to the institution where they were employed (no personal fees) and with no relation to the work reported in this paper. The other authors declare no conflict of interest.

Additional information

Publisher’s Note

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

Supplementary Information

Additional file 1: Appendix 1-2. Appendix 1.

Categorization of exposure groups according to their Anatomic Therapeutic Codes (ATC-codes). Appendix 2. Comorbidities included in the modified Charlson Comorbidity Index (CCI).

Additional file 2: Tables S1-S3. Table S1.

Crude and adjusted Cox regression models regarding the association of systemic estrogen treatment with/without progestogens in relation to death due to COVID-19, laboratory-confirmed SARS-CoV-2 infection or outpatient visits/inpatient hospitalizations with/without ICU admission. Table S2. Crude and adjusted Cox regression models regarding the association of sex steroid treatment in relation to all-cause mortality. Table S3. Crude and adjusted Cox regression models stratified by timing in the COVID-19 pandemic outbreak in relation to the first (1st January 2020–31st August 2020) or second wave of COVID-19 pandemic (1st September 2020–31st December 2020).

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Elenis, E., Kallner, H.K., Karalexi, M.A. et al. Estrogen-modulating treatment among mid-life women and COVID-19 morbidity and mortality: a multiregister nationwide matched cohort study in Sweden. BMC Med 22, 84 (2024). https://doi.org/10.1186/s12916-024-03297-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12916-024-03297-z

Keywords