Introduction

A consistent feature of the ongoing coronavirus disease 2019 (COVID-19) pandemic, caused by the novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)1, is the male bias towards severe disease2,3,4,5,6. Despite this, there are few analyses addressing whether this is a global rather than regional phenomenon. There are limited data that indicate whether the bias towards increased mortality in males is due to an increased proportion of infections in males, or a true representation of more severe disease. This gap in the literature highlights that sex remains an under-appreciated variable when interrogating outcomes in infectious diseases.

To address whether the reported sex-bias is validated by large-scale statistical analysis of global data, we have collected available case data from 90 reports including 46 different countries and 44 US states totalling 3,111,714 infected cases, and present a meta-analysis to investigate sex as a risk factor for SARS-CoV-2 infection, and COVID-19 morbidity and mortality. We demonstrate that while there is no difference in the proportion of males and females infected with SARS-CoV-2, males face higher odds of both intensive therapy unit (ITU) admission and death compared to females. The confirmation of this sex disparity with global data has important implications for the continuing public health response to this pandemic.

Results

Source selection and inclusion

To interrogate the sex-bias in the global COVID-19 pandemic, 107 reports were gathered from across the world, from 1st January 2020 until the 1st June 2020 (Fig. 1). Reports were included if they contained data on the total number of infections by sex, and the severity of disease as measured by admission to ITU and death (n = 92). Most reports (n = 88) were government websites, one report was from the British Medical Journal Global Health Blog website, and a minority comprised published research articles (n = 3), from which data were captured from full text. Importantly, sex disaggregated data for the United States of America (USA) were only available at individual state levels. Data from individual states were therefore analysed as separate regions to account for differences in how the data may have been collected. Reports were excluded if they did not report the overall number of infections by sex (one ITU case series and three mortality case series), or had small case numbers (less than five). Of the remaining 102 reports, multiple early reports originated from China, which were carefully examined for duplication. Following the exclusion of nine reports due to possible duplication, and one report from Canada with a large percentage of cases with unknown sex, a total of 92 reports from 47 different countries were analysed further (Supplementary Data 1)7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,79,80,81,82,83,84,85,86,87,88,89,90,91,92,93,94,95,96,97,98. The 92 reports included three reports from China, the largest of which included data on confirmed cases by sex and mortality by sex, but not ITU admission by sex9. For this reason, the other two Chinese reports10,11 were included in the analysis of ITU admission by sex, but were excluded from the analyses of confirmed cases and mortality by sex as it was likely these two reports overlapped with the larger report of Chinese cases9.

Fig. 1: Study selection.
figure 1

A total of 107 reports were found. Five reports were excluded as they did not report the total number of infections by sex (one intensive therapy unit (ITU) admission only case series and three mortality case series) and one report was excluded as it contained less than five cases. Of the 102 remaining reports, 9 were subsequently excluded due to possible duplication and 1 was excluded due to a large proportion of cases with unknown sex, yielding a total of 92 reports contributing to the analysis. These included three reports from China, which contributed differentially to the analysis. The largest report from China was used for analyses of confirmed cases and mortality by sex with the two other reports from China excluded from those analyses but used for the analysis of ITU admission by sex. This resulted in 90 reports of confirmed cases by sex, eight reports with data on ITU admission by sex and 70 reports with data on mortality by sex. Note that where totals for males and females do not add to the reported totals, this is because sex was unknown for some cases in the original source data.

Final reports and cases contributing to analysis: Ninety reports described the sex of 3,111,714 infected cases in 46 countries and 44 US states. Eight reports, representing 341,571 cases described 12,067 ITU admissions by sex. Seventy-one reports included mortality by sex, and after exclusion of the one possible overlapping report from China, 70 reports describing the sex of 2,751,115 COVID-19 cases and 214,361 related deaths were included in the final analysis. Counts of confirmed male cases and male deaths were calculated for 23 and 13 sources, respectively as these sources reported percentages and not counts. Thirty-four sources included a small proportion of cases of unknown sex and 17 sources included a small proportion of deaths of unknown sex (Supplementary Data 1).

