Skip to main content
  • Research article
  • Open access
  • Published:

Ethnic differences in COVID-19 mortality in the second and third waves of the pandemic in England during the vaccine rollout: a retrospective, population-based cohort study

Abstract

Background

Ethnic minority groups in England have been disproportionately affected by the COVID-19 pandemic and have lower vaccination rates than the White British population. We examined whether ethnic differences in COVID-19 mortality in England have continued since the vaccine rollout and to what extent differences in vaccination rates contributed to excess COVID-19 mortality after accounting for other risk factors.

Methods

We conducted a retrospective, population-based cohort study of 28.8 million adults aged 30–100 years in England. Self-reported ethnicity was obtained from the 2011 Census. The outcome was death involving COVID-19 during the second (8 December 2020 to 12 June 2021) and third wave (13 June 2021 to 1 December 2021). We calculated hazard ratios (HRs) for death involving COVID-19, sequentially adjusting for age, residence type, geographical factors, sociodemographic characteristics, pre-pandemic health, and vaccination status.

Results

Age-adjusted HRs of death involving COVID-19 were elevated for most ethnic minority groups during both waves, particularly for groups with lowest vaccination rates (Bangladeshi, Pakistani, Black African, and Black Caribbean). HRs were attenuated after adjusting for geographical factors, sociodemographic characteristics, and pre-pandemic health. Further adjusting for vaccination status substantially reduced residual HRs for Black African, Black Caribbean, and Pakistani groups in the third wave. Fully adjusted HRs only remained elevated for the Bangladeshi group (men: 2.19 [95% CI 1.72–2.78]; women: 2.12 [1.58–2.86]) and Pakistani men (1.24 [1.06–1.46]).

Conclusions

Lower COVID-19 vaccination uptake in several ethnic minority groups may drive some of the differences in COVID-19 mortality compared to White British. Public health strategies to increase vaccination uptake in ethnic minority groups would help reduce inequalities in COVID-19 mortality, which have remained substantial since the start of the vaccination campaign.

Peer Review reports

Background

The disproportionate impact of the coronavirus pandemic on ethnic minority groups has been widely reported [1,2,3,4]. In England, rates of hospitalisation for COVID-19, admission to intensive care, and death were higher among ethnic minority groups during the first and second waves of the pandemic [5]. However, compared with the first wave, excess COVID-19 mortality was reduced in the second wave among the Black African and Black Caribbean groups but increased among the Pakistani and Bangladeshi groups [6, 7]. Moreover, adjustments for geography, socioeconomic factors, and pre-existing health conditions accounted for a large proportion of the elevated COVID-19 mortality risk observed in the first and second waves. However, some residual risk remained unexplained, most notably for South Asian and Black African groups.

The UK began its coronavirus vaccination programme on 8 December 2020, starting with those most likely to experience severe outcomes (people aged 70 years and over and those with underlying health conditions) and those working in health and social care roles [8]. Survey data indicates that rates of vaccine hesitancy in the UK are highest among people from Black ethnic groups [9]. Consistent with this, differences in vaccination rates by ethnic group were evident early in the rollout of the vaccination programme and these differences widened over time [10, 11]; by 12 December 2021, 3.7% of White British adults aged 50 or over had not received any dose of a COVID-19 vaccine, compared with 26.2% from the Black Caribbean group and 17.4% of the Black African group [12].

This study used population-level linked administrative data for England to investigate whether inequalities in deaths involving COVID-19 by ethnic group have continued into the third wave. We also explored the extent to which elevated COVID-19 mortality in some ethnic groups can be explained by differences in age, residence type, geographical factors, sociodemographic characteristics, pre-pandemic health, and, for the first time, vaccination rates.

Methods

Data sources

We conducted a retrospective, population-based cohort study using data from the Office for National Statistics (ONS) Public Health Data Asset (PHDA). The ONS PHDA is a linked dataset combining the 2011 Census, mortality records, the General Practice Extraction Service (GPES) Data for Pandemic Planning and Research (GDPPR), Hospital Episode Statistics (HES), and vaccination data from the National Immunisation Management System (NIMS).

Person-level datasets were created from the HES and GDPPR record-level datasets by stacking and deduplicating on NHS number and date of birth. Records with blank or invalid NHS numbers or dates of birth were dropped as these could not be linked to the 2011 Census.

Data linkage

To obtain NHS numbers for the 2011 Census, we linked the 2011 Census to the 2011–2013 NHS Patient Registers using deterministic and probabilistic matching, with an overall linkage rate of 94.6% (see [13] for a detailed description of the linkage methodology and quality evaluation). Rates of linkage failure were higher for most ethnic minority groups than the White British group (Additional file 1: Table S1). Compared with the White British group, the unadjusted odds of linkage failure were highest for the ‘other’ ethnic group (odds ratio = 5.81, 95% confidence interval = 5.78 to 5.84), followed by the mixed ethnic group (4.44, 4.40 to 4.47) and the Chinese group (4.11, 4.07 to 4.16). Rates of linkage failure were also higher in men, younger age groups, people living in more deprived areas, and varied by region; the highest rate of linkage failure was observed for people living in London (9.2%). Once these factors were adjusted for, the odds of linkage failure in ethnic minority groups were substantially reduced. The adjusted odds ratio of linkage failure was below one for the Indian, Bangladeshi, and Pakistani ethnic groups but remained above one for the Black African (1.76, 1.74 to1.77) and the Black Caribbean (1.30, 1.28 to 1.31) ethnic groups.

Further linkage to deaths registrations data, GDPPR, HES, and NIMS data was performed deterministically based on NHS number. 86.2% of deaths recorded in deaths registrations data that occurred in England between 8 December 2020 and 1 December 2021 among people aged 30 to 100 years were included in the analysis. 79.0% of people aged 30 to 100 years who had received at least one dose of a COVID-19 vaccine during the study period according to NIMS records as of 4 January 2022 were linked to the PHDA. The unlinked deaths and vaccination records reflect people not included in our study population (see the ‘Generation of the study cohort’ section).

Generation of the study cohort

