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

Development and implementation of a Type I-C CRISPR-based programmable repression system for Neisseria gonorrhoeae

Clustered regularly interspaced short palindromic repeats (CRISPR) are prokaryotic adaptive immune systems regularly utilized as DNA-editing tools. While Neisseria gonorrhoeae does not have an endogenous CRISPR, the commensal species Neisseria lactamica encodes a functional Type I-C CRISPR-Cas system. We have established an isopropyl β-d-1-thiogalactopyranoside added (IPTG)-inducible, CRISPR interference (CRISPRi) platform based on the N. lactamica Type I-C CRISPR missing the Cas3 nuclease to allow locus-specific transcriptional repression. As proof of principle, we targeted a non-phase-variable version of the opaD gene. We show that CRISPRi can downregulate opaD gene and protein expression, resulting in bacterial inability to stimulate neutrophil oxidative responses and to bind to an N-terminal fragment of CEACAM1. Importantly, we used CRISPRi to effectively knockdown all the transcripts of all 11 opa genes using a five-spacer CRISPR array, allowing control of the entire phase-variable opa family in strain FA1090. We also report that repression is reversible following IPTG removal. Finally, we showed that the Type I-C CRISPRi system can conditionally reduce the expression of two essential genes. This CRISPRi system will allow the interrogation of every Gc gene, essential and non-essential, to study physiology and pathogenesis and aid in antimicrobial development.IMPORTANCEClustered regularly interspaced short palindromic repeats (CRISPR)-Cas systems have proven instrumental in genetically manipulating many eukaryotic and prokaryotic organisms. Despite its usefulness, a CRISPR system had yet to be developed for use in Neisseria gonorrhoeae (Gc), a bacterium that is the main etiological agent of gonorrhea infection. Here, we developed a programmable and IPTG-inducible Type I-C CRISPR interference (CRISPRi) system derived from the commensal species Neisseria lactamica as a gene repression system in Gc. As opposed to generating genetic knockouts, the Type I-C CRISPRi system allows us to block transcription of specific genes without generating deletions in the DNA. We explored the properties of this system and found that a minimal spacer array is sufficient for gene repression while also facilitating efficient spacer reprogramming. Importantly, we also show that we can use CRISPRi to knockdown genes that are essential to Gc that cannot normally be knocked out under laboratory settings. Gc encodes ~800 essential genes, many of which have no predicted function. We predict that this Type I-C CRISPRi system can be used to help categorize gene functions and perhaps contribute to the development of novel therapeutics for gonorrhea.

Effect on the Resistome of Dual vs Monotherapy for the Treatment of Neisseria gonorrhoeae: Results From a Randomized Controlled Trial (ResistAZM Trial)

Background

No randomized controlled trial (RCT) has compared the impact on the resistome of ceftriaxone (CRO) plus azithromycin (AZM) vs CRO for the treatment of Neisseria gonorrhoea (NG).

Methods

This was an open-label, single-center, RCT comparing the effect on the resistome of CRO plus AZM vs CRO for the treatment of NG. Men who have sex with men (MSM) with genital, anorectal, or pharyngeal NG infection were randomized into the CRO/AZM and CRO arms. Oral rinse and anorectal samples were taken for culture and resistome profiling at 2 visits (baseline and day 14). The primary outcome was the ratio of mean macrolide resistance determinants in anorectal samples from day 14 between arms.

Results

Twenty individuals were randomized into the CRO/AZM arm and 22 into the CRO arm. We found no significant difference in the mean macrolide resistance determinants in the day 14 anorectal samples between arms (ratio, 1.05; 95% CI, 0.55-1.83; P = .102). The prevalence of baseline macrolide resistance was high (CRO/AZM arm = 95.00%; CRO arm = 90.91%).

Conclusions

We could not demonstrate a significant effect of dual CRO/AZM therapy on the resistome compared with CRO alone, likely due to a high baseline resistance to AZM. Interventions to prevent the emergence of antimicrobial resistance in MSM are needed.

Substituted salicylic acid analogs offer improved potency against multidrug-resistant Neisseria gonorrhoeae and good selectivity against commensal vaginal bacteria

