A hybrid transmission model for Plasmodium vivax accounting for superinfection, immunity and the hypnozoite reservoir

Malaria is a vector-borne disease that exacts a grave toll in the Global South. The epidemiology of Plasmodium vivax, the most geographically expansive agent of human malaria, is characterised by the accrual of a reservoir of dormant parasites known as hypnozoites. Relapses, arising from hypnozoite activation events, comprise the majority of the blood-stage infection burden, with implications for the acquisition of immunity and the distribution of superinfection. Here, we construct a novel model for the transmission of P. vivax that concurrently accounts for the accrual of the hypnozoite reservoir, (blood-stage) superinfection and the acquisition of immunity. We begin by using an infinite-server queueing network model to characterise the within-host dynamics as a function of mosquito-to-human transmission intensity, extending our previous model to capture a discretised immunity level. To model transmission-blocking and antidisease immunity, we allow for geometric decay in the respective probabilities of successful human-to-mosquito transmission and symptomatic blood-stage infection as a function of this immunity level. Under a hybrid approximation-whereby probabilistic within-host distributions are cast as expected population-level proportions-we couple host and vector dynamics to recover a deterministic compartmental model in line with Ross-Macdonald theory. We then perform a steady-state analysis for this compartmental model, informed by the (analytic) distributions derived at the within-host level. To characterise transient dynamics, we derive a reduced system of integrodifferential equations, likewise informed by our within-host queueing network, allowing us to recover population-level distributions for various quantities of epidemiological interest. In capturing the interplay between hypnozoite accrual, superinfection and acquired immunity-and providing, to the best of our knowledge, the most complete population-level distributions for a range of epidemiological values-our model provides insights into important, but poorly understood, epidemiological features of P. vivax.

How immune dynamics shape multi-season epidemics: a continuous-discrete model in one dimensional antigenic space

We extend a previously published model for the dynamics of a single strain of an influenza-like infection. The model incorporates a waning acquired immunity to infection and punctuated antigenic drift of the virus, employing a set of coupled integral equations within a season and a discrete map between seasons. The long term behaviour of the model is demonstrated by examples where immunity to infection depends on the time since a host was last infected, and where immunity depends on the number of times that a host has been infected. The first scenario leads to complicated dynamics in some regions of parameter space, and to regions of parameter space with more than one attractor. The second scenario leads to a stable fixed point, corresponding to an identical epidemic each season. We also examine the model with both paradigms in combination, almost always but not exclusively observing a stable fixed point or periodic solution. Adding stochastic perturbations to the between season map fails to destroy the model's qualitative dynamics. Our results suggest that if the level of host immunity depends on the elapsed time since the last infection then the epidemiological dynamics may be unpredictable.

Characterisation of Plasmodium vivax lactate dehydrogenase dynamics in P. vivax infections

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.

Estimating measures to reduce the transmission of SARS-CoV-2 in Australia to guide a 'National Plan' to reopening

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.

Mathematical models of Plasmodium vivax transmission: A scoping review

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 focus on refining how key aspects of P. vivax dynamics are modelled, including spatial heterogeneity in exposure risk and heterogeneity in susceptibility to infection, the accumulation of hypnozoite variation, the interaction between P. falciparum and P. vivax, acquisition of immunity, and recovery under superinfection.

Key Challenges for Respiratory Virus Surveillance while Transitioning out of Acute Phase of COVID-19 Pandemic

To support the ongoing management of viral respiratory diseases while transitioning out of the acute phase of the COVID-19 pandemic, many countries are moving toward an integrated model of surveillance for SARS-CoV-2, influenza virus, and other respiratory pathogens. Although many surveillance approaches catalyzed by the COVID-19 pandemic provide novel epidemiologic insight, continuing them as implemented during the pandemic is unlikely to be feasible for nonemergency surveillance, and many have already been scaled back. Furthermore, given anticipated cocirculation of SARS-CoV-2 and influenza virus, surveillance activities in place before the pandemic require review and adjustment to ensure their ongoing value for public health. In this report, we highlight key challenges for the development of integrated models of surveillance. We discuss the relative strengths and limitations of different surveillance practices and studies as well as their contribution to epidemiologic assessment, forecasting, and public health decision-making.

Improving estimates of waning immunity rates in stochastic SIRS models with a hierarchical framework

