Search

Abstract We propose a method to estimate the household secondary attack rate (hSAR) of COVID-19 in the United Kingdom based on activity on the social media platform X, formerly known as Twitter. Conventional methods of hSAR estimation are resource intensive, requiring regular contact tracing of COVID-19 cases. Our proposed framework provides a complementary method that does not rely on conventional contact tracing or laboratory involvement, including the collection, processing, and analysis of biological samples. We use a text classifier to identify reports of people tweeting about themselves and/or members of their household having COVID-19 infections. A probabilistic analysis is then performed to estimate the hSAR based on the number of self or household, and self and household tweets of COVID-19 infection. The analysis includes adjustments for a reluctance of Twitter users to tweet about household members, and the possibility that the secondary infection was not acquired within the household. Experimental results for the UK, both monthly and weekly, are reported for the period from January 2020 to February 2022. Our results agree with previously reported hSAR estimates, varying with the primary variants of concern, e.g. delta and omicron. The serial interval (SI) is based on the time between the two tweets that indicate a primary and secondary infection. Experimental results, though larger than the consensus, are qualitatively similar. The estimation of hSAR and SI using social media data constitutes a new tool that may help in characterizing, forecasting and managing outbreaks and pandemics in a faster, affordable, and more efficient manner.

Mathematical modelling has played an important role in offering informed advice during the COVID-19 pandemic. In England, a cross government and academia collaboration generated medium-term projections (MTPs) of possible epidemic trajectories over the future 4–6 weeks from a collection of epidemiological models. In this article, we outline this collaborative modelling approach and evaluate the accuracy of the combined and individual model projections against the data over the period November 2021–December 2022 when various Omicron subvariants were spreading across England. Using a number of statistical methods, we quantify the predictive performance of the model projections for both the combined and individual MTPs, by evaluating the point and probabilistic accuracy. Our results illustrate that the combined MTPs, produced from an ensemble of heterogeneous epidemiological models, were a closer fit to the data than the individual models during the periods of epidemic growth or decline, with the 90% confidence intervals widest around the epidemic peaks. We also show that the combined MTPs increase the robustness and reduce the biases associated with a single model projection. Learning from our experience of ensemble modelling during the COVID-19 epidemic, our findings highlight the importance of developing cross-institutional multi-model infectious disease hubs for future outbreak control.

[This corrects the article DOI: 10.1371/journal.ppat.1011461.].

Abstract Background Predicting the likely size of future SARS-CoV-2 waves is necessary for public health planning. In England, voluntary “plan B” mitigation measures were introduced in December 2021 including increased home working and face coverings in shops but stopped short of restrictions on social contacts. The impact of voluntary risk mitigation behaviours on future SARS-CoV-2 burden is unknown. Methods We developed a rapid online survey of risk mitigation behaviours ahead of the winter 2021 festive period and deployed in two longitudinal cohort studies in the UK (Avon Longitudinal Study of Parents and Children (ALSPAC) and TwinsUK/COVID Symptom Study (CSS) Biobank) in December 2021. Using an individual-based, probabilistic model of COVID-19 transmission between social contacts with SARS-CoV-2 Omicron variant parameters and realistic vaccine coverage in England, we predicted the potential impact of the SARS-CoV-2 Omicron wave in England in terms of the effective reproduction number and cumulative infections, hospital admissions and deaths. Using survey results, we estimated in real-time the impact of voluntary risk mitigation behaviours on the Omicron wave in England, if implemented for the entire epidemic wave. Results Over 95% of survey respondents (NALSPAC = 2686 and NTwins = 6155) reported some risk mitigation behaviours, with vaccination and using home testing kits reported most frequently. Less than half of those respondents reported that their behaviour was due to “plan B”. We estimate that without risk mitigation behaviours, the Omicron variant is consistent with an effective reproduction number between 2.5 and 3.5. Due to the reduced vaccine effectiveness against infection with the Omicron variant, our modelled estimates suggest that between 55% and 60% of the English population could be infected during the current wave, translating into between 12,000 and 46,000 cumulative deaths, depending on assumptions about severity and vaccine effectiveness. The actual number of deaths was 15,208 (26 November 2021–1 March 2022). We estimate that voluntary risk reduction measures could reduce the effective reproduction number to between 1.8 and 2.2 and reduce the cumulative number of deaths by up to 24%. Conclusions Predicting future infection burden is affected by uncertainty in disease severity and vaccine effectiveness estimates. In addition to biological uncertainty, we show that voluntary measures substantially reduce the projected impact of the SARS-CoV-2 Omicron variant but that voluntary measures alone would be unlikely to completely control transmission.

