Evidence of Recent Genomic Evolution in Gonococcal Strains With Decreased Susceptibility to Cephalosporins or Azithromycin in the United States, 2014-2016

Emergence of Neisseria gonorrhoeae Clone with Reduced Susceptibility to Sitafloxacin in China: An In Vitro and Genomic Study

Drug-resistant Neisseria gonorrhoeae poses an urgent threat to public health. Recently, sitafloxacin, a new-generation fluoroquinolone, has shown high in vitro activity against drug-resistant N. gonorrhoeae. However, data on its effectiveness in clinical isolates remains limited. In this study, we collected 507 N. gonorrhoeae isolates from 21 hospitals in Shanghai, China, during 2020 and 2021. Antimicrobial susceptibility testing revealed that sitafloxacin minimum inhibitory concentrations (MICs) exhibited a bimodal distribution, ranging from <0.004 to 2 mg/L. The MIC50 and MIC90 for sitafloxacin were 0.125 mg/L and 0.5 mg/L, respectively, which are 32 and 16 times lower than those for ciprofloxacin (4 mg/L and 8 mg/L, respectively). Sitafloxacin demonstrated high in vitro activity against isolates resistant to either ceftriaxone, azithromycin, or both. Notably, among the isolates with reduced sitafloxacin susceptibility (MIC ≥ MIC90), 83.7% (36/43) were identified as sequence type (ST) 8123. Further phylogenetic analysis showed that ST8123 has evolved into two subclades, designated as subclade-I and subclade-II. A majority of the isolates (80%, 36/45) within subclade-I exhibited reduced susceptibility to sitafloxacin. In contrast, all isolates from subclade-II were found to be susceptible to sitafloxacin. Subsequent genomic investigations revealed that the GyrA-S91F, D95Y, and ParC-S87N mutations, which were exclusively found in ST8123 subclade-I, might be linked to reduced sitafloxacin susceptibility. Our study reveals that sitafloxacin is a promising antibiotic for combating drug-resistant N. gonorrhoeae. However, caution is advised in the clinical application of sitafloxacin for treating N. gonorrhoeae infections due to the emergence of a clone exhibiting reduced susceptibility.

Genomic analysis of 1710 surveillance-based Neisseria gonorrhoeae isolates from the USA in 2019 identifies predominant strain types and chromosomal antimicrobial-resistance determinants

This study characterized high-quality whole-genome sequences of a sentinel, surveillance-based collection of 1710 Neisseria gonorrhoeae (GC) isolates from 2019 collected in the USA as part of the Gonococcal Isolate Surveillance Project (GISP). It aims to provide a detailed report of strain diversity, phylogenetic relationships and resistance determinant profiles associated with reduced susceptibilities to antibiotics of concern. The 1710 isolates represented 164 multilocus sequence types and 21 predominant phylogenetic clades. Common genomic determinants defined most strains' phenotypic, reduced susceptibility to current and historic antibiotics (e.g. bla _TEM plasmid for penicillin, tetM plasmid for tetracycline, gyrA for ciprofloxacin, 23S rRNA and/or mosaic mtr operon for azithromycin, and mosaic penA for cefixime and ceftriaxone). The most predominant phylogenetic clade accounted for 21 % of the isolates, included a majority of the isolates with low-level elevated MICs to azithromycin (2.0 µg ml^-1), carried a mosaic mtr operon and variants in PorB, and showed expansion with respect to data previously reported from 2018. The second largest clade predominantly carried the GyrA S91F variant, was largely ciprofloxacin resistant (MIC ≥1.0 µg ml^-1), and showed significant expansion with respect to 2018. Overall, a low proportion of isolates had medium- to high-level elevated MIC to azithromycin ((≥4.0 µg ml^-1), based on C2611T or A2059G 23S rRNA variants). One isolate carried the penA 60.001 allele resulting in elevated MICs to cefixime and ceftriaxone of 1.0 µg ml^-1. This high-resolution snapshot of genetic profiles of 1710 GC sequences, through a comparison with 2018 data (1479 GC sequences) within the sentinel system, highlights change in proportions and expansion of select GC strains and the associated genetic mechanisms of resistance. The knowledge gained through molecular surveillance may support rapid identification of outbreaks of concern. Continued monitoring may inform public health responses to limit the development and spread of antibiotic-resistant gonorrhoea.

Repeated Emergence of Variant TetR Family Regulator, FarR, and Increased Resistance to Antimicrobial Unsaturated Fatty Acid among Clonal Complex 5 Methicillin-Resistant Staphylococcus aureus

Resistance-nodulation-division (RND) superfamily efflux pumps promote antibiotic resistance in Gram-negative pathogens, but their role in Gram-positive pathogens, including methicillin-resistant Staphylococcus aureus (MRSA) is undocumented. However, recent in vitro selections for resistance of S. aureus to an antimicrobial fatty acid, linoleic acid, and an antibiotic, rhodomyrtone, identified H121Y and C116R substitution variants, respectively, in a TetR family regulator, FarR, promoting increased expression of the RND pump FarE. Hypothesizing that in vivo selection pressures have also promoted the emergence of FarR variants, we searched available genome data and found that strains with FarR^H121Y from human and bovine hosts have emerged sporadically in clonal complexes (CCs) CC1, CC30, CC8, CC22, and CC97, whereas multiple FarR variants have occurred within CC5 hospital-associated (HA)-MRSA. Of these, FarR^E160G and FarR^E93EE were exclusive to CC5, while FarR^C116Y, FarR^P165L, and FarR^G166D also occurred in nonrelated CCs, primarily from bovine hosts. Within CC5, FarR^C116Y and FarR^G166D strains were polyphyletic, each exhibiting two emergence events. FarR^C116Y and FarR^E160G were individually sufficient to confer increased expression of FarE and enhanced resistance to linoleic acid (LA). Isolates with FarR^E93EE were most closely related to S. aureus N315 MRSA and exhibited increased resistance independently of FarR^E93EE. Accumulation of pseudogenes and additional polymorphisms in FarR^E93EE strains contributed to a multiresistance phenotype which included fosfomycin and fusidic acid resistance in addition to increased linoleic acid resistance. These findings underscore the remarkable adaptive capacity of CC5 MRSA, which includes the polyphyletic USA100 lineage of HA-MRSA that is endemic in the Western hemisphere and known for the acquisition of multiple resistance phenotypes.

