Reproducibility of COVID-era infectious disease models.

Infectious disease modelling has been prominent throughout the COVID-19 pandemic, helping to understand the virus' transmission dynamics and inform response policies. Given their potential importance and translational impact, we evaluated the computational reproducibility of infectious disease modelling articles from the COVID era. We found that four out of 100 randomly sampled studies released between January 2020 and August 2022 could be completely computationally reproduced using the resources provided (e.g., code, data, instructions) whilst a further eight were partially reproducible. For the 100 most highly cited articles from the same period we found that 11 were completely reproducible with a further 22 partially reproducible. Reflecting on our experience, we discuss common issues affecting computational reproducibility and how these might be addressed.

Replicating superspreader dynamics with compartmental models.

Infectious disease outbreaks often exhibit superspreader dynamics, where most infected people generate no, or few secondary cases, and only a small fraction of individuals are responsible for a large proportion of transmission. Although capturing this heterogeneity is critical for estimating outbreak risk and the effectiveness of group-specific interventions, it is typically neglected in compartmental models of infectious disease transmission-which constitute the most common transmission dynamic modeling framework. In this study we propose different classes of compartmental epidemic models that incorporate transmission heterogeneity, fit them to a number of real outbreak datasets, and benchmark their performance against the canonical superspreader model (i.e., the negative binomial branching process model). We find that properly constructed compartmental models can capably reproduce observed superspreader dynamics and we provide the pathogen-specific parameter settings required to do so. As a consequence, we also show that compartmental models parameterized according to a binary clinical classification have limited support.

Quantifying the impact of Wolbachia releases on dengue infection in Townsville, Australia.

From October 2014 to February 2019, local authorities in Townsville, North Queensland, Australia continually introduced Wolbachia-infected mosquitoes to control seasonal outbreaks of dengue infection. In this study, we develop a mathematical modelling framework to estimate the effectiveness of this intervention as well as the relative dengue transmission rates of Wolbachia-infected and wild-type mosquitoes. We find that the transmission rate of Wolbachia-infected mosquitoes is reduced approximately by a factor of 20 relative to the uninfected wild-type population. In addition, the Townsville Wolbachia release program led to a 65% reduction in predicted dengue incidence during the release period and over 95% reduction in the 24 months that followed. Finally, to investigate the potential impact of other Wolbachia release programs, we use our estimates of relative transmissibility to calculate the relationship between the reproductive number of dengue and the proportion of Wolbachia-infected mosquitoes in the vector population.

A Systematic Review of Mathematical Models of Dengue Transmission and Vector Control: 2010-2020.

Vector control methods are considered effective in averting dengue transmission. However, several factors may modify their impact. Of these controls, chemical methods, in the long run, may increase mosquitoes' resistance to chemicides, thereby decreasing control efficacy. The biological methods, which may be self-sustaining and very effective, could be hampered by seasonality or heatwaves (resulting in, e.g., loss of Wolbachia infection). The environmental methods that could be more effective than the chemical methods are under-investigated. In this study, a systematic review is conducted to explore the present understanding of the effectiveness of vector control approaches via dengue transmission models.

Modelling the ecological dynamics of mosquito populations with multiple co-circulating Wolbachia strains.

Wolbachia intracellular bacteria successfully reduce the transmissibility of arthropod-borne viruses (arboviruses) when introduced into virus-carrying vectors such as mosquitoes. Despite the progress made by introducing Wolbachia bacteria into the Aedes aegypti wild-type population to control arboviral infections, reports suggest that heat-induced loss-of-Wolbachia-infection as a result of climate change may reverse these gains. Novel, supplemental Wolbachia strains that are more resilient to increased temperatures may circumvent these concerns, and could potentially act synergistically with existing variants. In this article, we model the ecological dynamics among three distinct mosquito (sub)populations: a wild-type population free of any Wolbachia infection; an invading population infected with a particular Wolbachia strain; and a second invading population infected with a distinct Wolbachia strain from that of the first invader. We explore how the range of possible characteristics of each Wolbachia strain impacts mosquito prevalence. Further, we analyse the differential system governing the mosquito populations and the Wolbachia infection dynamics by computing the full set of basic and invasive reproduction numbers and use these to establish stability of identified equilibria. Our results show that releasing mosquitoes with two different strains of Wolbachia did not increase their prevalence, compared with a single-strain Wolbachia-infected mosquito introduction and only delayed Wolbachia dominance.

Phylodynamic signatures in the emergence of community-associated MRSA.

