Modeling the relative influence of socio-demographic variables on post-acute COVID-19 quality of life: an application to settings in Europe, Asia, Africa, and South America

Introduction

Long-term COVID-19 sequelae have resulted in a pressing public health crisis since early 2020. A diagnosis of long COVID is defined by the World Health Organization as unexplainable symptoms which persist at least three months after an infection, and occur over the span of two or more months.¹ The widespread presence and impacts of this condition have been immense: a multinational study found that nearly half of individuals who were previously infected with SARS-CoV-2 went on to experience long-term symptoms around four months post infection²; further, an estimated 59% of previously infected subjects reported a reduced quality of life (QoL), and total disability-adjusted life years (DALYs) for hospitalized patients have been estimated to be as high as 642.8 DALYs/1000 individuals.^3,4 Prior work has focused on identifying a myriad of clinical risk factors for long-term COVID-19 consequences, including co-infections and pre-existing conditions, a number of laboratory measures, severe acute infection, vaccination, age, and sex.^5–14 Conditions that have been consistently identified as key correlates of long-term sequelae include obesity, asthma and other pulmonary diseases, chronic cardiac disease, diabetes, tuberculosis (TB), liver disease, lung disease, and reactivated Epstein-Barr infection.^5–8,10

Beyond clinical factors, social vulnerabilities are often critical determinants of differential disease burden overall, with such inequities attributed to broader challenges in access to disease prevention and management services and an array of health-limiting exposures, including but not limited to food and housing insecurity, financial discrimination, and air pollution.^15–19

While there have been efforts to examine social factors potentially linked to long-term symptoms of COVID-19, findings on these relationships have been somewhat mixed.^{6,8,10,12,14,20–26} For instance, a 2021 study in the United Kingdom (UK) found that living in high-deprivation settings was associated with both higher and lower odds of symptom persistence, depending on the measure of deprivation index used, and a 2021 study in Michigan (USA) found that lower income was both significantly associated and not associated with long COVID symptoms’ prevalence, depending on the post-illness duration considered.^8,21 It is also important to highlight that many of these studies rely on self-reported binary measures of COVID recovery or continued symptomatology, some of which may be subjective classifications, with potential between-group differences in the tendency to report such experiences.

Given this context, we aimed to complement existing efforts centered on uncovering disparities in long-term COVID outcomes, leveraging a large dataset from a prospective, observational, multinational study of hospitalized and non-hospitalized COVID patients with post-infection follow-up data, focusing on Norway, the United Kingdom (UK), India, Brazil, Russia, and South Africa. Specifically, we formally assessed the relationship between a diverse group of biological and social exposures, and our long-term long COVID QoL measure, reasoning that factors like socio-economic status (SES) would matter as much or more than commonly considered comorbidities, as has been recently illustrated for related outcomes, such as “healthy aging”.²⁷ We further evaluated the extent to which clinical intermediates contribute to any observed disparities, hypothesizing that they may only partially explain these effects. Our analysis applies a similar mediation-centered lens to that of Vahidy et al.²⁵ and Lu et al.²⁸ However, rather than focusing on evaluating whether various mediators can independently explain the effects of a given social factor on long COVID risk²⁵ or on how much social factors mediate disparities in all-cause mortality²⁸ we measured the degree to which the effects of social variables on long COVID-associated QoL cannot be explained by comorbidities.

Methods

This study uses data from the International Severe Acute Respiratory and emerging Infection Consortium’s (ISARIC) multi-cohort consortium.²⁹ This prospective study across 76 countries collected demographic and illness-associated data during acute SARS-CoV-2 infection, with a subset of sites assessing participants 3 months (+/-1) following infection and 3-6 months thereafter. Recruitment was aimed at clinical settings. Complete details on the study design and recruitment procedures can be found in the published follow-up protocol.³⁰

We selected representative locations defined by varying pandemic experiences and demographic/socio-cultural contexts, with data available on indicators of socio-economic status (educational attainment and employment status), resulting in the following individual cohorts: Norway (n=1672) and the UK (n=1064). We combined data from India, Brazil, Russia, and South Africa into one cohort representing upper and lower middle-income countries (MICs), as defined by the World Bank³¹. Sample sizes were insufficient (n<1000) for each of these locations individually (n=264, n=125, n=57, for India, South Africa, and Brazil, respectively), with the exception of Russia (n=1155). Our study incorporates information on two outcome measures capturing long COVID experiences, self-reported continued symptomatology and post-illness quality of life, thereby reducing the influence of between-group reporting patterns that may exist within either of these measures.

Long-term health utility values were obtained using standard quality of life-adjustment estimation procedures, based on subjects’ responses to the EQ-5D-5L survey included in the ISARIC follow-up case report forms.³² The purpose of the EQ-5D-5L survey is to elicit self-reported rankings on five dimensions of health (mobility, self-care, usual activities, pain/discomfort, and anxiety/depression), representing the intensity of problems experienced along that dimension. All possible rankings for each of the five dimensions were assigned weights, which quantify preferences for different health states. Participants were asked to report these five rankings based on their experiences prior to their COVID infection and in the present, yielding two sets of rankings. Weights were sourced from studies conducted in the country in question or in a related country (see Supplementary Appendix Table S1).