Meta-analyses of infection, ITU and mortality risks: The proportion of male cases with COVID-19 in these reports was exactly half at 0.50 (95% confidence interval (CI) = 0.48,0.51; p = 0.56; n = 3,111,714), demonstrating that males and females have similar numbers of infections (Fig. 2). Male sex was associated with increased odds of ITU admission (odds ratio (OR) = 2.84; 95% CI = 2.06, 3.92; p = 1.86 × 10−10; n = 341,571) (Fig. 3), and increased odds of death (OR = 1.39; 95% CI = 1.31,1.47; p = 5.00 × 10−30; n = 2,751,115) (Fig. 4).

Fig. 2: There is no observed sex difference in the proportion of people with COVID-19.
figure 2

The table summarises the number of confirmed male COVID-19 cases and the total number of COVID-19 cases in 90 reports. The forest plot illustrates the estimated proportion of male cases for each report (grey boxes), with 95% confidence intervals (CI; horizontal black lines). The estimated pooled proportion of male cases (dark grey diamond) was 0.5 (95% CI = 0.48,0.51). A two-sided test confirmed the estimated pooled proportion was not significantly different from 0.5 (p = 0.56), indicating no difference between the proportions of male and female infections. Meta-analysis used a random effects model, which accounted for variance across reports and used the indicated weights for each report.

Fig. 3: Male sex is associated with a significantly increased risk of ITU admission within COVID-19 patients.
figure 3

The table summarises the number of ITU admissions and total number of confirmed COVID-19 cases for each sex for n = 8 reports with complete data on SARS-CoV-2 infections and ITU admissions in males and females. The forest plot represents the estimated odds ratio (OR) for the association of ITU admission with male sex for each report (grey boxes), with 95% confidence intervals (CI; horizontal black lines). The estimated pooled OR (dark grey diamond) was 2.84 (95% CI = 2.06, 3.92). A two-sided test confirmed the estimated pooled OR was significantly different from 1 (p = 1.86e-10). Meta-analysis used a random effects model with individual reports weighted using the indicated weights.

Fig. 4: Male sex is associated with a significantly increased risk of mortality within COVID-19 patients.
figure 4

The table summarises the number of deaths and total number of confirmed COVID-19 cases for each sex for n = 70 reports with complete data on infection and mortality in males and females. The forest plot illustrates the estimated odds ratio (OR) for the association of death with male sex for each report (grey boxes) with 95% confidence intervals (CI; horizontal black lines). The estimated pooled OR (dark grey diamond) was 1.39 (95% CI = 1.31, 1.47). A two-sided test confirmed the estimated pooled OR was significantly different from 1 (p = 5.00e-30). Meta-analysis used a random effects model with individual reports weighted using the indicated weights.

Funnel plots and sensitivity analyses indicated that the results for the estimated proportion of male cases and the estimated OR for ITU admission in men were unlikely to be influenced by reporting bias; however, the estimated OR for mortality in men may be an underestimate (Supplementary Information). Taken together these data highlight that, whilst there is no difference in the proportion of COVID-19 cases between sexes, men have a higher risk of ITU admission (OR 2.84) and death (OR 1.39).

Discussion

Here we present a large-scale global statistical analysis to show that whilst males and females are at equivalent risk of infection, male sex is associated with the development of severe disease as measured by ITU admission (OR = 2.84; 95% CI = 2.06, 3.92; p = 1.86 × 10−10) and death (OR = 1.39; 95% CI = 1.31,1.47; p = 5.00 × 10−30). Despite this notable feature of the pandemic, sex is still not routinely reported in all available regional data.

Sex differences in the prevalence and outcomes of infectious diseases occur at all ages, with an overall higher burden of bacterial, viral, fungal and parasitic infections in human males99,100,101,102. Previous coronavirus outbreaks have demonstrated the same sex bias. The Hong Kong SARS-CoV-1 epidemic showed an age-adjusted relative mortality risk ratio of 1.62 (95% CI = 1.21, 2.16) for males103. During the same outbreak in Singapore, male sex was associated with an odds ratio of 3.10 (95% CI = 1.64, 5.87; p ≤ 0.001) for ITU admission or death104. The Saudi Arabian MERS outbreak in 2013 - 2014 exhibited a case fatality rate of 52% in men and 23% in women105. These data suggest that, whilst socio-economic factors may be influencing some aspects of the pandemic, fundamental differences in the immune response between males and females are likely to be a driving factor behind the significant sex-bias observed in the COVID-19 pandemic.