Canary in the Coal Mine: How Resistance Surveillance in Commensals Could Help Curb the Spread of AMR in Pathogenic Neisseria

Antimicrobial resistance (AMR) is widespread within Neisseria gonorrhoeae populations. Recent work has highlighted the importance of commensal Neisseria (cN) as a source of AMR for their pathogenic relatives through horizontal gene transfer (HGT) of AMR alleles, such as mosaic penicillin binding protein 2 (penA), multiple transferable efflux pump (mtr), and DNA gyrase subunit A (gyrA) which impact beta-lactam, azithromycin, and ciprofloxacin susceptibility, respectively. However, nonpathogenic commensal species are rarely characterized. Here, we propose that surveillance of the universally carried commensal Neisseria may play the role of the "canary in the coal mine," and reveal circulating known and novel antimicrobial resistance determinants transferable to pathogenic Neisseria. We summarize the current understanding of commensal Neisseria as an AMR reservoir, and call to increase research on commensal Neisseria species, through expanding established gonococcal surveillance programs to include the collection, isolation, antimicrobial resistance phenotyping, and whole-genome sequencing (WGS) of commensal isolates. This will help combat AMR in the pathogenic Neisseria by: (i) determining the contemporary AMR profile of commensal Neisseria, (ii) correlating AMR phenotypes with known and novel genetic determinants, (iii) qualifying and quantifying horizontal gene transfer (HGT) for AMR determinants, and (iv) expanding commensal Neisseria genomic databases, perhaps leading to the identification of new drug and vaccine targets. The proposed modification to established Neisseria collection protocols could transform our ability to address AMR N. gonorrhoeae, while requiring minor modifications to current surveillance practices. IMPORTANCE Contemporary increases in the prevalence of antimicrobial resistance (AMR) in Neisseria gonorrhoeae populations is a direct threat to global public health and the effective treatment of gonorrhea. Substantial effort and financial support are being spent on identifying resistance mechanisms circulating within the gonococcal population. However, these surveys often overlook a known source of resistance for gonococci-the commensal Neisseria. Commensal Neisseria and pathogenic Neisseria frequently share DNA through horizontal gene transfer, which has played a large role in rendering antibiotic therapies ineffective in pathogenic Neisseria populations. Here, we propose the expansion of established gonococcal surveillance programs to integrate a collection, AMR profiling, and genomic sequencing pipeline for commensal species. This proposed expansion will enhance the field's ability to identify resistance in and from nonpathogenic reservoirs and anticipate AMR trends in pathogenic Neisseria.

Molecular Mechanisms of Drug Resistance and Epidemiology of Multidrug-Resistant Variants of Neisseria gonorrhoeae

The paper presents various issues related to the increasing drug resistance of Neisseria gonorrhoeae and the occurrence and spread of multidrug-resistant clones. One of the most important is the incidence and evolution of resistance mechanisms of N. gonorrhoeae to beta-lactam antibiotics. Chromosomal resistance to penicillins and oxyimino-cephalosporins and plasmid resistance to penicillins are discussed. Chromosomal resistance is associated with the presence of mutations in the PBP2 protein, containing mosaic variants and nonmosaic amino acid substitutions in the transpeptidase domain, and their correlation with mutations in the mtrR gene and its promoter regions (the MtrCDE membrane pump repressor) and in several other genes, which together determine reduced sensitivity or resistance to ceftriaxone and cefixime. Plasmid resistance to penicillins results from the production of beta-lactamases. There are different types of beta-lactamases as well as penicillinase plasmids. In addition to resistance to beta-lactam antibiotics, the paper covers the mechanisms and occurrence of resistance to macrolides (azithromycin), fluoroquinolones and some other antibiotics. Moreover, the most important epidemiological types of multidrug-resistant N. gonorrhoeae, prevalent in specific years and regions, are discussed. Epidemiological types are defined as sequence types, clonal complexes and genogroups obtained by various typing systems such as NG-STAR, NG-MAST and MLST. New perspectives on the treatment of N. gonorrhoeae infections are also presented, including new drugs active against multidrug-resistant strains.

First National Genomic Epidemiological Study of Neisseria gonorrhoeae Strains Spreading Across Sweden in 2016

