To Örebro University

The sexually transmitted pathogen Neisseria gonorrhoeae is regarded as being on the way to becoming an untreatable superbug. Despite its clinical importance, little is known about its emergence and evolution, and how this corresponds with the introduction of antimicrobials. We present a genome-based phylogeographical analysis of 419 gonococcal isolates from across the globe. Results indicate that modern gonococci originated in Europe or Africa, possibly as late as the sixteenth century and subsequently disseminated globally. We provide evidence that the modern gonococcal population has been shaped by antimicrobial treatment of sexually transmitted infections as well as other infections, leading to the emergence of two major lineages with different evolutionary strategies. The well-described multidrug-resistant lineage is associated with high rates of homologous recombination and infection in high-risk sexual networks. A second, multisusceptible lineage is more associated with heterosexual networks, with potential implications for infection control.

BACKGROUND: Multidrug-resistant Neisseria gonorrhoeae strains are prevalent, threatening gonorrhoea treatment globally, and understanding of emergence, evolution, and spread of antimicrobial resistance (AMR) in gonococci remains limited. We describe the genomic evolution of gonococci and their AMR, related to the introduction of antimicrobial therapies, examining isolates from 1928 (preantibiotic era) to 2013 in Denmark. This is, to our knowledge, the oldest gonococcal collection globally.

METHODS: Lyophilised isolates were revived and examined using Etest (18 antimicrobials) and whole-genome sequencing (WGS). Quality-assured genome sequences were obtained for 191 viable and 40 non-viable isolates and analysed with multiple phylogenomic approaches.

RESULTS: Gonococcal AMR, including an accumulation of multiple AMR determinants, started to emerge particularly in the 1950s-1970s. By the twenty-first century, resistance to most antimicrobials was common. Despite that some AMR determinants affect many physiological functions and fitness, AMR determinants were mainly selected by the use/misuse of gonorrhoea therapeutic antimicrobials. Most AMR developed in strains belonging to one multidrug-resistant (MDR) clade with close to three times higher genomic mutation rate. Modern N. gonorrhoeae was inferred to have emerged in the late-1500s and its genome became increasingly conserved over time.

CONCLUSIONS: WGS of gonococci from 1928 to 2013 showed that no AMR determinants, except penB, were in detectable frequency before the introduction of gonorrhoea therapeutic antimicrobials. The modern gonococcus is substantially younger than previously hypothesized and has been evolving into a more clonal species, driven by the use/misuse of antimicrobials. The MDR gonococcal clade should be further investigated for early detection of strains with predispositions to develop and maintain MDR and for initiation of public health interventions.

BACKGROUND: Regular quality-assured whole-genome sequencing linked to antimicrobial resistance (AMR) and patient metadata is imperative to elucidate the shifting gonorrhoea epidemiology, both nationally and internationally. We aimed to examine the gonococcal population in the European Economic Area (EEA) in 2020, elucidate emerging and disappearing gonococcal lineages associated with AMR and patient metadata, compare with 2013 and 2018 whole-genome sequencing data, and explain changes in gonococcal AMR and gonorrhoea epidemiology.

METHODS: In this retrospective genomic surveillance study, we analysed consecutive gonococcal isolates that were collected in EEA countries through the European Gonococcal Antimicrobial Surveillance Programme (Euro-GASP) in 2020, and made comparisons with Euro-GASP data from 2013 and 2018. All isolates had linked AMR data (based on minimum inhibitory concentration determination) and patient metadata. We performed whole-genome sequencing and molecular typing and AMR determinants were derived from quality-checked whole-genome sequencing data. Links between genomic lineages, AMR, and patient metadata were examined.

FINDINGS: 1932 gonococcal isolates collected in 2020 in 21 EEA countries were included. The majority (81·2%, 147 of 181 isolates) of azithromycin resistance (present in 9·4%, 181 of 1932) was explained by the continued expansion of the Neisseria gonorrhoeae sequence typing for antimicrobial resistance (NG-STAR) clonal complexes (CCs) 63, 168, and 213 (with mtrD/mtrR promoter mosaic 2) and the novel NG-STAR CC1031 (semi-mosaic mtrD variant 13), associated with men who have sex with men and anorectal or oropharyngeal infections. The declining cefixime resistance (0·5%, nine of 1932) and negligible ceftriaxone resistance (0·1%, one of 1932) was largely because of the progressive disappearance of NG-STAR CC90 (with mosaic penA allele), which was predominant in 2013. No known resistance determinants for novel antimicrobials (zoliflodacin, gepotidacin, and lefamulin) were found.

INTERPRETATION: Azithromycin-resistant clones, mainly with mtrD mosaic or semi-mosaic variants, appear to be stabilising at a relatively high level in the EEA. This mostly low-level azithromycin resistance might threaten the recommended ceftriaxone-azithromycin therapy, but the negligible ceftriaxone resistance is encouraging. The decreased genomic population diversity and increased clonality could be explained in part by the COVID-19 pandemic resulting in lower importation of novel strains into Europe.