Drug-resistant Neisseria gonorrhoeae represents a major threat to public health; without new effective antibiotics, untreatable gonococcal infections loom as a real possibility. In a previous drug-repurposing study, we reported that salicylic acid had good potency against azithromycin-resistant N. gonorrhoeae. We now report that the anti-gonococcal activity in this scaffold is easily lost by inopportune substitution, but that select substituted naphthyl analogs (3b, 3o and 3p) have superior activity to salicylic acid itself. Furthermore, these compounds retained potency against multiple ceftriaxone- and azithromycin-resistant strains, exhibited rapid bactericidal activity against N. gonorrhoeae, and showed high tolerability to mammalian cells (CC50 > 128 µg/mL). Promisingly, these compounds also show very weak growth inhibition of commensal vaginal bacteria.

Sustained Transmission of Neisseria gonorrhoeae Strains with High-Level Azithromycin Resistance (MIC ≥ 256 μg/mL) in Argentina, 2018 to 2022

Azithromycin combined with ceftriaxone is the recommended dual therapy for uncomplicated gonorrhea in many countries. Nevertheless, the increasing prevalence of azithromycin resistance compromises the effectiveness of this treatment strategy. From 2018 to 2022, we collected 13 gonococcal isolates with high-level azithromycin resistance (MIC ≥ 256 μg/mL) across Argentina. Whole-genome sequencing revealed that these isolates were mainly represented by the internationally spreading Neisseria gonorrhoeae multi-antigen sequence typing (NG-MAST) genogroup G12302, containing the 23S rRNA A2059G mutation (in all four alleles) together with mosaic mtrD and mtrR promoter 2 loci. This information is important to develop targeted public health policies to control the spread of azithromycin-resistant N. gonorrhoeae in Argentina and internationally. IMPORTANCE Azithromycin resistance in Neisseria gonorrhoeae has been increasing in numerous populations worldwide, which is of concern, as azithromycin is part of the recommended dual treatment in many countries. Here, we report 13 N. gonorrhoeae isolates with high-level azithromycin resistance (MIC ≥ 256 μg/mL). This study observed that high-level azithromycin-resistant gonococcal strains have shown sustained transmission in Argentina and are related to the successful international clone NG-MAST G12302. Genomic surveillance together with real-time tracing and data-sharing networks will be crucial in controlling the spread of azithromycin resistance in gonococcus.

Peptide barcode of multidrug-resistant strains of Neisseria gonorrhoeae isolated from patients in Thailand

The emergence of multidrug-resistant strains of Neisseria gonorrhoeae constitutes a serious threat to public health. The present study aimed to investigate peptidome-based biomarkers of multidrug-resistant N. gonorrhoeae, using matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS) and liquid chromatography tandem mass spectrometry (LC-MS). The peptide barcode database of multidrug resistant N. gonorrhoeae was generated from the whole-cell peptides of 93 N. gonorrhoeae isolated from patients in Thailand. The dendrogram of 93 independent isolates of antibiotic-resistant N. gonorrhoeae revealed five distinct clusters including azithromycin resistance group (AZ), ciprofloxacin resistance group (C), ciprofloxacin and penicillin resistance group (CP), ciprofloxacin and tetracycline resistance group (CT), ciprofloxacin, penicillin and tetracycline resistance group (CPT). The peptidomes of all clusters were comparatively analyzed using a high-performance liquid chromatography-mass spectrometry method (LC-MS). Nine peptides derived from 9 proteins were highly expressed in AZ (p value < 0.05). These peptides also played a crucial role in numerous pathways and showed a strong relationship with the antibiotic resistances. In conclusion, this study showed a rapid screening of antibiotic-resistant N. gonorrhoeae using MALDI-TOF MS. Additionally, potential specific peptidome-based biomarker candidates for AZ, C, CP, CT and CPT-resistant N. gonorrhoeae were identified.

High-level in vitro resistance to gentamicin acquired in a stepwise manner in Neisseria gonorrhoeae

OBJECTIVES

Gentamicin is used in several alternative treatments for gonorrhoea. Verified clinical Neisseria gonorrhoeae isolates with gentamicin resistance are mainly lacking and understanding the mechanisms for gonococcal gentamicin resistance is imperative. We selected gentamicin resistance in gonococci in vitro, identified the novel gentamicin-resistance mutations, and examined the biofitness of a high-level gentamicin-resistant mutant.

METHODS

