A dynamic power-law sexual network model of gonorrhoea outbreaks

Review of mathematical models of Neisseria gonorrhoeae vaccine impact: Implications for vaccine development

An effective prophylactic vaccine for prevention of Neisseria gonorrhoeae infection would have a major impact on sexual and reproductive health worldwide. Interest in developing gonorrhoea vaccines is growing due to the reported high rates of N. gonorrhoeae infections globally, and the threat of antimicrobial resistance. Several gonorrhoea vaccine candidates are currently under evaluation and various mathematical models have been used to assess the potential population-level impact a gonorrhoea vaccine may have once available. Here we review key aspects of gonorrhoea vaccine mathematical modelling studies, including model structures, populations considered, and assumptions used as well as vaccine characteristics and implementation scenarios investigated. The predicted vaccine impact varied between studies, ranging from as little as ∼17 % reduction in N. gonorrhoeae prevalence after 30 years up to 100 % reduction after 5 years. However, all studies predicted that even a partially effective gonorrhoea vaccine could have a substantial impact in reducing N. gonorrhoeae prevalence or incidence, particularly when high coverage is achieved within either important risk groups or the overall sexually active population. As expected, higher vaccine efficacy against acquisition of N. gonorrhoeae and longer duration of protection were linked to greater reductions in infections. A vaccine that alters onward transmission could also substantially reduce infections. Several gaps and research needs have been identified by researchers in the field and via this narrative literature review. For example, future modelling to inform gonorrhoea vaccine development and implementation should consider additional populations that are at high risk of N. gonorrhoeae infection, especially in low- and middle-income settings, as well as the impact of vaccination on the potential adverse sexual and reproductive health outcomes of infection. In addition, more detailed and robust epidemiological, biological, and behavioural data is needed to enable more accurate and robust modelling of gonorrhoea vaccine impact to inform future scientific and public health decision-making.

Underdetected dispersal and extensive local transmission drove the 2022 mpox epidemic

The World Health Organization declared mpox a public health emergency of international concern in July 2022. To investigate global mpox transmission and population-level changes associated with controlling spread, we built phylogeographic and phylodynamic models to analyze MPXV genomes from five global regions together with air traffic and epidemiological data. Our models reveal community transmission prior to detection, changes in case reporting throughout the epidemic, and a large degree of transmission heterogeneity. We find that viral introductions played a limited role in prolonging spread after initial dissemination, suggesting that travel bans would have had only a minor impact. We find that mpox transmission in North America began declining before more than 10% of high-risk individuals in the USA had vaccine-induced immunity. Our findings highlight the importance of broader routine specimen screening surveillance for emerging infectious diseases and of joint integration of genomic and epidemiological information for early outbreak control.

Different time scales used for sexual partner surveys pose a challenge in modelling dynamics of sexually transmitted infections

Mathematical models for sexually transmitted infections (STIs) are parameterised by empirical data on sexual behaviour (e.g. the number of partners over a given period) obtained from surveys. However, the time window for reporting sexual partnerships may vary between surveys and how data for different windows can be translated from one to another remains an open question. To highlight this issue, we compared the distributions of the number of sexual partners over one year and four weeks from the British National Surveys of Sexual Attitudes and Lifestyles. The results show that simple linear rescaling did not render the one-year and four-week partner distributions aligned. Parameterising STI models using survey-based sexual encounter rates without considering the implication of the reporting window used can lead to misleading results.

Early underdetected dissemination across countries followed by extensive local transmission propelled the 2022 mpox epidemic

The World Health Organization declared mpox a public health emergency of international concern in July 2022. To investigate global mpox transmission and population-level changes associated with controlling spread, we built phylogeographic and phylodynamic models to analyze MPXV genomes from five global regions together with air traffic and epidemiological data. Our models reveal community transmission prior to detection, changes in case-reporting throughout the epidemic, and a large degree of transmission heterogeneity. We find that viral introductions played a limited role in prolonging spread after initial dissemination, suggesting that travel bans would have had only a minor impact. We find that mpox transmission in North America began declining before more than 10% of high-risk individuals in the USA had vaccine-induced immunity. Our findings highlight the importance of broader routine specimen screening surveillance for emerging infectious diseases and of joint integration of genomic and epidemiological information for early outbreak control.

The role of vaccination and public awareness in forecasts of Mpox incidence in the United Kingdom

