In Fig 5, the stacked bar heights normalization is incorrect.The height of each stacked bar should be 1.0 instead of 0.8.Additionally, there is an error in the legend.The shade red should be associated with "Recovered Severe Infected Case (R S )" and yellow with "Recovered Mild Infected Case (R M )".Please see the correct Fig 5 and captions here.In S1 Text, there are errors in the legend of the Figures D, L, P, and T. The shade red should correspond to "Recovered severe infected case (R S )" and yellow to "Recovered mild infected case (R M )".Please view the correct S1 Text below.

BACKGROUND: Respiratory pathogens inflict a substantial burden on public health and the economy. Although the severity of symptoms caused by these pathogens can vary from asymptomatic to fatal, the factors that determine symptom severity are not fully understood. Correlations in symptoms between infector-infectee pairs, for which evidence is accumulating, can generate large-scale clusters of severe infections that could be devastating to those most at risk, whilst also conceivably leading to chains of mild or asymptomatic infections that generate widespread immunity with minimal cost to public health. Although this effect could be harnessed to amplify the impact of interventions that reduce symptom severity, the mechanistic representation of symptom propagation within mathematical and health economic modelling of respiratory diseases is understudied. METHODS AND FINDINGS: We propose a novel framework for incorporating different levels of symptom propagation into models of infectious disease transmission via a single parameter, α. Varying α tunes the model from having no symptom propagation (α = 0, as typically assumed) to one where symptoms always propagate (α = 1). For parameters corresponding to three respiratory pathogens-seasonal influenza, pandemic influenza and SARS-CoV-2-we explored how symptom propagation impacted the relative epidemiological and health-economic performance of three interventions, conceptualised as vaccines with different actions: symptom-attenuating (labelled SA), infection-blocking (IB) and infection-blocking admitting only mild breakthrough infections (IB_MB). In the absence of interventions, with fixed underlying epidemiological parameters, stronger symptom propagation increased the proportion of cases that were severe. For SA and IB_MB, interventions were more effective at reducing prevalence (all infections and severe cases) for higher strengths of symptom propagation. For IB, symptom propagation had no impact on effectiveness, and for seasonal influenza this intervention type was more effective than SA at reducing severe infections for all strengths of symptom propagation. For pandemic influenza and SARS-CoV-2, at low intervention uptake, SA was more effective than IB for all levels of symptom propagation; for high uptake, SA only became more effective under strong symptom propagation. Health economic assessments found that, for SA-type interventions, the amount one could spend on control whilst maintaining a cost-effective intervention (termed threshold unit intervention cost) was very sensitive to the strength of symptom propagation. CONCLUSIONS: Overall, the preferred intervention type depended on the combination of the strength of symptom propagation and uptake. Given the importance of determining robust public health responses, we highlight the need to gather further data on symptom propagation, with our modelling framework acting as a template for future analysis.

The rapid development of vaccines against SARS-CoV-2 altered the course of the COVID-19 pandemic. In most countries, vaccinations were initially targeted at high-risk populations, including older individuals and healthcare workers. Now, despite substantial infection- and vaccine-induced immunity in host populations worldwide, waning immunity and the emergence of novel variants continue to cause significant waves of infection and disease. Policy makers must determine how to deploy booster vaccinations, particularly when constraints in vaccine supply, delivery and cost mean that booster vaccines cannot be administered to everyone. A key question is therefore whether older individuals should again be prioritised for vaccination, or whether alternative strategies (e.g. offering booster vaccines to the individuals who have most contacts with others and therefore drive infection) can instead offer indirect protection to older individuals. Here, we use mathematical modelling to address this question, considering SARS-CoV-2 transmission in a range of countries with different socio-economic backgrounds. We show that the population structures of different countries can have a pronounced effect on the impact of booster vaccination, even when identical booster vaccination targeting strategies are adopted. However, under the assumed transmission model, prioritising older individuals for booster vaccination consistently leads to the most favourable public health outcomes in every setting considered. This remains true for a range of assumptions about booster vaccine supply and timing, and for different assumed policy objectives of booster vaccination.