The increasing transmission and antimicrobial resistance (AMR) in Neisseria gonorrhoeae is a global health concern with worrying trends of decreasing susceptibility to also the last-line extended-spectrum cephalosporin (ESC) ceftriaxone. A dramatic increase of reported gonorrhea cases has been observed in Sweden from 2016 and onward. The aim of the present study was to comprehensively investigate the genomic epidemiology of all cultured N. gonorrhoeae isolates in Sweden during 2016, in conjunction with phenotypic AMR and clinical and epidemiological data of patients. In total, 1279 isolates were examined. Etest and whole-genome sequencing (WGS) were performed, and epidemiological data obtained from the Public Health Agency of Sweden. Overall, 51.1%, 1.7%, and 1.3% resistance to ciprofloxacin, cefixime, and azithromycin, respectively, was found. No isolates were resistant to ceftriaxone, however, 9.3% of isolates showed a decreased susceptibility to ceftriaxone and 10.5% to cefixime. In total, 44 penA alleles were found of which six were mosaic (n = 92). Using the typing schemes of MLST, NG-MAST, and NG-STAR; 133, 422, and 280 sequence types, respectively, and 93 NG-STAR clonal complexes were found. The phylogenomic analysis revealed two main lineages (A and B) with lineage A divided into two main sublineages (A1 and A2). Resistance and decreased susceptibility to ESCs and azithromycin and associated AMR determinants, such as mosaic penA and mosaic mtrD, were predominantly found in sublineage A2. Resistance to cefixime and azithromycin was more prevalent among heterosexuals and MSM, respectively, and both were predominantly spread through domestic transmission. Continuous surveillance of the spread and evolution of N. gonorrhoeae, including phenotypic AMR testing and WGS, is essential for enhanced knowledge regarding the dynamic evolution of N. gonorrhoeae and gonorrhea epidemiology.

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.

Antimicrobial resistance prediction in Neisseria gonorrhoeae: Current status and future prospects

INTRODUCTION

Several nucleic acid amplification tests (NAATs), mostly real-time PCRs, to detect antimicrobial resistance (AMR) determinants and predict AMR in Neisseria gonorrhoeae are promising, and some may be ready to apply at the point-of-care (POC), but important limitations remain with most NAATs. Next-generation sequencing (NGS) can overcome many of these limitations.Areas covered: Recent advances, with main focus on publications since 2017, in the development and use of NAATs and NGS to predict gonococcal AMR for surveillance and clinical use, and pros and cons of these tests as well as future perspectives for appropriate use of molecular AMR prediction for N. gonorrhoeae.Expert Commentary: NAATs and/or NGS for AMR prediction should supplement culture-based AMR surveillance, which will remain because it detects also AMR due to unknown AMR determinants, and translation into POC tests is imperative for the end-goal of individualized treatment, sparing ceftriaxone±azithromycin. Several challenges for direct testing of clinical, especially pharyngeal, specimens and for accurate prediction of cephalosporins and azithromycin resistance, especially using NAATs, remain. The choice of AMR prediction assay needs to carefully consider the intended use of the assay; limitations intrinsic to the AMR prediction technology, algorithms and specific to chosen methodology; specimen types analyzed; and cost-effectiveness.

Demographic and Epidemiological Characteristics Associated With Reduced Antimicrobial Susceptibility to Neisseria gonorrhoeae in the United States, Strengthening the US Response to Resistant Gonorrhea, 2018 to 2019

BACKGROUND

Jurisdictions participating in Strengthening the US Response to Resistant Gonorrhea (SURRG) implemented specimen collection for culture and antimicrobial susceptibility testing from a sample of persons of all genders (at multiple anatomic sites) attending sexually transmitted disease clinics and community clinics. We describe the percentage and characteristics of patients whose isolates demonstrated reduced susceptibility (RS) to azithromycin, ceftriaxone, or cefixime.

METHODS

We included patients from clinics that participated in SURRG whose isolates underwent antimicrobial susceptibility testing by Etest. We defined RS as azithromycin minimum inhibitory concentrations (MICs) ≥2 μg/mL (AZM-RS), ceftriaxone MICs ≥0.125 μg/mL (CRO-RS), or cefixime MICs ≥0.25 μg/mL (CFX-RS). Patients with repeated infections appeared >1 time in the data. We calculated the frequency and percentage of patients with an isolate demonstrating RS by epidemiological characteristics.

RESULTS

During the period 2018-2019, 10,013 patients from 8 jurisdictions provided 10,735 isolates. Among 10,013 patients, 11.0% (n = 1099) had ≥1 isolate with AZM-RS (range by jurisdiction, 2.5%-18.0%). Approximately 11.3% of 8771 of patients visiting sexually transmitted disease clinics and approximately 8.8% of 1242 patients visiting community clinics had an AZM-RS isolate. Nearly 6% of 1013 females had an AZM-RS isolate; among males, the percents of patients with an AZM-RS isolate were 17.7% among 4177 men who have sex only with men and 6.1% among 3581 men who have sex only with women. Few (0.4%) patients had isolates with CFX-RS (n = 40) or CRO-RS (n = 43).

CONCLUSIONS

Although infections with reduced cephalosporin susceptibility were rare, AZM-RS infections were prevalent in this sample of patients in multiple jurisdictions and across gender and gender of sex partner categories.

Prediction of ceftriaxone MIC in Neisseria gonorrhoeae using DNA microarray technology and regression analysis

BACKGROUND

Decreased susceptibility of Neisseria gonorrhoeae to extended-spectrum cephalosporins is a major concern. Elucidation of the phenotypic and genetic characteristics of such isolates is a priority task.

METHODS

We developed a method for predicting the N. gonorrhoeae ceftriaxone susceptibility level (MICcro) by identifying genetic determinants of resistance using low-density hydrogel microarrays and a regression equation. A training dataset, containing 5631 isolates from the Pathogenwatch database and 181 isolates obtained in the Russian Federation during 2018-19, was used to build a regression model. The regression equation was tested on 14 WHO reference strains. Ceftriaxone resistance determinants for the 448 evaluated clinical isolates collected in Russia were identified using microarray analysis, and MICcro values were calculated using the regression equation and compared with those measured by the serial dilution method.

