Reducing disease burden in an influenza pandemic by targeted delivery of neuraminidase inhibitors: mathematical models in the Australian context
BMC Infectious Diseases volume 16, Article number: 552 (2016)
Many nations maintain stockpiles of neuraminidase inhibitor (NAI) antiviral agents for use in influenza pandemics to reduce transmission and mitigate the course of clinical infection. Pandemic preparedness plans include the use of these stockpiles to deliver proportionate responses, informed by emerging evidence of clinical impact. Recent uncertainty about the effectiveness of NAIs has prompted these nations to reconsider the role of NAIs in pandemic response, with implications for pandemic planning and for NAI stockpile size.
We combined a dynamic model of influenza epidemiology with a model of the clinical care pathways in the Australian health care system to identify effective NAI strategies for reducing morbidity and mortality in pandemic events, and the stockpile requirements for these strategies. The models were informed by a 2015 assessment of NAI effectiveness against susceptibility, pathogenicity, and transmission of influenza.
Liberal distribution of NAIs for early treatment in outpatient settings yielded the greatest benefits in all of the considered scenarios. Restriction of community-based treatment to risk groups was effective in those groups, but failed to prevent the large proportion of cases arising from lower risk individuals who comprise the majority of the population.
These targeted strategies are only effective if they can be deployed within the constraints of existing health care infrastructure. This finding highlights the critical importance of identifying optimal models of care delivery for effective emergency health care response.
Many developed nations maintain stockpiles of neuraminidase inhibitor (NAI) antiviral agents for use in the event of an influenza pandemic, and have developed management plans for using these stockpiles to deliver a proportionate response, informed by emerging evidence of likely clinical impact. In light of recent conflicting messages concerning the effectiveness of NAIs for reducing transmission and mitigating the course of clinical infection of pandemic influenza [1–4], many of these nations are now re-evaluating the best use of NAIs in pandemic response, with implications for future stockpile size.
The effectiveness of antiviral interventions on the transmission of pandemic influenza and the resulting burden on health care settings has previously been studied in the Australian context [5–8]. These studies have shown that NAIs are likely to be effective in constraining transmission of a pandemic virus only in a relatively small proportion of low transmissibility, high severity scenarios (where infection is highly visible to the health care system and the basic reproduction number is not much greater than 1). Accordingly, the focus has shifted from using NAIs for containment (e.g., limiting transmission until vaccines become available) to using them to mitigate complications and population impact. In these scenarios, where pharmaceutical interventions are likely unable to reduce transmission, evidence supports some degree of clinical effectiveness against complications and death  and indicates that early administration of NAIs to reduce hospitalisations, severe outcomes and death is an appropriate strategy . The benefits of such treatment appear greatest in individuals with underlying risk conditions .
Here, we use an updated assessment of parameter estimates of NAI effectiveness against susceptibility, pathogenicity, and transmission of influenza, to inform a dynamic model of influenza epidemiology in combination with a model of clinical-care pathways through the health care system. NAIs are distributed within the consultation and admission constraints of the Australian health care system, as estimated from public health care reports. This model is used to identify effective NAI strategies for reducing morbidity and mortality, and to determine the stockpile requirements (within plausible bounds) to achieve these goals. Our focus is on responding to the first pandemic wave, prior to the availability of a definitive vaccine intervention. The outcomes of this analysis are then evaluated in the global context of recommendations for updated pandemic preparedness plans in other high-income countries.
In accordance with the Australian Health Management Plan for Pandemic Influenza (AHMPPI) , we assume that all identified cases are provided with treatment and post-exposure prophylaxis (PEP) is provided to all identifiable contacts during the first four weeks of the pandemic (the Initial Action Stage). During the subsequent Targeted Action Stage, treatment and prophylaxis recommendations are rationalised and targeted for maximum effect for the remainder of the pandemic, based on impact assessment and the observed epidemiology.
For the purposes of this study, the Australian population was stratified into five distinct risk groups, whose sizes were informed by June 2014 demographic data from the Australian Bureau of Statistics : Young children aged 0–4 years comprise 6.5 % of the population, for whom NAI treatment is assumed to confer no beneficial effects [10, 13, 14]. Elderly aged 66+ years comprise 13.7 % of the population, for whom NAI treatment is assumed to confer the same benefits as per the general population. High Risk aged 5–65 years comprise 10 % of this age group (8 % of the total population) and have greater risks of requiring hospitalisation, being admitted to ICUs, and of death due to infection; early NAI treatment is assumed to confer a greater benefit for this group than for the general population. Health Care workers comprise 325,000 of the 5–65 age group (informed by expert advice and 2011 census data ); they have no significant risk factors  but are separated from the general population in order to estimate the impact of pandemic scenarios on the health care work force. General population comprise the remainder of the 5–65 age group and have no significant risk factors.
Each group can be targeted independently for treatment and/or post-exposure prophylaxis; the effects of these interventions on subsequent transmission and on case severity differ by risk group. The model framework used to investigate the effects of targeted NAI strategies combined a mechanistic compartment-model of infection (stratified by risk group) with a finite-capacity compartment model of clinical pathways, subject to likely health care capacities and mean lengths of stay for inpatient settings. We now describe these two models in turn.
The infection model has been described and analysed in previous studies [6, 8]. It is based on a classic susceptible-exposed-infectious-recovered (SEIR) paradigm. All individuals are assumed fully susceptible (S) at the outset of the epidemic, and vulnerable to acquiring infection (E) upon contact with an infectious (I) case. Once recovered (R), individuals are assumed to be fully resistant to reinfection. All simulations commence with 100 infections in the population, distributed between the E and I classes in proportion to the mean duration of the latent and infectious periods.
The model incorporates a dynamic “contact” label, applied to a fixed number of individuals drawn from the whole population each time a new infectious case appears. We define these contacts, based on the findings of sociological studies, as those people who have been sufficiently close to an infected individual to conceivably contract infection. Only contacts of an infectious case may proceed to the exposed and infectious classes, however the majority of contacts escape unscathed, returning to their original state within 72 hrs of exposure. The capture of “contacts” in the model framework allows simulation of targeted post-exposure antiviral prophylaxis (PEP) .
In modelling the delivery of antiviral agents to the population, the model accounts for (1) a drop in efficacy due to delays in distribution (e.g., due to the requirement for laboratory confirmation during the Initial Action Stage and (2) constraints on the maximum rate of delivery of antivirals to the population (related to health sector capacity). Both of these factors have previously been shown to dramatically influence the expected outcome from an antiviral intervention and to modify stockpile usage . We assumed that the maximum antiviral delivery ranged from a low estimate of 103 packets per day to an aspirational target of 105 packets per day and drew samples from a log-uniform distribution.
Clinical pathways model
Some proportion of infections will require hospitalisation (“severe cases”) while some proportion of infections will not require hospitalisation and may present to outpatient settings (“mild cases”). The proportion of mild cases that present to hospital EDs rather than to GP clinics in Australia was estimated to be 20 %, based on expert consultation. It is further assumed that a fraction of the cases that will ultimately become severe will present early in their clinical course to an outpatient setting and, should they receive early treatment, their risk of subsequent hospitalisation is reduced . Hospitalised cases have a risk of ICU admission (which varies by risk group and the provision of treatment) and the ICU admissions have a risk of death (which also varies by risk group and the provision of treatment) as shown in Fig. 1.
National consultation and admission capacities for each health care setting were informed by public reports of Australian health care infrastructure, under the assumption that in a worst-case scenario up to 50 % of total capacity in each health care setting could possibly be devoted to influenza patients (Table 1).
Patients are admitted to general wards with a mean length of stay of 5 days, and are admitted to ICUs with a mean length of stay of 10 days. Therefore, it is the prevalence of cases requiring hospitalisation that determines the available ward and ICU bed capacities for new admissions. Admissions are preferentially allocated by strata, with priority given to health care workers and high-risk adults.
In the event that there is insufficient capacity to admit a newly-presenting case, the following hierarchy of case priorities and overflows are applied:
ICU admissions are preferentially allocated in the following order: (1) health care workers; (2) high-risk adults; (3) children; (4) elderly; and (5) general adult population.
Any cases that cannot be admitted to an ICU are considered for admission to a general ward, subject to the same order of preferential allocation.
In the situation that there is insufficient capacity to admit all cases that require hospitalisation, these cases are assumed to instead present to hospital EDs.
A fixed proportion of the mild cases present to hospital EDs, subject to the same order of preference as for ICU and ward admissions.
Presentations that cannot receive consultation at an ED are assumed to present at GP clinics.
All remaining presentations present at GP clinics, subject to the same order of preference.
Presentations that cannot receive a GP consultation are unable to receive antiviral treatment, on the grounds that there was no capacity to consult with these patients.
The Australian Health Management Plan for Pandemic Influenza (AHMPPI) defines pandemic impact levels based on the clinical severity of the disease and on the transmissibility of the virus between humans, and characterises both qualities using a “Low”, “Moderate”, “High” scale . Consistent with previous modelling studies and as used in the AHMPPI, we used the pandemic scenarios defined in Table 2. The classification of past influenza pandemics according to these scenario definitions are shown in Fig. 2.
For each scenario, model uncertainties (e.g., epidemic time-course, effectiveness of NAIs) were accounted for by using Latin hypercube sampling (LHS) to randomly select model parameter combinations for 10,000 simulations. We report outcome measures in terms of their median, 5 th and 95 th percentiles, unless stated otherwise.
For the first four weeks of the epidemic—corresponding to the Initial Action Stage as defined in the AHMPPI—all identified cases are provided with treatment and post-exposure prophylaxis (PEP) is provided to all identifiable contacts. This intensive response has previously been shown necessary to have any chance of significantly reducing infection in those scenarios where such an outcome is possible. This outcome is only achievable in the particular scenario where the disease is both visible (i.e., high severity) and has limited transmissibility (i.e., R0<1.25 ) and an intensive combined treatment and prophylaxis recommendation is initiated early in the outbreak. Data collected during this stage will allow the outbreak to be classified according to the pandemic scenarios defined above.
The Initial Action Stage is followed by a Targeted Action Stage (for the remaining duration of the epidemic) where treatment and prophylaxis recommendations are revised and targeted for maximum effect as part of a proportionate response for the relevant pandemic scenario. Four alternative targeted strategies for this stage are considered here: Rx All/PEP Eld, HR Treatment of all identified cases regardless of risk stratum or setting of care, and provision of prophylaxis for individuals in the “Elderly” and “High-Risk” strata. Rx All Treatment of all identified cases regardless of risk stratum or setting of care, and no recommendation for prophylaxis. Rx At-Risk, Hosp Treatment of all identified cases in the “Children”, “Elderly”, “High-Risk” and “HCW” strata, and of all cases in hospital and ICU settings. Rx Hosp Treatment of all cases in hospital and ICU settings.
Effects of antiviral treatment and prophylaxis
The effects of NAI treatment on clinical outcomes for severe cases are listed in Table 3. We assumed that, with the provision of NAIs, high-risk adults were twice as likely to require ICU admission as other adults (0.25 vs 0.125) and were three times more likely to die due to severe infection than other adults . Using the risk ratios for total influenza-related complications reported by Falagas et al. (0.74 for otherwise healthy patients and 0.37 for high-risk patients ), we calculated the counter-factual risks of ICU admission and death in the high-risk and general adult populations (refer to the “No Rx” rows of Table 3). On the understanding that NAI treatment confers negligible benefits to children [10, 13, 14], we did not reduce either risk for the children stratum.
A recent population cohort data linkage study of the 2009 H1N1 pandemic in British Columbia reported that the hazard ratio of all-cause hospitalisation, given early NAI treatment, was approximately 0.84 in the general population and 0.52 among those with co-morbidity . We assumed that the relative risk of hospitalisation given early NAI treatment was identical to these hazard ratios: 0.52 for the high-risk stratum, 0.84 for all other adult strata, and 1.0 for children (assuming that NAI treatment confers no benefits to children). In the absence of early treatment, all severe cases require hospitalisation.
These assumptions are also consistent with a recent world-wide meta-analysis of patients hospitalised with A(H1N1)pdm09, which reported that early NAI treatment, in comparison to late NAI treatment, was associated with significant reductions in mortality and likelihood of requiring ventilatory support .
For each of the seven pandemic scenarios we performed 10,000 simulations with different combinations of parameter values, meaning that no single epidemic can characterise any of these scenarios. This design allows us to account for uncertainty in the precise nature of the epidemic itself, in the early (imprecise) estimates of transmissibility and severity obtained during initial action, in the effectiveness of NAIs to reduce susceptibility, pathogenicity, and transmission, in logistic capacities that limit NAI distribution, and in population compliance. An illustration of the variety of epidemics that were considered in this study are shown in Fig. 3, which depicts the median, 5 th and 95 th percentile epidemic curves, as measured by Clinical Attack Rate (CAR). The epidemic duration for each scenario is reported in Table 4.
Figures 4, 5, 6, 7, and 8 present an overview of the key health care outcomes for each of the pandemic scenarios and for each targeted antiviral strategy for the Targeted Action Stage. Here, we identify the key findings for each of the following pandemic scenarios: #4 (moderate severity, medium transmissibility), #6 (high severity, low transmissibility) and #7 (high severity, high transmissibility).
Moderate severity, medium transmissibility
Medium transmission scenarios are not controllable by antiviral interventions, regardless of the proportion of infections that are identifiable. In these scenarios antiviral interventions have only a marginal effect on the prevalence of infection, but can yield reductions in case severity, hospital admissions, and deaths.
Across the range of antiviral strategies for the Targeted Action Stage, the median peak ICU occupancy ranges from 119–164 (out of 1000 beds) and the median peak ward occupancy ranges from 317–420 (out of 27,600 beds), indicating that there is more than sufficient hospital capacity to care for all cases that might require hospitalisation. This holds true even for the 95th percentiles (540–723 ICU beds and 1412–1824 ward beds).
While some of the more conservative antiviral strategies cause ED capacity to be overwhelmed (in a small fraction of the simulations) by 3–7 % for a period of up to one week, this is a direct result of assuming that 1 in 5 mild presentations occurs at an ED rather than at a GP clinic. We note here that this assumption was based on expert consultation and it is unclear how to validate it against available public health data. The peak burden on GPs, on the other hand, only consumes up to 22 % of the available consultation capacity. The conclusion to draw from these observations is that there is more than sufficient capacity for all mild cases to receive consultation, but that the ED consultation capacity can be temporarily overwhelmed if a moderate proportion of these cases elect to present to EDs rather than to GP clinics.
The more liberal antiviral strategies produce median reductions of 400 fewer deaths, 370 fewer ICU admissions, and 700 fewer hospitalisations, given median stockpile usage of 650,000 to 1.3 million doses (95th percentile usage is less than 3 million doses). More conservative strategies produce median reductions of 330–380 fewer deaths, 235–330 fewer ICU admissions, and 15–427 few hospitalisations, given median stockpile usage of 60,000 to 300,000 doses (95th percentile usage is less than 600,000 doses).
High severity, low transmissibility
Low transmission scenarios are more readily controllable if a sufficient proportion of infections are identifiable and can be provided with antivirals for treatment or prophylaxis . This visibility substantially reduces the clinical attack rate (CAR) and can, in ideal circumstances, lead to successful mitigation of the epidemic.
For this scenario, the absence of antiviral interventions results in median CARs of 6.8 % in the High-Risk population and 4.9 % for the rest of the population. In contrast, the application of a 4 week Initial Action Stage followed by a Targeted Action Stage where antivirals are only provided to hospitalised cases results in median CARs of 5.5 % in the High-Risk population and 4.1 % for the rest of the population. As the antiviral strategy for the Targeted Action Stage becomes more liberal, the median CARs are reduced even further, to minimums of 0.9 % in the High-Risk population and 0.7 % for the rest of the population. The most liberal strategy (“Rx All/PEP Eld, HR”) substantially improves the relative risk of presentation (median relative risk: 0.376), hospital admission (median relative risk: 0.320), ICU admission (median relative risk: 0.223) and death (median relative risk: 0.146); see the summary tables in Additional file 1 for further details. This control is achieved with a median stockpile usage of fewer than one million doses for the most liberal antiviral strategies, where the majority of the doses are used for post-exposure prophylaxis.
Overall hospital capacity is never exceeded. ICU capacity is exceeded in a small number of simulations, but the mean and median number of days for which this occurs is zero. With even the most conservative antiviral strategy during the Targeted Action Stage, the mean and median peak ICU occupancy is 577 beds (out of 1000 beds), and the mean and median peak ward occupancy is 1394 beds (out of 27,600 beds).
The clearest indicators of the effect that antiviral interventions have on the epidemic burden are the reductions in deaths, ICU admissions and hospitalisations that are achieved, in comparison to the same scenario with no antiviral interventions. When the most conservative strategy is used during the Targeted Action Stage, and antivirals are only provided for treatment to hospitalised cases, the median reductions are 4,780 fewer deaths, 4,095 fewer ICU admissions and 5,656 fewer hospital admissions. These outcomes improve as more liberal antiviral strategies are used during the Targeted Action Stage, to median reductions of 6,437 fewer deaths, 7,793 fewer ICU admissions and 24,566 fewer hospital admissions.
High severity, high transmissibility
High transmission scenarios are not controllable by antiviral interventions, regardless of the proportion of infections that are identifiable. In these scenarios antiviral interventions cannot affect the prevalence of infection, but can yield reductions in case severity, hospital admissions, and deaths. Compared with the medium transmission scenarios, antiviral interventions will achieve greater absolute reductions in burden, but the higher transmission means that these are smaller fractional reductions in burden.
In this worst-case scenario, ICU bed capacity is exceeded for 6 weeks and ED consultation capacity is exceeded for 3–4 weeks, regardless of the antiviral strategy for the Targeted Action Stage (i.e., essentially for the duration of the pandemic). The more liberal antiviral strategies produce median reductions of 13,000 fewer deaths, 12,400 fewer ICU admissions, and 25,000 fewer hospitalisations, given median stockpile usage of 1.3 to 2.2 million doses (95th percentile usage is 5 million doses). More conservative strategies produce median reductions of 12,300–13,700 fewer deaths, 8,800–12,000 fewer ICU admissions, and 2,600–18,200 few hospitalisations, given median stockpile usage of 600,000 to 1 million doses (95th percentile usage is 2.2 million doses).
All choices of antiviral strategy for the Targeted Action Stage reduce deaths by approximately 50 % and also reduce ICU and hospital admissions, but while the proportional reduction in clinical burden is similar to the medium transmission and moderate severity scenario (above), the absolute burden remains substantially higher. In this case, the most important conclusions to draw are that antiviral strategies can substantially reduce the number of deaths in even the most severe pandemic scenarios where hospital bed capacities are substantially overwhelmed for many weeks, assuming that antivirals can continue to be delivered in a timely and effective manner (as they are through general practice in the model).
Stockpile consumption for the intensive response delivered in the Initial Action Stage is shown in Table 5 for each pandemic scenario. The salient detail is that delivering this response consumes a small number of doses, relative to the size of a national stockpile, even for the 95 th percentile across all scenarios.
Overall stockpile consumption for each scenario and each targeted antiviral strategy is shown in Fig. 6. Recall that the (imposed) maximum distribution rate varied from 103 to 105 packets per day. In the worst-case scenario (high severity, high transmissibility) and with large maximum distribution rates (e.g., at least 9×104 packets per day), the maximum rate was reached in two-thirds of the simulations with liberal treatment of cases in the community (the “Rx All/PEP Eld, HR” and “Rx All” strategies).
It is clearly evident that the most liberal strategy (“Rx All/PEP Eld, HR”) consumes a substantially greater number of treatment courses than all other strategies, due to the provision of post-exposure prophylaxis to all elderly and high-risk contacts, but generally confers no significant benefit over the more conservative “Rx All” strategy. The targeted treatment strategies (“Rx At-Risk, Hosp” and “Rx Hosp”) consume even fewer treatment courses than the “Rx All” strategy but, as identified above, at the cost of preventing substantially fewer hospital admissions.
NAI treatment strategies that allow liberal distribution of antivirals for early treatment in outpatient settings yielded the greatest benefits, evidenced by reductions in hospitalisations, critical care requirements and deaths, for all of the considered pandemic scenarios. Restriction of community-based treatment to risk groups is effective in those groups, but fails to prevent the large proportion of cases arising from lower risk individuals who comprise the majority of the population. In even the most severe scenarios, median stockpile consumption for treatment was 1.3 million doses, sufficient to cover 6.5 % of the population (95 th percentile: 3.3 million doses, sufficient to cover 16.5 % of the population).
In high severity scenarios, we predict that capacity constraints within ICUs and hospitals will be exceeded, placing greater pressure on community-based health services. NAI treatment strategies can substantially reduce the number of deaths in even the most severe pandemic scenarios where hospital bed capacities are substantially overwhelmed for many weeks, assuming that antivirals can continue to be delivered in a timely and effective manner (as they are through general practice in the model). While we have assumed that general practitioners will be able to deliver effective health care to the resulting “overflow” population, more consideration of models of care delivery in such situations is needed to ensure access to needed services.
Study strengths and weaknesses
Strengths include our use of LHS sampling and broad parameter ranges to consider a wide range of epidemics, logistic capacities, and compliance, for each of the pandemic scenarios; the spread of the results for each scenario (reported as 5th and 95th percentiles) is therefore indicative of the impact and outcomes, without being strongly tied to specific values of model parameters.
The major weaknesses are: (a) the assumption that up to 50 % of actual health care capacities can be devoted to treating influenza patients (approaching or exceeding these capacities therefore represents an “apocalpytic” pandemic, with broader ramifications for society as a whole); (b) that the health care workforce does not suffer from fatigue, absenteeism, or depletion due to illness (any of which would compromise intervention delivery, and consultation and bed capacities); (c) the role of emergency departments as “gateways” to the hospital system is not explicitly configured, including the potential for overwhelmed EDs to constrain capacity for hospital admissions; and (d) multiple presentations for a single case is only considered for severe cases who receive (ineffective) early treatment and subsequently require hospital admission, thus excluding the possibility of repeat presentations in general despite evidence that this may be a common occurrence .
Critically, the effectiveness of NAIs (see Table 3) is based on available evidence, which is mostly observational since there are no clinical trial data in hospitalised patients and in high-risk groups . The recent review by the Academy of Medical Sciences and the Wellcome Trust states that “there is a lack of evidence to guide decisions on NAI treatment for high-risk groups and children,” but also notes that “the steering group does not support the assumption that observational data are invariably of less use than data from RCTs” . The interpretation of the available evidence, however, remains a subject of some controversy.
We have not accounted for social distancing measures such as school closures in this study. Despite the lack of substantial quantitative evidence to inform mathematical models and optimal implementation , and concerns about societal costs [23, 24], such methods have been widely used in past epidemics (e.g., in the 2009 pandemic in Australia  and other countries , and in Australia during the 1918–19 pandemic ). Evidence from 2009 in the state of Victoria, Australia, is that compliance with behavioural and pharmaceutical recommendations was high, but were unlikely to have substantially altered the course of the epidemic . Due to the timing and circumstances, this likely reflects a “best case” estimate of public compliance during a moderate to severe influenza pandemic .
These limitations highlight the critical importance of considering tailored models of care delivery and how appropriate communication may influence health care seeking behaviour in target populations.
The timing of the pandemic waves with respect to season is another factor that can impact the basic reproduction number, causing substantial differences in the pandemic impact experienced in different regions of the world and confounding the applicability of overseas pathogen assessments for the local context. This reinforces the importance of having local or region-specific real-time surveillance and R0 estimation protocols in place.
We have also restricted this study to considering responses to a single pandemic wave. Multiple waves were observed in 2009 in the UK [28, 29] and Australia [30, 31], in 1918–19 in the UK [32–34] and Australia , and are typical of all pandemics in the 20th century . Accounting for secondary and tertiary waves is an integral part of pandemic preparedness, which we have not addressed here; implementation of a strain-specific vaccine is identified as the definitive measure to reduce morbidity and mortality in subsequent pandemic waves in an Australian pandemic response . In subsequent waves there may be greater knowledge about the pathogen, but external factors such as seasonality, vaccine availability , social behaviour, and changing host immunity [33, 34] may act as confounders and greatly affect the impact of these subsequent waves. All of these issues are clearly of real importance, but are beyond the scope of this study.
Meaning and implications
The results of this study suggest that the optimal use of an antiviral stockpile in the event of an influenza pandemic is to provide treatment to as many cases as possible, in both inpatient and outpatient settings. In the unlikely scenario that the pandemic strain exhibits both high clinical severity and low transmissibility, exhaustive contact-tracing and provision of post-exposure prophylaxis may be able to sufficiently reduce transmission as to mitigate the pandemic. In all other scenarios, this outcome is not achievable and the provision of post-exposure prophylaxis both confers no benefits to the population and substantially increases stockpile consumption. This finding highlights the critical importance of intensive early data gathering to inform impact assessment (e.g., based on results of “first few hundred” studies ), enabling reorientation of public health efforts for proportionate and effective response.
The key issue identified in this study is the importance of understanding how different models of care may enable the delivery effective interventions in the event of a pandemic, without overwhelming existing day-to-day and surge capacities. We have shown that a modest NAI stockpile permits liberal (early) treatment in both outpatient and inpatient settings, which can substantially reduce hospitalisations, critical care requirements and deaths in the event of an influenza pandemic. However, these outcomes are only achievable if liberal treatment strategies can be effectively deployed within the constraints of the existing health care infrastructure.
Australian Health Management Plan for Pandemic Influenza
Clinical attack rate
Health care worker
High Risk population stratum, aged 5-65 years and possessing underlying risk conditions
Intensive care unit
Latin hypercube sampling
Randomised controlled trial
- Rx All/PEP Eld, HR:
Targeted NAI strategy: treatment of all identified cases regardless of risk stratum or setting of care, and provision of prophylaxis for individuals in the “Elderly” and “High-Risk” strata
- Rx All:
Targeted NAI strategy: treatment of all identified cases regardless of risk stratum or setting of care, and no recommendation for prophylaxis
- Rx At-Risk, Hosp:
Targeted NAI strategy: treatment of all identified cases in the “Children”, “Elderly”, “High-Risk” and “HCW” strata, and of all cases in hospital and ICU settings
- Rx Hosp:
Targeted NAI strategy: treatment of all cases in hospital and ICU settings
Susceptible-exposed-infectious-recovered compartment model
Wang K, Shun-Shin M, Gill P, Perera R, Harnden A. Neuraminidase inhibitors for preventing and treating influenza in children (published trials only). Cochrane Database Syst Rev. 2012; 4:002744. doi:http://dx.doi.org/10.1002/14651858.cd002744.pub4.
Ebell MH, Call M, Shinholser J. Effectiveness of oseltamivir in adults: a meta-analysis of published and unpublished clinical trials. Fam Pract. 2012; 30(2):125–33. doi:http://dx.doi.org/10.1093/fampra/cms059.
Jefferson T, Jones MA, Doshi P, Del Mar CB, Hama R, Thompson MJ, Spencer EA, Onakpoya IJ, Mahtani KR, Nunan D, et al.Neuraminidase inhibitors for preventing and treating influenza in adults and children. Cochrane Database Syst Rev. 2014; 4:008965. doi:http://dx.doi.org/10.1002/14651858.cd008965.pub4.
Dobson J, Whitley RJ, Pocock S, Monto AS. Oseltamivir treatment for influenza in adults: a meta-analysis of randomised controlled trials. Lancet. 2015; 385(9979):1729–1737. doi:http://dx.doi.org/10.1016/s0140-6736(14)62449-1.
McVernon J, McCaw JM, Nolan TM. Modelling strategic use of the national antiviral stockpile during the CONTAIN and SUSTAIN phases of an Australian pandemic influenza response. Aust N Z J Public Health. 2010; 34(2):113–9. doi:http://dx.doi.org/10.1111/j.1753-6405.2010.00493.x.
Moss R, McCaw JM, McVernon J. Diagnosis and antiviral intervention strategies for mitigating an influenza epidemic. PLoS ONE. 2011; 6(2):14505. doi:http://dx.doi.org/10.1371/journal.pone.0014505.
McCaw JM, Moss R, McVernon J. A decision support tool for evaluating the impact of a diagnostic-capacity and antiviral-delivery constrained intervention strategy on an influenza pandemic. Influenza Other Respir Viruses. 2011; 5:212–5. doi:http://dx.doi.org/10.1111/j.1750-2659.2011.00209.x.
Dafilis MP, Moss R, McVernon J, McCaw J. Drivers and consequences of influenza antiviral resistant-strain emergence in a capacity-constrained pandemic response. Epidemics. 2012; 4(4):219–26. doi:http://dx.doi.org/10.1016/j.epidem.2012.12.002.
The Academy of Medical Sciences. Wellcome Trust: Use of Neuraminidase Inhibitors in Influenza. London: The Academy of Medical Sciences and Wellcome Trust; 2015.
Marra F, Chong M, Henry B, Patrick DM, Kendall P. Effectiveness of neuraminidase inhibitors in preventing hospitalization during the H1N1 influenza pandemic in British Columbia, Canada. J Antimicrob Chemother. 2013; 69(5):1397–1406. doi:http://dx.doi.org/10.1093/jac/dkt496.
Department of Health: Australian health management plan for pandemic influenza. Technical report, Australian Government 2014. http://www.health.gov.au/internet/main/publishing.nsf/Content/ohp-ahmppi.htm. Accessed 29 Sept 2016.
Australian Bureau of Statistics: Australian Demographic Statistics, Jun 2014. Technical Report 3101.0, Australian Bureau of Statistics 2014. http://www.abs.gov.au/AUSSTATS/abs@.nsf/DetailsPage/3101.0Jun%202014. Accessed 29 Sept 2016.
Muthuri SG, Venkatesan S, Myles PR, Leonardi-Bee J, Al Khuwaitir TSA, Al Mamun A, Anovadiya AP, Azziz-Baumgartner E, Báez C, Bassetti M, et al. Effectiveness of neuraminidase inhibitors in reducing mortality in patients admitted to hospital with influenza A H1N1pdm09 virus infection: a meta-analysis of individual participant data. Lancet Respir Med. 2014; 2(5):395–404. doi:http://dx.doi.org/10.1016/s2213-2600(14)70041-4.
Muthuri SG, Venkatesan S, Myles PR, Leonardi-Bee J, Lim WS, Mamun AA, Anovadiya AP, Araújo WN, Azziz-Baumgartner E, Báez C, et al.Impact of neuraminidase inhibitors on influenza A(H1N1)pdm09-related pneumonia: an IPD meta-analysis. Influenza Other Respir Viruses. 2015. doi:http://dx.doi.org/10.1111/irv.12363. Accepted, unedited article published online and citable.
Australian Bureau of Statistics: Australian Social Trends: Doctors and Nurses. Technical report, Australian Bureau of Statistics 2013. http://www.abs.gov.au/AUSSTATS/abs@.nsf/Lookup/4102.0Main+Features20April+2013. Accessed 29 Sept 2016.
Higa F, Tateyama M, Tomishima M, Mukatake S, Yamashiro T, Owan T, Haranaga S, Fujita J. Role of neuraminidase inhibitor chemoprophylaxis in controlling nosocomial influenza: an observational study. Influenza Other Respir Viruses. 2011; 6(4):299–303. doi:http://dx.doi.org/10.1111/j.1750-2659.2011.00311.x.
McCaw JM, McVernon J. Prophylaxis or treatment? optimal use of an antiviral stockpile during an influenza pandemic. Math Biosci. 2007; 209(2):336–60. doi:http://dx.doi.org/10.1016/j.mbs.2007.02.003.
McCaw JM, Glass K, Mercer GN, McVernon J. Pandemic controllability: a concept to guide a proportionate and flexible operational response to future influenza pandemics. J Public Health. 2014; 36(1):5–12. doi:http://dx.doi.org/10.1093/pubmed/fdt058.
Van Kerkhove MD, Vandemaele KAH, Shinde V, Jaramillo-Gutierrez G, Koukounari A, Donnelly CA, Carlino LO, Owen R, Paterson B, Pelletier L, et al.Risk factors for severe outcomes following 2009 influenza A (H1N1) infection: A global pooled analysis. PLoS Med. 2011; 8(7):1001053. doi:http://dx.doi.org/10.1371/journal.pmed.1001053.
Falagas ME, Koletsi PK, Vouloumanou EK, Rafailidis PI, Kapaskelis AM, Rello J. Effectiveness and safety of neuraminidase inhibitors in reducing influenza complications: a meta-analysis of randomized controlled trials. J Antimicrob Chemother. 2010; 65(7):1330–1346. doi:http://dx.doi.org/10.1093/jac/dkq158.
Pandemic Influenza Preparedness Team: The National Pandemic Flu Service: an evaluation. Technical Report 15652, UK Department of Health 2011. https://www.gov.uk/government/uploads/system/uploads/attachment_data/file/215679/dh_125338.pdf. Accessed 29 Sept 2016.
Cauchemez S, Ferguson NM, Wachtel C, Tegnell A, Saour G, Duncan B, Nicoll A. Closure of schools during an influenza pandemic. Lancet Infect Dis. 2009; 9(8):473–81. doi:http://dx.doi.org/10.1016/s1473-3099(09)70176-8.
Verikios G, McCaw JM, McVernon J, Harris AH. H1N1 influenza and the Australian macroeconomy. J Asia Pac Econ. 2012; 17(1):22–51. doi:http://dx.doi.org/10.1080/13547860.2012.639999.
Brown ST, Tai JH, Bailey RR, Cooley PC, Wheaton WD, Potter MA, Voorhees RE, LeJeune M, Grefenstette JJ, Burke DS, et al.Would school closure for the 2009 H1N1 influenza epidemic have been worth the cost?: a computational simulation of Pennsylvania. BMC Public Health. 2011; 11(1):353. doi:http://dx.doi.org/10.1186/1471-2458-11-353.
McVernon J, Mason K, Petrony S, Nathan P, LaMontagne AD, Bentley R, Fielding J, Studdert DM, Kavanagh A. Recommendations for and compliance with social restrictions during implementation of school closures in the early phase of the influenza A (H1N1) 2009 outbreak in Melbourne, Australia. BMC Infect Dis. 2011; 11(1):257. doi:http://dx.doi.org/10.1186/1471-2334-11-257.
Leung GM, Nicoll A. Reflections on pandemic (H1N1) 2009 and the international response. PLoS Med. 2010; 7(10):1000346. doi:http://dx.doi.org/10.1371/journal.pmed.1000346.
Caley P, Philp DJ, McCracken K. Quantifying social distancing arising from pandemic influenza. J R Soc Interface. 2008; 5(23):631–9. doi:http://dx.doi.org/10.1098/rsif.2007.1197.
Hardelid P, Andrews N, Hoschler K, Stanford E, Baguelin M, Waight P, Zambon M, Miller E. Assessment of baseline age-specific antibody prevalence and incidence of infection to novel influenza AH1N1 2009. Health Technol Assess. 2010; 14(55):3. doi:http://dx.doi.org/10.3310/hta14550-03.
Mytton OT, Rutter PD, Donaldson LJ. Influenza A(H1N1)pdm09 in England, 2009 to 2011: a greater burden of severe illness in the year after the pandemic than in the pandemic year. Euro Surveill. 2012; 17(14):20139.
McVernon J, Laurie K, Nolan T, Owen R, Irving D, Capper H, Hyland C, Faddy H, Carolan L, Barr I, Kelso A. Seroprevalence of 2009 pandemic influenza A(H1N1) virus in Australian blood donors, October–December 2009. Euro Surveill. 2010; 15(40):19678.
McVernon J, Laurie K, Faddy H, Irving D, Nolan T, Barr I, Kelso A. Seroprevalence of antibody to influenza A(H1N1)pdm09 attributed to vaccination or infection, before and after the second year=2010 pandemic wave in Australia. Influenza Other Respir Viruses. 2013; 8(2):194–200. doi:http://dx.doi.org/10.1111/irv.12225.
Chowell G, Bettencourt LMA, Johnson N, Alonso WJ, Viboud C. The 1918–1919 influenza pandemic in England and Wales: spatial patterns in transmissibility and mortality impact. Proc R Soc B. 2008; 275(1634):501–9. doi:http://dx.doi.org/10.1098/rspb.2007.1477.
Mathews JD, McBryde ES, McVernon J, Pallaghy PK, McCaw JM. Prior immunity helps to explain wave-like behaviour of pandemic influenza in 1918–9. BMC Infect Dis. 2010; 10(1):128. doi:http://dx.doi.org/10.1186/1471-2334-10-128.
Bolton KJ, McCaw JM, McVernon J, Mathews JD. The influence of changing host immunity on 1918–19 pandemic dynamics. Epidemics. 2014; 8:18–27. doi:http://dx.doi.org/10.1016/j.epidem.2014.07.004.
Miller MA, Viboud C, Balinska M, Simonsen L. The signature features of influenza pandemics — implications for policy. N Engl J Med. 2009; 360(25):2595–598. doi:http://dx.doi.org/10.1056/nejmp0903906.
McLean E, Pebody RG, Campbell C, Chamberland M, Hawkins C, Nguyen-Van-Tam JS, Oliver I, Smith GE, Ihekweazu C, Bracebridge S, Maguire H, Harris R, Kafatos G, White PJ, Wynne-Evans E, Green J, Myers R, Underwood A, Dallman T, Wreghitt T, Zambon M, Ellis J, Phin N, Smyth B, McMenamin J, Watson JM. Pandemic (H1N1) 2009 influenza in the UK: clinical and epidemiological findings from the first few hundred (FF100) cases. Epidemiol Infect. 2010; 138(11):1531–1541. doi:http://dx.doi.org/10.1017/s0950268810001366.
ANZICS: ANZICS 2014 annual report. Technical report, ANZICS 2014. http://www.anzics.com.au/Downloads/ANZICS%20Annual%20Report%202014.pdf. Accessed 29 Sept 2016.
AIHW: Australia’s hospitals 2009-10 at a glance. Technical report, AIHW 2011. http://www.aihw.gov.au/publication-detail/?id=10737418683. Accessed 29 Sept 2016.
AIHW: Australian hospital statistics 2011-12: Emergency department care. Health Services Series HSE 126. AIHW. 2012. http://www.aihw.gov.au/publication-detail/?id=10737423042. Accessed 29 Sept 2016.
AHWAC: Health workforce planning and models of care in emergency departments. Technical Report 2006.X. AHWAC. 2006. http://docplayer.net/15513519-Health-workforce-planning-and-models-of-care-in-emergencydepartments.html. Accessed 29 Sept 2016.
AIHW: Australia’s health 2012, Australia’s Health Series AUS 156. AIHW. 2012. http://www.aihw.gov.au/publication-detail/?id=10737422172. Accessed 29 Sept 2016.
The Melbourne WHO Collaborating Centre for Reference and Research on Influenza is supported by the Australian Government Department of Health.
This work was funded by a solicited tender, “Review current evidence on the use of neuraminidase inhibitors held in the National Medical Stockpile, in a pandemic,” from the Office of Health Protection, Australian Government Department of Health.
Jodie McVernon was supported by an Australian Government National Health and Medical Research Council Career Development Award (CDF1061321). James M McCaw was supported by an Australian Research Council Future Fellowship (FT11010025).
Availability of data and materials
The model description is provided in Additional file 1 (“Model description and additional results tables”). The source code used to perform all of the simulations and generate the results presented in this manuscript is published online (“Disease transmission and clinical pathways models”, doi:http://dx.doi.org/10.4225/49/57BA7BD0CC3D2).
RM, JMM, ACC, ACH and JM contributed to the experimental design. RM adapted the simulation models and performed the simulations. RM, JMM, ACC, ACH and JM contributed to the analysis of the simulation results. RM, JMM, ACC, ACH and JM contributed to authoring the manuscript. All authors have read and approved the final version of this manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Model description and additional results tables. This document presents the model equations, the distributions from which the model parameters were sampled, and tables of simulation results for each pandemic scenario and for each targeted NAI strategy. (PDF 103 kb)
About this article
Cite this article
Moss, R., McCaw, J.M., Cheng, A.C. et al. Reducing disease burden in an influenza pandemic by targeted delivery of neuraminidase inhibitors: mathematical models in the Australian context. BMC Infect Dis 16, 552 (2016). https://doi.org/10.1186/s12879-016-1866-7