Abstract
To successfully mitigate the extraordinary devastation caused by the Coronavirus disease 2019 (COVID-19) pandemic, it is crucial to identify important risk factors for this disease. One such neglected health determinant is the sex of the patient. This is an essential clinical characteristic, as it can factor into a patient’s clinical management and preventative measures. Some clinical studies have shown disparities in the proportion between males and females that have more severe clinical outcomes or, subsequently, die from this disease. However, this association has not been unequivocally established. Thus, the purpose of this investigation was to examine the association between male sex and COVID-19 severity. We systematically reviewed the literature, identified non-randomized studies that matched predetermined selection criteria, and performed a meta-analysis to evaluate the proportion of males among four disease severity categories. Appropriate assessment strategies were implemented to assess and minimize potential biases. The results of this meta-analysis indicated that males constituted a significantly higher proportion of those who had adverse clinical outcomes and died from COVID-19. As the coronavirus spread from the East to the West, male sex remained a consistent risk factor. Our results support the establishment of the male sex as an important risk factor for this disease. Early identification and appropriate medical care for males with lab-confirmed COVID-19 may substantially change the course of clinical prognosis, resulting in greater numbers of lives saved.
1 Introduction
Males and females have distinct biological, immunological, and endocrine differences that result in different disease processes and outcomes. Sex-specific differential gene expression and molecular-level variation have been reported to influence blood pressure, cardiovascular health, and kidney function (Convertino, 1998; Reckelhoff, 2001; Fischer et al., 2002; Kang and Miller, 2002; Sandberg and Ji, 2003; Hilliard et al., 2013). Females, in general, have a heightened capability to activate a greater and more robust immune response, offering protection against many infectious disease processes, but may predispose them to an array of autoimmune diseases (Fairweather and Rose, 2004; Pennell et al., 2012; Rubtsova et al., 2015; Klein and Flanagan, 2016; vom Steeg and Klein, 2016; Jaillon et al., 2019). Males and females also express immunological dimorphisms. Females have two X chromosomes in comparison to the XY in males. The random transcriptional inactivation of X chromosomes in females may also help offset certain mutation-related dysregulation of the immune system (Taneja, 2018). Differences in endocrine system regulation in females compared to males significantly affect disease processes including respiratory, cardiovascular, and renal disease (Hilliard et al., 2013; Blenck et al., 2016; Channappanavar et al., 2017; Palmisano et al., 2018; Vermillion et al., 2018; Wensveen et al., 2019). As nations across the world navigate their way through the Coronavirus disease 2019 (COVID-19) pandemic, clinical, research, and public health experts have observed that this disease does not affect all individuals alike.
Since the beginning of 2020, the world’s healthcare professionals have tirelessly attempted to mitigate the impact of the COVID-19 pandemic. With over 3.5 million confirmed cases and 243,000 deaths worldwide as of May 5th, 2020, a post-COVID-19 pandemic era is not within the near foreseeable future (WHO, 2020). The United States, one of the epicenters for the disease, has documented over 1.1 million confirmed cases and 68,000 deaths related to COVID-19 (CDC, 2020; WHO, 2020). Many recent studies have highlighted certain risk factors that cause specific populations to be disproportionately susceptible to the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection. Currently known risk factors for severe clinical outcomes of COVID-19 include: advanced age (65 years and older), chronic lung diseases, immunocompromised status, and other comorbidities such as hypertension, diabetes, and/or cardiovascular disease (Emami et al., 2020; Grasselli et al., 2020; Remuzzi and Remuzzi, 2020; Shahid et al., 2020; Zheng et al., 2020).
Observations in COVID-19 patient data involving clinical characteristics highlight specific disparities in males and females. A recent case-series study looking at COVID-19 and SARS patients showed that while males and females had the similar disease prevalence, males with COVID-19 were at higher risk for worse clinical outcomes and death (Jin et al., 2020). In this study, as the patient age and the documented comorbidities (i.e., cardiovascular diseases, diabetes, chronic lung diseases, or hypertension) increased, the risk of severity and mortality in both COVID-19 and SARS patients increased. However, the mortality rate in males was 2.4 times that of their age-matched female counterparts (70.3% and 29.7%, respectively).
Furthermore, a nationwide COVID-19 surveillance study conducted in Italy indicated that male mortality rates related to COVID-19 were disproportionately higher than that of female patients with a ratio as much as 4 to 1 (Remuzzi and Remuzzi, 2020). Other systematic reviews performed to characterize clinical features or risk factors for COVID-19, have also identified the sex-specific disparities in disease severity and mortality (Li et al., 2020a; Zheng et al., 2020). However, the clinical importance of male sex as a risk factor for COVID-19 has mainly been overlooked or explained as a potential confounder to other environmental factors such as smoking or tobacco product usage (Cai, 2020). While various studies have made observations of the sex-specific disparities of COVID-19, this specific relationship has not been adequately established. The sex-specific disease severity is an important clinical consideration as it affects all patient populations. Recognition of male sex as a risk factor for COVID-19 will impact both preventative measures and clinical patient management protocols.
The goal of this systematic review and meta-analysis is to identify whether males are more susceptible to COVID-19, severe forms of the disease, or mortality related to COVID-19. To address this question, we systematically reviewed the literature using the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) guidelines. We performed a meta-analysis of the selected study populations comparing male and female COVID-19 patients. This review incorporated three online databases and research studies published between December 15th, 2019, and April 16th, 2020. We characterized the influence of sex as a risk factor for COVID-19 measuring the following clinical outcomes: all lab-confirmed cases, severe cases, critically ill cases, and mortality.
2. Methods
2.1 Literature Search and Research Study Selection
We performed a comprehensive systematic research literature search of three online databases, PubMed (LitCOVID), Embase (OVID), and Web of Science (WoS), from December 15th, 2019, to April 16th, 2020. We identified all research articles related to COVID-19 that contained any sex-specific patient or clinical characterizations. The search terms and keywords used to identify research studies for the meta-analysis were: COVID-19, male, female, men, women, sex, and gender. We reviewed references of review, perspectives, systematic reviews, and meta-analysis articles of the include articles to ensure comprehensiveness of our search. All our search results were evaluated using the PRISMA statement. We reviewed the abstracts and tables of each of the articles to identify the presence of sex-specific (male and female) COVID-19 case numbers.
2.2 Eligibility criteria
The inclusion criteria for research article selection was as stated below. Study population: patients with lab-confirmed COVID-19 diagnosis. Study design: case series or cross-sectional study that did not exclude any lab-confirmed COVID-19 patients. Outcomes measure: at least one outcome reported with male to female ratio among lab-confirmed clinical cases, severe cases, critical cases, and mortality. Research study: only peer-reviewed research publications were included. Commentary articles, perspectives, review articles, and surveillance reports were excluded. The following case definitions were used in this study. All cases were lab-confirmed COVID-19 patients. Severe cases were defined as having at least one of the following clinical findings: (a) breathing rate ≥30/min, (b) oxygen saturation (SpO2) ≤ 93% at rest, or (c) ratio of the partial pressure of arterial oxygen (PaO2) to the fraction of inspired oxygen (FiO2) ≤300 mmHg. The severe case definition followed the American Thoracic Society guidelines for community-acquired pneumonia (Metlay et al., 2019). Critical cases were defined as: (a) received mechanical ventilation; (b) clinically diagnosed with shock symptoms, (c) received care in the intensive care unit (ICU) or (d) transfer to a higher level of medical care.
2.3 Data Extraction and Quality Assessment
All articles identified through the keyword search from the online databases were organized into an Excel® spreadsheet. Following the removal of duplicates, articles were subjected to evaluation, and five investigators did data extraction. Research studies were screened using the abstract and any tabulated clinical characteristics of COVID-19 patients. Directly after that, research articles were again screened to identify any discrepancies by an independent investigator. The screened articles were assessed against the study selection criteria by two independent investigators, and any differences in selected articles were revisited, and a definitive determination was made.
2.4 Selected Study Bias Risk Assessment
A bias risk assessment was conducted on studies included in the meta-analysis utilizing the methodological index for non-randomized studies (Minors) criteria at the study level (Slim et al., 2003). Each of the selected articles was scored with 0 (not reported), 1 (reported but inadequate), or 2 (reported and adequate). The highest score possible was 16 for non-comparative studies according to Minors guidelines.
2.5 Statistical Analysis of Selected Data Sets
Statistical analysis was conducted using R (R_Core_Team, 2019) with the meta-analysis packages meta (Balduzzi et al., 2019) and dmetar (Harrer et al., 2019). The principal summary measures of the meta-analysis were proportions of males in four different COVID-19 categories. The four groups were: (a) all confirmed COVID-19 cases, (b) severe cases of COVID-19 as defined in Section 2.3, (c) critically ill cases of COVID-19 as defined in Section 2.3, and (d) deaths associated with COVID-19. Agresti-Coull confidence intervals were used for individual studies. Studies were combined using the inverse variance method on the raw proportions with the DerSimonian-Laird estimator for the between-groups variance estimator (τ2) and the Jackson method for combined confidence intervals. Transformations of raw proportions were calculated for the combined estimates (log, logit, arcsin, and Freeman-Tukey double arcsin), but the results were so similar they are not shown. The proportion of variation in treatment effects was estimated with I2. To assess bias across studies, funnel plots were constructed for each of the four different categories, and Egger’s bias test conducted. In order to determine if there were region-specific differences among populations in Asian and Western countries, we sub-divided the COVID-19 critically ill patient populations into these two regions and analyzed them.
2.6 Clinical Outcomes Median Age Calculation
To combine the ages, in 20 of the articles, the median age of patients was given, along with sometimes interquartile range, sometimes min and max. In 10 of the articles, mean and standard deviation (SD) were presented. In one article (Easom et al., 2020), the mean age was given without SD. We used linear regression on the other 10 (mean, SD) pairs to estimate the SD to be 14.5 years. To combine the ages, we chose to convert means to medians because there would be fewer unknown statistics to estimate, and typical disease distributions are skewed. To convert, we fit a negative binomial distribution to the mean and SD using the method of moments. With the complete list of medians, we used R’s metamedian (McGrath et al., 2020) package to obtain summarized confidence intervals for each of the four categories.
3. Results
3.1 Research Study Selection and Quality Assessment
We identified 786 research articles that matched our search terms. After the duplicated were removed, 414 unique research articles were screened. Following the screening process, 353 articles with incomplete data were excluded. We then identified 61 research articles with sex-specific case numbers and reviewed full-length articles to assess their eligibility for our study according to the selection criteria. Thirty articles did not fit the selection criteria and were excluded from the meta-analysis. Reasons for exclusion were: not a primary research study (a surveillance report or perspective), did not include consecutive patients or did not meet with the case or severity definitions. The 31 research articles eligible for this meta-analysis were used for qualitative synthesis and quantitative analysis (Figure 1). The 31 eligible articles were subjected to a bias assessment using the Minors criteria at the study level (Slim et al., 2003). All 31 selected articles scored between 12 and 14 points, with 16 being the highest for non-randomized controlled studies (Table 1). The relatively high scores indicated that we were likely not introducing any significant systemic biases.
3.2 Study Population Demographics
Within our selected studies, 7556 lab-confirmed COVID-19 cases were identified. Of these 31 studies, 24 were from various cities in China and included a sample of 5629 lab-confirmed cases. Two studies were from South Korea and Singapore, which included a sample of 72 lab-confirmed cases. The other five studies were from Europe and North America, having a sample of 1855 lab-confirmed cases (Figure 2 and Table 2). Most of the early studies came from China with study periods from December 11th, 2019, to February 24th, 2020. Most of the later studies came from other countries with study periods from January 23rd to April 5th, 2020 (Figure 3). These patterns reflect the movement of epicenters for COVID-19 from the East to the West.
3.3 Meta-Analysis and Bias Assessment
The principal quantitative results are contained in the forest plots shown on the left side of Figures 4 and 5. The individual confidence intervals are shown, by study, with the combined proportion for each group and confidence interval at the bottom. A random-effects model was used for the combined proportion to check for heterogeneity (τ2= between-group variation and I2= proportion of total variation in the estimates of treatment effects due to heterogeneity). The heterogeneity statistics (τ2 and I2) are shown at the bottom left of the forest plots.
3.4 Sex-Specific COVID-19 Case Distribution
A total of 23 studies with 5408 lab-confirmed COVID-19 cases were analyzed (Table 3). Our results from the randomized effects model meta-analysis showed that in the sex-distribution of all COVID-19 cases, males accounted for 53% (95% CI [0.51, 0.55]) (Figure 4a). Female patients made up 47% of all COVID-19 cases. There is moderate heterogeneity between the set of overall population proportions (I2 = 64%, τ = 0.05). A funnel plot was drawn to assess the publication bias (Figure 4b). The publication bias test results: Egger’s test (p = 0.88) indicated that there was no publication bias.
3.5 Sex-Specific COVID-19 Severe Case Distribution
A total of 8 studies with 985 severe COVID-19 cases were analyzed (Table 4). Our results from the randomized effects model meta-analysis showed that in the sex-distribution of all COVID-19 severe cases, males accounted for 56% (95% CI [0.53, 0.59]) (Figure 5a). Female patients made up 44% of all COVID-19 severe cases. There is no heterogeneity for the severe population proportions (I2 = 0%, τ = 0.0). A funnel plot was drawn to assess the publication bias (Figure 5b). The publication bias test results: Egger’s test (p-value = 0.40) indicated that there was no publication bias.
3.6 Sex-Specific COVID-19 Critically Ill Case Distribution
A total of 9 studies with a total of 2025 critical COVID-19 cases were analyzed (Table 5). Our results from the randomized effects model meta-analysis showed that in the sex-distribution of all COVID-19 critically ill cases, males accounted for 71% (95% CI [0.63, 0.79]) (Figure 5c). Female patients made up 29% of all COVID-19 critical cases. There was strong heterogeneity between the critical population proportions (I2 = 83%, τ = 0.10). A funnel plot was drawn to assess the publication bias (Figure 5d). The publication bias test results: Egger’s test (p-value = 0.02) indicated that there could be some publication bias introduced by the Grasselli et al. (2020) study.
3.7 Sex-Specific COVID-19 Mortality Distribution
A total of 6 studies with a total of 399 mortalities related to COVID-19 cases were analyzed (Table 6). Our results from the randomized effects model meta-analysis showed that in the sexdistribution of all COVID-19 mortalities, males accounted for 69% (95% CI [0.63, 0.75]) (Figure 5e). Female patients made up 31% of all COVID-19 mortalities. The heterogeneity for the mortality population proportions is low (I2 = 34%, τ = 0.04). A funnel plot was drawn to assess the publication bias (Figure 5f). The publication bias test results: Egger’s test (p = 0.26) indicated that there was no observable publication bias.
3.8 Sex-Specific COVID-19 Distribution in Asia and the West
Sex-specific differences in clinical outcomes of COVID-19 cases in China were thought to be related to cultural and social differences in males and females (Cai, 2020). We investigated if our study results hold in different regions of the world. COVID-19 critically ill patient data sets were divided into two groups: Asia and West, and subgroup analyses were performed.
A total of 5 studies from Asia, with a total of 238 critical COVID-19 cases were analyzed. Our results from the randomized effects model meta-analysis showed that in the sex-distribution of COVID-19 critically ill cases from Asia, males accounted for 71% (95% CI [0.61, 0.81]) (Figure 6a). Female patients made up 29% of all COVID-19 critical cases in Asia. There was moderate heterogeneity between the critical population proportions (I2 = 64%, τ = 0.0082). A funnel plot was drawn to assess the publication bias in studies from Asia (Figure 6b). The publication bias test results: Egger’s test (p-value = 0.26) indicated that there was no observable publication bias.
A total of 4 studies from Western regions with a total of 1787 critical COVID-19 cases were analyzed. Our results from the randomized effects model meta-analysis showed that in the sex-distribution of COVID-19 critically ill cases from the West, males accounted for 70% (95% CI [0.59, 0.82]) (Figure 6c). Female patients made up 30% of all COVID-19 critical cases in the West. There was strong heterogeneity between the critical population proportions (I2 = 86%, τ = 0.0103). A funnel plot was drawn to assess the publication bias in studies from the West (Figure 6b). The publication bias test results: Egger’s test (p-value = 0.04) indicated that there could be some publication bias introduced by the Grasselli et al. (2020) study, as indicated previously. This comparative subgroup analysis of Asia and the West indicated that there was no geography-specific difference in the proportion of critically ill COVID-19 male patients. However, indicated by the moderate to strong heterogeneity observed, there are likely variations in male proportion between different studies and regions.
3.9 Disease severity stratification and age distribution
When extracting male and female proportions for each of the four COVID-19 disease severity categories, we obtained the age distributions of the cases stated as a mean ± SD or median and interquartile range (IQR). Using a skewed distribution assumption, the ages were aggregated as medians with 95% confidence intervals. The median age for all COVID-19 cases was 50, severe cases was 61, critically ill cases was 63, and mortality was 70 (Figure 7). A Kruskall-Wallis ranked-sum test conducted on the medians showed that age was significantly different between the COVID-19 disease severity groups (chi-squared = 24.07, df = 3, p-value = < 0.0001). Our data confirm that advanced age is a risk factor for more severe clinical outcomes and mortality related to COVID-19.
4. Discussion
In our systematic review and meta-analysis, we set forth to address the question of whether male sex is a risk factor for COVID-19 susceptibility, severe forms of the disease, or mortality related to COVID-19. Systematically reviewing all literature from December 15th, 2019, to April 16th, 2020, we selected 31 research studies that met our selection criteria and performed a meta-analysis on COVID-19 clinical outcomes. Our quality assessment measures indicated low heterogeneity in terms of a single-arm meta-analysis, and the sensitivity analysis showed that there was minimal publication bias. As of the time of completing this manuscript, there were no randomized controlled trials with COVID-19 patients that could address this particular question. The use of non-randomized studies for the meta-analysis is a limitation of this study. However, Abraham et al. (2010) suggested that, in the absence of randomized, controlled trials, that a well-designed meta-analysis is used non-randomized controlled trials can also present a high level of evidence (Abraham et al., 2010).
The four clinical outcome categories (overall, severe, critical, mortality) exhibited different levels of heterogeneity in our random-effects models. The explanation for these differences is most likely the region of the studies done within each category. The 23 overall studies exhibited 64% heterogeneity with one from Singapore and one from Great Britain. The eight severe studies exhibited 0% heterogeneity, all being from China. The nine critical studies exhibited 83% heterogeneity, with five from China, three from Europe, and one from the United States. The six mortality studies exhibited 35% heterogeneity, with five from China and one from Korea. The use of a randomized effects model for our meta-analysis takes into account these heterogeneities observed between different studies and regions. Based on the random-effects models shown, there appears to be a difference in the proportions of males with COVID-19 between at least some of the studies or regions. Due to the study designs, their sampling methods, and limited regions included in this study, it is neither possible nor wise to be more specific. This is a potential avenue for further research.
Our meta-analysis showed that while males accounted for 53% of all COVID-19 cases, males accounted for an increasing proportion of severe cases (56%), critically ill cases (71%), and mortalities (69%) compared to their counterpart. While similar male to female disproportions was observed among a few other studies looking at clinical characteristics of COVID-19, our study provides a comprehensive synthesis of data available across different world regions. This study helps establish male sex as a risk factor for COVID-19 clinical outcomes and shows that it is consistent in Asia and Western regions.
This study results do not come as a surprise. Several studies conducted on the two previous coronavirus epidemics, SARS CoV-1in 2002-2003 and MERS in 2012-2013, showed similar patterns with a male predominance toward greater severity and mortality risks. Studies on mortality rates during the MERS-CoV epidemic showed the male sex to be a risk factor (Matsuyama et al., 2016; Nam et al., 2017; Park et al., 2018). Epidemiological studies with SARS-CoV-1 showed similar patterns (Karlberg et al., 2004). To further support previous epidemiological observations, in controlled mouse model experiments, SARS-CoV-1 has displayed infectious dose-dependent higher mortality rates in male mice compared to female mice (Channappanavar et al., 2017). The mounting amount of evidence showing differences among males and female clinical outcomes to coronavirus infections highlights the importance of patient sex in determining the COVID-19 prognosis.
From a clinical standpoint, this information is very pertinent to the practice of patient care. As COVID-19 clinical outcomes are strongly associated with male sex, this can help guide preventative and treatment strategies. Male patients will likely warrant more aggressive inpatient care measures, and especially those that have other COVID-19 risk factors such as advanced age or underlying comorbidities. Susceptible males with other known risk factors may need to take extra precautions to help prevent SARS-CoV-2 infection. Infected males can be encouraged to obtain medical care at an earlier stage of the disease. In cases that require hospitalization, physicians can be more aggressive in their medical management of male patients.
In addition to preventative and COVID-19 treatment measures, this presents a unique clinical opportunity to address male and female differences at the molecular level, immunological response, and endocrine function (Sandberg and Ji, 2003; vom Steeg and Klein, 2016; Taneja, 2018). For example, SARS-CoV-2 binds to the Angiotensin-converting enzyme 2 (ACE2) receptors and use it as a mechanism for host cell entry (Hoffmann et al., 2020). Males have been shown to express more ACE2 receptors within the renin-angiotensin-aldosterone system (RAAS) (Komukai et al., 2010). This is likely to play an essential role in the severity of this disease observed in males (Komukai et al., 2010). Differences in male and female immunological responses will also be a clinically significant factor that can be appropriately modulated to better serve COVID-19 patients (Pennell et al., 2012; Jaillon et al., 2019). Besides sex-specific differences in immunological responses, hormonal regulation and the role of estrogen and testosterone in priming the ACE2 receptor sensitivity could hold the key to better explain the higher COVID-19 severity and mortality rates observed in males (Mishra et al., 2016; Bukowska et al., 2017). In an age of personalized medicine, if molecular level of differences in the disease processes of SARS-CoV-2 infection can be characterized, clinicians can use targeted therapy using to promote health equality and help save more patient lives.
Data Availability
All datasets presented in this study are included in the article or supplementary files.
Funding
The work of R.G. and T.G. was supported in part by the Discovery Institute and the Peter & Carla Roth Family.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Author Contributions
TG led the systematic review, helped prepare the tables and figure, and aided in writing and editing. BMP aided in the analytical evaluation of curated articles, writing, and editing. JA, AB, and DR research students analyzed data and joined in discussions. JW performed statistical analyses; drafted statistical sections. RSG conceptualized the problem and aided in writing and editing.
Supplementary Material
Supplementary Materials are provided as a separate file.
Acknowledgments
We acknowledge Dr. Siva Somasundaram, Ms. Avishka Jayasekera, and Dr. John P. Walsh of the University of Southern California for their collaborative support. We also thank Dr. Jeffrey S. Wang, Infection Disease Specialist at Kaiser Permanente, Anaheim, California, for his clinical insights.