Search

The role of direct cell-to-cell spread in viral infections-where virions spread between host and susceptible cells without needing to be secreted into the extracellular environment-has come to be understood as essential to the dynamics of medically significant viruses like hepatitis C and influenza. Recent work in both the experimental and mathematical modelling literature has attempted to quantify the prevalence of cell-to-cell infection compared to the conventional free virus route using a variety of methods and experimental data. However, estimates are subject to significant uncertainty and moreover rely on data collected by inhibiting one mode of infection by either chemical or physical factors, which may influence the other mode of infection to an extent which is difficult to quantify. In this work, we conduct a simulation-estimation study to probe the practical identifiability of the proportion of cell-to-cell infection, using two standard mathematical models and synthetic data that would likely be realistic to obtain in the laboratory. We show that this quantity cannot be estimated using non-spatial data alone, and that the collection of data which describes the spatial structure of the infection is necessary to infer the proportion of cell-to-cell infection. Our results provide guidance for the design of relevant experiments and mathematical tools for accurately inferring the prevalence of cell-to-cell infection in in vitro and in vivo contexts.

A characteristic of malaria in all its forms is the potential for superinfection (that is, multiple concurrent blood-stage infections). An additional characteristic of Plasmodium vivax malaria is a reservoir of latent parasites (hypnozoites) within the host liver, which activate to cause (blood-stage) relapses. Here, we present a model of hypnozoite accrual and superinfection for P. vivax. To couple host and vector dynamics for a homogeneously-mixing population, we construct a density-dependent Markov population process with countably many types, for which disease extinction is shown to occur almost surely. We also establish a functional law of large numbers, taking the form of an infinite-dimensional system of ordinary differential equations that can also be recovered by coupling expected host and vector dynamics (i.e. a hybrid approximation) or through a standard compartment modelling approach. Recognising that the subset of these equations that model the infection status of the human hosts has precisely the same form as the Kolmogorov forward equations for a Markovian network of infinite server queues with an inhomogeneous batch arrival process, we use physical insight into the evolution of the latter process to write down a time-dependent multivariate generating function for the solution. We use this characterisation to collapse the infinite-compartment model into a single integrodifferential equation (IDE) governing the intensity of mosquito-to-human transmission. Through a steady state analysis, we recover a threshold phenomenon for this IDE in terms of a parameter [Formula: see text] expressible in terms of the primitives of the model, with the disease-free equilibrium shown to be uniformly asymptotically stable if [Formula: see text] and an endemic equilibrium solution emerging if [Formula: see text]. Our work provides a theoretical basis to explore the epidemiology of P. vivax, and introduces a strategy for constructing tractable population-level models of

The availability of COVID-19 vaccines promised a reduction in the severity of disease and relief from the strict public health and social measures (PHSMs) imposed in many countries to limit spread and burden of COVID-19. We were asked to define vaccine coverage thresholds for Australia’s transition to easing restrictions and reopening international borders. Using evidence of vaccine effectiveness against the then-circulating Delta variant, we used a mathematical model to determine coverage targets. The absence of any COVID-19 infections in many sub-national jurisdictions in Australia posed particular methodological challenges. We used a novel metric called Transmission Potential (TP) as a proxy measure of the population-level effective reproduction number. We estimated TP of the Delta variant under a range of PHSMs, test-trace-isolate-quarantine (TTIQ) efficiencies, vaccination coverage thresholds, and age-based vaccine allocation strategies. We found that high coverage across all ages (>_70%) combined with ongoing TTIQ and minimal PHSMs was sufficient to avoid lockdowns. At lesser coverage (<_60%) rapid case escalation risked overwhelming of the health sector or the need to reimpose stricter restrictions. Maintaining low case numbers was most beneficial for health and the economy, and at higher coverage levels (>_80%) further easing of restrictions was deemed possible. These results directly informed easing of COVID-19 restrictions in Australia.

Disease surveillance data was critical in supporting public health decisions throughout the coronavirus disease 2019 (COVID-19) pandemic. At the same time, the unprecedented circumstances of the pandemic revealed many shortcomings of surveillance systems for viral respiratory pathogens. Strengthening of surveillance systems was identified as a priority for the recently established Australian Centre for Disease Control, which represents a critical opportunity to review pre-pandemic and pandemic surveillance practices, and to decide on future priorities, during both pandemic and inter-pandemic periods. On 20 October 2022, we ran a workshop with experts from the academic and government sectors who had contributed to the COVID-19 response in Australia on 'The role of surveillance in epidemic response', at the University of New South Wales, Sydney, Australia. Following the workshop, we developed five recommendations to strengthen respiratory virus surveillance systems in Australia, which we present here. Our recommendations are not intended to be exhaustive. We instead chose to focus on data types that are highly valuable yet typically overlooked by surveillance planners. Three of the recommendations focus on data collection activities that support the monitoring and prediction of disease impact and the effectiveness of interventions (what to measure) and two focus on surveillance methods and capabilities (how to measure). Implementation of our recommendations would enable more robust, timely, and impactful epidemic analysis.

Plasmodium vivax is one of the most geographically widespread malaria parasites in the world, primarily found across South-East Asia, Latin America, and parts of Africa. One of the significant characteristics of the P. vivax parasite is its ability to remain dormant in the human liver as hypnozoites and subsequently reactivate after the initial infection (i.e. relapse infections). Mathematical modelling approaches have been widely applied to understand P. vivax dynamics and predict the impact of intervention outcomes. Models that capture P. vivax dynamics differ from those that capture P. falciparum dynamics, as they must account for relapses caused by the activation of hypnozoites. In this article, we provide a scoping review of mathematical models that capture P. vivax transmission dynamics published between January 1988 and May 2023. The primary objective of this work is to provide a comprehensive summary of the mathematical models and techniques used to model P. vivax dynamics. In doing so, we aim to assist researchers working on mathematical epidemiology, disease transmission, and other aspects of P. vivax malaria by highlighting best practices in currently published models and highlighting where further model development is required. We categorise P. vivax models according to whether a deterministic or agent-based approach was used. We provide an overview of the different strategies used to incorporate the parasite's biology, use of multiple scales (within-host and population-level), superinfection, immunity, and treatment interventions. In most of the published literature, the rationale for different modelling approaches was driven by the research question at hand. Some models focus on the parasites' complicated biology, while others incorporate simplified assumptions to avoid model complexity. Overall, the existing literature on mathematical models for P. vivax encompasses various aspects of the parasite's dynamics. We recommend that future research should f