Vector control is a vital tool utilised by malaria control and elimination programmes worldwide, and as such it is important that we can accurately quantify the expected public health impact of these methods. There are very few previous models that consider vector-control-induced changes in the age-structure of the vector population and the resulting impact on transmission. We analytically derive the steady-state solution of a novel age-structured deterministic compartmental model describing the mosquito feeding cycle, with mosquito age represented discretely by parity-the number of cycles (or successful bloodmeals) completed. Our key model output comprises an explicit, analytically tractable solution that can be used to directly quantify key transmission statistics, such as the effective reproductive ratio under control, Rc, and investigate the age-structured impact of vector control. Application of this model reinforces current knowledge that adult-acting interventions, such as indoor residual spraying of insecticides (IRS) or long-lasting insecticidal nets (LLINs), can be highly effective at reducing transmission, due to the dual effects of repelling and killing mosquitoes. We also demonstrate how larval measures can be implemented in addition to adult-acting measures to reduce Rc and mitigate the impact of waning insecticidal efficacy, as well as how mid-ranges of LLIN coverage are likely to experience the largest effect of reduced net integrity on transmission. We conclude that whilst well-maintained adult-acting vector control measures are substantially more effective than larval-based interventions, incorporating larval control in existing LLIN or IRS programmes could substantially reduce transmission and help mitigate any waning effects of adult-acting measures.

[This corrects the article DOI: 10.1371/journal.pcbi.1007096.].

We discuss the invasion of the Omicron BA.1 variant into England as a paradigm for real-time model fitting and projection. Here we use a mixture of simple SIR-type models, analysis of the early data and a more complex age-structure model fit to the outbreak to understand the dynamics. In particular, we highlight that early data shows that the invading Omicron variant had a substantial growth advantage over the resident Delta variant. However, early data does not allow us to reliably infer other key epidemiological parameters-such as generation time and severity-which influence the expected peak hospital numbers. With more complete epidemic data from January 2022 are we able to capture the true scale of the epidemic in terms of both infections and hospital admissions, driven by different infection characteristics of Omicron compared to Delta and a substantial shift in estimated precautionary behaviour during December. This work highlights the challenges of real time forecasting, in a rapidly changing environment with limited information on the variant's epidemiological characteristics.

The rapid detection of outbreaks is a key step in the effective control and containment of infectious diseases. In particular, the identification of cases which might be epidemiologically linked is crucial in directing outbreak-containment efforts and shaping the intervention of public health authorities. Often this requires the detection of clusters of cases whose numbers exceed those expected by a background of sporadic cases. Quantifying exceedances rapidly is particularly challenging when only few cases are typically reported in a precise location and time. To address such important public health concerns, we present a general method which can detect spatio-temporal deviations from a Poisson point process and estimate the odds of an isolate being part of a cluster. This method can be applied to diseases where detailed geographical information is available. In addition, we propose an approach to explicitly take account of delays in microbial typing. As a case study, we considered invasive group A Streptococcus infection events as recorded and typed by Public Health England from 2015 to 2020.

The spread of infection amongst livestock depends not only on the traits of the pathogen and the livestock themselves, but also on the veterinary health behaviours of farmers and how this impacts their implementation of disease control measures. Controls that are costly may make it beneficial for individuals to rely on the protection offered by others, though that may be sub-optimal for the population. Failing to account for socio-behavioural properties may produce a substantial layer of bias in infectious disease models. We investigated the role of heterogeneity in vaccine response across a population of farmers on epidemic outbreaks amongst livestock, caused by pathogens with differential speed of spread over spatial landscapes of farms for two counties in England (Cumbria and Devon). Under different compositions of three vaccine behaviour groups (precautionary, reactionary, non-vaccination), we evaluated from population- and individual-level perspectives the optimum threshold distance to premises with notified infection that would trigger responsive vaccination by the reactionary vaccination group. We demonstrate a divergence between population and individual perspectives in the optimal scale of reactive voluntary vaccination response. In general, minimising the population-level perspective cost requires a broader reactive uptake of the intervention, whilst optimising the outcome for the average individual increased the likelihood of larger scale disease outbreaks. When the relative cost of vaccination was low and the majority of premises had undergone precautionary vaccination, then adopting a perspective that optimised the outcome for an individual gave a broader spatial extent of reactive response compared to a perspective wanting to optimise outcomes for everyone in the population. Under our assumed epidemiological context, the findings identify livestock disease intervention receptiveness and cost combinations where one would expect strong disagreement between the intervention stringency that is best from the perspective of a stakeholder responsible for supporting the livestock industry compared to a sole livestock owner. Were such discord anticipated and achieving a consensus view across perspectives desired, the findings may also inform those managing veterinary health policy the requisite reduction in intervention cost and/or the required extent of nurturing beneficial community attitudes towards interventions.