Of the 41,880,933 people enumerated at the 2011 Census in England and Wales who would be aged 30–100 in 2020, we excluded 354,036 people (0.9%) who were short-term residents (i.e. people who were enumerated at the 2011 Census but did not intend on staying in the country for at least 12 months), 2,257,221 people (5.4%) who could not be linked deterministically or probabilistically to the NHS Patient register, and 4,360,949 individuals (10.4%) who had died between the Census and 8 December 2020 (the start of the vaccination campaign in the UK). An additional 6,092,707 people (14.5%) were not linked to English primary care records because they either did not live in England in 2019 (the Census included people living in England and Wales) or were not registered with a GP practice in England that was participating in GDPPR (which collects data from 6535 GP practices, covering all GP system suppliers and 97.5% of open and active practices in England [14]). The final study population included 28,816,020 people aged 30–100 years in 2020 (representing 80% of the mid-year 2020 population estimate for England) (Fig. 1) [15].

Fig. 1
figure 1

Sample selection and number of participants

We restricted our analysis to people aged 30 to 100 in 2020 because most sociodemographic factors were drawn from the 2011 Census, which may not represent people’s circumstances at the beginning of the pandemic; younger people were thought particularly likely to have changed their circumstances, as evidenced by a greater proportion of younger people having different postcodes at the 2011 Census to the most recent postcode recorded in their GP records (Additional file 1: Fig. S1). In addition, very few deaths occurred in people aged below 30 years; official figures show that out of the 84,449 deaths involving COVID-19 in England and Wales in 2020, only 127 (0.2%) were among people less than 30 years old [16].

Exposure

The exposure was self-reported ethnic group, retrieved from the 2011 Census. We used a 10-category ethnic group classification (White British [White English/Welsh/Scottish/Northern Irish/British], Bangladeshi, Black African, Black Caribbean, Chinese, Indian, mixed [White and Asian, White and Black African, White and Black Caribbean, and other mixed], Pakistani, White other [White Irish, White Gypsy or Irish Traveller, and other White], and other [other Asian, Arab, other Black, and any other ethnic groups]).

Covariates

Ethnic differences in the risk of death involving COVID-19 across could be mediated by factors linked to the risk of infection (such as communal living, residing in an area with high infection rates, socioeconomic and demographic factors, and occupation) and factors associated the risk of death if infected (health status and vaccination status). These factors may fall on the causal pathway between ethnicity and COVID-19 mortality (Additional file 1: Fig. S2).

The following covariates were included from the 2011 Census data: age, residence type (private household, care home, or other communal establishment), household tenure, National Statistics Socio-economic Classification (NS-SEC), highest qualification, household size, household deprivation, family status, household composition, key worker in household, and key worker type (education and childcare, food and necessity goods, health and social care, key public services, national and local Government, public safety and national security, transport, utilities and communication; derived according to Census returns based on the 2010 Standard Occupational Classification (SOC) and the 2007 Standard Industrial Classification of Economic Activities [17]) (Additional file 1: Table S2). Body mass index (BMI; classified as underweight, normal weight, overweight, obese, unknown) and pre-existing health conditions (using the same Systematized Nomenclature of Medicine Clinical Terms [SNOMED-CT] codes as the QCovid2 risk prediction model [18]) were included as covariates from the GDPPR data. The QCovid risk prediction model was used in the UK to identify clinically extremely vulnerable individuals who should shield during the pandemic. The model has been previously shown to predict risk of COVID-19 hospitalisation and mortality in three independent datasets [19,20,21]. The number of admissions to, and number of days spent in, admitted patient care during the three years prior to the pandemic were included as covariates from the HES Admitted Patient Care data.

The following covariates were included from other data sources: vaccination status (from NIMS data); region and Rural Urban classification (from the National Statistics Postcode Lookup), population density of the Lower layer Super Output Area (from mid-2019 population estimates), and Index of Multiple Deprivation (IMD) (from the English Indices of Deprivation, 2019) derived from postcodes in GDPPR data; occupational exposure to disease and proximity to others for individuals and the maximum score among all individuals in each household (from the Occupational Information Network database, which collects a range of information about individuals’ working conditions and day-to-day tasks of their job. To calculate the proximity and exposure measures, the questions asked were as follows: (i) How physically close to other people are you when you perform your current job? (ii) How often does your current job require that you be exposed to diseases or infection? Scores ranging from 0 to 100 were calculated for these questions based on 2011 Census data on occupation [22]) and care-home residence status (from the 2019 NHS Patient Register).

Outcome

The outcome was death involving COVID-19, i.e. COVID-19 International Classification of Diseases 10 code of U07.1 (COVID-19, virus identified), U07.2 (COVID-19, virus not identified), or U09.9 (post-COVID condition, unspecified) in part I or II of the death certificate, occurring between 8 December 2020 and 1 December 2021.

Statistical analysis

Age-standardised vaccination rates by ethnic group were calculated for each wave based on the number of vaccine doses received by the end of the period (12 June 2021 for wave 2 and 1 December 2021 for wave 3). Crude rates (%) of the number of people in each ethnic group who were unvaccinated, single-vaccinated, double-vaccinated, or triple-vaccinated were age-standardised using the 2013 European Standardised Population [23].

We calculated age-standardised mortality rates (ASMRs) by ethnic group as deaths per 100,000 person-years at risk to examine the absolute risk of death involving COVID-19, standardised to the 2013 European Standardised Population [23]. ASMRs were calculated separately for each of the waves of the pandemic that occurred during the vaccine rollout (wave 2: 8 December 2020 to 12 June 2021; wave 3: 13 June 2021 to 1 December 2021). This analysis therefore excludes any deaths occurring early in the second wave, which is estimated to have started in September 2020 [24].

As the pandemic was ongoing at the end of the study period, the data were subject to right-censoring. We therefore used Cox proportional hazards models to assess whether differences in the risk of mortality involving COVID-19 by ethnic group could be accounted for by location, sociodemographic factors, pre-pandemic health, and vaccination status. Separate models were fitted for the second and third waves. The index date for start of follow-up time was 8 December 2020 for wave two (the start of the vaccination programme in the UK) and 13 June 2021 for wave three. End of follow-up was date of death for those who died or the end of the wave period for those who were still alive at the end of the period: 12 June 2021 for wave two and 1 December 2021 for wave three (see Additional file 1: Table S3 for mean follow-up times by ethnic group). For computational efficiency, we included all individuals who died of any cause during the analysis period and a random sample (selected by simple random sampling without replacement) of those who did not, with sampling rates of 1% for the White British ethnic group and 10% for every other ethnic group; case weights equal to the inverse probability of selection were included in the analysis, following previously published methods [6, 13]. The White British group was used as the reference category in all models.

The baseline model (model 1) only included adjustment for single year of age as a confounding variable, included as a second-order polynomial. We then introduced potential mediating factors sequentially, starting with factors associated with the risk of exposure to SARS-CoV-2 and then factors associated the risk of death if infected. Model 2 included additional adjustment for type of residence (private household, care home, or other communal establishments). In model 3, we included additional adjustment for geographical factors (region, Rural Urban classification and local population density). In model 4, we adjusted for socioeconomic and demographic factors that are likely to be linked to risk of infection (NS-SEC, highest qualification, IMD decile, household characteristics [tenure of the household, household deprivation, household size, family status, household composition, and key worker in the household], key worker type, individual and household exposure to disease, and individual and household proximity to others). We then adjusted for factors associated with the risk of death if infected. In model 5, additional adjustment was made for health status (pre-existing health conditions, BMI, and number of admissions to hospital and days spent in hospital over the previous 3 years). For all health variables, a binary interaction indicator was included, allowing the effects to vary depending on whether the individual was aged 70 years and older or younger than 70 years. In model 6, vaccination status (unvaccinated, one dose or two doses for wave two plus three doses for wave three) was included as a time-varying covariate, based on the date of vaccination plus 14 days. Therefore, an individual was classified as single-vaccinated 14 days after they received their first vaccine dose, double-vaccinated 14 days after they received their second dose, and triple-vaccinated 14 days after they received their third dose.

Missing Census data were imputed using nearest-neighbour donor imputation, the methodology employed by the Office for National Statistics across all 2011 Census variables [25]. Ethnicity was imputed in 3.0% of 2011 Census records due to item non-response. Individuals with missing data for BMI were placed into an unknown category. Health conditions were derived based on prescription and diagnosis codes, with the sample restricted to people who were registered with a GP practice in England that was participating in GDPPR. Therefore, there were no missing values.

All statistical analyses were stratified by sex and conducted using R, version 3.5. Cox proportional hazards models were implemented using the survival package (version 2.41-3) [26].

Results

Characteristics of the study population

The study population included 28,816,020 adults aged 30–100 years (mean age 55.8 years, SD 15.6) in England, 13,568,656 (47.1%) of whom were male (Table 1).

Table 1 Demographic and medical characteristics for the study cohort

446,093 deaths were recorded during follow-up, of which 66,558 were involving COVID-19 (54,770 in wave 2 and 11,788 in wave 3). Age-adjusted rates of receiving at least one COVID-19 vaccine dose by the end of the second wave (12 June 2021) were lowest among the Black Caribbean (41.0% unvaccinated) and Black African (31.5% unvaccinated) groups and highest among White British (10.1% unvaccinated) and Indian (12.0% unvaccinated) groups (Fig. 2). Ethnic inequalities in vaccination coverage continued into the third wave with the rollout of third/booster doses. As of 1 December 2021, age-adjusted rates of people receiving third/booster doses were lowest among the Pakistani group (24.6%) and highest for Indian (48.3%) and White British (46.1%) groups.

Fig. 2
figure 2

Age-standardised vaccination rates by ethnic group and wave of the pandemic. Percentage of people in each ethnic group that were unvaccinated or had received one or two doses of a COVID-19 vaccine by 12 June 2021 (left panel). Percentage of people in each ethnic group that were unvaccinated or had received one, two or three doses by 1 December 2021 (right panel)

Differences in COVID-19 mortality: age-standardised mortality rates

Overall, the ASMRs for death involving COVID-19 were higher than the White British group for all ethnic minority groups, except for the Chinese group and women in the White other group. Rates were highest for the Bangladeshi group and lowest for the Chinese group, men from the White British group, and women from the White other group (Table 2). Further disaggregating the ethnic groups showed that, compared with the White British group, ASMRs were elevated for the ‘other Asian’, ‘other Black’, and ‘any other’ ethnic groups and men from the Arab and White Gypsy or Irish Traveller groups (Additional file 1: Table S4). Breaking down the ASMRs by wave of the pandemic revealed that during both the second and third waves, the rate of death involving COVID-19 was consistently highest for the Bangladeshi group. Excess COVID-19 mortality relative to the White British group increased from the second to third wave for women from the Pakistani group and, to a lesser extent, women from the Bangladeshi and Black Caribbean group.

Table 2 Number of deaths and ASMRs (per 100,000 person-years) for deaths involving COVID-19, stratified by sex, ethnic group, and wave of the pandemic

Understanding the differences in COVID-19 mortality between ethnic groups

Second wave

Age-adjusted hazard ratios (HRs) were calculated to examine the rate of death involving COVID-19 for ethnic minority groups relative to the White British group, holding the effect of age constant to account for the younger age distribution of ethnic minority groups. Compared with the White British group, rates of death involving COVID-19 were higher for all ethnic minority groups (except for the Chinese group and women from the ‘White other’ group) during the second wave (Fig. 3; model 1). HRs were highest for the Bangladeshi group, at 5.03 (95% confidence interval [CI] 4.51 to 5.60) for men and 4.48 (3.91 to 5.13) for women. HRs were substantially reduced after adjustment for geographical factors, socioeconomic status, and pre-existing health conditions (model 5) but remained elevated for the Bangladeshi (men: 2.27, 2.01 to 2.56; women: 2.14, 1.84 to 2.50), Pakistani (men: 1.76, 1.61 to 1.91; women: 1.52, 1.35 to 1.70), Indian (men: 1.53, 1.43 to 1.65; women: 1.30, 1.19 to 1.42), and ‘other’ (men: 1.27, 1.16 to 1.39; women: 1.31, 1.18 to 1.45) groups and for men from the Black African (1.55, 1.37 to 1.76) and Black Caribbean groups (1.13, 1.03 to 1.24). After further adjustment for vaccination status (model 6), the rate of death involving COVID-19 for men from the Black Caribbean and White British groups were similar. However, excess risk remained for the Bangladeshi (men: 2.20, 1.95 to 2.47; women: 2.08, 1.78 to 2.43), Pakistani (men: 1.62, 1.49 to 1.77; women: 1.41, 1.26 to 1.58), Indian (men: 1.54, 1.44 to 1.66; women: 1.32, 1.21 to 1.45), and ‘other’ (men: 1.23, 1.12 to 1.34; women: 1.28, 1.15 to 1.42) groups and for men from the Black African group (1.48, 1.30 to 1.67).

Fig. 3
figure 3

HRs for death involving COVID-19 by ethnic group during the second wave of the pandemic (8 December 2020 to 12 June 2021), relative to the White British group, stratified by sex. Results obtained from Cox proportional hazards regression models adjusted for the following: model 1—age; model 2—age plus residence type (private household, care home or other communal establishment); model 3—age and residence type plus geographical factors (region, Rural Urban classification and population density); model 4—age, residence type, and geography, plus sociodemographic factors (highest qualification, IMD decile, NS-SEC, household characteristics [tenure of the household, household deprivation, household size, family status, household composition and key worker in household], key worker type, individual and household exposure to disease, and individual and household proximity to others); model 5—age, residence type, geography, and sociodemographic factors, plus health status (pre-existing health conditions, BMI and hospital admissions over the previous 3 years); and model 6—age, residence type, geography, sociodemographic factors, and health status plus vaccination status (unvaccinated, one dose or two doses). Error bars represent 95% CIs

Third wave

During the third wave, age-adjusted HRs were elevated for all ethnic minority groups except for the Chinese group, men from the ‘mixed’ group, and women from the ‘White other’ group (Fig. 4, model 1). HRs continued to be highest for the Bangladeshi group (men: 4.43, 3.54 to 5.53; women: 5.23, 4.02 to 6.80). After adjusting for geographical factors, socioeconomic status, and pre-existing health conditions (model 5), HRs remained elevated for Bangladeshi (men: 2.49, 1.96 to 3.17; women: 2.17, 1.61 to 2.93), Pakistani (men: 1.71, 1.46 to 2.01; women: 1.62, 1.32 to 1.99), Black Caribbean (men: 1.70, 1.43 to 2.04; women: 2.12, 1.76 to 2.55), Black African (men: 1.40, 1.06 to 1.85; women: 1.80, 1.38 to 2.35), and ‘other’ groups (men: 1.33, 1.11 to 1.61; women: 1.46, 1.17 to 1.83) and men from the ‘White other’ group (1.15, 1.01 to 1.31). After additional adjustment for vaccination status (model 6), HRs remained elevated for the Bangladeshi group (men: 2.19, 1.72 to 2.78; women: 2.12, 1.58 to 2.86) and men from the Pakistani group (1.24, 1.06 to 1.46), whereas rates of death involving COVID-19 for all other groups were similar to the White British group.

Fig. 4
figure 4

HRs for death involving COVID-19 by ethnic group during the third wave of the pandemic (13 June 2021 to 1 December 2021), relative to the White British group, stratified by sex. Results obtained from Cox proportional hazards regression models adjusted for the following: model 1—age; model 2—age plus residence type (private household, care home or other communal establishment); model 3—age and residence type plus geographical factors (region, Rural Urban classification and population density); model 4—age, residence type, and geography, plus sociodemographic factors (highest qualification, IMD decile, NS-SEC, household characteristics [tenure of the household, household deprivation, household size, family status, household composition and key worker in household], key worker type, individual and household exposure to disease, and individual and household proximity to others); model 5—age, residence type, geography, and sociodemographic factors, plus health status (pre-existing health conditions, BMI and hospital admissions over the previous 3 years); and model 6—age, residence type, geography, sociodemographic factors, and health status plus vaccination status (unvaccinated, one dose, two doses or third/booster dose). Error bars represent 95% CIs

Discussion

Our findings demonstrate that throughout the rollout of the vaccine programme, most ethnic minority groups have continued to experience greater rates of death involving COVID-19 compared with the White British group. Although the patterns of excess COVID-19 mortality risk by ethnic group have changed over the course of the pandemic, the Bangladeshi, Black African, Black Caribbean, and Pakistani groups remained the groups with highest rate of COVID-19 mortality in the third wave. We also found that adjusting for vaccination status explained some of the remaining increased risk of COVID-19 mortality that is not explained by other factors, particularly during the third wave for the Black African, Black Caribbean, and Pakistani groups. After adjusting for geographical factors, sociodemographic characteristics, pre-pandemic health, and vaccination status, the risk of COVID-19 mortality was similar to the White British group for all ethnic groups except the Bangladeshi group and men from the Pakistani group.

Previous analyses showed that differences in location, measures of disadvantage, occupation, living arrangements, and certain pre-existing health conditions explain a large proportion (but not all) of the excess COVID-19 mortality risk observed in some ethnic groups [6, 7, 13]. Based on these findings, we hypothesised that ethnic differences in COVID-19 mortality could be mediated by factors associated with the risk of exposure to SARS-CoV-2 (e.g. geographical factors, socioeconomic status, living arrangements) and the risk of severe outcomes once infected (e.g. pre-existing health conditions, vaccination status) (Additional file 1: Fig. S2). We found that adjustment for geographical factors was associated with the largest reduction in the HRs for most ethnic groups during the second wave. SARS-CoV-2 case rates in England were higher in more densely populated urban areas than in less densely populated rural areas during the second wave [27]. Therefore, the reduction in risk of COVID-19 mortality in ethnic minority groups following adjustment for geographical factors likely reflects increased risk of infection. Additional adjustment for socioeconomic and demographic characteristics further reduced the HRs for most ethnic groups, most notably for Bangladeshi and Pakistani groups. This suggests that socioeconomic inequalities also contribute to increased rates of COVID-19 mortality in ethnic minority groups, which is consistent with evidence from other countries [2, 4, 28]. Further adjusting for pre-existing health conditions had a modest effect on the HRs for most ethnic groups, with the largest reductions for Bangladeshi and Pakistani men during the second wave. The smaller reduction in HRs for other ethnic groups possibly indicates that differences in the prevalence of comorbidities associated with severe COVID-19 outcomes were already captured by adjusting for socioeconomic inequalities that are also associated with health status.

Our results are consistent with studies that investigated ethnic inequalities in SARS-CoV-2 infection and COVID-19 mortality in the first two waves of the pandemic. There is ample evidence that most ethnic minority groups were disproportionately affected in the first two waves of the pandemic [5,6,7, 13, 29]. This study shows that the excess rate of death involving COVID-19 observed among ethnic minority groups early in the pandemic has continued throughout the rollout of the vaccine programme until late 2021. This is despite SARS-CoV-2 case rates being higher among the White British population from summer 2021 onwards [30, 31]. However, among people aged 65 years and over, cumulative case rates between March 2020 and October 2021 were higher for Black and South Asian groups than the White British group [32], suggesting that the age distribution of cases may partly explain the continued higher rates of COVID-19 mortality in these groups.

In addition, our findings suggest that lower vaccination rates (especially among Black Caribbean, Black African and Pakistani groups) contributes to explaining why some ethnic groups are more likely to experience more severe COVID-19 outcomes once infected. However, residual unexplained risk remained in the Bangladeshi group and men from the Pakistani group, even after full adjustment. People from Pakistani and Bangladeshi groups are more likely to reside in deprived areas, in large households and in multigenerational families [33]. Living in large, overcrowded, multi-generational households is associated with increased risk of SARS-CoV-2 infection [34, 35], and there is some evidence that living in a multi-generational household explains some of the differences in mortality [36]. Differences in occupational exposure may also account for some of the differences in mortality between groups, as a higher proportion of men from Pakistani and Bangladeshi groups work in key worker roles, such as healthcare workers, taxi drivers, shopkeepers, and proprietors than any other ethnic group [37], and these occupations have been found to be at elevated risk of COVID-19 mortality [38]. Whilst our study adjusted for a range of sociodemographic factors, including household composition and occupational exposure, these variables were retrieved from the 2011 Census, which may not reflect the situation in 2020 accurately. Consequently, adjustment for these factors might have been incomplete, possibly contributing towards the residual association.

The main strength of our study derives from using the ONS Public Health Data Asset, a nationwide large-scale population-wide data source combining the 2011 Census, mortality records, the General Practice Extraction Service (GPES) Data for Pandemic Planning and Research (GDPPR), Hospital Episode Statistics (HES), and vaccination data from the National Immunisation Management System (NIMS). Unlike studies based only on electronic health records, our study relies on self-identified ethnicity, limiting the potential for exposure misclassification bias. The PHDA also contains both detailed sociodemographic characteristics, such as household composition, housing quality, and occupational exposure, and extensive information on pre-pandemic health based on primary care and hospital records. To our knowledge, our study is the first to use nationally representative population-based linked data to examine the association between ethnicity and COVID-19 mortality in the third wave of the pandemic and explore the role of differences in vaccination uptake as a potential additional explanatory factor for the differences in COVID-19 mortality.

The main limitation is that most sociodemographic characteristics included in our models reflect the situations of individuals as they were in 2011, not necessarily those at the start of the COVID-19 pandemic. To mitigate this, we excluded people aged less than 30 years old, whose circumstances are the most likely to have changed since the Census. We also updated place of residence based on information from primary care records. As a result, information on area deprivation, rural/urban classification, region, and care home residence were up-to-date at the beginning of the pandemic. In addition, measurement error is likely to be smaller for the people at greater risk, since the sociodemographic factors are less likely to have changed for older people than younger people. However, some measurement error may reduce the explanatory power of the sociodemographic factors and pre-existing conditions included in the model, thus reducing their effect on the hazard ratios. For example, people who have retired since the 2011 Census will have incorrect data for current occupation. We have no reason to believe that there was substantial misdiagnosis of COVID-19 where it was mentioned on death certificates, especially in the latter part of the pandemic when testing was widely available.

Another limitation is that the study population is limited to people enumerated at the 2011 Census and therefore did not include people who were living in England but did not respond to the Census (estimated to be 5% of households) [39]. We also excluded people who did not link to the 2011 to 2013 National Health Service (NHS) Patient Registers. Whilst this does not affect the estimated association between ethnicity and COVID-19 mortality in our sample, it may affect the external validity of our study. Since the rate of linkage failure was higher for ethnic minority groups (with the highest rates of failure for the ‘other’, mixed and Chinese ethnic groups), the most likely result is an underestimation of the association between ethnicity and COVID-19 mortality in the population living in England at the beginning of the rollout of the COVID-19 vaccination programme, because a larger proportion of deaths occurring in ethnic minority groups will have been excluded. We sought to remove people who had emigrated since the 2011 Census from the population at-risk by restricting our study population to people registered with the NHS at the start of the pandemic. However, the study population also did not include people who immigrated or were born between 2011 and 2020. As a result, it did not fully represent the population at risk. Rates of immigration are higher for ethnic minority groups [40], which may mean that the denominator in our study population was less representative for some groups. In addition, deaths occurring among recent migrants would not have been included in our analyses. However, migrants overall tend to be young and the risk of COVID-19 mortality has been shown to be lower for young people [41]. Finally, single imputation using nearest-neighbour donor input for missing Census data may have inflated the statistical precision of Census variables.

Our study shows that adjusting for vaccination status eliminates the remaining elevated risk of COVID-19 mortality that is not explained by geographical factors, socioeconomic and demographic characteristics, and pre-existing health conditions in the third wave for Black African and Black Caribbean groups and reduces the remaining excess risk substantially for the Pakistani group (particularly for Pakistani women) compared to the White British group. The reduction in risk was more modest in the second wave, possibly reflecting widening ethnic inequalities in vaccination rates over time [10] and lower uptake of third vaccine doses during the third wave in some ethnic groups. As of 1 December 2021, the age-standardised rate of receiving three vaccine doses was 25.8% for the Black Caribbean group and 27.8% for the Black African group, compared to 46.1% for the White British group. Our results suggest that increasing vaccination uptake in ethnic minority groups could substantially reduce the inequalities in COVID-19 mortality. Strategies to increase vaccine uptake include improving trust in the efficacy and safety of vaccines by engaging with local community networks with lower vaccine uptake using appropriately tailored communication strategies. Identifying barriers to access, providing vaccinations in more convenient settings, and improving vaccine uptake in deprived areas could also increase vaccination rates in ethnic groups with low uptake [42, 43]. Future studies could also explore the potential impact that prioritisation for vaccination based on age, health conditions, and occupation may have had on long-term health outcomes (e.g. post-acute COVID-19 syndrome) in ethnic minority groups.

Conclusions

The elevated rate of death involving COVID-19 observed among ethnic minority groups in the first phase of the pandemic has continued throughout the rollout of the vaccine programme until late 2021. Whilst much of these differences can be explained by geography, socioeconomic factors, and pre-existing health conditions, differences in COVID-19 vaccine coverage are also a key driver of the elevated risk of COVID-19 death, particularly among the Black African, Black Caribbean, and Pakistani groups. Increasing vaccine uptake in ethnic minority groups would help reduce the inequalities in COVID-19 mortality.

Availability of data and materials

Deidentified participant data and a data dictionary are available via the Secure Research Service for Accredited Researchers. For information on the Secure Research Service for Accredited Researchers see: https://www.ons.gov.uk/aboutus/whatwedo/statistics/requestingstatistics/approvedresearcherscheme. Analytical code is available on request to the corresponding author. As part of the Office for National Statistics wider research strategy, we are working towards making analytical code openly available.

Abbreviations

ASMR:

Age-standardised mortality rate

BMI:

Body mass index

CI:

Confidence interval

GDPPR:

General Practice Extraction Service Data for Pandemic Planning and Research

GPES:

General Practice Extraction Service

HES:

Hospital Episode Statistics

HR:

Hazard ratio

IMD:

Index of Multiple Deprivation

NHS:

National Health Service

NIMS:

National Immunisation Management Service

NS-SEC:

National Statistics Socio-economic Classification

ONS:

Office for National Statistics

PHDA:

Public Health Data Asset

SNOMED-CT:

Systematized Nomenclature of Medicine Clinical Terms

SOC:

Standard Occupational Classification

References

  1. Sze S, Pan D, Nevill CR, Gray LJ, Martin CA, Nazareth J, et al. Ethnicity and clinical outcomes in COVID-19: a systematic review and meta-analysis. EClinicalMedicine. 2020;29:100630.

    Article  Google Scholar 

  2. Magesh S, John D, Li WT, Li Y, Mattingly-app A, Jain S, et al. Disparities in COVID-19 outcomes by race, ethnicity, and socioeconomic status. JAMA Netw Open. 2021;4(11):e2134147.

    Article  Google Scholar 

  3. Gold JAW, Rossen LM, Ahmad FB, Sutton P, Li Z, Salvatore PP, et al. Race, ethnicity, and age trends in persons who died from COVID-19 - United States, May-August 2020. MMWR Morb Mortal Wkly Rep. 2020;69(42):1517–21.

    Article  CAS  Google Scholar 

  4. Khanijahani A, Iezadi S, Gholipour K, Azami-Aghdash S, Naghibi D. A systematic review of racial/ethnic and socioeconomic disparities in COVID-19. Int J Equity Health. 2021;20:248.

    Article  Google Scholar 

  5. Mathur R, Rentsch CT, Morton CE, Hulme WJ, Schultze A, MacKenna B, 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:1711–24.

    Article  CAS  Google Scholar 

  6. Nafilyan V, Islam N, Mathur R, Ayoubkhani D, Banerjee A, Glickman M, et al. Ethnic differences in COVID-19 mortality during the first two waves of the Coronavirus Pandemic: a nationwide cohort study of 29 million adults in England. Eur J Epidemiol. 2021;36(6):605–17.

    Article  CAS  Google Scholar 

  7. Office for National Statistics. Updating ethnic contrasts in deaths involving the coronavirus (COVID-19), England: 24 January 2020 to 31 March 2021. 2021. https://www.ons.gov.uk/peoplepopulationandcommunity/birthsdeathsandmarriages/deaths/articles/updatingethniccontrastsindeathsinvolvingthecoronaviruscovid19englandandwales/24january2020to31march2021. Accessed 25 Jan 2022.

  8. Department of Health and Social Care. Priority groups for coronavirus (COVID-19) vaccination: advice from the JCVI, 30 December 2020. 2020. https://www.gov.uk/government/publications/priority-groups-for-coronavirus-covid-19-vaccination-advice-from-the-jcvi-30-december-2020. Accessed 25 Jan 2022.

  9. Office for National Statistics. Coronavirus and vaccine hesitancy, Great Britain: 9 August 2021. 2021. https://www.ons.gov.uk/peoplepopulationandcommunity/healthandsocialcare/healthandwellbeing/bulletins/coronavirusandvaccinehesitancygreatbritain/9august2021. Accessed 25 Jan 2022.

  10. Dolby T, Finning K, Baker A, Dowd L, Khunti K, Razieh C, et al. Monitoring sociodemographic inequality in COVID-19 vaccination coverage in England: a national linked data study. J Epidemiol Community Health. 2022;76:646–52.

    Article  Google Scholar 

  11. Nafilyan V, Dolby T, Razieh C, Gaughan CH, Morgan J, Ayoubkhani D, et al. Sociodemographic inequality in COVID-19 vaccination coverage among elderly adults in England: a national linked data study. BMJ Open. 2021;11(7):e053402.

    Article  Google Scholar 

  12. Office for National Statistics. Coronavirus and vaccination rates in people aged 50 years and over by socio-demographic characteristic, England: 8 December 2020 to 12 December 2021. 2021. https://www.ons.gov.uk/peoplepopulationandcommunity/healthandsocialcare/healthinequalities/datasets/coronavirusandvaccinationratesinpeopleaged50yearsandoverbysociodemographiccharacteristicengland. Accessed 25 Jan 2022.

  13. Ayoubkhani D, Nafilyan V, White C, Goldblatt P, Gaughan C, Blackwell L, et al. Ethnic-minority groups in England and Wales—factors associated with the size and timing of elevated COVID-19 mortality: a retrospective cohort study linking census and death records. Int J Epidemiol. 2021;49(6):1951–62.

    Article  Google Scholar 

  14. NHS Digital. General Practice Extraction Service (GPES) Data for pandemic planning and research: a guide for analysts. 2022. https://digital.nhs.uk/coronavirus/gpes-data-for-pandemic-planning-and-research/guide-for-analysts-and-users-of-the-data#data-coverage-management-information. Accessed 1 Aug 2022.

  15. Office for National Statistics. Estimates of the population for the UK, England and Wales, Scotland and Northern Ireland: Mid-2020 edition. 2021. https://www.ons.gov.uk/peoplepopulationandcommunity/populationandmigration/populationestimates/datasets/populationestimatesforukenglandandwalesscotlandandnorthernireland. Accessed 05 Apr 2022.

  16. Office for National Statistics. Deaths registered weekly in England and Wales, 2020 edition. 2020. https://www.ons.gov.uk/peoplepopulationandcommunity/birthsdeathsandmarriages/deaths/datasets/weeklyprovisionalfiguresondeathsregisteredinenglandandwales. Accessed 25 Jan 2022.

  17. Office for National Statistics. Coronavirus and key workers in the UK. 2020. https://www.ons.gov.uk/employmentandlabourmarket/peopleinwork/earningsandworkinghours/articles/coronavirusandkeyworkersintheuk/2020-05-15. Accessed 26 Aug 2022.

  18. Hippisley-Cox J, Coupland CA, Mehta N, Keogh RH, Diaz-Ordaz K, Khunti K, et al. Risk prediction of COVID-19 related death and hospital admission in adults after COVID-19 vaccination: national prospective cohort study. BMJ. 2021;374:n2244.

    Article  Google Scholar 

  19. Nafilyan V, Humberstone B, Mehta N, Diamond I, Coupland C, Lorenzi L, et al. An external validation of the QCovid risk prediction algorithm for risk of mortality from COVID-19 in adults: a national validation cohort study in England. Lancet Digit Health. 2021;3(7):e425–33.

    Article  Google Scholar 

  20. Simpson CR, Robertson C, Kerr S, Shi T, Vasileiou E, Moore E, et al. External validation of the QCovid risk prediction algorithm for risk of COVID-19 hospitalisation and mortality in adults: national validation cohort study in Scotland. Thorax. 2022;77:497–504.

    Article  Google Scholar 

  21. Lyons J, Nafilyan V, Akbari A, Davies G, Griffiths R, Harrison E, et al. Validating the QCOVID risk prediction algorithm for risk of mortality from COVID-19 in the adult population in Wales, UK. Int J Popul Data Sci. 2022;5(4):1697.

    Google Scholar 

  22. Office for National Statistics. Which occupations have the highest potential exposure to the coronavirus (COVID-19)? 2020. https://www.ons.gov.uk/employmentandlabourmarket/peopleinwork/employmentandemployeetypes/articles/whichoccupationshavethehighestpotentialexposuretothecoronaviruscovid19/2020-05-11. Accessed 25 Jan 2022.

  23. Pace M, Lanzieri G, Glickman M, Grande E, Zupanic T, Wojtyniak B, et al. Revision of the European standard population: report of Eurostat’s task force. Luxembourg: Publication's Office of the European Union; 2013.

    Google Scholar 

  24. Office for National Statistics. Coronavirus (COVID-19) Infection Survey technical article: waves and lags of COVID-19 in England. 2021. https://www.ons.gov.uk/peoplepopulationandcommunity/healthandsocialcare/conditionsanddiseases/articles/coronaviruscovid19infectionsurveytechnicalarticle/wavesandlagsofcovid19inenglandjune2021. Accessed 07 Feb 2022.

  25. Office for National Statistics. 2011 Census item edit and imputation process. 2012. https://webarchive.nationalarchives.gov.uk/ukgwa/20160105160709/http:/www.ons.gov.uk/ons/guide-method/census/2011/census-data/2011-census-user-guide/quality-and-methods/quality/quality-measures/response-and-imputation-rates/item-edit-and-imputation-proces. Accessed 25 Jan 2022.

  26. Thernau T. A package for survival analysis in S. R package version 2.41-3. 2017. https://CRAN.R-project.org/package=survival. Accessed 25 Jan 2022.

  27. Larsen T, Bosworth ML, Ayoubkhani D, Schofield R, Ali R, Khunti K, et al. Disparities in SARS-CoV-2 case rates by ethnicity, religion, measures of socio-economic position, English proficiency, and self-reported disability: cohort study of 39 million people in England during the Alpha and Delta waves. medRxiv. 2022. https://doi.org/10.1101/2022.03.02.22271762.

  28. Yao L, Aleya L, Howard SC, Cao Y, Wang C-Y, Day SW, et al. Variations of COVID-19 mortality are affected by economic disparities across countries. Sci Total Environ. 2022;832:154770.

    Article  CAS  Google Scholar 

  29. Niedzwiedz CL, O’Donnell CA, Jani BD, Demou E, Ho FK, Celis-Morales C, et al. Ethnic and socioeconomic differences in SARS-CoV-2 infection: prospective cohort study using UK Biobank. BMC Med. 2020;18:160.

    Article  CAS  Google Scholar 

  30. Office for National Statistics. Coronavirus (COVID-19) case rates by socio-demographic characteristics, England: 1 September 2020 to 25 July 2021. 2021. https://www.ons.gov.uk/peoplepopulationandcommunity/healthandsocialcare/conditionsanddiseases/bulletins/coronaviruscovid19caseratesbysociodemographiccharacteristicsengland/1september2020to25july2021. Accessed 25 Jan 2022.

  31. Office for Health Improvements and Disparities. COVID-19 Health Inequalities Monitoring for England (CHIME) tool. 2021. https://analytics.phe.gov.uk/apps/chime/. Accessed 25 Jan 2022.

  32. Race Disparity Unit. Final report on progress to address COVID-19 health inequalities. 2021. https://www.gov.uk/government/publications/final-report-on-progress-to-address-covid-19-health-inequalities. Accessed 07 Feb 2022.

  33. Office for National Statistics. Why have Black and South Asian people been hit hardest by COVID-19? 2020. https://www.ons.gov.uk/peoplepopulationandcommunity/healthandsocialcare/conditionsanddiseases/articles/whyhaveblackandsouthasianpeoplebeenhithardestbycovid19/2020-12-14. Accessed 07 Feb 2022.

  34. Aldridge RW, Pineo H, Fragaszy E, Eyre MT, Kovar J, Nguyen V, et al. Household overcrowding and risk of SARS-CoV-2: analysis of the Virus Watch prospective community cohort study in England and Wales. Wellcome Open Res. 2021;6:347.

    Article  Google Scholar 

  35. Pritchard E, Jones J, Vihta K, Stoesser N, Matthews PC, Eyre DW, et al. Monitoring populations at increased risk for SARS-CoV-2 infection in the community using population-level demographic and behavioural surveillance. Lancet Reg Health Eur. 2022;13:100282.

    Article  Google Scholar 

  36. Nafilyan V, Islam N, Ayoubkhani D, Gilles C, Katikireddi SV, Mathur R, et al. Ethnicity, household composition and COVID-19 mortality: a national linked data study. J R Soc Med. 2021;114(4):182–211.

    Article  Google Scholar 

  37. Selden TM, Berdahl TA. COVID-19 and racial/ethnic disparities in health risk, employment, and household composition. Health Aff (Millwood). 2020;39(9):1624–32.

    Article  Google Scholar 

  38. Nafilyan V, Pawelek P, Ayoubkhani D, Rhodes S, Pembrey L, Matz M, et al. Occupation and COVID-19 mortality in England: a national linked data study of 14.3 million adults. Occup Environ Med. 2022;79:433–41.

    Article  Google Scholar 

  39. Office for National Statistics. 2011 Census Statistics for England and Wales: March 2011 QMI. 2013. https://www.ons.gov.uk/peoplepopulationandcommunity/populationandmigration/populationestimates/methodologies/2011censusstatisticsforenglandandwalesmarch2011qmi. Accessed 07 Feb 2022.

  40. Office for National Statistics. Migration by ethnic group by age. 2015. https://www.nomisweb.co.uk/census/2011/ukmig003. Accessed 05 Apr 2022.

  41. Williamson EJ, Walker AJ, Bhaskaran K, Bacon S, Bates C, Morton CE, et al. Factors associated with COVID-19-related death using OpenSAFELY. Nature. 2020;584(7821):430–6.

    Article  CAS  Google Scholar 

  42. Ethnicity sub-group of the Scientific Advisory Group for Emergencies (SAGE). Factors influencing COVID-19 vaccine uptake among minority ethnic groups. 2021. https://www.gov.uk/government/publications/factors-influencing-covid-19-vaccine-uptake-among-minority-ethnic-groups-17-december-2020. Accessed 05 Apr 2022.

  43. Schmidt H, Weintraub R, Williams MA, Miller K, Buttenheim A, Sadecki E, et al. Equitable allocation of COVID-19 vaccines in the United States. Nat Med. 2021;27(7):1298–307.

    Article  CAS  Google Scholar 

Download references

Acknowledgments

None.

Funding

The study received no external funding.

Author information

Authors and Affiliations

Authors

Contributions

MLB and VN conceptualised and designed the study. JM, MLB and TA prepared the study data. MLB and TA performed the statistical analysis. MLB, TL, and LL quality assured the underlying data and results. MLB, TA, TL, LL, JM, RA, PG, NI, KK, VR, DA, NB, MG, and VN contributed to the interpretation of the findings. MLB and VN wrote the original draft. MLB, TA, TL, LL, JM, RA, PG, NI, KK, VR, DA, NB, MG, and VN contributed to the review and editing of the manuscript and approved the final manuscript.

Corresponding author

Correspondence to Vahé Nafilyan.

Ethics declarations

Ethics approval and consent to participate

Ethical approval was obtained from the National Statistician’s Data Ethics Advisory Committee (NSDEC(20)12). This study involving secondary use of administrative datasets. Therefore, informed consent was not required.

Consent for publication

Not applicable.

Competing interests

All authors have completed the ICMJE uniform disclosure form at www.icmje.org/coi_disclosure.pdf and declare: no support from any organisation for the submitted work, no financial relationships with any organisations that might have an interest in the submitted work in the previous 3 years; KK is a member of the UK Scientific Advisory Group for Emergencies (SAGE) and chair of the ethnicity subgroup of SAGE.

Additional information

Publisher’s Note

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

Supplementary Information

Additional file 1: Fig. S1.

Percentage of people with different postcodes in GP records and 2011 Census by age group. Fig. S2. Directed acyclic graph of the hypothesised relationship between ethnicity and COVID-19 mortality. Table S1. Rates and odds ratios of linkage failure between the 2011 Census and 2011 to 2013 NHS Patient Registers for people living in England at the 2011 Census who were aged 30-100 years in 2020, by sex, age group, ethnicity, region and IMD decile. Table S2. Coding and source of variables included in the analysis. Table S3. Weighted mean follow-up time in days by ethnic groups for wave two and wave three. Table S4. Age-standardised mortality rates (ASMRs) of death involving COVID-19 (8 December 2020 to 1 December 2021) for disaggregated ‘White other’ and ‘Other’ ethnic groups.

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

Bosworth, M.L., Ahmed, T., Larsen, T. et al. Ethnic differences in COVID-19 mortality in the second and third waves of the pandemic in England during the vaccine rollout: a retrospective, population-based cohort study. BMC Med 21, 13 (2023). https://doi.org/10.1186/s12916-022-02704-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12916-022-02704-7

Keywords