Low- and high-level gentamicin resistance was selected in WHO X (gentamicin MIC = 4 mg/L) on gentamicin-gradient agar plates. Selected mutants were whole-genome sequenced. Potential gentamicin-resistance fusA mutations were transformed into WT strains to verify their impact on gentamicin MICs. The biofitness of high-level gentamicin-resistant mutants was examined using a competitive assay in a hollow-fibre infection model.

RESULTS

WHO X mutants with gentamicin MICs of up to 128 mg/L were selected. Primarily selected fusA mutations were further investigated, and fusAR635L and fusAM520I + R635L were particularly interesting. Different mutations in fusA and ubiM were found in low-level gentamicin-resistant mutants, while fusAM520I was associated with high-level gentamicin resistance. Protein structure predictions showed that fusAM520I is located in domain IV of the elongation factor-G (EF-G). The high-level gentamicin-resistant WHO X mutant was outcompeted by the gentamicin-susceptible WHO X parental strain, suggesting lower biofitness.

CONCLUSIONS

We describe the first high-level gentamicin-resistant gonococcal isolate (MIC = 128 mg/L), which was selected in vitro through experimental evolution. The most substantial increases of the gentamicin MICs were caused by mutations in fusA (G1560A and G1904T encoding EF-G M520I and R635L, respectively) and ubiM (D186N). The high-level gentamicin-resistant N. gonorrhoeae mutant showed impaired biofitness.

Whole genome analysis and antimicrobial resistance of Neisseria gonorrhoeae isolates from Ghana

Introduction

Gonorrhoea is a major public health concern. With the global emergence and spread of resistance to last-line antibiotic treatment options, gonorrhoea threatens to be untreatable in the future. Therefore, this study performed whole genome characterization of Neisseria gonorrhoeae collected in Ghana to identify lineages of circulating strains as well as their phenotypic and genotypic antimicrobial resistance (AMR) profiles.

Methods

Whole genome sequencing (WGS) was performed on 56 isolates using both the Oxford Nanopore MinION and Illumina MiSeq sequencing platforms. The Comprehensive Antimicrobial Resistance Database (CARD) and PUBMLST.org/neisseria databases were used to catalogue chromosomal and plasmid genes implicated in AMR. The core genome multi-locus sequence typing (cgMLST) approach was used for comparative genomics analysis.

Results and Discussion

In vitro resistance measured by the E-test method revealed 100%, 91.0% and 85.7% resistance to tetracycline, penicillin and ciprofloxacin, respectively. A total of 22 sequence types (STs) were identified by multilocus sequence typing (MLST), with ST-14422 (n = 10), ST-1927 (n = 8) and ST-11210 (n = 7) being the most prevalent. Six novel STs were also identified (ST-15634, 15636-15639 and 15641). All isolates harboured chromosomal AMR determinants that confer resistance to beta-lactam antimicrobials and tetracycline. A single cefixime-resistant strain, that belongs to N. gonorrhoeae multiantigen sequence type (NG-MAST) ST1407, a type associated with widespread cephalosporin resistance was identified. Neisseria gonorrhoeae Sequence Typing for Antimicrobial Resistance (NG-STAR), identified 29 unique sequence types, with ST-464 (n = 8) and the novel ST-3366 (n = 8) being the most prevalent. Notably, 20 of the 29 STs were novel, indicative of the unique nature of molecular AMR determinants in the Ghanaian strains. Plasmids were highly prevalent: pTetM and pblaTEM were found in 96% and 92% of isolates, respectively. The TEM-135 allele, which is an amino acid change away from producing a stable extended-spectrum β-lactamase that could result in complete cephalosporin resistance, was identified in 28.5% of the isolates. Using WGS, we characterized N. gonorrhoeae strains from Ghana, giving a snapshot of the current state of gonococcal AMR in the country and highlighting the need for constant genomic surveillance.

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.

The Prevalence of Undiagnosed Salmonella enterica Serovar Typhi in Healthy School-Aged Children in Osun State, Nigeria