The widespread, and in many countries unprecedented, use of non-pharmaceutical interventions (NPIs) during the COVID-19 pandemic has highlighted the need for mathematical models which can estimate the impact of these measures while accounting for the highly heterogeneous risk profile of COVID-19. Models accounting either for age structure or the household structure necessary to explicitly model many NPIs are commonly used in infectious disease modelling, but models incorporating both levels of structure present substantial computational and mathematical challenges due to their high dimensionality. Here we present a modelling framework for the spread of an epidemic that includes explicit representation of age structure and household structure. Our model is formulated in terms of tractable systems of ordinary differential equations for which we provide an open-source Python implementation. Such tractability leads to significant benefits for model calibration, exhaustive evaluation of possible parameter values, and interpretability of results. We demonstrate the flexibility of our model through four policy case studies, where we quantify the likely benefits of the following measures which were either considered or implemented in the UK during the current COVID-19 pandemic: control of within- and between-household mixing through NPIs; formation of support bubbles during lockdown periods; out-of-household isolation (OOHI); and temporary relaxation of NPIs during holiday periods. Our ordinary differential equation formulation and associated analysis demonstrate that multiple dimensions of risk stratification and social structure can be incorporated into infectious disease models without sacrificing mathematical tractability. This model and its software implementation expand the range of tools available to infectious disease policy analysts.

Rapid testing strategies that replace the isolation of close contacts through the use of lateral flow device tests (LFTs) have been suggested as a way of controlling SARS-CoV-2 transmission within schools that maintain low levels of pupil absences. We developed an individual-based model of a secondary school formed of exclusive year group bubbles (five year groups, with 200 pupils per year) to assess the likely impact of strategies using LFTs in secondary schools over the course of a seven-week half-term on transmission, absences, and testing volume, compared to a policy of isolating year group bubbles upon a pupil returning a positive polymerase chain reaction (PCR) test. We also considered the sensitivity of results to levels of participation in rapid testing and underlying model assumptions. While repeated testing of year group bubbles following case detection is less effective at reducing infections than a policy of isolating year group bubbles, strategies involving twice weekly mass testing can reduce infections to lower levels than would occur under year group isolation. By combining regular testing with serial contact testing or isolation, infection levels can be reduced further still. At high levels of pupil participation in lateral flow testing, strategies replacing the isolation of year group bubbles with testing substantially reduce absences, but require a high volume of testing. Our results highlight the conflict between the goals of minimising within-school transmission, minimising absences and minimising testing burden. While rapid testing strategies can reduce school transmission and absences, they may lead to a large number of daily tests.

Gambiense human African trypanosomiasis (gHAT, sleeping sickness) is one of several neglected tropical diseases (NTDs) where there is evidence of asymptomatic human infection but there is uncertainty of the role it plays in transmission and maintenance. To explore possible consequences of asymptomatic infections, particularly in the context of elimination of transmission-a goal set to be achieved by 2030-we propose a novel dynamic transmission model to account for the asymptomatic population. This extends an established framework, basing infection progression on a number of experimental and observation gHAT studies. Asymptomatic gHAT infections include those in people with blood-dwelling trypanosomes, but no discernible symptoms, or those with parasites only detectable in skin. Given current protocols, asymptomatic infection with blood parasites may be diagnosed and treated, based on observable parasitaemia, in contrast to many other diseases for which treatment (and/or diagnosis) may be based on symptomatic infection. We construct a model in which exposed people can either progress to either asymptomatic skin-only parasite infection, which would not be diagnosed through active screening algorithms, or blood-parasite infection, which is likely to be diagnosed if tested. We add extra parameters to the baseline model including different self-cure, recovery, transmission and detection rates for skin-only or blood infections. Performing sensitivity analysis suggests all the new parameters introduced in the asymptomatic model can impact the infection dynamics substantially. Among them, the proportion of exposures resulting in initial skin or blood infection appears the most influential parameter. For some plausible parameterisations, an initial fall in infection prevalence due to interventions could subsequently stagnate even under continued screening due to the formation of a new, lower endemic equilibrium. Excluding this scenario, our results still highlight the possibility for asymptomatic infection to slow down progress towards elimination of transmission. Location-specific model fitting will be needed to determine if and where this could pose a threat.

