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.

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.

On the probability of strain invasion in endemic settings: Accounting for individual heterogeneity and control in multi-strain dynamics.

Pathogen evolution is an imminent threat to global health that has warranted, and duly received, considerable attention within the medical, microbiological and modelling communities. Outbreaks of new pathogens are often ignited by the emergence and transmission of mutant variants descended from wild-type strains circulating in the community. In this work we investigate the stochastic dynamics of the emergence of a novel disease strain, introduced into a population in which it must compete with an existing endemic strain. In analogy with past work on single-strain epidemic outbreaks, we apply a branching process approximation to calculate the probability that the new strain becomes established. As expected, a critical determinant of the survival prospects of any invading strain is the magnitude of its reproduction number relative to that of the background endemic strain. Whilst in most circumstances this ratio must exceed unity in order for invasion to be viable, we show that differential control scenarios can lead to less-fit novel strains invading populations hosting a fitter endemic one. This analysis and the accompanying findings will inform our understanding of the mechanisms that have led to past instances of successful strain invasion, and provide valuable lessons for thwarting future drug-resistant strain incursions.

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.

Cost-effectiveness of 3 months of weekly rifapentine and isoniazid compared with other standard treatment regimens for latent tuberculosis infection: a decision analysis study.

Background: Latent tuberculosis infection (LTBI) is a critical driver of the global burden of active TB, and therefore LTBI treatment is key for TB elimination. Treatment regimens for LTBI include self-administered daily isoniazid for 6 (6H) or 9 (9H) months, self-administered daily rifampicin plus isoniazid for 3 months (3RH), self-administered daily rifampicin for 4 months (4R) and weekly rifapentine plus isoniazid for 3 months self-administered (3HP-SAT) or administered by a healthcare worker as directly observed therapy (3HP-DOT). Data on the relative cost-effectiveness of these regimens are needed to assist policymakers and clinicians in selecting an LTBI regimen. Objectives: To evaluate the cost-effectiveness of all regimens for treating LTBI. Methods: We developed a Markov model to investigate the cost-effectiveness of 3HP-DOT, 3HP-SAT, 4R, 3RH, 9H and 6H for LTBI treatment in a cohort of 10000 adults with LTBI. Cost-effectiveness was evaluated from a health system perspective over a 20 year time horizon. Results: Compared with no preventive treatment, 3HP-DOT, 3HP-SAT, 4R, 3RH, 9H and 6H prevented 496, 470, 442, 418, 370 and 276 additional cases of active TB per 10000 patients, respectively. All regimens reduced costs and increased QALYs compared with no preventive treatment. 3HP was more cost-effective under DOT than under SAT at a cost of US$27948 per QALY gained. Conclusions: Three months of weekly rifapentine plus isoniazid is more cost-effective than other regimens. Greater recognition of the benefits of short-course regimens can contribute to the scale-up of prevention and achieving the 'End TB' targets.

Global stability properties of a class of renewal epidemic models.

We investigate the global dynamics of a general Kermack-McKendrick-type epidemic model formulated in terms of a system of renewal equations. Specifically, we consider a renewal model for which both the force of infection and the infected removal rates are arbitrary functions of the infection age, [Formula: see text], and use the direct Lyapunov method to establish the global asymptotic stability of the equilibrium solutions. In particular, we show that the basic reproduction number, [Formula: see text], represents a sharp threshold parameter such that for [Formula: see text], the infection-free equilibrium is globally asymptotically stable; whereas the endemic equilibrium becomes globally asymptotically stable when [Formula: see text], i.e. when it exists.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.