Community-associated, methicillin-resistant Staphylococcus aureus (MRSA) lineages have emerged in many geographically distinct regions around the world during the past 30 y. Here, we apply consistent phylodynamic methods across multiple community-associated MRSA lineages to describe and contrast their patterns of emergence and dissemination. We generated whole-genome sequencing data for the Australian sequence type (ST) ST93-MRSA-IV from remote communities in Far North Queensland and Papua New Guinea, and the Bengal Bay ST772-MRSA-V clone from metropolitan communities in Pakistan. Increases in the effective reproduction number (Re) and sustained transmission (Re > 1) coincided with spread of progenitor methicillin-susceptible S. aureus (MSSA) in remote northern Australian populations, dissemination of the ST93-MRSA-IV genotype into population centers on the Australian East Coast, and subsequent importation into the highlands of Papua New Guinea and Far North Queensland. Applying the same phylodynamic methods to existing lineage datasets, we identified common signatures of epidemic growth in the emergence and epidemiological trajectory of community-associated S. aureus lineages from America, Asia, Australasia, and Europe. Surges in Re were observed at the divergence of antibiotic-resistant strains, coinciding with their establishment in regional population centers. Epidemic growth was also observed among drug-resistant MSSA clades in Africa and northern Australia. Our data suggest that the emergence of community-associated MRSA in the late 20th century was driven by a combination of antibiotic-resistant genotypes and host epidemiology, leading to abrupt changes in lineage-wide transmission dynamics and sustained transmission in regional population centers.

Mathematical analysis of a two-strain tuberculosis model in Bangladesh.

Tuberculosis (TB) is an airborne infectious disease that causes millions of deaths worldwide each year (1.2 million people died in 2019). Alarmingly, several strains of the causative agent, Mycobacterium tuberculosis (MTB)-including drug-susceptible (DS) and drug-resistant (DR) variants-already circulate throughout most developing and developed countries, particularly in Bangladesh, with totally drug-resistant strains starting to emerge. In this study we develop a two-strain DS and DR TB transmission model and perform an analysis of the system properties and solutions. Both analytical and numerical results show that the prevalence of drug-resistant infection increases with an increasing drug use through amplification. Both analytic results and numerical simulations suggest that if the basic reproduction numbers of both DS ([Formula: see text]) and DR ([Formula: see text]) TB are less than one, i.e. [Formula: see text] the disease-free equilibrium is asymptotically stable, meaning that the disease naturally dies out. Furthermore, if [Formula: see text], then DS TB dies out but DR TB persists in the population, and if [Formula: see text] both DS TB and DR TB persist in the population. Further, sensitivity analysis of the model parameters found that the transmission rate of both strains had the greatest influence on DS and DR TB prevalence. We also investigated the effect of treatment rates and amplification on both DS and DR TB prevalence; results indicate that inadequate or inappropriate treatment makes co-existence more likely and increases the relative abundance of DR TB infections.

Estimating the global burden of Epstein-Barr virus-related cancers.

BACKGROUND: More than 90% of the adult population globally is chronically infected by the Epstein-Barr virus (EBV). It is well established that EBV is associated with a number of malignancies, and advances in knowledge of EBV-related malignancies are being made every year. Several studies have analysed the global epidemiology and geographic distribution of EBV-related cancers. However, most have only described a single cancer type or subtype in isolation or limited their study to the three or four most common EBV-related cancers. This review will present an overview on the spectrum of cancers linked to EBV based on observations of associations and proportions in the published literature while also using these observations to estimate the incidence and mortality burden of some of these cancers. METHOD: We have reviewed the literature on defining features, distribution and outcomes across six cancers with a relatively large EBV-related case burden: Nasopharyngeal carcinoma (NPC), Gastric carcinoma (GC), Hodgkin lymphoma (HL), Burkitt lymphoma (BL), Diffuse large B-cell lymphoma (DLBCL) and Extranodal NK/T-cell lymphoma, Nasal type (ENKTL-NT). We retrieved published region-specific EBV-related case proportions for NPC, GC, HL and BL and performed meta-analyses on pooled region-specific studies of EBV-related case proportions for DLBCL and ENKTL-NT. We match these pooled proportions with their respective regional incidence and mortality numbers retrieved from a publicly available cancer database. Additionally, we also reviewed the literature on several other less common EBV-related cancers to summarize their key characteristics herein. CONCLUSION: We estimated that EBV-related cases from these six cancers accounted for 239,700-357,900 new cases and 137,900-208,700 deaths in 2020. This review highlights the significant global impact of EBV-related cancers and extends the spectrum of disease that could benefit from an EBV-specific therapeutic.

