Antimicrobial-resistant Neisseria gonorrhoeae in Europe in 2020 compared with in 2013 and 2018: a retrospective genomic surveillance study.

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. FUNDING: European Centre for Disease Prevention and Control and Örebro University Hospital.

pyngoST: fast, simultaneous and accurate multiple sequence typing of Neisseria gonorrhoeae genome collections.

Extensive gonococcal surveillance has been performed using molecular typing at global, regional, national and local levels. The three main genotyping schemes for this pathogen, multi-locus sequence typing (MLST), Neisseria gonorrhoeae multi-antigen sequence typing (NG-MAST) and N. gonorrhoeae sequence typing for antimicrobial resistance (NG-STAR), allow inter-laboratory and inter-study comparability and reproducibility and provide an approximation to the gonococcal population structure. With whole-genome sequencing (WGS), we obtain a substantially higher and more accurate discrimination between strains compared to previous molecular typing schemes. However, WGS remains unavailable or not affordable in many laboratories, and thus bioinformatic tools that allow the integration of data among laboratories with and without access to WGS are imperative for a joint effort to increase our understanding of global pathogen threats. Here, we present pyngoST, a command-line Python tool for fast, simultaneous and accurate sequence typing of N. gonorrhoeae from WGS assemblies. pyngoST integrates MLST, NG-MAST and NG-STAR, and can also designate NG-STAR clonal complexes, NG-MAST genogroups and penA mosaicism, facilitating multiple sequence typing from large WGS assembly collections. Exact and closest matches for existing alleles and sequence types are reported. The implementation of a fast multi-pattern searching algorithm allows pyngoST to be rapid and report results on 500 WGS assemblies in under 1 min. The mapping of typing results on a core genome tree of 2375 gonococcal genomes revealed that NG-STAR is the scheme that best represents the population structure of this pathogen, emphasizing the role of antimicrobial use and antimicrobial resistance as a driver of gonococcal evolution. This article contains data hosted by Microreact.

Genomics for public health and international surveillance of antimicrobial resistance.

Historically, epidemiological investigation and surveillance for bacterial antimicrobial resistance (AMR) has relied on low-resolution isolate-based phenotypic analyses undertaken at local and national reference laboratories. Genomic sequencing has the potential to provide a far more high-resolution picture of AMR evolution and transmission, and is already beginning to revolutionise how public health surveillance networks monitor and tackle bacterial AMR. However, the routine integration of genomics in surveillance pipelines still has considerable barriers to overcome. In 2022, a workshop series and online consultation brought together international experts in AMR and pathogen genomics to assess the status of genomic applications for AMR surveillance in a range of settings. Here we focus on discussions around the use of genomics for public health and international AMR surveillance, noting the potential advantages of, and barriers to, implementation, and proposing recommendations from the working group to help to drive the adoption of genomics in public health AMR surveillance. These recommendations include the need to build capacity for genome sequencing and analysis, harmonising and standardising surveillance systems, developing equitable data sharing and governance frameworks, and strengthening interactions and relationships among stakeholders at multiple levels.

GyrB in silico mining in 27 151 global gonococcal genomes from 1928-2021 combined with zoliflodacin in vitro testing of 71 international gonococcal isolates with different GyrB, ParC and ParE substitutions confirms high susceptibility.