In this study, we evaluated the impact of viral variant, in addition to other variables, on within-host viral burden, by analysing cycle threshold (Ct) values derived from nose and throat swabs, collected as part of the UK COVID-19 Infection Survey. Because viral burden distributions determined from community survey data can be biased due to the impact of variant epidemiology on the time-since-infection of samples, we developed a method to explicitly adjust observed Ct value distributions to account for the expected bias. By analysing the adjusted Ct values using partial least squares regression, we found that among unvaccinated individuals with no known prior exposure, viral burden was 44% lower among Alpha variant infections, compared to those with the predecessor strain, B.1.177. Vaccination reduced viral burden by 67%, and among vaccinated individuals, viral burden was 286% higher among Delta variant, compared to Alpha variant, infections. In addition, viral burden increased by 17% for every 10-year age increment of the infected individual. In summary, within-host viral burden increases with age, is reduced by vaccination, and is influenced by the interplay of vaccination status and viral variant.

Abstract Some social settings such as households and workplaces, have been identified as high risk for SARS-CoV-2 transmission. Identifying and quantifying the importance of these settings is critical for designing interventions. A tightly-knit religious community in the UK experienced a very large COVID-19 epidemic in 2020, reaching 64.3% seroprevalence within 10 months, and we surveyed this community both for serological status and individual-level attendance at particular settings. Using these data, and a network model of people and places represented as a stochastic graph rewriting system, we estimated the relative contribution of transmission in households, schools and religious institutions to the epidemic, and the relative risk of infection in each of these settings. All congregate settings were important for transmission, with some such as primary schools and places of worship having a higher share of transmission than others. We found that the model needed a higher general-community transmission rate for women (3.3-fold), and lower susceptibility to infection in children to recreate the observed serological data. The precise share of transmission in each place was related to assumptions about the internal structure of those places. Identification of key settings of transmission can allow public health interventions to be targeted at these locations.

ObjectiveWe aimed to use mathematical models of SARS-COV-2 to assess the potential efficacy of non-pharmaceutical interventions on transmission in the parcel delivery and logistics sector.MethodsWe devloped a network-based model of workplace contacts based on data and consultations from companies in the parcel delivery and logistics sectors. We used these in stochastic simulations of disease transmission to predict the probability of workplace outbreaks in this settings. Individuals in the model have different viral load trajectories based on SARS-CoV-2 in-host dynamics, which couple to their infectiousness and test positive probability over time, in order to determine the impact of testing and isolation measures.ResultsThe baseline model (without any interventions) showed different workplace infection rates for staff in different job roles. Based on our assumptions of contact patterns in the parcel delivery work setting we found that when a delivery driver was the index case, on average they infect only 0.14 other employees, while for warehouse and office workers this went up to 0.65 and 2.24 respectively. In the LIDD setting this was predicted to be 1.40, 0.98, and 1.34 respectively. Nonetheless, the vast majority of simulations resulted in 0 secondary cases among customers (even without contact-free delivery). Our results showed that a combination of social distancing, office staff working from home, and fixed driver pairings (all interventions carried out by the companies we consulted) reduce the risk of workplace outbreaks by 3-4 times.ConclusionThis work suggests that, without interventions, significant transmission could have occured in these workplaces, but that these posed minimal risk to customers. We found that identifying and isolating regular close-contacts of infectious individuals (i.e. house-share, carpools, or delivery pairs) is an efficient measure for stopping workplace outbreaks. Regular testing can make these isolation measures even more effective but also increases the number of staff isolating at one time. It is therefore more efficient to use these isolation measures in addition to social distancing and contact reduction interventions, rather than instead of, as these reduce both transmission and the number of people needing to isolate at one time.

