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

Simulation-based inference of the time-dependent reproduction number from temporally aggregated and under-reported disease incidence time series data

During infectious disease outbreaks, the time-dependent reproduction number ([Formula: see text]) can be estimated to monitor pathogen transmission. In previous work, we developed a simulation-based method for estimating [Formula: see text] from temporally aggregated disease incidence data (e.g. weekly case reports). While that approach is straightforward to use, it assumes implicitly that all cases are reported and the computation can be slow when applied to large datasets. In this article, we extend our previous approach and develop a computationally efficient simulation-based method for estimating [Formula: see text] in real-time accounting for both temporal aggregation of incidence data and under-reporting (with a fixed reporting probability per case). Using simulated data, we show that failing to consider stochastic under-reporting can lead to inappropriately precise estimates, including scenarios in which the true [Formula: see text] value lies outside inferred credible intervals more often than expected. We then apply our approach to data from the 2018 to 2020 Ebola outbreak in the Democratic Republic of the Congo (DRC), again exploring the effects of case under-reporting. Finally, we show how our method can be extended to account for temporal variations in reporting. Given information about the level of case reporting, our framework can be used to estimate [Formula: see text] during future outbreaks with under-reported and temporally aggregated case data.This article is part of the theme issue 'Uncertainty quantification for healthcare and biological systems (Part 2)'.

Quantifying age-specific household contacts in Aotearoa New Zealand for infectious disease modelling

Accounting for population age structure and age-specific contact patterns is crucial for accurate modelling of human infectious disease dynamics and impact. A common approach is to use contact matrices, which estimate the number of contacts between individuals of different ages. These contact matrices are frequently based on data collected from populations with very different demographic and socio-economic characteristics from the population of interest. Here we use a comprehensive household composition dataset based on Aotearoa New Zealand census and administrative data to construct a household contact matrix and a synthetic population that can be used for modelling. We investigate the behaviour of a compartment-based and an agent-based epidemic model parametrized using these data, compared with a commonly used contact matrix that was constructed by projecting international data onto New Zealand's population. We find that using the New Zealand household data, either in a compartment-based model or in an agent-based model, leads to lower attack rates in older age groups compared with using the projected contact matrix. This difference becomes larger when household transmission is more dominant relative to non-household transmission. We provide electronic versions of the synthetic population and household contact matrix for other researchers to use in infectious disease models.

Structured methods for parameter inference and uncertainty quantification for mechanistic models in the life sciences

Parameter inference and uncertainty quantification are important steps when relating mathematical models to real-world observations and when estimating uncertainty in model predictions. However, methods for doing this can be computationally expensive, particularly when the number of unknown model parameters is large. The aim of this study is to develop and test an efficient profile likelihood-based method, which takes advantage of the structure of the mathematical model being used. We do this by identifying specific parameters that affect model output in a known way, such as a linear scaling. We illustrate the method by applying it to three toy models from different areas of the life sciences: (i) a predator-prey model from ecology; (ii) a compartment-based epidemic model from health sciences; and (iii) an advection-diffusion reaction model describing the transport of dissolved solutes from environmental science. We show that the new method produces results of comparable accuracy to existing profile likelihood methods but with substantially fewer evaluations of the forward model. We conclude that our method could provide a much more efficient approach to parameter inference for models where a structured approach is feasible. Computer code to apply the new method to user-supplied models and data is provided via a publicly accessible repository.

The timing of local SARS-Cov-2 outbreaks and vaccination coverage during the Delta wave in Melbourne

OBJECTIVE

This article presents a longitudinal analysis of COVID-19 infection and vaccination coverage in Melbourne metropolitan local government areas (LGAs) during the 2021 Delta wave.

METHODS

COVID-19 vaccination and infection data from 12 July to 27 November 2021 were sourced from government websites. Summary statistics and associated 95% confidence intervals (95% CI) were compared by LGA ranked according to socioeconomic status: total "burden" (total infections per thousand), "peak" (highest weekly infection rate), "lag" (interval between peak and 70% double vaccination).

RESULTS

LGAs in the bottom five deciles for social advantage experienced higher infection rates (39.0 per thousand [95% CI: 38.5, 39.5] vs. 14.8 [14.7, 14.9]), and had lower two-dose vaccination coverage (23.8% [23.6, 23.9] vs. 32.7% [32.6, 32.7]) compared with LGAs in the top five deciles. LGAs that achieved 70% coverage two weeks or more after the infection peak experienced nearly twice the total infection burden (27.7 per 1000 [27.3, 28.0] compared with 14.9 [14.7, 15.0]) than LGAs with a shorter lag.

CONCLUSIONS

Exposure and transmission risk factors cluster within disadvantaged LGAs. The potential for large local outbreaks is heightened if vaccination uptake trails in these communities.

IMPLICATIONS FOR PUBLIC HEALTH

In a pandemic, decision-makers must prioritise disease control and harm reduction interventions for at-risk LGAs.

Jointly estimating epidemiological dynamics of Covid-19 from case and wastewater data in Aotearoa New Zealand

BACKGROUND