RESULTS

The regression equation for calculating MICcro values included 20 chromosomal resistance determinants. The greatest contributions to the increase in MICcro were shown to be PBP2: Ala-501→Pro, Ala-311→Val, Gly-545→Ser substitutions, Asp(345-346) insertion; and PorB: Gly-120→Arg substitution. The substitutions PBP2: Ala-501→Thr/Val, PorB: Gly-120→Asn/Asp/Lys and PBP1: Leu-421→Pro had weaker effects. For 94.4% of the isolates in the evaluation set, the predicted MICcro was within one doubling dilution of the experimentally determined MICcro. No ceftriaxone-resistant isolates were identified in the analysed samples from Russia, and no interpretative errors were detected in the MICcro calculations.

CONCLUSIONS

The developed strategy for predicting ceftriaxone MIC can be used for the continuous surveillance of known and emerging resistant N. gonorrhoeae isolates.

Genomic Epidemiology of Azithromycin-Nonsusceptible Neisseria gonorrhoeae, Argentina, 2005-2019

Azithromycin-nonsusceptible Neisseria gonorrhoeae strains are an emerging global public health threat. During 2015–2018, the prevalence of azithromycin-nonsusceptible gonococcal infection increased significantly in Argentina. To investigate the genomic epidemiology and resistance mechanisms of these strains, we sequenced 96 nonsusceptible isolates collected in Argentina during 2005–2019. Phylogenomic analysis revealed 2 main clades, which were characterized by a limited geographic distribution, circulating during January 2015–November 2019. These clades included the internationally spreading multilocus sequence types (STs) 1580 and 9363. The ST1580 isolates, which had MICs of 2–4 μg/mL, had mutations in the 23S rRNA. The ST9363 isolates, which had MICs of 2–4 or >256 μg/mL, had mutations in the 23S rRNA, a mosaic mtr locus, or both. Identifying the geographic dissemination and characteristics of these predominant clones will guide public health policies to control the spread of azithromycin-nonsusceptible N. gonorrhoeae in Argentina.

The Accuracy of Molecular Detection Targeting the Mutation C2611T for Detecting Moderate-Level Azithromycin Resistance in Neisseria gonorrhoeae: A Systematic Review and Meta-Analysis

BACKGROUND

Neisseria gonorrhoeae (N. gonorrhoeae) is now recognized as a commonly reported sexually transmitted pathogen, and the increasing drug resistance of N. gonorrhoeae has become a serious public health problem. The accuracy of molecular detection for detecting moderate-level azithromycin resistance is not well-established. We summarized the data from studies of the N. gonorrhoeae 23S rRNA mutation at position 2611 with azithromycin resistance to determine the relationship between the mutation and resistance.

METHODS AND FINDINGS

In this systematic review and meta-analysis, two researchers independently searched six databases for studies with data for the azithromycin minimum inhibitory concentrations (MICs) and the 23S rRNA mutation C2611T of each N. gonorrhoeae isolate. Since the breakpoint of moderate-level resistance to azithromycin (ML-AzmR) was not determined, we divided the moderate level into two groups according to the range of MICs (moderate resistance limited to 2-128 mg/L or 4-128 mg/L) for data extraction. A random-effects model was used to calculate the pooled sensitivity rate, the specificity rate, the pooled positive likelihood ratio (PLR), the negative likelihood ratio (NLR), and the diagnostic odds ratio (DOR). Meta-regression analyses by detection method, isolates sampling (a random sample or not), location, and sample size were performed to explore the possible causes of heterogeneity. The potential publication bias of the included studies was conducted by the Deeks' test. We included 20 studies in our study: 20 studies have data of N. gonorrhoeae with MICs between 2 and 128 mg/L with mutation or without mutation at position 2611(4759 samples), and 14 studies have data of N. gonorrhoeae with MICs between 4 and 128 mg/L (3367 samples). In the group with the moderate level of 2-128 mg/L, the pooled sensitivity rate of the molecular assays was determined to be 71.9% (95% CI, 67.6-74%), the pooled specificity rate was 98.7% (95% CI, 98.2-99.0%), and the DOR ranged from 55.0 to 351.3 (mean, 139.1). In the 4-128 mg/L group, the pooled sensitivity rate was 91.9% (95% CI, 88.9-94.2%), the pooled specificity rate was 95.9% (95% CI, 95.1-96.6%), and the DOR ranged from 41.9 to 364.1 (mean, 123.6).

CONCLUSION

Through this meta-analysis, we found that the C2611T mutation of 23S rRNA is valuable for the molecular diagnostic of moderate-level azithromycin resistance (ML-AzmR) in N. gonorrhoeae, especially when the moderate level is set at 4-128 mg/L. This rapid molecular detection method can be used for the rapid identification of ML-AzmR isolates in the clinic.

Genomic Analysis of the Predominant Strains and Antimicrobial Resistance Determinants Within 1479 Neisseria gonorrhoeae Isolates From the US Gonococcal Isolate Surveillance Project in 2018

BACKGROUND

The prevalence of Neisseria gonorrhoeae (GC) isolates with elevated minimum inhibitory concentrations to various antibiotics continues to rise in the United States and globally. Genomic analysis provides a powerful tool for surveillance of circulating strains, antimicrobial resistance determinants, and understanding of transmission through a population.

METHODS