The COVID-19 pandemic has had high mortality rates in the elderly and frail worldwide, particularly in care homes. This is driven by the difficulty of isolating care homes from the wider community, the large population sizes within care facilities (relative to typical households), and the age/frailty of the residents. To quantify the mortality risk posed by disease, the case fatality risk (CFR) is an important tool. This quantifies the proportion of cases that result in death. Throughout the pandemic, CFR amongst care home residents in England has been monitored closely. To estimate CFR, we apply both novel and existing methods to data on deaths in care homes, collected by Public Health England and the Care Quality Commission. We compare these different methods, evaluating their relative strengths and weaknesses. Using these methods, we estimate temporal trends in the instantaneous CFR (at both daily and weekly resolutions) and the overall CFR across the whole of England, and dis-aggregated at regional level. We also investigate how the CFR varies based on age and on the type of care required, dis-aggregating by whether care homes include nursing staff and by age of residents. This work has contributed to the summary of measures used for monitoring the UK epidemic.

Understanding the potential impacts of new variants of SARS-CoV-2 is important for pandemic planning. Here, the authors develop a model incorporating hypothetical new variants with varying transmissibility and immune evasion properties, and use it to project possible future epidemic waves in the UK.

Abstract Background Predicting hospital length of stay (LoS) for patients with COVID-19 infection is essential to ensure that adequate bed capacity can be provided without unnecessarily restricting care for patients with other conditions. Here, we demonstrate the utility of three complementary methods for predicting LoS using UK national- and hospital-level data. Method On a national scale, relevant patients were identified from the COVID-19 Hospitalisation in England Surveillance System (CHESS) reports. An Accelerated Failure Time (AFT) survival model and a truncation corrected method (TC), both with underlying Weibull distributions, were fitted to the data to estimate LoS from hospital admission date to an outcome (death or discharge) and from hospital admission date to Intensive Care Unit (ICU) admission date. In a second approach we fit a multi-state (MS) survival model to data directly from the Manchester University NHS Foundation Trust (MFT). We develop a planning tool that uses LoS estimates from these models to predict bed occupancy. Results All methods produced similar overall estimates of LoS for overall hospital stay, given a patient is not admitted to ICU (8.4, 9.1 and 8.0 days for AFT, TC and MS, respectively). Estimates differ more significantly between the local and national level when considering ICU. National estimates for ICU LoS from AFT and TC were 12.4 and 13.4 days, whereas in local data the MS method produced estimates of 18.9 days. Conclusions Given the complexity and partiality of different data sources and the rapidly evolving nature of the COVID-19 pandemic, it is most appropriate to use multiple analysis methods on multiple datasets. The AFT method accounts for censored cases, but does not allow for simultaneous consideration of different outcomes. The TC method does not include censored cases, instead correcting for truncation in the data, but does consider these different outcomes. The MS method can model complex pathways to different outcomes whilst accounting for censoring, but cannot handle non-random case missingness. Overall, we conclude that data-driven modelling approaches of LoS using these methods is useful in epidemic planning and management, and should be considered for widespread adoption throughout healthcare systems internationally where similar data resources exist.

The first year of the COVID-19 pandemic put considerable strain on healthcare systems worldwide. In order to predict the effect of the local epidemic on hospital capacity in England, we used a variety of data streams to inform the construction and parameterisation of a hospital progression model, EpiBeds, which was coupled to a model of the generalised epidemic. In this model, individuals progress through different pathways (e.g. may recover, die, or progress to intensive care and recover or die) and data from a partially complete patient-pathway line-list was used to provide initial estimates of the mean duration that individuals spend in the different hospital compartments. We then fitted EpiBeds using complete data on hospital occupancy and hospital deaths, enabling estimation of the proportion of individuals that follow the different clinical pathways, the reproduction number of the generalised epidemic, and to make short-term predictions of hospital bed demand. The construction of EpiBeds makes it straightforward to adapt to different patient pathways and settings beyond England. As part of the UK response to the pandemic, EpiBeds provided weekly forecasts to the NHS for hospital bed occupancy and admissions in England, Wales, Scotland, and Northern Ireland at national and regional scales.

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.

