Article Text

Original research
Does comorbidity burden explain the higher COVID-19 mortality risk among men? A retrospective cross-sectional analysis of a well-defined cohort of patients in Bronx, New York
  1. Aastha Vasa1,
  2. Maya Kini2,
  3. Joel Neugarten1,
  4. Eran Bellin3,
  5. Ladan Golestaneh1
  1. 1Department of Medicine, Montefiore Medical Center/Albert Einstein College of Medicine, Bronx, New York, USA
  2. 2Rye Country Day School, Rye, New York, USA
  3. 3Epidemiology and Population Health, Montefiore Medical Center, Yonkers, New York, USA
  1. Correspondence to Aastha Vasa; aastha.vasa{at}einsteinmed.edu

Abstract

Objectives Men have a higher mortality rate and more severe COVID-19 infection than women. The mechanism for this is unclear. We hypothesise that innate sex differences, rather than comorbidity burden, drive higher male mortality.

Design Retrospective cohort.

Setting Montefiore Health System (MHS) in Bronx, New York, USA.

Participants A cohort population of 364 992 patients at MHS between 1 January 2018 and 1 January 2020 was defined, from which individuals hospitalised during the pre-COVID period (1 January 2020–15 February 2020) (n=5856) and individuals hospitalised during the COVID-19 surge (1 March 2020–15 April 2020) (n=4793) were examined for outcomes. A subcohort with confirmed COVID-19+ hospitalisation was also examined (n=1742).

Primary and secondary outcome measures Hospitalisation and in-hospital mortality.

Results Men were older, had more comorbidities, lower body mass index and were more likely to smoke. Unadjusted logistic regression showed a higher odds of death in hospitalised men than women during both the pre-COVID-19 and COVID-19 periods (pre-COVID-19, OR: 1.66 vs COVID-19 OR: 1.98). After adjustment for relevant clinical and demographic factors, the higher risk of male death attenuated towards the null in the pre-COVID-19 period (OR 1.36, 95% CI 1.05 to 1.76) but remained significantly higher in the COVID-19 period (OR 2.02; 95% CI 1.73 to 2.34).

In the subcohort of COVID-19+ hospitalised patients, men had 1.37 higher odds of in-hospital death (95% CI 1.09 to 1.72), which was not altered by adjustment for comorbidity (OR remained at 1.38 (95% CI 1.08 to 1.76)) but was attenuated with addition of initial pulse oximetry on presentation (OR 1.26, 95% CI 0.99 to 1.62).

Conclusions Higher male mortality risk during the COVID-19 period despite adjustment for comorbidity supports the role of innate physiological susceptibility to COVID-19 death. Attenuation of higher male risk towards the null after adjustment for severity of lung disease in hospitalised COVID-19+ patients further supports the role of higher severity of COVID-19 pneumonia in men.

  • public health
  • public health
  • epidemiology
  • epidemiology
  • respiratory infections
  • COVID-19

Data availability statement

All data relevant to the study are included in the article or uploaded as online supplemental information. All data relevant to this study are provided upon request to the senior author.

http://creativecommons.org/licenses/by-nc/4.0/

This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/.

Statistics from Altmetric.com

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

STRENGTHS AND LIMITATIONS OF THIS STUDY

  • Using Looking Glass, we built a database of patients (n=364 992) who received outpatient care between 1 January 2018 and 1 January 2020 and was defined as the denominator patient population.

  • Those patient with hospitalisations and in-hospital deaths during the pre-COVID-19 and COVID-19 periods were examined as a proportion of the denominator population allowing for a robust assessment of risk.

  • To elucidate the role of baseline comorbidity in higher risk of COVID-19 mortality in men, we used mixed-effects logistic regression to examine the association between sex and hospitalisation/death (outcomes) with and without adjustment for patient comorbidity.

  • The analysis did not include outpatient deaths which limited the scalability of the findings and introduced sampling bias.

Introduction

