Skip to main content
Download PDF

Ethnic disparities in COVID-19 outcomes: a multinational cohort study of 20 million individuals from England and Canada

BMC Public Health volume 23, Article number: 399 (2023) Cite this article

Abstract

Background

Heterogeneous studies have demonstrated ethnic inequalities in the risk of SARS-CoV-2 infection and adverse COVID-19 outcomes. This study evaluates the association between ethnicity and COVID-19 outcomes in two large population-based cohorts from England and Canada and investigates potential explanatory factors for ethnic patterning of severe outcomes.

Methods

We identified adults aged 18 to 99 years in the QResearch primary care (England) and Ontario (Canada) healthcare administrative population-based datasets (start of follow-up: 24th and 25th Jan 2020 in England and Canada, respectively; end of follow-up: 31st Oct and 30th Sept 2020, respectively). We harmonised the definitions and the design of two cohorts to investigate associations between ethnicity and COVID-19-related death, hospitalisation, and intensive care (ICU) admission, adjusted for confounders, and combined the estimates obtained from survival analyses. We calculated the ‘percentage of excess risk mediated’ by these risk factors in the QResearch cohort.

Results

There were 9.83 million adults in the QResearch cohort (11,597 deaths; 21,917 hospitalisations; 2932 ICU admissions) and 10.27 million adults in the Ontario cohort (951 deaths; 5132 hospitalisations; 1191 ICU admissions). Compared to the general population, pooled random-effects estimates showed that South Asian ethnicity was associated with an increased risk of COVID-19 death (hazard ratio: 1.63, 95% CI: 1.09-2.44), hospitalisation (1.53; 1.32-1.76), and ICU admission (1.67; 1.23-2.28). Associations with ethnic groups were consistent across levels of deprivation. In QResearch, sociodemographic, lifestyle, and clinical factors accounted for 42.9% (South Asian) and 39.4% (Black) of the excess risk of COVID-19 death.

Conclusion

International population-level analyses demonstrate clear ethnic inequalities in COVID-19 risks. Policymakers should be cognisant of the increased risks in some ethnic populations and design equitable health policy as the pandemic continues.

Peer Review reports

Introduction

There have been almost 450 million SARS-CoV-2 infections and 6 million deaths (as of March 2022) worldwide since the novel coronavirus disease 2019 (COVID-19) pandemic emerged [1]. Several studies have demonstrated ethnic inequalities in the risk of infection and adverse outcomes, which has generated substantial concern [2,3,4,5]. In the United Kingdom (UK), compared to White individuals, men from all backgrounds other than Chinese, and women from any ethnic group other than Bangladeshi, Chinese or mixed ethnicity, had an increased risk of COVID-19 mortality when accounting for differences in demographics, socioeconomic status, and comorbidities [6]. Notably, in the UK Black African men and women were over 2 times as likely to die from COVID-19 than those of White ethnicity [6]. Other large-scale epidemiological analyses from the UK demonstrated that those from South Asian, Black, and ‘Mixed’ ethnic groups had increased rates of COVID-19 death compared to the White group [2]. However, the evidence from other health systems such, as the United States, is conflicting [7], and data on COVID-19 cases and mortality in Canada by ethnicity are more limited [8, 9].

The mechanisms driving the inequalities are unclear but have been posited to be related to a complex and interrelated patterning of multiple factors, including medical factors – such as comorbidities and medication use – as well as social determinants, including cultural, behavioural, and occupational factors, and structural inequalities [10]. The presence of comorbidities has been associated with both a higher risk of SARS-CoV-2 infection and worse outcomes in individuals with COVID-19 [11, 12], while medications (i.e., certain glucose-lowering [13] or immune-modifying drugs [14]) have been linked to either an increased or reduced risk of COVID-19 outcomes. Among the social factors, poor living and working conditions, low income, health literacy, poverty, or exposure to air pollution have all been associated with COVID-19 infectivity and mortality [15,16,17]. Robustly ascertaining the comparative contributory effects has been difficult to untangle, and one study [18] has sought to quantify potential mediators, rather than reporting ‘overall’ effects [2, 4, 19,20,21].

Establishing a nuanced understanding of ethnic inequalities in COVID-19-related outcomes is needed to reduce the burden of COVID-19 and may permit rapid public health interventions should modifiable factors be identified. Here, we carried out two observational studies with harmonised designs to reduce the bias due to heterogeneous definitions of exposures, outcomes, and confounders, in the UK (England) and Ontario to quantify the associations between ethnicity and COVID-19 severity and explore potential modifiable and non-modifiable explanatory factors. We then sought to synthesise these cohort-level estimates using a meta-analysis.

Methods

Data sources and study population

QResearch

QResearch database (version 45) comprises individuals registered across 1321 general practices covering 18% of the English population with linkages of primary care data to hospitalisation, intensive care (ICU) admission, and mortality data. For this study, we included 9,828,099 adults aged 18 to 99 years contributing to the QResearch database with at least 12 months of continuous prior registration. The study period ran from the date of the first confirmed SARS-CoV-2 infection in England (24th January 2020, start of follow-up) until 31st October 2020 (end of follow-up), the occurrence of outcome, or death, whichever occurred earlier.