Identification of those at greatest risk of death due to the substantial threat of COVID-19 can benefit from novel approaches to epidemiology that leverage large datasets and complex machine-learning models, provide data-driven intelligence, and guide decisions such as intensive-care unit admission (ICUA). The objective of this study is two-fold, one substantive and one methodological: substantively to evaluate the association of demographic and health records with two related, yet different, outcomes of severe COVID-19 (viz., death and ICUA); methodologically to compare interpretations based on logistic regression and on gradient-boosted decision tree (GBDT) predictions interpreted by means of the Shapley impacts of covariates. Very different association of some factors, e.g., obesity and chronic respiratory diseases, with death and ICUA may guide review of practice. Shapley explanation of GBDTs identified varying effects of some factors among patients, thus emphasising the importance of individual patient assessment. The results of this study are also relevant for the evaluation of complex automated clinical decision systems, which should optimise prediction scores whilst remaining interpretable to clinicians and mitigating potential biases.

Gambiense human African trypanosomiasis (gHAT) is a virulent disease declining in burden but still endemic in West and Central Africa. Although it is targeted for elimination of transmission by 2030, there remain numerous questions about the drivers of infection and how these vary geographically. In this study we focus on the Democratic Republic of Congo (DRC), which accounted for 84% of the global case burden in 2016, to explore changes in transmission across the country and elucidate factors which may have contributed to the persistence of disease or success of interventions in different regions. We present a Bayesian fitting methodology, applied to 168 endemic health zones (∼100,000 population size), which allows for calibration of a mechanistic gHAT model to case data (from the World Health Organization HAT Atlas) in an adaptive and automated framework. It was found that the model needed to capture improvements in passive detection to match observed trends in the data within former Bandundu and Bas Congo provinces indicating these regions have substantially reduced time to detection. Health zones in these provinces generally had longer burn-in periods during fitting due to additional model parameters. Posterior probability distributions were found for a range of fitted parameters in each health zone; these included the basic reproduction number estimates for pre-1998 (R0) which was inferred to be between 1 and 1.14, in line with previous gHAT estimates, with higher median values typically in health zones with more case reporting in the 2000s. Previously, it was not clear whether a fall in active case finding in the period contributed to the declining case numbers. The modelling here accounts for variable screening and suggests that underlying transmission has also reduced greatly-on average 96% in former Equateur, 93% in former Bas Congo and 89% in former Bandundu-Equateur and Bandundu having had the highest case burdens in 2000. This analysis also sets out a framework to enable future predictions for the country.

As part of a concerted pandemic response to protect public health, businesses can enact non-pharmaceutical controls to minimise exposure to pathogens in workplaces and premises open to the public. Amendments to working practices can lead to the amount, duration and/or proximity of interactions being changed, ultimately altering the dynamics of disease spread. These modifications could be specific to the type of business being operated. We use a data-driven approach to parameterise an individual-based network model for transmission of SARS-CoV-2 amongst the working population, stratified into work sectors. The network is comprised of layered contacts to consider the risk of spread in multiple encounter settings (workplaces, households, social and other). We analyse several interventions targeted towards working practices: mandating a fraction of the population to work from home; using temporally asynchronous work patterns; and introducing measures to create 'COVID-secure' workplaces. We also assess the general role of adherence to (or effectiveness of) isolation and test and trace measures and demonstrate the impact of all these interventions across a variety of relevant metrics. The progress of the epidemic can be significantly hindered by instructing a significant proportion of the workforce to work from home. Furthermore, if required to be present at the workplace, asynchronous work patterns can help to reduce infections when compared with scenarios where all workers work on the same days, particularly for longer working weeks. When assessing COVID-secure workplace measures, we found that smaller work teams and a greater reduction in transmission risk reduced the probability of large, prolonged outbreaks. Finally, following isolation guidance and engaging with contact tracing without other measures is an effective tool to curb transmission, but is highly sensitive to adherence levels. In the absence of sufficient adherence to non-pharmaceutical interventions, our results indicate a high likelihood of SARS-CoV-2 spreading widely throughout a worker population. Given the heterogeneity of demographic attributes across worker roles, in addition to the individual nature of controls such as contact tracing, we demonstrate the utility of a network model approach to investigate workplace-targeted intervention strategies and the role of test, trace and isolation in tackling disease spread.