The use of data has been essential throughout the unfolding COVID-19 pandemic. We have needed it to populate our models, inform our understanding, and shape our responses to the disease. However, data has not always been easy to find and access, it has varied in quality and coverage, been difficult to reuse or repurpose. This paper reviews these and other challenges and recommends steps to develop a data ecosystem better able to deal with future pandemics by better supporting preparedness, prevention, detection and response.

New disease challenges, societal demands and better or novel types of data, drive innovations in the structure, formulation and analysis of epidemic models. Innovations in modelling can lead to new insights into epidemic processes and better use of available data, yielding improved disease control and stimulating collection of better data and new data types. Here we identify key challenges for the structure, formulation, analysis and use of mathematical models of pathogen transmission relevant to current and future pandemics.

Models of emerging epidemics can be exceedingly helpful in planning the response, but early on model selection is a difficult task. Here, the authors explore the joint contribution of age stratification and household structure on epidemic spread, and provides a rule of thumb to guide model choice.

During an infectious disease outbreak, biases in the data and complexities of the underlying dynamics pose significant challenges in mathematically modelling the outbreak and designing policy. Motivated by the ongoing response to COVID-19, we provide a toolkit of statistical and mathematical models beyond the simple SIR-type differential equation models for analysing the early stages of an outbreak and assessing interventions. In particular, we focus on parameter estimation in the presence of known biases in the data, and the effect of non-pharmaceutical interventions in enclosed subpopulations, such as households and care homes. We illustrate these methods by applying them to the COVID-19 pandemic.

Background: After the declaration of COVID-19 pandemic on March 11th, 2020, local transmission chains starting in different countries including Canada are forcing governments to take decisions on public health interventions to mitigate the spread of the epidemic. Methods: We conduct data-driven and model-free estimations for the growth rates of the COVID-19 epidemics in Italy and Canada, by fitting an exponential curve to the daily reported cases. We use these estimates to predict epidemic trends in Canada under different scenarios of public health interventions. Results: In Italy, the initial growth rate (0.22) has reduced to 0.1 two weeks after the lockdown of the country on March 8th, 2020. This corresponds to an increase of the doubling time from about 3.15 to almost 7 days. In comparison, the growth rate in Canada has increased from 0.13 between March 1st and 13th, to 0.25 between March 13th to 22nd. This current growth rate corresponds to a doubling time of 2.7 days, and therefore, unless further public health interventions are escalated in Canada, we project 15,000 cases by March 31st. However, the case number may be reduced to 4000 if escalated public health interventions could instantly reduce the growth rate to 0.1, the same level achieved in Italy. Interpretation: Prompt and farsighted interventions are critical to counteract the very rapid initial growth of the COVID-19 epidemic in Canada. Mitigation plans must take into account the delayed effect of interventions by up to 2-weeks and the short doubling time of 3–4 days.

The estimation of parameters and model structure for informing infectious disease response has become a focal point of the recent pandemic. However, it has also highlighted a plethora of challenges remaining in the fast and robust extraction of information using data and models to help inform policy. In this paper, we identify and discuss four broad challenges in the estimation paradigm relating to infectious disease modelling, namely the Uncertainty Quantification framework, data challenges in estimation, model-based inference and prediction, and expert judgement. We also postulate priorities in estimation methodology to facilitate preparation for future pandemics.

Mathematical modelling and statistical inference provide a framework to evaluate different non-pharmaceutical and pharmaceutical interventions for the control of epidemics that has been widely used during the COVID-19 pandemic. In this paper, lessons learned from this and previous epidemics are used to highlight the challenges for future pandemic control. We consider the availability and use of data, as well as the need for correct parameterisation and calibration for different model frameworks. We discuss challenges that arise in describing and distinguishing between different interventions, within different modelling structures, and allowing both within and between host dynamics. We also highlight challenges in modelling the health economic and political aspects of interventions. Given the diversity of these challenges, a broad variety of interdisciplinary expertise is needed to address them, combining mathematical knowledge with biological and social insights, and including health economics and communication skills. Addressing these challenges for the future requires strong cross-disciplinary collaboration together with close communication between scientists and policy makers.