Typhoid fever remains a significant public health concern due to cases of mis-/overdiagnosis. Asymptomatic carriers play a role in the transmission and persistence of typhoid fever, especially among children, where limited data exist in Nigeria and other endemic countries. We aim to elucidate the burden of typhoid fever among healthy school-aged children using the best surveillance tool(s). In a semi-urban/urban state (Osun), 120 healthy school-aged children under 15 years were enrolled. Whole blood and fecal samples were obtained from consenting children. ELISA targeting the antigen lipopolysaccharide (LPS) and anti-LPS antibodies of Salmonella Typhi, culture, polymerase chain reaction (PCR), and next-generation sequencing (NGS) were used to analyze the samples. At least one of the immunological markers was detected in 65.8% of children, with 40.8%, 37.5%, and 39% of children testing positive for IgM, IgG, and antigen, respectively. Culture, PCR, and NGS assays did not detect the presence of Salmonella Typhi in the isolates. This study demonstrates a high seroprevalence of Salmonella Typhi in these healthy children but no carriage, indicating the inability to sustain transmission. We also demonstrate that using a single technique is insufficient for typhoid fever surveillance in healthy children living in endemic areas.

Machine learning models for Neisseria gonorrhoeae antimicrobial susceptibility tests

Neisseria gonorrhoeae is an urgent public health threat due to the emergence of antibiotic resistance. As most isolates in the United States are susceptible to at least one antibiotic, rapid molecular antimicrobial susceptibility tests (ASTs) would offer the opportunity to tailor antibiotic therapy, thereby expanding treatment options. With genome sequence and antibiotic resistance phenotype data for nearly 20,000 clinical N. gonorrhoeae isolates now available, there is an opportunity to use statistical methods to develop sequence-based diagnostics that predict antibiotic susceptibility from genotype. N. gonorrhoeae, therefore, provides a useful example illustrating how to apply machine learning models to aid in the design of sequence-based ASTs. We present an overview of this framework, which begins with establishing the assay technology, the performance criteria, the population in which the diagnostic will be used, and the clinical goals, and extends to the choices that must be made to arrive at a set of features with the desired properties for predicting susceptibility phenotype from genotype. While we focus on the example of N. gonorrhoeae, the framework generalizes to other organisms for which large-scale genotype and antibiotic resistance data can be combined to aid in diagnostics development.

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.

Etiological Surveillance of Male Urethritis Syndrome in South Africa: 2019 to 2020

BACKGROUND

In South Africa, male urethritis syndrome (MUS) is the most common sexually transmitted infection (STI) syndrome in men. We determined the distribution of STI etiologies and the susceptibility profiles of Neisseria gonorrhoeae isolates from men presenting with MUS to 3 sentinel surveillance health care facilities. Secondary objectives were to determine the seroprevalence of coinfections (HIV, syphilis, herpes simplex virus 2).

METHODS

Consecutive, consenting men with symptomatic urethral discharge were enrolled between January 1, 2019, and December 31, 2020. Genital discharge swab and blood specimens were collected and transported to a central STI reference laboratory in Johannesburg, South Africa.

RESULTS

Among 769 men enrolled, N. gonorrhoeae was the commonest cause of MUS (674 [87.8%]; 95% confidence interval [CI], 85.2%-89.9%), followed by Chlamydia trachomatis (161 [21.0%]; 95% CI, 18.2%-24.0%). Of 542 cultivable N. gonorrhoeae isolates, all were susceptible to ceftriaxone (modal minimum inhibitory concentration, 0.004 mg/L) and azithromycin (modal minimum inhibitory concentration, 0.128 mg/L). Seroprevalence rates of HIV, syphilis, and HSV-2 were 21.4% (95% CI, 18.5%-24.5%), 2.3%, and 50.1%, respectively. Condom use at last sexual encounter was reported by only 7%, less than 50% had been medically circumcised, and only 66.7% (58 of 87) who self-reported an HIV-positive status were adherent on antiretroviral drugs.

CONCLUSIONS

Neisseria gonorrhoeae and C. trachomatis were the predominant causes of MUS. Currently recommended dual ceftriaxone and azithromycin therapy are appropriate for MUS syndromic management; however, surveillance must be maintained to timeously detect emerging and increasing gonococcal resistance. Clinic-based interventions must be intensified in men seeing sexual health care to reduce the community transmission and burden of STI and HIV.

Transcriptional regulation of the mtrCDE efflux pump operon: importance for Neisseria gonorrhoeae antimicrobial resistance