Numerous studies have shown that men have a higher risk of COVID-19-related death as compared with women.1–14 Men also have a more severe disease course, with higher odds of in-hospital adverse outcomes such as ICU admission, need for ventilatory support, and higher odds of acute kidney and cardiac injury.3 Mediators of the risk difference between the sexes have been attributed to differences in comorbidities, risk behaviours (such as smoking and not seeking preventive healthcare), exposure rate (more men work in occupations that force them to leave their homes and have contact with other people) and biological differences in susceptibility.2 8 For example, preliminary data from Gebhard et al describe an association between comorbidities, risk behaviours (ie, smoking) and COVID-19 disease severity.11 Another study from a cohort of patients in Bronx, New York, showed that obesity, male sex and older age were associated with greater in-hospital mortality.15 Other researchers argue that inherent biological differences in ACE expression, androgen-mediated regulation of transmembrane serine protease 2 (TMPRSS2), as well as sex differences in immune responses, are the actual mediators.5 7–9 16 17

To address the question of what is mediating the higher male risk for COVID-19-related mortality, cohort studies should adjust for comorbidity and viral exposure to test the strength of association between male sex and COVID-19 outcome, with the understanding that a demonstration of the attenuation of this strength supports the role of these explanatory variables in imparting higher male risk. In this study, we first establish the disproportionate higher risk of in-hospital male death in one large health system in the Bronx during a period impacted by the COVID-19 surge of 2020 as compared with a pre-COVID time period. Because of the comprehensive nature of our clinical database and because this was still a period where vaccinations were not available, we are uniquely able to evaluate the role of comorbidity and smoking behaviour as an explanation for the higher risk of in-hospital COVID-19-related death in men. We then examine the role of severity of SARS-CoV-2-related pneumonia as a plausible explanation for any higher male risk found in a subcohort of patients hospitalised with confirmed COVID-19 illness.

Materials and methods

To test the association of sex with COVID-19-related outcomes, we: (1) test the relative risk of male hospitalisation and in-hospital mortality, as compared with female, during the initial COVID-19 surge of New York (1 March 2020–15 April 2020) and compare this risk to a pre-COVID-19 time period (1 January 2020–15 February 2020)18; (2) then compare the odds of male versus female death in a subset of patients hospitalised with a confirmed positive COVID-19 test during the COVID-19 surge time period and (3) provide a descriptive analysis of the city/state and national level impact of COVID-19 death by sex using Vital Statistics/Centers for Disease Control and Prevention (CDC) data.19 20

Study populations

  1. To minimise sampling bias, a baseline population of 364 992 Bronx addressed adult patients (>18 years) receiving care at Bronx Montefiore Health System (BMHS) between 1 January 2018 and 1 January 2020 was defined. From this population, we extracted variables and outcomes data on individuals hospitalised during the pre-COVID-19 period (1 January 2020–15 February 2020) and during the COVID-19 surge (1 March 2020–15 April 2020). Clinical Looking Glass (CLG),21 Montefiore’s clinical/claims database, was used to collect individual demographic and co-morbidity data: including sex, age (categorised by deciles), race/ethnicity (self-identified and categorised into racial/ethnic categories), diabetes mellitus status (based on International Classification of Diseases (ICD) coding and/or evidence for Hemoglobin A1c (HgbA1c) level >6.5 within 2 years prior to entry into cohort), hypertension (based on ICD coding and/or evidence for systolic blood pressure> 150 mm Hg within 2 years prior to entry into cohort), asthma (based on ICD coding), Charlson score (comorbidity score calculated by CLG based on ICD-10 coding for diagnoses required for calculation of the score within 6 months prior to entry into cohort22 23), any history of smoking and body mass index (BMI, based on latest BMI recorded at any type of clinical visit at Montefiore). Individual addresses were linked to the American Community Survey24 to define census tract level sociodemographic variables (% of households living under the poverty line, completed high school, had health insurance, used public transportation to commute to work and had access to internet). We also defined mean household size by census tract.

  2. The subset of patients, drawn from the baseline population (see above), who were hospitalised during the COVID-19 period (n=1742) with evidence for COVID-19 testing and a positivity of COVID-19 test was analysed separately.

  3. Finally, to test the generalisability of the findings of our cohort to a larger city, state and national population, the New York Department of Health’s Vital Statistics registry was accessed for total population and death estimates to build proportions of all-cause male/female death rates in four time periods (2015–2018).19 The time periods include 4 years of data prior to the COVID-19 surge, and proportion of city-wide death during COVID-19, were similarly reported. Using 2020 Census data to derive baseline breakdown of sex in the denominator population, proportions were calculated to determine the rates of COVID-19 cases, admissions and hospitalisations for men versus women.25 Furthermore, CDC data were accessed to evaluate the proportion of COVID-19-related as compared with COVID-19-unrelated deaths by sex in New York State and across the country as of 29 March 2022.