Efforts to suppress transmission of SARS-CoV-2 in the UK have seen non-pharmaceutical interventions being invoked. The most severe measures to date include all restaurants, pubs and cafes being ordered to close on 20th March, followed by a "stay at home" order on the 23rd March and the closure of all non-essential retail outlets for an indefinite period. Government agencies are presently analysing how best to develop an exit strategy from these measures and to determine how the epidemic may progress once measures are lifted. Mathematical models are currently providing short and long term forecasts regarding the future course of the COVID-19 outbreak in the UK to support evidence-based policymaking. We present a deterministic, age-structured transmission model that uses real-time data on confirmed cases requiring hospital care and mortality to provide up-to-date predictions on epidemic spread in ten regions of the UK. The model captures a range of age-dependent heterogeneities, reduced transmission from asymptomatic infections and produces a good fit to the key epidemic features over time. We simulated a suite of scenarios to assess the impact of differing approaches to relaxing social distancing measures from 7th May 2020 on the estimated number of patients requiring inpatient and critical care treatment, and deaths. With regard to future epidemic outcomes, we investigated the impact of reducing compliance, ongoing shielding of elder age groups, reapplying stringent social distancing measures using region based triggers and the role of asymptomatic transmission. We find that significant relaxation of social distancing measures from 7th May onwards can lead to a rapid resurgence of COVID-19 disease and the health system being quickly overwhelmed by a sizeable, second epidemic wave. In all considered age-shielding based strategies, we projected serious demand on critical care resources during the course of the pandemic. The reintroduction and release of strict measures on a regional basis, based on ICU bed occupancy, results in a long epidemic tail, until the second half of 2021, but ensures that the health service is protected by reintroducing social distancing measures for all individuals in a region when required. Our work confirms the effectiveness of stringent non-pharmaceutical measures in March 2020 to suppress the epidemic. It also provides strong evidence to support the need for a cautious, measured approach to relaxation of lockdown measures, to protect the most vulnerable members of society and support the health service through subduing demand on hospital beds, in particular bed occupancy in intensive care units.

The COVID-19 outbreak has highlighted our vulnerability to novel infections. Faced with this threat and no effective treatment, in line with many other countries, the UK adopted enforced social distancing (lockdown) to reduce transmission-successfully reducing the reproductive number R below one. However, given the large pool of susceptible individuals that remain, complete relaxation of controls is likely to generate a substantial further outbreak. Vaccination remains the only foreseeable means of both containing the infection and returning to normal interactions and behaviour. Here, we consider the optimal targeting of vaccination within the UK, with the aim of minimising future deaths or quality adjusted life year (QALY) losses. We show that, for a range of assumptions on the action and efficacy of the vaccine, targeting older age groups first is optimal and may be sufficient to stem the epidemic if the vaccine prevents transmission as well as disease.

We explore the spatial and temporal spread of the novel SARS-CoV-2 virus under containment measures in three European countries based on fits to data of the early outbreak. Using data from Spain and Italy, we estimate an age dependent infection fatality ratio for SARS-CoV-2, as well as risks of hospitalization and intensive care admission. We use them in a model that simulates the dynamics of the virus using an age structured, spatially detailed agent based approach, that explicitly incorporates governmental interventions and changes in mobility and contact patterns occurred during the COVID-19 outbreak in each country. Our simulations reproduce several of the features of its spatio-temporal spread in the three countries studied. They show that containment measures combined with high density are responsible for the containment of cases within densely populated areas, and that spread to less densely populated areas occurred during the late stages of the first wave. The capability to reproduce observed features of the spatio-temporal dynamics of SARS-CoV-2 makes this model a potential candidate for forecasting the dynamics of SARS-CoV-2 in other settings, and we recommend its application in low and lower-middle income countries which remain understudied.

