HPVsim: An agent-based model of HPV transmission and cervical disease.

In 2020, the WHO launched its first global strategy to accelerate the elimination of cervical cancer, outlining an ambitious set of targets for countries to achieve over the next decade. At the same time, new tools, technologies, and strategies are in the pipeline that may improve screening performance, expand the reach of prophylactic vaccines, and prevent the acquisition, persistence and progression of oncogenic HPV. Detailed mechanistic modelling can help identify the combinations of current and future strategies to combat cervical cancer. Open-source modelling tools are needed to shift the capacity for such evaluations in-country. Here, we introduce the Human papillomavirus simulator (HPVsim), a new open-source software package for creating flexible agent-based models parameterised with country-specific vital dynamics, structured sexual networks, and co-transmitting HPV genotypes. HPVsim includes a novel methodology for modelling cervical disease progression, designed to be readily adaptable to new forms of screening. The software itself is implemented in Python, has built-in tools for simulating commonly-used interventions, includes a comprehensive set of tests and documentation, and runs quickly (seconds to minutes) on a laptop. Performance is greatly enhanced by HPVsim's multiscale modelling functionality. HPVsim is open source under the MIT License and available via both the Python Package Index (via pip install) and GitHub (hpvsim.org).

COVID-19 epidemic modelling for policy decision support in Victoria, Australia 2020-2021.

BACKGROUND: Policy responses to COVID-19 in Victoria, Australia over 2020-2021 have been supported by evidence generated through mathematical modelling. This study describes the design, key findings, and process for policy translation of a series of modelling studies conducted for the Victorian Department of Health COVID-19 response team during this period. METHODS: An agent-based model, Covasim, was used to simulate the impact of policy interventions on COVID-19 outbreaks and epidemic waves. The model was continually adapted to enable scenario analysis of settings or policies being considered at the time (e.g. elimination of community transmission versus disease control). Model scenarios were co-designed with government, to fill evidence gaps prior to key decisions. RESULTS: Understanding outbreak risk following incursions was critical to eliminating community COVID-19 transmission. Analyses showed risk depended on whether the first detected case was the index case, a primary contact of the index case, or a 'mystery case'. There were benefits of early lockdown on first case detection and gradual easing of restrictions to minimise resurgence risk from undetected cases. As vaccination coverage increased and the focus shifted to controlling rather than eliminating community transmission, understanding health system demand was critical. Analyses showed that vaccines alone could not protect health systems and need to be complemented with other public health measures. CONCLUSIONS: Model evidence offered the greatest value when decisions needed to be made pre-emptively, or for questions that could not be answered with empiric data and data analysis alone. Co-designing scenarios with policy-makers ensured relevance and increased policy translation.

Optima TB: A tool to help optimally allocate tuberculosis spending.

Approximately 85% of tuberculosis (TB) related deaths occur in low- and middle-income countries where health resources are scarce. Effective priority setting is required to maximise the impact of limited budgets. The Optima TB tool has been developed to support analytical capacity and inform evidence-based priority setting processes for TB health benefits package design. This paper outlines the Optima TB framework and how it was applied in Belarus, an upper-middle income country in Eastern Europe with a relatively high burden of TB. Optima TB is a population-based disease transmission model, with programmatic cost functions and an optimisation algorithm. Modelled populations include age-differentiated general populations and higher-risk populations such as people living with HIV. Populations and prospective interventions are defined in consultation with local stakeholders. In partnership with the latter, demographic, epidemiological, programmatic, as well as cost and spending data for these populations and interventions are then collated. An optimisation analysis of TB spending was conducted in Belarus, using program objectives and constraints defined in collaboration with local stakeholders, which included experts, decision makers, funders and organisations involved in service delivery, support and technical assistance. These analyses show that it is possible to improve health impact by redistributing current TB spending in Belarus. Specifically, shifting funding from inpatient- to outpatient-focused care models, and from mass screening to active case finding strategies, could reduce TB prevalence and mortality by up to 45% and 50%, respectively, by 2035. In addition, an optimised allocation of TB spending could lead to a reduction in drug-resistant TB infections by 40% over this period. This would support progress towards national TB targets without additional financial resources. The case study in Belarus demonstrates how reallocations of spending across existing and new interventions could have a substantial impact on TB outcomes. This highlights the potential for Optima TB and similar modelling tools to support evidence-based priority setting.

Covasim: An agent-based model of COVID-19 dynamics and interventions.

