Opportunities to strengthen respiratory virus surveillance systems in Australia: lessons learned from the COVID-19 response.

Abstract: 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.

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.

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.

Estimating the impact of test-trace-isolate-quarantine systems on SARS-CoV-2 transmission in Australia.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.