Beginning in May 2022, Mpox virus spread rapidly in high-income countries through close human-to-human contact primarily amongst communities of gay, bisexual and men who have sex with men (GBMSM). Behavioural change arising from increased knowledge and health warnings may have reduced the rate of transmission and modified Vaccinia-based vaccination is likely to be an effective longer-term intervention. We investigate the UK epidemic presenting 26-week projections using a stochastic discrete-population transmission model which includes GBMSM status, rate of formation of new sexual partnerships, and clique partitioning of the population. The Mpox cases peaked in mid-July; our analysis is that the decline was due to decreased transmission rate per infected individual and infection-induced immunity among GBMSM, especially those with the highest rate of new partners. Vaccination did not cause Mpox incidence to turn over, however, we predict that a rebound in cases due to behaviour reversion was prevented by high-risk group-targeted vaccination.

Estimating the fitness cost and benefit of antimicrobial resistance from pathogen genomic data

Increasing levels of antibiotic resistance in many bacterial pathogen populations are a major threat to public health. Resistance to an antibiotic provides a fitness benefit when the bacteria are exposed to this antibiotic, but resistance also often comes at a cost to the resistant pathogen relative to susceptible counterparts. We lack a good understanding of these benefits and costs of resistance for many bacterial pathogens and antibiotics, but estimating them could lead to better use of antibiotics in a way that reduces or prevents the spread of resistance. Here, we propose a new model for the joint epidemiology of susceptible and resistant variants, which includes explicit parameters for the cost and benefit of resistance. We show how Bayesian inference can be performed under this model using phylogenetic data from susceptible and resistant lineages and that by combining data from both we are able to disentangle and estimate the resistance cost and benefit parameters separately. We applied our inferential methodology to several simulated datasets to demonstrate good scalability and accuracy. We analysed a dataset of Neisseria gonorrhoeae genomes collected between 2000 and 2013 in the USA. We found that two unrelated lineages resistant to fluoroquinolones shared similar epidemic dynamics and resistance parameters. Fluoroquinolones were abandoned for the treatment of gonorrhoea due to increasing levels of resistance, but our results suggest that they could be used to treat a minority of around 10% of cases without causing resistance to grow again.

Power Law for Estimating Underdetection of Foodborne Disease Outbreaks, United States

We fit a power law distribution to US foodborne disease outbreaks to assess underdetection and underreporting. We predicted that 788 fewer than expected small outbreaks were identified annually during 1998-2017 and 365 fewer during 2018-2019, after whole-genome sequencing was implemented. Power law can help assess effectiveness of public health interventions.

Heavy-tailed sexual contact networks and monkeypox epidemiology in the global outbreak, 2022

The outbreak of monkeypox across non-endemic regions confirmed in May 2022 shows epidemiological features distinct from previously imported outbreaks, most notably its observed growth and predominance amongst men who have sex with men (MSM). We use a transmission model fitted to empirical sexual partnership data to show that the heavy-tailed sexual partnership distribution, in which a handful of individuals have disproportionately many partners, can explain the sustained growth of monkeypox among MSM despite the absence of such patterns previously. We suggest that the basic reproduction number (R0) for monkeypox over the MSM sexual network may be substantially above 1, which poses challenges to outbreak containment. Ensuring support and tailored messaging to facilitate prevention and early detection among MSM with high numbers of partners is warranted.

Non-selective distribution of infectious disease prevention may outperform risk-based targeting

Epidemic control often requires optimal distribution of available vaccines and prophylactic tools, to protect from infection those susceptible. Well-established theory recommends prioritizing those at the highest risk of exposure. But the risk is hard to estimate, especially for diseases involving stigma and marginalization. We address this conundrum by proving that one should target those at high risk only if the infection-averting efficacy of prevention is above a critical value, which we derive analytically. We apply this to the distribution of pre-exposure prophylaxis (PrEP) of the Human Immunodeficiency Virus (HIV) among men-having-sex-with-men (MSM), a population particularly vulnerable to HIV. PrEP is effective in averting infections, but its global scale-up has been slow, showing the need to revisit distribution strategies, currently risk-based. Using data from MSM communities in 58 countries, we find that non-selective PrEP distribution often outperforms risk-based, showing that a logistically simpler strategy is also more effective. Our theory may help design more feasible and successful prevention.

Pre-exposure prophylaxis (PrEP) is an effective HIV prevention measure but identifying those most at risk to target for treatment is challenging. Here, the authors demonstrate that non-selective PrEP distribution outperforms targeted strategies when use is not consistent, and/or prevalence of untreated HIV is high.