The COVID-19 pandemic has created an urgent need for models that can project epidemic trends, explore intervention scenarios, and estimate resource needs. Here we describe the methodology of Covasim (COVID-19 Agent-based Simulator), an open-source model developed to help address these questions. Covasim includes country-specific demographic information on age structure and population size; realistic transmission networks in different social layers, including households, schools, workplaces, long-term care facilities, and communities; age-specific disease outcomes; and intrahost viral dynamics, including viral-load-based transmissibility. Covasim also supports an extensive set of interventions, including non-pharmaceutical interventions, such as physical distancing and protective equipment; pharmaceutical interventions, including vaccination; and testing interventions, such as symptomatic and asymptomatic testing, isolation, contact tracing, and quarantine. These interventions can incorporate the effects of delays, loss-to-follow-up, micro-targeting, and other factors. Implemented in pure Python, Covasim has been designed with equal emphasis on performance, ease of use, and flexibility: realistic and highly customized scenarios can be run on a standard laptop in under a minute. In collaboration with local health agencies and policymakers, Covasim has already been applied to examine epidemic dynamics and inform policy decisions in more than a dozen countries in Africa, Asia-Pacific, Europe, and North America.

Preventing a cluster from becoming a new wave in settings with zero community COVID-19 cases.

BACKGROUND: In settings with zero community transmission, any new SARS-CoV-2 outbreaks are likely to be the result of random incursions. The level of restrictions in place at the time of the incursion is likely to considerably affect possible outbreak trajectories, but the probability that a large outbreak eventuates is not known. METHODS: We used an agent-based model to investigate the relationship between ongoing restrictions and behavioural factors, and the probability of an incursion causing an outbreak and the resulting growth rate. We applied our model to the state of Victoria, Australia, which has reached zero community transmission as of November 2020. RESULTS: We found that a future incursion has a 45% probability of causing an outbreak (defined as a 7-day average of > 5 new cases per day within 60 days) if no restrictions were in place, decreasing to 23% with a mandatory masks policy, density restrictions on venues such as restaurants, and if employees worked from home where possible. A drop in community symptomatic testing rates was associated with up to a 10-percentage point increase in outbreak probability, highlighting the importance of maintaining high testing rates as part of a suppression strategy. CONCLUSIONS: Because the chance of an incursion occurring is closely related to border controls, outbreak risk management strategies require an integrated approaching spanning border controls, ongoing restrictions, and plans for response. Each individual restriction or control strategy reduces the risk of an outbreak. They can be traded off against each other, but if too many are removed there is a danger of accumulating an unsafe level of risk. The outbreak probabilities estimated in this study are of particular relevance in assessing the downstream risks associated with increased international travel.

Modelling herd immunity requirements in Queensland: impact of vaccination effectiveness, hesitancy and variants of SARS-CoV-2.

Long-term control of SARS-CoV-2 outbreaks depends on the widespread coverage of effective vaccines. In Australia, two-dose vaccination coverage of above 90% of the adult population was achieved. However, between August 2020 and August 2021, hesitancy fluctuated dramatically. This raised the question of whether settings with low naturally derived immunity, such as Queensland where less than [Formula: see text] of the population is known to have been infected in 2020, could have achieved herd immunity against 2021's variants of concern. To address this question, we used the agent-based model Covasim. We simulated outbreak scenarios (with the Alpha, Delta and Omicron variants) and assumed ongoing interventions (testing, tracing, isolation and quarantine). We modelled vaccination using two approaches with different levels of realism. Hesitancy was modelled using Australian survey data. We found that with a vaccine effectiveness against infection of 80%, it was possible to control outbreaks of Alpha, but not Delta or Omicron. With 90% effectiveness, Delta outbreaks may have been preventable, but not Omicron outbreaks. We also estimated that a decrease in hesitancy from 20% to 14% reduced the number of infections, hospitalizations and deaths by over 30%. Overall, we demonstrate that while herd immunity may not be attainable, modest reductions in hesitancy and increases in vaccine uptake may greatly improve health outcomes. This article is part of the theme issue 'Technical challenges of modelling real-life epidemics and examples of overcoming these'.

A framework for considering the utility of models when facing tough decisions in public health: a guideline for policy-makers.