Ontario

The second data source is the population-level healthcare administrative data in Ontario, Canada’s most populous and most ethnically diverse province. These data include the entire population of Ontario (currently 14.5 million, representing nearly 40% of the Canadian population) and are linked to sociodemographic information, hospital, and ICU admissions; in this investigation, 10,273,496 people aged over 18 years were included. The study period ran from the 25th January 2020 (start of follow-up) to 30th September 2020 (end of follow-up), the occurrence of outcome, or death, whichever occurred earlier.

Ethnicity and COVID-19 outcome: pooled analysis

In the first analysis, we explored the association between self-reported ethnicity and COVID-19 related death, hospitalisation, and ICU admission: these outcomes were slightly different in the QResearch and Ontario cohorts as based on country-specific definitions. For QResearch, outcomes included: (a) COVID-19 death, defined as either confirmed or suspected COVID-19 on death certificate, or a death from any cause with a confirmed positive SARS-CoV-2 test in the immediately preceding 28 days; (b) Hospitalisation due to COVID-19, defined as an admission with confirmed or suspected COVID-19 (as per ICD-10 codes U07.1 and U07.2), or new hospitalisation with a positive SARS-CoV-2 test in the immediately preceding 14 days; (c) ICU admission due to COVID-19, defined as admission to ICU with confirmed or suspected SARS-CoV-2 test in the preceding 28 days. In the Ontario database, outcomes were defined as: (a) COVID-19 death, defined as any death with a confirmed positive SARS-CoV-2 test in the immediately preceding 28 days; (b) Hospitalisation due to COVID-19, defined as an admission with confirmed or suspected COVID-19 (as per ICD-10 codes U07.1 and U07.2), or with a positive SARS-CoV-2 test between 28 days prior to and 14 days after the admission date; (c) ICU admission due to COVID-19, defined as a hospital admission that included ICU stay with confirmed or suspected COVID-19 (as per ICD-10 codes), or with a positive SARS-CoV-2 test between 28 days prior to and 14 days after the admission date.

We utilised a 3-level ethnicity classification comprised South Asian (Indian, Bangladeshi, Pakistani), Chinese, and ‘General Population’ (all other ethnic groups, 87.5% White in this cohort), based on the UK Office for National Statistics Census ethnic classification. In the Ontario linked healthcare administrative database, ethnicity was ascertained based on surnames, using lists that have been previously validated in this population to identify the two largest ethnic groups in Canada: South Asian and Chinese [22]. The positive predictive values for this approach to identifying ethnicity, when compared to self-reported ethnicity, are high: 89.3% for South Asians and 91.9% for Chinese; specificity 99.7% for both. People whose surnames were not on either list were labelled as ‘General population’ (all other ethnic groups, approximately 80% White). ‘General Population’ was used as the reference category for analyses.

The analyses were adjusted for demographic, clinical, and lifestyle factors (Supplementary Material; Tables S1 and S2); estimates of the associations between ethnicity and each of the three outcomes obtained in the QResearch and Ontario cohorts were combined in a two-stage random-effects meta-analysis. Further details on the definitions of the population and confounders are reported in the Supplementary Material.

Percentage of excess risk mediated by risk factors

The contribution of possible ‘risk factor’ classes to the increased relative risks in different ethnic groups was quantified in the QResearch data as the ‘percentage of excess risk mediated’ (PERM) [23]. By evaluating the change in the magnitude of the exposure-outcome association in models with different confounders, this analysis helps clarify the extent to which a confounder (or a set of confounders) accounts for the association between ethnicity and COVID-19 outcome. For the PERM analyses, we defined 5-level ethnic groups as Mixed ethnicity, South Asian, Black, and ‘Other’ ethnic groups; hazard ratios (HRs), relative to White, were estimated separately for each of the three outcomes and the following set of confounders: ‘minimally adjusted’ model (age, sex, and region); household and social factors; comorbidities; lifestyle factors (including BMI); and ‘maximally adjusted’ model.

Statistical analyses

Country-specific baseline socio-demographic and clinical characteristics were summarised using descriptive statistics by COVID-19 related hospitalisation, ICU admission, and mortality. In QResearch, survival analyses to evaluate the adjusted association of 3-level ethnicity with COVID-19 outcomes, accounting for clustering of practices (robust standard error), were performed with the Royston-Parmar model [24, 25]. We performed multiple imputation to replace missing values for ethnicity (20.2% missing), deprivation (0.6%), BMI (16.7%), and smoking status (5.2%) using chained equations under the missing at random assumption. These variables were modelled following a multinomial logistic model for ethnicity, ordinal logistic model for smoking and alcohol, and truncated regression for BMI. Results from five imputations were pooled using Rubin’s rules [26]. Complete case analyses and time-varying associations were performed as sensitivity analyses. Cox survival regressions were conducted to evaluate the adjusted associations of ethnicity with COVID-19 outcomes in the Ontario database; as the only missing data were a small number of people missing deprivation data (0.2%), we used complete cases regressions. The proportional hazards assumption was checked in both survival analyses by plotting log-log plot. HRs from maximally adjusted models were used as the common measure of association across QResearch and Ontario cohorts and combined with the DerSimonian-Laird random-effects method in a two-stage meta-analysis; heterogeneity was assessed with I 2.