This review focuses on the mechanisms of transcriptional control of an important multidrug efflux pump system (MtrCDE) possessed by Neisseria gonorrhoeae, the aetiological agent of the sexually transmitted infection termed gonorrhoea. The mtrCDE operon that encodes this tripartite protein efflux pump is subject to both cis- and trans-acting transcriptional factors that negatively or positively influence expression. Critically, levels of MtrCDE can influence levels of gonococcal susceptibility to classical antibiotics, host-derived antimicrobials and various biocides. The regulatory systems that control mtrCDE can have profound influences on the capacity of gonococci to resist current and past antibiotic therapy regimens as well as virulence. The emergence, mechanisms of action and clinical significance of the transcriptional regulatory systems that impact mtrCDE expression in gonococci are reviewed here with the aim of linking bacterial antimicrobial resistance with multidrug efflux capability.

Phylogenomic comparison of Neisseria gonorrhoeae causing disseminated gonococcal infections and uncomplicated gonorrhea in Georgia, United States

Disseminated gonococcal infection (DGI) is a rare complication caused by the systemic dissemination of Neisseria gonorrhoeae (Ng) to normally sterile anatomical sites. Little is known about the genetic diversity of DGI gonococcal strains and how they relate to other gonococcal strains causing uncomplicated mucosal infections. We used whole genome sequencing to characterize DGI isolates (n = 30) collected from a surveillance system in Georgia (GA), USA during 2017-2020 to understand phylogenetic clustering among DGI as well as uncomplicated uro-and-extragenital gonococcal (UGI) isolates (n = 110) collected in Fulton County, GA during 2017-2019. We also investigated the presence or absence of genetic markers related to antimicrobial resistance (AMR) as well as surveyed the genomes for putative virulence genetic factors associated with normal human-serum (NHS) resistance that might facilitate DGI. We found that DGI strains demonstrated significant genetic variability similar to the population structure of isolates causing UGI, with sporadic incidences of geographically clustered DGI strains. DGI isolates contained various AMR markers and genetic mechanisms associated with NHS resistance. DGI isolates had a higher frequency of the porB1A allele compared with UGI (67% vs. 9%, p &lt; 0.0001); however, no single NHS resistance marker was found in all DGI isolates. Continued DGI surveillance with genome-based characterization of DGI isolates is necessary to better understand specific factors that promote systemic dissemination.

Controversies and evidence on Chlamydia testing and treatment in asymptomatic women and men who have sex with men: a narrative review

Background

Chlamydia trachomatis (CT) is the most common bacterial sexually transmitted infection (STI) worldwide. CT is mainly asymptomatic. Test-and-treat strategies are widely implemented to prevent transmission and complications. Strategies are not without controversy in asymptomatic women and men who have sex with men (MSM). Concerns are emerging to test and treat asymptomatic persons for urogenital CT (‘Controversy 1’) and pharyngeal or rectal CT (‘Controversy 2’), whereby testing symptomatic persons is not under debate. Opposed views in CT treatment involve using azithromycin versus doxycycline (‘Controversy 3’). The objective of this review is to provide coverage of these public health and clinical controversies by reviewing the current scientific evidence.

Methods

A literature search was performed using PubMed for relevant publications between 2018 and September 2021, and iterative retrieval of additional relevant publications.

Results

Controversy 1. In women, the majority of asymptomatic CT are at the urogenital site, and detections mostly include viable CT. CT easily transmits to a partner and potentially also between the vaginal and rectal areas; the clinical impact of urogenital CT is established, although risks for adverse outcomes are uncertain. Wide-scale testing in asymptomatic women has not resulted in reduced prevalence. In MSM, evidence for the clinical impact of asymptomatic urogenital CT is lacking. Controversy 2. Rectal CT is common in women diagnosed with urogenital CT, but the clinical impact of asymptomatic rectal CT is uncertain. In MSM, rectal CT is common, and most CT infections are at the rectal site, yet the risk of longer term complications is unknown. In both sexes, pharyngeal CT is uncommon and has no documented clinical impact. Controversy 3. In the treatment of rectal CT, doxycycline has superior effectiveness to azithromycin. Evidence has also accumulated on the harms of test-and-treat strategies.

Conclusions

Current practices vary widely, from widescale test-and-treat approaches to more individual patient- and partner-level case management. Choosing which asymptomatic people to test at what anatomic site, and whether to test or not, requires an urgent (re-)definition of the goals of testing and treating asymptomatic persons. Treatment guidelines are shifting toward universal doxycycline use, and clinical practice now faces the challenge of implementation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12879-022-07171-2.