Neisseria gonorrhoeae isolates collected from the US Gonococcal Isolate Surveillance Project in 2018 (n = 1479) were sequenced and characterized. Whole-genome sequencing was used to identify sequence types, antimicrobial resistance profiles, and phylogenetic relationships across demographic and geographic populations.

RESULTS

Genetic characterization identified that (1) 80% of the GC isolates were represented in 33 multilocus sequence types, (2) isolates clustered in 23 major phylogenetic clusters with select phenotypic and demographic prevalence, and (3) common antimicrobial resistance determinants associated with low-level or high-level decreased susceptibility or resistance to relevant antibiotics.

CONCLUSIONS

Characterization of this 2018 Gonococcal Isolate Surveillance Project genomic data set, which is the largest US whole-genome sequence data set to date, sets the basis for future prospective studies, and establishes a genomic baseline of GC populations for local and national monitoring.

Emergence and dissemination of three mild outbreaks of Neisseria gonorrhoeae with high-level resistance to azithromycin in Barcelona, 2016-18

BACKGROUND

Neisseria gonorrhoeae (NG) isolates with high-level azithromycin resistance (HL-AziR) have emerged worldwide in recent decades, threatening the sustainability of current dual-antimicrobial therapy.

OBJECTIVES

This study aimed to characterize the first 16 NG isolates with HL-AziR in Barcelona between 2016 and 2018.

METHODS

WGS was used to identify the mechanisms of antimicrobial resistance, to establish the MLST ST, NG multiantigen sequence typing (NG-MAST) ST and NG sequence typing for antimicrobial resistance (NG-STAR) ST and to identify the clonal relatedness of the isolates with other closely related NG previously described in other countries based on a whole-genome SNP analysis approach. The sociodemographic characteristics of the patients included in the study were collected by comprehensive review of their medical records.

RESULTS

Twelve out of 16 HL-AziR isolates belonged to the MLST ST7823/NG-MAST ST5309 genotype and 4 to MLST ST9363/NG-MAST ST3935. All presented the A2059G mutation in all four alleles of the 23S rRNA gene. MLST ST7823/NG-MAST ST5309 isolates were only identified in men who have sex with women and MLST ST9363/NG-MAST ST3935 were found in MSM. Phylogenomic analysis revealed the presence of three transmission clusters of three different NG strains independently associated with sexual behaviour.

CONCLUSIONS

Our findings support the first appearance of three mild outbreaks of NG with HL-AziR in Spain. These results highlight the continuous capacity of NG to develop antimicrobial resistance and spread among sexual networks. The enhanced resolution of WGS provides valuable information for outbreak investigation, complementing the implementation of public health measures focused on the prevention and dissemination of MDR NG.

Bioinformatics tools used for whole-genome sequencing analysis of Neisseria gonorrhoeae: a literature review

Whole-genome sequencing (WGS) data are well established for the investigation of gonococcal transmission, antimicrobial resistance prediction, population structure determination and population dynamics. A variety of bioinformatics tools, repositories, services and platforms have been applied to manage and analyze Neisseria gonorrhoeae WGS datasets. This review provides an overview of the various bioinformatics approaches and resources used in 105 published studies (as of 30 April 2021). The challenges in the analysis of N. gonorrhoeae WGS datasets, as well as future bioinformatics requirements, are also discussed.

Utility of MALDI-TOF MS for differentiation of Neisseria gonorrhoeae isolates with dissimilar azithromycin susceptibility profiles

BACKGROUND

Antibiotic-resistant gonorrhoea has been a chronic public health burden since the mid-1930s. Recent emergence of isolates resistant to the current recommended antibiotics for gonorrhoea further magnifies the threat of untreatable gonorrhoea. The lack of new, effective antibiotics highlights the need for better understanding of the population structure of Neisseria gonorrhoeae in order to provide greater insight on how to curtail the spread of antimicrobial-resistant N. gonorrhoeae.

OBJECTIVES

To explore a potential application of MALDI-TOF MS to differentiate N. gonorrhoeae displaying different levels of susceptibility to the antibiotic azithromycin.

METHODS

We conducted MALDI-TOF MS using the Bruker Biotyper on 392 N. gonorrhoeae isolates collected through the Gonococcal Isolate Surveillance Project (GISP) and/or the Strengthening the United States Response to Resistant Gonorrhea (SURRG) project. The MALDI-TOF MS spectra were visually analysed to assess the presence of distinctive peak(s). Statistical analysis was performed to assess the relationship between gonococcal isolates with the distinct protein peak and antibiotic susceptibility.

RESULTS

In this study, we were able to differentiate N. gonorrhoeae isolates into two distinct subpopulations using MALDI-TOF MS. Isolates were distinguished by the presence or absence of a spectral peak at 11 300 Da. Notably, these two groups exhibited different levels of susceptibility to azithromycin.

CONCLUSIONS

We have shown that in addition to its ability to identify N. gonorrhoeae, MALDI-TOF MS could also be used to differentiate gonococcal isolates with different levels of susceptibility to azithromycin.

Evidence of Horizontal Gene Transfer of 50S Ribosomal Genes rplB, rplD, and rplY in Neisseria gonorrhoeae