For infectious disease prevention, policy-makers are typically required to base policy decisions in light of operational and monetary restrictions, prohibiting implementation of all candidate interventions. To inform the evidence-base underpinning policy decision making, mathematical and health economic modelling can be a valuable constituent. Applied to England, this study aims to identify the optimal target age groups when extending a seasonal influenza vaccination programme of at-risk individuals to those individuals at low risk of developing complications following infection. To perform this analysis, we utilise an age- and strain-structured transmission model that includes immunity propagation mechanisms which link prior season epidemiological outcomes to immunity at the beginning of the following season. Making use of surveillance data from the past decade in conjunction with our dynamic model, we simulate transmission dynamics of seasonal influenza in England from 2012 to 2018. We infer that modified susceptibility due to natural infection in the previous influenza season is the only immunity propagation mechanism to deliver a non-negligible impact on the transmission dynamics. Further, we discerned case ascertainment to be higher for young infants compared to adults under 65 years old, and uncovered a decrease in case ascertainment as age increased from 65 to 85 years of age. Our health economic appraisal sweeps vaccination age space to determine threshold vaccine dose prices achieving cost-effectiveness under differing paired strategies. In particular, we model offering vaccination to all those low-risk individuals younger than a given age (but no younger than two years old) and all low-risk individuals older than a given age, while maintaining vaccination of at-risk individuals of any age. All posited strategies were deemed cost-effective. In general, the addition of low-risk vaccination programmes whose coverage encompassed children and young adults (aged 20 and below) were highly cost-effective. The inclusion of elder age-groups to the low-risk programme typically lessened the cost-effectiveness. Notably, elderly-centric programmes vaccinating from 65-75 years and above had the least permitted expense per vaccine.

The 2019-2020 pandemic of atypical pneumonia (COVID-19) caused by the virus SARS-CoV-2 has spread globally and has the potential to infect large numbers of people in every country. Estimating the country-specific basic reproductive ratio is a vital first step in public-health planning. The basic reproductive ratio (R0) is determined by both the nature of pathogen and the network of human contacts through which the disease can spread, which is itself dependent on population age structure and household composition. Here we introduce a transmission model combining age-stratified contact frequencies with age-dependent susceptibility, probability of clinical symptoms, and transmission from asymptomatic (or mild) cases, which we use to estimate the country-specific basic reproductive ratio of COVID-19 for 152 countries. Using early outbreak data from China and a synthetic contact matrix, we estimate an age-stratified transmission structure which can then be extrapolated to 151 other countries for which synthetic contact matrices also exist. This defines a set of country-specific transmission structures from which we can calculate the basic reproductive ratio for each country. Our predicted R0 is critically sensitive to the intensity of transmission from asymptomatic cases; with low asymptomatic transmission the highest values are predicted across Eastern Europe and Japan and the lowest across Africa, Central America and South-Western Asia. This pattern is largely driven by the ratio of children to older adults in each country and the observed propensity of clinical cases in the elderly. If asymptomatic cases have comparable transmission to detected cases, the pattern is reversed. Our results demonstrate the importance of age-specific heterogeneities going beyond contact structure to the spread of COVID-19. These heterogeneities give COVID-19 the capacity to spread particularly quickly in countries with older populations, and that intensive control measures are likely to be necessary to impede its progress in these countries.