Plasmodium vivax lactate dehydrogenase (PvLDH) is an essential enzyme in the glycolytic pathway of P. vivax. It is widely used as a diagnostic biomarker and a measure of total-body parasite biomass in vivax malaria. However, the dynamics of PvLDH remains poorly understood. Here, we developed mathematical models that capture parasite and matrix PvLDH dynamics in ex vivo culture and the human host. We estimated key biological parameters characterising in vivo PvLDH dynamics based on longitudinal data of parasitemia and PvLDH concentration collected from P. vivax-infected humans, with the estimates informed by the ex vivo data as prior knowledge in a Bayesian hierarchical framework. We found that the in vivo accumulation rate of intraerythrocytic PvLDH peaks at 10-20 h post-invasion (late ring stage) with a median estimate of intraerythrocytic PvLDH mass at the end of the life cycle to be 9.4 × 10-3ng. We also found that the median estimate of in vivo PvLDH half-life was approximately 21.9 h. Our findings provide a foundation with which to advance our quantitative understanding of P. vivax biology and will facilitate the improvement of PvLDH-based diagnostic tools.

BACKGROUND: Australian states and territories used test-trace-isolate-quarantine (TTIQ) systems extensively in their response to the COVID-19 pandemic in 2020-2021. We report on an analysis of Australian case data to estimate the impact of test-trace-isolate-quarantine systems on SARS-CoV-2 transmission. METHODS: Our analysis uses a novel mathematical modelling framework and detailed surveillance data on COVID-19 cases including dates of infection and dates of isolation. First, we directly translate an empirical distribution of times from infection to isolation into reductions in potential for onward transmission during periods of relatively low caseloads (tens to hundreds of reported cases per day). We then apply a simulation approach, validated against case data, to assess the impact of case-initiated contact tracing on transmission during a period of relatively higher caseloads and system stress (up to thousands of cases per day). RESULTS: We estimate that under relatively low caseloads in the state of New South Wales (tens of cases per day), TTIQ contributed to a 54% reduction in transmission. Under higher caseloads in the state of Victoria (hundreds of cases per day), TTIQ contributed to a 42% reduction in transmission. Our results also suggest that case-initiated contact tracing can support timely quarantine in times of system stress (thousands of cases per day). CONCLUSION: Contact tracing systems for COVID-19 in Australia were highly effective and adaptable in supporting the national suppression strategy from 2020-21, prior to the emergence of the Omicron variant in November 2021. TTIQ systems were critical to the maintenance of the strong suppression strategy and were more effective when caseloads were (relatively) low.

Since the emergence of SARS-CoV-2 in 2019 through to mid-2021, much of the Australian population lived in a COVID-19-free environment. This followed the broadly successful implementation of a strong suppression strategy, including international border closures. With the availability of COVID-19 vaccines in early 2021, the national government sought to transition from a state of minimal incidence and strong suppression activities to one of high vaccine coverage and reduced restrictions but with still-manageable transmission. This transition is articulated in the national 're-opening' plan released in July 2021. Here, we report on the dynamic modelling study that directly informed policies within the national re-opening plan including the identification of priority age groups for vaccination, target vaccine coverage thresholds and the anticipated requirements for continued public health measures-assuming circulation of the Delta SARS-CoV-2 variant. Our findings demonstrated that adult vaccine coverage needed to be at least 60% to minimize public health and clinical impacts following the establishment of community transmission. They also supported the need for continued application of test-trace-isolate-quarantine and social measures during the vaccine roll-out phase and beyond.

There has been an increasing recognition of the utility of models of the spatial dynamics of viral spread within tissues. Multicellular models, where cells are represented as discrete regions of space coupled to a virus density surface, are a popular approach to capture these dynamics. Conventionally, such models are simulated by discretising the viral surface and depending on the rate of viral diffusion and other considerations, a finer or coarser discretisation may be used. The impact that this choice may have on the behaviour of the system has not been studied. Here we demonstrate that under realistic parameter regimes - where viral diffusion is small enough to support the formation of familiar ring-shaped infection plaques - the choice of spatial discretisation of the viral surface can qualitatively change key model outcomes including the time scale of infection. Importantly, we show that the choice between implementing viral spread as a cell-scale process, or as a high-resolution converged PDE can generate distinct model outcomes, which raises important conceptual questions about the strength of assumptions underpinning the spatial structure of the model. We investigate the mechanisms driving these discretisation artefacts, the impacts they may have on model predictions, and provide guidance on the design and implementation of spatial and especially multicellular models of viral dynamics. We obtain our results using the simplest TIV construct for the viral dynamics, and therefore anticipate that the important effects we describe will also influence model predictions in more complex models of virus-cell-immune system interactions. This analysis will aid in the construction of models for robust and biologically realistic modelling and inference.

The ability for vaccines to protect against infectious diseases varies among individuals, but computational models employed to inform policy typically do not account for this variation. Here we examine this issue: we implement a model of vaccine efficacy developed in the context of SARS-CoV-2 in order to evaluate the general implications of modelling correlates of protection on the individual level. Due to high levels of variation in immune response, the distributions of individual-level protection emerging from this model tend to be highly dispersed, and are often bimodal. We describe the specification of the model, provide an intuitive parameterisation, and comment on its general robustness. We show that the model can be viewed as an intermediate between the typical approaches that consider the mode of vaccine action to be either "all-or-nothing" or "leaky". Our view based on this analysis is that individual variation in correlates of protection is an important consideration that may be crucial to designing and implementing models for estimating population-level impacts of vaccination programs.

As of January 2021, Australia had effectively controlled local transmission of COVID-19 despite a steady influx of imported cases and several local, but contained, outbreaks in 2020. Throughout 2020, state and territory public health responses were informed by weekly situational reports that included an ensemble forecast of daily COVID-19 cases for each jurisdiction. We present here an analysis of one forecasting model included in this ensemble across the variety of scenarios experienced by each jurisdiction from May to October 2020. We examine how successfully the forecasts characterised future case incidence, subject to variations in data timeliness and completeness, showcase how we adapted these forecasts to support decisions of public health priority in rapidly-evolving situations, evaluate the impact of key model features on forecast skill, and demonstrate how to assess forecast skill in real-time before the ground truth is known. Conditioning the model on the most recent, but incomplete, data improved the forecast skill, emphasising the importance of developing strong quantitative models of surveillance system characteristics, such as ascertainment delay distributions. Forecast skill was highest when there were at least 10 reported cases per day, the circumstances in which authorities were most in need of forecasts to aid in planning and response.