Test-and-treat is a key strategy in the control of Chlamydia trachomatis (CT). We discuss recent controversies and present scientific evidence regarding urogenital, rectal, and pharyngeal CT test-and-treat strategies in women and in men who have sex with men (MSM). This should inform best practices for the prevention and management of the most common bacterial sexually transmitted infection (STI) worldwide.

Gonococcal Clinical Strains Bearing a Common gdhR Single Nucleotide Polymorphism That Results in Enhanced Expression of the Virulence Gene lctP Frequently Possess a mtrR Promoter Mutation That Decreases Antibiotic Susceptibility

GdhR is a transcriptional repressor of the virulence factor gene lctP, which encodes a unique l-lactate permease that has been linked to pathogenesis of Neisseria gonorrhoeae, and loss of gdhR can confer increased fitness of gonococci in a female mouse model of lower genital tract infection. In this work, we identified a single nucleotide polymorphism (SNP) in gdhR, which is often present in both recent and historical gonococcal clinical strains and results in a proline (P)-to-serine (S) change at amino acid position 6 (P6S) of GdhR. This mutation (gdhR6) was found to reduce GdhR transcriptional repression at lctP in gonococcal strains containing the mutant protein compared to wild-type GdhR. By using purified recombinant proteins and in vitro DNA-binding and cross-linking experiments, we found that gdhR6 impairs the DNA-binding activity of GdhR at lctP without an apparent effect on protein oligomerization. By analyzing a panel of U.S. (from 2017 to 2018) and Danish (1928 to 2013) clinical isolates, we observed a statistical association between gdhR6 and the previously described adenine deletion in the promoter of mtrR (mtrR-P A-del), encoding the repressor (MtrR) of the mtrCDE operon that encodes the MtrCDE multidrug efflux pump that can export antibiotics, host antimicrobials, and biocides. The frequent association of gdhR6 with the mtrR promoter mutation in these clinical isolates suggests that it has persisted in this genetic background to enhance lctP expression, thereby promoting virulence.

Antimicrobial susceptibility of commensal Neisseria in a general population and men who have sex with men in Belgium

Non-pathogenic Neisseria are a reservoir of antimicrobial resistance genes for pathogenic Neisseria meningitidis and Neisseria gonorrhoeae. Men who have sex with men (MSM) are at risk of co-colonization with resistant non-pathogenic and pathogenic Neisseria. We assessed if the antimicrobial susceptibility of non-pathogenic Neisseria among MSM differs from a general population and if antimicrobial exposure impacts susceptibility. We recruited 96 participants at our center in Belgium: 32 employees, 32 MSM who did not use antibiotics in the previous 6 months, and 32 MSM who did. Oropharyngeal Neisseria were cultured and identified with MALDI-TOF–MS. Minimum inhibitory concentrations for azithromycin, ceftriaxone and ciprofloxacin were determined using E-tests^® and compared between groups with non-parametric tests. Non-pathogenic Neisseria from employees as well as MSM were remarkably resistant. Those from MSM were significantly less susceptible than employees to azithromycin and ciprofloxacin (p < 0.0001, p < 0.001), but not ceftriaxone (p = 0.3). Susceptibility did not differ significantly according to recent antimicrobial exposure in MSM. Surveilling antimicrobial susceptibility of non-pathogenic Neisseria may be a sensitive way to assess impact of antimicrobial exposure in a population. The high levels of antimicrobial resistance in this survey indicate that novel resistance determinants may be readily available for future transfer from non-pathogenic to pathogenic Neisseria.

Increasing azithromycin resistance in Neisseria gonorrhoeae due to NG-MAST 12302 clonal spread in Canada, 2015-2018

Objectives: Azithromycin resistant (AZIR) gonorrhea has been steadily increasing in Canada over the past decade which is cause for alarm as azithromycin (AZI) has been part of the combination therapy recommended by the Canadian Guidelines on Sexually Transmitted Infections (CGSTI) since 2012.