Seasonal influenza poses serious problems for global public health, being a significant contributor to morbidity and mortality. In England, there has been a long-standing national vaccination programme, with vaccination of at-risk groups and children offering partial protection against infection. Transmission models have been a fundamental component of analysis, informing the efficient use of limited resources. However, these models generally treat each season and each strain circulating within that season in isolation. Here, we amalgamate multiple data sources to calibrate a susceptible-latent-infected-recovered type transmission model for seasonal influenza, incorporating the four main strains and mechanisms linking prior season epidemiological outcomes to immunity at the beginning of the following season. Data pertaining to nine influenza seasons, starting with the 2009/10 season, informed our estimates for epidemiological processes, virological sample positivity, vaccine uptake and efficacy attributes, and general practitioner influenza-like-illness consultations as reported by the Royal College of General Practitioners (RCGP) Research and Surveillance Centre (RSC). We performed parameter inference via approximate Bayesian computation to assess strain transmissibility, dependence of present season influenza immunity on prior protection, and variability in the influenza case ascertainment across seasons. This produced reasonable agreement between model and data on the annual strain composition. Parameter fits indicated that the propagation of immunity from one season to the next is weaker if vaccine derived, compared to natural immunity from infection. Projecting the dynamics forward in time suggests that while historic immunity plays an important role in determining annual strain composition, the variability in vaccine efficacy hampers our ability to make long-term predictions.

In the event of a new infectious disease outbreak, mathematical and simulation models are commonly used to inform policy by evaluating which control strategies will minimize the impact of the epidemic. In the early stages of such outbreaks, substantial parameter uncertainty may limit the ability of models to provide accurate predictions, and policymakers do not have the luxury of waiting for data to alleviate this state of uncertainty. For policymakers, however, it is the selection of the optimal control intervention in the face of uncertainty, rather than accuracy of model predictions, that is the measure of success that counts. We simulate the process of real-time decision-making by fitting an epidemic model to observed, spatially-explicit, infection data at weekly intervals throughout two historical outbreaks of foot-and-mouth disease, UK in 2001 and Miyazaki, Japan in 2010, and compare forward simulations of the impact of switching to an alternative control intervention at the time point in question. These are compared to policy recommendations generated in hindsight using data from the entire outbreak, thereby comparing the best we could have done at the time with the best we could have done in retrospect. Our results show that the control policy that would have been chosen using all the data is also identified from an early stage in an outbreak using only the available data, despite high variability in projections of epidemic size. Critically, we find that it is an improved understanding of the locations of infected farms, rather than improved estimates of transmission parameters, that drives improved prediction of the relative performance of control interventions. However, the ability to estimate undetected infectious premises is a function of uncertainty in the transmission parameters. Here, we demonstrate the need for both real-time model fitting and generating projections to evaluate alternative control interventions throughout an outbreak. Our results highlight the use of using models at outbreak onset to inform policy and the importance of state-dependent interventions that adapt in response to additional information throughout an outbreak.

Numerical models for simulating outbreaks of infectious diseases are powerful tools for informing surveillance and control strategy decisions. However, large-scale spatially explicit models can be limited by the amount of computational resources they require, which poses a problem when multiple scenarios need to be explored to provide policy recommendations. We introduce an easily implemented method that can reduce computation time in a standard Susceptible-Exposed-Infectious-Removed (SEIR) model without introducing any further approximations or truncations. It is based on a hierarchical infection process that operates on entire groups of spatially related nodes (cells in a grid) in order to efficiently filter out large volumes of susceptible nodes that would otherwise have required expensive calculations. After the filtering of the cells, only a subset of the nodes that were originally at risk are then evaluated for actual infection. The increase in efficiency is sensitive to the exact configuration of the grid, and we describe a simple method to find an estimate of the optimal configuration of a given landscape as well as a method to partition the landscape into a grid configuration. To investigate its efficiency, we compare the introduced methods to other algorithms and evaluate computation time, focusing on simulated outbreaks of foot-and-mouth disease (FMD) on the farm population of the USA, the UK and Sweden, as well as on three randomly generated populations with varying degree of clustering. The introduced method provided up to 500 times faster calculations than pairwise computation, and consistently performed as well or better than other available methods. This enables large scale, spatially explicit simulations such as for the entire continental USA without sacrificing realism or predictive power.