Outcome variable

Proportion of patients who had an in-hospital death in either time period.

Exposure variable

Patient sex.

Statistical analysis

Stata V.15.0 was used for analyses. Bivariate analysis was done to compare distribution between male and female of demographic (race and age), clinical risk (history of diabetes, hypertension, asthma, smoking status, Charlson Comorbidity Score and BMI) and community-level variables using t-test for normally distributed variables and Wilcoxon rank-sum test for non-normally distributed variables. In the subset of hospitalised COVID+patients, initial respiratory rate, white ell count, systolic blood pressure, respiratory rate and pulseoximetry at initial hospital encounter were also included in sex-category-based bivariate analyses.

Mixed effects logit modelling (random effects presented by individual patients in duplicated baseline populations, each representing a different time period) tested the probability of male versus female hospitalisation rate and death in and between the pre-COVID-19 and COVID-19 time periods, while adjusting for all demographic and clinical risk variables that were significantly different by sex in bivariate analyses.26 Interaction between sex and COVID-19 versus pre-COVID-19 time period with respect to hospitalisation and mortality outcome was tested as an additive ‘lincom’ postestimation test.27 We then stratified logistic regression models by time period while adjusting for similar confounding variables.

For the hospitalised COVID+ subpopulation, logistic regression modelling, using stepwise forward addition techniques (at a significance for inclusion of p<0.05), tested the association of sex and mortality while adjusting for demographic and comorbid confounding variables, and those clinical variables on presentation to the hospital.

Simple proportions were used to test the sex-based all-cause mortality differences between sexes of data from New York City during the 2015–2018 and the city-wide mortality during the COVID-19 surge.

Sensitivity analysis

The BMI variable had ~18% missing observations and was excluded from both regression models. In sensitivity analyses, we included BMI into the multivariate adjusted models to examine any difference in the coefficient of interest with and without its inclusion, and any change in the accuracy of the drawn inferences.

Because of the inherent difficulty with interpreting differences in community-level attributes between men and woman (units in census data defined as households and not individuals), they were left out of the final regression models. In sensitivity analyses, we included community level attributes (in bulk, see table 1) and examined any difference in coefficient of interest with and without their inclusion.

Table 1

Bivariate comparison of demographic, clinical and community-level variables in men versus women

Patient and public involvement

Patients and the public were not involved in the design, conduct, outcome measures or recruitment of this study.

Results

The mean age of the baseline population was 49.9 years old, 7.8% were white, 28.5% were black and 37.8% were Hispanic. 36.5% of the baseline population was male, with 16.4% diagnosed with diabetes, 44.9% with hypertension and 11.8% with asthma (table 1).