Objectives Importations of novel variants of concern (VOC), particularly B.1.617.2, have become the impetus behind recent outbreaks of SARS-CoV-2. Concerns around the impact on vaccine effectiveness, transmissibility and severity are now driving the public health response to these variants. This paper analyses the patterns of growth in hospitalisations and confirmed cases for novel VOCs by age groups, geography and ethnicity in the context of changing behaviour, non-pharmaceutical interventions (NPIs) and the UK vaccination programme. We seek to highlight where strategies have been effective and periods that have facilitated the establishment of new variants.Design We have algorithmically linked the most complete testing and hospitalisation data in England to create a data set of confirmed infections and hospitalisations by SARS-CoV-2 genomic variant. We have used these linked data sets to analyse temporal, geographic and demographic distinctions.Setting and participants The setting is England from October 2020 to July 2021. Participants included all COVID-19 tests that included RT-PCR CT gene target data or underwent sequencing and hospitalisations that could be linked to these tests.Methods To calculate the instantaneous growth rate for VOCs we have developed a generalised additive model fit to multiple splines and varying day of the week effects. We have further modelled the instantaneous reproduction number Rt for the B.1.1.7 and B.1.617.2 variants and included a doubly interval censored model to temporally adjust the confirmed variant cases.Results We observed a clear replacement of the predominant B.1.1.7 by the B.1.617.2 variant without observing sustained exponential growth in other novel variants. Modelled exponential growth of RT PCR gene target triple-positive cases was initially detected in the youngest age groups, although we now observe across all ages a very small doubling time of 10.7 (95% CI 9.1 to 13.2) days and 8 (95% CI 6.9 to 9.1) days for cases and hospitalisations, respectively. We observe that growth in RT PCR gene target triple-positive cases was first detected in the Indian ethnicity group in late February, with a peak of 0.06 (95% CI 0.07 to 0.05) in the instantaneous growth rate, but is now maintained by the white ethnicity groups, observing a doubling time of 6.8 (95% CI 4.9 to 11) days. Rt analysis indicates a reproduction number advantage of 0.45 for B.1.617.2 relative to B.1.1.7, with the Rt value peaking at 1.85 for B.1.617.2.Conclusions Our results illustrate a clear transmission advantage for the B.1.617.2 variant and the growth in hospitalisations illustrates that this variant is able to maintain exponential growth within age groups that are largely doubly vaccinated. There are concerning signs of intermittent growth in the B.1.351 variant, reaching a 28-day doubling time peak in March 2021, although this variant is presently not showing any evidence of a transmission advantage over B.1.617.2. Step 1b of the UK national lockdown easing was sufficient to precipitate exponential growth in B.1.617.2 cases for most regions and younger adult age groups. The final stages of NPI easing appeared to have a negligible impact on the growth of B.1.617.2 with every region experiencing sustained exponential growth from step 2. Nonetheless, early targeted local NPIs appeared to markedly reduced growth of B.1.617.2. Later localised interventions, at a time of higher prevalence and greater geographic dispersion of this variant, appeared to have a negligible impact on growth.

We propose a deterministic model capturing essential features of contact tracing as part of public health non-pharmaceutical interventions to mitigate an outbreak of an infectious disease. By incorporating a mechanistic formulation of the processes at the individual level, we obtain an integral equation (delayed in calendar time and advanced in time since infection) for the probability that an infected individual is detected and isolated at any point in time. This is then coupled with a renewal equation for the total incidence to form a closed system describing the transmission dynamics involving contact tracing. We define and calculate basic and effective reproduction numbers in terms of pathogen characteristics and contact tracing implementation constraints. When applied to the case of SARS-CoV-2, our results show that only combinations of diagnosis of symptomatic infections and contact tracing that are almost perfect in terms of speed and coverage can attain control, unless additional measures to reduce overall community transmission are in place. Under constraints on the testing or tracing capacity, a temporary interruption of contact tracing may, depending on the overall growth rate and prevalence of the infection, lead to an irreversible loss of control even when the epidemic was previously contained.