Foot-and-mouth disease outbreaks in non-endemic countries can lead to large economic costs and livestock losses but the use of vaccination has been contentious, partly due to uncertainty about emergency FMD vaccination. Value of information methods can be applied to disease outbreak problems such as FMD in order to investigate the performance improvement from resolving uncertainties. Here we calculate the expected value of resolving uncertainty about vaccine efficacy, time delay to immunity after vaccination and daily vaccination capacity for a hypothetical FMD outbreak in the UK. If it were possible to resolve all uncertainty prior to the introduction of control, we could expect savings of £55 million in outbreak cost, 221,900 livestock culled and 4.3 days of outbreak duration. All vaccination strategies were found to be preferable to a culling only strategy. However, the optimal vaccination radius was found to be highly dependent upon vaccination capacity for all management objectives. We calculate that by resolving the uncertainty surrounding vaccination capacity we would expect to return over 85% of the above savings, regardless of management objective. It may be possible to resolve uncertainty about daily vaccination capacity before an outbreak, and this would enable decision makers to select the optimal control action via careful contingency planning.

Epidemiological modelling has a vital role to play in policy planning and prediction for the control of vectors, and hence the subsequent control of vector-borne diseases. To decide between competing policies requires models that can generate accurate predictions, which in turn requires accurate knowledge of vector natural histories. Here we highlight the importance of the distribution of times between life-history events, using short-lived midge species as an example. In particular we focus on the distribution of the extrinsic incubation period (EIP) which determines the time between infection and becoming infectious, and the distribution of the length of the gonotrophic cycle which determines the time between successful bites. We show how different assumptions for these periods can radically change the basic reproductive ratio (R0) of an infection and additionally the impact of vector control on the infection. These findings highlight the need for detailed entomological data, based on laboratory experiments and field data, to correctly construct the next-generation of policy-informing models.

Network-based infectious disease models have been highly effective in elucidating the role of contact structure in the spread of infection. As such, pair- and neighbourhood-based approximation models have played a key role in linking findings from network simulations to standard (random-mixing) results. Recently, for SIR-type infections (that produce one epidemic in a closed population) on locally tree-like networks, these approximations have been shown to be exact. However, network models are ideally suited for Sexually Transmitted Infections (STIs) due to the greater level of detail available for sexual contact networks, and these diseases often possess SIS-type dynamics. Here, we consider the accuracy of three systematic approximations that can be applied to arbitrary disease dynamics, including SIS behaviour. We focus in particular on low degree networks, in which the small number of neighbours causes build-up of local correlations between the state of adjacent nodes that are challenging to capture. By examining how and when these approximation models converge to simulation results, we generate insights into the role of network structure in the infection dynamics of SIS-type infections.

Prediction and control of the spread of infectious disease in human populations benefits greatly from our growing capacity to quantify human movement behavior. Here we develop a mathematical model for non-transmissible infections contracted from a localized environmental source, informed by a detailed description of movement patterns of the population of Great Britain. The model is applied to outbreaks of Legionnaires' disease, a potentially life-threatening form of pneumonia caused by the bacteria Legionella pneumophilia. We use case-report data from three recent outbreaks that have occurred in Great Britain where the source has already been identified by public health agencies. We first demonstrate that the amount of individual-level heterogeneity incorporated in the movement data greatly influences our ability to predict the source location. The most accurate predictions were obtained using reported travel histories to describe movements of infected individuals, but using detailed simulation models to estimate movement patterns offers an effective fast alternative. Secondly, once the source is identified, we show that our model can be used to accurately determine the population likely to have been exposed to the pathogen, and hence predict the residential locations of infected individuals. The results give rise to an effective control strategy that can be implemented rapidly in response to an outbreak.

The tracing of potentially infectious contacts has become an important part of the control strategy for many infectious diseases, from early cases of novel infections to endemic sexually transmitted infections. Here, we make use of mathematical models to consider the case of partner notification for sexually transmitted infection, however these models are sufficiently simple to allow more general conclusions to be drawn. We show that, when contact network structure is considered in addition to contact tracing, standard "mass action" models are generally inadequate. To consider the impact of mutual contacts (specifically clustering) we develop an improvement to existing pairwise network models, which we use to demonstrate that ceteris paribus, clustering improves the efficacy of contact tracing for a large region of parameter space. This result is sometimes reversed, however, for the case of highly effective contact tracing. We also develop stochastic simulations for comparison, using simple re-wiring methods that allow the generation of appropriate comparator networks. In this way we contribute to the general theory of network-based interventions against infectious disease.