Epidemiological and clinical evidence indicates that humans infected with the 1918 pandemic H1N1 influenza virus and highly pathogenic avian H5N1 influenza viruses often displayed severe lung pathology. High viral load and extensive infiltration of macrophages are the hallmarks of highly pathogenic (HP) influenza viral infections. However, it remains unclear what biological mechanisms primarily determine the observed difference in the kinetics of viral load and macrophages between HP and low pathogenic (LP) viral infections, and how the mechanistic differences are associated with viral pathogenicity. In this study, we develop a mathematical model of viral dynamics that includes the dynamics of different macrophage populations and interferon. We fit the model to in vivo kinetic data of viral load and macrophage level from BALB/c mice infected with an HP or LP strain of H1N1/H5N1 virus to estimate model parameters using Bayesian inference. Our primary finding is that HP viruses have a higher viral infection rate, a lower interferon production rate and a lower macrophage recruitment rate compared to LP viruses, which are strongly associated with more severe tissue damage (quantified by a higher percentage of epithelial cell loss). We also quantify the relative contribution of macrophages to viral clearance and find that macrophages do not play a dominant role in the direct clearance of free viruses although their role in mediating immune responses such as interferon production is crucial. Our work provides new insight into the mechanisms that convey the observed difference in viral and macrophage kinetics between HP and LP infections and establishes an improved model-fitting framework to enhance the analysis of new data on viral pathogenicity.

Plasmodium vivax is the most geographically widespread malaria-causing parasite resulting in significant associated global morbidity and mortality. One of the factors driving this widespread phenomenon is the ability of the parasites to remain dormant in the liver. Known as 'hypnozoites', they reside in the liver following an initial exposure, before activating later to cause further infections, referred to as 'relapses'. As around 79-96% of infections are attributed to relapses from activating hypnozoites, we expect it will be highly impactful to apply treatment to target the hypnozoite reservoir (i.e. the collection of dormant parasites) to eliminate P. vivax. Treatment with radical cure, for example tafenoquine or primaquine, to target the hypnozoite reservoir is a potential tool to control and/or eliminate P. vivax. We have developed a deterministic multiscale mathematical model as a system of integro-differential equations that captures the complex dynamics of P. vivax hypnozoites and the effect of hypnozoite relapse on disease transmission. Here, we use our multiscale model to study the anticipated effect of radical cure treatment administered via a mass drug administration (MDA) program. We implement multiple rounds of MDA with a fixed interval between rounds, starting from different steady-state disease prevalences. We then construct an optimisation model with three different objective functions motivated on a public health basis to obtain the optimal MDA interval. We also incorporate mosquito seasonality in our model to study its effect on the optimal treatment regime. We find that the effect of MDA interventions is temporary and depends on the pre-intervention disease prevalence (and choice of model parameters) as well as the number of MDA rounds under consideration. The optimal interval between MDA rounds also depends on the objective (combinations of expected intervention outcomes). We find radical cure alone may not be enough to lead to P. vivax elimina

Against a backdrop of widespread global transmission, a number of countries have successfully brought large outbreaks of COVID-19 under control and maintained near-elimination status. A key element of epidemic response is the tracking of disease transmissibility in near real-time. During major outbreaks, the effective reproduction number can be estimated from a time-series of case, hospitalisation or death counts. In low or zero incidence settings, knowing the potential for the virus to spread is a response priority. Absence of case data means that this potential cannot be estimated directly. We present a semi-mechanistic modelling framework that draws on time-series of both behavioural data and case data (when disease activity is present) to estimate the transmissibility of SARS-CoV-2 from periods of high to low - or zero - case incidence, with a coherent transition in interpretation across the changing epidemiological situations. Of note, during periods of epidemic activity, our analysis recovers the effective reproduction number, while during periods of low - or zero - case incidence, it provides an estimate of transmission risk. This enables tracking and planning of progress towards the control of large outbreaks, maintenance of virus suppression, and monitoring the risk posed by re-introduction of the virus. We demonstrate the value of our methods by reporting on their use throughout 2020 in Australia, where they have become a central component of the national COVID-19 response.

Plasmodium falciparum and P. vivax are the two most common causes of malaria. While the majority of deaths and severe morbidity are due to P. falciparum, P. vivax poses a greater challenge to eliminating malaria outside of Africa due to its ability to form latent liver stage parasites (hypnozoites), which can cause relapsing episodes within an individual patient. In areas where P. falciparum and P. vivax are co-endemic, individuals can carry parasites of both species simultaneously. These mixed infections complicate dynamics in several ways: treatment of mixed infections will simultaneously affect both species, P. falciparum can mask the detection of P. vivax, and it has been hypothesised that clearing P. falciparum may trigger a relapse of dormant P. vivax. When mixed infections are treated for only blood-stage parasites, patients are at risk of relapse infections due to P. vivax hypnozoites. We present a stochastic mathematical model that captures interactions between P. falciparum and P. vivax, and incorporates both standard schizonticidal treatment (which targets blood-stage parasites) and radical cure treatment (which additionally targets liver-stage parasites). We apply this model via a hypothetical simulation study to assess the implications of different treatment coverages of radical cure for mixed and P. vivax infections and a "unified radical cure" treatment strategy where P. falciparum, P. vivax, and mixed infections all receive radical cure after screening glucose-6-phosphate dehydrogenase (G6PD) normal. In addition, we investigated the impact of mass drug administration (MDA) of blood-stage treatment. We find that a unified radical cure strategy leads to a substantially lower incidence of malaria cases and deaths overall. MDA with schizonticidal treatment was found to decrease P. falciparum with little effect on P. vivax. We perform a univariate sensitivity analysis to highlight important model parameters.