OBJECTIVES: Antimicrobial resistance (AMR) in Neisseria gonorrhoeae is a global threat and novel treatment alternatives are imperative. Herein, susceptibility to the novel antimicrobial zoliflodacin, currently in a global Phase 3 randomized controlled clinical trial for gonorrhoea treatment, was investigated by screening for zoliflodacin GyrB target mutations in publicly available gonococcal genomes and, where feasible, determination of the associated zoliflodacin MIC. METHODS: The European Nucleotide Archive was queried using the search term 'Taxon: 485'. DNA sequences from 27 151 gonococcal isolates were analysed and gyrB, gyrA, parC and parE alleles characterized. RESULTS: GyrB amino acid alterations were rare (97.0% of isolates had a wild-type GyrB sequence). GyrB V470L (2.7% of isolates) was the most prevalent alteration, followed by S467N (0.12%), N. meningitidis GyrB (0.092%), V470I (0.059%), Q468R/P (0.015%), A466T (0.0074%), L425I + L465I (0.0037%), L465I (0.0037%), G482S (0.0037%) and D429V (0.0037%). Only one isolate (0.0037%) carried a substitution in a resistance-associated GyrB codon (D429V), resulting in a zoliflodacin MIC of 8 mg/L. None of the other detected gyrB, gyrA, parC or parE mutations caused a zoliflodacin MIC outside the wild-type MIC distribution. CONCLUSIONS: The zoliflodacin target GyrB was highly conserved among 27 151 global gonococcal isolates cultured in 1928-2021. The single zoliflodacin-resistant clinical isolate (0.0037%) was cultured from a male patient in Japan in 2000. Evidently, this strain has not clonally expanded nor has the gyrB zoliflodacin-resistance mutation disseminated through horizontal gene transfer to other strains. Phenotypic and genomic surveillance, including gyrB mutations, of zoliflodacin susceptibility are imperative.

Europe-wide expansion and eradication of multidrug-resistant Neisseria gonorrhoeae lineages: a genomic surveillance study.

BACKGROUND: Genomic surveillance using quality-assured whole-genome sequencing (WGS) together with epidemiological and antimicrobial resistance (AMR) data is essential to characterise the circulating Neisseria gonorrhoeae lineages and their association to patient groups (defined by demographic and epidemiological factors). In 2013, the European gonococcal population was characterised genomically for the first time. We describe the European gonococcal population in 2018 and identify emerging or vanishing lineages associated with AMR and epidemiological characteristics of patients, to elucidate recent changes in AMR and gonorrhoea epidemiology in Europe. METHODS: We did WGS on 2375 gonococcal isolates from 2018 (mainly Sept 1-Nov 30) in 26 EU and EEA countries. Molecular typing and AMR determinants were extracted from quality-checked genomic data. Association analyses identified links between genomic lineages, AMR, and epidemiological data. FINDINGS: Azithromycin-resistant N gonorrhoeae (8·0% [191/2375] in 2018) is rising in Europe due to the introduction or emergence and subsequent expansion of a novel N gonorrhoeae multi-antigen sequence typing (NG-MAST) genogroup, G12302 (132 [5·6%] of 2375; N gonorrhoeae sequence typing for antimicrobial resistance [NG-STAR] clonal complex [CC]168/63), carrying a mosaic mtrR promoter and mtrD sequence and found in 24 countries in 2018. CC63 was associated with pharyngeal infections in men who have sex with men. Susceptibility to ceftriaxone and cefixime is increasing, as the resistance-associated lineage, NG-MAST G1407 (51 [2·1%] of 2375), is progressively vanishing since 2009-10. INTERPRETATION: Enhanced gonococcal AMR surveillance is imperative worldwide. WGS, linked to epidemiological and AMR data, is essential to elucidate the dynamics in gonorrhoea epidemiology and gonococcal populations as well as to predict AMR. When feasible, WGS should supplement the national and international AMR surveillance programmes to elucidate AMR changes over time. In the EU and EEA, increasing low-level azithromycin resistance could threaten the recommended ceftriaxone-azithromycin dual therapy, and an evidence-based clinical azithromycin resistance breakpoint is needed. Nevertheless, increasing ceftriaxone susceptibility, declining cefixime resistance, and absence of known resistance mutations for new treatments (zoliflodacin, gepotidacin) are promising. FUNDING: European Centre for Disease Prevention and Control, Centre for Genomic Pathogen Surveillance, Örebro University Hospital, Wellcome.

Evolutionary Processes in the Emergence and Recent Spread of the Syphilis Agent, Treponema pallidum.