Men were slightly older (mean age of 50.0 years for men vs 49.9 years for women; p=0.001), had more comorbidities (20.1% men had a >3 Charlson score vs 19.2% women; p<0.001), more had diabetes (17.8% vs 16.5%; p<0.001) and hypertension (47.2% vs 43.6%; p<0.001) and had a lower BMI than women (median of 28.0 for men vs median of 29.2 for women; p<0.001). More men smoked than women (36.4% vs 24.0%; p<0.001), but a lower proportion had asthma (8.4% vs 13.8%; p<0.001). There was no difference between men and women with respect to census tract average household size or percentage of households that completed high school. More men had an active internet subscription (mean of 73.2% for men vs 72.8% for women; p<0.001), fewer men came from households that lived under the poverty line (mean of 23.5% for men vs 24.2% for women; p<0.001) and fewer men used public transportation to commute to work (mean of 59.1% men vs 59.7% women; p<0.001) (table 1).

During the pre-COVID-19 period, 5856/364 992 individuals were hospitalised (1.6%), while during the COVID-19 period 4793/364 992 (1.3%) were hospitalised. In unadjusted and adjusted models, men had a higher odds of hospitalisation in both time periods (table 2). The interaction between male sex and time period with respect to odds of hospitalisation was statistically significant (p<0.001) suggesting a higher risk of male hospitalisations during COVID-19 as compared with pre-COVID-19 time period. Furthermore, the odds of pre-COVID-19 male hospitalisation attenuated towards the null after adjustment for relevant comorbidities including Charlson Comorbidity Score, but remained robust without any attenuation during the COVID-19 period (table 2).

Table 2

Odds of male hospitalisation and death during the pre-COVID-19 and COVID-19 time period

Two hundred and fifty-four individuals from the pre-COVID-19 time period and 748 individuals from the COVID-19 period had an in-hospital death event. During the pre-COVID-19 time period, men had greater in-hospital deaths, with a 0.03% difference in proportion who died between men and women. During the COVID-19 time period, the difference between male and female proportions of in-hospital deaths was 0.15%. Stratified analyses showed men were had a positive and significant odds of death during both pre-COVID-19 and COVID-19 time periods (OR pre-COVID 1.66 (95% CI 1.30 to 2.12) vs OR COVID-19 1.98 (95% CI 1.71 to 2.28)) as compared with women. Addition of confounding demographic and clinical variables attenuated the odds of male death in the pre-COVID-19 period but not in the COVID-19 period (table 2).

