Advertisement
Original Article| Volume 158, P61-69, April 2022

Sex- and gender-related differences linked to SARS-CoV-2 infection among the participants in the web-based EPICOVID19 survey: the hormonal hypothesis

      Highlights

      • COVID-19 is characterized by marked sex and gender differences.
      • The role of sex hormones in SARS-CoV-2 infection is largely unknown.
      • Female sex and younger age were negatively associated with the severity of the disease in this study.
      • Older females seem to benefit from hormonal therapy against the SARS-CoV-2 infection.

      Abstract

      Objective To investigate sex- and gender-based differences linked to SARS-COV-2 infection and to explore the role of hormonal therapy (HT) in females.
      Study design Data from the self-administered, cross-sectional, web-based EPICOVID19 survey of 198,822 adults living in Italy who completed an online questionnaire during the first wave of the epidemic in Italy (April-May 2020) were analyzed.
      Main outcomes measures Multivariate binary logistic and multinomial regression models were respectively used to estimate the odds ratios (ORs) with 95% confidence intervals (CIs) for positive nasopharyngeal swab (NPS) test results and severe SARS-CoV-2 infection.
      Results The data from 6,873 participants (mean age 47.9 ± 14.1 years, 65.8% females) who had a known result from an NPS test were analyzed. According to the multivariate analysis, females had lower odds of a positive result from the NPS test (aOR 0.75, 95%CI 0.66–0.85) and of having a severe infection (aOR 0.46, 95%CI 0.37–0.57) than did their male counterparts. These differences were greater with decreasing age in both sexes. In addition, females aged ≥60 years receiving HT (N = 2,153, 47.6%) had a 46% lower probability of having a positive NPS test (aOR 0.54, 95%CI 0.36–0.80) than their same-aged peers who had never used HT; there were no differences in the younger age groups with respect to HT status.
      Conclusion Female sex was associated with an age-dependent lower risk of having a severe SARS-CoV-2 infection than their male counterparts. Age seemed to modify the relationship between HT status and infection: while the two were not related among younger participants, it was negative in the older ones. Future prospective studies are needed to elucidate the potential protective role sex hormones may play.
      Trial registration ClinicalTrials.gov NCT04471701.

      Keywords

      Abbreviations:

      coronavirus disease (COVID-19), severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), angiotensin-converting enzyme 2 (ACE 2), Italian national epidemiological survey on COVID-19 (EPICOVID19), nasopharyngeal swab (NPS), EU general data protection regulation (EU-GDPR), hormonal therapy (HT), standard deviation (SD), odds ratio (OR), 95% confidence interval (95%CI)

      Introduction

      Since its onset, the coronavirus disease (COVID-19) pandemic caused by the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) has been characterized by marked sex differences [
      • Guan W.J.
      • Ni Z.Y.
      • Hu Y.
      • et al.
      Clinical characteristics of coronavirus disease 2019 in China.
      ]. Although epidemiological evidence collected early on indicated that males have a higher infection rate than females [
      • Stall N.M.
      • Wu W.
      • Lapointe-Shaw L.
      • Fisman D.N.
      • Giannakeas V.
      • Hillmer M.P.
      • Rochon P.A.
      Sex- and age-specific differences in COVID19 testing, cases, and outcomes: a population-wide study in Ontario, Canada.
      ,
      • Zhao S.
      • Cao P.
      • Chong M.K.C.
      • Gao D.
      • Lou Y.
      • Ran J.
      • Wang K.
      • Wang W.
      • Yang L.
      • He D.
      • Wang M.H.
      COVID-19 and gender-specific difference: analysis of public surveillance data in Hong Kong and Shenzhen, China, from January 10 to February 15, 2020.
      ,
      • Abate B.B.
      • Kassie A.M.
      • Kassaw M.W.
      • Aragie T.G.
      • Masresha S.A.
      Sex difference in coronavirus disease (COVID-19): a systematic review and meta-analysis.
      ,
      • Vahidy F.S.
      • Pan A.P.
      • Ahnstedt H.
      • Munshi Y.
      • Choi H.A.
      • Tiruneh Y.
      • Nasir K.
      • Kash B.A.
      • Andrieni J.D.
      • McCullough L.D.
      Sex differences in susceptibility, severity, and outcomes of coronavirus disease 2019: cross-sectional analysis from a diverse US metropolitan area.
      ], sex-disaggregated data collected by the Global Health 50/50 research initiative showed that there was no substantial difference in the male:female ratio for SARS-CoV-2 []. Some studies have however shown that COVID-19 male patients in the most affected age group (60 years old and older) have a higher risk of requiring intensive care treatment, worse outcomes, and mortality with respect to their female counterparts [
      • Gebhard C.
      • Regitz-Zagrosek V.
      • Neuhauser H.K.
      • Morgan R.
      • Klein S.L.
      Impact of sex and gender on COVID-19 outcomes in Europe.
      ,
      • Peckham H.
      • de Gruijter N.M.
      • Raine C.
      • et al.
      Male sex identified by global COVID-19 meta-analysis as a risk factor for death and ITU admission.
      ]. This sex disparity is not entirely surprising since other studies have already demonstrated that males of all age groups are more susceptible than females to other respiratory tract infections (e.g. severe acute respiratory syndrome [SARS] and the Middle East respiratory syndrome [MERS]) [
      • Channappanavar R.
      • Fett C.
      • Mack M.
      • Ten Eyck P.P.
      • Meyerholz D.K.
      • Perlman S.
      Sex-based differences in susceptibility to severe acute respiratory syndrome coronavirus infection.
      ,
      • Alghamdi I.G.
      • Hussain I.I.
      • Almalki S.S.
      • Alghamdi M.S.
      • Alghamdi M.M.
      • El-Sheemy M.A.
      The pattern of middle east respiratory syndrome coronavirus in Saudi Arabia: a descriptive epidemiological analysis of data from the Saudi Ministry of Health.
      ].
      Although the mechanisms underlying these differences are not fully understood, it is probable that they involve an interplay between social, behavioral, and biological factors. Gender-related factors, including socioeconomic status, lifestyles (e.g. smoking habit and alcohol drinking), personal hygiene patterns (e.g. handwashing), healthcare-seeking behavior and access to medical assistance, which in turn can affect the risk of developing diseases, may only partially explain sex differences linked to the SARS-CoV-2 infection [
      • Galasso V.
      • Pons V.
      • Profeta P.
      • Becher M.
      • Brouard S.
      • Foucault M.
      Gender differences in COVID-19 attitudes and behavior: panel evidence from eight countries.
      ,
      • Cai H.
      Sex difference and smoking predisposition in patients with Covid-19.
      ].
      From a biological standpoint, it is known that females present an enhanced immune reactivity that makes them both more vulnerable to developing autoimmune diseases as well as more predisposed to mounting an effective immunity to viral infection [
      • Klein S.L.
      • Flanagan K.L.
      Sex differences in immune responses.
      ]. Genetics and sex hormones may influence both the expression of viral receptors and the differential regulation of immune responses [
      • Klein S.L.
      • Flanagan K.L.
      Sex differences in immune responses.
      ]. Indeed, since they seem to be able to modulate the immune and inflammatory responses and the expression of the Angiotensin-converting enzyme (ACE)-2 gene, which binds the SARS-CoV-2 viral spike protein [
      • Bukowska A.
      • Spiller L.
      • Wolke C.
      • Lendeckel U.
      • Weinert S.
      • Hoffmann J.
      • Bornfleth P.
      • Kutschka I.
      • Gardemann A.
      • Isermann B.
      • Goette A.
      Protective regulation of the ACE2/ACE gene expression by estrogen in human atrial tissue from elderly men.
      ], female hormones may have a protective effect against COVID-19 disease [
      • Pinna G.
      Sex and COVID-19: a protective role for reproductive steroids.
      ]. Conversely, androgens may predispose males to a more severe COVID-19 progression [
      • Giagulli V.A.
      • Guastamacchia E.
      • Magrone T.
      • Jirillo E.
      • Lisco G.
      • De Pergola G.
      • Triggiani V.
      Worse progression of COVID-19 in men: is testosterone a key factor?.
      ]. To the best of our knowledge, only a few studies have investigated the sex- and gender differences linked to the SARS-CoV-2 infection [
      • Stall N.M.
      • Wu W.
      • Lapointe-Shaw L.
      • Fisman D.N.
      • Giannakeas V.
      • Hillmer M.P.
      • Rochon P.A.
      Sex- and age-specific differences in COVID19 testing, cases, and outcomes: a population-wide study in Ontario, Canada.
      ,
      • Zhao S.
      • Cao P.
      • Chong M.K.C.
      • Gao D.
      • Lou Y.
      • Ran J.
      • Wang K.
      • Wang W.
      • Yang L.
      • He D.
      • Wang M.H.
      COVID-19 and gender-specific difference: analysis of public surveillance data in Hong Kong and Shenzhen, China, from January 10 to February 15, 2020.
      ,
      • Vahidy F.S.
      • Pan A.P.
      • Ahnstedt H.
      • Munshi Y.
      • Choi H.A.
      • Tiruneh Y.
      • Nasir K.
      • Kash B.A.
      • Andrieni J.D.
      • McCullough L.D.
      Sex differences in susceptibility, severity, and outcomes of coronavirus disease 2019: cross-sectional analysis from a diverse US metropolitan area.
      ] and the potential therapeutic role of sex hormones [

      R. Costeira, K.A. Lee, B. Murray, C. Christiansen et al., Estrogen and COVID-19 symptoms: associations in women from the COVID symptom study. medRxiv preprint doi: 10.1101/2020.07.30.20164921.

      ,
      • Ding T.
      • Zhang J.
      • Wang T.
      • Cui P.
      • Chen Z.
      • Jiang J.
      • Zhou Su
      • Dai J.
      • Wang Bo
      • Yuan S.
      • Ma W.
      • Ma L.
      • Rong Y.
      • Chang J.
      • Miao X.
      • Ma X.
      • Wang S.
      Potential influence of menstrual status and sex hormones on female severe acute respiratory syndrome coronavirus 2 infection: a cross-sectional multicenter study in Wuhan, China.
      ,
      • Lee J.H.
      • Kim Y.C.
      • Cho S.H.
      • Lee J.
      • You S.C.
      • Song Y.G.
      • Won Y.B.
      • Choi Y.S.
      • Park Y.S.
      Effect of sex hormones on coronavirus disease 2019: an analysis of 5,061 laboratory-confirmed cases in South Korea.
      ,
      • Seeland U.
      • Coluzzi F.
      • Simmaco M.
      • Mura C.
      • Bourne P.E.
      • Heiland M.
      • Preissner R.
      • Preissner S.
      Evidence for treatment with estradiol for women with SARS-CoV-2 infection.
      ].
      Considering the health system burden due to medical care and disabling sequelae of COVID-19 [
      • de Sir A
      • Andrenelli E
      • Negrini F
      • Lazzarini SG
      • Patrini M
      • Ceravolo MG
      International Multiproessional Steering Committee of Cochranef Rehabilitation REH-COVER action.
      Rehabilitation and COVID-19: the Cochrane Rehabilitation 2020 rapid living systematic review.
      ], understanding the extent to which sex hormones underlie the sex and gender differences in the severity of the coronavirus disease could have important clinical and public health implications. Data collected by the large web-based Italian National Epidemiological Survey on COVID-19 (EPICOVID19) during the first early wave of the pandemic were thus analyzed with the intent to: (i) examine if, after adjusting for social, clinical, and behavioral factors, males have a higher probability of having a positive SARS-CoV-2 test result and of developing a severe infection compared to females; and (ii) evaluate the role of hormonal therapy (HT) usage in the female participants.

      2. Material and methods

      2.1 Study design, setting and population

      The EPICOVID19 survey consists of a self-selected convenience sample of 198,822 males and females aged 18–100 living in Italy during the first lockdown who filled in a web-based questionnaire between April and May 2020. The study's methodology has been described in detail elsewhere [
      • Adorni F.
      • Prinelli F.
      • Bianchi F.
      • Giacomelli A.
      • Pagani G.
      • Bernacchia D.
      • Rusconi S.
      • Maggi S.
      • Trevisan C.
      • Noale M.
      • Molinaro S.
      • Bastiani L.
      • Fortunato L.
      • Jesuthasan N.
      • Sojic A.
      • Pettenati C.
      • Tavio M.
      • Andreoni M.
      • Mastroianni C.
      • Antonelli Incalzi R.
      • Galli M.
      Self-reported symptoms of SARS-CoV-2 infection in a nonhospitalized population in Italy: cross-sectional study of the EPICOVID19 web-based survey.
      ]. The inclusion criteria were: age ≥18 years; access to devices connected to internet; and providing online consent to participate in the study. Out of the 198,822 participants who filled out the web questionnaire, two subsets were identified for the purposes of the current study: 1) Sample A = the male and female participants who underwent a nasopharyngeal swab (NPS) test with a known result (N = 6,873); 2) Sample B = the female participants who underwent a NPS test with a known result (N = 4,521) (Fig. 1).
      Fig. 1
      Fig. 1Flow-chart of the study samples A and B.

      2.2 Data collection and exposures definition

      As described elsewhere, all the participants of the EPICOVID19 survey were asked to complete an anonymous 38-item questionnaire [
      • Adorni F.
      • Prinelli F.
      • Bianchi F.
      • Giacomelli A.
      • Pagani G.
      • Bernacchia D.
      • Rusconi S.
      • Maggi S.
      • Trevisan C.
      • Noale M.
      • Molinaro S.
      • Bastiani L.
      • Fortunato L.
      • Jesuthasan N.
      • Sojic A.
      • Pettenati C.
      • Tavio M.
      • Andreoni M.
      • Mastroianni C.
      • Antonelli Incalzi R.
      • Galli M.
      Self-reported symptoms of SARS-CoV-2 infection in a nonhospitalized population in Italy: cross-sectional study of the EPICOVID19 web-based survey.
      ]. The questionnaire was designed to gain information about the participant's life and in particular about his/her age, sex, educational level (illiterate or primary school, middle or high school, and university or postgraduate degree), and employment status (employed, student, unemployed, retired, or other). The questionnaire also asked: if the participant was a healthcare professional, about the participant's residence area (in the northern, central, southern or island regions of Italy), if the participant was living with other individuals at high risk of infection, had self-reported diseases, was taking any medicines, had a flu vaccination during the autumn of 2019 and/or an anti-pneumococcal vaccination over the last 12 months. Other items enquired about how much physical activity the participant were engaged in (never or less than 10 min/week, 10 min to two hours and half/week, and more than two hours and half/week), if the participant had a smoking habit (classified as a non-, former, or current smoker status), had any contacts with confirmed/suspected COVID-19 cases, his/her self-perceived health status (bad, adequate, or good), if she/he had exhibited healthcare-seeking behavior (had contacted an emergency number or his/her general practitioner), had any self-reported SARS-CoV-2 related symptoms, pneumonia, had been hospitalized for confirmed/suspected COVID-19, or had performed a NPS test result (positive vs negative) (Annex 1).

      2.3 Exposure to hormonal therapy

      Some items specifically addressing the female participants inquired if they were currently taking or had in the past (for more or less than 5 years) taken any form of HT, including hormonal contraceptives and/or hormone replacement therapy. The responses were divided into three groups: (1) never-HT users (reference group) vs - HT users; (2) past HT users vs current HT users; (3) participants who used HT less than 5 years vs participants who used HT more than 5 years.

      2.4 Main outcomes

      The study's primary outcome measures were: the results of NPS molecular testing and SARS-CoV-2 infection severity identified by combining information on the NPS test and responses regarding symptoms and hospitalization for COVID-19. Each participant was classified as having: (i) a negative NPS test; (ii) an asymptomatic or mild infection (a positive NPS test without symptoms or with at least one COVID-19-like symptom); or, (iii) severe infection (positive NPS test with pneumonia and/or hospitalization for COVID-19) .

      2.5 Statistical analysis

      Age as a continuous variable was summarized in the descriptive analyses as means and standard deviation (SD); the categorical variables were presented as counts and percentages. The participants’ characteristics according to sex (Sample A) and HT usage (Sample B) were compared using the t-test for age and the chi-squared test for all other categorical variables. Multivariate logistic and multinomial regression models were fitted to assess the relationships between sex (Sample A) and HT usage (Sample B) with positive versus negative NPS test results and the infection severity, respectively. Odds ratio (OR) and 95% confidence interval (95% CI) were estimated. The potential confounders of the two models were selected on the basis of theoretical knowledge and empirical criteria (P-value<0.05 in univariate analysis). Interaction terms were included in the model to investigate if age was a moderator of the effects of sex or HT usage on the NPS test and infection severity. The Wald tests were used to assess the age × sex (and HT interactions). When heterogeneity was present, stratum-specific estimates were evaluated. Healthcare professionals were stratified in a sensitivity analysis to exclude any bias due to the selective inclusion of this category in the survey. In oder to evaluate the effect of HT usage among post-menopausal females, a supplementary analysis focusing only on females aged 50 years or older has been performed. Sample size calculation is available as Supplementary material (S1). All statistical analyses were performed using Stata 15.0 version (StataCorp LP, College Station, Texas, USA) and a two-sided P-value <0.05 was considered statistically significant.

      3. Results

      3.1 Characteristics of the study samples

      The distributions of the participants’ characteristics according to sex (Sample A) and by HT usage (Sample B) are shown in Table 1. The mean age of the sample population that underwent NPS testing was 47.9 ± 14.1 years; 65.8% were females, 24.8% had a positive result and 7.1% developed a severe disease (sample A). With respect to the males, the females were significantly younger and were more likely to have self-reported headaches, heart palpitations, gastrointestinal disturbances, conjunctivitis, and sore throat/rhinorrhoea. The females had a lower rate of NPS positive test results, hospitalizations, and severe COVID-19 infection compared with males. Nearly half (47.6%) of the females in Sample B were HT users. The never-HT users were older, had more comorbidities, and were more likely to have a positive NPS test result with respect to their HT user counterparts.
      Table 1Socio-demographic and behavioral characteristics of the study participants according to sex and HT usage in Sample A (N = 6873) and Sample B (females, N = 4521).
      Sample A (N = 6,873)Sample B (N = 4,521)
      CharacteristicsOverall N = 6,873Males N = 2,352 (34.2%)Females N = 4,521 (65.8%)P-valueNever-HT users N = 2,368 (52.4%)HT users N = 2,153 (47.6%)P-value
      Socio-demographic a
      Age, years (mean, SD)47.914.149.514.147.114.1< 0.00149.915.344.111.9< 0.001
      Age classes< 0.001< 0.001
      < 503,69353.71,12747.92,56656.81,19050.21,37663.9
      50–591,85627.059625.31,26027.968929.157126.5
      ≥601,32419.362926.769515.448920.72069.6
      Low educational level5307.71516.43798.40.01529112.3884.1<0.001
      Retired4666.82159.12515.60.0012008.5512.4<0.001
      Healthcare professionals3,47450.698742.02,48755.0< 0.0011,24352.51,24457.8< 0.001
      Northern area of residence5,16975.21,72273.23,44776.2< 0.00131,80276.11,64576.40.904
      Personal characteristics
      Sedentary habits2,09430.561025.91,48432.8< 0.00188737.559727.7< 0.001
      Current smokers1,06815.534614.772217.9< 0.00134914.737317.3<0.001
      Good self-perceived health status5,51880.3194182.53,57779.1< 0.0011,80276.11,77582.4< 0.001
      Co-habitants at risk
      Older people or anyone with immunocompromising or chronic disease conditions.
      1,37720.040217.197521.6< 0.00159825.337717.5< 0.001
      Contact with COVID-19 cases*
      Suspected/confirmed.
      4,87470.91,60668.33,26872.30.0011,72973.01,53971.50.250
      Healthcare seeking
      Contact the emergency number or the general practitioner.
      2,97743.31,04244.31,93542.80.2331,02743.490842.20.417
      Medical conditions
      Lung diseases5538.12068.83477.70.1171978.31507.00.088
      Heart diseases and/or use of drugs3064.51295.51773.90.0032179.2803.7<0.001
      Hypertension and/or use of drugs1,20217.558725.061513.6< 0.00140217.02139.9< 0.001
      Oncological diseases2.223.2793.41433.20.663773.3663.10.721
      Liver diseases610.9241.0370.80.397311.360.3<0.001
      Renal diseases781.1361.5420.90.025311.3110.50.005
      Metabolic diseases and/or use of drugs3755.51355.72405.30.45527111.41527.1<0.001
      Depression/anxiety and/or use of drugs80111.72209.458112.9< 0.00137615.92059.5< 0.001
      Immune system diseases6709.81094.656112.4< 0.00130512.925611.90.313
      Flu shot vaccination2,31333.787637.21,43731.8< 0.00177332.666430.80.194
      Anti-pneumococcal vaccination3314.81436.11884.2< 0.0011174.9713.30.006
      Oncological drugs721.1231.0491.10.682170.7321.50.013
      Corticosteroids1602.3592.51012.20.474502.1512.40.559
      Thyroid drugs5688.3753.249310.9< 0.00126811.322510.50.350
      Anti-inflammatory drugs4336.31114.73227.1< 0.0011797.61436.60.231
      COVID-19 related variables
      No symptoms1,79226.162026.41,17225.90.69565127.552124.20.012
      Fever1,89827.676232.41,13625.1< 0.00159925.353724.90.784
      Headache2,56237.371530.41,84740.9< 0.00189437.895344.3<0.001
      Muscle/bone pain2,38134.679233.71,58935.20.22379733.779236.80.028
      Olfactory and taste disorders1,45021.149721.195321.10.96047520.147822.20.078
      Shortness of breath1,03615.137315.966314.70.18931913.534416.00.017
      Chest pain96514.030813.165714.50.10430913.134816.20.003
      Heart palpitations87612.82058.767114.8< 0.00131613.335516.50.003
      Gastrointestinal disturbances1,92928.158825.01,34129.7< 0.00164327.269832.4<0.001
      Conjunctivitis82112.025410.856712.50.03530012.726712.40.786
      Sore throat/rhinorrhoea2,53136.875332.0177839.3< 0.00188637.489241.40.006
      Cough2,37134.583935.7153233.90.14075331.877936.20.002
      Pneumonia5578.128912.32685.9< 0.0011375.81316.10.671
      Hospitalized for COVID-195287.728712.22415.3< 0.0011385.81034.80.119
      NPS test positive result1,70224.867728.8102522.7< 0.00159525.143020.0< 0.001
      Infection severity
      Asymptomatic or mild infection (positive NPS test without symptoms or with at least one COVID-19 like symptom) and severe infection (positive NPS test with pneumonia and/or hospitalization for COVID-19).
      < 0.001< 0.001
      Negative NPS test5,17175.2167571.23,49677.3177374.91,72380.0
      Asymptomatic or mild1,21417.72224.980317.847119.933215.4
      Severe4887.126611.32224.91245.2984.6
      ° Older people or anyone with immunocompromising or chronic disease conditions.
      low asterisk Suspected/confirmed.
      Contact the emergency number or the general practitioner.
      § Asymptomatic or mild infection (positive NPS test without symptoms or with at least one COVID-19 like symptom) and severe infection (positive NPS test with pneumonia and/or hospitalization for COVID-19).

      3.2 Association analyses

      Table 2 shows the logistic regression results considering a positive NPS test as the outcome for the whole sample (Sample A). With respect to their male counterparts, being female was inversely associated with the odds of having a positive NPS test (OR 0.75, 95%CI 0.66–0.85). Since a significant interaction between sex and age was found, six levels of indicator variables of males/females aged ≥60, males/females 50–59, and males/females <50 were created. With respect to the males aged ≥60, a pattern across age classes was noted in both sexes; the males younger than 50, the females between 50 and 59, and the females older than 50 showed a significantly lower probability of testing positive by 39%, 44%, and 52%, respectively. Table 3 reports the association between sex and the probability of getting a severe infection. Overall, females had a lower probability of having a severe infection (aOR 0.46, 95%CI 0.37–0.57) with respect to their male counterparts. Considering the older males with no infection as the reference group, data analysis uncovered that the males aged over 50, the females 60 and older, the females between 50 and 59, and the females younger than 50 had, respectively, a significant risk reduction of 66%, 42%, 67%, and 85%.
      Table 2Odds ratios of positive NPS molecular test according to sex and age groups in Sample A (N = 6873).
      TotalNegativePositiveAdjusted
      Sample A°N = 6,873N = 5,171%=75.2N = 1,702%=24.8aOR (95% CI)
      Males235234.2167532.467739.81 (ref.)
      Females4,52165.83,49667.61,02560.20.75 (0.66–0.85)
      SEX X AGE
      Males aged≥606299,33917.623814.01 (ref.)
      Males aged 50–595968,84118.018510.90.87 (0.66–1.14)
      Males aged <501,12716,687316.925414.90.61 (0.47–0.79)
      Females aged≥6069510,24288.326715.70.86 (0.67–1.11)
      Females aged 50–591,26018,697718.928316.60.56 (0.44–0.72)
      Females aged <502,56637,8209140.447527.90.48 (0.38–0.61)
      Males reference category.
      Model adjusted for age°, education, employment status, area of residence, healthcare professionals, physical activity, smoking status, living with at risk co-habitants, contact with COVID-19 cases, heart diseases, hypertension, depression, renal diseases, immune system disorders, flu and anti-pneumococcal vaccine, thyroid drugs, anti-inflammatory drugs, and self-perceived health status.
      Table 3Odds ratios of infection severity according to sex and age groups in Sample A (N = 6873).
      Sample A°Negative NPS test N = 5,171 (75.2)Asymptomatic or mild infection N = 1,214 (17.7)Severe infection N = 488 (7.1)
      N (%)OR(95% CI)N (%)OR(95% CI)N (%)
      Males1,675 (34.4)1 (ref.)

      411 (33.9)1 (ref.)

      266 (54.5)
      Females3,496 (67.6)0.920.79–1.06803 (66.1)0.460.37–0.57222 (45.5)
      SEX X AGE
      Males aged≥60391 (7.6)

      1 (ref.)

      120 (9.9)1 (ref.)

      118 (24.2)
      Males aged 50–59411 (8.0)0.960.69–1.34103 (8.5)0.740.51–1.0882 (16.8)
      Males aged <50873 (16.9)0.850.63–1.16188 (15.5)0.340.23–0.5066 (13.5)
      Females aged≥60428 (8.3)1.150.87–1.53187 (15.4)0.580.41–0.8280 (16.4)
      Females aged 50–59977 (18.9)0.770.57–1.04208 (17.3)0.330.22–0.4875 (15.4)
      Females aged <502,091 (40.4)0.770.58–1.02408 (33.6)0.150.10–0.2267 (13.7)
      Males and negative NPS test as reference category. Asymptomatic or mild infection (positive NPS test without symptoms or with at least one COVID-19 like symptom) and severe infection (positive NPS test with pneumonia and/or hospitalization for COVID-19).
      Model adjusted for age°, education, employment status, area of residence, healthcare professionals, physical activity, smoking status, living with at risk co-habitants, contact with COVID-19 cases, heart diseases, hypertension, depression, renal diseases, immune system disorders, flu and anti-pneumococcal vaccine, thyroid drugs, anti-inflammatory drugs, and self-perceived health status.
      Tables 4 and 5 show that the aORs of having a positive NPS test and of developing a severe infection in the Sample B population significantly decrease with decreasing age. Although no association between HT usage and SARS-CoV-2 infection was found, a statistically significant interaction was observed between age and HT (P-value <0.001). When age class was combined with HT, it was found that, with respect to the never-HT users who were 60 or older, the HT users of the same age class had a 46% reduced odds of receiving a positive result (aOR 0.54). Similar results were found for the younger age classes, irrespective of HT usage. With respect to the never-HT users who were 60 or older, the HT users of the same age class had a lower probability of having an asymptomatic or mild infection (aOR 0.37, 95%CI 0.22–0.61); both the never- and the HT users in the younger age groups (50–59 and <50 years) showed a lower odds of severity and there was no association with severe infection severity. A similar pattern of reduced association in the females who were 60 and over was observed when the HT status and duration were considered (Fig. 2 and Supplementary S23). The associations between the female sex and a positive NPS test result and severity were less pronounced in the healthcare professionals (Supplementary S4). Supplementary material S5–S8 shows the associations between HT usage in females older than 50 years (post-menopausal). With respect to the never-HT users the association between HT users and NPS status and infection severity, even considering the HT status and duration, persisted although the effect was in part attenuated. This is partly due to the differences that the oldest age classes may have when compared to the youngest one, especially with respect to SARS-CoV-2 infection.
      Table 4Adjusted odds ratios of positive NPS test according to age and HT usage in Sample B (females, N = 4521).
      Model BTotal N = 4,521Negative N = 3,496%=77.3Positive N = 1,025%=22.7Adjusted OR (95%CI)
      AGE
      Females aged≥6026715,442812.226726.11 (ref)
      Females aged 50–6028327,997728.028327.60.69 (0.53–0.90)
      Females aged <5047556,8209159.847546.30.62 (0.48–0.80)
      HORMONAL THERAPY°N%N%
      Never-HT users59552,4177350.759558.11 (ref)
      HT users43047,6172349.343042.00.95 (0.81–1.10)
      AGE X HORMONAL THERAPY
      Females never-HT users aged≥6021910,82707.721921.41 (ref)
      Females HT users aged ≥60484,61584.5484.70.54 (0.36–0.80)
      Females never-HT users aged 50–6015915,253015.215915.50.57 (0.42–0.77)
      Females HT users aged 50–6012412,644712.812412.10.57 (0.41–0.79)
      Females never-HT users aged <5021726,397327.821721.20.50 (0.37–0.67)
      Females HT users aged <5025830,4111832.025825.20.52 (0.38–0.69)
      Never-HT users as reference category.
      Model adjusted for age°, education, employment status, area of residence, healthcare professionals, physical activity, smoking status, living with at risk co-habitants, contact with COVID-19 cases, heart diseases, depression, liver and metabolic diseases, flu and anti-pneumococcal vaccine, anti-inflammatory and oncological drugs.
      Table 5Adjusted odds ratios of infection severity according to age and HT usage in Sample B (females, N = 4521).
      Negative NPS test N = 3,496 (77.3)Asymptomatic or mild infection N = 803 (17.8)Severe infection N = 222 (4.9)
      N (%)OR(95% CI)N (%)OR(95% CI)N (%)
      AGE
      Females aged≥60428 (12.2)1 (ref.)187 (23.3)1 (ref)80 (36.0)
      Females aged 50–60977 (28.0)0.750.56–1.01208 (25.9)0.530.35–0.8275 (33.8)
      Females aged <502,091 (59.8)0.800.60–1.06408 (50.8)0.240.16–0.3867 (30.2)
      HORMONAL THERAPY°
      Never-HT users1,773 (50.7)1 (ref.)471 (58.7)1 (ref)124 (55.9)
      HT users1,723 (49.3)0.890.76–1.06332 (41.3)1.190.88–1.6298 (44.1)
      AGE X HORMONAL THERAPY
      Females never-HT users aged≥60270 (7.7)1 (ref.)163 (20.3)1 (ref.)56 (25.2)
      Females HT users aged ≥60158 (4.5)0.370.22–0.6124 (3.0)1.040.58–1.8624 (10.8)
      Females never-HT users aged 50–60530 (15.2)0.590.42–0.84121 (15.1)0.470.28–0.8138 (17.1)
      Females HT users aged 50–60447 (12.8)0.550.38–0.7987 (10.8)0.630.36–1.0937 (16.7)
      Females never-HT users aged <50973 (27.8)0.590.43–0.82187 (23.3)0.240.14–0.4330 (13.5)
      Females HT users aged <501,118 (32.0)0.610.44–0.85221 (27.5)0.250.15–0.4437 (16.7)
      Negative NPS test and never-HT users as reference category. Asymptomatic or mild infection (positive NPS test without symptoms or with at least one COVID-19 like symptom) and severe infection (positive NPS test with pneumonia and/or hospitalization for COVID-19).
      Model adjusted for age°, education, employment status, area of residence, healthcare professionals, physical activity, smoking status, living area, living with at risk co-habitants, contact with COVID-19 cases, heart diseases, depression, liver and metabolic diseases, flu and anti-pneumococcal vaccine, anti-inflammatory and oncological drugs.
      Fig. 2
      Fig. 2Model adjusted for age, education, employment status, area of residence, healthcare professionals, physical activity, smoking status, living area, living with at risk co-habitants, contact with COVID-19 cases, heart diseases, depression, liver and metabolic diseases, flu and anti-pneumococcal vaccine, anti-inflammatory and oncological drugs.

      4. Discussion

      4.1 Main findings

      The data were collected from a large web-based survey of an adult population during the first wave of the COVID-19 pandemic in Italy. After adjusting for several socio-demographic, clinical, and behavioral factors, data analysis showed that the female sex was associated with a 25% lower probability of a positive NPS test result and a 46% lower risk of developing a severe infection; the strength of these associations increased as age decreased. These findings are in line with other epidemiological data gathered from the Italian population performed between March and August 2020 showing that males were more likely than females to test positive [
      • Stall N.M.
      • Wu W.
      • Lapointe-Shaw L.
      • Fisman D.N.
      • Giannakeas V.
      • Hillmer M.P.
      • Rochon P.A.
      Sex- and age-specific differences in COVID19 testing, cases, and outcomes: a population-wide study in Ontario, Canada.
      ,
      • Zhao S.
      • Cao P.
      • Chong M.K.C.
      • Gao D.
      • Lou Y.
      • Ran J.
      • Wang K.
      • Wang W.
      • Yang L.
      • He D.
      • Wang M.H.
      COVID-19 and gender-specific difference: analysis of public surveillance data in Hong Kong and Shenzhen, China, from January 10 to February 15, 2020.
      ,
      • Vahidy F.S.
      • Pan A.P.
      • Ahnstedt H.
      • Munshi Y.
      • Choi H.A.
      • Tiruneh Y.
      • Nasir K.
      • Kash B.A.
      • Andrieni J.D.
      • McCullough L.D.
      Sex differences in susceptibility, severity, and outcomes of coronavirus disease 2019: cross-sectional analysis from a diverse US metropolitan area.
      ] and to manifest severe disease leading to an increased risk of COVID-19-related hospitalization, intensive care unit admission, and death [
      • Stall N.M.
      • Wu W.
      • Lapointe-Shaw L.
      • Fisman D.N.
      • Giannakeas V.
      • Hillmer M.P.
      • Rochon P.A.
      Sex- and age-specific differences in COVID19 testing, cases, and outcomes: a population-wide study in Ontario, Canada.
      ,
      • Vahidy F.S.
      • Pan A.P.
      • Ahnstedt H.
      • Munshi Y.
      • Choi H.A.
      • Tiruneh Y.
      • Nasir K.
      • Kash B.A.
      • Andrieni J.D.
      • McCullough L.D.
      Sex differences in susceptibility, severity, and outcomes of coronavirus disease 2019: cross-sectional analysis from a diverse US metropolitan area.
      ,
      • Peckham H.
      • de Gruijter N.M.
      • Raine C.
      • et al.
      Male sex identified by global COVID-19 meta-analysis as a risk factor for death and ITU admission.
      ]. A systematic review and meta-analysis of 57 studies performed between December 2019 and April 2020 uncovered a pooled prevalence of COVID-19 confirmed cases in males and females of 55% and 45%, respectively [
      • Abate B.B.
      • Kassie A.M.
      • Kassaw M.W.
      • Aragie T.G.
      • Masresha S.A.
      Sex difference in coronavirus disease (COVID-19): a systematic review and meta-analysis.
      ]; these findings were similar to those reported by the Italian National Institute of Health (ISS) at the end of March 2020 [

      Istituto Superiore di Sanità. Epidemia Covid-19. Aggiornamento Nazionale. 2020. Disponibile all'indirizzo: https://www.epicentro.iss.it/coronavirus/bollettino/Infografica_31marzo%20ITA.pdf.

      ]. The Global 50/50 initiative [] and some reviews and meta-analyses focusing on COVID-19 cases worldwide published after that date, however, reported no difference in the male:female ratio of individuals infected with SARS-CoV-2, although the former face a higher risk of hospitalization, intensive care unit admission, and death with respect to their female counterparts [
      • Peckham H.
      • de Gruijter N.M.
      • Raine C.
      • et al.
      Male sex identified by global COVID-19 meta-analysis as a risk factor for death and ITU admission.
      ,
      • Pérez-López F.R.
      • Tajada M.
      • Savirón-Cornudella R.
      • Sánchez-Prieto M.
      • Chedraui P.
      • Terán E.
      Coronavirus disease 2019 and gender-related mortality in European countries: a meta-analysis.
      ,
      • Rozenberga S.
      • Vandrommea J.
      • Martin C.
      Are we equal in adversity? Does Covid-19 affect women and men differently?.
      ]. A seroprevalence systematic review and meta-analysis likewise reported no substantial sex difference in individuals infected with SARS-CoV-2 [
      • Rostami A.
      • Sepidarkish M.
      • Leeflang Mariska M.G.
      • Riahi S.M.
      • Shiadeh M.N.
      • Esfandyari S.
      • Mokdad A.H.
      • Hotez P.J.
      • Gasser R.B
      SARS-CoV-2 seroprevalence worldwide: a systematic review and meta-analysis.
      ]. Our cross-sectional data were collected during the early stages of the epidemic in Italy (April 2020) when only individuals experiencing symptoms severe enough to require medical attention underwent NPS testing. The fact that there were more males with positive NPS results may have distorted our results just as the fact there were more female than male participants that underwent the NPS tests possibly because of the high percentage of female healthcare professionals workers.
      The male participants of our survey were older and more frequently reported comorbidities and severe COVID-19 symptoms such as fever and pneumonia leading to worse disease progression. The females instead were more likely to be paucisymptomatic, presenting with atypical symptoms characterized by sore throat/rhinorrhea, gastrointestinal disturbances, headache, conjunctivitis, and palpitations that were associated with less severe outcomes, as previously documented [
      • Lechien J.R.
      • Chiesa-Estomba C.M.
      • Place S.
      • Van Laethem Y.
      • Cabaraux P.
      • Mat Q.
      • Huet K.
      • Plzak J.
      • Horoi M.
      • Hans S.
      • Barillari M.R
      • Cammaroto G.
      • Fakhry N.
      • Martiny D.
      • Ayad T.
      • Jouffe L.
      • Hopkins C.
      • Saussez S.
      COVID-19 task force of YO-IFOS
      Clinical and epidemiological characteristics of 1420 European patients with mild-to-moderate coronavirus disease 2019.
      ]. As the age-dependent sex disparities in our population persisted after controlling for comorbidities and high-risk behaviors, they seemed to pointing in the direction of possible biological-related explanations. Sex differences in immune responses throughout the lifecourse are influenced by both age and reproductive status. In fact, it has been reported that females tend to be less susceptible to infections than their male counterparts, and it has been hypothesized that sex steroids contribute to the differential regulation of immune responses between sexes [
      • Klein S.L.
      • Flanagan K.L.
      Sex differences in immune responses.
      ]. Both oestrogens and progesterone act by suppressing the production of proinflammatory cytokines associated with the COVID-19 cytokine storm and by enhancing the anti-inflammatory cytokines [
      • Mauvais-Jarvis F.
      • Klein S.L.
      • Levin E.R.
      Estradiol, progesterone, immunomodulation, and COVID-19 outcomes.
      ]. Furthermore, it has been reported that ACE2, which is located on the X chromosome, was significantly downregulated after binding viral Spike protein with consequent reduced ACE2 expression in the lung resulting in severe acute respiratory failure [
      • Hanff T.C.
      • Harhay M.O.
      • Brown T.S.
      • Cohen J.B.
      • Mohareb A.M.
      Is there an association between COVID-19 mortality and the renin-angiotensin system? A call for epidemiologic investigations.
      ]. High estrogen concentrations might up-regulate ACE2 expression leading to an over-expression of ACE2 in females protecting them against viral entry. Low oestrogens levels in males may instead contribute to higher disease susceptibility and death rates [
      • Foresta C.
      • Rocca M.S.
      • Di Nisio A.
      Gender susceptibility to COVID-19: a review of the putative role of sex hormones and X chromosome.
      ]. Androgens could also promote the transcription of the TMPRSS2 gene facilitating viral entry into the cells [
      • Giagulli V.A.
      • Guastamacchia E.
      • Magrone T.
      • Jirillo E.
      • Lisco G.
      • De Pergola G.
      • Triggiani V.
      Worse progression of COVID-19 in men: is testosterone a key factor?.
      ] and decreasing the antibody response to viral infections [
      • Klein S.L.
      • Flanagan K.L.
      Sex differences in immune responses.
      ]. Interestingly, patients with prostate cancer undergoing treatment with androgen deprivation therapy were less likely to develop severe COVID-19 with respect to their non-treated counterparts. Moreover, the low levels of androgens in females may reduce TMPRSS2 expression, further protecting them against the SARS-CoV-2 infection [
      • Montopoli M.
      • Zumerle S.
      • Vettor R.
      • Rugge M.
      • Zorzi M.
      • Catapano C.V.
      • Carbone G.M.
      • Cavalli A.
      • Pagano F.
      • Ragazzi E.
      • Prayer-Galetti T.
      • Alimonti A.
      Androgen-deprivation therapies for prostate cancer and risk of infection by SARS-CoV-2: a population-based study (n= 4532).
      ].
      Our analysis also uncovered that the older females (with natural estrogen deficiency) who were currently receiving or had received HT in the past had a significantly lower odds of getting the infection with respect to their same age counterparts who had never used HT (aOR 0.54). The odds was similar for younger females irrespective of HT use, possibly because of their higher endogenous estrogen levels. These associations seemed to persist independently of the current or past usage of HT or its duration. As females age, the hormone's potential protective effect is attenuated because menopause causes a drastic decline in the natural estrogen levels and affects B and T cells causing post-menopausal women to be more prone to chronic and infectious diseases. In addition, post-menopausal females are reported to have higher levels of pro-inflammatory cytokines that seem to be reduced by HT [
      • Straub R.H.
      The complex role of estrogens in inflammation.
      ].
      Our preliminary data support these hypotheses and agree with the findings of other studies. A study based on an analysis of electronic health records from a large international COVID-19 patient cohort (N = 68,466) recently reported that although estradiol treatment in pre-menopausal females had no effect, fatality risk for post-menopausal females receiving estradiol therapy was reduced by more than 50% [
      • Seeland U.
      • Coluzzi F.
      • Simmaco M.
      • Mura C.
      • Bourne P.E.
      • Heiland M.
      • Preissner R.
      • Preissner S.
      Evidence for treatment with estradiol for women with SARS-CoV-2 infection.
      ]. Menopause was found to be an independent risk factor for COVID-19 in a cross-sectional study investigating 1902 COVID-19 female patients, while anti-Mullerian hormone and estradiol appeared to be potential protective factors [
      • Ding T.
      • Zhang J.
      • Wang T.
      • Cui P.
      • Chen Z.
      • Jiang J.
      • Zhou Su
      • Dai J.
      • Wang Bo
      • Yuan S.
      • Ma W.
      • Ma L.
      • Rong Y.
      • Chang J.
      • Miao X.
      • Ma X.
      • Wang S.
      Potential influence of menstrual status and sex hormones on female severe acute respiratory syndrome coronavirus 2 infection: a cross-sectional multicenter study in Wuhan, China.
      ]. A study investigating 152,637 female users of a COVID Symptom Tracker Application in the UK reported that postmenopausal femlaes had a higher rate of predicted COVID while pre-menopausal females taking the contraceptive pill had a significantly lower rate of predicted COVID-19 and hospitalizations with respect to their post-menopausal counterparts [

      R. Costeira, K.A. Lee, B. Murray, C. Christiansen et al., Estrogen and COVID-19 symptoms: associations in women from the COVID symptom study. medRxiv preprint doi: 10.1101/2020.07.30.20164921.

      ]. The authors of a study carried out in South Korea did not find any association between females who had been taking HT over the past year and morbidity and clinical outcomes of COVID-19 [
      • Lee J.H.
      • Kim Y.C.
      • Cho S.H.
      • Lee J.
      • You S.C.
      • Song Y.G.
      • Won Y.B.
      • Choi Y.S.
      • Park Y.S.
      Effect of sex hormones on coronavirus disease 2019: an analysis of 5,061 laboratory-confirmed cases in South Korea.
      ]. Although it is impossible to exclude that post-menopausal females are more susceptible to SARS-COV-2 infection due to age-related factors and comorbidities rather than to lower estrogen levels, our data suggest that HT may play a protective role against COVID-19 in older females. At present, several clinical trials investigating the potential effect of hormonal therapies (i.e. the selective estrogen receptor modulator, raloxifene) on the infection are currently underway [
      • Hong S.
      • Chang J.
      • Jeong K.
      • Lee W.
      Raloxifene as a treatment option for viral infections.
      ].

      4.2 Limitations and strengths

      These findings need to be interpreted cautiously since the sample was self-selected and cannot be considered entirely representative of the Italian population. The study's generalizability is also limited because restricted to relatively young, female, highly educated, healthcare workers, health-conscious individuals, while a low percentage was made up of participants severely ill with COVID-19 . Furthermore, this study suffers from sampling bias because of the exclusion of individuals who were hospitalized for severe infection or have died because of COVID-19. Other investigators have pointed out that a non-random sampling based on the availability of testing may have led to restrict the analysis to individuals who have been tested for active SARS-CoV 2-infection with severe symptoms [
      • Griffith G.J.
      • Morris T.T.
      • Tudball M.J.
      • et al.
      Collider bias undermines our understanding of COVID-19 disease risk and severity.
      ]. The cross-sectional nature of the study did not enable cause-and-effect inference. The study may be subjected to the “healthy women effect”, females taking HT were more educated, had higher socioeconomic status, a healthier lifestyle, and a lower risk of CVD and metabolic diseases, all factors known to be associated with COVID-19, than their non-HT counterpart, although the models have been adjusted for these variables.
      The study may also be prone to a recall bias because data regarding homone-related exposure were self-reported and the questions regarding HT were non-standardized therefore, misunderstanding and internal errors due to the lack of qualified assessment might have affected the study results. It was moreover impossible to analyze oral contraceptive use and hormone replacement therapy separately; in addition, no specific drug type description was collected. Lastly, although we controlled for several potential confounders, other conditions not considered may have affected our results, therefore unmeasured or residual confounders cannot completely ruled out. The study does nevertheless boast a large sample, and its data, which was collected directly from the general population, reflected the geographical distribution of SARS-CoV-2 infection in the country over the observation period [
      • Adorni F.
      • Prinelli F.
      • Bianchi F.
      • Giacomelli A.
      • Pagani G.
      • Bernacchia D.
      • Rusconi S.
      • Maggi S.
      • Trevisan C.
      • Noale M.
      • Molinaro S.
      • Bastiani L.
      • Fortunato L.
      • Jesuthasan N.
      • Sojic A.
      • Pettenati C.
      • Tavio M.
      • Andreoni M.
      • Mastroianni C.
      • Antonelli Incalzi R.
      • Galli M.
      Self-reported symptoms of SARS-CoV-2 infection in a nonhospitalized population in Italy: cross-sectional study of the EPICOVID19 web-based survey.
      ]. In addition, information on the length of time in years that hormonal therapies were being or had been taken made it possible to perform a time-response analysis. To our knowledge, EPICOVID19 is the largest web survey performed in Italy during the first wave of the pandemic; it collected exhaustive information on socio-demographic, medical, behavioral, and environmental factors, which made it possible, to control for many potential confounding factors.

      5. Conclusion

      These preliminary findings based on data from an Italian adult population collected during the first wave of the pandemic uncovered a sex bias for the SARS-CoV-2 infection rate favoring females. Females over 60 years who were receiving or had received hormonal therapy had a lower probability of having a positive NPS test compared with non-user same-age counterparts. Future prospective studies are warranted to explore how hormonal therapy can affect infection vulnerability and long term COVID19-related clinical outcomes.

      Contributors

      Federica Prinelli was project administrator and contributed to the study concept and design, performed the statistical analysis, and drafted the original manuscript.
      Caterina Trevisan contributed to the study concept and design and to the preparation of the manuscript.
      Marianna Noale contributed to the study concept and design and to the preparation of the manuscript.
      Michela Franchini contributed to the editing and revision of the manuscript and to the interpretation of the results.
      Andrea Giacomelli contributed to the editing and revision of the manuscript and to the interpretation of the results.
      Liliana Cori contributed to the editing and revision of the manuscript and to the interpretation of the results.
      Nithiya Jesuthasan contributed to the editing and revision of the manuscript and to the interpretation of the results.
      Raffaele Antonelli Incalzi contributed to the editing and revision of the manuscript and to the interpretation of the results.
      Stefania Maggi contributed to the study concept and design and to the preparation of the manuscript, and supervised the study.
      Fulvio Adorni was project administrator and contributed to the study concept and design and to the preparation of the manuscript, and supervised the study.

      Funding

      This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.

      Ethical approval

      The Ethics Committee of the Istituto Nazionale per le Malattie Infettive IRCCS Lazzaro Spallanzani (protocol No. 70, 12/4/2020) approved the EPICOVID19 study protocol. The participants were requested to give their informed consent when they first accessed the web-based platform. Participation was voluntary, and no compensation was expected for respondents. The planning, conduct, and reporting of the study were in line with the Declaration of Helsinki, as revised in 2013. All data were handled and stored following the EU General Data Protection Regulation (EU-GDPR) 2016/679, and data transfer was safeguarded by encrypting and decrypting data and password protection. The study has been registered in ClinicalTrials.Gov (NCT04471701).

      Provenance and peer review

      This article was not commissioned and was externally peer reviewed.

      Research data (data sharing and collaboration)

      There are no linked research data sets for this paper. Data will be made available on reasonable request to the corresponding author.

      Declaration of competing interest

      The authors declare that they have no competing interests.

      Acknowledgments

      The authors would like to thank all the participants who took part in this study and made it possible and all the collaborators of the EPICOVID19 Working Group°. The authors would like to thank Linda Inverso for editing the English version of the manuscript.

      Member of the EPICOVID19 Working Group° (in alphabetical order)

      • Adorni Fulvio, National Research Council, Institute of Biomedical Technologies, Epidemiology Unit, Via Fratelli Cervi 93, 20090 Segrate (MI), Italy. [email protected]
      • Andreoni Massimo, Infectious Diseases Clinic, Department of System Medicine, Tor Vergata University of Rome, 00133 Rome, Italy. [email protected]
      • Antonelli Incalzi Raffaele, Unit of Geriatrics, Department of Medicine, Biomedical Campus of Rome, via Alvaro del Portillo, 21, 00128 Rome, Italy. [email protected]
      • Bastiani Luca, National Research Council, Institute of Clinical Physiology, Via G. Moruzzi 1, 56124 Pisa (PI), Italy. [email protected]
      • Bianchi Fabrizio, National Research Council, Institute of Clinical Physiology, Via G. Moruzzi 1, 56124 Pisa (PI), Italy. [email protected]
      • Di Bari Mauro, Geriatric Intensive Care Medicine, University of Florence and Azienda Ospedaliero-Universitaria Careggi, Viale Peraccini 18, 50139 Florence, Italy. [email protected]
      • Fortunato Loredana, National Research Council, Institute of Clinical Physiology, Via G. Moruzzi 1, 56124 Pisa (PI), Italy. [email protected]
      • Galli Massimo, Infectious Diseases Unit, Department of Biomedical and Clinical Sciences L. Sacco, Università di Milano, ASST Fatebenefratelli Sacco, 20157 Milan, Italy. [email protected]
      • Giacomelli Andrea, Infectious Diseases Unit, Department of Biomedical and Clinical Sciences L. Sacco, Università di Milano, ASST Fatebenefratelli Sacco, 20157 Milan, Italy. [email protected]
      • Jesuthasan Nithiya, National Research Council, Institute of Biomedical Technologies, Epidemiology Unit, Via Fratelli Cervi 93, 20090 Segrate (MI), Italy. [email protected]
      • Maggi Stefania, National Research Council-Neuroscience Institute, Aging Branch, Via Vincenzo Maria Gallucci 16, 35128 Padova, Italy. [email protected]
      • Mastroianni Claudio, Public Health and Infectious Disease Department, “Sapienza” University, Piazzale Aldo Moro 1, 00185, Rome, Italy. [email protected]
      • Molinaro Sabrina, National Research Council, Institute of Clinical Physiology, Via G. Moruzzi 1, 56124 Pisa (PI), Italy. [email protected]
      • Noale Marianna, National Research Council-Neuroscience Institute, Aging Branch, Via Vincenzo Maria Gallucci 16, 35128 Padova, Italy. [email protected]
      • Pagani Gabriele, Infectious Diseases Unit, Department of Biomedical and Clinical Sciences L. Sacco, Università di Milano, ASST Fatebenefratelli Sacco, 20157 Milan, Italy. [email protected]
      • Pedone Claudio, Unit of Geriatrics, Department of Medicine, Biomedical Campus of Rome, via Alvaro del Portillo, 21, 00128 Rome, Italy. [email protected]
      • Pettenati Carla, National Research Council, Institute of Biomedical Technologies, Via Fratelli Cervi 93, 20090 Segrate (MI), Italy. [email protected]
      • Prinelli Federica, National Research Council, Institute of Biomedical Technologies, Epidemiology Unit, Via Fratelli Cervi 93, 20090 Segrate (MI), Italy. [email protected]
      • Rusconi Stefano, Infectious Diseases Unit, Department of Biomedical and Clinical Sciences L. Sacco, Università di Milano, ASST Fatebenefratelli Sacco, 20157 Milan, Italy. [email protected]
      • Sojic Aleksandra, National Research Council, Institute of Biomedical Technologies, Epidemiology Unit, Via Fratelli Cervi 93, 20090 Segrate (MI), Italy. [email protected]
      • Tavio Marcello, Division of Infectious Diseases, Azienda Ospedaliero Universitaria Ospedali Riuniti, Via Conca 71, Torrette, Ancona, Italy. [email protected]
      • Trevisan Caterina, Geriatric Unit, Department of Medicine (DIMED), University of Padova, Via
      • Giustiniani 2, 35128 Padova, Italy; National Research Council-Neuroscience Institute, Aging Branch, Via Vincenzo Maria Gallucci 16, 35128 Padova, Italy. [email protected]

      Appendix. Supplementary materials

      References

        • Guan W.J.
        • Ni Z.Y.
        • Hu Y.
        • et al.
        Clinical characteristics of coronavirus disease 2019 in China.
        N. Engl. J. Med. 2020; 382: 1708-1720https://doi.org/10.1056/NEJMoa2002032
        • Stall N.M.
        • Wu W.
        • Lapointe-Shaw L.
        • Fisman D.N.
        • Giannakeas V.
        • Hillmer M.P.
        • Rochon P.A.
        Sex- and age-specific differences in COVID19 testing, cases, and outcomes: a population-wide study in Ontario, Canada.
        J. Am. Geriatr. Soc. 2020; 68 (OctEpub 2020 Aug 15. PMID:32743827): 2188-2191https://doi.org/10.1111/jgs.16761
        • Zhao S.
        • Cao P.
        • Chong M.K.C.
        • Gao D.
        • Lou Y.
        • Ran J.
        • Wang K.
        • Wang W.
        • Yang L.
        • He D.
        • Wang M.H.
        COVID-19 and gender-specific difference: analysis of public surveillance data in Hong Kong and Shenzhen, China, from January 10 to February 15, 2020.
        Infect. Control Hosp. Epidemiol. 2020; 41 (JunPMID:32146921; PMCID: PMC7113032): 1-2https://doi.org/10.1017/ice.2020.64
        • Abate B.B.
        • Kassie A.M.
        • Kassaw M.W.
        • Aragie T.G.
        • Masresha S.A.
        Sex difference in coronavirus disease (COVID-19): a systematic review and meta-analysis.
        BMJ Open. 2020; 10 (Oct 6PMID:33028563; PMCID: PMC7539579)e040129https://doi.org/10.1136/bmjopen-2020-040129
        • Vahidy F.S.
        • Pan A.P.
        • Ahnstedt H.
        • Munshi Y.
        • Choi H.A.
        • Tiruneh Y.
        • Nasir K.
        • Kash B.A.
        • Andrieni J.D.
        • McCullough L.D.
        Sex differences in susceptibility, severity, and outcomes of coronavirus disease 2019: cross-sectional analysis from a diverse US metropolitan area.
        PLoS One. 2021; 16 (Jan 13PMID:33439908; PMCID: PMC7806140)e0245556https://doi.org/10.1371/journal.pone.0245556
      1. The sex, gender and COVID-19 project, Men, sex gender and Covid-19. https://globalhealth5050.org/the-sex-gender-and-covid-19-project/men-sex-gender-and-covid-19/.

        • Gebhard C.
        • Regitz-Zagrosek V.
        • Neuhauser H.K.
        • Morgan R.
        • Klein S.L.
        Impact of sex and gender on COVID-19 outcomes in Europe.
        Biol. Sex Differ. 2020; 11 (May 25PMID:32450906; PMCID: PMC7247289): 29https://doi.org/10.1186/s13293-020-00304-9
        • Peckham H.
        • de Gruijter N.M.
        • Raine C.
        • et al.
        Male sex identified by global COVID-19 meta-analysis as a risk factor for death and ITU admission.
        Nat. Commun. 2020; 11: 6317https://doi.org/10.1038/s41467-020-19741-6
        • Channappanavar R.
        • Fett C.
        • Mack M.
        • Ten Eyck P.P.
        • Meyerholz D.K.
        • Perlman S.
        Sex-based differences in susceptibility to severe acute respiratory syndrome coronavirus infection.
        J. Immunol. 2017; 198 (May 15Epub 2017 Apr 3. PMID:28373583; PMCID: PMC5450662): 4046-4053https://doi.org/10.4049/jimmunol.1601896
        • Alghamdi I.G.
        • Hussain I.I.
        • Almalki S.S.
        • Alghamdi M.S.
        • Alghamdi M.M.
        • El-Sheemy M.A.
        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 (Aug 20PMID:25187734PMCID: PMC4149400): 417-423https://doi.org/10.2147/ijgm.s67061
        • Galasso V.
        • Pons V.
        • Profeta P.
        • Becher M.
        • Brouard S.
        • Foucault M.
        Gender differences in COVID-19 attitudes and behavior: panel evidence from eight countries.
        Proc. Natl. Acad. Sci. U. S. A. 2020; 117 (Nov 3Epub 2020 Oct 15. PMID:33060298; PMCID: PMC7959517): 27285-27291https://doi.org/10.1073/pnas.2012520117
        • Cai H.
        Sex difference and smoking predisposition in patients with Covid-19.
        Lancet Respir. Med. 2020; 8 (AprEpub 2020 Mar 11. Erratum in: Lancet Respir Med. 2020 Apr;8(4):e26. PMID:32171067; PMCID: PMC7103991): e20https://doi.org/10.1016/S2213-2600(20)30117-X
        • Klein S.L.
        • Flanagan K.L.
        Sex differences in immune responses.
        Nat. Rev. Immunol. 2016; 16 (OctEpub 2016 Aug 22. PMID:27546235): 626-638https://doi.org/10.1038/nri.2016.90
        • Bukowska A.
        • Spiller L.
        • Wolke C.
        • Lendeckel U.
        • Weinert S.
        • Hoffmann J.
        • Bornfleth P.
        • Kutschka I.
        • Gardemann A.
        • Isermann B.
        • Goette A.
        Protective regulation of the ACE2/ACE gene expression by estrogen in human atrial tissue from elderly men.
        Exp. Biol. Med. (Maywood). 2017; 242 (AugEpub 2017 Jun 29. PMID:28661206; PMCID: PMC5544171): 1412-1423https://doi.org/10.1177/1535370217718808
        • Pinna G.
        Sex and COVID-19: a protective role for reproductive steroids.
        Trends Endocrinol. Metab. 2021; 32 (JanEpub 2020 Nov 9. PMID:33229187; PMCID: PMC7649655): 3-6https://doi.org/10.1016/j.tem.2020.11.004
        • Giagulli V.A.
        • Guastamacchia E.
        • Magrone T.
        • Jirillo E.
        • Lisco G.
        • De Pergola G.
        • Triggiani V.
        Worse progression of COVID-19 in men: is testosterone a key factor?.
        Andrology. 2021; 9 (JanEpub 2020 Jun 28. PMID:32524732; PMCID: PMC7307026): 53-64https://doi.org/10.1111/andr.12836
      2. R. Costeira, K.A. Lee, B. Murray, C. Christiansen et al., Estrogen and COVID-19 symptoms: associations in women from the COVID symptom study. medRxiv preprint doi: 10.1101/2020.07.30.20164921.

        • Ding T.
        • Zhang J.
        • Wang T.
        • Cui P.
        • Chen Z.
        • Jiang J.
        • Zhou Su
        • Dai J.
        • Wang Bo
        • Yuan S.
        • Ma W.
        • Ma L.
        • Rong Y.
        • Chang J.
        • Miao X.
        • Ma X.
        • Wang S.
        Potential influence of menstrual status and sex hormones on female severe acute respiratory syndrome coronavirus 2 infection: a cross-sectional multicenter study in Wuhan, China.
        Clin. Infect. Dis. 2020; (Jul): ciaa1022https://doi.org/10.1093/cid/ciaa1022
        • Lee J.H.
        • Kim Y.C.
        • Cho S.H.
        • Lee J.
        • You S.C.
        • Song Y.G.
        • Won Y.B.
        • Choi Y.S.
        • Park Y.S.
        Effect of sex hormones on coronavirus disease 2019: an analysis of 5,061 laboratory-confirmed cases in South Korea.
        Menopause. 2020; 27 (DecPMID:33003134; PMCID: PMC7709921): 1376-1381https://doi.org/10.1097/GME.0000000000001657
        • Seeland U.
        • Coluzzi F.
        • Simmaco M.
        • Mura C.
        • Bourne P.E.
        • Heiland M.
        • Preissner R.
        • Preissner S.
        Evidence for treatment with estradiol for women with SARS-CoV-2 infection.
        BMC Med. 2020; 18 (Nov 25PMID:33234138; PMCID: PMC7685778): 369https://doi.org/10.1186/s12916-020-01851-z
        • de Sir A
        • Andrenelli E
        • Negrini F
        • Lazzarini SG
        • Patrini M
        • Ceravolo MG
        • International Multiproessional Steering Committee of Cochranef Rehabilitation REH-COVER action.
        Rehabilitation and COVID-19: the Cochrane Rehabilitation 2020 rapid living systematic review.
        Eur J Phys Rehabil Med. 2020; https://doi.org/10.23736/S1973-9087.20.06614-9
        • Adorni F.
        • Prinelli F.
        • Bianchi F.
        • Giacomelli A.
        • Pagani G.
        • Bernacchia D.
        • Rusconi S.
        • Maggi S.
        • Trevisan C.
        • Noale M.
        • Molinaro S.
        • Bastiani L.
        • Fortunato L.
        • Jesuthasan N.
        • Sojic A.
        • Pettenati C.
        • Tavio M.
        • Andreoni M.
        • Mastroianni C.
        • Antonelli Incalzi R.
        • Galli M.
        Self-reported symptoms of SARS-CoV-2 infection in a nonhospitalized population in Italy: cross-sectional study of the EPICOVID19 web-based survey.
        JMIR Public Health Surveill. 2020; 6 (Sep 18PMID:32650305PMCID: 7505691): e21866https://doi.org/10.2196/21866
      3. Istituto Superiore di Sanità. Epidemia Covid-19. Aggiornamento Nazionale. 2020. Disponibile all'indirizzo: https://www.epicentro.iss.it/coronavirus/bollettino/Infografica_31marzo%20ITA.pdf.

        • Pérez-López F.R.
        • Tajada M.
        • Savirón-Cornudella R.
        • Sánchez-Prieto M.
        • Chedraui P.
        • Terán E.
        Coronavirus disease 2019 and gender-related mortality in European countries: a meta-analysis.
        Maturitas. 2020; 141 (NovEpub 2020 Jun 23PMID:33036704; PMCID: PMC7309755): 59-62https://doi.org/10.1016/j.maturitas.2020.06.017
        • Rozenberga S.
        • Vandrommea J.
        • Martin C.
        Are we equal in adversity? Does Covid-19 affect women and men differently?.
        Maturitas. 2020; 138: 62-68https://doi.org/10.1016/j.maturitas.2020.05.009
        • Rostami A.
        • Sepidarkish M.
        • Leeflang Mariska M.G.
        • Riahi S.M.
        • Shiadeh M.N.
        • Esfandyari S.
        • Mokdad A.H.
        • Hotez P.J.
        • Gasser R.B
        SARS-CoV-2 seroprevalence worldwide: a systematic review and meta-analysis.
        Clin. Microbiol. Infect. 2021; 27 (Epub 2020 Oct 24. PMID:33228974; PMCID: PMC7584920): 331e340https://doi.org/10.1016/j.cmi.2020.10.020
        • Lechien J.R.
        • Chiesa-Estomba C.M.
        • Place S.
        • Van Laethem Y.
        • Cabaraux P.
        • Mat Q.
        • Huet K.
        • Plzak J.
        • Horoi M.
        • Hans S.
        • Barillari M.R
        • Cammaroto G.
        • Fakhry N.
        • Martiny D.
        • Ayad T.
        • Jouffe L.
        • Hopkins C.
        • Saussez S.
        • COVID-19 task force of YO-IFOS
        Clinical and epidemiological characteristics of 1420 European patients with mild-to-moderate coronavirus disease 2019.
        J. Intern. Med. 2020; 288 (SepEpub 2020 Jun 17. PMID:32352202; PMCID: PMC7267446): 335-344https://doi.org/10.1111/joim.13089
        • Mauvais-Jarvis F.
        • Klein S.L.
        • Levin E.R.
        Estradiol, progesterone, immunomodulation, and COVID-19 outcomes.
        Endocrinology. 2020; 161 (SepPMID:32730568; PMCID: PMC7438701): bqaa127https://doi.org/10.1210/endocr/bqaa127
        • Hanff T.C.
        • Harhay M.O.
        • Brown T.S.
        • Cohen J.B.
        • Mohareb A.M.
        Is there an association between COVID-19 mortality and the renin-angiotensin system? A call for epidemiologic investigations.
        Clin. Infect. Dis. 2020; 71 (Jul 28PMID:32215613; PMCID: PMC7184340): 870-874https://doi.org/10.1093/cid/ciaa329
        • Foresta C.
        • Rocca M.S.
        • Di Nisio A.
        Gender susceptibility to COVID-19: a review of the putative role of sex hormones and X chromosome.
        J. Endocrinol. Invest. 2021; 44 (MayEpub 2020 Sep 16. PMID:32936429; PMCID: PMC7492232): 951-956https://doi.org/10.1007/s40618-020-01383-6
        • Montopoli M.
        • Zumerle S.
        • Vettor R.
        • Rugge M.
        • Zorzi M.
        • Catapano C.V.
        • Carbone G.M.
        • Cavalli A.
        • Pagano F.
        • Ragazzi E.
        • Prayer-Galetti T.
        • Alimonti A.
        Androgen-deprivation therapies for prostate cancer and risk of infection by SARS-CoV-2: a population-based study (n= 4532).
        Ann. Oncol. 2020; 31 (AugEpub 2020 May 6. PMID:32387456; PMCID: PMC7202813): 1040-1045https://doi.org/10.1016/j.annonc.2020.04.479
        • Straub R.H.
        The complex role of estrogens in inflammation.
        Endocr. Rev. 2007; 28: 521-574
        • Hong S.
        • Chang J.
        • Jeong K.
        • Lee W.
        Raloxifene as a treatment option for viral infections.
        J. Microbiol. 2021; 59 (FebEpub 2021 Feb 1. PMID:33527314; PMCID: PMC7849956): 124-131https://doi.org/10.1007/s12275-021-0617-7
        • Griffith G.J.
        • Morris T.T.
        • Tudball M.J.
        • et al.
        Collider bias undermines our understanding of COVID-19 disease risk and severity.
        Nat. Commun. 2020; 11 (PMID:33184277; PMCID: PMC7665028): 5749https://doi.org/10.1038/s41467-020-19478-2