As most disease causing pathogens require transmission from an infectious individual to a susceptible individual, continued persistence of the pathogen within the population requires the replenishment of susceptibles through births, immigration, or waning immunity. Consider the introduction of an unknown infectious disease into a fully susceptible population where it is not known how long immunity is conferred once an individual recovers from infection. If, initially, the prevalence of disease increases (that is, the infection takes off), the number of infectives will usually decrease to a low level after the first major outbreak. During this post-outbreak period, the disease dynamics may be influenced by stochastic effects and there is a non-zero probability that the epidemic will die out. Die out in this period following the first major outbreak is known as an epidemic fade-out. If the disease does not die out, the susceptible population may be replenished by the waning of immunity, and a second wave may start. In this study, we investigate if the rate of waning immunity (and other epidemiological parameters) can be reliably estimated from multiple outbreak data, in which some outbreaks display epidemic fade-out and others do not. We generated synthetic outbreak data from independent simulations of stochastic SIRS models in multiple communities. Some outbreaks faded-out and some did not. We conducted Bayesian parameter estimation under two alternative approaches: independently on each outbreak and under a hierarchical framework. When conducting independent estimation, the waning immunity rate was poorly estimated and biased towards zero when an epidemic fade-out was observed. However, under a hierarchical approach, we obtained more accurate and precise posterior estimates for the rate of waning immunity and other epidemiological parameters. The greatest improvement in estimates was obtained for those communities in which epidemic fade-out was observed. Our findings demonstrate the feasibility and value of adopting a Bayesian hierarchical approach for parameter inference for stochastic epidemic models.

Superinfection and the hypnozoite reservoir for Plasmodium vivax: a general framework

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 malarial superinfection that can be generalised to allow for greater biological realism in a number of directions.

Seasonality as a driver of pH1N12009 influenza vaccination campaign impact

Although the most recent respiratory virus pandemic was triggered by a Coronavirus, sustained and elevated prevalence of highly pathogenic avian influenza viruses able to infect mammalian hosts highlight the continued threat of pandemics of influenza A virus (IAV) to global health. Retrospective analysis of pandemic outcomes, including comparative investigation of intervention efficacy in different regions, provide important contributions to the evidence base for future pandemic planning. The swine-origin IAV pandemic of 2009 exhibited regional variation in onset, infection dynamics and annual infection attack rates (IARs). For example, the UK experienced three severe peaks of infection over two influenza seasons, whilst Australia experienced a single severe wave. We adopt a seasonally forced 2-subtype model for the transmission of pH1N12009 and seasonal H3N2 to examine the role vaccination campaigns may play in explaining differences in pandemic trajectories in temperate regions. Our model differentiates between the nature of vaccine- and infection-acquired immunity. In particular, we assume that immunity triggered by infection elicits heterologous cross-protection against viral shedding in addition to long-lasting neutralising antibody, whereas vaccination induces imperfect reduction in susceptibility. We employ an Approximate Bayesian Computation (ABC) framework to calibrate the model using data for pH1N12009 seroprevalence, relative subtype dominance, and annual IARs for Australia and the UK. Heterologous cross-protection substantially suppressed the pandemic IAR over the posterior, with the strength of protection against onward transmission inversely correlated with the initial reproduction number. We show that IAV pandemic timing relative to the usual seasonal influenza cycle influenced the size of the initial waves of pH1N12009 in temperate regions and the impact of vaccination campaigns.

How effective were Australian Quarantine Stations in mitigating transmission aboard ships during the influenza pandemic of 1918-19?

The influenza pandemic of 1918-19 was the most devastating pandemic of the 20th century. It killed an estimated 50-100 million people worldwide. In late 1918, when the severity of the disease was apparent, the Australian Quarantine Service was established. Vessels returning from overseas and inter-state were intercepted, and people were examined for signs of illness and quarantined. Some of these vessels carried the infection throughout their voyage and cases were prevalent by the time the ship arrived at a Quarantine Station. We study four outbreaks that took place on board the Medic, Boonah, Devon, and Manuka in late 1918. These ships had returned from overseas and some of them were carrying troops that served in the First World War. By analysing these outbreaks under a stochastic Bayesian hierarchical modeling framework, we estimate the transmission rates among crew and passengers aboard these ships. Furthermore, we ask whether the removal of infectious, convalescent, and healthy individuals after arriving at a Quarantine Station in Australia was an effective public health response.

Individual variation in vaccine immune response can produce bimodal distributions of protection

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.

Choice of spatial discretisation influences the progression of viral infection within multicellular tissues

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.

A model for malaria treatment evaluation in the presence of multiple species

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.

COVID-19 vaccine coverage targets to inform reopening plans in a low incidence setting

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.

Understanding the impact of disease and vaccine mechanisms on the importance of optimal vaccine allocation