Deterministic epidemic models that allow for replenishment of susceptibles typically display damped oscillatory behaviour. If the population is initially fully susceptible, once an epidemic takes off a distinct trough will exist between the first and second waves of infection. Epidemic dynamics are, however, influenced by stochastic effects, particularly when the prevalence is low. At the beginning of an epidemic, stochastic die-out is possible and well characterised through use of a branching process approximation. Conditional on an epidemic taking off, stochastic extinction is highly unlikely during the first epidemic wave, but the probability of extinction increases again as the wave declines. Extinction during this period, prior to a potential second wave of infection, is defined as 'epidemic fade-out'. We consider a set of observed epidemics, each distinct and having evolved independently, in which some display fade-out and some do not. While fade-out is necessarily a stochastic phenomenon, the probability of fade-out will depend on the model parameters associated with each epidemic. Accordingly, we ask whether time-series data for the epidemics contain sufficient information to identify the key driver(s) of different outcomes-fade-out or otherwise-across the sub-populations supporting each epidemic. We apply a Bayesian hierarchical modelling framework to synthetic data from an SIRS model of epidemic dynamics and demonstrate that we can (1) identify when the sub-population specific model parameters supporting each epidemic have significant variability and (2) estimate the probability of epidemic fade-out for each sub-population. We demonstrate that a hierarchical analysis can provide precise estimates of the probability of fade-out than is possible if considering each epidemic in isolation. Our methods may be applied to both epidemiological and other biological data to identify where differences in outcome-fade-out or recurrent infection/waves are purely due to

Scientific knowledge and advances are a cornerstone of modern society. They improve our understanding of the world we live in and help us navigate global challenges including emerging infectious diseases, climate change and the biodiversity crisis. However, there is a perpetual challenge in translating scientific insight into policy. Many articles explain how to better bridge the gap through improved communication and engagement, but we believe that communication and engagement are only one part of the puzzle. There is a fundamental tension between science and policy because scientific endeavors are rightfully grounded in discovery, but policymakers formulate problems in terms of objectives, actions and outcomes. Decision science provides a solution by framing scientific questions in a way that is beneficial to policy development, facilitating scientists’ contribution to public discussion and policy. At its core, decision science is a field that aims to pinpoint evidence-based management strategies by focussing on those objectives, actions, and outcomes defined through the policy process. The importance of scientific discovery here is in linking actions to outcomes, helping decision-makers determine which actions best meet their objectives. In this paper we explain how problems can be formulated through the structured decision-making process. We give our vision for what decision science may grow to be, describing current gaps in methodology and application. By better understanding and engaging with the decision-making processes, scientists can have greater impact and make stronger contributions to important societal problems.

During the early stages of an emerging disease outbreak, governments are required to make critical decisions on how to respond, despite limited data being available to inform these decisions. Analytical risk assessment is a valuable approach to guide decision-making on travel restrictions and border measures during the early phase of an outbreak. Here we describe a rapid risk assessment framework that was developed in February 2020 to support time-critical decisions on the risk of SARS-CoV-2 importation into Australia. We briefly describe the context in which our framework was developed, the framework itself, and provide an example of the type of decision support provided to the Australian government. We then report a critical evaluation of the modelling choices made in February 2020, assessing the impact of our assumptions on estimated rates of importation, and provide a summary of "lessons learned". The framework presented and evaluated here provides a flexible approach to rapid assessment of importation risk, of relevance to current and future pandemic scenarios.

Malaria is caused by Plasmodium parasites which are transmitted to humans by the bite of an infected Anopheles mosquito. Plasmodium vivax is distinct from other malaria species in its ability to remain dormant in the liver (as hypnozoites) and activate later to cause further infections (referred to as relapses). Mathematical models to describe the transmission dynamics of P. vivax have been developed, but most of them fail to capture realistic dynamics of hypnozoites. Models that do capture the complexity tend to involve many governing equations, making them difficult to extend to incorporate other important factors for P. vivax, such as treatment status, age and pregnancy. In this paper, we have developed a multiscale model (a system of integro-differential equations) that involves a minimal set of equations at the population scale, with an embedded within-host model that can capture the dynamics of the hypnozoite reservoir. In this way, we can gain key insights into dynamics of P. vivax transmission with a minimum number of equations at the population scale, making this framework readily scalable to incorporate more complexity. We performed a sensitivity analysis of our multiscale model over key parameters and found that prevalence of P. vivax blood-stage infection increases with both bite rate and number of mosquitoes but decreases with hypnozoite death rate. Since our mathematical model captures the complex dynamics of P. vivax and the hypnozoite reservoir, it has the potential to become a key tool to inform elimination strategies for P. vivax.

Malaria remains a major public health burden in South-East Asia and Africa. Mathematical models of within-host infection dynamics and drug action, developed in support of malaria elimination initiatives, have significantly advanced our understanding of the dynamics of infection and supported development of effective drug-treatment regimens. However, the mathematical models supporting these initiatives are predominately based on deterministic dynamics and therefore cannot capture stochastic phenomena such as extinction (no parasitized red blood cells) following treatment, with potential consequences for our interpretation of data sets in which recrudescence is observed. Here we develop a stochastic within-host infection model to study the growth, decline and possible stochastic extinction of parasitized red blood cells in malaria-infected human volunteers. We show that stochastic extinction can occur when the inoculation size is small or when the number of parasitized red blood cells reduces significantly after an antimalarial treatment. We further show that the drug related parameters, such as the maximum killing rate and half-maximum effective concentration, are the primary factors determining the probability of stochastic extinction following treatment, highlighting the importance of highly-efficacious antimalarials in increasing the probability of cure for the treatment of malaria patients.

Malaria is a mosquito-borne disease with a devastating global impact. Plasmodium vivax is a major cause of human malaria beyond sub-Saharan Africa. Relapsing infections, driven by a reservoir of liver-stage parasites known as hypnozoites, present unique challenges for the control of P. vivax malaria. Following indeterminate dormancy periods, hypnozoites may activate to trigger relapses. Clearance of the hypnozoite reservoir through drug treatment (radical cure) has been proposed as a potential tool for the elimination of P. vivax malaria. Here, we introduce a stochastic, within-host model to jointly characterise hypnozoite and infection dynamics for an individual in a general transmission setting, allowing for radical cure. We begin by extending an existing activation-clearance model for a single hypnozoite, adapted to both short- and long-latency strains, to include drug treatment. We then embed this activation-clearance model in an epidemiological framework accounting for repeated mosquito inoculation and the administration of radical cure. By constructing an open network of infinite server queues, we derive analytic expressions for several quantities of epidemiological significance, including the size of the hypnozoite reservoir; the relapse rate; the relative contribution of relapses to the infection burden; the distribution of multiple infections; the cumulative number of recurrences over time, and the time to first recurrence following drug treatment. We derive from first principles the functional dependence between within-host and transmission parameters and patterns of blood- and liver-stage infection, whilst allowing for treatment under a mass drug administration regime. To yield population-level insights, our analytic within-host distributions can be embedded in multiscale models. Our work thus contributes to the epidemiological understanding of the effects of radical cure on P. vivax malaria.