Global emergence and dissemination of Neisseria gonorrhoeae ST-9363 isolates with reduced susceptibility to azithromycin

Neisseria gonorrhoeae multi-locus sequence type (ST) 9363 core-genogroup isolates have been associated with reduced azithromycin susceptibility (AZMrs) and show evidence of clonal expansion in the U.S. Here we analyze a global collection of ST-9363 core-genogroup genomes to shed light on the emergence and dissemination of this strain. The global population structure of ST-9363 core-genogroup falls into three lineages: Basal, European, and North American; with 32 clades within all lineages. Although, ST-9363 core-genogroup is inferred to have originated from Asia in the mid-19th century; we estimate the three modern lineages emerged from Europe in the late 1970s to early 1980s. The European lineage appears to have emerged and expanded from around 1986 to 1998, spreading into North America and Oceania in the mid-2000s with multiple introductions, along with multiple secondary reintroductions into Europe. Our results suggest two separate acquisition events of mosaic mtrR and mtrR promoter alleles: first during 2009-2011 and again during the 2012-2013 time, facilitating the clonal expansion of this core-genogroup with AZMrs in the U.S. By tracking phylodynamic evolutionary trajectories of clades that share distinct demography as well as population-based genomic statistics, we demonstrate how recombination and selective pressures in the mtrCDE efflux operon granted a fitness advantage to establish ST-9363 as a successful gonococcal lineage in the U.S. and elsewhere. Although it is difficult to pinpoint the exact timing and emergence of this young core-genogroup, it remains critically important to continue monitoring it, as it could acquire additional resistance markers.

Genomic Epidemiology Analysis of Infectious Disease Outbreaks Using TransPhylo

Comparing the pathogen genomes from several cases of an infectious disease has the potential to help us understand and control outbreaks. Many methods exist to reconstruct a phylogeny from such genomes, which represents how the genomes are related to one another. However, such a phylogeny is not directly informative about transmission events between individuals. TransPhylo is a software tool implemented as an R package designed to bridge the gap between pathogen phylogenies and transmission trees. TransPhylo is based on a combined model of transmission between hosts and pathogen evolution within each host. It can simulate both phylogenies and transmission trees jointly under this combined model. TransPhylo can also reconstruct a transmission tree based on a dated phylogeny, by exploring the space of transmission trees compatible with the phylogeny. A transmission tree can be represented as a coloring of a phylogeny where each color represents a different host of the pathogen, and TransPhylo provides convenient ways to plot these colorings and explore the results. This article presents the basic protocols that can be used to make the most of TransPhylo. © 2021 The Authors. Basic Protocol 1: First steps with TransPhylo Basic Protocol 2: Simulation of outbreak data Basic Protocol 3: Inference of transmission Basic Protocol 4: Exploring the results of inference.

Sudden emergence of a Neisseria gonorrhoeae clade with reduced susceptibility to extended-spectrum cephalosporins, Norway

Neisseria gonorrhoeae multilocus sequence type (ST)-7827 emerged in a dramatic fashion in Norway in the period 2016-2018. Here, we aim to shed light on the provenance and expansion of this ST. ST-7827 was found to be polyphyletic, but the majority of members belonged to a monophyletic clade we termed PopPUNK cluster 7827 (PC-7827). In Norway, both PC-7827 and ST-7827 isolates were almost exclusively isolated from men. Phylogeographical analyses demonstrated an Asian origin of the genogroup, with multiple inferred exports to Europe and the USA. The genogroup was uniformly resistant to fluoroquinolones, and associated with reduced susceptibility to both azithromycin and the extended-spectrum cephalosporins (ESCs) cefixime and ceftriaxone. From a genetic background including the penA allele 13.001, associated with reduced ESC susceptibility, we identified repeated events of acquisition of porB alleles associated with further reduction in ceftriaxone susceptibility. Transmission of the strain was significantly reduced in Norway in 2019, but our results indicate the existence of a recently established global reservoir. The worrisome drug-resistance profile and rapid emergence of PC-7827 calls for close monitoring of the situation.

Assessment of the Potential of Vaccination to Combat Antibiotic Resistance in Gonorrhea: A Modeling Analysis to Determine Preferred Product Characteristics

BACKGROUND