Horizontal gene transfer (HGT) in the penA and multidrug efflux pump genes has been shown to play a key role in the genesis of antimicrobial resistance in Neisseria gonorrhoeae. In this study, we evaluated if there was evidence of HGT in the genes coding for the ribosomal proteins in the Neisseria genus. We did this in a collection of 11,659 isolates of Neisseria, including N. gonorrhoeae and commensal Neisseria species (N. cinerea, N. elongata, N. flavescens, N. mucosa, N. polysaccharea, and N. subflava). Comparative genomic analyses identified HGT events in three genes: rplB, rplD, and rplY coding for ribosomal proteins L2, L4 and L25, respectively. Recombination events were predicted in N. gonorrhoeae and N. cinerea, N. subflava, and N. lactamica were identified as likely progenitors. In total, 2,337, 2,355, and 1,127 isolates possessed L2, L4, and L25 HGT events. Strong associations were found between HGT in L2/L4 and the C2597T 23S rRNA mutation that confers reduced susceptibility to macrolides. Whilst previous studies have found evidence of HGT of entire genes coding for ribosomal proteins in other bacterial species, this is the first study to find evidence of HGT-mediated chimerization of ribosomal proteins.

First Case of High-Level Azithromycin-Resistant Neisseria gonorrhoeae in North Carolina

We report on the first high-level azithromycin-resistant Neisseria gonorrhoeae isolate (minimum inhibitory concentration, ≥256 μg/mL) in North Carolina isolated from a pharyngeal swab of a 33-year-old HIV-negative man who has sex with men. In addition, the isolate was found to be susceptible to cefixime, ceftriaxone, and penicillin and resistant to tetracycline. By whole-genome sequencing, the strain was assigned as MLST ST9363, NG-MAST ST5035, and a novel NG-STAR sequence type, ST1993.

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.

Draft Genome Sequences of Three Penicillin-Resistant Neisseria gonorrhoeae Strains Isolated in Cincinnati, Ohio, in 1994

Here, we report the draft genome sequences of three penicillin-resistant Neisseria gonorrhoeae isolates. We include associated data on MICs and genetic relationships to other N. gonorrhoeae strains collected from across the United States. Resistance mutations known to contribute to reduced penicillin susceptibility are annotated in each genome.

Here, we report the draft genome sequences of three penicillin-resistant Neisseria gonorrhoeae isolates. We include associated data on MICs and genetic relationships to other N. gonorrhoeae strains collected from across the United States. Resistance mutations known to contribute to reduced penicillin susceptibility are annotated in each genome.

The mosaic mtr locus as major genetic determinant of azithromycin resistance of Neisseria gonorrhoeae, Germany, 2018

Within the German Gonococcal Resistance Network (GORENET) Neisseria gonorrhoeae (NG) sample collection, azithromycin-resistant NG isolates increased from 4.3% in 2016 to 9.2% in 2018. We aim to understand this observed increase using whole genome sequencing of NG isolates combined with epidemiological and clinical data. Reduced susceptibility to azithromycin in 2018 was predominately clonal (NG-MAST G12302) and could mainly be accounted to the recently described mosaic-like mtr locus. Our data suggest that, together with horizontal gene transfer of resistance determinants and well-established point mutations, international spread of resistant lineages plays a major role regarding azithromycin resistance in Germany.

Cyr SB, Pham CD, Kirkcaldy RD. Sustained Transmission of Neisseria gonorrhoeae with High-Level Resistance to Azithromycin, in Indianapolis, Indiana, 2017–2018

Background

Since 2014, Neisseria gonorrhoeae azithromycin (AZM) susceptibility has declined in the United States, but high-level AZM resistance (HL-AZMR) has been infrequent and sporadic. We describe a cluster of 14 N. gonorrhoeae isolates with HL-AZMR identified in Indianapolis over 13 months.

Methods

N. gonorrhoeae culture specimens (genital and extragenital) were collected from attendees of the Bell Flower Clinic. Isolates underwent antimicrobial susceptibility testing (AST) using Etest. AZM minimum inhibitory concentrations ≥256 µg/mL were classified as HL-AZMR. Local disease intervention specialists interviewed patients whose isolates demonstrated HL-AZMR and conducted partner services. Relatedness of isolates was investigated by genomic analyses.

Results

During 2017–2018, AST was performed in 1016 N. gonorrhoeae isolates collected at the Bell Flower Clinic. Fourteen isolates (1.4%) from 12 men collected over 13 months demonstrated HL-AZMR; all were cephalosporin susceptible. Of the 12 men, 9 were white and reported male sex partners. Nine of the men were able to be retested; all were cured with 250-mg ceftriaxone plus 1-g AZM. Two men named each other as partners; no other partners in common were reported. Genomic analysis demonstrated close relatedness of the HL-AZMR isolates and a novel combination of a mosaic-mtrR promoter along with 23S ribosomal RNA mutations that appear to have emerged from circulating strains.

Conclusions

The close genetic relatedness with limited epidemiologic linkages between patients highlights the challenges of gonorrhea partner investigations and suggests undetected local transmission. Local AST, rapid public health action, and epidemiologic investigations combined with genomic analysis provides a multipronged approach to understanding an outbreak of sexually transmitted disease.

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.

The distribution and spread of susceptible and resistant Neisseria gonorrhoeae across demographic groups in a major metropolitan center

Background

Genomic epidemiology studies of gonorrhea in the United States have primarily focused on national surveillance for antibiotic resistance, and patterns of local transmission between demographic groups of resistant and susceptible strains are unknown.

Methods

We analyzed a convenience sample of genome sequences, antibiotic susceptibility, and patient data from 897 gonococcal isolates cultured at the New York City (NYC) Public Health Laboratory from NYC Department of Health and Mental Hygiene (DOHMH) Sexual Health Clinic (SHC) patients, primarily in 2012–2013. We reconstructed the gonococcal phylogeny, defined transmission clusters using a 10 nonrecombinant single nucleotide polymorphism threshold, tested for clustering of demographic groups, and placed NYC isolates in a global phylogenetic context.