The COVID-19 pandemic has brought the combined disciplines of public health, infectious disease and policy modelling squarely into the spotlight. Never before have decisions regarding public health measures and their impacts been such a topic of international deliberation, from the level of individuals and communities through to global leaders. Nor have models-developed at rapid pace and often in the absence of complete information-ever been so central to the decision-making process. However, after nearly 3 years of experience with modelling, policy-makers need to be more confident about which models will be most helpful to support them when taking public health decisions, and modellers need to better understand the factors that will lead to successful model adoption and utilization. We present a three-stage framework for achieving these ends.

Modelling the impact of relaxing COVID-19 control measures during a period of low viral transmission.

OBJECTIVES: To assess the risks associated with relaxing coronavirus disease 2019 (COVID-19)-related physical distancing restrictions and lockdown policies during a period of low viral transmission. DESIGN: Network-based viral transmission risks in households, schools, workplaces, and a variety of community spaces and activities were simulated in an agent-based model, Covasim. SETTING: The model was calibrated for a baseline scenario reflecting the epidemiological and policy environment in Victoria during March-May 2020, a period of low community viral transmission. INTERVENTION: Policy changes for easing COVID-19-related restrictions from May 2020 were simulated in the context of interventions that included testing, contact tracing (including with a smartphone app), and quarantine. MAIN OUTCOME MEASURE: Increase in detected COVID-19 cases following relaxation of restrictions. RESULTS: Policy changes that facilitate contact of individuals with large numbers of unknown people (eg, opening bars, increased public transport use) were associated with the greatest risk of COVID-19 case numbers increasing; changes leading to smaller, structured gatherings with known contacts (eg, small social gatherings, opening schools) were associated with lower risks. In our model, the rise in case numbers following some policy changes was notable only two months after their implementation. CONCLUSIONS: Removing several COVID-19-related restrictions within a short period of time should be undertaken with care, as the consequences may not be apparent for more than two months. Our findings support continuation of work from home policies (to reduce public transport use) and strategies that mitigate the risk associated with re-opening of social venues.

Diabetes care cascade in Ukraine: an analysis of breakpoints and opportunities for improved diabetes outcomes.

BACKGROUND: Diabetes is one of the leading causes of poor health and high care costs in Ukraine. To prevent diabetes complications and alleviate the financial burden of diabetes care on patients, the Ukrainian government reimburses diabetes medication and provides glucose monitoring, but there are significant gaps in the care continuum. We estimate the costs of providing diabetes care and the most cost-effective ways to address these gaps in the Poltava region of Ukraine. METHODS: We gathered data on the unit costs of diabetes interventions in Poltava and estimated expenditure on diabetes care. We estimated the optimal combination of facility-based and outreach screening and investigated how additional funding could best be allocated to improve glucose control outcomes. RESULTS: Of the ~ 40,000 adults in diabetes care, only ~ 25% achieved sustained glucose control. Monitoring costs were higher for those who did not: by 10% for patients receiving non-pharmacological treatment, by 61% for insulin patients, and twice as high for patients prescribed oral treatment. Initiatives to improve treatment adherence (e.g. medication copayment schemes, enhanced adherence counseling) would address barriers along the care continuum and we estimate such expenditures may be recouped by reductions in patient monitoring costs. Improvements in case detection are also needed, with only around two-thirds of estimated cases having been diagnosed. Outreach screening campaigns could play a significant role: depending on how well-targeted and scalable such campaigns are, we estimate that 10-46% of all screening could be conducted via outreach, at a cost per positive patient identified of US$7.12-9.63. CONCLUSIONS: Investments to improve case detection and treatment adherence are the most efficient interventions for improved diabetes control in Poltava. Quantitative tools provide essential decision support for targeting investment to close the gaps in care.

How do spatially distinct frequency specific MEG networks emerge from one underlying structural connectome? The role of the structural eigenmodes.

Functional networks obtained from magnetoencephalography (MEG) from different frequency bands show distinct spatial patterns. It remains to be elucidated how distinct spatial patterns in MEG networks emerge given a single underlying structural network. Recent work has suggested that the eigenmodes of the structural network might serve as a basis set for functional network patterns in the case of functional MRI. Here, we take this notion further in the context of frequency band specific MEG networks. We show that a selected set of eigenmodes of the structural network can predict different frequency band specific networks in the resting state, ranging from delta (1-4 Hz) to the high gamma band (40-70 Hz). These predictions outperform predictions based from surrogate data, suggesting a genuine relationship between eigenmodes of the structural network and frequency specific MEG networks. We then show that the relevant set of eigenmodes can be excited in a network of neural mass models using linear stability analysis only by including delays. Excitation of an eigenmode in this context refers to a dynamic instability of a network steady state to a spatial pattern with a corresponding coherent temporal oscillation. Simulations verify the results from linear stability analysis and suggest that theta, alpha and beta band networks emerge very near to the bifurcation. The delta and gamma bands in the resting state emerges further away from the bifurcation. These results show for the first time how delayed interactions can excite the relevant set of eigenmodes that give rise to frequency specific functional connectivity patterns.