Gonorrhea incidence is increasing rapidly in many countries, while antibiotic resistance is making treatment more difficult. Combined with evidence that two meningococcal vaccines are likely partially protective against gonorrhea, this has renewed interest in a gonococcal vaccine, and several candidates are in development. Key questions are how protective and long-lasting a vaccine needs to be, and how to target it. We assessed vaccination's potential impact and the feasibility of achieving the World Health Organization's (WHO) target of reducing gonorrhea incidence by 90% during 2018-2030, by comparing realistic vaccination strategies under a range of scenarios of vaccine efficacy and duration of protection, and emergence of extensively-resistant gonorrhea.

METHODS

We developed a stochastic transmission-dynamic model, incorporating asymptomatic and symptomatic infection and heterogeneous sexual behavior in men who have sex with men (MSM). We used data from England, which has a comprehensive, consistent nationwide surveillance system. Using particle Markov chain Monte Carlo methods, we fitted to gonorrhea incidence in 2008-2017, then used Bayesian forecasting to examine an extensive range of scenarios.

RESULTS

Even in the worst-case scenario of untreatable infection emerging, the WHO target is achievable if all MSM attending sexual health clinics receive a vaccine offering ≥ 52% protection for ≥ 6 years. A vaccine conferring 31% protection (as estimated for MeNZB) for 2-4 years could reduce incidence in 2030 by 45% in the worst-case scenario, and by 75% if > 70% of resistant gonorrhea remains treatable.

CONCLUSIONS

Even a partially-protective vaccine, delivered through a realistic targeting strategy, could substantially reduce gonorrhea incidence, despite antibiotic resistance.

Improving pairwise approximations for network models with susceptible-infected-susceptible dynamics

Network models of disease spread play an important role in elucidating the impact of long-lasting infectious contacts on the dynamics of epidemics. Moment-closure approximation is a common method of generating low-dimensional deterministic models of epidemics on networks, which has found particular success for diseases with susceptible-infected-recovered (SIR) dynamics. However, the effect of network structure is arguably more important for sexually transmitted infections, where epidemiologically relevant contacts are comparatively rare and longstanding, and which are in general modelled via the susceptible-infected-susceptible (SIS)-paradigm. In this paper, we introduce an improvement to the standard pairwise approximation for network models with SIS-dynamics for two different network structures: the isolated open triple (three connected individuals in a line) and the k-regular network. This improvement is achieved by tracking the rate of change of errors between triple values and their standard pairwise approximation. For the isolated open triple, this improved pairwise model is exact, while for k-regular networks a closure is made at the level of triples to obtain a closed set of equations. This improved pairwise approximation provides an insight into the errors introduced by the standard pairwise approximation, and more closely matches both higher-order moment-closure approximations and explicit stochastic simulations with only a modest increase in dimensionality to the standard pairwise approximation.

Agent-based modelling study of antimicrobial-resistant Neisseria gonorrhoeae transmission in men who have sex with men: towards individualised diagnosis and treatment

Background Antimicrobial-resistant (AMR) gonorrhoea is a global public health threat. Discriminatory point-of-care tests (POCT) to detect drug sensitivity are under development, enabling individualised resistance-guided therapy.

METHODS

An individual-based dynamic transmission model of gonorrhoea infection in MSM living in London has been developed, incorporating ciprofloxacin-sensitive and resistant strains. The time-dependent sexual contact network is captured by periodically restructuring active connections to reflect the transience of contacts. Different strategies to improve treatment selection were explored, including discriminatory POCT and selecting partner treatment based on either the index case or partner susceptibility. Outcomes included population prevalence of gonorrhoea and drug dose counts.

RESULTS

It is shown that using POCT to detect ciprofloxacin-sensitive infections could result in a large decrease in ceftriaxone doses (by 70% compared with the reference case in the simulations of this study). It also suggests that ceftriaxone use can be reduced with existing technologies, albeit to a lesser degree; either using index case sensitivity profiles to direct treatment of partners, or testing notified partners with strain discriminatory laboratory tests before treatment, reduced ceftriaxone use in our model (by 27% and 47% respectively).

CONCLUSIONS

POCT to detect ciprofloxacin-sensitive gonorrhoea are likely to dramatically reduce reliance on ceftriaxone, but requires the implementation of new technology. In the meantime, the proportion of unnecessary ceftriaxone treatment by testing partners before treatment could be reduced significantly. Alternatively, index case sensitivity profiles could be used to select effective treatments for partners.