Results

The NYC gonococcal phylogeny reflected global diversity with isolates from 22/23 of the prevalent global lineages (96%). Isolates clustered on the phylogeny by patient sexual behavior (P < .001) and race/ethnicity (P < .001). Minimum inhibitory concentrations were higher across antibiotics in isolates from men who have sex with men compared to heterosexuals (P < .001) and white heterosexuals compared to black heterosexuals (P < .01). In our dataset, all large transmission clusters (≥10 samples) of N. gonorrhoeae were susceptible to ciprofloxacin, ceftriaxone, and azithromycin, and comprised isolates from patients across demographic groups.

Conclusions

All large transmission clusters were susceptible to gonorrhea therapies, suggesting that resistance to empiric therapy was not a main driver of spread, even as risk for resistance varied across demographic groups. Further study of local transmission networks is needed to identify drivers of transmission.

Azithromycin susceptibility of Neisseria gonorrhoeae in the USA in 2017: a genomic analysis of surveillance data

Background

The number of cases of gonorrhoea in the USA and worldwide caused by Neisseria gonorrhoeae is increasing (555 608 reported US cases in 2017, and 87 million cases worldwide in 2016). Many countries report declining in vitro susceptibility of azithromycin, which is a concern because azithromycin and ceftriaxone are the recommended dual treatment in many countries. We aimed to identify strain types associated with decreased susceptibility to azithromycin.

Methods

We did a genomic analysis of N gonorrhoeae isolates obtained by the US Gonococcal Isolate Surveillance Project. Isolates were whole-genome sequenced based on decreased susceptibility to azithromycin (minimal inhibitory concentration [MIC] ≥2 μg/mL, using agar dilution antibiotic susceptibility testing) and geographical representation. Bioinformatic analyses established genomic diversity, strain population dynamics, and antimicrobial resistance profiles.

Findings

410 isolates were sorted into more than 20 unique phylogenetic clades. One predominant persistent clade (consisting of 97 isolates) included the most isolates with azithromycin MICs of 2 μg/mL or higher (61 of 97 [63%] vs 59 of 311 [19%]; p<0·0001) and carried a mosaic mtr (multiple transferable resistance) locus (68 of 97 [70%] vs two of 313 [1%]; p<0·0001). Of the remaining 313 isolates, 57 (18%) had decreased susceptibility to azithromycin (MIC ≥4 μg/mL), which was attributed to 23S rRNA variants (56 of 57 [98%]) and formed phylogenetically diverse clades, showing various levels of clonal expansion.

Interpretation

Reduced azithromycin susceptibility was associated with expanding and persistent clades harbouring two well described resistance mechanisms, mosaic mtr locus and 23S rRNA variants. Understanding the role of recombination, particularly within the mtr locus, on the fitness and expansion of strains with decreased susceptibility has important implications for the public health response to minimise gonorrhoea transmission.

Funding

US Centers for Disease Control and Prevention (CDC), CDC Combating Antibiotic Resistant Bacteria initiative, Oak Ridge Institute for Science Education, US Department of Energy/CDC/Emory University, National Institutes of Health, and Biomedical Laboratory Research and Development Service of the US Department of Veterans Affairs.

Transcriptional control of the gonococcal ompA gene by the MisR/MisS two-component regulatory system

Neisseria gonorrhoeae, the causative agent of gonorrhea, is an exclusive human pathogen whose growing antibiotic resistance is causing worldwide concern. The increasing rise of antibiotic resistance expressed by gonococci highlights the need to find alternative approaches to current gonorrhea treatment such as vaccine development or novel therapeutics. The gonococcal OmpA protein was previously identified as a potential vaccine candidate due to its conservation and stable expression amongst strains of Neisseria gonorrhoeae. However, factors that might modulate levels of OmpA and therefore potential vaccine efficacy are unknown. Earlier work indicated that ompA is part of the MisR/MisS regulon and suggested that it was a MisR-activated gene. Herein, we confirmed MisR/MisS regulation of ompA and report that the MisR response regulator can bind upstream of the ompA translational start codon. Further, we describe the contribution of a DNA sequence upstream of the ompA promoter that is critical for MisR activation of ompA transcription. Our results provide a framework for understanding the transcription of gonococcal ompA through a regulatory system known to be important for survival of gonococci during experimental infection.

Genomic Characterization of Neisseria gonorrhoeae Strains from 2016 U.S. Sentinel Surveillance Displaying Reduced Susceptibility to Azithromycin

In 2016, the proportion of Neisseria gonorrhoeae isolates with reduced susceptibility to azithromycin rose to 3.6%. A phylogenetic analysis of 334 N. gonorrhoeae isolates collected in 2016 revealed a single, geographically diverse lineage of isolates with MICs of 2 to 16 μg/ml that carried a mosaic-like mtr locus, whereas the majority of isolates with MICs of ≥16 μg/ml appeared sporadically and carried 23S rRNA mutations. Continued molecular surveillance of N. gonorrhoeae isolates will identify new resistance mechanisms.

Genomic epidemiology and population structure of Neisseria gonorrhoeae in Norway, 2016-2017