Timely and informed public health responses to infectious diseases such as COVID-19 necessitate reliable information about infection dynamics. The case ascertainment rate (CAR), the proportion of infections that are reported as cases, is typically much less than one and varies with testing practices and behaviours, making reported cases unreliable as the sole source of data. The concentration of viral RNA in wastewater samples provides an alternate measure of infection prevalence that is not affected by clinical testing, healthcare-seeking behaviour or access to care.

METHODS

We construct a state-space model with observed data of levels of SARS-CoV-2 in wastewater and reported case incidence and estimate the hidden states of the effective reproduction number, R, and CAR using sequential Monte Carlo methods.

RESULTS

We analyse data from 1 January 2022 to 31 March 2023 from Aotearoa New Zealand. Our model estimates that R peaks at 2.76 (95% CrI 2.20, 3.83) around 18 February 2022 and the CAR peaks around 12 March 2022. We calculate that New Zealand's second Omicron wave in July 2022 is similar in size to the first, despite fewer reported cases. We estimate that the CAR in the BA.5 Omicron wave in July 2022 is approximately 50% lower than in the BA.1/BA.2 Omicron wave in March 2022.

CONCLUSIONS

Estimating R, CAR, and cumulative number of infections provides useful information for planning public health responses and understanding the state of immunity in the population. This model is a useful disease surveillance tool, improving situational awareness of infectious disease dynamics in real-time.

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.

Is SARS-CoV-2 elimination or mitigation best? Regional and disease characteristics determine the recommended strategy

Public health responses to the COVID-19 pandemic varied across the world. Some countries (e.g. mainland China, New Zealand and Taiwan) implemented elimination strategies involving strict travel measures and periods of rigorous non-pharmaceutical interventions (NPIs) in the community, aiming to achieve periods with no disease spread; while others (e.g. many European countries and the USA) implemented mitigation strategies involving less strict NPIs for prolonged periods, aiming to limit community spread. Travel measures and community NPIs have high economic and social costs, and there is a need for guidelines that evaluate the appropriateness of an elimination or mitigation strategy in regional contexts. To guide decisions, we identify key criteria and provide indicators and visualizations to help answer each question. Considerations include determining whether disease elimination is: (1) necessary to ensure healthcare provision; (2) feasible from an epidemiological point of view and (3) cost-effective when considering, in particular, the economic costs of travel measures and treating infections. We discuss our recommendations by considering the regional and economic variability of Canadian provinces and territories, and the epidemiological characteristics of different SARS-CoV-2 variants. While elimination may be a preferable strategy for regions with limited healthcare capacity, low travel volumes, and few ports of entry, mitigation may be more feasible in large urban areas with dense infrastructure, strong economies, and with high connectivity to other regions.

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.

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.

Near-term forecasting of Covid-19 cases and hospitalisations in Aotearoa New Zealand

Near-term forecasting of infectious disease incidence and consequent demand for acute healthcare services can support capacity planning and public health responses. Despite well-developed scenario modelling to support the Covid-19 response, Aotearoa New Zealand lacks advanced infectious disease forecasting capacity. We develop a model using Aotearoa New Zealand's unique Covid-19 data streams to predict reported Covid-19 cases, hospital admissions and hospital occupancy. The method combines a semi-mechanistic model for disease transmission to predict cases with Gaussian process regression models to predict the fraction of reported cases that will require hospital treatment. We evaluate forecast performance against out-of-sample data over the period from 2 October 2022 to 23 July 2023. Our results show that forecast performance is reasonably good over a 1-3 week time horizon, although generally deteriorates as the time horizon is lengthened. The model has been operationalised to provide weekly national and regional forecasts in real-time. This study is an important step towards development of more sophisticated situational awareness and infectious disease forecasting tools in Aotearoa New Zealand.

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.

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.

Influence of setting-dependent contacts and protective behaviours on asymptomatic SARS-CoV-2 infection amongst members of a UK university

We survey 62 users of a university asymptomatic SARS-CoV-2 testing service on details of their activities, protective behaviours and contacts in the 7 days prior to receiving a positive or negative SARS-CoV-2 PCR test result in the period October 2020-March 2021. The resulting data set is novel in capturing very detailed social contact history linked to asymptomatic disease status during a period of significant restriction on social activities. We use this data to explore 3 questions: (i) Did participation in university activities enhance infection risk? (ii) How do contact definitions rank in their ability to explain test outcome during periods of social restrictions? (iii) Do patterns in the protective behaviours help explain discrepancies between the explanatory performance of different contact measures? We classify activities into settings and use Bayesian logistic regression to model test outcome, computing posterior model probabilities to compare the performance of models adopting different contact definitions. Associations between protective behaviours, participant characteristics and setting are explored at the level of individual activities using multiple correspondence analysis (MCA). We find that participation in air travel or non-university work activities was associated with a positive asymptomatic SARS-CoV-2 PCR test, in contrast to participation in research and teaching settings. Intriguingly, logistic regression models with binary measures of contact in a setting performed better than more traditional contact numbers or person contact hours (PCH). The MCA indicates that patterns of protective behaviours vary between setting, in a manner which may help explain the preference for any participation as a contact measure. We conclude that linked PCR testing and social contact data can in principle be used to test the utility of contact definitions, and the investigation of contact definitions in larger linked studies is warranted to ensure contact data can capture environmental and social factors influencing transmission risk.

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.