Sustaining effective COVID-19 control in Malaysia through large-scale vaccination.

INTRODUCTION: As of 3rd June 2021, Malaysia is experiencing a resurgence of COVID-19 cases. In response, the federal government has implemented various non-pharmaceutical interventions (NPIs) under a series of Movement Control Orders and, more recently, a vaccination campaign to regain epidemic control. In this study, we assessed the potential for the vaccination campaign to control the epidemic in Malaysia and four high-burden regions of interest, under various public health response scenarios. METHODS: A modified susceptible-exposed-infectious-recovered compartmental model was developed that included two sequential incubation and infectious periods, with stratification by clinical state. The model was further stratified by age and incorporated population mobility to capture NPIs and micro-distancing (behaviour changes not captured through population mobility). Emerging variants of concern (VoC) were included as an additional strain competing with the existing wild-type strain. Several scenarios that included different vaccination strategies (i.e. vaccines that reduce disease severity and/or prevent infection, vaccination coverage) and mobility restrictions were implemented. RESULTS: The national model and the regional models all fit well to notification data but underestimated ICU occupancy and deaths in recent weeks, which may be attributable to increased severity of VoC or saturation of case detection. However, the true case detection proportion showed wide credible intervals, highlighting incomplete understanding of the true epidemic size. The scenario projections suggested that under current vaccination rates complete relaxation of all NPIs would trigger a major epidemic. The results emphasise the importance of micro-distancing, maintaining mobility restrictions during vaccination roll-out and accelerating the pace of vaccination for future control. Malaysia is particularly susceptible to a major COVID-19 resurgence resulting from its limited population immunity due to the country's historical success in maintaining control throughout much of 2020.

Understanding how Victoria, Australia gained control of its second COVID-19 wave.

During 2020, Victoria was the Australian state hardest hit by COVID-19, but was successful in controlling its second wave through aggressive policy interventions. We calibrated a detailed compartmental model of Victoria's second wave to multiple geographically-structured epidemic time-series indicators. We achieved a good fit overall and for individual health services through a combination of time-varying processes, including case detection, population mobility, school closures, physical distancing and face covering usage. Estimates of the risk of death in those aged ≥75 and of hospitalisation were higher than international estimates, reflecting concentration of cases in high-risk settings. We estimated significant effects for each of the calibrated time-varying processes, with estimates for the individual-level effect of physical distancing of 37.4% (95%CrI 7.2-56.4%) and of face coverings of 45.9% (95%CrI 32.9-55.6%). That the multi-faceted interventions led to the dramatic reversal in the epidemic trajectory is supported by our results, with face coverings likely particularly important.

Modelling direct and herd protection effects of vaccination against the SARS-CoV-2 Delta variant in Australia.

OBJECTIVES: To analyse the outcomes of COVID-19 vaccination by vaccine type, age group eligibility, vaccination strategy, and population coverage. DESIGN: Epidemiologic modelling to assess the final size of a COVID-19 epidemic in Australia, with vaccination program (Pfizer, AstraZeneca, mixed), vaccination strategy (vulnerable first, transmitters first, untargeted), age group eligibility threshold (5 or 15 years), population coverage, and pre-vaccination effective reproduction number ( Reffv¯ ) for the SARS-CoV-2 Delta variant as factors. MAIN OUTCOME MEASURES: Numbers of SARS-CoV-2 infections; cumulative hospitalisations, deaths, and years of life lost. RESULTS: Assuming Reffv¯ = 5, the current mixed vaccination program (vaccinating people aged 60 or more with the AstraZeneca vaccine and people under 60 with the Pfizer vaccine) will not achieve herd protection unless population vaccination coverage reaches 85% by lowering the vaccination eligibility age to 5 years. At Reffv¯ = 3, the mixed program could achieve herd protection at 60-70% population coverage and without vaccinating 5-15-year-old children. At Reffv¯ = 7, herd protection is unlikely to be achieved with currently available vaccines, but they would still reduce the number of COVID-19-related deaths by 85%. CONCLUSION: Vaccinating vulnerable people first is the optimal policy when population vaccination coverage is low, but vaccinating more socially active people becomes more important as the Reffv¯ declines and vaccination coverage increases. Assuming the most plausible Reffv¯ of 5, vaccinating more than 85% of the population, including children, would be needed to achieve herd protection. Even without herd protection, vaccines are highly effective in reducing the number of deaths.

Vaccines and variants: Modelling insights into emerging issues in COVID-19 epidemiology.