Sex differences in both the innate and adaptive immune system have been previously reported and may account for the female advantage in COVID-19. Within the adaptive immune system, females have higher numbers of CD4+ T cells106,107,108,109,110,111, more robust CD8+ T cell cytotoxic activity112, and increased B cell production of immunoglobulin compared to males106,113. Women report more severe local and systemic side effects, and produce higher antibody titres in response to the trivalent inactivated seasonal influenza vaccination (TIV)100, as well as to most other pathogen vaccines114. More specifically, females achieve equivalent protective antibody titres to males at half the dose of TIV115, with serum testosterone levels inversely correlating with TIV antibody titres116. Female B cells also produce more antigen-specific IgG in response to TIV117. These findings imply that females have an increased capacity to mount humoral immune responses compared to males, and together with the data from this meta-analysis, may have important implications for the development of vaccination strategies for COVID-19.

Females produce more type 1 interferon (IFN), a potent anti-viral cytokine, upon toll-like receptor 7 sensing of viral RNA than males118,119,120,121,122,123,124, which is important for the early response in COVID-19125. Increased IFN production by females is associated with both sex hormone concentration and the number of X chromosomes present119,121. The X chromosome contains many immune-related genes126, as evidenced by the existence of many X-linked immunodeficiency disorders127. Furthermore X-encoded immune genes may be variably expressed on both alleles in immune cells in females, increasing immune response diversity128,129. Oestradiol offers an advantage against infectious disease by augmenting T cell responses112,130,131,132, increasing antibody production, somatic hyper-mutation and class switching133. Oestradiol also increases abundance of neutrophils134; and monocyte/macrophage cytokine production135.

An early report in pre-print suggested that the declining oestrogen levels in post-menopausal women may be associated with increased inflammatory cytokine production following infection with SARS-CoV-2136. This suggests a potential protective effect of oestradiol against the development of hyperinflammatory immune responses associated with mortality in COVID-19137. In contrast, the male sex hormone testosterone suppresses the immune system: hypo-androgenism is associated with increased inflammatory cytokines, antibody titres, CD4/CD8 ratios, and natural killer cells, and a decrease in regulatory T cells101,138. Interestingly, testosterone-deprivation therapy for prostate cancer has been associated with improved outcomes for COVID-19, suggesting that suppression of the immune response by testosterone, as well as the protective effect of oestrogen, may underlie the observed sex-bias139.

Age-related changes in the immune system are also different between sexes140 and there is a marked association between morbidity/mortality and advanced age in COVID-19. For example, males show an age-related decline in B cells and a trend towards accelerated immune ageing141,142. This may further contribute to the sex bias seen in COVID-19.

Other biological factors may influence the sex-bias observed in this study. Expression of angiotensin converting enzyme 2 (ACE2) receptors – which facilitate SARS-CoV-2 viral entry and human to human transmission143 – is different between the sexes144,145. Oestradiol may influence ACE2 expression146, and the gene for ACE2 is located in the X chromosome147, which may render it susceptible to escaping X-inactivation in women.

Sex-based differences in co-morbidities that are associated with severe COVID-19 may also drive some of the differences observed in this study. However, due to the nature of these high level, publicly available summary data, metadata including age, ethnicity and comorbidities for individual cases are not available. The lack of adjustment for these factors limits our ability to accurately predict the role of sex in disease severity. Notably, there are no marked sex differences in the proportions of adults globally with hypertension (33% of women vs 36% of men)148 or diabetes (9% of women vs 9.6% of men)149, the most common reported comorbidities in hospitalised COVID-19 patients150. Once more data become available, future studies can adjust for additional factors using techniques such as mediation analysis.

Gender-based socio-cultural and behavioural differences could contribute to the sex difference seen in COVID-19 disease severity. Men are more likely to smoke151, although smoking has not emerged as a clear risk factor for severe disease144. Men are less likely to wash their hands with soap after entering a restroom152, and in many cultures, men may be more likely to leave the house and enter crowded areas. Unequal access to healthcare and testing between sexes may skew towards a male bias in infection rates. The data, however, show no difference in the numbers of infected cases between sexes overall, so gender differences in hygiene behaviours and testing are unlikely to explain the sex disparity in disease severity. Regional gender differences in health-seeking behaviours and access to care may predispose men towards access to hospital and ITU admission153,154,155,156. However, the ubiquitous nature of the sex-bias in these data argues for a true biological difference in the response to SARS-CoV-2 between sexes.

