Dynamics of Neisseria gonorrhoeae transmission among female sex workers and clients: A mathematical modeling study

Abstract

BackgroundThis study aimed to examine the transmission dynamics of Neisseria gonorrhoeae (NG) in heterosexual sex work networks (HSWNs) and the impact of variation in sexual behavior and interventions on NG epidemiology. MethodsThe study employed an individual-based mathematical model to simulate NG transmission dynamics in sexual networks involving female sex workers (FSWs) and their clients, primarily focusing on the Middle East and North Africa region. A deterministic model was also used to describe NG transmission from clients to their spouses. ResultsNG epidemiology in HSWNs displays two distinct patterns. In the common low-partner-number HSWNs, a significant proportion of NG incidence occurs among FSWs, with NG prevalence 13 times higher among FSWs than clients, and three times higher among clients than their spouses. Interventions substantially reduce incidence. Increasing condom use from 10 % to 50 % lowers NG prevalence among FSWs, clients, and their spouses from 12.2 % to 6.4 %, 1.2 % to 0.5 %, and 0.4 % to 0.2 %, respectively. Increasing symptomatic treatment coverage among FSWs from 0 % to 100 % decreases prevalence from 10.6 % to 4.5 %, 0.8 % to 0.4 %, and 0.3 % to 0.1 %, respectively. Increasing asymptomatic treatment coverage among FSWs from 0 % to 50 % decreases prevalence from 8.2 % to 0.4 %, 0.6 % to 0.1 %, and 0.2 % to 0.0 %, respectively, with very low prevalence when coverage exceeds 50 %. In high-partner-number HSWNs, prevalence among FSWs saturates at a high level, and the vast majority of incidence occurs among clients and their spouses, with a limited impact of incremental increases in interventions. ConclusionNG epidemiology in HSWNs is typically a "fragile epidemiology" that is responsive to a range of interventions even if the interventions are incremental, partially efficacious, and only applied to FSWs.