National influenza pandemic plans have evolved substantially over recent decades, as has the scientific research that underpins the advice contained within them. While the knowledge generated by many research activities has been directly incorporated into the current generation of pandemic plans, scientists and policymakers are yet to capitalise fully on the potential for near real-time analytics to formally contribute to epidemic decision-making. Theoretical studies demonstrate that it is now possible to make robust estimates of pandemic impact in the earliest stages of a pandemic using first few hundred household cohort (FFX) studies and algorithms designed specifically for analysing FFX data. Pandemic plans already recognise the importance of both situational awareness i.e., knowing pandemic impact and its key drivers, and the need for pandemic special studies and related analytic methods for estimating these drivers. An important next step is considering how information from these situational assessment activities can be integrated into the decision-making processes articulated in pandemic planning documents. Here we introduce a decision support tool that directly uses outputs from FFX algorithms to present recommendations on response options, including a quantification of uncertainty, to decision makers. We illustrate this approach using response information from within the Australian influenza pandemic plan.

Most deaths from severe falciparum malaria occur within 24 h of presentation to a hospital. Intravenous (i.v.) artesunate is the first-line treatment for severe falciparum malaria, but its efficacy may be compromised by delayed parasitological responses. In patients with severe malaria, the life-saving benefit of the artemisinin derivatives is their ability to clear circulating parasites rapidly, before they can sequester and obstruct the microcirculation. To evaluate the dosing of i.v. artesunate for the treatment of artemisinin-sensitive and reduced ring stage sensitivity to artemisinin severe falciparum malaria infections, Bayesian pharmacokinetic-pharmacodynamic modeling of data from 94 patients with severe malaria (80 children from Africa and 14 adults from Southeast Asia) was performed. Assuming that delayed parasite clearance reflects a loss of ring stage sensitivity to artemisinin derivatives, the median (95% credible interval) percentage of patients clearing ≥99% of parasites within 24 h (PC24≥99%) for standard (2.4 mg/kg body weight i.v. artesunate at 0 and 12 h) and simplified (4 mg/kg i.v. artesunate at 0 h) regimens was 65% (52.5% to 74.5%) versus 44% (25% to 61.5%) for adults, 62% (51.5% to 74.5%) versus 39% (20.5% to 58.5%) for larger children (≥20 kg), and 60% (48.5% to 70%) versus 36% (20% to 53.5%) for smaller children (<20 kg). The upper limit of the credible intervals for all regimens was below a PC24≥99% of 80%, a threshold achieved on average in clinical studies of severe falciparum malaria infections. In severe falciparum malaria caused by parasites with reduced ring stage susceptibility to artemisinin, parasite clearance is predicted to be slower with both the currently recommended and proposed simplified i.v. artesunate dosing regimens.

MUC1 belongs to the family of cell surface (cs-) mucins. Experimental evidence indicates that its presence reduces in vivo influenza viral infection severity. However, the mechanisms by which MUC1 influences viral dynamics and the host immune response are not yet well understood, limiting our ability to predict the efficacy of potential treatments that target MUC1. To address this limitation, we use available in vivo kinetic data for both virus and macrophage populations in wildtype and MUC1 knockout mice. We apply two mathematical models of within-host influenza dynamics to this data. The models differ in how they categorise the mechanisms of viral control. Both models provide evidence that MUC1 reduces the susceptibility of epithelial cells to influenza virus and regulates macrophage recruitment. Furthermore, we predict and compare some key infection-related quantities between the two mice groups. We find that MUC1 significantly reduces the basic reproduction number of viral replication as well as the number of cumulative macrophages but has little impact on the cumulative viral load. Our analyses suggest that the viral replication rate in the early stages of infection influences the kinetics of the host immune response, with consequences for infection outcomes, such as severity. We also show that MUC1 plays a strong anti-inflammatory role in the regulation of the host immune response. This study improves our understanding of the dynamic role of MUC1 against influenza infection and may support the development of novel antiviral treatments and immunomodulators that target MUC1.

BACKGROUND: Allocation of scarce resources during a pandemic extends to the allocation of vaccines when they eventually become available. We describe a framework for priority vaccine allocation that employed a cross-disciplinary approach, guided by ethical considerations and informed by local risk assessment. METHODS: Published and grey literature was reviewed, and augmented by consultation with key informants, to collate past experience, existing guidelines and emerging strategies for pandemic vaccine deployment. Identified ethical issues and decision-making processes were also included. Concurrently, simulation modelling studies estimated the likely impacts of alternative vaccine allocation approaches. Assembled evidence was presented to a workshop of national experts in pandemic preparedness, vaccine strategy, implementation and ethics. All of this evidence was then used to generate a proposed ethical framework for vaccine priorities best suited to the Australian context. FINDINGS: Published and emerging guidance for priority pandemic vaccine distribution differed widely with respect to strategic objectives, specification of target groups, and explicit discussion of ethical considerations and decision-making processes. Flexibility in response was universally emphasised, informed by real-time assessment of the pandemic impact level, and identification of disproportionately affected groups. Model outputs aided identification of vaccine approaches most likely to achieve overarching goals in pandemics of varying transmissibility and severity. Pandemic response aims deemed most relevant for an Australian framework were: creating and maintaining trust, promoting equity, and reducing harmful outcomes. INTERPRETATION: Defining clear and ethically-defendable objectives for pandemic response in context aids development of flexible and adaptive decision support frameworks and facilitates clear communication and engagement activities.