NFTsim: Theory and Simulation of Multiscale Neural Field Dynamics.

A user ready, portable, documented software package, NFTsim, is presented to facilitate numerical simulations of a wide range of brain systems using continuum neural field modeling. NFTsim enables users to simulate key aspects of brain activity at multiple scales. At the microscopic scale, it incorporates characteristics of local interactions between cells, neurotransmitter effects, synaptodendritic delays and feedbacks. At the mesoscopic scale, it incorporates information about medium to large scale axonal ranges of fibers, which are essential to model dissipative wave transmission and to produce synchronous oscillations and associated cross-correlation patterns as observed in local field potential recordings of active tissue. At the scale of the whole brain, NFTsim allows for the inclusion of long range pathways, such as thalamocortical projections, when generating macroscopic activity fields. The multiscale nature of the neural activity produced by NFTsim has the potential to enable the modeling of resulting quantities measurable via various neuroimaging techniques. In this work, we give a comprehensive description of the design and implementation of the software. Due to its modularity and flexibility, NFTsim enables the systematic study of an unlimited number of neural systems with multiple neural populations under a unified framework and allows for direct comparison with analytic and experimental predictions. The code is written in C++ and bundled with Matlab routines for a rapid quantitative analysis and visualization of the outputs. The output of NFTsim is stored in plain text file enabling users to select from a broad range of tools for offline analysis. This software enables a wide and convenient use of powerful physiologically-based neural field approaches to brain modeling. NFTsim is distributed under the Apache 2.0 license.

A biophysical model of dynamic balancing of excitation and inhibition in fast oscillatory large-scale networks.

Over long timescales, neuronal dynamics can be robust to quite large perturbations, such as changes in white matter connectivity and grey matter structure through processes including learning, aging, development and certain disease processes. One possible explanation is that robust dynamics are facilitated by homeostatic mechanisms that can dynamically rebalance brain networks. In this study, we simulate a cortical brain network using the Wilson-Cowan neural mass model with conduction delays and noise, and use inhibitory synaptic plasticity (ISP) to dynamically achieve a spatially local balance between excitation and inhibition. Using MEG data from 55 subjects we find that ISP enables us to simultaneously achieve high correlation with multiple measures of functional connectivity, including amplitude envelope correlation and phase locking. Further, we find that ISP successfully achieves local E/I balance, and can consistently predict the functional connectivity computed from real MEG data, for a much wider range of model parameters than is possible with a model without ISP.

Discovering dynamic brain networks from big data in rest and task.

Brain activity is a dynamic combination of the responses to sensory inputs and its own spontaneous processing. Consequently, such brain activity is continuously changing whether or not one is focusing on an externally imposed task. Previously, we have introduced an analysis method that allows us, using Hidden Markov Models (HMM), to model task or rest brain activity as a dynamic sequence of distinct brain networks, overcoming many of the limitations posed by sliding window approaches. Here, we present an advance that enables the HMM to handle very large amounts of data, making possible the inference of very reproducible and interpretable dynamic brain networks in a range of different datasets, including task, rest, MEG and fMRI, with potentially thousands of subjects. We anticipate that the generation of large and publicly available datasets from initiatives such as the Human Connectome Project and UK Biobank, in combination with computational methods that can work at this scale, will bring a breakthrough in our understanding of brain function in both health and disease.

A unified model of melatonin, 6-sulfatoxymelatonin, and sleep dynamics.

A biophysical model of the key aspects of melatonin synthesis and excretion has been developed, which is able to predict experimental dynamics of melatonin in plasma and saliva, and of its urinary metabolite 6-sulfatoxymelatonin (aMT6s). This new model is coupled to an established model of arousal dynamics, which predicts sleep and circadian dynamics based on light exposure and times of wakefulness. The combined model thus predicts melatonin levels over the sleep-wake/dark-light cycle and enables prediction of melatonin-based circadian phase markers, such as dim light melatonin onset (DLMO) and aMT6s acrophase under conditions of normal sleep and circadian misalignment. The model is calibrated and tested against group average data from 10 published experimental studies and is found to reproduce quantitatively the key dynamics of melatonin and aMT6s, including the timing of release and amplitude, as well as response to controlled lighting and shift work.