These large-scale data demonstrate that although there is no sex difference in the proportion of people infected with SARS-CoV-2, males are at a significantly higher risk of severe disease and death than females. Previous reports describe fundamental differences between sexes in the immune response to infection, which include a more robust antiviral innate interferon response and increased adaptive immunity towards viral antigens in females. In people infected with SARS-CoV-2 these differences are likely to lead to more effective viral control in females, which may contribute to the relatively lower risk of developing severe disease. Although further studies are needed, these data have implications for the clinical management of COVID-19 and highlight the importance of considering sex as a variable in fundamental and clinical research.

Methods

Search strategy

An online search of government websites and published literature was performed by multiple researchers working remotely for regional data reports on COVID-19 cases that included sex as a variable from 01 January 2020 up until 01 June 2020 (Search terms: SARS-CoV-2/COVID-19/case/sex/country/data/death/ICU/ITU). In order to ensure unbiased representation from as many regions as possible, a cross check was done using the list of countries reporting data on ‘Worldometer’157, and an attempt was made to include as many regions reporting sex data as possible. Reports were translated using Google translate if they were not in English.

Data selection, extraction and synthesis

Reports were included if they contained sex as a variable in data describing case number, ITU admission, or mortality. Data were entered directly by individual researchers into an online structured data extraction table. For some sources, counts of male confirmed cases or male deaths were not provided, but percentages of male cases or male deaths were provided instead. To include these sources and avoid biases that might be introduced by their exclusion, we calculated counts of male confirmed cases and male deaths from the reported percentages with rounding to the nearest integer. We acknowledge that this approach assumes that the reported percentages are reflective of the true percentages. For some sources, data included confirmed cases and deaths of unknown sex. For these sources, the reported totals were used where the proportion of unknown sex was small (see Supplementary Data 1). This approach was preferred to excluding cases of unknown sex in order to avoid bias. The estimates represent the proportion of known male infections and ORs for mortality associated with known male sex, and will differ slightly from what the true values would be if the sex had been reported for all cases. Data were available at the level of country or regional summary data representing distinct individuals for each report, but not at the level of covariates for all individuals within a study. Consequently, covariates such as lifestyle, comorbidities, testing method and case type (hospital vs. community) could not be controlled for. Data represent unique individuals at a single time snapshot. All data are provided in Supplementary Data 17,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,79,80,81,82,83,84,85,86,87,88,89,90,91,92,93,94,95,96,97,98 and https://doi.org/10.25375/uct.12952151.

Data analysis

Meta-analysis was performed to estimate an overall proportion of male infected cases with 95% CI. This meta-analysis was two-sided and tested the null hypothesis that the proportion of male infected cases was 0.5. For this analysis, the classic inverse variance method for estimation of single proportions and standard errors was used, which uses logit-transformed proportions. The inverse variance method accounts for differing sample sizes of individual studies by weighting studies by the variance of their estimates, such that small studies with large variance have less weighting, and large studies with small variance have more weighting. A sensitivity analysis was also performed, in which the generalised linear mixed model (GLMM) method was used to estimate an overall proportion. This yielded identical results, indicating the differing assumptions of these different methods were inconsequential for these data.

Meta-analyses were also performed to estimate ORs with 95% CI associated with male sex for ITU admission and death, based on pooled average effect measures that were weighted according to the size and precision of each report. Random effects models were estimated and are reported rather than fixed effects models, since these do not assume uniformity across reports and account for variance between reports. The Mantel-Haenszel and DerSimonian-Laird methods were used to calculate the fixed effects and random effects estimates, respectively. Similar to the inverse variance weighting method, individual studies are weighted according to size and variance, and estimates were almost identical when the inverse variance weighting method was used. Both of these meta-analyses for the ITU admission and death outcomes were two-sided and tested the null hypothesis that the estimated OR was 1. Reports that did not contain the data required to calculate ORs were automatically excluded from meta-analyses.

Meta-analyses were performed using R version 3.6.1 and the “meta” package version 4.11-0158. All R code and input data are available at https://github.com/claire-deakin/covid19-sex-bias and at https://doi.org/10.25375/uct.12952151.

Reporting Summary

Further information on research design is available in the Nature Research Reporting Summary linked to this article.