Vaccination is an important epidemic intervention strategy. However, it is generally unclear how the outcomes of different vaccine strategies change depending on population characteristics, vaccine mechanisms and allocation objective. In this paper we develop a conceptual mathematical model to simulate strategies for pre-epidemic vaccination. We extend the SEIR model to incorporate a range of vaccine mechanisms and disease characteristics. We then compare the outcomes of optimal and suboptimal vaccination strategies for three public health objectives (total infections, total symptomatic infections and total deaths) using numerical optimisation. Our comparison shows that the difference in outcomes between vaccinating optimally and suboptimally depends on vaccine mechanisms, disease characteristics, and objective considered. Our modelling finds vaccines that impact transmission produce better outcomes as transmission is reduced for all strategies. For vaccines that impact the likelihood of symptomatic disease or dying due to infection, the improvement in outcome as we decrease these variables is dependent on the strategy implemented. Through a principled model-based process, this work highlights the importance of designing effective vaccine allocation strategies. We conclude that efficient allocation of resources can be just as crucial to the success of a vaccination strategy as the vaccine effectiveness and/or amount of vaccines available.

Forecasting COVID-19 activity in Australia to support pandemic response: May to October 2020

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.

Optimal Interruption of P. vivax Malaria Transmission Using Mass Drug Administration

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 elimination under our mathematical model (and choice of model parameters) since the prevalence of infection eventually returns to pre-MDA levels.

Enhanced viral infectivity and reduced interferon production are associated with high pathogenicity for influenza viruses

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.

A modelling approach to estimate the transmissibility of SARS-CoV-2 during periods of high, low, and zero case incidence

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.

Real-time analysis of hospital length of stay in a mixed SARS-CoV-2 Omicron and Delta epidemic in New South Wales, Australia

BACKGROUND

The distribution of the duration that clinical cases of COVID-19 occupy hospital beds (the 'length of stay') is a key factor in determining how incident caseloads translate into health system burden. Robust estimation of length of stay in real-time requires the use of survival methods that can account for right-censoring induced by yet unobserved events in patient progression (e.g. discharge, death). In this study, we estimate in real-time the length of stay distributions of hospitalised COVID-19 cases in New South Wales, Australia, comparing estimates between a period where Delta was the dominant variant and a subsequent period where Omicron was dominant.

METHODS

Using data on the hospital stays of 19,574 individuals who tested positive to COVID-19 prior to admission, we performed a competing-risk survival analysis of COVID-19 clinical progression.

RESULTS

During the mixed Omicron-Delta epidemic, we found that the mean length of stay for individuals who were discharged directly from ward without an ICU stay was, for age groups 0-39, 40-69 and 70 +, respectively, 2.16 (95% CI: 2.12-2.21), 3.93 (95% CI: 3.78-4.07) and 7.61 days (95% CI: 7.31-8.01), compared to 3.60 (95% CI: 3.48-3.81), 5.78 (95% CI: 5.59-5.99) and 12.31 days (95% CI: 11.75-12.95) across the preceding Delta epidemic (1 July 2021-15 December 2021). We also considered data on the stays of individuals within the Hunter New England Local Health District, where it was reported that Omicron was the only circulating variant, and found mean ward-to-discharge length of stays of 2.05 (95% CI: 1.80-2.30), 2.92 (95% CI: 2.50-3.67) and 6.02 days (95% CI: 4.91-7.01) for the same age groups.

CONCLUSIONS

Hospital length of stay was substantially reduced across all clinical pathways during a mixed Omicron-Delta epidemic compared to a prior Delta epidemic, contributing to a lessened health system burden despite a greatly increased infection burden. Our results demonstrate the utility of survival analysis in producing real-time estimates of hospital length of stay for assisting in situational assessment and planning of the COVID-19 response.

A Multiscale Mathematical Model of Plasmodium Vivax Transmission

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.

From Climate Change to Pandemics: Decision Science Can Help Scientists Have Impact

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.

Rapid assessment of the risk of SARS-CoV-2 importation to Australia in February 2020

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.

Development of an influenza pandemic decision support tool linking situational analytics to national response policy

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.

Development and Validation of an In Silico Decision Tool To Guide Optimization of Intravenous Artesunate Dosing Regimens for Severe Falciparum Malaria Patients

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.

Modelling the Effect of MUC1 on Influenza Virus Infection Kinetics and Macrophage Dynamics

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.

Estimation of the force of infection and infectious period of skin sores in remote Australian communities using interval-censored data

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.

Coronavirus Disease Model to Inform Transmission-Reducing Measures and Health System Preparedness, Australia

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.

Early analysis of the Australian COVID-19 epidemic

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).

Key questions for modelling COVID-19 exit strategies

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.

Influencing public health policy with data-informed mathematical models of infectious diseases: Recent developments and new challenges

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.

Coordinating the real-time use of global influenza activity data for better public health planning

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.

Infectious disease pandemic planning and response: Incorporating decision analysis

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

Modeling the dynamics of Plasmodium falciparum gametocytes in humans during malaria infection

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.

Within-host modeling of blood-stage malaria

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.

Anatomy of a seasonal influenza epidemic forecast

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.

Accounting for Healthcare-Seeking Behaviours and Testing Practices in Real-Time Influenza Forecasts

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.

Sequential infection experiments for quantifying innate and adaptive immunity during influenza infection

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.