We subset each country cohort to only subjects reporting, at any follow-up survey, at least one long COVID-associated symptom that was not present prior to illness. Utility scores were computed following standard practice and then time-transformed (Supplementary Appendix: Methods). Similar to Sandmann et al.¹⁴, we used a measure of quality-adjusted life days or QALDs, as opposed to quality-adjusted life years (QALYs), given the limited variability across subjects expected for the latter measure, where utility scores would be simply scaled by one year, i.e., the assumed duration, based on other study findings.^33,34 We additionally chose to focus on quality of life at least three months following infection, or quality of life in the ‘present’ as recorded in each follow-up survey. We did not consider QALD differences post-versus prior to illness, given that issues with recalling experiences several months in the past are likely to result in erroneous measures of pre-COVID QALDs.

We incorporated data on age, socio-demographic variables (female sex at birth and socio-economic status proxies), a set of clinical comorbidities and risk behaviors, COVID-19 severity, antiviral treatment, and vaccination status (Supplementary Appendix: Methods). For all countries included in our analysis, variables directly reflecting socio-economic status, such as income/wealth or job category are not available, so we generated a series of associated indicators depending on the data available. For Norway, we used quintiles of educational attainment (years). For the UK and the combined MIC cohort, we used information on pre-illness employment status (Supplementary Appendix Table S2).

For a subset of variables with incomplete data, we used multiple imputation by chained equations (MICE) to assign values to missing entries (Supplementary Appendix: Methods). We cannot make a definitive conclusion on the missingness mechanisms for sex at birth, i.e., we have no reason to believe that females compared to males are less likely to self-report their sex, which would result in a missing not at random (MNAR) mechanism. Thus, the use of MICE, which assumes a missing (completely) at random mechanism, appears justified. However, vaccination status and socio-economic indicators of educational attainment and employment status are more likely to be MNAR due to social stigmas and fear of judgment. Nonetheless, a relatively low fraction of most of these values were missing, with the exception of the vaccination status indicator in the UK and combined MIC cohorts and antiviral treatment in the Norwegian and combined MIC cohorts (Supplementary Appendix: Methods). Consequently, to avoid imposing significant bias, we chose not to impute these specific variables and excluded them from their respective variable sets. Finally, as the distribution of sex at birth in our sampled populations non-negligibly deviated from that of the underlying population in the individual cohorts, i.e., Norway and the UK, we conducted sensitivity analyses where subjects were assigned post-stratification weights according to the raking method (Supplementary Appendix: Methods).^35–37

