Changes in severity of 2009 pandemic A/H1N1 influenza in England: a Bayesian evidence synthesis
BMJ 2011; 343 doi: https://doi.org/10.1136/bmj.d5408 (Published 08 September 2011) Cite this as: BMJ 2011;343:d5408- A M Presanis, investigator statistician1,
- R G Pebody, consultant epidemiologist2,
- B J Paterson, research fellow in epidemiology23,
- B D M Tom, investigator scientist1,
- P J Birrell, career development fellow in biostatistics1,
- A Charlett, Head of statistics, modelling and economics department2,
- M Lipsitch, professor of epidemiology4,
- D De Angelis, senior statistician12
- 1Medical Research Council Biostatistics Unit, Institute of Public Health, University Forvie Site, Cambridge CB2 0SR, UK
- 2Health Protection Agency Health Protection Services, London NW9 5EQ, UK
- 3Hunter Medical Research Institute, University of Newcastle, Newcastle, Australia
- 4Centre for Communicable Disease Dynamics, Departments of Epidemiology and Immunology and Infectious Diseases, Harvard School of Public Health, 677 Huntington Avenue, Boston, MA 02115, USA
- Correspondence to: A M Presanis anne.presanis{at}mrc-bsu.cam.ac.uk
- Accepted 27 June 2011
Abstract
Objective To assess the impact of the 2009 A/H1N1 influenza pandemic in England during the two waves of activity up to end of February 2010 by estimating the probabilities of cases leading to severe events and the proportion of the population infected.
Design A Bayesian evidence synthesis of all available relevant surveillance data in England to estimate severity of the pandemic.
Data sources All available surveillance systems relevant to the pandemic 2009 A/H1N1 influenza outbreak in England from June 2009 to February 2010. Pre-existing influenza surveillance systems, including estimated numbers of symptomatic cases based on consultations to the health service for influenza-like illness and cross sectional population serological surveys, as well as systems set up in response to the pandemic, including follow-up of laboratory confirmed cases up to end of June 2009 (FF100 and Fluzone databases), retrospective and prospective follow-up of confirmed hospitalised cases, and reported deaths associated with pandemic 2009 A/H1N1 influenza.
Main outcome measures Age specific and wave specific probabilities of infection and symptomatic infection resulting in hospitalisation, intensive care admission, and death, as well as infection attack rates (both symptomatic and total). The probabilities of intensive care admission and death given hospitalisation over time are also estimated to evaluate potential changes in severity across waves.
Results In the summer wave of A/H1N1 influenza, 0.54% (95% credible interval 0.33% to 0.82%) of the estimated 606 100 (419 300 to 886 300) symptomatic cases were hospitalised, 0.05% (0.03% to 0.08%) entered intensive care, and 0.015% (0.010% to 0.022%) died. These correspond to 3200 (2300 to 4700) hospital admissions, 310 (200 to 480) intensive care admissions, and 90 (80 to 110) deaths in the summer wave. In the second wave, 0.55% (0.28% to 0.89%) of the 1 352 000 (829 900 to 2 806 000) estimated symptomatic cases were hospitalised, 0.10% (0.05% to 0.16%) were admitted to intensive care, and 0.025% (0.013% to 0.040%) died. These correspond to 7500 (5900 to 9700) hospitalisations, 1340 (1030 to 1790) admissions to intensive care, and 240 (310 to 380) deaths. Just over a third (35% (26% to 45%)) of infections were estimated to be symptomatic. The estimated probabilities of infections resulting in severe events were therefore 0.19% (0.12% to 0.29%), 0.02% (0.01% to 0.03%), and 0.005% (0.004% to 0.008%) in the summer wave for hospitalisation, intensive care admission, and death respectively. The corresponding second wave probabilities are 0.19% (0.10% to 0.32%), 0.03% (0.02% to 0.06%), and 0.009% (0.004% to 0.014%). An estimated 30% (20% to 43%) of hospitalisations were detected in surveillance systems in the summer, compared with 20% (15% to 25%) in the second wave. Across the two waves, a mid-estimate of 11.2% (7.4% to 18.9%) of the population of England were infected, rising to 29.5% (16.9% to 64.1%) in 5-14 year olds. Sensitivity analyses to the evidence included suggest this infection attack rate could be as low as 5.9% (4.2% to 8.7%) or as high as 28.4% (26.0% to 30.8%). In terms of the probability that an infection leads to death in the second wave, these correspond, respectively, to a high estimate of 0.017% (0.011% to 0.024%) and a low estimate of 0.0027% (0.0024% to 0.0031%).
Conclusions This study suggests a mild pandemic, characterised by case and infection severity ratios increasing between waves. Results suggest low ascertainment rates, highlighting the importance of systems enabling early robust estimation of severity, to inform optimal public health responses, particularly in light of the apparent resurgence of the 2009 A/H1N1 strain in the 2010-11 influenza season.
Introduction
Since the first confirmed cases of pandemic influenza A/H1N1 were reported in April 2009, 18 449 deaths have been notified worldwide.1 Reported numbers of laboratory confirmed cases and deaths are of limited value in determining the severity of a disease because of the difficulties in identifying, testing, confirming, and reporting cases, particularly during a pandemic.2 3 4 5 6 7 8 For policy makers, however, understanding severity is crucial to determine appropriate public health responses. Severity estimates are needed early on in an outbreak, when timely (but perhaps not the most robust) estimates are required, but also when an epidemic has run its course, to quantify robustly both the severity and total burden of disease and to assess whether public health responses and surveillance systems were adequate in the midst of the epidemic. An assessment of severity is a key element in understanding the epidemiology of the 2009 pandemic and for planning for future pandemics, particularly if more severe or more transmissible influenza strains emerge. Furthermore, if unusual patterns of age specific mortality continue for several years, as has happened in past pandemics9 and as seems to be the case with the 2009 pandemic,10 knowledge of the severity and burden of the strain is important for planning healthcare resource allocation as well as for understanding if severity is changing over time.
A measure of the severity of an infection is its case severity, the probability that an infected individual develops severe disease. Specifically, the case fatality ratio (the probability that an infection leads to death, estimated by taking the ratio of deaths to cases) is often used, together with the case hospitalisation ratio and case intensive care admission ratio (table 1⇓). A case may be defined either as any infection or a symptomatic infection, meaning a febrile influenza-like illness. Although influenza typically causes symptomatic infection, a substantial proportion of infected individuals will be asymptomatic or have only mild symptoms. These individuals are unlikely to be identified unless serological testing is undertaken on a population basis.11 12 Estimation is therefore required to obtain a quantification of severity. In this paper, we estimate the symptomatic case fatality, case intensive care admission, and case hospitalisation ratios (the probabilities of symptomatic infections leading to severe events (table 1⇓)) for pandemic A/H1N1 in England, from multiple sources of data available from various surveillance systems, and derive the corresponding infection severity ratios (the probabilities of all infections leading to severe events) using data from population sero-incidence surveys on the infection attack rate.13
Severity estimation is potentially complicated by two key problems affecting the observed numbers of severe events and cases—censoring and ascertainment bias.2 3 4 6 7 8 Censoring happens in the midst of an epidemic, when some severe events resulting from infections to date have yet to occur, leading to underestimation of the number of severe events. Ascertainment bias refers more generally to undercounting of cases because of surveillance systems not capturing all cases for multiple reasons. To estimate case severity ratios correctly, methods accounting for these biases are required. Previous attempts to estimate the symptomatic case fatality ratio of pandemic A/H1N1 in England3 4 5 6 7 8 9 10 11 12 13 14 15 have accounted for censoring but not all ascertainment biases and, in particular, have not accounted for all the uncertainty inherent in the data. Crucially, they have not used all the available information in the estimation process. To address these limitations, we adapted a Bayesian approach to severity estimation,8 synthesising all available relevant data and prior information on biases to derive estimates of the infection and case severity ratios, the infection attack rate, and the symptomatic attack rate in England. For a retrospective analysis, as presented here, censoring is no longer an issue, so we did not need to account for it, but our analytical approach is easily adapted to a mid-pandemic situation where censoring does induce bias.8 In accounting for ascertainment bias, we assessed the proportion of infections detected in both pre-existing surveillance systems monitoring seasonal influenza and systems specifically set up in response to the pandemic A/H1N1 outbreak.
As England experienced two waves of pandemic A/H1N1 infection in 2009-10, we investigated changes in severity across the waves, comparing the periods 1 June to 31 August 2009 (summer) and 1 September 2009 to 28 February 2010 (autumn-winter). In absolute terms, more deaths and hospitalisations were observed from September to November than in the three summer months: it is important to understand whether this reflects greater numbers of infections in the autumn, better ascertainment, or a real increase in severity. The answer has important implications for surveillance and public health responses in future pandemics.
Methods
Sources of information
Several influenza surveillance systems at the Health Protection Agency provide data on pandemic A/H1N1 cases.
Confirmed cases—The FF10016 and FluZone17 databases comprise detailed information on each of the first few thousand confirmed cases of pandemic A/H1N1 influenza, including dates of severe outcomes. From these, we estimated the delay from symptom onset to hospitalisation using a parametric mixture survival model,18 giving an estimate (assumed unbiased) of the proportion of confirmed cases hospitalised, which we used indirectly in our Bayesian analysis. Since the survival model was fitted to data from early on in the epidemic, on confirmed cases that were followed up only until the end of June 2009, this is the only analysis where we have accounted for censoring.
Symptomatic cases—During the pandemic, the Health Protection Agency fitted a regression model to data on the number of consultations for influenza-like illness (in primary care and through the National Pandemic Flu Service) and to data on the positivity rate, to produce weekly age specific and region specific estimates19 20 of the number of symptomatic cases accessing healthcare. These were divided by estimates of the proportion of influenza-like illness cases in patients who consulted healthcare services, to obtain estimates of the total number of symptomatic cases. These are thought to underestimate the true symptomatic attack rate, given other evidence on the infection attack rate.13
Seroincidence—Repeated cross sectional serological surveys were undertaken as part of the annual collection of residual sera for the Health Protection Agency sero-epidemiology programme from patients accessing healthcare in England. Data from these surveys, taken before the pandemic (a 2008 baseline), after the first wave of infections in August-September 2009,13 and after the second wave21 are published, providing evidence (though highly uncertain because of small sample sizes, low power, and potential biases in the sampled population) on the age specific infection attack rate during both waves.
Hospitalisation—A web based surveillance system in England was established by the Health Protection Agency and Department of Health after the end of the first wave in August 200922 to ascertain and collect data on all confirmed cases of pandemic A/H1N1 influenza hospitalised in acute NHS trusts in England. Of these, 970 cases, confirmed during the first wave, were retrospectively added to the database and 1306 cases were added prospectively during the second wave. Ascertainment of cases, however, was incomplete, with only 129 (77%) of the 168 acute NHS hospital trusts in England participating. This database provided information (from a line listing at 9 July 2010) on the number of hospitalisations of confirmed cases and on the ratios of intensive care admissions and deaths to hospitalisations. We considered the observed hospitalisations to reflect only a subset of all pandemic A/H1N1 hospitalisations but assumed the observed ratios of intensive care admission to hospitalisation and death to hospitalisation provided unbiased information on the probabilities of severe events given hospitalisation.18
Death—As at 9 July 2010, 380 deaths among people with confirmed pandemic A/H1N1 infection or mention of influenza on the death certificate had been reported to the Health Protection Agency or the chief medical officer.14 15 Among individuals with symptom onset from 1 June to 31 August, 79 deaths were observed in the summer wave, and 301 deaths were observed in the second wave among individuals with symptom onset 1 September to 28 February. A capture-recapture analysis suggests that deaths were under-ascertained by approximately 10%.23
Estimation approach
We adapted a Bayesian approach to estimate the attack rates, infection severity ratio, and symptomatic case severity ratio.8 We considered five severity levels (infection, symptomatic infection, hospitalisation, intensive care admission, and death (fig 1⇓)), denoted g=INF, S, H, I, and D; and seven age groups, <1, 1-4, 5-14, 15-24, 25-44, 45-64, and ≥65 years.
We defined Na,t,g to be the number of infections in age group a, wave t, and severity level g (table 2⇓), and defined ca,t,g\k to be the probability that pandemic A/H1N1 cases in age group a and wave t at severity level k progress to level g (namely, the ratio of cases at level g to those at level k if cases at g are a subset of those at k (such as if all deaths occur in hospital)). Then the infection and symptomatic case severity ratios were expressed as products of the component probabilities of progressing through each successive level of the severity pyramid (fig 1⇓, table 2⇓). For example, the case fatality ratio was defined as the product of the probability of hospitalised cases dying, the probability of symptomatic cases being hospitalised, and the probability of an infection being symptomatic: ca,t,D\INF=ca,t,D\Hca,t,H\Sca,t,S\INF. These components were estimated by synthesising the available information on infection at each severity level and on ascertainment probabilities, da,t,g, the probabilities that severe events were captured by surveillance systems (table 2⇓). Further combining with population sizes, we also obtained estimates of the infection and symptomatic attack rates.
To do this, we combined the observed surveillance data with our knowledge of each probability before observation, summarised by a “prior” distribution, to obtain an updated distribution (the “posterior”) which formally summarised our final knowledge of the quantities of interest. This distribution, accounting for both imperfect detection and all uncertainty inherent in the data as well as uncertainty about some model assumptions through the introduction of informative prior distributions,18 was summarised by its centiles—the median and the 2.5 and 97.5 centiles (denoted the 95% credible interval). The Bayesian analysis was carried out in OpenBUGS.24
Statistical model
Figure 2⇓ illustrates the relation between the data and the quantities to be estimated. Full model details are given by Presanis et al18 but are briefly described here, with table 2⇑ summarising the model parameters. The Bayesian evidence synthesis combined both direct (solid lines, fig 2⇓) and indirect information on each parameter, where indirect evidence might take the form of data on other quantities indirectly influencing the parameter of interest (such as the number of deaths informing the number of hospitalisations) or model assumptions.
Data from the sero-incidence study13 were used to inform the infection attack rate in the summer wave (ca,Summer,INF\POP). We considered the Health Protection Agency estimates of symptomatic infection19 20 to be biased downwards, providing information on a quantity (Na,t,SB) that was a proportion of the true number of symptomatic cases (Na,t,S, where da,t,S denotes this proportion). This proportion was estimated in the first wave of influenza from the information on the infection attack rate in combination with prior information on the proportion (cS\INF) of infections that are symptomatic (table 1⇑). This proportion symptomatic was assumed to be equal across age groups25 and waves. By assuming the age specific bias in the Health Protection Agency estimates, da,t,S, was similar, but not necessarily equal, in the two waves (table 1⇑), we estimated the number of symptomatic infections Na,t,S in the autumn-winter period.18 Combined again with the prior information on the proportion of infections which are symptomatic, we obtained an estimate of Na,t,INF in the autumn-winter wave.
The detected number of hospitalised cases22 and deaths14 15 were considered a subset of the true number of cases, Na.t.g, for g=H,D respectively, with corresponding detection probability, dt,g, assumed equal across age groups. Similarly, the observed number of intensive care admissions and deaths from the subset of hospitalisations with non-missing data on outcomes were a proportion of the observed number in the subset, denoted ca,t,g\H for g=I,D, respectively (table 2⇑).
Each age specific symptomatic case hospitalisation ratio, ca,t,H\S, was constrained to be less than the corresponding age specific confirmed case hospitalisation ratio. These confirmed case hospitalisation ratios were themselves informed by the age specific estimates given by a parametric mixture survival model.18 After accounting for censoring, the estimated confirmed case hospitalisation ratio averaged over age was 1.20% (95% confidence interval 0.97% to 1.49%).
Importantly, we assumed the detection probabilities to be less than 100% for both hospitalisations and deaths. We gave dt,H a vague prior distribution (table 2⇑) that accounts for under-ascertainment. We chose a prior distribution for dt,D with mean 90%, ranging between 80% and 97%, based on a capture-recapture study,23 to reflect under-ascertainment due to test sensitivity or other reasons for failing to appear in surveillance systems. We allowed both dt,H and dt,D to vary by wave, but not age.
Results
Table 3⇓ shows that in the summer wave of influenza, 0.54% (95% credible interval 0.33% to 0.82%) of the estimated 606 100 (419 300 to 886 300) symptomatic cases were hospitalised, corresponding to 540 (330 to 820) out of every 100 000 symptomatic cases and to a total of 3200 (2300 to 4700) hospitalisations. Of the symptomatic cases in this first wave, 0.05% (0.03% to 0.08%) entered intensive care (corresponding to 310 (200 to 480) intensive care admissions) and 0.015% (0.010% to 0.022%) died (90 (80 to 110) deaths).
In the second wave, 0.55% (0.28% to 0.89%) of the 1 352 000 (829 900 to 2 806 000) estimated symptomatic cases were hospitalised, corresponding to 7500 (5900 to 9700) admissions. Of the symptomatic cases in the second wave, 0.10% (0.05% to 0.16%) entered intensive care (1340 (1030 to 1790) admissions) and 0.025% (0.013% to 0.040%) died (340 (310 to 380) deaths).
Just over a third (35% (26% to 45%)) of infections were estimated to be symptomatic. The estimated probabilities of infections resulting in severe events were therefore 0.19% (0.12% to 0.29%), 0.02% (0.01% to 0.03%), and 0.005% (0.004% to 0.008%) in the summer wave for hospitalisation, intensive care, and death, respectively. The corresponding second wave proportions were 0.19% (0.10% to 0.32%), 0.03% (0.02% to 0.06%), and 0.009% (0.004% to 0.014%) (table 3⇑). Age specific infection and case severity ratios are given in figure 3⇓. The case fatality ratios, for both symptomatic cases and all cases, differed substantially across waves (table 3⇑), with 93.5% posterior probability that these were greater in the autumn-winter wave than in the summer. The age specific differences over time were less pronounced (fig 3⇓). The probabilities of intensive care admission and death among hospitalised patients also increase over time (fig 4⇓). We estimated a cumulative infection attack rate across the two waves of 11.2% (7.4% to 18.9%), with age specific and wave specific estimates given in table 4⇓.
The observed number of hospitalisations in the autumn-winter wave was 1481, representing an estimated 20% (15% to 25%) of the estimated 7500 (5900 to 9700) hospitalised cases (table 3⇑), a substantial bias. Surveillance of suspected (rather than confirmed) cases of pandemic A/H1N1 influenza indicated that 17 519 of these were hospitalised in the autumn-winter wave17: comparison with our estimate of the number of hospitalised cases (table 3⇑) therefore suggests between 34% and 55% of suspected cases were true cases. An estimated 30% (20% to 43%) of hospitalisations were detected in surveillance systems in the summer, compared with 20% (15% to 25%) in the second wave. The bias in the Health Protection Agency case estimates, averaged across age, was substantial but uncertain, with the case estimates Na,t,SB representing 40% (24% to 59%) of the true number of symptomatic cases, Na,t,S, in the summer and 36% (19% to 57%) of them in the autumn-winter.
Sensitivity analyses
Several sensitivity analyses to the choice of denominator data and to the prior distribution of the infection attack rate were performed to further assess the uncertainty in the attack rates.18 Differences due to the choice of prior distribution of the infection attack rate were small relative to the differences due to choice of denominator data, so we concentrated on sensitivity to the denominator. Specifically, four models with different denominators were considered: (1) Health Protection Agency case estimates only, assumed unbiased; (2) the model presented in the main results section, with both the case estimates (assumed biased downwards) and the summer sero-incidence data, and the second wave attack rates obtained by assuming a similar bias in the Health Protection Agency case estimates in the two waves; (3) using the sero-incidence data from both waves,13 21 assuming these are unbiased and that the Health Protection Agency case estimates are biased downwards; and (4) using the sero-incidence data from both waves,13 21 assuming these are biased upwards and the Health Protection Agency case estimates are biased downwards.
The estimates of the infection attack rate in both waves across the four different choices of denominator may be as low as 5.9% (4.2% to 8.7%) (model 1) or as high as 28.4% (26.0% to 30.8%) (model 3), with the corresponding estimates for 5-14 year olds varying between 17.9% (11.5% to 28.4%) and 58.8% (52.8% to 64.5%). The estimates of the case fatality ratio in the second wave were therefore highest in model 1 at 0.017% (0.011% to 0.024%) and lowest in model 3 at 0.0027% (0.0024% to 0.0031%).
Discussion
This study suggests a mild pandemic, with case and infection severity ratios highest in children and older adults and increasing over time. Infection attack rates were highest among school age children. The results suggest under-ascertainment of severe cases through routine systems and considerable uncertainty in the denominators of symptomatic and all infection, highlighting the limits of current surveillance. Robust systems are essential for future early estimation of severity.
Variation by age and wave of influenza
The case severity ratio estimates show a clear U shape to the age distribution (fig 3⇑), in contrast to the estimated infection and symptomatic attack rates (table 4⇑), which have an inverted U shape, with the peak in school age children. These distributions indicate that proportionately fewer adults were infected than children, with 5-14 year olds having an infection attack rate of 29.5% (16.9% to 64.1%), whereas young and old cases were most likely to be severely affected. Among those hospitalised, the age distribution was again different, with the probability of severe outcomes increasing with age (fig 4⇓), possibly due to increased risk of comorbidities.22 Although children and older adults were more likely than other adults to be admitted to hospital, once there, children were less likely than adults to enter intensive care or to die.
Our estimates of the number of hospitalisations and corresponding detection probability suggest a large bias in the observed hospitalisations due to under-ascertainment of confirmed cases. However, the estimated number was still substantially smaller than the observed suspected number of hospitalisations, suggesting between 34% and 55% of suspected hospitalisations in the autumn-winter wave were true pandemic A/H1N1 influenza cases. We estimated an increase in infection, case, and hospitalisation severity ratios between waves. The age specific increases were slightly less evident than the overall increases, because of the smaller sample sizes and hence greater uncertainty.
The higher severity among hospitalised cases older than 25 years estimated in the second wave (fig 4⇑) may indicate a change in healthcare seeking behaviour or hospital admission policy, or may represent a real increase in severity. Initially, when less was known about risk groups and outcomes, individuals may have been more likely to seek healthcare and hospitals more cautious in admitting suspected cases than during the second wave,26 27 so that those observed in hospital in the second wave may have had more severe illness on average than those hospitalised in the first wave. The difference between waves was less pronounced in children and young adults, possibly because the same cautious hospital admission behaviour may have applied to this population throughout.
Uncertainties and assumptions
Although behaviour change might explain the increase in probabilities of severe outcomes after hospitalisation, it might not be the only contributory factor to the increase in symptomatic case severity ratios. It is possible that a true increase in severity may have played a part in the estimated increase if, for example, other effect modifying factors such as bacterial super-infection were more common in the autumn-winter wave than in the summer, as was observed.28 29
If behaviour change had played a large part in the estimated increase in severity, then we might not have accounted adequately for under-ascertainment in the denominator of symptomatic cases, despite finding a substantial though uncertain bias in the Health Protection Agency estimates, or we may have overestimated the case hospitalisation ratio in the second wave. Alternatively, the observed ratios of intensive care admission to hospitalisation and death to hospitalisation ratios (assumed to provide unbiased information) may actually have been biased.
Other unaccounted factors might also have influenced apparent changes in severity, such as any potential shift in the age distribution of incident cases,30 the use of antiviral drugs as prophylaxis or treatment, and the vaccination campaign. These last two factors might be postulated to reduce rather than increase severity.17 However, vaccination uptake increased only towards the second half of the second wave,17 so it might not have had a large effect on severity. Also, the per population rate of collection of antivirals via the National Pandemic Flu Service was lower in the second wave than in the first,17 so if the attack rate was larger in the second than in the first wave, the effect of antivirals on reducing severity might have been greater in the first wave than in the second.
The estimated infection attack rate across all ages in the summer wave of 3.4% (2.4% to 4.8%), and the corresponding rate in 5-14 year olds of 10.0% (5.9% to 14.9%), are somewhat smaller than those estimated from the sero-incidence data alone in London and the West Midlands, but larger than those estimated from the sero-incidence data alone in other regions.13 Thus there may have been regional differences13 that were smoothed over by our estimates averaged across England. A regional analysis may be beneficial, although the consequent small sample sizes would entail large uncertainty in such estimates.
The differences also highlight the uncertainties surrounding the estimation of attack rates, as evidenced by the relatively wide credible intervals, particularly for the age specific rates (table 4⇑). Results from our sensitivity analyses show that the attack rate estimates are highly dependent on the model and data used, further emphasising the uncertainties surrounding the denominators. The sensitivity of results to data used to inform denominators18 suggests further work is needed in synthesising all available data. It is clear that the evidence on denominators provided by the sero-incidence data and by the Health Protection Agency estimates conflict, the extent to which either source may be biased is less apparent.
The Health Protection Agency case estimates relied on broad assumptions about the proportion of individuals with influenza-like illness who contacted healthcare services,20 without allowing for a likely substantial change in this proportion over the course of the epidemic.31 Further, the Health Protection Agency case estimates represent numbers symptomatic where symptomatic is interpreted as presentation with febrile influenza-like illness due to pandemic A/H1N1, and may therefore miss individuals with milder symptoms due to pandemic A/H1N1.20
On the other hand, as Miller and colleagues13 and Hardelid et al21 point out, their cumulative incidence estimates rely on the definition of “positivity” used, as well as on an assumption that the study population is representative of the population of interest. There is concern that the serology samples taken during the pandemic might have come from a population at greater risk of infection than average, and therefore more likely to have been vaccinated in the second wave, possibly introducing bias.13 21 32 The authors also point out that small sample sizes imply limited power to detect changes in prevalence, particularly in older age groups. Care should therefore be taken not to over-interpret the estimates from our sensitivity analyses.
The model on which our main results are based (model 2 in our sensitivity analyses) incorporates uncertainty in the estimates that approximately covers the range, by age group, of the uncertainty in the estimates from model 1 and model 3, the two models which assume unbiased Health Protection Agency case estimates and unbiased sero-incidence data, respectively. Our estimates of the proportion of cases that were symptomatic (ranging from 30% (21% to 40%) to 39% (30% to 49%) in the sensitivity analyses) are broadly comparable with estimates from studies of seasonal influenza11 12 and other studies of pandemic A/H1N1.25 33 34 35 A further sensitivity analysis18 to the choice of prior information on the proportion of infections that are symptomatic (mean prior value 63% (95% credible interval 50% to 76%)12) gives a posterior estimate of 55% (40% to 69%), with estimates of symptomatic case severity and the infection attack rate mildly sensitive to this choice. However, the difference in estimates is outweighed by the uncertainties in the denominators of symptomatic infection and all infection.
Our estimates rely on some assumptions of representativeness, such as that observations from hospital surveillance are representative of all English hospitals, not just those participating in the scheme, and that reporting practices did not change between the two waves of influenza. Finally, by using the sero-incidence data to infer the infection attack rate, we are implicitly assuming that by “infection” we mean infection with detectable antibody response, since we did not adjust for test sensitivity. Given the small sample sizes, concerns about the representativeness of the serology data, and the uncertainty in the denominators, we believe any uncertainty due to test sensitivity will be relatively small.
Comparison with other studies
Our analysis gives definitive estimates of the case severity ratios in England, lying between the US estimates based on medical attendances in Milwaukee and hospitalisations in New York (approach 1 in paper by Presanis et al8; symptomatic case fatality ratio 0.048% (0.026% to 0.096%)) and the estimates based on self reported influenza-like illness in New York (approach 2 in paper by Presanis et al8; symptomatic case fatality ratio 0.007% (0.005% to 0.009%)). Our estimates are somewhat lower than the early estimates of the confirmed case fatality ratio provided by Garske et al,3 which ranged from 0.11% to 1.47% overall and from 0.13% to 0.41% in the UK, accounting for censoring. These were based, however, only on confirmed cases as a denominator, so were expected to be larger. Nishiura et al4 likewise considered confirmed cases as a denominator, accounting for censoring, and estimating between 0.16% and 4.48% of these in the US and Canada to be fatal.
Our estimates are more comparable with other estimates where infections or symptomatic infections were used as a denominator: Wilson and Baker6 estimated a range for the case fatality ratio of 0.0004% to 0.06%, and Baker et al7 estimated a symptomatic case fatality ratio of 0.005%, somewhat lower than ours. Hadler et al36 estimated a range for the symptomatic case fatality ratio in New York between 0.0054% and 0.0086%, whereas Donaldson et al14 reported a symptomatic case fatality ratio for England of 0.026% (range 0.011% to 0.066%), and Pebody et al15 reported a ratio of 0.04% (range 0.02% to 0.10%). Wu et al37 estimated a case fatality ratio of 0.0109% (0.0041% to 0.0377%).
Conclusions and policy implications
While we have been careful to outline underlying assumptions and possible uncertainties in our analysis, we have nevertheless estimated severity in a robust statistical framework, systematically accounting for possible biases. Identification of biases is possible only in an evidence synthesis framework, through the “triangulation” of multiple data sources: each source on its own provides a (potentially biased) view of only one aspect of the severity of an epidemic. The credible intervals reported fully reflect the uncertainty in the observed data, the estimation process, and some (though not all) model assumptions.
Finally, our study suggests a mild pandemic, characterised by case and infection severity ratios increasing between the two influenza waves, while the process of synthesis and reconciliation of the data available from different sources has highlighted the importance of transparent design of routine influenza surveillance to derive robust estimates of key measures of severity. Multiple information streams—and an established framework to interpret these data as quickly as possible—are critically important, particularly in light of the apparent resurgence10 of the pandemic A/H1N1 strain in the 2010-11 influenza season.
What is already known on this topic
Initial estimates of the severity of the 2009 A/H1N1 influenza pandemic were highly uncertain, with later estimates suggesting a relatively mild pandemic but with an age distribution of infections different from that for seasonal influenza
Previous studies of the severity of the 2009 A/H1N1 influenza pandemic in England have not fully accounted for biases, have not made comprehensive use of all available relevant evidence to quantify the uncertainty in these severity measures, and have not provided an assessment of changes in severity nor of the adequacy of routine surveillance systems
What this paper adds
This study is the first to provide a complete overview of the severity of the first two waves of the 2009 pandemic A/H1N1 influenza outbreak in England, synthesising all available relevant influenza surveillance data
The analysis fully accounts for biases and under-ascertainment of confirmed cases in each surveillance system and accounts robustly for uncertainty in the data
Notes
Cite this as: BMJ 2011;343:d5408
Footnotes
We thank all colleagues throughout the Health Protection Agency who contributed to data collection during the pandemic, Professor E Miller and her team for providing the sero-incidence data, and the medical staff of the chief medical officer who collected data on fatal cases.
Contributors: AMP performed the main statistical analysis, developed the model, and wrote the paper. RGP, BJP, and AC collated the data, contributed to model development, and helped revise the paper. BDMT and PJB contributed to statistical analysis, to model development, and to critical revisions of the paper. ML and DDA conceived and planned the project and contributed to model development and to critical revisions of the paper. RGP, ML, and DDA are the guarantors for the study. All authors gave final approval of the published paper and had full access to the data used in the study.
Funding: AMP, PJB, and DDA were funded by the UK Medical Research Council (grant MC_US_A030_0016.01). BDMT was also funded by the council (grant MC_US_A030_0022.01). DDA, RGP, BJP, and AC were funded by the UK Health Protection Agency. BJP was funded as a research fellow by the Hunter Medical Research Institute, Australia. ML was supported by the Models of Infectious Disease Agent Study Consortium from the National Institute of General Medical Sciences (award No U54GM088558). The content is the responsibility of the authors and does not necessarily represent the views of the National Institute of General Medical Sciences or the National Institutes of Health. The funders had no role in study design; the collection, analysis, and interpretation of data; in the writing of the report; or in the decision to publish the article.
Competing interests: All authors have completed the Unified Competing Interest form at www.icmje.org/coi_disclosure.pdf (available on request from the corresponding author) and declare: ML received consulting income in 2007-8 from the Avian/Pandemic Flu Registry, a project of Outcome Sciences (Cambridge, MA), which was sponsored by Roche, and has received consulting income or honorariums from Novartis and Pfizer. All other authors declare no support for the submitted work from anyone other than their employer; no financial relationships with any organisations that might have an interest in the submitted work in the previous three years; no other relationships or activities that could appear to have influenced the submitted work.
Ethical approval: None required.
Data sharing: No additional data available.
This is an open-access article distributed under the terms of the Creative Commons Attribution Non-commercial License, which permits use, distribution, and reproduction in any medium, provided the original work is properly cited, the use is non commercial and is otherwise in compliance with the license. See: http://creativecommons.org/licenses/by-nc/2.0/ and http://creativecommons.org/licenses/by-nc/2.0/legalcode.