In QResearch, we applied the following formula to estimate PERM using the HR across imputed datasets was:

$$\textrm{PERM}=100\frac{\ \left[\textrm{HR}\left(\textrm{age},\textrm{sex},\textrm{region}\right)\kern0.5em -\kern0.5em \textrm{HR}\left(\textrm{age},\textrm{sex},\textrm{region}+\textrm{risk}\ \textrm{factor}\ \textrm{group}\right)\right]}{\left[\textrm{HR}\left(\textrm{age},\textrm{sex},\textrm{region}\right)-1\right]\kern0.5em }$$

The PERM was also calculated for ‘maximal adjustment’ in each of the non-White ethnic groups to assess the extent to which inequalities were potentially attributable to the large set of measured adjustment factors.

All p-values are two sided and nominal statistical significance was considered at p < 0.05. We used Stata v.17 for the QResearch statistical analyses and SAS v.9.4 for the Ontario analyses. We followed current guidance for conducting and reporting observational studies using routinely collected health data (RECORD checklist in the Supplementary Material).

Patient and public involvement reporting

Two public representatives advised on interest and appropriateness of the research questions, were involved in writing the protocol for the wider study, and input on lay-summaries describing the planned study.

Results

Study populations

In the QResearch cohort, there were 9,828,099 individuals; during follow-up, 11,597 COVID-19 deaths, 21,917 hospitalisations and 2932 ICU admissions occurred; in the Ontario cohort, corresponding figures were 10,273,496 individuals, 951 COVID-19 deaths, 5132 hospitalisations, and 1191 ICU admissions (Table 1). Ethnicity data and their classifications are summarised in Table 1 and characteristics stratified by ethnicity are provided in Tables S1 -S2.

Table 1 Baseline sociodemographic, clinical characteristics and outcomes in the QResearch and Ontario cohorts
Full size table

Cohort studies and meta-analyses

In QResearch, South Asian ethnicity was associated with increased rates of COVID-19 mortality (HR: 1.35; 95% CI: 1.20, 1.51; Fig. 1 and S1), hospitalisation (1.63; 1.51, 1.75; Fig. 1 and S2), and ICU admission (1.93; 1.67, 2.25; Fig. 1 and S3) compared to the general population group; corresponding estimates in Ontario were 2.04 (1.56, 2.68) for mortality, 1.41 (1.24, 1.59) for hospitalisation, and 1.41 (1.10, 1.79) for ICU admission. In the same maximally adjusted models, in QResearch there was no evidence of increased rates of COVID-19 mortality (HR: 1.12; 0.75, 1.66), hospitalisation (0.86; 0.67, 1.11), or ICU admission (1.20; 0.68, 2.11) in Chinese ethnic group compared to the general population group, whilst in Ontario the HRs were 0.92 (0.67, 1.25) for mortality, 0.79 (0.69, 0.91) for hospitalisation, and 1.29 (1.02, 1.63) for ICU admission. For all three outcomes, the direction of associations was similar for most of the confounders available in both the QResearch and Ontario cohorts, indicating an increased risk associated with the presence of medical conditions and a progressively higher risk in older people and larger households (Fig. S1-S3). In the QResearch cohort, complete case estimations were largely similar to those of the main analyses using multiple imputation (Fig. S4); time-varying associations by ethnic groups are presented in Fig. S5.

Fig. 1
figure 1

Cohort-level meta-analysis of individual participant data from QResearch and Ontario. Estimates and number of events and participants are shown following multiple imputation in QResearch cohort and for complete-case analysis in Ontario cohort. The reference ethnic group is “general population”, including: (1) people not South Asian and Chinese in Ontario (approximately 80% White); (2) White, Other Asian, Black African, Black Caribbean, and Other in QResearch

Full size image

Combining estimates for South Asian ethnicity across QResearch and Ontario cohorts resulted in a random-effects HR of 1.63 (1.09, 2.44) for COVID-19 related mortality, with considerable heterogeneity between the two estimates (I 2 86.9%; Fig. 1). Corresponding estimates for hospitalisation and ICU admission were 1.53 (1.32, 1.76) and 1.67 (1.23, 2.28), with considerable heterogeneity: I 2 75.4% and I2 74.9%, respectively. The pooled random-effects HRs comparing Chinese ethnicity to the general population were 0.99 (0.77, 1.26) for mortality, 0.81 (0.72, 0.91) for hospitalisation, and 1.28 (1.03, 1.58) for ICU admission; there was no evidence of heterogeneity for all three outcomes (I2 0%; Fig. 1). There was no clear trend in the mortality, hospitalisation, or ICU admission HRs comparing ethnic groups across levels of deprivation (Fig. S6 ).