Method: Neisseria gonorrhoeae (NG) with AZI MICs ≥ 1 mg/L collected between 2015 and 2018 as part of the Gonococcal Antimicrobial Surveillance Program-Canada underwent antimicrobial susceptibility testing, molecular typing and whole genome sequencing. Regional, demographic and clinical isolation site comparisons were made to aid in our understanding of AZI susceptibility trending.

Results: 3,447 NG with AZI MICs ≥ 1 mg/L were identified in Canada, increasing from 6.3% in 2015 to 26.5% of isolates in 2018. Central Canada had the highest proportion rising from 9.2% in 2015 to 31.2% in 2018. 273 different NG-MAST sequence types were identified among these isolates with ST-12302 the most prevalent (50.9%). Whole genome sequencing identified the Neisseria lactamica -like mosaic mtr locus as the mechanism of AZIR in isolates of ST-12302 and isolates genetically similar (differ by ≤ 5 base pairs) designated as the ST-12302 genogroup, accounting for 65.2% of study isolateswhich were originally identified in central Canada but spread to other regions by 2018.

Conclusion: Genomic analysis indicated that AZIR in Canadian NG expanded rapidly due to clonal spread of the ST-12302 genogroup. The rapid expansion of this AZIR clonal group in all regions of Canada is of concern. CGSTI are currently under review to address the increase in AZIR in Canada.

Design, synthesis and biological evaluation of P2-modified proline analogues targeting the HtrA serine protease in Chlamydia

High temperature requirement A (HtrA) serine proteases have emerged as a novel class of antibacterial target, which are crucial in protein quality control and are involved in the pathogenesis of a wide array of bacterial infections. Previously, we demonstrated that HtrA in Chlamydia is essential for bacterial survival, replication and virulence. Here, we report a new series of proline (P2)-modified inhibitors of Chlamydia trachomatis HtrA (CtHtrA) developed by proline ring expansion and Cγ-substitutions. The structure-based drug optimization process was guided by molecular modelling and in vitro pharmacological evaluation of inhibitory potency, selectivity and cytotoxicity. Compound 25 from the first-generation 4-substituted proline analogues increased antiCtHtrA potency and selectivity over human neutrophil elastase (HNE) by approximately 6- and 12-fold, respectively, relative to the peptidic lead compound 1. Based on this compound, second-generation substituted proline residues containing 1,2,3-triazole moieties were synthesized by regioselective azide-alkyne click chemistry. Compound 49 demonstrated significantly improved antichlamydial activity in whole cell assays, diminishing the bacterial infectious progeny below the detection limit at the lowest dose tested. Compound 49 resulted in approximately 9- and 22-fold improvement in the inhibitory potency and selectivity relative to 1, respectively. To date, compound 49 is the most potent HtrA inhibitor developed against Chlamydia spp.

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.

Repurposing FDA-approved sulphonamide carbonic anhydrase inhibitors for treatment of Neisseria gonorrhoeae

Neisseria gonorrhoeae is a high-priority pathogen of concern due to the growing prevalence of resistance development against approved antibiotics. Herein, we report the anti-gonococcal activity of ethoxzolamide, the FDA-approved human carbonic anhydrase inhibitor. Ethoxzolamide displayed an MIC_50, against a panel of N. gonorrhoeae isolates, of 0.125 µg/mL, 16-fold more potent than acetazolamide, although both molecules exhibited almost similar potency against the gonococcal carbonic anhydrase enzyme (NgCA) in vitro. Acetazolamide displayed an inhibition constant (K_i) versus NgCA of 74 nM, while Ethoxzolamide^’s K_i was estimated to 94 nM. Therefore, the increased anti-gonococcal potency of ethoxzolamide was attributed to its increased permeability in N. gonorrhoeae as compared to that of acetazolamide. Both drugs demonstrated bacteriostatic activity against N. gonorrhoeae, exhibited post-antibiotic effects up to 10 hours, and resistance was not observed against both. Taken together, these results indicate that acetazolamide and ethoxzolamide warrant further investigation for translation into effective anti-N. gonorrhoeae agents.

Development and application of Cas13a-based diagnostic assay for Neisseria gonorrhoeae detection and azithromycin resistance identification

BACKGROUND

Gonorrhoea, caused by Neisseria gonorrhoeae, has spread worldwide. Strains resistant to most antibiotics, including ceftriaxone and azithromycin, have emerged to an alarming level. Rapid testing for N. gonorrhoeae and its antimicrobial resistance will therefore contribute to clinical decision making for early diagnosis and rational drug use.