Inferring the natural history of HPV from global cancer registries: insights from a multi-country calibration.

Human papillomavirus (HPV) is the cause of almost all cases of cervical cancer, a disease that kills some 340,000 women per year. The timeline from initial infection with HPV to the onset of invasive cervical cancer spans decades, and observational studies of this process are limited to settings in which treatment of precancerous lesions was withheld or inadequate. Such studies have been critical for understanding the natural history of HPV. Modeling can shed additional insight on the natural history of HPV, especially across geographical settings with varying prevalence of factors known to affect the host-side immune response to HPV, such as HIV and tobacco use. In this study, we create models for the 30 most populous countries in Sub-Saharan Africa, each with country-specific demographic, and behavioral inputs. We found that it was not possible to fit the data if we assumed that the natural history parameters were exactly identical for all countries, even after accounting for demographic and behavioral differences, but that we could achieve a good fit with the addition of a single immunocompetence parameter for each country. Our results indicate that variation in host immune responses may play a role in explaining the differences in the burden of cervical cancer between countries, which in turn implies a greater need for more geographically diverse data collection to understand the natural history of HPV.

A Framework for Assessing the Impact of Outbreak Response Immunization Programs.

The impact of outbreak response immunization (ORI) can be estimated by comparing observed outcomes to modelled counterfactual scenarios without ORI, but the most appropriate metrics depend on stakeholder needs and data availability. This study developed a framework for using mathematical models to assess the impact of ORI for vaccine-preventable diseases. Framework development involved (1) the assessment of impact metrics based on stakeholder interviews and literature reviews determining data availability and capacity to capture as model outcomes; (2) mapping investment in ORI elements to model parameters to define scenarios; (3) developing a system for engaging stakeholders and formulating model questions, performing analyses, and interpreting results; and (4) example applications for different settings and pathogens. The metrics identified as most useful were health impacts, economic impacts, and the risk of severe outbreaks. Scenario categories included investment in the response scale, response speed, and vaccine targeting. The framework defines four phases: (1) problem framing and data sourcing (identification of stakeholder needs, metrics, and scenarios); (2) model choice; (3) model implementation; and (4) interpretation and communication. The use of the framework is demonstrated by application to two outbreaks, measles in Papua New Guinea and Ebola in the Democratic Republic of the Congo. The framework is a systematic way to engage with stakeholders and ensure that an analysis is fit for purpose, makes the best use of available data, and uses suitable modelling methodology.

COVID-19 outbreaks in residential aged care facilities: an agent-based modeling study.

Introduction: A disproportionate number of COVID-19 deaths occur in Residential Aged Care Facilities (RACFs), where better evidence is needed to target COVID-19 interventions to prevent mortality. This study used an agent-based model to assess the role of community prevalence, vaccination strategies, and non-pharmaceutical interventions (NPIs) on COVID-19 outcomes in RACFs in Victoria, Australia. Methods: The model simulated outbreaks in RACFs over time, and was calibrated to distributions for outbreak size, outbreak duration, and case fatality rate in Victorian RACFs over 2022. The number of incursions to RACFs per day were estimated to fit total deaths and diagnoses over time and community prevalence.Total infections, diagnoses, and deaths in RACFs were estimated over July 2023-June 2024 under scenarios of different: community epidemic wave assumptions (magnitude and frequency); RACF vaccination strategies (6-monthly, 12-monthly, no further vaccines); additional non-pharmaceutical interventions (10, 25, 50% efficacy); and reduction in incursions (30% or 60%). Results: Total RACF outcomes were proportional to cumulative community infections and incursion rates, suggesting potential for strategic visitation/staff policies or community-based interventions to reduce deaths. Recency of vaccination when epidemic waves occurred was critical; compared with 6-monthly boosters, 12-monthly boosters had approximately 1.2 times more deaths and no further boosters had approximately 1.6 times more deaths over July 2023-June 2024. Additional NPIs, even with only 10-25% efficacy, could lead to a 13-31% reduction in deaths in RACFs. Conclusion: Future community epidemic wave patterns are unknown but will be major drivers of outcomes in RACFs. Maintaining high coverage of recent vaccination, minimizing incursions, and increasing NPIs can have a major impact on cumulative infections and deaths.

Admission screening testing of patients and staff N95 masks are cost-effective in reducing COVID-19 hospital acquired infections.