Percentage of excess risk mediated by risk factor classes (QResearch)

The percentage of excess risk mediated by separate groups of potential attributable factors across the entirety of follow-up in QResearch is reported in Table 2. We estimated that approximately 20-30% of the excess risk of COVID-19-related hospitalisation in non-White ethnic groups may be mediated by household size/status and deprivation; and that differences in comorbidity prevalence may mediate up to approximately 20% of excess risk (in South Asian). For COVID-19-related ICU admission, adjustment for comorbidities accounted for up to approximately 30% of the excess risk, whereas maximal adjustment accounted for up to approximately 40% of the excess risk (in Black ethnic group). Differences in smoking habits and BMI did not appear to mediate any degree of excess risk of COVID-19-related death in any non-White ethnic group. Maximal adjustment accounted for 42.9% (South Asian) and 39.4% (Black) of the excess risks of death. Therefore, the majority of excess risk in non-White groups may not be accounted for the range of sociodemographic, lifestyle, and comorbidity factors considered in this analysis.

Table 2 Percentage of excess risk mediated in COVID-19-related outcomes in ethnic minority groups (QResearch cohort)
Full size table

Discussion

In this international study of population-level healthcare databases covering over 20 million individuals, we showed that adults of South Asian background had a 63% increased risk of COVID-19 mortality, 53% increased risk of COVID-19-related hospital admission, and 67% increased risk of ICU admission overall compared to the general population in England and Ontario. This compares to 28% of increased risk of ICU admission in Chinese, with no evidence of increased mortality and hospitalisation risks. In England, sociodemographic, lifestyle, and clinical factors accounted for approximately 40% of excess risks of COVID-19 death.

Our results are consistent with other UK population-level analyses derived from data using combinations of different IT systems, which also reported similar estimates of risk in non-White ethnic groups [2]. In this respect, it is important to note that the risks of COVID-19 outcomes estimated in QResearch across ethnic groups, and combined with the results from Ontario, should be considered in view of some variations in the magnitude of associations between ethnicity and COVID-19 outcomes both between waves and within the same wave; more importantly, the public health implications of these variations are primarily determined by the country- and region-specific change in the absolute risk of each outcome over time [27].

Whilst it is increasingly established in the literature that non-white ethnicity is associated with increased risk of severe COVID-19 outcomes, the degree to which modifiable and other factors may contribute to this risk in different ethnic groups is poorly understood. Some ethnic communities may be disadvantaged as living in poorer socioeconomic environments where the risk of infection and worse outcomes is higher, including overcrowded multigenerational houses or occupations with a high degree of public contact [2, 18]; at the same time, biological factors have been suggested to play a role as well, such as an unfavourable metabolic-inflammatory milieu (i.e., obesity, multimorbidity) [11, 20, 28]. In our investigation, rather than reporting summary effect estimates after full or serial adjustment, our approach in the QResearch also included assessment of relative contribution of potential attributable factors and suggests that there may be heterogeneity in the mechanistic underpinnings the increased risks in different ethnic groups. Our study found that the sociodemographic, lifestyle, and clinical factors considered in this investigation accounted for approximately 40% of excess risks of COVID-19 death. Hence, further research should investigate whether other factors, not captured in our data, may explain the proportion of excess risks in some ethnic groups and possible causal pathways in the COVID-19 syndemic [29]. It is possible, in fact, that ethnic differences are at least in part the epiphenomenon of a complex network of other risk factors associated with a higher risk of COVID-19 outcomes, including overcrowding and occupation [30].