Pathogens have evolved a variety of life-history strategies. An important strategy consists of successful transmission by an infected host before the appearance of symptoms, that is, while the host is still partially or fully asymptomatic. During this initial stage of infection, it is possible for another pathogen to superinfect an already infected host and replace the previously infecting pathogen. Here, we study the effect of superinfection during the first stage of an infection on the evolutionary dynamics of the degree to which the host is asymptomatic (host latency) in that same stage. We find that superinfection can lead to major differences in evolutionary behaviour. Most strikingly, the duration of immunity following infection can significantly influence pathogen evolutionary dynamics, whereas without superinfection the outcomes are independent of host immunity. For example, changes in host immunity can drive evolutionary transitions from a fully symptomatic to a fully asymptomatic first infection stage. Additionally, if superinfection relative to susceptible infection is strong enough, evolution can lead to a unique strategy of latency that corresponds to a local fitness minimum, and is therefore invasible by nearby mutants. Thus, this strategy is a branching point, and can lead to coexistence of pathogens with different latencies. Furthermore, in this new framework with superinfection, we also find that there can exist two interior singular strategies. Overall, new evolutionary outcomes can cascade from superinfection.

Background: ​ During the coronavirus disease 2019 (COVID-19) lockdown, contact clustering in social bubbles may allow extending contacts beyond the household at minimal additional risk and hence has been considered as part of modified lockdown policy or a gradual lockdown exit strategy. We estimated the impact of such strategies on epidemic and mortality risk using the UK as a case study. Methods: ​ We used an individual based model for a synthetic population similar to the UK, stratified into transmission risks from the community, within the household and from other households in the same social bubble. The base case considers a situation where non-essential shops and schools are closed, the secondary household attack rate is 20% and the initial reproduction number is 0.8. We simulate social bubble strategies (where two households form an exclusive pair) for households including children, for single occupancy households, and for all households. We test the sensitivity of results to a range of alternative model assumptions and parameters. Results: Clustering contacts outside the household into exclusive bubbles is an effective strategy of increasing contacts while limiting the associated increase in epidemic risk. In the base case, social bubbles reduced fatalities by 42% compared to an unclustered increase of contacts. We find that if all households were to form social bubbles the reproduction number would likely increase to above the epidemic threshold of R=1. Strategies allowing households with young children or single occupancy households to form social bubbles increased the reproduction number by less than 11%. The corresponding increase in mortality is proportional to the increase in the epidemic risk but is focussed in older adults irrespective of inclusion in social bubbles. Conclusions: ​ If managed appropriately, social bubbles can be an effective way of extending contacts beyond the household while limiting the increase in epidemic risk.

BACKGROUND:Visceral leishmaniasis (VL) is characterised by a high degree of spatial clustering at all scales, and this feature remains even with successful control measures. VL is targeted for elimination as a public health problem in the Indian subcontinent by 2020, and incidence has been falling rapidly since 2011. Current control is based on early diagnosis and treatment of clinical cases, and blanket indoor residual spraying of insecticide (IRS) in endemic villages to kill the sandfly vectors. Spatially targeting active case detection and/or IRS to higher risk areas would greatly reduce costs of control, but its effectiveness as a control strategy is unknown. The effectiveness depends on two key unknowns: how quickly transmission risk decreases with distance from a VL case and how much asymptomatically infected individuals contribute to transmission. METHODOLOGY/PRINCIPAL FINDINGS:To estimate these key parameters, a spatiotemporal transmission model for VL was developed and fitted to geo-located epidemiological data on 2494 individuals from a highly endemic village in Mymensingh, Bangladesh. A Bayesian inference framework that could account for the unknown infection times of the VL cases, and missing symptom onset and recovery times, was developed to perform the parameter estimation. The parameter estimates obtained suggest that, in a highly endemic setting, VL risk decreases relatively quickly with distance from a case-halving within 90m-and that VL cases contribute significantly more to transmission than asymptomatic individuals. CONCLUSIONS/SIGNIFICANCE:These results suggest that spatially-targeted interventions may be effective for limiting transmission. However, the extent to which spatial transmission patterns and the asymptomatic contribution vary with VL endemicity and over time is uncertain. In any event, interventions would need to be performed promptly and in a large radius (≥300m) around a new case to reduce transmission risk.