Health planners from global to local levels must anticipate year-to-year and week-to-week variation in seasonal influenza activity when planning for and responding to epidemics to mitigate their impact. To help with this, countries routinely collect incidence of mild and severe respiratory illness and virologic data on circulating subtypes and use these data for situational awareness, burden of disease estimates and severity assessments. Advanced analytics and modelling are increasingly used to aid planning and response activities by describing key features of influenza activity for a given location and generating forecasts that can be translated to useful actions such as enhanced risk communications, and informing clinical supply chains. Here, we describe the formation of the Influenza Incidence Analytics Group (IIAG), a coordinated global effort to apply advanced analytics and modelling to public influenza data, both epidemiological and virologic, in real-time and thus provide additional insights to countries who provide routine surveillance data to WHO. Our objectives are to systematically increase the value of data to health planners by applying advanced analytics and forecasting and for results to be immediately reproducible and deployable using an open repository of data and code. We expect the resources we develop and the associated community to provide an attractive option for the open analysis of key epidemiological data during seasonal epidemics and the early stages of an influenza pandemic.

The ability of health systems to cope with coronavirus disease (COVID-19) cases is of major concern. In preparation, we used clinical pathway models to estimate healthcare requirements for COVID-19 patients in the context of broader public health measures in Australia. An age- and risk-stratified transmission model of COVID-19 demonstrated that an unmitigated epidemic would dramatically exceed the capacity of the health system of Australia over a prolonged period. Case isolation and contact quarantine alone are insufficient to constrain healthcare needs within feasible levels of expansion of health sector capacity. Overlaid social restrictions must be applied over the course of the epidemic to ensure systems do not become overwhelmed and essential health sector functions, including care of COVID-19 patients, can be maintained. Attention to the full pathway of clinical care is needed, along with ongoing strengthening of capacity.

Freya Shearer and co-authors discuss the use of decision analysis in planning for infectious disease pandemics.

Prevalence of impetigo (skin sores) remains high in remote Australian Aboriginal communities, Fiji, and other areas of socio-economic disadvantage. Skin sore infections, driven primarily in these settings by Group A Streptococcus (GAS) contribute substantially to the disease burden in these areas. Despite this, estimates for the force of infection, infectious period and basic reproductive ratio—all necessary for the construction of dynamic transmission models—have not been obtained. By utilising three datasets each containing longitudinal infection information on individuals, we estimate each of these epidemiologically important parameters. With an eye to future study design, we also quantify the optimal sampling intervals for obtaining information about these parameters. We verify the estimation method through a simulation estimation study, and test each dataset to ensure suitability to the estimation method. We find that the force of infection differs by population prevalence, and the infectious period is estimated to be between 12 and 20 days. We also find that optimal sampling interval depends on setting, with an optimal sampling interval between 9 and 11 days in a high prevalence setting, and 21 and 27 days for a lower prevalence setting. These estimates unlock future model-based investigations on the transmission dynamics of skin sores.

ABSTRACT

Background

Methods

Findings

Interpretation

Modern data and computational resources, coupled with algorithmic and theoretical advances to exploit these, allow disease dynamic models to be parameterised with increasing detail and accuracy. While this enhances models' usefulness in prediction and policy, major challenges remain. In particular, lack of identifiability of a model's parameters may limit the usefulness of the model. While lack of parameter identifiability may be resolved through incorporation into an inference procedure of prior knowledge, formulating such knowledge is often difficult. Furthermore, there are practical challenges associated with acquiring data of sufficient quantity and quality. Here, we discuss recent progress on these issues.

As of 1 May 2020, there had been 6808 confirmed cases of COVID-19 in Australia. Of these, 98 had died from the disease. The epidemic had been in decline since mid-March, with 308 cases confirmed nationally since 14 April. This suggests that the collective actions of the Australian public and government authorities in response to COVID-19 were sufficiently early and assiduous to avert a public health crisis - for now. Analysing factors that contribute to individual country experiences of COVID-19, such as the intensity and timing of public health interventions, will assist in the next stage of response planning globally. We describe how the epidemic and public health response unfolded in Australia up to 13 April. We estimate that the effective reproduction number was likely below one in each Australian state since mid-March and forecast that clinical demand would remain below capacity thresholds over the forecast period (from mid-to-late April).

Combinations of intense non-pharmaceutical interventions (lockdowns) were introduced worldwide to reduce SARS-CoV-2 transmission. Many governments have begun to implement exit strategies that relax restrictions while attempting to control the risk of a surge in cases. Mathematical modelling has played a central role in guiding interventions, but the challenge of designing optimal exit strategies in the face of ongoing transmission is unprecedented. Here, we report discussions from the Isaac Newton Institute 'Models for an exit strategy' workshop (11-15 May 2020). A diverse community of modellers who are providing evidence to governments worldwide were asked to identify the main questions that, if answered, would allow for more accurate predictions of the effects of different exit strategies. Based on these questions, we propose a roadmap to facilitate the development of reliable models to guide exit strategies. This roadmap requires a global collaborative effort from the scientific community and policymakers, and has three parts: (i) improve estimation of key epidemiological parameters; (ii) understand sources of heterogeneity in populations; and (iii) focus on requirements for data collection, particularly in low-to-middle-income countries. This will provide important information for planning exit strategies that balance socio-economic benefits with public health.

Bayesian methods have been used to predict the timing of infectious disease epidemics in various settings and for many infectious diseases, including seasonal influenza. But integrating these techniques into public health practice remains an ongoing challenge, and requires close collaboration between modellers, epidemiologists, and public health staff. During the 2016 and 2017 Australian influenza seasons, weekly seasonal influenza forecasts were produced for cities in the three states with the largest populations: Victoria, New South Wales and Queensland. Forecast results were presented to Health Department disease surveillance units in these jurisdictions, who provided feedback about the plausibility and public health utility of these predictions. In earlier studies we found that delays in reporting and processing of surveillance data substantially limited forecast performance, and that incorporating climatic effects on transmission improved forecast performance. In this study of the 2016 and 2017 seasons, we sought to refine the forecasting method to account for delays in receiving the data, and used meteorological data from past years to modulate the force of infection. We demonstrate how these refinements improved the forecast’s predictive capacity, and use the 2017 influenza season to highlight challenges in accounting for population and clinician behaviour changes in response to a severe season.