Our study analysed in greater detail the differential effects of deprivation within ethnic groups, as well as the relative contributions of different factors to the increased risks in non-White groups, given the suggested interplay between ethnicity and deprivation on the risk of COVID-19 outcomes [31]. Our results also expand and clarify the evidence base regarding ethnic inequalities in COVID-19 outcomes in several ways. First, in contrast to evidence generated using data only from those attending hospitals or registered with providers within fragmented healthcare systems investigating the role of sociodemographic and clinical characteristics on the risk of outcomes across ethnic groups [27, 32], our population-level approach examined the relevant risk trajectories and avoided conditioning on positive tests or other intermediates [33]. Second, much of the available evidence about ethnicity and COVID-19 related outcomes is highly heterogeneous in terms of study designs, population, definitions of outcomes/exposures, confounders adjusted for (if any), and settings (geographical and healthcare system). This negatively affected individual study interpretation but also limited the cohesive synthesis of evidence via meta-analytical approaches due to significant within- and between-study heterogeneity [4, 34]. We explicitly sought to harmonise analytical approaches to facilitate pooling of robust estimates from multiple geographical units, namely different nations (England and Canada). Other key strengths of our study include the use of two large, population-level and representative healthcare databases without selection bias, which possess individual-level linkages across the healthcare network enabling accurate ascertainment of exposures, confounders, and outcomes. Our flexible harmonisation of definitions and analytical approaches facilitated cohort-level meta-analysis of results from both main study databases; we also used the Royston-Parmar survival model which allowed us to explore whether the association between ethnicity and COVID-19 related outcomes changed across the first and second wave. Lastly, we investigated the possible mediation role of some factors in explaining the increased risk observed across ethnic groups in UK. In this respect, it should be noted that different methods exist to investigate mediation (including the possibility to account for intermediate confounding) [35]; furthermore, while the difference between a confounder and mediator is well-known, the same factors may be considered mediators in some context and confounders in others, or even in the same context by different investigators [36], further highlighting the complex interactions among multiple factors in determining the health status. Moreover, some potential mediators have not been included in our analyses (i.e., education, employment status, income). As such, our PERM results should be considered explorative and no definitive causal inference can be derived from them: it is plausible that the comparative causal role of these factors would be different in heterogeneous healthcare systems and societies. Our study has also some limitations, including the inability to further disaggregate ethnicity into more granular groups in Ontario; lack of recorded other information that may be relevant to disease risks (such as occupation, which is relevant to SARS-CoV-2 exposure, and detailed household composition) [37]; the risk of residual confounding, which affects every observational analysis and hampers a conclusive causal interpretation; missing data, which were addressed assuming a missing at random mechanism, yet previous evidence would indicate that ethnicity could be missing not at random: [38] however, the complete-case analysis for the latter scenario [39] resulted in estimates virtually identical to those obtained using multiple imputation; and the potential variations in the ascertainment of COVID-19 infections over time, between countries, and among ethnic groups [40]. Furthermore, the contribution of potential attributable factors was explored only in the QResearch cohort as several of these factors were not available in the Ontario administrative data.

Evidence from large-scale cohort studies in England and Canada and from meta-analyses provide robust evidence of ethnic inequalities in COVID-19 outcomes. Not only do these persist despite accounting for potential sociodemographic and clinical confounders but the risks in individual ethnic groups have varied during the pandemic. The currently unexplainable proportion of excess risks in non-White groups requires careful consideration of economic, healthcare system, and other factors to guide public health strategy to protect everyone as the pandemic progresses globally.

Availability of data and materials

Public access to the databases is closed. To guarantee the confidentiality of personal and health information, only the authors have had access to the data during the study in accordance with the relevant license agreements.

Ontario: The dataset from this study is held securely in coded form at ICES. While legal data sharing agreements between ICES and data providers (e.g., healthcare organizations and government) prohibit ICES from making the dataset publicly available, access may be granted to those who meet pre-specified criteria for confidential access, available at https://www.ices.on.ca/DAS. QResearch: Access to QResearch data is according to the information on the QResearch website (www.qresearch.org).