This study presents the nationwide epidemiology of Neisseria gonorrhoeae, using whole-genome sequencing of all culture-positive cases, which comprise roughly 40 % of all cases of gonorrhea reported in Norway from 2016 to 2017. Isolates were assigned to sequence types and Bayesian analysis clusters and variation in genes coding for antibiotic resistance was linked to phenotypic resistance data. The study also included isolates taken from the same patients from different anatomical sites at one or more time points. Comparing these isolates allows for observation of patterns of infections, i.e. multiple reinfections of genetically related clones vs. reinfections of genetically distant clones, and quantification of the genomic variation of closely related isolates from samples taken from a patient within the same day. Demographically, the patients in the study could be split into two groups; one group of patients from the capital with a high proportion of men who have sex with men (MSM), and another consisting of young adults with transmission primarily between males and females from outside the capital. Some clusters of N. gonorrhoeae were restricted to one of these two demographic groups. Pairwise comparison of multiple isolates from the same patients revealed that most were reinfected with different clones. Observations of frequent reinfections in patients is a concern and should be taken into account in the development of improved information and treatment guidelines.

Expanding U.S. Laboratory Capacity for Neisseria gonorrhoeae Antimicrobial Susceptibility Testing and Whole-Genome Sequencing through the CDC's Antibiotic Resistance Laboratory Network

U.S. gonorrhea rates are rising, and antibiotic-resistant Neisseria gonorrhoeae (AR-Ng) is an urgent public health threat. Since implementation of nucleic acid amplification tests for N. gonorrhoeae identification, the capacity for culturing N. gonorrhoeae in the United States has declined, along with the ability to perform culture-based antimicrobial susceptibility testing (AST). Yet AST is critical for detecting and monitoring AR-Ng. In 2016, the CDC established the Antibiotic Resistance Laboratory Network (AR Lab Network) to shore up the national capacity for detecting several resistance threats including N. gonorrhoeae AR-Ng testing, a subactivity of the CDC's AR Lab Network, is performed in a tiered network of approximately 35 local laboratories, four regional laboratories (state public health laboratories in Maryland, Tennessee, Texas, and Washington), and the CDC's national reference laboratory. Local laboratories receive specimens from approximately 60 clinics associated with the Gonococcal Isolate Surveillance Project (GISP), enhanced GISP (eGISP), and the program Strengthening the U.S. Response to Resistant Gonorrhea (SURRG). They isolate and ship up to 20,000 isolates to regional laboratories for culture-based agar dilution AST with seven antibiotics and for whole-genome sequencing of up to 5,000 isolates. The CDC further examines concerning isolates and monitors genetic AR markers. During 2017 and 2018, the network tested 8,214 and 8,628 N. gonorrhoeae isolates, respectively, and the CDC received 531 and 646 concerning isolates and 605 and 3,159 sequences, respectively. In summary, the AR Lab Network supported the laboratory capacity for N. gonorrhoeae AST and associated genetic marker detection, expanding preexisting notification and analysis systems for resistance detection. Continued, robust AST and genomic capacity can help inform national public health monitoring and intervention.

Bridging of Neisseria gonorrhoeae lineages across sexual networks in the HIV pre-exposure prophylaxis era

Whole genome sequencing (WGS) has been used to investigate transmission of Neisseria gonorrhoeae, but to date, most studies have not combined genomic data with detailed information on sexual behaviour to define the extent of transmission across population risk groups (bridging). Here, through combined epidemiological and genomic analysis of 2,186N. gonorrhoeae isolates from Australia, we show widespread transmission of N. gonorrhoeae within and between population groups. We describe distinct transmission clusters associated with men who have sex with men (MSM) and heterosexuals, and men who have sex with men and women (MSMW) are identified as a possible bridging population between these groups. Further, the study identifies transmission of N. gonorrhoeae between HIV-positive and HIV-negative individuals receiving pre-exposure prophylaxis (PrEP). Our data highlight several groups that can be targeted for interventions aimed at improving gonorrhoea control, including returning travellers, sex workers, and PrEP users.

Here, Williamson et al. combine epidemiological and genomic analysis of 2,186 Neisseria gonorrhoeae isolates from Australia and show that men who have sex with men and women are a possible ‘bridging’ population between men who have sex with men and heterosexuals.

Antibacterial evaluation of plants extracts against ampicillin-resistant Escherichia coli (E. coli) by microcalorimetry and principal component analysis

Antibiotics abuse has caused increased bacterial resistance, which severely limits the application of antibiotics to the treatment of bacterial infections. Therefore, it is urgent to develop new antibacterial drugs through other sources. Dracontomelon dao (Blanco) Merr. & Rolfe (Ren Mianzi in Chinese) is a traditional medicinal material derived from Anacardiaceae with a long history of treating various infectious diseases, such as decubitus and skin ulcers. Recent research has indicated that different extracts from the leaves of D. dao, especially the ethyl acetate (EtOAc) fraction containing flavonoids and phenolic acids, exhibit potent antibacterial activities. In this research, the combined anti-drug-resistant bacterial activities of these active ingredients were investigated. Six samples (S1-S6) were obtained from the EtOAc fraction of D. dao leaves. Microcalorimetric measurements and principal component analysis were performed on the in vitro samples. The results showed that all six samples had notable antibacterial activities. Specifically, sample S6 exhibited a prominent antibacterial effect, with an IC50 value of 84.3 μg mL-1, which was significantly lower than that of other samples. The relative contents of main flavonoids and phenolic acids in S6 sample were confirmed by UPLC/Q-TOF-MS. In conclusion, sample S6 from the EtOAc fraction of D. dao leaves could be used as a potential antimicrobial resource in the treatment of infectious diseases. This work provides an insight into the effect of traditional Chinese medicine on drug-resistant bacteria. Moreover, the purification and characterization of the chemical compounds from the sample S6 deserve further analysis.