Sensitivity analysis showed that addition of BMI to the logit model altered the adjusted OR of male death in the pre-COVID-19 (from OR 1.37 (95% CI 1.05 to 1.77) to OR 1.29 (95% CI 0.99 to 1.67); LR test of model with BMI nested vs without: p<0.001; n=3 19 753) and in the COVID-19 (from OR 1.99 (95% CI 1.71 to 2.31) to OR 1.99 (95% CI 1.71 to 2.31); LR test p=0.94). The addition of BMI to the logit model testing for interaction showed no change in the p value for interaction between sex and time period with respect to mortality (p for interaction 0.24 with vs 0.25 without BMI in the model; LR test p=0.6) Addition of community-level variables including percentage of households with high school education, living under poverty line, using public transport to commute to work, access to internet and mean household size, did not alter the OR of male death in the pre-COVID-19 period (OR 1.35 (95% CI 1.04 to 1.77) and in the COVID-19 period (OR 1.99 (95% CI 1.7 to 2.32)) significantly.

In the sample of patients who tested positive for COVID-19 and were hospitalised (n=1742/364 992; 0.48%), 849 were male (48.7%), median age was 64.9 years, median Charlson score was 3 (IQR 1–4) with 737 individuals (42.3%) diagnosed with diabetes, 1467 (84.2%) diagnosed with hypertension and 398 (22.8%) diagnosed with asthma. Mean BMI was 30.2 (±7.5) among patients in the sample. Men had higher comorbidity with 464 (53.7%), as compared with 436 (48.8%), categorised in the highest Charlson score category (table 2). Men also had higher rates of smoking than women (table 3). Men had 1.37 odds of in-hospital death (95% CI 1.09 to 1.72) as compared with women after adjusting for race/ethnicity and age, which was not altered by the addition of clinical comorbidity (OR remained at 1.38 (95% CI 1.08 to 1.76)). However, addition of pulse oximetry on presentation to the hospital attenuated the odds of male death towards the null (1.26 (95% CI 0.99 to 1.62)).

Table 3

Demographic and clinical variables by sex among those hospitalised with a positive COVID-19 test

City-wide male versus female death rates showed a difference of 0.02% between male and female mortality for 2015 and a 0.04% difference in 2016, 2017 and 2018 (online supplemental table 1). City-wide data from the department of health on COVID-19 cases, hospitalisations and deaths from 29 February 2020–1 June 2020 was analysed.28 As shown in online supplemental table 2, there were 0.4% more men who tested positive, 0.2% more men who were hospitalised and 0.11% more men who died from COVID-19 than women. This suggests a relative increase of male mortality during COVID-19 as compared with prior years, portraying the differential impact of COVID-19 on men at a population level.

State-wide and national CDC data on COVID-19 and all-cause deaths for people ≥18 years old and people ≥65 years old was analysed.20 In New York State and nationally, more men died from COVID-19 than women. This difference was greater than the sex differences in all-cause mortality suggesting an even greater male vulnerability to COVID-19 (online supplemental table 3).

Discussion

In this analysis of an ambulatory population of the BMHS we showed a higher proportion of men than women were hospitalised during both pre-COVID-19 and COVID-19 time periods. The higher male risk for hospitalisation was driven by higher comorbidity in the pre-COVID-19 period but not during the COVID-19 period. Men were also at disproportionately higher risk of death during both time periods, and even though their greater risk of death attenuated towards the null in the pre-COVID-19 time period with adjustment for relevant comorbidities it remained significantly elevated compared with females during the COVID-19 period. These data show that the higher comorbidity burden and higher prevalence of smoking among males does not explain the higher risk of COVID-19-related mortality in men. Thus, even though comorbidity at least partially explains the higher male risk for death in the general population, it fails to do so during the COVID-19 surge of 2020. This suggests the unique vulnerability of men to the virus. In our subgroup analysis of those hospitalised patients with an established positive COVID-19 test, we showed that the higher risk of male death was attenuated towards the null only when adjusting for pulseoximetry readings on hospitalisation, suggesting that men were more critically ill with COVID-19 pneumonia than women and this may have been an explanation for their higher death rate.

Using city-wide data from the department of health, and state-wide and national data from the CDC, a higher risk of male COVID-19 cases, hospitalisations and deaths was also observed.19 20 The greater ARDS and respiratory failure rates in males were attributed to a dysregulated immune response in other studies.8 Furthermore, the higher male risk for mortality may be related to sexual dimorphism in ACE2 and TMPRSS2 expression and regulation.8 16 29 To enter a host cell, SARS-CoV-2 binds to the ACE2 surface receptor via its spike protein. Androgen mediated upregulation of TMPRSS2 has been hypothesised to contribute to increased male disease severity with COVID-19.7–9 16 In fact, men on androgen deprivation therapy for prostate cancer have a significant reduction in COVID-19 infection compared with men with prostate cancer not on this therapy.7 Sex-based differences in immune responses have been theorised as another explanation for the increased male mortality with COVID-19. Sex hormones and the X chromosome (which encodes for various immune regulatory genes) have a significant influence on the immune response.9 Kelada et al report that female sex hormones are immunostimulatory.5 Oestrogen contributes to B cell development and induces production of cytokines in response to viral entry while facilitating reduction of viral load.5 High levels of oestrogen and progesterone in females have been shown to suppress pro-inflammatory cytokine production, and stimulate cell-mediated immunity and antibody production all of which contribute to faster viral clearance.1 When hormone replacement therapy was given to postmenopausal women, a decreased mortality from COVID-19 was observed.30 Animal studies demonstrated greater disease severity in female mice after a ovariectomy or administration of an oestrogen receptor antagonist.31 These findings suggest a protective effect of female sex that is not explained by comorbidities.1 Sex differences in vaccine responses have also been demonstrated, with women mounting stronger responses.17 However, our data predate COVID-19 vaccines so that vaccine responsiveness and or vaccine seeking behaviours do not contribute to our findings.

Higher exposure to virus is another potential explanation for the higher male risks observed. Men are more likely to break social distancing rules and less likely to wear masks or perform proper hand hygiene than women.2 5 Also, men smoke more than women which has clear implications on cardiovascular and pulmonary health.32 Such differences in health behaviour may be a result of cultural and societal constructs of what ‘masculine’ behaviour entails.2 To examine the potential role of men’s social behaviour in putting them at risk of higher exposure to, and of a higher inoculum dose when exposed, examining outcomes in a controlled environment, such as nursing homes, where male and female residents have similar exposures, may be instructive. Nursing home data show that male residents die at a greater rate than female residents from COVID-19.33–37 Heras et al reported the mortality rate for COVID-positive nursing home male residents was 42% compared with 12.9% for COVID-19-positive female residents.35 Another study from an Italian nursing home found that the male COVID-19 mortality rate was 46%, while for females, it was 41%, a significant difference.37 Similar findings were present in another study from multiple US nursing homes, where the male COVID-19 mortality rate was also reported as being higher.36 Thus, the role of exposure to virus as an explanation for why men have worse COVID-19 outcomes than women is not compelling.

Limitations of this study include lack of data on deaths occurring outside of the hospital in the pre COVID-19 and COVID-19 time periods, which would shed light on the role of differential self-referrals to the hospital in the observed outcomes. Furthermore, hospital and population data may provide incomplete documentation of clinical and demographic variables. Unfortunately, New York City’s Office of Vital Statistics has not yet completed death data for 2020, which would have provided total death counts in the city by sex. Moreover, our study was an observational analysis and as such any inference is not conclusive. Residual confounding is also a limitation as there are factors that we did not adjust for (eg, environmental factors). Although our database lacked data on male vs female exposures to the virus, we have used the nursing home data in the previous paragraph to illustrate how exposure rates are playing a lesser role than once thought. Lastly, our data may not apply to the current wave of the pandemic. A review of the literature found that reduction in death with newer variants of the COVID-19 virus, and with the protective effects of vaccination and boosters seemed to favour men more than women. Notwithstanding the enhanced mortality benefit of late among men, they remained at higher risk of mortality than women in most reports.38–40

In conclusion, we showed that comorbidity burden and risk behaviours (smoking) are not the main drivers of the higher risk of adverse COVID-19 outcomes in men, but severity of COVID-19 related pneumonia is likely a significant driver. Future studies should look more closely at sex-based differences in physiologic, hormonal and immune mediators in response to COVID-19 infection.

Data availability statement

All data relevant to the study are included in the article or uploaded as online supplemental information. All data relevant to this study are provided upon request to the senior author.

Ethics statements

Patient consent for publication

Ethics approval

Our study was approved by the Albert Einstein College of Medicine Institutional Review Board (IRB number: 2020-12286).

Acknowledgments

We would like to acknowledge Dr. Yungtai Lo for his assistance with some of the statistical analyses.

References

Supplementary materials

  • Supplementary Data

    This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.

Footnotes

  • Contributors Study conception, design, analysis, interpretation of results and drafting of manuscript done by AV, LG. Analysis, interpretation of results by EB and JN. Writing of manuscript, presentation of data by MK. All authors reviewed and approved the final manuscript. LG acts as the guarantor of this study and accepts full responsibility for the work and/or conduct of the study, had access to the data, and controlled the decision to publish.

  • Funding The authors have not declared a specific grant for this research from any funding agency in the public, commercial or not-for-profit sectors.

  • Competing interests None declared.

  • Patient and public involvement Patients and/or the public were not involved in the design, or conduct, or reporting, or dissemination plans of this research.

  • Provenance and peer review Not commissioned; externally peer reviewed.

  • Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise.