INTRODUCTION
Evidence has been published since the 2009 influenza pandemic which indicates that ethnic minority populations experienced worse health outcomes following pandemic influenza infection compared to other groups [1, Reference Sachedina and Donaldson2]. This has also been observed during previous influenza pandemics [Reference Groom3–Reference La Ruche5]. The higher impact of pandemic influenza in minority populations has been related to a range of potential factors including higher rates of underlying chronic disease (increasing risk of influenza-related complications [Reference Hutchins4], in particular secondary bacterial infections) and barriers to accessing healthcare. The explanations underlying the latter have been related to socioeconomic status and cultural, educational and language factors [Reference Groom3, Reference Hutchins4].
The 2009 H1N1 influenza pandemic is now recognized to have been due to a relatively mild novel influenza virus compared to historical influenza pandemics in the 20th century, in particular the 1918 Spanish influenza pandemic [Reference Bautista6, Reference McLean7]. However, severe and fatal cases did still occur due to infection with the 2009 pandemic virus [influenza A(H1N1)pdm09]; including in previously healthy individuals [Reference Sachedina and Donaldson2, Reference Pebody8, Reference Donaldson9]. Higher rates of influenza A(H1N1)pdm09-associated morbidity and mortality were noted particularly in various ethnic minority or indigenous populations including in the USA [10], New Zealand [Reference La Ruche5, Reference Bandaranayake11], Australia [Reference La Ruche5], Canada [Reference La Ruche5, Reference Zarychanski12] and in children in England [Reference Sachedina and Donaldson2]. The latter study [Reference Sachedina and Donaldson2] highlighted that children in particular ethnic groups were at elevated risk of death. Deprivation has also been identified as a risk factor for fatal outcome during the 2009 influenza pandemic [Reference Rutter13]. The first post-pandemic influenza season in England, 2010/2011, was dominated by circulation of A(H1N1)pdm09, with intense activity mainly in young adults with many deaths and ICU admissions reported [14–Reference Green16]. The role of ethnicity and deprivation as risk factors for severe disease during this first post-pandemic season has not been fully explored.
The present study aims to provide a more detailed account of the relationship between ethnicity, deprivation and influenza mortality in England during the 2009/2010 pandemic period and to extend the analysis to the immediate post-pandemic season in 2010/2011, which was dominated by circulation of influenza A(H1N1)pdm09, and to determine if there were any differences between these two time periods.
METHODS
Study population and case definition
Individual-level mortality surveillance was established from the start of the 2009/2010 pandemic to ascertain, confirm and report fatal cases related to influenza A(H1N1)pdm09 infection in England [Reference Pebody8]. A confirmed fatal case was defined as a resident in England who died between 27 April 2009 (the date of the first pandemic cases reported from the UK [17]) and 30 April 2010, and between 1 September 2010 and 31 August 2011, where influenza A(H1N1)pdm09 infection was laboratory-confirmed and/or recorded on the death certificate as described previously [Reference Pebody8]. As ethnicity status was not recorded for fatal cases in this study, patients were classified as ‘from an ethnic minority (i.e. non-White), or ‘non-ethnic’ (i.e. White) based on available surname and first name using Onomap software [Reference Lakha18] (http://www.onomap.org/Index.aspx). The Onomap software assigned each study subject into one of the UK 2001 census ethnic groups as follows (a) White – British; (b) White – Irish; (c) White – any other White background; (d) Asian or Asian British – Indian; (e) Asian or Asian British – Pakistani; (f) Asian or Asian British – Bangladeshi; (g) Asian or Asian British – any other Asian background; (h) Black or Black British – African; (i) Chinese; (j) Other ethnic group – any other ethnic group and (k) Unclassified. The patients were then grouped into White for groups (a)–(c), non-White for groups (d)–(j) and Unknown for group (k) patients.
Using patients' home postcode, place of residence was classified into rural/suburban and urban areas as defined by the Office for National Statistics (ONS) [19]. The deprivation status of their geographical area of place of residence was defined using the Index of Multiple Deprivation (IMD 2007) [Reference Noble20]. The deprivation score can be obtained for each postcode by data linkage between IMD 2007 and postcodes. The IMD 2007 consists of seven dimensions of deprivation (income; employment; health deprivation and disability; education, skills and training; barriers to housing and services; crime and disorder; living environment). These are weighted and combined to create an overall index of multiple deprivation.
Additional information was available on a subset of cases on underlying medical condition, antiviral use and influenza vaccine uptake through surveys to patients' family physicians. Due to incomplete data reporting and missing data, these variables were only analysed using univariable analysis (χ 2 test).
Population data
Population denominator data were obtained from ONS (mid-2009 population) for England and used to calculate cumulative mortality rates. England was divided into the nine government office regions at that time. Population data by 362 local authorities (LAs), ethnicity, sex and age group (0–15, 16–64, ⩾65 years) were used in conjunction with data for the whole of the England by ethnicity and age in 5-year age groups. These were combined on the assumption that the age distributions of the 5-year age groupings were constant within the broader groups across LAs. The resulting data were aggregated to regions to provide a region, age, gender and ethnicity-specific denominator for analysis. IMD and residence area type were available at a much finer geographical level, the lower super output area (LSOA), of which there are 32 378 areas in England to enable analysis of the relationship with deprivation.
Analysis
Poisson regression models were used to estimate the mortality risk in relation to ethnicity, age, sex, region, socioeconomic status (using IMD 2007 in quartiles), and residence area type (rural/suburban vs. urban). As no data source exists that provides a full population breakdown according to all of the variables of interest, analyses were conducted using two models based on two sets of variables, for which denominator data were available. These were: model A (ethnicity, age, sex and region, with interaction terms between age and ethnicity); and model B (IMD, residence area type, and region). IMD must be analysed at a fine geographical level due to high local variability. Results are presented as adjusted incident rate ratios (IRRs) with 95% confidence intervals (CIs) from Poisson regression analysis.
RESULTS
A total of 389 fatal cases associated with influenza A(H1N1)pdm09 infection were reported between 1 June 2009 and 30 April 2010 in England. Two hundred and ninety-one (74·8%) were from both laboratory-confirmed cases and death certificates; 51 (13·1%) were from laboratory-confirmed cases only and 47 (12·1%) from death certificates only. Out of these 389 influenza deaths, 211 (54·2%) were male. For 348 (89·4%) of these deaths full name information was available. Of these, for 342 cases ethnicity information could be derived using Onomap software, with ethnicity status unclassified for six deaths, and for five cases other information was not available. Therefore 337 cases were included in the analysis. Sixty-seven (19·9%) of these 337 cases were classified as being from an ethnic minority (i.e. non-White), a significantly higher proportion than observed for the whole population in England (12·5%) (P < 0·0001).
Table 1a shows the crude and adjusted cumulative mortality rates by ethnicity from model A for the 2009/2010 pandemic period. We found a significant interaction effect between age group and ethnicity, but not for other covariates; the Pearson goodness-of-fit test for the final model provided a P value of 0·889 (no evidence of lack of fit). Within the White population, adjusting for gender and geography, there was only a significantly increased mortality risk in those aged 45–64 years compared to those aged 25–44 years in 2009/2010, with a significantly decreased mortality risk in those aged 5–14 years (Table 1a ). The overall adjusted mortality risk of non-White vs. White ethnicity in 2009/2010 was significantly higher with an IRR of 1·84 (95% CI 1·39-2·45). Stratifying by age group, non-White populations had significantly increased risk for those aged 1–4, 5–14, 45–64 and ⩾65 years compared to those aged 25–44 years in the White population in 2009/2010. There was no evidence of significant differences in mortality by region or gender in 2009/2010 after adjusting for ethnicity and age.
Rates of influenza A(H1N1)pdm09-associated mortality per million according to ethnicity and age, sex and region during (a) 2009/2010 and (b) 2010/2011 seasons in England
IRR, Incident rate ratio; CI, confidence interval; ref., reference category.
During the 2010/2011 season there were 463 deaths associated with influenza A(H1N1)pdm09 infection in England; 291 (95·5%) were from both laboratory-confirmed cases and death certificates; one (0·2%) was from a laboratory-confirmed case only and 20 (4·3%) from death certificates only. For 15 cases information on name was incomplete. Therefore, in total, 448 cases were used in the analyses, giving a crude mortality rate of 8·7/million in 2010/2011 compared to 6·5/million in 2009/2010. Forty-nine (10·9%) of these fatal cases were classified as from an ethnic minority, non-statistically lower than the national average (P = 0·195).
The adjusted cumulative mortality risk changed significantly between the two periods, with significantly higher mortality in 2010/2011 compared to 2009/2010 (IRR 1·48, 95% CI 1·27-1·72, adjusted for age, region and gender) for White populations and a non-significant decrease in mortality risk for non-White populations (adjusted IRR 0·73, 95% CI 0·51-1·06).
Table 1b shows the results from model A for 2010/2011. The adjusted mortality risk in the White population was again significantly elevated in those aged 45–64 years compared to those aged 25–44 years. However, there were significant decreases in risk for those aged 1–4, 5–14 and 15–24 years compared to those aged 25–44 years in 2010–2011. The risk profile changed in 2010/2011 with no significantly increased risk found in non-White ethnic populations by age compared to those aged 25–44 years of White ethnicity, except a borderline significant increase in those aged 1–4 years. Unlike 2009/2010, there was evidence of significant regional differences in 2010/2011, with higher mortality rates in the northern and central parts of the country compared to South-East England. The Pearson goodness-of-fit test for the final model provided a P value of 0·998.
We also examined mortality risk by individual ethnic groups. Denominator data for these groups are only available from ONS using broader age groups of <15, 15–59 and ⩾60 years for females and <15, 15–64 and ⩾65 years for males. Results are shown in Table 2(a, b) for the two periods. Crude mortality rates were significantly higher for those of Pakistani and Indian ethnicity compared to the White population in 2009/2010, and after adjusting for region and gender, the mortality risk was still significantly higher in Pakistani and Indian groups with the greatest risk observed in those of Pakistani ethnicity (Table 2a ). In general, there was little power to test for age differences across individual ethnicities (interactions), but there were significantly higher rates in Pakistanis aged <15 years in 2009/2010 (P = 0·001) with 55% in this age group; and significantly higher rates in Indians aged >60 years in 2010/2011 (P = 0·018) with 31% in this age group; otherwise, the bulk of cases in non-White groups (~61%) was in the 15–59 years age group for non-White ethnicities. Overall, the adjusted mortality risk was not significantly different during the 2010/2011 season in different ethnic groups compared to those of White ethnicity (Table 2b ). A non-significantly decreased mortality risk was seen in those of Black ethnicity compared to the White population in both periods.
Univariable and multivariable results between influenza A(H1N1)pdm09-associated mortality risk and age, sex and ethnicity group during the (a) 2009/2010 and (b) 2010/2011 seasons in England
IRR, Incident rate ratio; CI, confidence interval; ref., reference category.
* Adjusted by region and gender.
Table 3(a, b) presents the results by IMD (quartiles) and area type (rural/suburban vs. urban). Mortality risk adjusted for age group and gender was significantly increased in patients living in areas with higher levels of deprivation in both 2009/2010 and 2010/2011, with a slightly higher risk in IMD quartile 3 and a significant, twofold increase for IMD quartile 4 (most deprived) compared to IMD quartile 1 (least deprived) for both periods. The effect of IMD was found to have a linear relationship with risk of death with an overall IRR of 1·29 (95% CI 1·20-1·38) per quartile increase in IMD. Patients living in rural or suburban areas had a lower risk of death than their urban counterparts for both study periods, although this difference was not significant after adjusting for deprivation levels.
Modelling results between influenza A(H1N1)pdm09-associated mortality risk and Index of Multiple Deprivation (IMD) quartiles and rural/suburban and urban area type during the (a) 2009/2010 and (b) 2010/2011 seasons in England
IRR, Incident rate ratio; CI, confidence interval; ref., reference category.
* Adjusted for age and gender.
During the 2009/2010 pandemic period, no difference was found in the proportion of fatal cases reported to have any underlying medical conditions between White and non-White patients (74·3%, 182/245 vs. 74·3%, 52/70, P = 1·0). Available data from the 2010/2011 period showed no significant difference in proportion of cases with any underlying medical conditions between White (64·9%, 237/365) and non-White cases (61·5%, 32/52, P = 0·632).
During the 2009/2010 pandemic period, data for reported influenza antiviral use between White and non-White cases did not show a significant difference (65·2%, 45/69 vs. 75·0%, 12/16, P = 0·453). In 2010/2011, there was no significant difference in the proportion of cases who received antivirals between White and non-White patients (33·6%, 37/110 vs. 30·8%, 4/13, P = 0·836). No influenza vaccine uptake data were available for analysis during the 2009/2010 period, and no significant difference between White and non-White cases was found in influenza vaccine uptake during 2010/2011 (20·5%, 45/220 vs. 25·0%, 6/24, P = 0·603).
DISCUSSION
Our study has shown that non-White ethnic minorities, in particular those of south Asian descent, experienced a significantly increased risk of fatal outcome following influenza A(H1N1)pdm09 infection during the 2009/2010 influenza pandemic period in England compared to the White population. This excess risk appeared to be focused in the youngest and oldest age groups; but there were no significant differences according to gender or region in non-White ethnicities. This elevated risk of death in those from ethnic subgroups was no longer observed in the immediate post-pandemic 2010/2011 season, although the overall mortality risk in the general population was higher in 2010/2011 compared to that seen in 2009/2010. In addition, in both periods, patients living in areas with a higher level of deprivation had an increased risk of mortality associated with influenza A(H1N1)pdm09 infection compared to those living in the most affluent areas.
Our findings are of public health importance in that the novel pandemic influenza A(H1N1)pdm09 virus was associated with increased impact in non-White ethnic populations, particularly those of Asian heritage, during the 2009/2010 pandemic period compared to the White population. No such increased relative impact was found during the following non-pandemic influenza season, when the same virus strain continued to circulate with even higher population impact. Available historical data suggest that an increased risk of adverse health outcome including higher mortality is likely to occur in specific ethnic subgroups during an influenza pandemic [Reference Hutchins4]. Earlier published data indeed showed the disproportionately higher impact on ethnic minority populations during the 2009/2010 pandemic [Reference Sachedina and Donaldson2, Reference La Ruche5, 10]. Previous studies on the impact of influenza have focused primarily on the impact during a pandemic period [1–Reference La Ruche5] but not during a post-pandemic period. The explanations for our findings of changes in risk over time are likely to be multi-factorial and complex. The elevated impact of pandemic influenza in specific ethnic sub-populations could be due to either increased rates of infection and/or infection-severity ratios compared to the White population. Disentangling these two components is challenging. Importation and initial spread of a novel respiratory virus depends upon a range of factors including migration and mixing patterns. The roles of genetic factors on human susceptibility to influenza have long been suspected [21] but conclusive evidence is not available [Reference Horby22]. In terms of case severity, a recent study seemed to suggest that a certain gene might be linked to increased severity of influenza in Chinese populations during the 2009/2010 pandemic period when comparing severe with mild cases [Reference Zhang23].
Other explanatory factors for elevated case severity may relate to poor underlying health status and access to preventive interventions. However, the proportion of cases with pre-existing underlying medical conditions during the 2009/2010 pandemic was similar for White and non-White patients. Available data did not show any significant differences in influenza vaccine uptake between the White and non-White cases. Although, data on influenza vaccination coverage among different ethnic populations have not been collected in England; an earlier study on pneumococcal vaccine uptake in England and Wales [Reference Pebody24] did show that the uptake rates for pneumococcal vaccine for the total population are substantially higher in areas where >99% of the population are White, which might indicate a similar situation for other vaccines, such as that for influenza. This highlights the importance of gathering influenza vaccine uptake data by ethnicity, which is now being implemented into the UK.
Differences in healthcare-seeking behaviour might be a possible contributory factor. The disproportionately poorer health outcomes experienced by ethnic minority populations during the pandemic period may be associated with limited access to healthcare facilities, and greater cultural, educational and linguistic obstacles in adopting pandemic influenza control measures, unfavourable environmental factors (such as household crowding and poor quality housing), behavioural differences in response to influenza, [Reference La Ruche5, Reference Bandaranayake25]. Available preliminary data, however, did not show any significant differences in influenza antiviral use between the White and non-White patients. Further work is required to better understand any differences.
Finally the strong association we found between higher deprivation and increased risk for severe disease might play a role in the higher mortality risk observed in the non-White ethnic minorities. An earlier study [Reference Rutter13] showed that the mortality risk significantly increased in the deprived areas during the 2009/2010 influenza pandemic. Our results demonstrate that this increased risk was maintained into the post-pandemic period in 2010/2011 influenza season. This indicates that deprived areas are more vulnerable to the impact of influenza infection. Socioeconomic deprivation may result in poorer health outcomes because of a range of potential explanatory factors, such as poorer access to or use of healthcare services, lower vaccine uptake, overcrowding and poorer nutrition.
Some assumptions had to be made in these analyses. Population data were not available for the full age, sex and ethnicity breakdown, and we have assumed that age distributions within the broad age groups available for LA data are constant. However, this should not impact adversely on our estimates of ethnicity and age-specific risks, as violation of this assumption would principally affect the estimation of regional differences. Our attempts to estimate all factors within one model were hampered by different levels of data granularity, and required us to examine risk factors within two different models. IMD in particular must be examined at a very local level as neighbourhoods with highly variable deprivation levels may be closely mixed, particularly in areas such as inner London [Reference Noble20]. In the absence of a fine-level breakdown of ethnicity and deprivation, it was not possible to disentangle the effects of inherent cultural differences between ethnic groups, and other factors such as deprivation and area type, which may well be correlated [26]. A potential limitation of this study is that the IMD is classified at a geographical level so there is the possibility of ecological bias, although IMD data is at the finest level (LSOA). Using the Onomap software package to classify patients into different ethnic groups may bring in some misclassification errors, although it is considered a standardized methodology for classification of ethnicity [Reference Lakha18]. This ethnicity classification method is based on the patient's name, and if the maternal side is non-White but the paternal side is White, the name of the mixed-race status of the next generation may not be correctly classified using this method, hence a degree of caution should be exercised when interpreting the results.
Another limitation of this study is that a lack of denominator data on underlying medical conditions in the population means that we could not account for their contributions to mortality within our models; and underlying medical conditions would certainly play a major role. Nevertheless, the shift in patterns from 2009/2010 to 2010/2011 shows that there was increased risk for non-White ethnicities in the pandemic period compared to the post-pandemic period; regardless of potential differences in risk due to underlying conditions, which would have been unlikely to change in the general population. In conclusion, an increased mortality impact was associated with influenza A(H1N1)pdm09 during the 2009/2010 pandemic period in non-White ethnic minorities (especially in those of certain south Asian groups): an observation that disappeared in the first post-pandemic season. This suggests these groups were among the first affected by this novel virus, rather than necessarily inherently higher case-severity ratios and reflect initial transmission patterns. In addition elevated risk was also seen for those residing in areas with higher levels of deprivation in England in both the pandemic and first post-pandemic period. These results emphasize the importance of ensuring that influenza control and prevention measures are in place for at-risk and vulnerable populations both for seasonal influenza and when a novel influenza virus emerges. Further studies are warranted to further investigate the reasons for these differences in impact.
ACKNOWLEDGEMENTS
We thank our data providers for reporting death data to us: Office for National Statistics, Health Protection teams across England, Chief Medical Officer's Office, and the Respiratory Virus Unit at PHE Colindale. We also thank the following colleagues for their help for providing data: Pauline Kaye, Rashmi Malkani, Paragi Patel, Erol Yusuf, and Osman Dar. We thank Estelle McLean for earlier data matching and cleaning process, and Matthew Bull and Sam Lattimore for postcode and ethnicity data linkage work.
DECLARATION OF INTEREST
None.