References

  1. Coronavirus cases: Reported Cases and Deaths by Country, Territory, or Conveyance. https://www.worldometers.info/coronavirus/ (Accessed 3 Aug 2021).

  2. Mathur R, Rentsch CT, Morton CE, et al. Ethnic differences in SARS-CoV-2 infection and COVID-19-related hospitalisation, intensive care unit admission, and death in 17 million adults in England: an observational cohort study using the OpenSAFELY platform. Lancet. 2021;397(10286):1711–24.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Rimmer A. Covid-19: disproportionate impact on ethnic minority healthcare workers will be explored by government. BMJ. 2020;369:m1562.

    Article  PubMed  Google Scholar 

  4. Sze S, Pan D, Nevill CR, et al. Ethnicity and clinical outcomes in COVID-19: a systematic review and meta-analysis. EClinicalMedicine. 2020;29:100630.

    Article  PubMed  Google Scholar 

  5. Magesh S, John D, Li WT, et al. Disparities in COVID-19 outcomes by race, ethnicity, and socioeconomic status: a systematic-review and Meta-analysis. JAMA Netw Open. 2021;4(11):e2134147.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Updating ethnic contrasts in deaths involving the coronavirus (COVID-19), England and Wales: deaths occurring 2 March to 28 July 2020. https://www.ons.gov.uk/peoplepopulationandcommunity/birthsdeathsandmarriages/deaths/articles/updatingethniccontrastsindeathsinvolvingthecoronaviruscovid19englandandwales/deathsoccurring2marchto28july2020 (Accessed 29 Apr 2021).

  7. Price-Haywood EG, Burton J, Fort D, Seoane L. Hospitalization and mortality among black patients and white patients with Covid-19. N Engl J Med. 2020;382(26):2534–43.

    Article  CAS  PubMed  Google Scholar 

  8. Choi KH, Denice P, Haan M, Zajacova A. Studying the social determinants of COVID-19 in a data vacuum. Can Rev Sociol. 2021;58(2):146–64.

    Article  PubMed  PubMed Central  Google Scholar 

  9. COVID-19 in Immigrants, Refugees and Other Newcomers in Ontario: Characteristics of Those Tested and Those Confirmed Positive, as of June 13, 2020. 2020. https://www.ices.on.ca/Publications/Atlases-and-Reports/2020/COVID-19-in-Immigrants-Refugees-and-Other-Newcomers-in-Ontario.

  10. Pareek M, Bangash MN, Pareek N, et al. Ethnicity and COVID-19: an urgent public health research priority. Lancet. 2020;395(10234):1421–2.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Chudasama YV, Zaccardi F, Gillies CL, et al. Patterns of multimorbidity and risk of severe SARS-CoV-2 infection: an observational study in the U.K. BMC Infect Dis. 2021;21(1):908.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Agrawal U, Azcoaga-Lorenzo A, Fagbamigbe AF, et al. Association between multimorbidity and mortality in a cohort of patients admitted to hospital with COVID-19 in Scotland. J R Soc Med. 2022;115(1):22–30.

    Article  PubMed  Google Scholar 

  13. Khunti K, Knighton P, Zaccardi F, et al. Prescription of glucose-lowering therapies and risk of COVID-19 mortality in people with type 2 diabetes: a nationwide observational study in England. Lancet Diab Endocrinol. 2021;9(5):293–303.

    Article  CAS  Google Scholar 

  14. MacKenna B, Kennedy NA, Mehrkar A, et al. Risk of severe COVID-19 outcomes associated with immune-mediated inflammatory diseases and immune-modifying therapies: a nationwide cohort study in the OpenSAFELY platform. Lancet Rheumatol. 2022;4(7):e490–506.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Paakkari L, Okan O. COVID-19: health literacy is an underestimated problem. Lancet Public Health. 2020;5(5):e249–e50.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Travaglio M, Yu Y, Popovic R, Selley L, Leal NS, Martins LM. Links between air pollution and COVID-19 in England. Environ Pollut. 2021;268(Pt A):115859.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Upshaw TL, Brown C, Smith R, Perri M, Ziegler C, Pinto AD. Social determinants of COVID-19 incidence and outcomes: a rapid review. PLoS One. 2021;16(3):e0248336.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Razieh C, Zaccardi F, Islam N, et al. Ethnic minorities and COVID-19: examining whether excess risk is mediated through deprivation. Eur J Pub Health. 2021;31(3):630–4.

    Article  Google Scholar 

  19. Amato M, Werba JP, Frigerio B, et al. Relationship between influenza vaccination coverage rate and COVID-19 outbreak: an Italian ecological study. Vaccines (Basel). 2020;8(3).

  20. Khunti K, Singh AK, Pareek M, Hanif W. Is ethnicity linked to incidence or outcomes of covid-19? BMJ. 2020;369:m1548.

    Article  PubMed  Google Scholar 

  21. Zanettini C, Omar M, Dinalankara W, et al. Influenza Vaccination and COVID19 Mortality in the USA. medRxiv. 2020.

  22. Shah BR, Chiu M, Amin S, Ramani M, Sadry S, Tu JV. Surname lists to identify south Asian and Chinese ethnicity from secondary data in Ontario, Canada: a validation study. BMC Med Res Methodol. 2010;10:42.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Elovainio M, Hakulinen C, Pulkki-Raback L, et al. Contribution of risk factors to excess mortality in isolated and lonely individuals: an analysis of data from the UK biobank cohort study. Lancet Public Health. 2017;2(6):e260–e6.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Lambert PC, Royston P. Further development of flexible parametric models for survival analysis. Stata J. 2009;9(2):265–90.

    Article  Google Scholar 

  25. Royston P, Parmar MK. Flexible parametric proportional-hazards and proportional-odds models for censored survival data, with application to prognostic modelling and estimation of treatment effects. Stat Med. 2002;21(15):2175–97.

    Article  PubMed  Google Scholar 

  26. Rubin D. Multiple imputation for nonresponse in surveys: John Wiley; 1987.

    Book  Google Scholar 

  27. Navaratnam AV, Gray WK, Day J, Wendon J, Briggs TWR. Patient factors and temporal trends associated with COVID-19 in-hospital mortality in England: an observational study using administrative data. Lancet Respir Med. 2021;9(4):397–406.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Yates T, Summerfield A, Razieh C, et al. A population-based cohort study of obesity, ethnicity and COVID-19 mortality in 12.6 million adults in England. Nat Commun. 2022;13(1):624.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Horton R. Offline: COVID-19 is not a pandemic. Lancet. 2020;396(10255):874.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Hawkins D. Differential occupational risk for COVID-19 and other infection exposure according to race and ethnicity. Am J Ind Med. 2020;63(9):817–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Razai MS, Kankam HKN, Majeed A, Esmail A, Williams DR. Mitigating ethnic disparities in covid-19 and beyond. BMJ. 2021;372:m4921.

    Article  PubMed  Google Scholar 

  32. Gray WK, Navaratnam AV, Day J, Wendon J, Briggs TWR. Changes in COVID-19 in-hospital mortality in hospitalised adults in England over the first seven months of the pandemic: an observational study using administrative data. Lancet Reg Health Eur. 2021;5:100104.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Griffith GJ, Morris TT, Tudball MJ, et al. Collider bias undermines our understanding of COVID-19 disease risk and severity. Nat Commun. 2020;11(1):5749.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Raharja A, Tamara A, Kok LT. Association between ethnicity and severe COVID-19 disease: a systematic review and Meta-analysis. J Racial Ethn Health Disparities. 2021;8(6):1563–72.

    Article  PubMed  Google Scholar 

  35. VanderWeele TJ. Mediation analysis: a Practitioner's guide. Annu Rev Public Health. 2016;37:17–32.

    Article  PubMed  Google Scholar 

  36. Babyak MA. Understanding confounding and mediation. Evid Based Ment Health. 2009;12(3):68–71.

    Article  PubMed  Google Scholar 

  37. Coronavirus (COVID-19) related deaths by occupation, England and Wales: deaths registered between 9 March and 28 December 2020. https://www.ons.gov.uk/peoplepopulationandcommunity/healthandsocialcare/causesofdeath/bulletins/coronaviruscovid19relateddeathsbyoccupationenglandandwales/deathsregisteredbetween9marchand28december2020 (accessed 13 May 2021).

  38. Mathur R, Bhaskaran K, Chaturvedi N, et al. Completeness and usability of ethnicity data in UK-based primary care and hospital databases. J Public Health (Oxf). 2014;36(4):684–92.

    Article  PubMed  Google Scholar 

  39. Caleyachetty R, Barber TM, Mohammed NI, et al. Ethnicity-specific BMI cutoffs for obesity based on type 2 diabetes risk in England: a population-based cohort study. Lancet Diabetes Endocrinol. 2021;9(7):419–26.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Rentsch CT, Kidwai-Khan F, Tate JP, et al. Patterns of COVID-19 testing and mortality by race and ethnicity among United States veterans: a nationwide cohort study. PLoS Med. 2020;17(9):e1003379.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