Renewed efforts to eliminate malaria have highlighted the potential to interrupt human-to-mosquito transmission - a process mediated by gametocyte kinetics in human hosts. Here we study the in vivo dynamics of Plasmodium falciparum gametocytes by establishing a framework which incorporates improved measurements of parasitemia, a novel gametocyte dynamics model and model fitting using Bayesian hierarchical inference. We found that the model provides an excellent fit to the clinical data from 17 volunteers infected with P. falciparum (3D7 strain) and reliably predicts observed gametocytemia. We estimated the sexual commitment rate and gametocyte sequestration time to be 0.54% (95% credible interval: 0.30-1.00%) per asexual replication cycle and 8.39 (6.54-10.59) days respectively. We used the data-calibrated model to investigate human-to-mosquito transmissibility, providing a method to link within-human host infection kinetics to epidemiological-scale infection and transmission patterns.

Laboratory models are often used to understand the interaction of related pathogens via host immunity. For example, recent experiments where ferrets were exposed to two influenza strains within a short period of time have shown how the effects of cross-immunity vary with the time between exposures and the specific strains used. On the other hand, studies of the workings of different arms of the immune response, and their relative importance, typically use experiments involving a single infection. However, inferring the relative importance of different immune components from this type of data is challenging. Using simulations and mathematical modelling, here we investigate whether the sequential infection experiment design can be used not only to determine immune components contributing to cross-protection, but also to gain insight into the immune response during a single infection. We show that virological data from sequential infection experiments can be used to accurately extract the timing and extent of cross-protection. Moreover, the broad immune components responsible for such cross-protection can be determined. Such data can also be used to infer the timing and strength of some immune components in controlling a primary infection, even in the absence of serological data. By contrast, single infection data cannot be used to reliably recover this information. Hence, sequential infection data enhances our understanding of the mechanisms underlying the control and resolution of infection, and generates new insight into how previous exposure influences the time course of a subsequent infection.

For diseases such as influenza, where the majority of infected persons experience mild (if any) symptoms, surveillance systems are sensitive to changes in healthcare-seeking and clinical decision-making behaviours. This presents a challenge when trying to interpret surveillance data in near-real-time (e.g., to provide public health decision-support). Australia experienced a particularly large and severe influenza season in 2017, perhaps in part due to: (a) mild cases being more likely to seek healthcare; and (b) clinicians being more likely to collect specimens for reverse transcription polymerase chain reaction (RT-PCR) influenza tests. In this study, we used weekly Flutracking surveillance data to estimate the probability that a person with influenza-like illness (ILI) would seek healthcare and have a specimen collected. We then used this estimated probability to calibrate near-real-time seasonal influenza forecasts at each week of the 2017 season, to see whether predictive skill could be improved. While the number of self-reported influenza tests in the weekly surveys are typically very low, we were able to detect a substantial change in healthcare seeking behaviour and clinician testing behaviour prior to the high epidemic peak. Adjusting for these changes in behaviour in the forecasting framework improved predictive skill. Our analysis demonstrates a unique value of community-level surveillance systems, such as Flutracking, when interpreting traditional surveillance data. These methods are also applicable beyond the Australian context, as similar community-level surveillance systems operate in other countries.

Antimalarial treatment currently relies on an artemisinin derivative and a longer-acting partner drug. With the emergence of resistance to the artemisinin derivatives and the potential pressure this exerts on the partner drugs, the impact of resistance to each drug on efficacy needs to be investigated. An in silico exploration of dihydroartemisinin-piperaquine and mefloquine-artesunate, two artemisinin-based combination therapies that are commonly used in Southeast Asia, was performed. The percentage of treatment failures was simulated from a within-host pharmacokinetic-pharmacodynamic (PKPD) model, assuming that parasites developed increasing levels of (i) artemisinin derivative resistance or (ii) concomitant resistance to both the artemisinin derivative and the partner drug. Because the exact nature of how resistant Plasmodium falciparum parasites respond to treatment is unknown, we examined the impact on treatment failure rates of artemisinin resistance that (i) reduced the maximal killing rate, (ii) increased the concentration of drug required for 50% killing, or (iii) shortened the window of parasite stages that were susceptible to artemisinin derivatives until the drugs had no effect on the ring stages. The loss of the ring-stage activity of the artemisinin derivative caused the greatest increase in the treatment failure rate, and this result held irrespective of whether partner drug resistance was assumed to be present or not. To capture the uncertainty regarding how artemisinin derivative and partner drug resistance affects the assumed concentration-killing effect relationship, a variety of changes to this relationship should be considered when using within-host PKPD models to simulate clinical outcomes to guide treatment strategies for resistant infections.

Malaria infection continues to be a major health problem worldwide and drug resistance in the major human parasite species, Plasmodium falciparum, is increasing in South East Asia. Control measures including novel drugs and vaccines are in development, and contributions to the rational design and optimal usage of these interventions are urgently needed. Infection involves the complex interaction of parasite dynamics, host immunity, and drug effects. The long life cycle (48 hours in the common human species) and synchronized replication cycle of the parasite population present significant challenges to modeling the dynamics of Plasmodium infection. Coupled with these, variation in immune recognition and drug action at different life cycle stages leads to further complexity. We review the development and progress of “within-host” models of Plasmodium infection, and how these have been applied to understanding and interpreting human infection and animal models of infection.

© 2018 Elsevier Inc. Infections with Sarcoptes scabiei, or scabies, remain common in many disadvantaged populations. Mass drug administration (MDA) has been used in such settings to achieve a rapid reduction in infection and transmission, with the goal of eliminating the public health burden of scabies. While prevalence has been observed to fall substantially following such an intervention, in some instances resurgence of infection to baseline levels has occurred over several years. To explore the biology underpinning this phenomenon, we have developed a theoretical model of scabies life-cycle and transmission dynamics in a homogeneously mixing population, and simulate the impact of mass drug treatment strategies acting on egg and mite life cycle stages (ovicidal) or mites alone (non-ovicidal). In order to investigate the dynamics of the system, we first define and calculate the optimal interval between treatment doses. We calculate the probability of eradication as a function of the number of optimally-timed successive treatment doses and the number of years over which a program is run. For the non-ovicidal intervention, we first show that at least two optimally-timed doses are required to achieve eradication. We then demonstrate that while more doses over a small number of years provides the highest chance of eradication, a similar outcome can be achieved with fewer doses delivered annually over a longer period of time. For the ovicidal intervention, we find that doses should be delivered as close together as possible. This work provides a platform for further research into optimal treatment strategies which may incorporate heterogeneity of transmission, and the interplay between MDA and enhancement of continuing scabies surveillance and treatment strategies.