Networks offer a fertile framework for studying the spread of infection in human and animal populations. However, owing to the inherent high-dimensionality of networks themselves, modelling transmission through networks is mathematically and computationally challenging. Even the simplest network epidemic models present unanswered questions. Attempts to improve the practical usefulness of network models by including realistic features of contact networks and of host–pathogen biology (e.g. waning immunity) have made some progress, but robust analytical results remain scarce. A more general theory is needed to understand the impact of network structure on the dynamics and control of infection. Here we identify a set of challenges that provide scope for active research in the field of network epidemic models.

This paper considers metapopulation models in the general sense, i.e. where the population is partitioned into sub-populations (groups, patches,...), irrespective of the biological interpretation they have, e.g. spatially segregated large sub-populations, small households or hosts themselves modelled as populations of pathogens. This framework has traditionally provided an attractive approach to incorporating more realistic contact structure into epidemic models, since it often preserves analytic tractability (in stochastic as well as deterministic models) but also captures the most salient structural inhomogeneity in contact patterns in many applied contexts. Despite the progress that has been made in both the theory and application of such metapopulation models, we present here several major challenges that remain for future work, focusing on models that, in contrast to agent-based ones, are amenable to mathematical analysis. The challenges range from clarifying the usefulness of systems of weakly-coupled large sub-populations in modelling the spread of specific diseases to developing a theory for endemic models with household structure. They include also developing inferential methods for data on the emerging phase of epidemics, extending metapopulation models to more complex forms of human social structure, developing metapopulation models to reflect spatial population structure, developing computationally efficient methods for calculating key epidemiological model quantities, and integrating within- and between-host dynamics in models.

Deterministic models have a long history of being applied to the study of infectious disease epidemiology. We highlight and discuss nine challenges in this area. The first two concern the endemic equilibrium and its stability. We indicate the need for models that describe multi-strain infections, infections with time-varying infectivity, and those where superinfection is possible. We then consider the need for advances in spatial epidemic models, and draw attention to the lack of models that explore the relationship between communicable and non-communicable diseases. The final two challenges concern the uses and limitations of deterministic models as approximations to stochastic systems.

The population dynamics of infectious disease is a mature field in terms of theory and to some extent, application. However for microparasites, the theory and application of models of the dynamics within a single infected host is still an open field. Further, connecting across the scales – from cellular to host level, to population level – has potential to vastly improve our understanding of pathogen dynamics and evolution. Here, we highlight seven challenges in the following areas: transmission bottlenecks, heterogeneity within host, dynamic fitness landscapes within hosts, making use of next-generation sequencing data, capturing superinfection and when and how to model more than two scales.

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.

The viral population of HIV-1, like many pathogens that cause systemic infection, is structured and differentiated within the body. The dynamics of cellular immune trafficking through the blood and within compartments of the body has also received wide attention. Despite these advances, mathematical models, which are widely used to interpret and predict viral and immune dynamics in infection, typically treat the infected host as a well-mixed homogeneous environment. Here, we present mathematical, analytical, and computational results that demonstrate that consideration of the spatial structure of the viral population within the host radically alters predictions of previous models. We study the dynamics of virus replication and cytotoxic T lymphocytes (CTLs) within a metapopulation of spatially segregated patches, representing T cell areas connected by circulating blood and lymph. The dynamics of the system depend critically on the interaction between CTLs and infected cells at the within-patch level. We show that for a wide range of parameters, the system admits an unexpected outcome called the shifting-mosaic steady state. In this state, the whole body's viral population is stable over time, but the equilibrium results from an underlying, highly dynamic process of local infection and clearance within T-cell centers. Notably, and in contrast to previous models, this new model can explain the large differences in set-point viral load (SPVL) observed between patients and their distribution, as well as the relatively low proportion of cells infected at any one time, and alters the predicted determinants of viral load variation.