Parts of this material are based on data and/or information compiled and provided by the Canadian Institute for Health Information (CIHI). The analyses, conclusions, opinions and statements expressed herein are solely those of the authors and do not reflect those of the MOH, the MLTC or CIHI; no endorsement is intended or should be inferred. We acknowledge the contribution of EMIS practices who contribute to QResearch and EMIS Health and the Universities of Nottingham and Oxford for expertise in establishing, developing, or supporting the QResearch database. This project involves data derived from patient level information collected by the NHS, as part of the care and support of cancer patients. The data are collated, maintained, and quality assured by the National Cancer Registration and Analysis Service, which is part of Public Health England (PHE). Access to the data was facilitated by the PHE Office for Data Release. The Hospital Episode Statistics data used in this analysis are reused by permission from NHS Digital, which retains the copyright in that data. We thank the Office for National Statistics (ONS) for providing the mortality data. NHS Digital, PHE, and the ONS bear no responsibility for the analysis or interpretation of the data. The investigators acknowledge the philanthropic support of the donors to the University of Oxford’s COVID-19 Research Response Fund. The views and opinions expressed by authors in this publication are those of the authors and do not necessarily reflect those of the UK National Institute for Health Research (NIHR) or the Department of Health and Social Care.

Funding

This study (UK) is jointly funded by UKRI and NIHR [COV0130 /MR/V027778/1].

SJG is supported by an MRC Epidemiology Unit programme: MC_UU_12015/4. The University of Cambridge has received salary support in respect of SJG from the NHS in the East of England through the Clinical Academic Reserve.

AKC is funded by a Clinical Research Fellowship from Cancer Research UK (C2195/A31310).

JHC has received grants from the NHS, National Institute for Health Research, Oxford, John Fell Oxford University Press Research Fund, Cancer Research UK (grant number C5255/A18085) through the Cancer Research UK Oxford Centre, and the Oxford Wellcome Institutional Strategic Support Fund (204826/Z/16/Z) during the conduct of the study.

KK and FZ are supported by the National Institute for Health Research (NIHR) Applied Research Collaboration East Midlands (ARC EM) and the NIHR Leicester Biomedical Research Centre (BRC). [NIHR200171]. HDM is a National Institute for Health Research funded Academic Clinical Lecturer.

This study (Canada) was supported by ICES, which is funded by an annual grant from the Ontario Ministry of.

Health (MOH) and the Ministry of Long-Term Care (MLTC).

Author information

Author notes

  1. Francesco Zaccardi, Pui San Tan and Baiju R. Shah are joint first authors.

  2. Hajira Dambha-Miller and Julia Hippisley-Cox are joint senior authors.