In some disease systems, the process of waning immunity can be subtle, involving a complex relationship between the duration of immunity-acquired either through natural infection or vaccination-and subsequent boosting of immunity through asymptomatic re-exposure. We present and analyse a model of infectious disease transmission where primary and secondary infections are distinguished to examine the interplay between infection and immunity. Additionally we allow the duration of infection-acquired immunity to differ from that of vaccine-acquired immunity to explore the impact on long-term disease patterns and prevalence of infection in the presence of immune boosting. Our model demonstrates that vaccination may induce cyclic behaviour, and the ability of vaccinations to reduce primary infections may not lead to decreased transmission. Where the boosting of vaccine-acquired immunity delays a primary infection, the driver of transmission largely remains primary infections. In contrast, if the immune boosting bypasses a primary infection, secondary infections become the main driver of transmission under a sufficiently long duration of immunity. Our results show that the epidemiological patterns of an infectious disease may change considerably when the duration of vaccine-acquired immunity differs from that of infection-acquired immunity. Our study highlights that for any particular disease and associated vaccine, a detailed understanding of the waning and boosting of immunity and how the duration of protection is influenced by infection prevalence are important as we seek to optimise vaccination strategies.

The first line treatment for uncomplicated falciparum malaria is artemisinin-based combination therapy (ACT), which consists of an artemisinin derivative coadministered with a longer-acting partner drug. However, the spread of Plasmodium falciparum resistant to both artemisinin and its partner drugs poses a major global threat to malaria control activities. Novel strategies are needed to retard and reverse the spread of these resistant parasites. One such strategy is triple artemisinin-based combination therapy (TACT). We developed a mechanistic within-host mathematical model to investigate the efficacy of a TACT (dihydroartemisinin-piperaquine-mefloquine [DHA-PPQ-MQ]) for use in South-East Asia, where DHA and PPQ resistance are now increasingly prevalent. Comprehensive model simulations were used to explore the degree to which the underlying resistance influences the parasitological outcomes. The effect of MQ dosing on the efficacy of TACT was quantified at various degrees of DHA and PPQ resistance. To incorporate interactions between drugs, a novel model is presented for the combined effect of DHA-PPQ-MQ, which illustrates how the interactions can influence treatment efficacy. When combined with a standard regimen of DHA and PPQ, the administration of three 6.7-mg/kg doses of MQ was sufficient to achieve parasitological efficacy greater than that currently recommended by World Health Organization (WHO) guidelines. As a result, three 8.3-mg/kg doses of MQ, the current WHO-recommended dosing regimen for MQ, combined with DHA-PPQ, has the potential to produce high cure rates in regions where resistance to DHA-PPQ has emerged.

OBJECTIVE: Recent studies have used Bayesian methods to predict timing of influenza epidemics many weeks in advance, but there is no documented evaluation of how such forecasts might support the day-to-day operations of public health staff. METHODS: During the 2015 influenza season in Melbourne, Australia, weekly forecasts were presented at Health Department surveillance unit meetings, where they were evaluated and updated in light of expert opinion to improve their accuracy and usefulness. RESULTS: Predictive capacity of the model was substantially limited by delays in reporting and processing arising from an unprecedented number of notifications, disproportionate to seasonal intensity. Adjustment of the predictive algorithm to account for these delays and increased reporting propensity improved both current situational awareness and forecasting accuracy. CONCLUSIONS: Collaborative engagement with public health practitioners in model development improved understanding of the context and limitations of emerging surveillance data. Incorporation of these insights in a quantitative model resulted in more robust estimates of disease activity for public health use. Implications for public health: In addition to predicting future disease trends, forecasting methods can quantify the impact of delays in data availability and variable reporting practice on the accuracy of current epidemic assessment. Such evidence supports investment in systems capacity.

Treatment options for influenza B virus infections are limited to neuraminidase inhibitors (NAIs), which block the neuraminidase (NA) glycoprotein on the virion surface. The development of NAI resistance would therefore result in a loss of antiviral treatment options for influenza B virus infections. This study characterized two contemporary influenza B viruses with known resistance-conferring NA amino acid substitutions, D197N and H273Y, detected during routine surveillance. The D197N and H273Y variants were characterized in vitro by assessing NA enzyme activity and affinity, as well as replication in cell culture compared to those of NAI-sensitive wild-type viruses. In vivo studies were also performed in ferrets to assess the replication and transmissibility of each variant. Mathematical models were used to analyze within-host and between-host fitness of variants relative to wild-type viruses. The data revealed that the H273Y variant had NA enzyme function similar to that of its wild type but had slightly reduced replication and transmission efficiency in vivo The D197N variant had impaired NA enzyme function, but there was no evidence of reduction in replication or transmission efficiency in ferrets. Our data suggest that the influenza B virus variant with the H273Y NA substitution had a more notable reduction in fitness compared to wild-type viruses than the influenza B variant with the D197N NA substitution. Although a D197N variant is yet to become widespread, it is the most commonly detected NAI-resistant influenza B virus in surveillance studies. Our results highlight the need to carefully monitor circulating viruses for the spread of influenza B viruses with the D197N NA substitution.

Falciparum malaria is a major parasitic disease causing widespread morbidity and mortality globally. Artemisinin derivatives—the most effective and widely-used antimalarials that have helped reduce the burden of malaria by 60% in some areas over the past decade—have recently been found to induce growth retardation of blood-stage Plasmodium falciparum when applied at clinically relevant concentrations. To date, no model has been designed to quantify the growth retardation effect and to predict the influence of this property on in vivo parasite killing. Here we introduce a mechanistic model of parasite growth from the ring to trophozoite stage of the parasite's life cycle, and by modelling the level of staining with an RNA-binding dye, we demonstrate that the model is able to reproduce fluorescence distribution data from in vitro experiments using the laboratory 3D7 strain. We quantify the dependence of growth retardation on drug concentration and identify the concentration threshold above which growth retardation is evident. We estimate that the parasite life cycle is prolonged by up to 10 hours. We illustrate that even such a relatively short delay in growth may significantly influence in vivo parasite dynamics, demonstrating the importance of considering growth retardation in the design of optimal artemisinin-based dosing regimens.