Mathematical modelling has played a pivotal role in understanding the epidemiology of and guiding public health responses to the ongoing coronavirus disease of 2019 (COVID-19) pandemic. Here, we review the role of epidemiological models in understanding evolving epidemic characteristics, including the effects of vaccination and Variants of Concern (VoC). We highlight ways in which models continue to provide important insights, including (1) calculating the herd immunity threshold and evaluating its limitations; (2) verifying that nascent vaccines can prevent severe disease, infection, and transmission but may be less efficacious against VoC; (3) determining optimal vaccine allocation strategies under efficacy and supply constraints; and (4) determining that VoC are more transmissible and lethal than previously circulating strains, and that immune escape may jeopardize vaccine-induced herd immunity. Finally, we explore how models can help us anticipate and prepare for future stages of COVID-19 epidemiology (and that of other diseases) through forecasts and scenario projections, given current uncertainties and data limitations.

Major amputation rates and outcomes for Aboriginal and Torres Strait Islander and non-Indigenous people in North Queensland Australia between 2000 and 2015.

BACKGROUND: This study estimated the incidence of major amputation for people in North Queensland, Australia, examined changes in amputation rates over time and investigated survival after major amputation. METHODS: This was a retrospective study of patients who underwent a major amputation above the ankle between 2000 and 2015. Major amputation rates and incidence rate ratios (IRR) were calculated using census data to define the at-risk population. Associations between risk factors and calendar year with major amputation were assessed using quasipoisson regression. Kaplan-Meier survival and Cox-proportional hazard analyses estimated the incidence of and risk factors for all-cause mortality. RESULTS: The annual incidence of major amputation was estimated to be greater in Aboriginal and Torres Strait Islanders than non-Indigenous people (IRR 2.75, 95 % CI 1.92 to 3.84). After adjusting for population growth, the annual incidence of major amputations did not change significantly over time for either groups. Aboriginal and Torres Strait Islander people were at greater risk of all-cause mortality after major amputation compared to non-Indigenous people, although this association was not significant after adjusting for other risk factors (hazard ratio 1.24, 95 % CI 0.82 to 1.90). CONCLUSIONS: The incidence of major amputation in North Queensland has not reduced over time, indicating the need for better preventative treatments, particularly in Aboriginal and Torres Strait Islander people.

Scenario analysis for programmatic tuberculosis control in Bangladesh: a mathematical modelling study.

Tuberculosis (TB) is a major public health problem in Bangladesh. Although the National TB control program of Bangladesh is implementing a comprehensive expansion of TB control strategies, logistical challenges exist, and there is significant uncertainty concerning the disease burden. Mathematical modelling of TB is considered one of the most effective ways to understand the dynamics of infection transmission and allows quantification of parameters in different settings, including Bangladesh. In this study, we present a two-strain mathematical modelling framework to explore the dynamics of drug-susceptible (DS) and multidrug-resistant (MDR) TB in Bangladesh. We calibrated the model using DS and MDR-TB annual incidence data from Bangladesh from years 2001 to 2015. Further, we performed a sensitivity analysis of the model parameters and found that the contact rate of both strains had the largest influence on the basic reproduction numbers [Formula: see text] and [Formula: see text] of DS and MDR-TB, respectively. Increasingly powerful intervention strategies were developed, with realistic impact and coverage determined with the help of local staff. We simulated for the period from 2020 to 2035. Here, we projected the DS and MDR-TB burden (as measured by the number of incident cases and mortality) under a range of intervention scenarios to determine which of these scenario is the most effective at reducing burden. Of the single-intervention strategies, enhanced case detection is the most effective and prompt in reducing DS and MDR-TB incidence and mortality in Bangladesh and that with GeneXpert testing was also highly effective in decreasing the burden of MDR-TB. Our findings also suggest combining additional interventions simultaneously leads to greater effectiveness, particularly for MDR-TB, which we estimate requires a modest investment to substantially reduce, whereas DS-TB requires a strong sustained investment.

A Review: Aedes-Borne Arboviral Infections, Controls and Wolbachia-Based Strategies.