Authors and Affiliations

  1. Leicester Real World Evidence Unit, Leicester Diabetes Centre, University of Leicester, Leicester, England

    Francesco Zaccardi & Kamlesh Khunti

  2. Nuffield Department of Primary Care Health Sciences, University of Oxford, Oxford, England

    Pui San Tan, Ash Kieran Clift, Martina Patone, Defne Saatci, Carol Coupland & Julia Hippisley-Cox

  3. Department of Medicine, University of Toronto; Division of Endocrinology, Sunnybrook Health Sciences Centre, Institute for Clinical Evaluative Sciences, Toronto, Canada

    Baiju R. Shah & Karl Everett

  4. Cancer Research UK Oxford Centre, Department of Oncology, University of Oxford, Oxford, England

    Ash Kieran Clift

  5. Division of Primary Care, School of Medicine, University of Nottingham, Nottingham, England

    Carol Coupland

  6. Primary Care Unit, School of Clinical Medicine, University of Cambridge, Cambridge, England

    Simon J. Griffin

  7. MRC Epidemiology Unit, School of Clinical Medicine, University of Cambridge, Cambridge, England

    Simon J. Griffin

  8. Primary Care Research Centre, University of Southampton, Southampton, England

    Hajira Dambha-Miller

Authors
  1. Francesco Zaccardi
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Pui San Tan
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Baiju R. Shah
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Karl Everett
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Ash Kieran Clift
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Martina Patone
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Defne Saatci
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Carol Coupland
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Simon J. Griffin
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. Kamlesh Khunti
    View author publications

    You can also search for this author in PubMed Google Scholar

  11. Hajira Dambha-Miller
    View author publications

    You can also search for this author in PubMed Google Scholar

  12. Julia Hippisley-Cox
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

Study conceptualisation was led by HDM and JHC. All authors contributed to development of research questions and formulated the study design. FZ, PST, BRS, KE and AKC performed the statistical analyses. All authors contributed to interpretation of results. Clinical code groups were developed and checked by JHC, AKC, PST (QResearch) and BRS, KE (Canada). FZ, PST, AKC and HDM wrote the first draft of the manuscript, which was critically revised for important intellectual content by all other authors. All authors approved the final version of the manuscript. The corresponding author attests that all listed authors meet authorship criteria and that no others meeting the criteria have been omitted. JHC is the guarantor.

Corresponding author

Correspondence to Julia Hippisley-Cox.

Ethics declarations

Ethics approval and consent to participate

This project used data routinely collected from clinical purpose and no experiments were conducted. The research has been conducted in accordance with the relevant following regulations:

Ontario - The consent was unnecessary according to national regulations. ICES is a prescribed entity under Ontario’s Personal Health Information Protection Act (PHIPA). According to the national regulation, section 45 of PHIPA authorizes ICES to collect personal health information, without consent, for the purpose of analysis or compiling statistical information with respect to the management of, evaluation or monitoring of, the allocation of resources to or planning for all or part of the health system (https://www.ontario.ca/laws/statute/04p03). Projects that use data collected by ICES under section 45 of PHIPA are exempt from REB (Research Ethics Board) review. The use of the data in this project was reviewed by ICES’ Privacy and Legal Office and approved under section 45. The data were anonymized before its use. Storage, handling and access to these data are governed by policies and procedures that are regularly reviewed and approved by the Information and Privacy Commissioner of Ontario.

QResearch - Access to data for this study has been granted ethics approval by East Midlands-Derby Research Ethics Committee (reference 18/EM/0400). This study has also been independently peer-reviewed by the QResearch Scientific board (project reference OX102). Consent for QResearch was not necessary according to national regulations as the database is anonymised and no patient consent is required: as this study used routinely collected data, no consent is needed but patients can opt out if they prefer their data not to be used. Data are anonymised in which no individuals can be identified.

Consent for publication

Not applicable. To preserve confidentiality, in tables cells with frequencies less than 5 have been suppressed.

Competing interests

KK is a Member of the Scientific Advisory Group for Emergencies (SAGE), Member of Independent SAGE, Director of the University of Leicester Centre for Black Minority Health and Trustee of the South Asian Health Foundation. JHC is a member of several SAGE committees and chair of the risk stratification subgroup of the NERVTAG. She is unpaid director of QResearch and founder and former medical director of ClinRisk Ltd. (outside the submitted work). PST reports past consultation with AstraZeneca and Duke-NUS outside the submitted work. All other authors have no competing interests to declare.

JHC is an unpaid director of QResearch, a not-for-profit organisation which is a partnership between the University of Oxford and EMIS Health who supply the QResearch database used for this work, and is a founder and shareholder of ClinRisk Ltd. and was its medical director until 31 May 2019; ClinRisk produces open and closed source software to implement clinical risk algorithms (outside this work) into clinical computer systems.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zaccardi, F., Tan, P.S., Shah, B.R. et al. Ethnic disparities in COVID-19 outcomes: a multinational cohort study of 20 million individuals from England and Canada. BMC Public Health 23, 399 (2023). https://doi.org/10.1186/s12889-023-15223-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12889-023-15223-8

Keywords

BMC Public Health

ISSN: 1471-2458

Contact us