Cardiorespiratory fitness and lifestyle on severe COVID-19 risk in 279,455 adults: a case control study
International Journal of Behavioral Nutrition and Physical Activity volume 18, Article number: 135 (2021)
The impact of cardiorespiratory fitness (CRF) and other lifestyle-related factors on severe COVID-19 risk is understudied. The present study aims to investigate lifestyle-related and socioeconomic factors as possible predictors of COVID-19, with special focus on CRF, and to further study whether these factors may attenuate obesity- and hypertension-related risks, as well as mediate associations between socioeconomic factors and severe COVID-19 risk.
Out of initially 407,131 participants who participated in nationwide occupational health service screening between 1992 and 2020, n = 857 cases (70% men, mean age 49.9 years) of severe COVID-19 were identified. CRF was estimated using a sub-maximum cycle test, and other lifestyle variables were self-reported. Analyses were performed including both unmatched, n = 278,598, and sex-and age-matched, n = 3426, controls. Severe COVID-19 included hospitalization, intensive care or death due to COVID-19.
Patients with more severe COVID-19 had significantly lower CRF, higher BMI, a greater presence of comorbidities and were more often daily smokers. In matched analyses, there was a graded decrease in odds for severe COVID-19 with each ml in CRF (OR = 0.98, 95% CI 0.970 to 0.998), and a two-fold increase in odds between the lowest and highest (< 32 vs. ≥ 46 ml·min−1·kg−1) CRF group. Higher BMI (per unit increase, OR = 1.09, 1.06 to 1.12), larger waist circumference (per cm, OR = 1.04, 1.02 to 1.06), daily smoking (OR = 0.60, 0.41 to 0.89) and high overall stress (OR = 1.36, 1.001 to 1.84) also remained significantly associated with severe COVID-19 risk. Obesity- and blood pressure-related risks were attenuated by adjustment for CRF and lifestyle variables. Mediation through CRF, BMI and smoking accounted for 9% to 54% of the associations between low education, low income and blue collar/low skilled occupations and severe COVID-19 risk. The results were consistent using either matched or unmatched controls.
Both lifestyle-related and socioeconomic factors were associated with risk of severe COVID-19. However, higher CRF attenuated the risk associated with obesity and high blood pressure, and mediated the risk associated with various socioeconomic factors. This emphasises the importance of interventions to maintain or increase CRF in the general population to strengthen the resilience to severe COVID-19, especially in high-risk individuals.
The coronavirus disease 2019 (COVID-19) has become a public health emergency worldwide. Among approximately one million laboratory confirmed COVID-19 cases in Sweden, over 57,000 have been hospitalised and more than 14,000 COVID-19 related deaths have been confirmed (up until May 15th, 2021). Mechanisms explaining a higher vulnerability to severe COVID-19 have been linked to inflammation characterized by increased levels of several pro-inflammatory cytokines and the inflammasome . In turn, this has resulted in an inter-individual variation in severity of COVID-19 infection, so that, for example, older age, male gender and one or more comorbidities have been associated with increased risk for hospitalization and mortality due to COVID-19 [2,3,4]. Also, lifestyle-related factors have been linked to COVID-19 severity. Overall/central obesity and hypertension were the first and most frequently reported factors found to be more prevalent in individuals who were hospitalized or died due to COVID-19 [5,6,7,8]. In later papers, physical inactivity has been linked to severe COVID-19 risk [9, 10] and in a small sample of men and women, lower cardiorespiratory fitness (CRF) has also been associated with a higher risk of hospitalisation for COVID-19 .
Only a few studies have investigated the importance of lifestyle factors on severe COVID-19 risk, and it is plausible that, based on previous knowledge, a healthy lifestyle before infection may reduce the risk of severe COVID-19. A positive impact on inflammation and the immune system is one possible mechanistic pathway [12, 13], as low-grade inflammation is considered to be a strong causal factor for chronic diseases such as cardiovascular disease and cancer . Also, the possible impact of a healthy lifestyle on other risk factors, such as overweight/obesity and hypertension [6, 15], may induce protection against severe COVID-19, and regular physical activity (PA) has been suggested as a protective non-pharmacological tool against COVID-19 [12, 16]. However, the research underpinning these assumptions is limited, as are previous studies looking at the importance of and the interaction between different lifestyle-related factors for COVID-19 severity.
Apart from the above predictors, lower socioeconomic status (assessed as, for example, educational level, income or area of residence) has been related to more severe COVID-19 [2, 17, 18]. The subsequent severity of the COVID-19 infection may, however, not only be explained by structural socioeconomic factors, but also by more unfavourable lifestyle habits and poorer health status before infection in individuals with lower socioeconomic status [19, 20]. If and how lifestyle mediates some of the associations seen between socioeconomic factors and severe COVID-19 has not yet been investigated.
The identified knowledge gaps above are addressed in the present study, with the main aim being to study a wide span of lifestyle-related and socioeconomic factors as potential predictors of severe COVID-19, and with special focus on CRF. Secondary aims are to study whether CRF may attenuate obesity- and hypertension-related risk of severe COVID-19, and whether lifestyle-related factors mediate the associations between socio-economic factors and severe COVID-19 risk. The hypotheses are that lifestyle-related (in particular CRF) and socioeconomic factors can predict severe COVID-19, and that variations in lifestyle-related factors mediate a large proportion of the risk of severe COVID-19 associated with socioeconomic factors.
The study is a nested case–control study based on data from the Health Profile Assessment (HPA) database (www.hpi.se). HPAs have been carried out in health services all around Sweden since the middle of the 1970s and is offered to all employees working for a company or an organization connected to occupational or health-related services. An HPA includes a questionnaire about lifestyle and health experiences, measurements of anthropometrics and blood pressure, estimations of maximal oxygen consumption (VO2max) from a submaximal cycle ergometer test, and a person-centred dialogue with a HPA coach.
In February 2021, a total of 407,131 HPAs between 1992 and 2020 were available in the database, and the database was linked to national registries with data on severe COVID-19 (defined as hospitalization, intensive care or death due to COVID-19) using the unique Swedish personal identity number. A total of 857 (0.2%) confirmed cases with severe COVID-19 were identified, including COVID-19 hospitalization (n = 547, 0.1%), intensive care (n = 172, 0.04%) and death (n = 138, 0.03%). Controls were recruited from the same HPA database. All deceased controls before 2019–12-31 according to the national cause of death registry were excluded. To minimize internal drop-out, only participants without severe COVID-19 and with valid data on sex, age, educational level, CRF, body mass index (BMI), exercise and smoking were eligible as controls (a total of n = 278,598). Eligible controls were in general more often women, older, and had higher CRF and lower BMI (see overview of included and excluded participants in Additional file 1). The study was approved by the ethics board at the Stockholm Ethics Review Board (Dnr 2020–02,727). Informed consent was obtained from the participants prior to participation in the HPA. It was not possible to involve participants or the public in the design, conduct, reporting or dissemination plans of our research, due to its retrospective design.
Estimation of VO2max
Measurement of CRF as actual maximal oxygen uptake (VO2max), using a graded test to exhaustion, is limited in the general population for numerous reasons including health risks in non-athletic populations and dependence on laboratory equipment and expertise. Therefore, CRF was assessed as estimated VO2max (estVO2max) from the standardized submaximal Åstrand cycle ergometer test in L·min−1 and also expressed in relative values (ml·min−1·kg−1) . Previous validation studies on adult population samples show small and non-significant mean differences on a group level (− 0.07 L·min−1 95% CI − 0.21 to 0.06) between estVO2max from the Åstrand protocol and direct measured VO2max during treadmill running with an absolute error and coefficient of variance similar to other submaximal tests (SEE = 0.48 L min−1, CV = 18.1%) . To minimize well-known errors with submaximal testing, participants were requested to refrain from vigorous activity the day before the test, consuming a heavy meal and smoking/using snuff three hours and one hour respectively before the test, as well as avoiding physiologic and emotional stress prior to the test. The participant cycled on a calibrated ergometer at an individually adapted submaximal work rate (aiming at a rate of perceived exertion of “Somewhat hard”, 13–14, on the Borg RPE scale) for 6 min to achieve a steady-state pulse assessed during the last minute of cycling. Using the steady- state pulse and the work rate, VO2max was estimated from a sex-specific nomogram, with corresponding age-correction factors .
Body mass was assessed in light-weight clothing using a calibrated scale and to the nearest 0.5 kg. Body height was assessed to the nearest 0.5 cm using a wall-mounted stadiometer. Body mass index (BMI) (kg·m−2) was subsequently calculated. Central obesity was assessed as waist circumference and measured to the nearest 0.5 cm with a tape measure at the midpoint between the top of the iliac crest and the lower margin of the last palpable rib in the mid axillary line after normal exhalation. Systolic and diastolic blood pressure (BP) were measured manually by the standard auscultation method in the seated position after 20 min of resting.
Self-reported and register data
Exercise, commute type, physical work situation, diet habits, alcohol habits, smoking, overall stress, and perceived health were self-reported (see Additional file 2). Highest educational attainment, occupation, income, civil status, and data on country of birth at the time for the HPA were obtained from Statistics Sweden by linking of the participants’ personal identity numbers. Educational attainment was collected from the Swedish education nomenclature 2000 and was categorised into three categories: Elementary school, High school/Vocational education, and University. Each occupation is labelled and defined by a four-digit code according to the Swedish Standard Classification of Occupation . In the present study, occupations were further aggregated according to the first digit into white-collar high-skilled (Major group 1–3), white-collar low-skilled (Major group 4–5), blue-collar high skilled (Major group 6–7) and blue-collar low-skilled (Major group 8–9). For a more detailed description see publication  by Väisänen et al. Total income from employment for the specific year was derived from the Income and Taxation Register. Cases of chronic disease before 2020 were ascertained through the Swedish national patient registry using ICD-coding; C00-D48 to define tumour disease, E10 to E14 for diabetes, I10 to I15 for hypertension, I26 to I28 + J44 to J45 for lung disease and I20 to I25 + I30 to I52 + I60 to I69 for cardiovascular disease.
Severe COVID-19 surveillance
The main outcome in the present study is severe COVID-19, which is defined as either hospital admission, admission to ICU and/or death due to COVID-19. Participants’ personal identity numbers were linked to national registers. Data on hospitalization was obtained from the Swedish National Patient Register, data on intensive care was obtained from the Swedish Intensive Care Registry, and data on death was obtained from the National Cause of Death Register.
Main analyses include cases with sex- and age-matched controls, as the risk of severe COVID-19 is strongly associated with male gender and higher age. Sensitivity analyses using unmatched controls are included in supplementary material (see Additional file 3). For matched analyses, each case was matched to four eligible controls out of the 278,598 eligible controls in the HPA database, with no tolerance in variation between sex or age (per year) between the case and the controls and without replacement of controls in the sampling. For unmatched analyses all eligible controls were used. To compare descriptive characteristics of cases and controls, paired t-test (continuous data), Cochran's Q test (nominal data) and Related-Samples Friedman's Two-Way Analysis of Variance by Ranks (categorical data) were used, and effect sizes as Cohen´s d is presented for continuous data. To compare descriptive characteristics between cases of severe COVID-19, chi-square test of independence with false discovery rate correction for multiple testing (categorical data) and ANCOVA (continuous data) were used. Logistic regression modelling was used to calculate odds ratio (OR) with 95% confidence intervals (95% CI) for different predictors of serve COVID-19. Three models were used and adjusted for an increasing number of variables (see under each table/figure). Model 1 included sex, age and performed year of HPA. Model 2 additionally adjusted for educational level, civil status and country of birth, and Model 3 also adjusted for CRF, BMI, number of chronic diseases, exercise habits, smoking and overall stress. Due to internal drop-out for variables included in Model 2 and 3, two Model 1’s are presented; one Model 1 with OR and 95% CI including all available individuals (labelled Non-complete data in the tables, referred to as Model 1-nc), and a second Model 1 including only individuals with complete data for all adjusting variables (labelled Complete data for all adjusting variables in the tables, referred to as Model 1-c). Further, BMI, waist circumference, blood pressure, estVO2max and income from employment were analyzed as continuous variables as well as after aggregation; BMI was aggregated into normal weight < 25, overweight 25–29.9, obesity 30–34.9 and severe obesity ≥ 35 kg·m−2; large waist circumference was defined as ≥ 88 cm for women and ≥ 102 cm for men, (both BMI and waist circumference were categorized according to recommendations by the world health organization ); high systolic and diastolic BP were defined as ≥ 140 mmHg and ≥ 90 mmHg, respectively; estVO2max was arbitrarily categorised based on multiples of one METs (3.5 ml·min−1·kg−1) into < 25 (very low), 25- < 32 (low), 32- < 46 (moderate) and ≥ 46 (high) ml·min−1·kg−1; and total income from employment into quartiles, percentile 25 = 281,143, percentile 50 = 362,718, percentile 75 = 479,764 Swedish crowns. All above analyses were performed using IBM SPSS (V.18.104.22.168) and Jamovi (The jamovi project (2021) Version 1.6. retrieved from https://www.jamovi.org). Marginal effects plots for severe COVID were calculated by setting the covariates at a mean (for continuous variables) or average (for factor variables) level while varying the focal variables, using R (R Core Team, 2021) and the packages Tidyverse  and ggeffects . Mplus version 8.6  was used to estimate Bayesian parallel mediation models linking socioeconomic indicators to severe COVID-19 via multiple mediators (Fig. 2). Separate models were estimated for each socioeconomic indicator. The highest socioeconomic category was used as the reference category in each model. CRF, BMI, exercise, and stress were treated as continuous variables whereas smoking was dichotomized into daily smoker or never/seldom smoking. We used the proportion of the total effect that is mediated as an effect size measure . The proportion mediated by each mediator was calculated by dividing the specific indirect effect by the total effect (a1b1) / (a1b1 + a2b2 + a3b3+ a4b4+ a5b5+ c’). The total proportion mediated was calculated by dividing the sum of the indirect effects by the total effect (a1b1 + a2b2 + a3b3+ a4b4+ a5b5) / (a1b1 + a2b2 + a3b3+ a4b4+ a5b5+ c’). Models were estimated using four Markov chain Monte Carlo chains and a minimum of 50,000 iterations. The first half of the iterations were discarded as burn-in and the remaining iterations were used to estimate the posterior distribution of the parameters. A stable potential scale reduction factor (PSFR) close to 1 was considered as evidence of chain convergence alongside inspection of trace plots and autocorrelation plots. Indirect effects were evaluated using 95% highest posterior density (HPD) credibility intervals . The credibility interval indicates the probability that the parameter lies between the lower and upper bound of the interval . If an interval did not include zero, the indirect effect was credible. The default non-informative prior specification in Mplus was used.
Characteristics of cases and controls
In the matched analyses, 857 cases of severe COVID-19 and 3426 matched controls were included (for one case, only two exact matched controls were identified). Mean age was 49.9 years (SD 10.7) and 70.4% (n = 603 cases and n = 2 412 controls) were men. The median year that the HPA was performed was 2012 (Q1 2008, Q3 2016) for controls and 2011 (Q1 2006, Q3 2016) for cases. In the unmatched analyses, the mean age for all eligible controls was significantly lower compared to cases (43.7 years (SD 11.6), p < 0.001), and with a significantly lower proportion of men (53.8%, p < 0.001) compared to the matched analyses. The median year that the HPI was performed for all eligible controls was similar to the matched controls 2012 (Q1 2007, Q3 2016).
There were several differences between cases and matched controls for established COVID-19 risk factors, such as cases having higher BMI, blood pressure and presence of comorbidities as well as greater waist circumferences (Table 1). Cases also demonstrated significantly lower estVO2max and more unfavourable exercise patterns. There were also several differences in terms of educational level, country of birth, dietary habits and self-rated health. Moreover, cases with more severe complications from COVID-19 (death vs intensive care or hospitalization, and intensive care vs hospitalization) had significantly lower estVO2max, higher BMI, greater presence of comorbidities and were more often daily smokers (see Table 2).
Impact of lifestyle-related characteristics
Four models were used to quantify independent associations between potential lifestyle-related predictors and severe COVID-19 in the matched analyses, where two Model 1’s (Model 1-nc and Model 1-c) enabled comparative analyses with non-complete and complete data for all adjusting variables in Model 2 and 3, see Table 3.
In terms of CRF, there was a graded increase in odds with lower compared to high levels, OR = 1.62 (95% CI, 1.00 to 2.62) for moderate CRF (32 to < 46 ml·min−1·kg−1) and an approximately two-fold increased odds for low (OR = 2.02, 1.22 to 3.35) and very low fitness (OR = 1.91, 1.09 to 3.34), respectively (Table 3). Similarly, being overweight was associated with two-fold increased odds compared to normal weight (OR = 1.98, 1.53 to 2.56), and obesity and severe obesity was associated with three-fold increased odds (OR = 2.94, 2.13 to 4.07 and OR = 2.98, 1.80 to 4.94 respectively). A larger WC was associated with higher odds in the fully adjusted model, OR = 1.75, 1.20 to 2.55. Presence of chronic disease had a graded increase for every additional diagnosis, OR = 1.88 (95% CI: 1.44 to 2.45) for one chronic disease, and OR = 4.55 (1.83 to 11.33) for 4 to 5 chronic diseases. Neither high systolic nor diastolic blood pressure remained significantly associated with severe COVID-19 after multi-adjustment. Reporting daily smoking (OR = 0.60, 0.41 to 0.89) as well as a high level of stress (OR = 1.36, 1.001 to 1.84) were significantly associated with severe COVID-19 in the fully adjusted model.
Figure 1 presents the predicted probability of severe COVID-19 according to continuous levels of CRF, overall and central obesity (BMI and WC), and systolic and diastolic blood pressure. All obesity and blood pressure measures were attenuated by adjustment for lifestyle variables and CRF, however, these associations with severe COVID-19 remained significant.
Impact of sociodemographic factors
Low education (elementary vs. university) predicted severe COVID-19 (OR = 1.81, 1.28 to 2.54) after multi-adjustment, as did being born outside Sweden vs. being born in Sweden (OR = 2.58, 1.97 to 3.38) (Table 4). No significant associations were seen for civil status, occupational groups, or income quartiles.
In additional sensitivity analyses using unmatched controls (see Additional file 3, Supplement Tables 1 and 2) the odds for severe COVID-19 were higher in men (Model 3, OR = 1.97, 1.62 to 2.40) and with increasing age (Model 3, per year OR = 1.02, 1.01 to 1.03). Individuals ≥ 70 years and 60 to 69 years had higher odds compared to those < 60 years. However, the odds were attenuated by additional adjustments for lifestyle-related factors (Model 2 and 3 adjustment). Further, lifestyle-related and sociodemographic predictors showed similar associations as in the matched analyses, with only marginal variations in both OR and CI, which did not alter the results or conclusions of the unmatched analyses.
The mediation analyses are summarized in Table 5 and Fig. 2. Indirect effects were observed through BMI, CRF and smoking, whereas no credible indirect effects were observed through exercise and stress. The proportion mediated ranged from 12 to 23% for BMI, 9% to 17% for CRF, and 24% to 54% for smoking. Compared to those in the highest socioeconomic category, lower socioeconomic status was related to an increased risk of severe COVID-19 through BMI and CRF, and a lower risk for severe COVID-19 through smoking. The findings were similar across all three socioeconomic indicators. The total proportion mediated across all five mediators ranged from 49 to 86%. Because the mediation models were inconsistent (i.e., they included both positive and negative effects on the dependent variable), we calculated the proportion mediated based on absolute values . Thus, these values represent the proportion of the absolute total effect that was mediated. Mediation analyses using an unmatched sample (N = 279,455) showed similar results (see Additional file 3, Supplement Table 3).
The main results of the present study include strong associations of several lifestyle-related risk factors, including CRF, overweight/obesity, perceived stress, and high blood pressure, with severe COVID-19, even after adjustments for sociodemographic factors and previous diseases. Among patients with severe COVID-19, those with more severe COVID-19 (death vs intensive care or hospitalization, and intensive care vs hospitalization) had lower CRF. In mutually adjusted analyses, higher CRF attenuated some of the risks related to both obesity and hypertension. Low educational level, low income as well as blue collar/low skilled occupations were associated with increased risk of severe COVID-19. However, these associations were, to a relatively large proportion, mediated by CRF, BMI and smoking. The results were consistent when using matched or unmatched controls.
Comparison with other studies
This is, to our knowledge, the first study investigating the association between a wide variation of lifestyle-related risk factors, including CRF, and severe COVID-19. The results are consistent with the few previous existing studies that have found evidence of associations between PA, overweight/obesity and smoking with severe COVID-19 [9, 10]. In 387,109 middle-aged men and women from the UK Biobank, there were 760 cases of hospitalization for COVID-19 . After multi-adjustment, participants reporting no regular PA had a 32% higher risk of hospitalization for COVID-19 compared to those reporting some PA (active but below guidelines) or meeting activity guidelines (≥ 150 min per week of moderate PA or 75 min per week of vigorous PA). Smoking, compared to not smoking, was associated with a 42% higher risk. Further, a lifestyle score was derived including both physical inactivity, smoking, heavy alcohol consumption, and overweight/obesity, which showed a dose-dependent increased risk of hospitalization for COVID-19 partly explained by C-reactive protein levels. Moreover, in 48,440 adult patients with a COVID-19 diagnosis, those who had been consistently inactive (0–10 min of PA per week) in the two years preceding COVID-19 infection, had a significantly higher odds of hospital admission (OR 2.26), admission to intensive care (OR 1.73) and death (OR 2.49) due to COVID-19 compared to patients reporting being inconsistently (10 to < 150 min per week) or consistently (≥ 150 min per week) moderately to vigorously physically active . Although the present study did not find a significant association between PA levels and severe COVID-19, the strong and consistent association of CRF and COVID-19 may be even more important. All previous studies have relied on self-reported PA, which is a subjective measure of recent PA levels containing well-known errors (recall-bias) that permit valid analyses on mainly aggregated PA levels . In the present analyses, CRF was included as a more objective measure of recent PA as well as an indicator of the status of the cardiorespiratory system. This showed a lower risk of severe COVID-19 per ml·min−1·kg−1 with a doubling of risk between the two lowest and the highest categories (< 32 ml·min−1·kg−1 and ≥ 46 ml·min−1·kg−1 respectively). This is similar to a previous report on all-cause mortality and CVD morbidity risk, where decreases of 2.3% and 2.6% per ml increase in estVO2max were seen . Only one previous study has studied the association between recent CRF and COVID-19. In a small sample of patients (n = 246) with positive tests for COVID-19, men (but not women) with lower CRF were more likely to be hospitalized than those with higher CRF . A study using data from military conscript (≈18 years of age) between 1968 and 2005 showed that high CRF at conscript was associated with lower odds of severe COVID-19 later in life .
In the fully adjusted analyses, both perceived stress and smoking remained significantly associated with severe COVID-19. Reporting high overall stress was associated with significantly higher OR (1.36) compared to low stress. This is partly supported by findings from the UK Biobank participants  where a 58% increased risk of hospitalization due to COVID-19 was found among individuals reporting high psychological distress. In contrast to the present results, the association did not remain after full adjustment with comorbidities, other lifestyle variables and socioeconomics. More surprisingly in the present study, smokers had a significantly lower OR (0.60) compared to non-smokers, which adds to equivocal results in the current literature [4, 9]. A hypothesis has been raised that nicotine may have beneficial effects on COVID-19 due to its interaction with the renin-angiotensin and effects on the immunomodulatory system , but further investigation of the mechanisms associated with these findings remains to be elucidated by better controlled studies.
Consistent with other publications [5, 6], both overweight and obesity were associated with a higher risk of severe COVID-19. This could partly be explained by a higher prevalence of metabolic risk factors and low-grade inflammation in overweight/obese individuals, as these have been identified as central mechanisms for a higher vulnerability to severe COVID-19 . Interestingly, a recent paper including over 17 million individuals found similar associations between commonly accepted risk factors (age, male sex, deprivation, obesity, and some comorbidities) for non-COVID (including CVD, cancer, dementia etc.) deaths and for COVID-19 deaths, suggesting that COVID-19 largely mirrors existing risks faced by patients . However in the present study, obesity-risk was at least partly attenuated by CRF. Attenuation by CRF were also seen for central obesity (waist circumference) and high systolic and diastolic blood pressure related risks. These findings are highly clinically relevant and in line with previous studies on cardiovascular disease risk and premature death, where “fat but fit” individuals had significantly better prognoses for cardiovascular outcomes and mortality compared to obese but unfit individuals [38, 39].
There are several suggested mechanisms for the beneficial effects of regular PA and higher CRF levels on both COVID-19 severity per se, as well as attenuation of the obesity- and hypertension-related risks [12, 13, 16, 40]. One is the lower prevalence of obesity and hypertension in more active individuals [6, 15]. Moreover, regular exercise induces a marked increase in several anti-inflammatory cytokines, counteracting the low-grade inflammatory state present in many chronic metabolic diseases (such as obesity and type 2 diabetes) [12, 13]. It also induces a natural immune-protection against more severe COVID-19 by reducing the so-called “cytokine storm” (peaking of pro-inflammatory cytokines including interleukin-6 and tumour necrosis factor-alpha) that ICU-patients with severe COVID-19 experience [41, 42]. Also, regular PA has shown a direct and positive effect on lung function, and the antibody concentration after vaccination is higher among regularly physically active individuals .
There were differences in sociodemographic factors between cases and controls in the present study, which is consistent with previous studies. Among 431,051 British adults, low levels of education, income and area deprivation doubled the risk of hospitalization due to COVID-19, with a 39% higher risk for those with occupations including personal service and sales compared to managers . Across 3135 US counties, the counties with a higher percentage of households with poor housing had a higher incidence of COVID-19, as well as mortality due to COVID-19 . These findings are supported by a large Swedish study, indicating that an educational level only up to elementary school, compared to higher educational levels, was associated with a higher risk for both intensive care and non-intensive care hospitalisation due to COVID-19 . Also, blue-collar workers were significantly less likely to work from home or to change commuting habits in relation to the COVID-19 pandemic, compared to white-collar workers . However, as health status prior to infection seems to heavily impact the severity of COVID-19, we hypothesized that the variation in health lifestyle factors would mediate some of the risk associated with socioeconomic factors. In the mediation analyses, lower socioeconomic status (indicated by education, income, and occupation) was related to an increased risk of severe COVID-19 through higher BMI and lower CRF, whereas lower socioeconomic status was related to a lower risk of severe COVID-19 through smoking. Similar mediation analyses have been performed for cardiovascular disease  and cancer morbidity and mortality , where modifiable factors including BMI and smoking explained between 42 and 46% of the association between low socioeconomic position and the outcomes. The proportion mediated in the current study ranged from 49 to 86%, indicating that the mediators accounted for a relatively large proportion of the association between socioeconomic factors and risk of severe COVID-19. The present indirect effects on severe COVID-19 risk through BMI and CRF highlight factors that could be targeted in interventions to strengthen the resilience for future severe infections.
Strengths and weaknesses of the study
A case–control study is not as powerful as other types of studies in confirming a causal relationship . However, the strengths of this study are the large cohort of different aged women and men with variations in socioeconomic gradients, and the available data on several lifestyle-related factors assessed by standardised methods. Another strength is the highly corresponding results obtained using either the sex- and age-matched controls or all eligible controls in the analyses. The mediation analyses are also a strength, as they highlight processes through which socioeconomic inequalities may influence disease risk. In Sweden, patients from both low and high socioeconomic status have similar access to healthcare, which strengthens the argument for the role of lifestyle factors, including CRF, in preventing severe COVID-19. Limitations of the study include self-reported data regarding lifestyle habits, which risks recall bias . However, questionnaires with categorical answer modes as used in the present study have been reported to provide superior validity compared to open answer modes for PA level . The study design explores associations over time, but does not give information about causality, in this case between lifestyle related and socioeconomic risk factors and severe COVID-19. Moreover, the clinical status of the cases and controls between the time of their HPA and the follow-up period were not monitored. There is a risk of reversed causality due to individuals with a better health status possibly having higher CRF, lower BMI and lower blood pressure. However, the size of the study population made it possible to identify the effect of low CRF, obesity and elevated blood pressure by adjusting for multiple potential confounders and thereby reducing the risk of reverse causality.
Higher CRF was associated with better resilience for severe COVID-19, which is of great clinical value, particularly for high-risk individuals with obesity and/or hypertension. Further, the mediation analyses included in the present paper add important initial evidence of modifiable factors mediating the associations between socioeconomic variables and severe COVID-19. This should shift the focus from structural factors, such as educational level or income per se, having direct effects on disease risk, to instead highlighting and targeting modifiable factors, including CRF and BMI, to increase resilience. This is particularly important as a decrease by 10% (4.2 ml·min−1·kg−1) in CRF has been reported over the last two decades in the Swedish working population . This has been confirmed in international data . During the same time period, the prevalence of obesity and severe obesity has increased by 153% and 86%, respectively . This calls for an urgent need to implement interventions, such as PA on prescription, to increase CRF, preferably specifically targeting high-risk individuals. Further analyses on how sex and age moderate the association between CRF and severe COVID-19 are needed, as are studies including objective measures (e.g. accelerometers) for assessment of PA patterns. Although the mediation analyses highlight processes through which socioeconomic inequalities may influence disease risk, given the correlational nature of the data, these findings need to be replicated in future studies using designs that allow for stronger causal conclusions.
Availability of data and materials
The datasets generated and/or analysed during the current study are not publicly available due to them being the property of the HPI Health Profile Institute, but are available from the corresponding author, email@example.com. Additional information regarding technical details, statistical code, and derived data are also available from the corresponding author.
Body mass index
Coronavirus disease 2019
Estimated maximal oxygen consumption
Health Profile Assessment
- HPD CI:
Highest posterior density credibility interval,
International Classification of Diseases
Maximal oxygen consumption
Ramos-Casals M, Brito-Zeron P, Mariette X. Systemic and organ-specific immune-related manifestations of COVID-19. Nat Rev Rheumatol. 2021;17(6):315–32.
Bergman J, Ballin M, Nordström A, Nordström P. Risk factors for COVID-19 diagnosis, hospitalization, and subsequent all-cause mortality in Sweden: a nationwide study. Eur J Epidemiol. 2021;36(3):287–98.
Li J, Huang DQ, Zou B, Yang H, Hui WZ, Rui F, et al. Epidemiology of COVID-19: A systematic review and meta-analysis of clinical characteristics, risk factors, and outcomes. J Med Virol. 2021;93(3):1449–58.
Williamson EJ, Walker AJ, Bhaskaran K, Bacon S, Bates C, Morton CE, et al. Factors associated with COVID-19-related death using OpenSAFELY. Nature. 2020;584(7821):430–6.
Hamer M, Gale CR, Kivimäki M, Batty GD. Overweight, obesity, and risk of hospitalization for COVID-19: A community-based cohort study of adults in the United Kingdom. Proc Natl Acad Sci USA. 2020;117(35):21011–3.
Huang Y, Lu Y, Huang YM, Wang M, Ling W, Sui Y, et al. Obesity in patients with COVID-19: a systematic review and meta-analysis. Metabolism. 2020;113:154378.
Fang X, Li S, Yu H, Wang P, Zhang Y, Chen Z, et al. Epidemiological, comorbidity factors with severity and prognosis of COVID-19: a systematic review and meta-analysis. Aging. 2020;12(13):12493–503.
Gao M, Piernas C, Astbury NM, Hippisley-Cox J, O'Rahilly S, Aveyard P, et al. Associations between body-mass index and COVID-19 severity in 6·9 million people in England: a prospective, community-based, cohort study. Lancet Diab Endocrinol. 2021;9(6):350–9.
Hamer M, Kivimäki M, Gale CR, Batty GD. Lifestyle risk factors, inflammatory mechanisms, and COVID-19 hospitalization: a community-based cohort study of 387,109 adults in UK. Brain Behav Immun. 2020;87:184–7.
Sallis R, Young DR, Tartof SY, Sallis JF, Sall J, Li Q, et al. Physical inactivity is associated with a higher risk for severe COVID-19 outcomes: a study in 48 440 adult patients. Br J Sports Med. 2021. Published Online First: 13 April 2021. https://doi.org/10.1136/bjsports-2021-104080.
Brawner CA, Ehrman JK, Bole S, Kerrigan DJ, Parikh SS, Lewis BK, et al. Inverse relationship of maximal exercise capacity to hospitalization secondary to coronavirus disease 2019. Mayo Clin Proc. 2021;96(1):32–9.
Zbinden-Foncea H, Francaux M, Deldicque L, Hawley JA. Does high cardiorespiratory fitness confer some protection against proinflammatory responses after infection by SARS-CoV-2? Obesity (Silver Spring, Md). 2020;28(8):1378–81.
Jakobsson J, Malm C, Furberg M, Ekelund U, Svensson M. Physical activity during the coronavirus (covid-19) pandemic: prevention of a decline in metabolic and immunological functions. Front Sports Act Living. 2020;2:57.
Kraus WE, Powell KE, Haskell WL, Janz KF, Campbell WW, Jakicic JM, et al. Physical activity, all-cause and cardiovascular mortality, and cardiovascular disease. Med Sci Sports Exerc. 2019;51(6):1270–81.
Pranata R, Lim MA, Huang I, Raharjo SB, Lukito AA. Hypertension is associated with increased mortality and severity of disease in COVID-19 pneumonia: A systematic review, meta-analysis and meta-regression. J Renin Angiotensin Aldosterone Syst. 2020;21(2):1470320320926899.
Nigro E, Polito R, Alfieri A, Mancini A, Imperlini E, Elce A, et al. Molecular mechanisms involved in the positive effects of physical activity on coping with COVID-19. Eur J Appl Physiol. 2020;120(12):2569–82.
Batty GD, Deary IJ, Luciano M, Altschul DM, Kivimäki M, Gale CR. Psychosocial factors and hospitalisations for COVID-19: prospective cohort study based on a community sample. Brain Behav Immun. 2020;89:569–78.
Ahmad K, Erqou S, Shah N, Nazir U, Morrison AR, Choudhary G, et al. Association of poor housing conditions with COVID-19 incidence and mortality across US counties. PloS one. 2020;15(11):e0241327.
Carter AR, Gill D, Davies NM, Taylor AE, Tillmann T, Vaucher J, et al. Understanding the consequences of education inequality on cardiovascular disease: mendelian randomisation study. Bmj. 2019;365:l1855.
Nejatinamini S, Godley J, Minaker LM, Sajobi TT, McCormack GR, Cooke MJ, et al. Quantifying the contribution of modifiable risk factors to socio-economic inequities in cancer morbidity and mortality: a nationally representative population-based cohort study. Int J Epidemiol. 2021. Online first. https://doi.org/10.1093/ije/dyab067.
Astrand I. Aerobic work capacity in men and women with special reference to age. Acta Physiol Scand Suppl. 1960;49(169):1–92.
Björkman F, Ekblom-Bak E, Ekblom Ö, Ekblom B. Validity of the revised Ekblom Bak cycle ergometer test in adults. Eur J Appl Physiol. 2016;116(9):1627–38.
Stastistics Sweden. Standard för svensk yrkesklassificering (SSYK) [Standard for Swedish occupational classification (SSYK)]. (In Swedish) [cited 2020 May 18]. [Available from: http://www.scb.se/dokumentation/klassifikationer-och-standarder/standard-for-svenskyrkesklassificering-ssyk/.
Väisänen D, Kallings LV, Andersson G, Wallin P, Hemmingsson E, Ekblom-Bak E. Lifestyle-associated health risk indicators across a wide range of occupational groups: a cross-sectional analysis in 72,855 workers. BMC Public Health. 2020;20(1):1656.
World Health Organization. Waist circumference and waist-hip ratio: report of a WHO expert consultation. Geneva. December 2008;8–11:2011.
Wickham H, Averick M, Bryan J, Chang W, McGown L, Francois R, et al. Welcome to the Tidyverse. The Journal of Open Source Software. 2019;https://doi.org/10.21105/joss.01686.
Lüdecke D. ggeffects: Tidy Data Frames of Marginal Effects from Regression Models. The Journal of Open Source Software. 2018;https://doi.org/10.21105/joss.00772.
Muthén LK, Muthén BO. Mplus User’s Guide. 8th ed. Los Angeles: Muthén & Muthén. 1998–2017.
MacKinnon DP. Introduction to statistical mediation analysis. 1st ed. New York: Taylor & Francis Group; 2008.
Miocevic M, MacKinnon DP, Levy R. Power in Bayesian Mediation Analysis for Small Sample Research. Struct Equ Modeling. 2017;24(5):666–83.
Yuan Y, MacKinnon DP. Bayesian mediation analysis. Psychol Methods. 2009;14(4):301–22.
Ekblom O, Ekblom-Bak E, Bolam KA, Ekblom B, Schmidt C, Soderberg S, et al. Concurrent and predictive validity of physical activity measurement items commonly used in clinical settings–data from SCAPIS pilot study. BMC Public Health. 2015;15:978.
Ekblom-Bak E, Ekblom B, Soderling J, Borjesson M, Blom V, Kallings LV, et al. Sex- and age-specific associations between cardiorespiratory fitness, CVD morbidity and all-cause mortality in 266.109 adults. Prev Med. 2019;127:105799.
Af Geijerstam A, Mehlig K, Borjesson M, Robertson J, Nyberg J, Adiels M, et al. Fitness, strength and severity of COVID-19: a prospective register study of 1 559 187 Swedish conscripts. BMJ Open. 2021;11(7):e051316.
Farsalinos K. BA, Niaura R. Smoking, vaping and hospitalization for COVID-19. Qeios. 2020. https://doi.org/10.32388/Z69O8A.13.
Popkin BM, Du S, Green WD, Beck MA, Algaith T, Herbst CH, et al. Individuals with obesity and COVID-19: A global perspective on the epidemiology and biological relationships. Obesity reviews : an official journal of the International Association for the Study of Obesity. 2020;21(11):e13128.
Bhaskaran K, Bacon SCJ, Evans SJW, Bates CJ, Rentsch CT, MacKenna B, et al. Factors associated with deaths due to COVID-19 versus other causes: population-based cohort analysis of UK primary care data and linked national death registrations within the OpenSAFELY platform. The Lancet Regional Health. 2021;6:https://doi.org/10.1016/j.lanepe.2021.100109.
Lavie CJ, Carbone S, Kachur S, OʼKeefe EL, Elagizi A. Effects of physical activity, exercise, and fitness on obesity-related morbidity and mortality. Curr Sports Med Rep. 2019;18(8):292–8.
Sanchez-Lastra MA, Ding D, Dalene KE, Ekelund U, Tarp J. Physical activity and mortality across levels of adiposity: a prospective cohort study from the UK Biobank. Mayo Clin Proc. 2021;96(1):105–19.
Chastin SFM, Abaraogu U, Bourgois JG, Dall PM, Darnborough J, Duncan E, et al. Effects of Regular Physical Activity on the Immune System, Vaccination and Risk of Community-Acquired Infectious Disease in the General Population: Systematic Review and Meta-Analysis. Sports Med (Auckland, NZ). 2021:51(8):1673–86.
Mehta P, McAuley DF, Brown M, Sanchez E, Tattersall RS, Manson JJ, et al. COVID-19: consider cytokine storm syndromes and immunosuppression. Lancet. 2020;395(10229):1033–4.
Nieman DC, Pence BD. Exercise immunology: Future directions. J Sport Health Sci. 2020;9(5):432–45.
Blom V, Lonn A, Ekblom B, Kallings LV, Vaisanen D, Hemmingsson E, et al. Lifestyle Habits and Mental Health in Light of the Two COVID-19 Pandemic Waves in Sweden, 2020. Int J Environ Res Public Health. 2021;18(6):3313.
Munnangi S, Boktor SW. Epidemiology Of Study Design. 1st ed. Treasure Island: StatPearls Publishing; 2021.
Ainsworth B, Cahalin L, Buman M, Ross R. The current state of physical activity assessment tools. Prog Cardiovasc Dis. 2015;57(4):387–95.
Olsson SJ, Ekblom O, Andersson E, Borjesson M, Kallings LV. Categorical answer modes provide superior validity to open answers when asking for level of physical activity: A cross-sectional study. Scandinavian journal of public health. 2016;44(1):70–6.
Ekblom-Bak E, Ekblom O, Andersson G, Wallin P, Soderling J, Hemmingsson E, et al. Decline in cardiorespiratory fitness in the Swedish working force between 1995 and 2017. Scand J Med Sci Sports. 2019;29(2):232–9.
Lamoureux NR, Fitzgerald JS, Norton KI, Sabato T, Tremblay MS, Tomkinson GR. Temporal trends in the cardiorespiratory fitness of 2,525,827 adults between 1967 and 2016: a systematic review. Sports Med. 2019;49(1):41–55.
Hemmingsson E, Ekblom O, Kallings LV, Andersson G, Wallin P, Soderling J, et al. Prevalence and time trends of overweight, obesity and severe obesity in 447,925 Swedish adults, 1995–2017. Scand J Public Health. 2021;49(4):377–83.
Thanks go to the Swedish Heart-Lung Foundation and the Swedish Military Forces Research Authority for funding this study. We also thank the Swedish Intensive Care Registry, a non-profit organisation, for sharing data on intensive care admission.
This research was funded by The Swedish Heart–Lung Foundation, grant number 20200564 and The Swedish Military Forces Research Authority (Grant # AF 922 0915). Author AS is supported with an international postdoc grant from The Swedish Research Council grant number 2017–00273. The study sponsors had no role in designing the study, the collection, analysis, or interpretation of data, writing the report, or in the decision to submit the article for publication. The researchers and their work have been independent from the funders and all authors, external and internal, had full access to all of the data (including statistical reports and tables) and take responsibility for the integrity of the data and the accuracy of the data analysis.
Ethics approval and consent to participate
The protocols used were approved by the institutional review boards of the institutions involved in this study and all participants provided informed consent prior to data collection. The study was approved by the Stockholm Ethics Review Board (Dnr 2020–02727), and adhered to the Declaration of Helsinki.
Consent for publication
Author GA (responsible for research and method) and PW (CEO and responsible for research and method) are employed by HPI Health Profile Institute. The remaining authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplement figure 1. Flow chart of included and excluded participants. Contains a flow chart of included and excluded cases and controls.
Questions from Health Profile Assessment. Contains the questions from the Health Profile Assessment for self-reported data in the study.
Supplement Tables. Contains supplement table 1 to 3.
About this article
Cite this article
Ekblom-Bak, E., Väisänen, D., Ekblom, B. et al. Cardiorespiratory fitness and lifestyle on severe COVID-19 risk in 279,455 adults: a case control study. Int J Behav Nutr Phys Act 18, 135 (2021). https://doi.org/10.1186/s12966-021-01198-5