Arthropod-borne viruses (Arboviruses) continue to generate significant health and economic burdens for people living in endemic regions. Of these viruses, some of the most important (e.g., dengue, Zika, chikungunya, and yellow fever virus), are transmitted mainly by Aedes mosquitoes. Over the years, viral infection control has targeted vector population reduction and inhibition of arboviral replication and transmission. This control includes the vector control methods which are classified into chemical, environmental, and biological methods. Some of these control methods may be largely experimental (both field and laboratory investigations) or widely practised. Perceptively, one of the biological methods of vector control, in particular, Wolbachia-based control, shows a promising control strategy for eradicating Aedes-borne arboviruses. This can either be through the artificial introduction of Wolbachia, a naturally present bacterium that impedes viral growth in mosquitoes into heterologous Aedes aegypti mosquito vectors (vectors that are not natural hosts of Wolbachia) thereby limiting arboviral transmission or via Aedes albopictus mosquitoes, which naturally harbour Wolbachia infection. These strategies are potentially undermined by the tendency of mosquitoes to lose Wolbachia infection in unfavourable weather conditions (e.g., high temperature) and the inhibitory competitive dynamics among co-circulating Wolbachia strains. The main objective of this review was to critically appraise published articles on vector control strategies and specifically highlight the use of Wolbachia-based control to suppress vector population growth or disrupt viral transmission. We retrieved studies on the control strategies for arboviral transmissions via arthropod vectors and discussed the use of Wolbachia control strategies for eradicating arboviral diseases to identify literature gaps that will be instrumental in developing models to estimate the impact of these control strategies and, in essence, the use of different Wolbachia strains and features.

Modeling the potential of wAu-Wolbachia strain invasion in mosquitoes to control Aedes-borne arboviral infections.

Arboviral infections such as dengue, Zika and chikungunya are fast spreading diseases that pose significant health problems globally. In order to control these infections, an intracellular bacterium called Wolbachia has been introduced into wild-type mosquito populations in the hopes of replacing the vector transmitting agent, Aedes aegypti with one that is incapable of transmission. In this study, we developed a Wolbachia transmission model for the novel wAu strain which possesses several favourable traits (e.g., enhanced viral blockage and maintenance at higher temperature) but not cyctoplasmic incompatibility (CI)-when a Wolbachia-infected male mosquito mates with an uninfected female mosquito, producing no viable offspring. This model describes the competitive dynamics between wAu-Wolbachia-infected and uninfected mosquitoes and the role of imperfect maternal transmission. By analysing the system via computing the basic reproduction number(s) and stability properties, the potential of the wAu strain as a viable strategy to control arboviral infections is established. The results of this work show that enhanced maintenance of Wolbachia infection at higher temperatures can overcome the lack of CI induction to support wAu-Wolbachia infected mosquito invasion. This study will support future arboviral control programs, that rely on the introduction of new Wolbachia variants.

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

Modelling insights into the COVID-19 pandemic.

Coronavirus disease 2019 (COVID-19) is a newly emerged infectious disease caused by the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) that was declared a pandemic by the World Health Organization on 11th March, 2020. Response to this ongoing pandemic requires extensive collaboration across the scientific community in an attempt to contain its impact and limit further transmission. Mathematical modelling has been at the forefront of these response efforts by: (1) providing initial estimates of the SARS-CoV-2 reproduction rate, R0 (of approximately 2-3); (2) updating these estimates following the implementation of various interventions (with significantly reduced, often sub-critical, transmission rates); (3) assessing the potential for global spread before significant case numbers had been reported internationally; and (4) quantifying the expected disease severity and burden of COVID-19, indicating that the likely true infection rate is often orders of magnitude greater than estimates based on confirmed case counts alone. In this review, we highlight the critical role played by mathematical modelling to understand COVID-19 thus far, the challenges posed by data availability and uncertainty, and the continuing utility of modelling-based approaches to guide decision making and inform the public health response. †Unless otherwise stated, all bracketed error margins correspond to the 95% credible interval (CrI) for reported estimates.

Role of modelling in COVID-19 policy development.

Models have played an important role in policy development to address the COVID-19 outbreak from its emergence in China to the current global pandemic. Early projections of international spread influenced travel restrictions and border closures. Model projections based on the virus's infectiousness demonstrated its pandemic potential, which guided the global response to and prepared countries for increases in hospitalisations and deaths. Tracking the impact of distancing and movement policies and behaviour changes has been critical in evaluating these decisions. Models have provided insights into the epidemiological differences between higher and lower income countries, as well as vulnerable population groups within countries to help design fit-for-purpose policies. Economic evaluation and policies have combined epidemic models and traditional economic models to address the economic consequences of COVID-19, which have informed policy calls for easing restrictions. Social contact and mobility models have allowed evaluation of the pathways to safely relax mobility restrictions and distancing measures. Finally, models can consider future end-game scenarios, including how suppression can be achieved and the impact of different vaccination strategies.