Investigating the Efficacy of Triple Artemisinin-Based Combination Therapies for Treating Plasmodium falciparum Malaria Patients Using Mathematical Modeling

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.

Characterization of Influenza B Virus Variants with Reduced Neuraminidase Inhibitor Susceptibility

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.

A biological model of scabies infection dynamics and treatment informs mass drug administration strategies to increase the likelihood of elimination

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.

Infection-acquired versus vaccine-acquired immunity in an SIRWS model

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.

Clonally diverse CD38+HLA-DR+CD8+ T cells persist during fatal H7N9 disease

Severe influenza A virus (IAV) infection is associated with immune dysfunction. Here, we show circulating CD8+ T-cell profiles from patients hospitalized with avian H7N9, seasonal IAV, and influenza vaccinees. Patient survival reflects an early, transient prevalence of highly activated CD38+HLA-DR+PD-1+ CD8+ T cells, whereas the prolonged persistence of this set is found in ultimately fatal cases. Single-cell T cell receptor (TCR)-αβ analyses of activated CD38+HLA-DR+CD8+ T cells show similar TCRαβ diversity but differential clonal expansion kinetics in surviving and fatal H7N9 patients. Delayed clonal expansion associated with an early dichotomy at a transcriptome level (as detected by single-cell RNAseq) is found in CD38+HLA-DR+CD8+ T cells from patients who succumbed to the disease, suggesting a divergent differentiation pathway of CD38+HLA-DR+CD8+ T cells from the outset during fatal disease. Our study proposes that effective expansion of cross-reactive influenza-specific TCRαβ clonotypes with appropriate transcriptome signatures is needed for early protection against severe influenza disease.

Evidence for Viral Interference and Cross-reactive Protective Immunity Between Influenza B Virus Lineages

Background

Two influenza B virus lineages, B/Victoria and B/Yamagata, cocirculate in the human population. While the lineages are serologically distinct, cross-reactive responses to both lineages have been detected. Viral interference describes the situation whereby infection with one virus limits infection and replication of a second virus. We investigated the potential for viral interference between the influenza B virus lineages.

Methods

Ferrets were infected and then challenged 3, 10, or 28 days later with pairs of influenza B/Victoria and B/Yamagata viruses.

Results

Viral interference occurred at challenge intervals of 3 and 10 days and occasionally at 28 days. At the longer interval, shedding of challenge virus was reduced, and this correlated with cross-reactive interferon γ responses from lymph nodes from virus-infected animals. Viruses from both lineages could prevent or significantly limit subsequent infection with a virus from the other lineage. Coinfections were rare, indicating the potential for reassortment between lineages is limited.

Conclusions

These data suggest that innate and cross-reactive immunity mediate viral interference and that this may contribute to the dominance of a specific influenza B virus lineage in any given influenza season. Furthermore, infection with one influenza B virus lineage may be beneficial in protecting against subsequent infection with either influenza B virus lineage.

Epidemic forecasts as a tool for public health: interpretation and (re)calibration

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.

A mechanistic model quantifies artemisinin-induced parasite growth retardation in blood-stage Plasmodium falciparum infection

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.

Influenza as a trigger for cardiovascular disease: An investigation of serotype, subtype and geographic location

BACKGROUND

Seasonal peaks of influenza and cardiovascular disease tend to coincide. Many excess deaths may be triggered by influenza, and the severity of this effect may vary with the virulence of the circulating influenza strain and host susceptibility. We aimed to explore the association between hospital admissions for influenza and/or pneumonia (IP) and acute myocardial infarction (AMI) or ischaemic heart disease (IHD) in Queensland, Australia, taking into account temporal and spatial variation of influenza virus type and subtype in 2007, 2008 and 2009.

METHODS

This ecological study at Statistical Subdivision level (SSD, n=38) used linked patient-level data. For each study year, Standardized Morbidity Ratios (SMRs) were calculated for hospital admissions with diagnoses of IP, AMI and IHD. We investigated the associations between IP and AMI or IHD using spatial autoregressive modelling, adjusting for socio-demographic factors.

RESULTS

Spatial autocorrelation was detected in SMRs, possibly reflecting underlying social and behavioural risk factors, but consistent with infectious disease spread. SMRs for IP were consistently predictive of SMRs for AMI and IHD when adjusted for socioeconomic status, population density and per cent Indigenous population (coefficient: 0.707, 95% confidence interval (CI): 0.318 - 1.096; 0.553, 0.222 - 0.884; 0.598, 0.307 - 0.888 and 1.017, 0.711 - 1.323; 0.650, 0.342 - 0.958; 1.031, 0.827 - 1.236) in 2007, 2008 and 2009, respectively.

CONCLUSIONS

This ecological study provides further evidence that severe respiratory infections may trigger the onset of cardiovascular events, implicating the influenza virus as a contributing factor.