To identify social predictors of long COVID QALDs, we implemented a series of random forest ensemble learners for each country, fit to all available clinical and demographic data, where variables were either treated individually (RF #1), grouped algorithmically via hierarchical clustering (RF #3), or pre-grouped based on subject matter knowledge (RF #2). We implemented a pre-grouped procedure, which incorporates subject matter knowledge, as an alternative to model-grouped approaches agnostic to such context.

For RF #1, we trained a simple random forest treating each variable as individual predictors over an ensemble of decision trees. For RF #2, we manually grouped comorbidities depending on similar clinical categorizations or latent mechanisms (Supplementary Appendix Table S3), applying multiple factor analysis³⁸ to the groups to provide a summary measure for each prior to implementing the random forest (Supplementary Appendix: Methods).For RF #3, we applied a combined “variable clustering/variable selection” method introduced by Chavent et al.³⁹, i.e., Clustering of Variables-Variable Selection Using Random Forest (CoV-VSURF), which enables related variables to be grouped prior to being used as inputs in training a random forest.³⁹

To estimate the natural direct effects (NDE) and natural indirect effects (NIE) of binary educational attainment and employment status categorizations or female sex on long COVID QALDs in each cohort, we applied a flexible semi-parametric statistical approach.^40,41 We note that the NDE and NIE are, under well-studied assumptions, interpreted causally, as defined by Pearl⁴² and Robins and Greenland⁴³, although our “exposures” of interest are not directly intervenable as they are social constructs.⁴⁴ For Norway, our binary proxy of SES was high/low educational attainment, where ‘high’ encompassed quintiles 3-5 and ‘low’ encompassed quintiles 1 and 2. For the UK and the combined MIC cohorts, our binary proxy of SES was high/low employment status. For the UK, ‘high’ encompassed full-time employment and ‘low’ encompassed all other employment status categories (furloughed, part-time employment, student, and the unemployed), with the exception of retirement. For the combined MIC cohort, ‘high’ encompassed full-time employment and ‘low’ encompassed all other employment status categories, where retirees are included due to notable disparities in representation of full-time employees vs all other categories in this cohort.

Estimates of the NDE and NIE are obtained using targeted maximum likelihood estimation.^40,45 Through this procedure, “nuisance parameters”^40,41,45,46 are first estimated via ensemble machine learning or super learning, then updated in such a way to diminish bias and inform a more correct characterization of uncertainty for the target parameter of interest (Supplementary Appendix: Methods).^40,41 In this context, nuisance parameters include the conditional probability of each of our binary SES proxies or female sex, dependent on confounders; the conditional probability of each of our binary SES proxies or female sex, dependent on confounders and mediators; and the conditional mean of long COVID QALDs dependent on each of our binary SES proxies or female sex, the confounders, and mediators.

To obtain these estimates, we incorporated a multi-model learner consisting of an assortment of parametric and nonparametric approaches (Supplementary Appendix: Methods). The causal structures assumed here are depicted in the Directed Acyclic Graphs (DAGs) presented in Supplementary Appendix Figures S1-2. Mediators were selected following a literature search of studies distinguishing key long COVID predictors^11,12, consistent with the data available for each country, resulting in the following final mediator sets: asthma, chronic cardiac disease (not hypertension), hypertension, chronic pulmonary disease (not asthma), type 1 diabetes (T1D), type 2 diabetes (T2D), diabetes (type not specified), psychological disorder, smoking, and vaccination status (Norway), asthma, chronic cardiac disease (not hypertension), hypertension, chronic pulmonary disease (not asthma), T2D, diabetes (type not specified), psychological disorder, ischemic heart disease, smoking, obesity, and antiviral treatment (UK) and asthma, chronic cardiac disease (not hypertension), hypertension, chronic pulmonary disease (not asthma), T1D, T2D, diabetes (type not specified), smoking and obesity (combined MIC cohort).

Here, the NDE broadly describes the effect of a given social variable on long COVID QALDs not operating through the included mediators while the NIE specifically describes the effect of that social variable on long COVID QALDs through the included mediators. We define the proportion non-mediated as the % contribution of NDE to the total effect (TE), reflecting the share of the TE of the social variable on long COVID QALDs that cannot be explained by the clinical intermediates, consistent with the social disparities conceptual framing provided in Bellavia et al., 2018.⁴⁷ For the combined MIC cohort, we additionally conducted country-specific analyses to assess whether the directionality of effects were generally consistent compared to their analogues in the overall analysis.

Data availability statement

The data that underpin this analysis are highly detailed clinical data on individuals hospitalised with COVID-19. Due to the sensitive nature of these data and the associated privacy concerns, they are available via a governed data access mechanism following review of a data access committee. Data can be requested via the IDDO COVID-19 Data Sharing Platform (http://www.iddo.org/covid-19). The Data Access Application, Terms of Access and details of the Data Access Committee are available on the website. Briefly, the requirements for access are a request from a qualified researcher working with a legal entity who have a health and/or research remit; a scientifically valid reason for data access which adheres to appropriate ethical principles.The full terms are at: https://www.iddo.org/document/covid-19-data-access-guidelines. A small subset of sites who contributed data to this analysis have not agreed to pooled data sharing as above. In the case of requiring access to these data, please contact the corresponding author in the first instance who will look to facilitate access.

All code (with the exception of code used to process the individual datasets) is publicly available at: https://github.com/goshgondar2018/social_long_covid.

Results

Combined MIC cohort (India, Brazil, Russia, and South Africa)

The mean age of participants was 58.4 years (SD: 14.4 years). Females outnumbered males, constituting 57% of the cohort. Hypertension was the most frequently reported comorbidity (58%). Estimated long COVID QALDs were greatest in this cohort (median=346; Interquartile range (IQR): 316-365), when compared to the Norway and UK cohorts.

View this table:

• View inline
• View popup

Table 1: Summary of demographic variables (excluding SES proxies) and common comorbidities in the final study populations for each cohort, post-missing data imputation. Note the combined middle-income country cohort consists of the following countries: India, Brazil, Russia, and South Africa

Employment status was markedly skewed towards full-employment (53%) and retirement (35%), with unemployment (7%) and part-time employment (2%) seeing less representation. Carers and students were least represented (2% and 0.3%). Long COVID QALDs were highest among students, full-time employees, and carers, and lowest among those in the retired and unemployed categories. Males reported marginally higher long COVID QALDs than females (p<0.0001).

Age dramatically eclipsed all other variables in terms of its predictive power for long COVID QALDs in this cohort, followed by employment status indicators, chronic neurological disorder, and hypertension. Country fell next in the rankings, indicating that underlying cross-country differences in long COVID QALD trends that are not tied to the variables under interest merit consideration, validating our adjustment by country in the subsequent mediation analyses as a source of confounding (Figure 1a). The acute COVID-19 severity indicator was found to be unimportant. RF #2 supported these findings. The cluster containing solely age led the rankings, followed by the principal components (PCs) of the cluster containing socio-demographic variables employment status and sex, and the PCs of the cluster containing psychological disorder and chronic neurological disorder (Figure 1c). Similarly, for RF #3, age, chronic cardiac disease (not hypertension), chronic neurological disorder, as well as dementia, employment status indicators, hypertension, rheumatological disorder, and sex led the set of most frequently selected variables (Figure 1b). However, there were no perceptible differences in variable selection across nearly all predictors, indicating that the CoV-VSURF procedure may not be well-suited for variable rankings in a combined MIC cohort.

• Download figure
• Open in new tab

Figure 1a.

Estimated variable importance measures, i.e. % increase in mean squared error or MSE, from individual random forest implementation (RF #1) for the combined middle-income country cohort.

• Download figure
• Open in new tab

Figure 1b.

Number of times (frequency) each variable appears in clusters selected for each CoV-VSURF run (RF #3) for the combined middle-income country cohort.

• Download figure
• Open in new tab

Figure 1c.

Estimated variable importance measures, i.e. % increase in mean squared error or MSE, from pre-grouped random forest implementation (RF #2) for the combined middle-income country cohort. Rows indicate cluster names (a full list of variables belonging to each cluster can be found in Supplementary Table S3) and corresponding principal components, if the cluster consists of multiple variables. PC1 denotes principal component 1 and PC2 denotes principal component 2.

While we cannot directly comment on the directionality of associations based on these results, the correlation matrix for all variables in the combined MIC cohort (Supplementary Figure S5) indicates that QALDs were independently negatively associated with age, the acute COVID-19 severity indicator, and all top predicting comorbidities, particularly for chronic neurological disorder, hypertension, and rheumatological disorder. Subjects with those comorbidities present reported lower long COVID QALDs.

Full-time employment was associated with higher long COVID QALDs compared to all other employment status categories. 8.95 (95% CI: 1.82, 16.1) more long COVID QALDs were expected among subjects self-reporting full-time employment compared to all other employment categories due to the direct effect of self-reported full-time employment on long COVID QALDs (the NDE). An additional 4.03 (0.82, 7.24) long COVID QALDs were expected among subjects self-reporting full-time employment compared to all other employment categories, due to the indirect effect of self-reported full-time employment on long COVID QALDs through the intermediates (the NIE). The proportion non-mediated was estimated to be 0.74 (0.46, 1), i.e., 74.2% of the effect of full-time employment versus all other employment status categories on long COVID QALDs could not be explained by the included mediators and must thus be attributed to broader structural factors or other mechanisms. Female sex was associated with lower COVID QALDs, with an estimated NDE of -5.24 (-10.8, 0.33), NIE of -0.85 (-2.59, 0.90), and proportion non-mediated of 0.84 (0.62, 1), indicating that a high fraction of the sex-at-birth effect cannot be tied to mediation by the considered clinical intermediates. We note that the upper and lower bounds of the confidence intervals (CIs) for the proportions mediated for all cohorts are not necessarily exactly aligned with what we expect based on the bounds of the CIs for the corresponding NDEs and NIEs. The relationship between these measures may not hold exactly for computational reasons, because CIs are estimated separately for each of the measures using cross-validation, which may introduce technical noise due to the randomness inherent in sample splitting. Finally, from the country-specific analyses, we observed the same directionality (positive), albeit at far reduced precision, for the NDE for the binary high SES indicators in all countries as in the combined analysis, but opposing directionality (negative) for the NIE for the binary high SES indicators in India and Brazil (Supplementary Figure S9).

Nonetheless, the TE for both these countries was positive overall, suggesting that, as in our main analysis, high SES is associated with increased long COVID QALDs, but that the mediators have an opposing intermediary effect (Supplementary Figure S9). Importantly, for all countries except South Africa, the direct effect comprises the majority of the total effect, although to differing degrees (Supplementary Figure S9). For the female sex ‘exposure’, we found consistent directionality for all NDEs and NIEs (negative) as in the combined analysis, with the exception of the estimated NDE for South Africa, and NDEs exceeding NIEs in magnitude (Supplementary Figure S9). For Brazil, NDEs and NIEs could not be estimated because of the large discrepancy between males and females (9 versus 84).

Norway

The Norwegian cohort was the youngest, with a mean age of 51.8 years (SD: 13.6 years). Most participants self-identified as white (93%), with non-white race/ethnic minorities notably underrepresented when compared to the approximate race/ethnicity distribution of Norway.⁴⁸ As a consequence of this imbalance and sparse numbers in some ethnic groups, we concluded that any assessment of the role of race/ethnicity in shaping long COVID QALDs would not be justified, omitting it from our social variables under consideration. Females constituted a majority of this cohort (68%) compared to males and are thus overrepresented in relation to Norway’s underlying population.⁴⁹

The most commonly reported comorbidity was asthma (22%). From the total cohort, 50% reported receiving at least one dose of any COVID-19 vaccine. The median (IQR) value for estimated long COVID QALDs was 345 (313-360).

There was no broadly consistent trend in long COVID QALDs across increasing quintiles of educational attainment, although the lowest mean QALDs occurred in the bottom two quintiles. We observed the greatest differences in long COVID QALDs between, in order of increasing magnitude, quintiles 3 and 1, quintiles 5 and 1, and quintiles 4 and 1. Estimated long COVID QALDs in males slightly exceeded those in females (p<0.0001).

Among the leading individual predictors of long COVID QALDs in Norway, anxiety/depression dominated the rankings, followed by educational attainment, rheumatological disorder, and age (Figure 2a). Sex fell lower in the rankings. For RF #2, we observed that the first and second PCs of the cluster containing all socio-demographic variables, i.e., educational attainment indicators and sex, ranked below the first and second PCs of the cluster containing psychological disorder and chronic neurologic disorder (Figure 2c). RF #3 largely corroborated these orderings, where psychological disorder, rheumatological disorder, chronic neurological disorder, and asthma were the most consistently selected variables within identified important clusters, followed by educational attainment (in years) and a dummy educational attainment indicator for quintile 5 (vs 1) (Figure 2b). Our sensitivity analyses applying the sex-based population corrections produced nearly identical findings for RF #1 and RF #3 (Supplementary Figure S12a and c). As expected, due to its greater intrinsic variability across model runs, the model-grouped CoV-VSURF procedure provided somewhat contrasting rankings, with variables like T2D and diabetes (type not specified) more frequently selected, and age and sex less frequently selected, than in the main analysis (Supplementary Figure S12b). However, the top-ranking variables – psychological disorder, rheumatological disorder, and chronic neurological disorder – were preserved, and educational attainment indicators retained a similarly high ranking (Supplementary Figure S12b).

• Download figure
• Open in new tab

Figure 2a.

Estimated variable importance measures, i.e. % increase in mean squared error or MSE, from individual random forest implementation (RF #1) for Norway.

• Download figure
• Open in new tab

Figure 2b.

Number of times (frequency) each variable appears in clusters selected for each CoV-VSURF run (RF #3) for Norway.

• Download figure
• Open in new tab

Figure 2c.

Estimated variable importance measures, i.e. % increase in mean squared error or MSE, from pre-grouped random forest implementation (RF #2) for Norway. Rows indicate cluster names (a full list of variables belonging to each cluster can be found in Supplementary Table S3) and corresponding principal components, if the cluster consists of multiple variables. PC1 denotes principal component 1 and PC2 denotes principal component 2.

QALDs were found to negatively correlate with all top comorbidities and moderately positively correlate with age and the vaccination status indicator (Supplementary Figure S6).

We estimated that falling in the top two quintiles of educational attainment was associated with 12.3 (6.49,18.2) additional long COVID QALDs, on average, via the NDE and 0.67 (-0.98, 2.32) additional long COVID QALDs, on average, via the NIE. The corresponding proportion non-mediated was estimated to be 0.90 (0.77, 1). We obtained consistent directionality in findings for pairwise comparisons of quintiles 3 and 1, 4 and 1, and 5 and 1, with the greatest proportion non-mediated observed for the quintile 5 vs 1 comparison. However, we note that such pairwise comparisons warrant multiple testing corrections before any inferential claims can be made. A clear negative association was also observed between female sex and long COVID QALDs, with an estimated NDE of -6.79 (-12.8, -0.72), NIE of -3.05 (-5.89, -0.22) and proportion non-mediated of 0.77 (0.46,1). Effect estimates for the NDE, NIE, and proportion non-mediated were nearly equivalent, but with reduced stability, in the sensitivity analysis applying population weights (Supplementary Figure 10).

UK

The UK cohort was skewed towards older adults (mean (SD): 59.0 (12.6) years). Nearly all participants self-identified as white (96%), representing an even greater deviation from the underlying race/ethnicity distribution in the UK.⁵⁰ Again, for this reason, we exclude race/ethnicity from our “exposure” variables/predictors for long COVID QALDs. Most were male (59%), which is not reflective of the underlying UK population, as females are the slight majority (51%).⁵¹

The most commonly reported comorbidity was hypertension (36%). The median (IQR) value for estimated long COVID QALDs was 295 (233, 342), the lowest of all three cohorts.

Employment status was markedly skewed towards full-employment (51%), retirement (30%), part-time employment (10%) and unemployment (7%). The least represented categories were students (0.6%) and the furloughed (0.5%). Estimated long COVID QALDs were greatest among participants who reported being furloughed, students or full-time employees and lowest among those in the unemployed and retired categories.

Employment status was the leading predictor for long COVID QALDs in the UK, followed by psychological disorder, age, employment status category, chronic neurological disorder, and rheumatological disorder (Figure 3a). Sex followed in the rankings, which, along with the acute COVID-19 severity indicator, fell among the top ten predictors (Figure 3a). RF #2 also revealed the prime predictive role of employment status and sex as a group, with the PCs of the socio-demographic variables in aggregate dominating other grouped variables, closely followed by the PCs of the group of mental health and neurological disorders (Figure 3c). Age alone ranked highly, even in comparison to grouped factors (Figure 3c). Findings from RF #3 aligned well with these results with age, chronic neurological disorder, employment status indicators, and psychological disorder the most commonly selected variables, followed by rheumatological disorder, across key clusters (Figure 3b). Finally, the population correction sensitivity analysis reported very similar rankings of variables for both RF #1 and RF #3 implementations (Supplementary Figure S13a and c). The same was true for RF #2, where all top ten most selected variables were the same in the main and sensitivity analysis, with the exception of obesity (Supplementary Figure S13b).

• Download figure
• Open in new tab

Figure 3a.

Estimated variable importance measures, i.e. % increase in mean squared error or MSE, from individual random forest implementation (RF #1) for the UK.

• Download figure
• Open in new tab

Figure 3b.

Number of times (frequency) each variable appears in clusters selected for each CoV-VSURF run (RF #3) for the UK.

• Download figure
• Open in new tab

Figure 3c.

Estimated variable importance measures, i.e. % increase in mean squared error or MSE, from pre-grouped random forest implementation (RF #2) for the UK. Rows indicate cluster names (a full list of variables belonging to each cluster can be found in Supplementary Table S3) and corresponding principal components, if the cluster consists of multiple variables. PC1 denotes principal component 1 and PC2 denotes principal component 2.

QALDs were negatively associated with all major predictive comorbidities, particularly for psychological disorder and chronic neurological disorder, the acute COVID-19 severity indicator, and acute treatment with antivirals (Supplementary Figure S7). In contrast, QALDs were modestly positively associated with age (Supplementary Figure S7).

Increased income/job stability, as proxied by employment status, was consistently associated with increased long COVID QALDs, irrespective of the binary designation used. We found that self-reported full-time employment compared to any other employment status category (excluding retirement) was associated with, on average, 31.7 (14.2, 49.3) higher expected long COVID QALDs, via the NDE, and 4.90 (-0.065, 9.86) higher expected long COVID QALDs, via the NIE. These results yield a proportion non-mediated of 0.86 (0.73, 0.996), i.e., 86% of the effects of full-time employment versus all other employment status categories on long COVID QALDs were unexplained by the mediators in question. We obtained an even stronger relationship between self-reported full-time employment versus unemployment status, with, on average, 79.5 (50.0, 109.0) increased long COVID QALDs among the full-time employed relative to the unemployed (NDE) and 9.50 (1.04, 18.0) increased long COVID QALDs among the full-time employed versus unemployed (NIE). The proportion non-mediated is thus slightly higher at 0.91 (0.83, 0.98), suggesting that around 91% of the effect of full-time employment versus unemployment on long COVID QALDs cannot be attributed as operating through the considered mediators.

Female sex was associated with lower expected long COVID QALDs, with an NDE and NIE of -24.2 (-37.8, -10.7) and -9.61 (-16.3, -2.95), respectively. The corresponding proportion non-mediated was the lowest observed among all contrasts: only 73% (52%, 94%) of the effect of female sex on long COVID QALDs were explained by the clinical intermediates. As in Norway, the population correction sensitivity analysis reproduced these findings, but with greater uncertainty (Supplementary Figure S11).

Discussion

In this study, we provided a quantitative assessment of the extent to which social factors, compared to commonly highlighted clinical conditions, contribute to varying experiences with long COVID.

The data provided compelling evidence for specific categories of pre-existing comorbidities, namely neurological, psychological, and rheumatological, being major predictors of long COVID-associated quality of life. Age also consistently ranked highly as a predictor across the multiple middle and high income countries represented in Europe, Asia, Africa, and South America, although the direction of this association differed by context. However, of note, our evaluated bivariate associations reflect the unadjusted role of a given factor of interest on long COVID QALDs and cannot be used to draw any meaningful conclusions about causality. For instance, the negative association observed between antiviral treatment and long COVID QALDs in the UK is likely to be an artifact of confounding by acute disease severity.

Importantly, we observed that socio-demographic variables educational attainment or employment status and sex at birth were generally as or more predictive of long COVID QALDs than the aforementioned comorbidities. Our mediation analyses further suggest that not only are indicators of societal disadvantage highly predictive of lower long COVID QALDs, but also that the impacts of these variables on long COVID QALDs could only be partially explained by key long COVID-predicting comorbidities. This general finding, i.e., that disparities are not solely attributable to underlying differences in comorbidity rates across various demographic groups, has been validated in other studies conducted over the course of the pandemic.^20,25,52

Our study benefited from the use of a sizable and diverse international cohort with data on long-term quality of life following SARS-CoV-2 infection, providing more exhaustive information on post-acute COVID-19 experiences beyond simply whether patients experience long COVID symptoms. Given the sufficiently large sample sizes for our selected cohorts, we were able to apply adaptable machine learning tools, including recent developments in causal machine learning.^40,46 The variable selection methods we used avoid strict and potentially biased modeling assumptions required for commonly used parametric regressions and accommodate inherent variable groupings. Additionally, the causal mediation approaches we applied extend traditional causal inference frameworks⁵³ through the integration of flexible, but simultaneously relatively precise, model-free regression algorithms ^40,46, providing a promising alternative, as aligned with previously outlined frameworks⁵⁴, to parametric, model-based approaches.

There are several important limitations of our analysis. First, as our outcome of interest was post-acute COVID-19 long-term quality of life, we were unable to examine the varying roles of the different social and clinical factors on changes in quality of life. For future work, it is thus imperative to collect information on quality of life measures at all stages of illness, not simply ex post facto, to avoid issues with recall by positioning readily implementable epidemic study protocols at the outset of an outbreak. Relatedly, we did not have information on subject-specific duration of long COVID, and instead assumed uniform duration for all participants consistent with examples in the literature on the duration of long COVID.^33,34 Additionally, while participant recruitment was extensive for each of the cohorts considered in this analysis, filtering each dataset to subjects with available demographic, comorbidity, and quality of life data led to sharply reduced sample sizes, especially for MICs. This led to reduced power and the need to combine multiple, highly varying, MIC countries into one cohort, despite potentially profoundly different relationships between each of the social factors and long COVID QALDs in each country. While we found evidence for generally comparable directions of association for the country-specific analyses compared to those of the full combined cohort, which would indicate that there may not be meaningful effect modification by country, some measures did vary in sign and magnitude, notably lower estimated proportions non-mediated, and all results were generally less stable. Consequently, country-specific analyses where sample sizes are sufficient are necessary for disentangling whether these represent real differences by country within the combined cohort, or possible artifacts of small samples. Finally, variables like vaccination status and antiviral treatment were sparsely recorded in some cohorts mainly due to early implementation of the follow-up studies before vaccination was available.

In the Norway cohort, we observed that vaccination status was a very weak predictor of long COVID QALDs, but we note that there could be possible type-specific, dose-specific, or duration-specific effects that we were unable to evaluate here given that the relevant data was nonexistent or insufficient. Thus, the prioritization of both data unification procedures that enable the synchronous collection of the range of variables considered, and the collection of data in less resource-rich/high-attention settings, is crucial.

Finally, it can be concluded that survivorship bias might affect our results, as only subjects who completed a follow-up survey at any follow-up interval can have their quality of life measures recorded. Those lost to follow up due to death from early myocardial infarction, vascular strokes, etc. are also likely to have a reduced quality of life prior to this event.⁵⁵ However, no participants in the Norway and UK cohorts died at any point during follow up and only 1.3% of subjects in the combined MIC cohort with follow-up information died by the conclusion of the study.

There were additional limitations related to the type of information we had at our disposal. First, the marked underrepresentation of race/ethnic minorities in the UK and Norway cohorts, and the lack of race/ethnicity data or relevance of this construct in some contexts, such as the Russian and Indian cohorts, prevented us from investigating the impacts of race/ethnicity. It is crucial that future studies seek to recruit more balanced cohorts in this respect, and, where applicable, address constraints with documenting such information. Importantly, in order for a meaningful analysis of race/ethnic risk factors to be conducted in settings where race/ethnicity is a relevant determinant of health, controlled definitions of minority groups must be tailored to and assessed independently in each context. Additionally, while we did have data available for socio-economic proxies, we were limited to specific variables that may not fully reflect participants’ levels of socio-economic deprivation and had to apply somewhat arbitrary aggregations of variables due to data and methodological constraints. For future work, it would be useful to both emphasize the collection of more proximate indicators of socio-economic status and consider extensions of the mediation approaches that allow for categorical exposures. Finally, we were limited to data on sex at birth, which does not capture important gender-based disparities that exist beyond this binary.⁵⁶

Despite these challenges, the central aim of this analysis was to outline a robust statistical and causal-analytic framework, applied to example case studies, for highlighting the contribution of social disparities to chronic ill-health. Our framework can be used both as a standardized comparison of a collective of diverse variables, grouping related factors when necessary, as determinants of worsened post-acute COVID-19 health and well-being and a way to distill the unique impacts of any social variable of interest on these outcomes. Our data highlights the multifactorial relationship between pre-existing risk factors and socio-economic factors and recovery from SARS-CoV-2 infection. As such, we demonstrate that accounting for social vulnerabilities when evaluating determinants of post-acute COVID-19 trajectories is essential and that studies focusing solely on clinical targets may not be sufficient. Consequently, approaches aiming to alleviate social disparities in long COVID recovery by exclusively targeting comorbidity prevention and management are likely to be restricted in their impact. Conversely, transformational societal interventions, that can address disease exposures and access to care, educational, and employment, and other social determinants of health, have the opportunity to lead to potentially more comprehensive benefits and improve overall well-being in marginalized communities.

Funding Statement

TFM acknowledges support from NIH Training Grant 2T32AI007535. LFR was funded by Universidad de La Sabana (MED-309-2021). MS has been funded (in part) by contracts 200-2016-91779 and cooperative agreement CDC-RFA-FT-23-0069 with the Centers for Disease Control and Prevention (CDC). The findings, conclusions, and views expressed are those of the author(s) and do not necessarily represent the official position of the CDC. MS was also partially supported by the National Institute of General Medical Sciences of the National Institutes of Health under award number R01GM130668.This work was made possible by the UK Foreign, Commonwealth and Development Office and Wellcome [215091/Z/18/Z, 222410/Z/21/Z, 225288/Z/22/Z and 220757/Z/20/Z]; the Bill & Melinda Gates Foundation [OPP1209135]; the philanthropic support of the donors to the University of Oxford’s COVID-19 Research Response Fund (0009109); grants from the National Institute for Health Research (NIHR; award CO-CIN-01/DH_/Department of Health/United Kingdom), the Medical Research Council (MRC; grant MC_PC_19059), and by the NIHR Health Protection Research Unit (HPRU) in Emerging and Zoonotic Infections at University of Liverpool in partnership with Public Health England (PHE), (award 200907), NIHR HPRU in Respiratory Infections at Imperial College London with PHE (award 200927), Liverpool Experimental Cancer Medicine Centre (grant C18616/A25153), NIHR Biomedical Research Centre at Imperial College London (award ISBRC-1215-20013), and NIHR Clinical Research Network providing infrastructure support; the Comprehensive Local Research Networks (CLRNs) of which PJMO is an NIHR Senior Investigator (NIHR201385); Cambridge NIHR Biomedical Research Centre (award NIHR203312); the Research Council of Norway grant no 312780, and a philanthropic donation from Vivaldi Invest A/S owned by Jon Stephenson von Tetzchner to the Norwegian SARS-CoV-2 study; the South Eastern Norway Health Authority and the Research Council of Norway.

Competing interests

MS has received institutional research funds from the Johnson and Johnson foundation and from Janssen global public health. MS also received institutional research funding from Pfizer.

ISARIC Clinical Characterization Group

Beatrice Alex, Eyvind W. Axelsen, Benjamin Bach, John Kenneth Baillie, Wendy S. Barclay, Joaquín Baruch, Husna Begum, Lucille Blumberg, Debby Bogaert, Fernando Augusto Bozza, Sonja Hjellegjerde Brunvoll, Polina Bugaeva, Aidan Burrell, Denis Butnaru, Roar Bævre-Jensen, Gail Carson, Meera Chand, Barbara Wanjiru Citarella, Sara Clohisey, Marie Connor, Graham S. Cooke, Andrew Dagens, John Arne Dahl, Jo Dalton, Ana da Silva Filipe, Emmanuelle Denis, Thushan de Silva, Pathik Dhangar, Annemarie B. Docherty, Christl A. Donnelly, Thomas Drake, Murray Dryden, Susanne Dudman, Jake Dunning, Anne Margarita Dyrhol-Riise, Linn Margrete Eggesbø, Merete Ellingjord-Dale, Cameron J. Fairfield, Tom Fletcher, Victor Fomin, Robert A. Fowler, Christophe Fraser, Linda Gail Skeie, Carrol Gamble, Michelle Girvan, Petr Glybochko, Christopher A. Green, William Greenhalf, Fiona Griffiths, Matthew Hall, Sophie Halpin, Bato Hammarström, Hayley Hardwick, Ewen M. Harrison, Janet Harrison, Lars Heggelund, Ross Hendry, Rupert Higgins, Antonia Ho, Jan Cato Holter, Peter Horby, Samreen Ijaz, Mette Stausland Istre, Clare Jackson, Waasila Jassat, Synne Jenum, Silje Bakken Jørgensen, Karl Trygve Kalleberg, Christiana Kartsonaki, Seán Keating, Sadie Kelly, Kalynn Kennon, Saye Khoo, Beathe Kiland Granerud, Anders Benjamin Kildal, Eyrun Floerecke Kjetland, Paul Klenerman, Gry Kloumann Bekken, Stephen R Knight, Andy Law, Jennifer Lee, Gary Leeming, Wei Shen Lim, Andreas Lind, Miles Lunn, Laura Marsh, John Marshall, Colin McArthur, Sarah E. McDonald, Kenneth A. McLean, Alexander J. Mentzer, Laura Merson, Alison M. Meynert, Sarah Moore, Shona C. Moore, Caroline Mudara, Daniel Munblit, Srinivas Murthy, Fredrik Müller, Karl Erik Müller, Nikita Nekliudov, Alistair D Nichol, Mahdad Noursadeghi, Anders Benteson Nygaard, Piero L. Olliaro, Wilna Oosthuyzen, Peter Openshaw, Massimo Palmarini, Carlo Palmieri, Prasan Kumar Panda, Rachael Parke, William A. Paxton, Frank Olav Pettersen, Riinu Pius, Georgios Pollakis, Mark G. Pritchard, Else Quist-Paulsen, Dag Henrik Reikvam, David L. Robertson, Amanda Rojek, Clark D. Russell, Aleksander Rygh Holten, Vanessa Sancho-Shimizu, Egle Saviciute, Janet T. Scott, Malcolm G. Semple, Catherine A. Shaw, Victoria Shaw, Louise Sigfrid, Mahendra Singh, Vegard Skogen, Sue Smith, Lene Bergendal Solberg, Tom Solomon, Shiranee Sriskandan, Trude Steinsvik, Birgitte Stiksrud, David Stuart, Charlotte Summers, Andrey Svistunov, Arne Søraas, Emma C. Thomson, Mathew Thorpe, Ryan S. Thwaites, Peter S Timashev, Kristian Tonby, Lance C.W. Turtle, Anders Tveita, Timothy M. Uyeki, Steve Webb, Jia Wei, Murray Wham, Maria Zambon.