The incidence of syphilis has risen worldwide in the last decade in spite of being an easily treated infection. The causative agent of this sexually transmitted disease is the bacterium Treponema pallidum subspecies pallidum (TPA), very closely related to subsp. pertenue (TPE) and endemicum (TEN), responsible for the human treponematoses yaws and bejel, respectively. Although much focus has been placed on the question of the spatial and temporary origins of TPA, the processes driving the evolution and epidemiological spread of TPA since its divergence from TPE and TEN are not well understood. Here, we investigate the effects of recombination and selection as forces of genetic diversity and differentiation acting during the evolution of T. pallidum subspecies. Using a custom-tailored procedure, named phylogenetic incongruence method, with 75 complete genome sequences, we found strong evidence for recombination among the T. pallidum subspecies, involving 12 genes and 21 events. In most cases, only one recombination event per gene was detected and all but one event corresponded to intersubspecies transfers, from TPE/TEN to TPA. We found a clear signal of natural selection acting on the recombinant genes, which is more intense in their recombinant regions. The phylogenetic location of the recombination events detected and the functional role of the genes with signals of positive selection suggest that these evolutionary processes had a key role in the evolution and recent expansion of the syphilis bacteria and significant implications for the selection of vaccine candidates and the design of a broadly protective syphilis vaccine.

A global resource for genomic predictions of antimicrobial resistance and surveillance of Salmonella Typhi at pathogenwatch.

As whole-genome sequencing capacity becomes increasingly decentralized, there is a growing opportunity for collaboration and the sharing of surveillance data within and between countries to inform typhoid control policies. This vision requires free, community-driven tools that facilitate access to genomic data for public health on a global scale. Here we present the Pathogenwatch scheme for Salmonella enterica serovar Typhi (S. Typhi), a web application enabling the rapid identification of genomic markers of antimicrobial resistance (AMR) and contextualization with public genomic data. We show that the clustering of S. Typhi genomes in Pathogenwatch is comparable to established bioinformatics methods, and that genomic predictions of AMR are highly concordant with phenotypic susceptibility data. We demonstrate the public health utility of Pathogenwatch with examples selected from >4,300 public genomes available in the application. Pathogenwatch provides an intuitive entry point to monitor of the emergence and spread of S. Typhi high risk clones.

A community-driven resource for genomic epidemiology and antimicrobial resistance prediction of Neisseria gonorrhoeae at Pathogenwatch.