BACKGROUND: COVID-19 outbreaks in acute care settings can have severe consequences for patients due to their underlying vulnerabilities, and can be costly due to additional patient bed days and the need to replace isolating staff. This study assessed the cost-effectiveness of clinical staff N95 masks and admission screening testing of patients to reduce COVID-19 hospital-acquired infections. METHODS: An agent-based model was calibrated to data on 178 outbreaks in acute care settings in Victoria, Australia between October 2021 and July 2023. Outbreaks were simulated under different combinations of staff masking (surgical, N95) and patient admission screening testing (none, RAT, PCR). For each scenario, average diagnoses, COVID-19 deaths, quality-adjusted life years (QALYs) from discharged patients, and costs (masks, testing, patient COVID-19 bed days, staff replacement costs while isolating) from acute COVID-19 were estimated over a 12-month period. FINDINGS: Compared to no admission screening testing and staff surgical masks, all scenarios were cost saving with health gains. Staff N95s + RAT admission screening of patients was the cheapest, saving A$78.4M [95%UI 44.4M-135.3M] and preventing 1,543 [1,070-2,146] deaths state-wide per annum. Both interventions were individually beneficial: staff N95s in isolation saved A$54.7M and 854 deaths state-wide per annum, while RAT admission screening of patients in isolation saved A$57.6M and 1,176 deaths state-wide per annum. INTERPRETATION: In acute care settings, staff N95 mask use and admission screening testing of patients can reduce hospital-acquired COVID-19 infections, COVID-19 deaths, and are cost-saving because of reduced patient bed days and staff replacement needs.

A global investment case for hepatitis B elimination: a modelling study.

Hepatitis B is estimated to cause 500 000-900 000 deaths globally each year. WHO has targets for elimination by 2030; however, progress has stalled due to multiple barriers, notably a paucity of global funding and insufficient evidence on the economic burden of disease. Using a dynamic mathematical model of hepatitis B transmission, disease progression, and mortality in the six WHO regions, we estimate the costs and benefits of reaching 90% vaccination, 90% diagnosis, and 80% treatment coverage by either 2030 (as targeted), 2040, or 2050. Without increased intervention coverage, hepatitis B mortality was estimated to cost US$784·35 billion (95% Crl 731·63-798·33 billion) globally in lost productivity over 2022-50. Achieving targets by 2030 averted 25·64 million infections (95% Crl 17·39-34·55 million) and 8·63 million hepatitis B-attributable deaths (95% Crl 7·12-9·74 million) over 2022-50. This achievement incurred an incremental cost of $2934·55 (95% Crl 2778·55-3173·52) per disability-adjusted life year averted by 2050 under a health systems perspective, and was cost-saving with a net economic benefit of $99·03 billion (95% Crl 78·66-108·96 billion) by 2050 from a societal perspective. Delayed achievement of intervention coverage targets had reduced health and economic benefits. These findings highlight that hepatitis B is an underappreciated cause of economic burden and show investment toward elimination will probably yield substantial returns.

Keeping kids in school: modelling school-based testing and quarantine strategies during the COVID-19 pandemic in Australia.

Background: In 2021, the Australian Government Department of Health commissioned a consortium of modelling groups to generate evidence assisting the transition from a goal of no community COVID-19 transmission to 'living with COVID-19', with adverse health and social consequences limited by vaccination and other measures. Due to the extended school closures over 2020-21, maximizing face-to-face teaching was a major objective during this transition. The consortium was tasked with informing school surveillance and contact management strategies to minimize infections and support this goal. Methods: Outcomes considered were infections and days of face-to-face teaching lost in the 45 days following an outbreak within an otherwise COVID-naïve school setting. A stochastic agent-based model of COVID-19 transmission was used to evaluate a 'test-to-stay' strategy using daily rapid antigen tests (RATs) for close contacts of a case for 7 days compared with home quarantine; and an asymptomatic surveillance strategy involving twice-weekly screening of all students and/or teachers using RATs. Findings: Test-to-stay had similar effectiveness for reducing school infections as extended home quarantine, without the associated days of face-to-face teaching lost. Asymptomatic screening was beneficial in reducing both infections and days of face-to-face teaching lost and was most beneficial when community prevalence was high. Interpretation: Use of RATs in school settings for surveillance and contact management can help to maximize face-to-face teaching and minimize outbreaks. This evidence supported the implementation of surveillance testing in schools in several Australian jurisdictions from January 2022.