METHODS

A Cas13a-based assay (specific high-sensitivity enzymatic reporter unlocking; SHERLOCK) was developed for N. gonorrhoeae detection (porA gene) and azithromycin resistance identification (A2059G, C2611T). Assays were evaluated for sensitivity with purified dsDNA and specificity with 17 non-gonococcal strains. Performance of SHERLOCK (porA) was compared with Roche Cobas 4800 using 43 urine samples. Identification of azithromycin resistance mutations (A2059G, C2611T) was evaluated using a total of 84 clinical isolates and 18 urine samples. Lateral flow was tested for this assay as a readout tool. Moreover, we directly assayed 27 urethral swabs from patients with urethritis to evaluate their status in terms of N. gonorrhoeae infection and azithromycin resistance.

RESULTS

The SHERLOCK assay was successfully developed with a sensitivity of 10 copies/reaction, except 100 copies/reaction for A2059G, and no cross-reaction with other species. Comparison of the SHERLOCK assay with the Cobas 4800 revealed 100% concordance within 18 positive and 25 negative urine samples. Of the 84 isolates, 21 strains with azithromycin resistance mutations were distinguished and further verified by sequencing and MIC determination. In addition, 62.96% (17/27) strains from swab samples were detected with no mutant strains confirmed by sequencing.

CONCLUSIONS

The SHERLOCK assay for rapid N. gonorrhoeae detection combined with azithromycin resistance testing is a promising method for application in clinical practice.

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.

Large Increase in Azithromycin-Resistant Neisseria gonorrhoeae in Northern Spain

The aim of this study was to characterize the evolution of gonorrhea in the general population by correlating epidemiological, genotypic, and antimicrobial resistance data of Neisseria gonorrhoeae isolates collected in northern Spain from 2014 to 2018. One hundred ninety-four strains underwent antimicrobial susceptibility testing and were genetically analyzed by N. gonorrhoeae multiantigen sequence typing. Increasing cases of gonococcal infections have been observed after 2015. Most occurred in male with urethritis. Sequence type (ST)-9972 and ST-1576, the predominant genotypes identified, have not been previously described as epidemic clones. Of great concern was the significant increase in azithromycin-resistant N. gonorrhoeae. More than 30% of these isolates were obtained from men who have sex with men (MSM). ST-12302 was the most prevalent clone among the azithromycin-resistant strains, and was also resistant to penicillin, ciprofloxacin, and tetracycline. This multidrug-resistant clone was exclusively isolated from MSM during 2018. The incidence rates of gonorrhea and azithromycin-resistant N. gonorrhoeae have significantly increased due to the emergence of new clones. ST-12302 has recently been recognized as an epidemic clone; therefore, its surveillance could be the key in controlling further dissemination of azithromycin resistance. These data highlight the need to perform local studies to update treatment guidelines and reinforce preventive measures against gonorrhea.

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.

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.

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.

Structures of Neisseria gonorrhoeae MtrR-operator complexes reveal molecular mechanisms of DNA recognition and antibiotic resistance-conferring clinical mutations

Mutations within the mtrR gene are commonly found amongst multidrug resistant clinical isolates of Neisseria gonorrhoeae, which has been labelled a superbug by the Centers for Disease Control and Prevention. These mutations appear to contribute to antibiotic resistance by interfering with the ability of MtrR to bind to and repress expression of its target genes, which include the mtrCDE multidrug efflux transporter genes and the rpoH oxidative stress response sigma factor gene. However, the DNA-recognition mechanism of MtrR and the consensus sequence within these operators to which MtrR binds has remained unknown. In this work, we report the crystal structures of MtrR bound to the mtrCDE and rpoH operators, which reveal a conserved, but degenerate, DNA consensus binding site 5'-MCRTRCRN4YGYAYGK-3'. We complement our structural data with a comprehensive mutational analysis of key MtrR-DNA contacts to reveal their importance for MtrR-DNA binding both in vitro and in vivo. Furthermore, we model and generate common clinical mutations of MtrR to provide plausible biochemical explanations for the contribution of these mutations to multidrug resistance in N. gonorrhoeae. Collectively, our findings unveil key biological mechanisms underlying the global stress responses of N. gonorrhoeae.

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.

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.