BACKGROUND: Antimicrobial-resistant (AMR) Neisseria gonorrhoeae is an urgent threat to public health, as strains resistant to at least one of the two last-line antibiotics used in empiric therapy of gonorrhoea, ceftriaxone and azithromycin, have spread internationally. Whole genome sequencing (WGS) data can be used to identify new AMR clones and transmission networks and inform the development of point-of-care tests for antimicrobial susceptibility, novel antimicrobials and vaccines. Community-driven tools that provide an easy access to and analysis of genomic and epidemiological data is the way forward for public health surveillance. METHODS: Here we present a public health-focussed scheme for genomic epidemiology of N. gonorrhoeae at Pathogenwatch ( https://pathogen.watch/ngonorrhoeae ). An international advisory group of experts in epidemiology, public health, genetics and genomics of N. gonorrhoeae was convened to inform on the utility of current and future analytics in the platform. We implement backwards compatibility with MLST, NG-MAST and NG-STAR typing schemes as well as an exhaustive library of genetic AMR determinants linked to a genotypic prediction of resistance to eight antibiotics. A collection of over 12,000 N. gonorrhoeae genome sequences from public archives has been quality-checked, assembled and made public together with available metadata for contextualization. RESULTS: AMR prediction from genome data revealed specificity values over 99% for azithromycin, ciprofloxacin and ceftriaxone and sensitivity values around 99% for benzylpenicillin and tetracycline. A case study using the Pathogenwatch collection of N. gonorrhoeae public genomes showed the global expansion of an azithromycin-resistant lineage carrying a mosaic mtr over at least the last 10 years, emphasising the power of Pathogenwatch to explore and evaluate genomic epidemiology questions of public health concern. CONCLUSIONS: The N. gonorrhoeae scheme in Pathogenwatch provides customised bioinformatic pipelines guided by expert opinion that can be adapted to public health agencies and departments with little expertise in bioinformatics and lower-resourced settings with internet connection but limited computational infrastructure. The advisory group will assess and identify ongoing public health needs in the field of gonorrhoea, particularly regarding gonococcal AMR, in order to further enhance utility with modified or new analytic methods.

Phylogenomics of Mycobacterium africanum reveals a new lineage and a complex evolutionary history.

Human tuberculosis (TB) is caused by members of the Mycobacterium tuberculosis complex (MTBC). The MTBC comprises several human-adapted lineages known as M. tuberculosis sensu stricto, as well as two lineages (L5 and L6) traditionally referred to as Mycobacterium africanum. Strains of L5 and L6 are largely limited to West Africa for reasons unknown, and little is known of their genomic diversity, phylogeography and evolution. Here, we analysed the genomes of 350 L5 and 320 L6 strains, isolated from patients from 21 African countries, plus 5 related genomes that had not been classified into any of the known MTBC lineages. Our population genomic and phylogeographical analyses showed that the unclassified genomes belonged to a new group that we propose to name MTBC lineage 9 (L9). While the most likely ancestral distribution of L9 was predicted to be East Africa, the most likely ancestral distribution for both L5 and L6 was the Eastern part of West Africa. Moreover, we found important differences between L5 and L6 strains with respect to their phylogeographical substructure and genetic diversity. Finally, we could not confirm the previous association of drug-resistance markers with lineage and sublineages. Instead, our results indicate that the association of drug resistance with lineage is most likely driven by sample bias or geography. In conclusion, our study sheds new light onto the genomic diversity and evolutionary history of M. africanum, and highlights the need to consider the particularities of each MTBC lineage for understanding the ecology and epidemiology of TB in Africa and globally.

High susceptibility to zoliflodacin and conserved target (GyrB) for zoliflodacin among 1209 consecutive clinical Neisseria gonorrhoeae isolates from 25 European countries, 2018.

OBJECTIVES: Novel antimicrobials for treatment of gonorrhoea are imperative. The first-in-class spiropyrimidinetrione zoliflodacin is promising and currently in an international Phase 3 randomized controlled clinical trial (RCT) for treatment of uncomplicated gonorrhoea. We evaluated the in vitro activity of and the genetic conservation of the target (GyrB) and other potential zoliflodacin resistance determinants among 1209 consecutive clinical Neisseria gonorrhoeae isolates obtained from 25 EU/European Economic Area (EEA) countries in 2018 and compared the activity of zoliflodacin with that of therapeutic antimicrobials currently used. METHODS: MICs of zoliflodacin, ceftriaxone, cefixime, azithromycin and ciprofloxacin were determined using an agar dilution technique for zoliflodacin or using MIC gradient strip tests or an agar dilution technique for the other antimicrobials. Genome sequences were available for 96.1% of isolates. RESULTS: Zoliflodacin modal MIC, MIC50, MIC90 and MIC range were 0.125, 0.125, 0.125 and ≤0.004-0.5 mg/L, respectively. The resistance was 49.9%, 6.7%, 1.6% and 0.2% to ciprofloxacin, azithromycin, cefixime and ceftriaxone, respectively. Zoliflodacin did not show any cross-resistance to other tested antimicrobials. GyrB was highly conserved and no zoliflodacin gyrB resistance mutations were found. No fluoroquinolone target GyrA or ParC resistance mutations or mutations causing overexpression of the MtrCDE efflux pump substantially affected the MICs of zoliflodacin. CONCLUSIONS: The in vitro susceptibility to zoliflodacin was high and the zoliflodacin target GyrB was conserved among EU/EEA gonococcal isolates in 2018. This study supports further clinical development of zoliflodacin. However, additional zoliflodacin data regarding particularly the treatment of pharyngeal gonorrhoea, pharmacokinetics/pharmacodynamics and resistance selection, including suppression, would be valuable.

Neisseria gonorrhoeae Sequence Typing for Antimicrobial Resistance (NG-STAR) clonal complexes are consistent with genomic phylogeny and provide simple nomenclature, rapid visualization and antimicrobial resistance (AMR) lineage predictions.

OBJECTIVES: Surveillance of antimicrobial resistance (AMR) in Neisseria gonorrhoeae, supported by molecular typing, ideally through genome sequencing, is imperative. We defined N. gonorrhoeae Sequence Typing for Antimicrobial Resistance (NG-STAR) clonal complexes (CCs) and validated their usefulness in gonococcal AMR surveillance. METHODS: All NG-STAR alleles and STs available in the public database (https://ngstar.canada.ca/) were analysed using PHYLOViZ 2.0 to define CCs according to the closest founder ST with ≥5 identical alleles and founding ST with the highest number of links. The published 2013 European gonococcal dataset (n = 1054), the 2016 WHO reference strain panel (n = 14) and N. gonorrhoeae isolates with ceftriaxone resistance determinant penA-60.001 (n = 7) from several countries were used for validation. RESULTS: The majority of the isolates (n = 1063) were designated to 71 CCs. The most common CC was CC90 (n = 194), followed by CC63 (n = 166), CC139 (n = 73), CC158 (n = 73) and CC127 (n = 62). CC90 included isolates belonging to the internationally spread MDR clone N. gonorrhoeae Multi-Antigen Sequence Typing (NG-MAST) G1407 (predominantly MLST ST1901). The ceftriaxone-resistant isolates with penA-60.001 (n = 7) belonged to CC73 or STs linking between CC90 and CC73 (ST233 and ST1133). Phylogenomic analysis revealed that NG-STAR CCs more appropriately correlated to phylogenomic AMR clusters compared with MLST STs, NG-MAST STs, NG-MAST genogroups and NG-STAR STs. CONCLUSIONS: NG-STAR CCs: are consistent with the gonococcal genome phylogeny; allow rapid visualizations with limited computational requirements; provide a simple, reproducible and portable nomenclature (for WGS and conventional Sanger sequencing data); and predict AMR lineages. Phenotypic AMR surveillance, supplemented with WGS, is imperative and NG-STAR CCs can effectively support this.

Core and Accessory Genome Analysis of Vibrio mimicus.

Vibrio mimicus is an emerging pathogen, mainly associated with contaminated seafood consumption. However, little is known about its evolution, biodiversity, and pathogenic potential. This study analyzes the pan-, core, and accessory genomes of nine V. mimicus strains. The core genome yielded 2424 genes in chromosome I (ChI) and 822 genes in chromosome II (ChII), with an accessory genome comprising an average of 10.9% of the whole genome for ChI and 29% for ChII. Core genome phylogenetic trees were obtained, and V. mimicus ATCC-33654 strain was the closest to the outgroup in both chromosomes. Additionally, a phylogenetic study of eight conserved genes (ftsZ, gapA, gyrB, topA, rpoA, recA, mreB, and pyrH), including Vibrio cholerae, Vibrio parilis, Vibrio metoecus, and Vibrio caribbenthicus, clearly showed clade differentiation. The main virulence genes found in ChI corresponded with type I secretion proteins, extracellular components, flagellar proteins, and potential regulators, while, in ChII, the main categories were type-I secretion proteins, chemotaxis proteins, and antibiotic resistance proteins. The accessory genome was characterized by the presence of mobile elements and toxin encoding genes in both chromosomes. Based on the genome atlas, it was possible to characterize differential regions between strains. The pan-genome of V. mimicus encompassed 3539 genes for ChI and 2355 genes for ChII. These results give us an insight into the virulence and gene content of V. mimicus, as well as constitute the first approach to its diversity.

Increased power from conditional bacterial genome-wide association identifies macrolide resistance mutations in Neisseria gonorrhoeae.

The emergence of resistance to azithromycin complicates treatment of Neisseria gonorrhoeae, the etiologic agent of gonorrhea. Substantial azithromycin resistance remains unexplained after accounting for known resistance mutations. Bacterial genome-wide association studies (GWAS) can identify novel resistance genes but must control for genetic confounders while maintaining power. Here, we show that compared to single-locus GWAS, conducting GWAS conditioned on known resistance mutations reduces the number of false positives and identifies a G70D mutation in the RplD 50S ribosomal protein L4 as significantly associated with increased azithromycin resistance (p-value = 1.08 × 10-11). We experimentally confirm our GWAS results and demonstrate that RplD G70D and other macrolide binding site mutations are prevalent (present in 5.42% of 4850 isolates) and widespread (identified in 21/65 countries across two decades). Overall, our findings demonstrate the utility of conditional associations for improving the performance of microbial GWAS and advance our understanding of the genetic basis of macrolide resistance.

Adaptation to the cervical environment is associated with increased antibiotic susceptibility in Neisseria gonorrhoeae.

Neisseria gonorrhoeae is an urgent public health threat due to rapidly increasing incidence and antibiotic resistance. In contrast with the trend of increasing resistance, clinical isolates that have reverted to susceptibility regularly appear, prompting questions about which pressures compete with antibiotics to shape gonococcal evolution. Here, we used genome-wide association to identify loss-of-function (LOF) mutations in the efflux pump mtrCDE operon as a mechanism of increased antibiotic susceptibility and demonstrate that these mutations are overrepresented in cervical relative to urethral isolates. This enrichment holds true for LOF mutations in another efflux pump, farAB, and in urogenitally-adapted versus typical N. meningitidis, providing evidence for a model in which expression of these pumps in the female urogenital tract incurs a fitness cost for pathogenic Neisseria. Overall, our findings highlight the impact of integrating microbial population genomics with host metadata and demonstrate how host environmental pressures can lead to increased antibiotic susceptibility.

Genomic and Phenotypic Variability in Neisseria gonorrhoeae Antimicrobial Susceptibility, England.

Antimicrobial resistance (AMR) in Neisseria gonorrhoeae is a global concern. Phylogenetic analyses resolve uncertainties regarding genetic relatedness of isolates with identical phenotypes and inform whether AMR is due to new mutations and clonal expansion or separate introductions by importation. We sequenced 1,277 isolates with associated epidemiologic and antimicrobial susceptibility data collected during 2013-2016 to investigate N. gonorrhoeae genomic variability in England. Comparing genetic markers and phenotypes for AMR, we identified 2 N. gonorrhoeae lineages with different antimicrobial susceptibility profiles and 3 clusters with elevated MICs for ceftriaxone, varying mutations in the penA allele, and different epidemiologic characteristics. Our results indicate N. gonorrhoeae with reduced antimicrobial susceptibility emerged independently and multiple times in different sexual networks in England, through new mutation or recombination events and by importation. Monitoring and control for AMR in N. gonorrhoeae should cover the entire population affected, rather than focusing on specific risk groups or locations.

Genomic evolution of Neisseria gonorrhoeae since the preantibiotic era (1928-2013): antimicrobial use/misuse selects for resistance and drives evolution.

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.

Phylogenomic analysis of Neisseria gonorrhoeae transmission to assess sexual mixing and HIV transmission risk in England: a cross-sectional, observational, whole-genome sequencing study.

BACKGROUND: Characterising sexual networks with transmission of sexually transmitted infections might allow identification of individuals at increased risk of infection. We aimed to investigate sexual mixing in Neisseria gonorrhoeae transmission networks between women, heterosexual men, and men who report sex with men (MSM), and between people with and without HIV. METHODS: In this cross-sectional observational study, we whole-genome sequenced N gonorrhoeae isolates from the archive of the Gonococcal Resistance to Antimicrobials Surveillance Programme (GRASP).w Isolates that varied by five single nucleotide polymorphisms or fewer were grouped into clusters that represented sexual networks with N gonorrhoeae transmission. Clusters were described by gender, sexual risk group, and HIV status. FINDINGS: We sequenced 1277 N gonorrhoeae isolates with linked clinical and sociodemographic data that were collected in five clinics in England during 2013-16 (July 1 to Sept 30 in 2013-15; July 1 to Sept 9 in 2016). The isolates grouped into 213 clusters. 30 (14%) clusters contained isolates from heterosexual men and MSM but no women and three (1%) clusters contained isolates from only women and MSM. 146 (69%) clusters comprised solely people with negative or unknown HIV status and seven (3%) comprised only HIV-positive people. 60 (28%) clusters comprised MSM with positive and negative or unknown HIV status. INTERPRETATION: N gonorrhoeae molecular data can provide information indicating risk of HIV or other sexually transmitted infections for some individuals for whom such risk might not be known from clinical history. These findings have implications for sexual health care, including offering testing, prevention advice, and preventive treatment, such as HIV pre-exposure prophylaxis. FUNDING: National Institute for Health Research Health Protection Research Unit; Wellcome; Public Health England.

Genetic variation regulates the activation and specificity of Restriction-Modification systems in Neisseria gonorrhoeae.

Restriction-Modification systems (RMS) are one of the main mechanisms of defence against foreign DNA invasion and can have an important role in the regulation of gene expression. The obligate human pathogen Neisseria gonorrhoeae carries one of the highest loads of RMS in its genome; between 13 to 15 of the three main types. Previous work has described their organization in the reference genome FA1090 and has inferred the associated methylated motifs. Here, we studied the structure of RMS and target methylated motifs in 25 gonococcal strains sequenced with Single Molecule Real-Time (SMRT) technology, which provides data on DNA modification. The results showed a variable picture of active RMS in different strains, with phase variation switching the activity of Type III RMS, and both the activity and specificity of a Type I RMS. Interestingly, the Dam methylase was found in place of the NgoAXI endonuclease in two of the strains, despite being previously thought to be absent in the gonococcus. We also identified the real methylation target of NgoAXII as 5'-GCAGA-3', different from that previously described. Results from this work give further insights into the diversity and dynamics of RMS and methylation patterns in N. gonorrhoeae.

Evaluation of parameters affecting performance and reliability of machine learning-based antibiotic susceptibility testing from whole genome sequencing data.

Prediction of antibiotic resistance phenotypes from whole genome sequencing data by machine learning methods has been proposed as a promising platform for the development of sequence-based diagnostics. However, there has been no systematic evaluation of factors that may influence performance of such models, how they might apply to and vary across clinical populations, and what the implications might be in the clinical setting. Here, we performed a meta-analysis of seven large Neisseria gonorrhoeae datasets, as well as Klebsiella pneumoniae and Acinetobacter baumannii datasets, with whole genome sequence data and antibiotic susceptibility phenotypes using set covering machine classification, random forest classification, and random forest regression models to predict resistance phenotypes from genotype. We demonstrate how model performance varies by drug, dataset, resistance metric, and species, reflecting the complexities of generating clinically relevant conclusions from machine learning-derived models. Our findings underscore the importance of incorporating relevant biological and epidemiological knowledge into model design and assessment and suggest that doing so can inform tailored modeling for individual drugs, pathogens, and clinical populations. We further suggest that continued comprehensive sampling and incorporation of up-to-date whole genome sequence data, resistance phenotypes, and treatment outcome data into model training will be crucial to the clinical utility and sustainability of machine learning-based molecular diagnostics.