Abstract

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.

Introduction

Gonorrhoea is an important public health concern as the second most common global bacterial sexually transmitted infection (STI) with 82.4 million new cases estimated among adults aged 15–49 years annually.1 The causative bacterium, Neisseria gonorrhoeae (gonococcus), has developed resistance to all introduced gonorrhoea therapies since the introduction of antibiotics.2–4 Current first-line treatment for uncomplicated gonorrhoea is ceftriaxone, an extended-spectrum cephalosporin (ESC), either in a high-dose monotherapy or given in combination with azithromycin.5–8 Since 2018, sporadic N. gonorrhoeae strains with resistance to ceftriaxone combined with high-level resistance to azithromycin causing treatment failures have been reported.9–11 Furthermore, international transmission of the ceftriaxone-resistant gonococcal FC428 clone has been verified since 2015.12–21 National and international antimicrobial surveillance reports have also documented increases in azithromycin resistance and expansions of azithromycin-resistant clones in the gonococcal population.22–25 This spread of N. gonorrhoeae strains resistant to ceftriaxone and/or azithromycin jeopardize the last remaining gonorrhoea therapies.

In 2017, the WHO published a list of 12 antimicrobial-resistant (AMR) pathogens including N. gonorrhoeae, which pose the greatest threat to human health,26 promoting research and development of new antimicrobials, but equally important, repurposing and optimized use of existing drugs. Aminoglycosides, such as gentamicin and tobramycin, in combination therapy or monotherapy remain active against several of the Gram-negative priority pathogens published by the WHO, e.g. Pseudomonas aeruginosa and Escherichia coli.27 Verified clinical N. gonorrhoeae isolates expressing in vitro resistance to gentamicin (MIC > 16 mg/L)28 have also been mainly lacking.28–30 Accordingly, several alternative treatment regimens for uncomplicated gonorrhoea recommended by the WHO, and in European, UK and US guidelines include gentamicin [240 mg intramuscularly (IM)] as an alternative treatment in combination with azithromycin single oral 2 g dose.5–8,22

In the advent of more widespread resistance to ceftriaxone, the use of gentamicin for treatment of gonorrhoea may rapidly increase. Consequently, it is imperative to understand how N. gonorrhoeae could potentially develop gentamicin resistance. In different bacterial species, the resistance determinants for gentamicin have involved enzymatic modifications,31 ineffective transport into the cell (although not fully understood),32 and rare ribosomal modifications causing resistance to streptomycin and spectinomycin (an aminocyclitol).33 A recent report described for the first time a resistance determinant for low-level resistance to gentamicin (MIC = 32 mg/L) in N. gonorrhoeae, i.e. SNPs in the fusA gene resulting in amino acid alterations, especially A563V, but also G564D and V651F in the elongation factor G (EF-G).34

Herein, we selected in vitro stepwise gentamicin-resistance determinants in N. gonorrhoeae, causing the gentamicin-susceptible WHO X reference strain35 to express high-level gentamicin resistance, identified the novel gentamicin-resistance mutations, and examined the biofitness of the WHO X mutant with high-level gentamicin resistance using competitive coculture in a dynamic hollow-fibre infection model (HFIM) for gonorrhoea.36

Materials and methods

N. gonorrhoeae isolate and antimicrobial susceptibility testing

The N. gonorrhoeae reference strain WHO X35 (gentamicin MIC of 4 mg/L) was cultured and subcultured once on a GCAGP agar plate (3.6% Difco GC Medium Base agar; BD Diagnostics, Sparks, MD, USA) supplemented with 1% haemoglobin (BD Diagnostics), 1% IsoVitaleX (BD Diagnostics) and 10% horse serum37 prior to 100 serial culture passages (subcultures) on gentamicin-gradient GCAGP agar plates. All incubations were performed in a 5% CO2-enriched humid atmosphere at 36 ± 1°C for 24 h. MICs (mg/L) of gentamicin, kanamycin, streptomycin, tobramycin, amikacin, ceftriaxone, cefixime, ciprofloxacin, ampicillin, tetracycline, erythromycin, spectinomycin and azithromycin were determined using Etest (bioMérieux, Marcy-l’Étoile, France), as previously described.37 The agar dilution method, in accordance with recommendations by CLSI,38 was used to further verify the MICs (mg/L) of gentamicin, kanamycin and amikacin. No gonococcal resistance breakpoints for gentamicin have been stated by EUCAST or CLSI. Consequently, in the present study a gentamicin-resistance breakpoint (MIC > 16 mg/L) suggested for N. gonorrhoeae in previous studies was used.28

In vitro selection for gentamicin resistance

GCAGP agar was prepared in two batches, of which one additionally included gentamicin. Gentamicin solution (10 mg/mL) in deionized water (Sigma–Aldrich, Saint Louis, USA) was used to prepare different gentamicin concentrations. Gentamicin-gradient agar plates were prepared as previously described.39 Briefly, the first layer of GCAGP agar including gentamicin (approximately 13 mL) was poured into a Petri dish placed on a slight slope and this layer was allowed to solidify prior to adding the second GCAGP agar layer without gentamicin (approximately 13 mL) when the Petri dish was placed on an even surface. The agar when solidified contained up to 8× the gentamicin MIC of WHO X.

The first culture passage was performed by suspending the WHO X reference strain in sterile saline solution to McFarland 0.5 and 100 µL of this bacterial solution was inoculated on the gentamicin-gradient GCAGP agar plate followed by incubation (36 ± 1°C, 5% CO2-enriched humid atmosphere, 24 h). A single colony growing at the highest gentamicin concentration was picked and subcultured on GCAGP agar plate without gentamicin (36 ± 1°C, 5% CO2-enriched humid atmosphere, 24 h). This type of culture passage was repeated and Etest for gentamicin was performed at passages 0, 2, 4, 6, 8, 10, 30, 40, 50, 60, 80 and 100, referred to as WHO X0–WHO X100. As a control, the WHO X reference strain was also serially cultured for 100 passages on a GCVIT agar plate [3.6% Difco GC Medium Base agar (BD) supplemented with 1% IsoVitaleX (BD)] without gentamicin.

WGS and analysis

DNA was extracted from 12 WHO X0–100 isolates using the automated QIAsymphony DSP Virus/Pathogen kit (QIAGEN, Hilden, Germany). Illumina DNA Prep (Illumina, Inc., San Diego, CA, USA) libraries were quality controlled using Qubit (Thermo Fisher Scientific, Waltham, MA, USA) and TapeStation (Agilent Technologies, Santa Clara, CA, USA), according to manufacturer instructions prior to library normalization. Illumina MiSeq reagent kit v2 (Illumina) was used to produce 250 paired-end reads with an average coverage of >100×.

All reads were quality controlled, trimmed (Phred quality score Q30), and assembled using a customized CLC Genomics Workbench (v22.0.1).3 Furthermore, a pipeline to generate a multiple sequence alignment, phylogenetic tree and additional quality control reports was performed using Nullarbor (https://github.com/tseemann/nullarbor) with WHO X (accession: NZ_LT592155.1) as reference with the following parameters; –ref NZ_LT592155.1.fasta –assembler spades –assembler-opt ‘–careful’ –taxoner kraken2. The multiple sequence alignment was used to compute pairwise distances between the WHO X reference genome35 and WHO X0–100 using MEGA (v11.0.10).40 The phylogenetic tree was visualized using Microreact.41 All raw sequence reads are available through the European Nucleotide Archive (project accession number: PRJEB56638).

Mutations found in the gentamicin-resistance selected isolates were also queried for in an in-house database with 33 306 publicly available gonococcal genomes (downloaded 31 January 2023) to investigate the occurrences in the gonococcal population, as previously described.42

Transformation of fusA

To examine whether the novel fusA mutations contributed to the high-level resistance to gentamicin in WHO X, the full-length fusA gene was PCR amplified using the primers fusA_F (CCACGACGACAGAACGTATT) and fusA_R (TTACGGGCTTCAGTTACAGC). The amplified product was transformed into WHO M35 and WHO X35 using spot transformation with 0.1 μg of purified fusA PCR product. Transformants with increased MICs of gentamicin were selected on gentamicin-containing GCVIT agar plates (16 mg/L). WGS, as described above, was performed to verify all fusA transformants.

Biofitness experiments using competitive coculture in an HFIM

To evaluate the biofitness of the WHO X mutant with high-level gentamicin resistance (WHO X100 with EF-G R635L and M520I substitutions) selected on the gentamicin-gradient GCAGP37 agar plates compared with the gentamicin-susceptible WHO X parent strain, competition experiments using coculture were performed in an HFIM, as earlier described.36,43 Briefly, bacteria were harvested from GCAGP agar plates and suspended in modified fastidious broth (mFB)36,43 to a quantity of ∼105–106 cfu/mL. Equal volumes (5 mL) of the suspensions of each strain were inoculated into the same HFIM cartridge. Aliquots (1 mL) were sampled at 24, 48, 72, 96, 120, 144 and 168 h, serially diluted in mFB and quantitatively plated on GCAGP agar plates and GCAGP agar plates containing 3×MIC of gentamicin, as previously described.36,43 Colony counts (log10 cfu/mL) were quantified after incubation for up to 72 h at 36°C in a humidified 5% CO2-enriched humid atmosphere using an automated colony counter (Scan 4000, Interscience, Saint-Nom-la-Bretèche, France). The competitive index (CI)36,43,44 was determined by dividing the ratio of the WHO X100 (EF-G R635L and M520I) mutant to WHO X parent strain at each timepoint with the ratio of the WHO X100 (EF-G R635L and M520I) mutant to WHO X parent strain in the initial inoculum.

Protein structure prediction

The protein structure of WHO X EF-G was predicted using the primary structure of fusA as input to AlphaFold (v2.1.0) (DeepMind, London, UK) executed using Colaboratory (Google Research, San Francisco, USA). It was analysed as a monomer and queried using the full database (https://github.com/deepmind/alphafold).45 Pfam was used to identify protein domains46 and CASTp was used to identify binding pockets.47 PyMOL was used to visualize the protein, amino acids and used for mutagenesis analysis (v2.0) (Schrödinger, LLC).

Results

Long-term in vitro gentamicin exposure causes high-level resistance in N. gonorrhoeae

The MICs of gentamicin for WHO X and several other aminoglycosides after serial culture passages (subcultures) on gentamicin-gradient agar plates are summarized in Table 1. As documented earlier,30 the gentamicin MIC determined using Etest is frequently approximately one gentamicin MIC doubling dilution lower compared with the agar dilution gentamicin MIC. The first 10 passages resulted in only minor increases in the antimicrobial MICs (one MIC doubling dilution) in the isolates WHO X0–10 (Table 1). For example, the Etest MIC of gentamicin only increased from 4 mg/L in WHO X0 to 8 mg/L in WHO X10. A more substantial increase in gentamicin MIC was documented after 30 passages (WHO X30) when low-level gentamicin resistance (MIC = 32 mg/L) and substantially increased MICs of also the additional aminoglycosides (MIC = 32 to >256 mg/L) were determined. The WHO X30 mutant appeared to be relatively stable in both growth and MICs of gentamicin and the additional aminoglycosides. Accordingly, no further increases in the MICs of gentamicin were recorded until after 100 serial culture passages (WHO X100), when high-level gentamicin resistance (MIC = 128 mg/L) and Etest MICs ≥ 256 mg/L of kanamycin, tobramycin and amikacin were identified. The MIC of streptomycin was also substantially increased, i.e. 32 mg/L. Notably, the WHO X MICs of ceftriaxone (4 mg/L), cefixime (8 mg/L), azithromycin (1 mg/L), ciprofloxacin (>32 mg/L), ampicillin (2 mg/L), tetracycline (2 mg/L), erythromycin (2 mg/L) and spectinomycin (16 mg/L) did not change during the 100 serial culture passages.

Table 1.

Changes in MICs (mg/L) of five aminoglycosides for N. gonorrhoeae reference strain WHO X35 during 100 serial culture passages (subcultures) on gentamicin-gradient agar plates

Etest MIC (agar dilution MIC), mg/La,b
PassageGentamicinKanamycinStreptomycinTobramycinAmikacin
WHO X04 (8)16 (16)8432 (128)
WHO X24 (8)16 (32)16864 (64)
WHO X44 (16)16 (32)168128 (128)
WHO X64 (16)32 (32)16864 (64)
WHO X84 (16)32 (32)16864 (64)
WHO X108 (16)32 (32)16864 (128)
WHO X3032 (64)128 (128)3232>256 (256)
WHO X4032 (64)>256 (128)3264>256 (256)
WHO X5032 (64)>256 (128)3264>256 (256)
WHO X6032 (32)>256 (128)3264>256 (256)
WHO X8032 (64)>256 (256)3264>256 (256)
WHO X100128 (128)>256 (512)32256>256 (256)
Etest MIC (agar dilution MIC), mg/La,b
PassageGentamicinKanamycinStreptomycinTobramycinAmikacin
WHO X04 (8)16 (16)8432 (128)
WHO X24 (8)16 (32)16864 (64)
WHO X44 (16)16 (32)168128 (128)
WHO X64 (16)32 (32)16864 (64)
WHO X84 (16)32 (32)16864 (64)
WHO X108 (16)32 (32)16864 (128)
WHO X3032 (64)128 (128)3232>256 (256)
WHO X4032 (64)>256 (128)3264>256 (256)
WHO X5032 (64)>256 (128)3264>256 (256)
WHO X6032 (32)>256 (128)3264>256 (256)
WHO X8032 (64)>256 (256)3264>256 (256)
WHO X100128 (128)>256 (512)32256>256 (256)

MICs (mg/L) were determined using both Etest and agar dilution method (in parentheses).

For in vitro selected mutations associated with increased aminoglycosides MICs, see Figure 1.

Table 1.

Changes in MICs (mg/L) of five aminoglycosides for N. gonorrhoeae reference strain WHO X35 during 100 serial culture passages (subcultures) on gentamicin-gradient agar plates

Etest MIC (agar dilution MIC), mg/La,b
PassageGentamicinKanamycinStreptomycinTobramycinAmikacin
WHO X04 (8)16 (16)8432 (128)
WHO X24 (8)16 (32)16864 (64)
WHO X44 (16)16 (32)168128 (128)
WHO X64 (16)32 (32)16864 (64)
WHO X84 (16)32 (32)16864 (64)
WHO X108 (16)32 (32)16864 (128)
WHO X3032 (64)128 (128)3232>256 (256)
WHO X4032 (64)>256 (128)3264>256 (256)
WHO X5032 (64)>256 (128)3264>256 (256)
WHO X6032 (32)>256 (128)3264>256 (256)
WHO X8032 (64)>256 (256)3264>256 (256)
WHO X100128 (128)>256 (512)32256>256 (256)
Etest MIC (agar dilution MIC), mg/La,b
PassageGentamicinKanamycinStreptomycinTobramycinAmikacin
WHO X04 (8)16 (16)8432 (128)
WHO X24 (8)16 (32)16864 (64)
WHO X44 (16)16 (32)168128 (128)
WHO X64 (16)32 (32)16864 (64)
WHO X84 (16)32 (32)16864 (64)
WHO X108 (16)32 (32)16864 (128)
WHO X3032 (64)128 (128)3232>256 (256)
WHO X4032 (64)>256 (128)3264>256 (256)
WHO X5032 (64)>256 (128)3264>256 (256)
WHO X6032 (32)>256 (128)3264>256 (256)
WHO X8032 (64)>256 (256)3264>256 (256)
WHO X100128 (128)>256 (512)32256>256 (256)

MICs (mg/L) were determined using both Etest and agar dilution method (in parentheses).

For in vitro selected mutations associated with increased aminoglycosides MICs, see Figure 1.

Mutations in several genes including fusA cause high-level gentamicin resistance

The results of the pairwise genome comparisons between WHO X and the 11 WHO X2–100 isolates from serial culture passages on gentamicin-gradient agar plates are summarized in Table 2. In total, 19 SNPs were found in the seven WHO X10–100 isolates. However, 9 of these 19 SNPs (all in opa and pilE genes) were also detected in the gentamicin-unexposed WHO X control strain, which also passed 100 serial culture passages. This gentamicin-unexposed strain did not show any increased MIC of gentamicin and, accordingly, these SNPs were considered as normal SNPs in repetitive regions of the gonococcal genome and excluded from further analysis. Interestingly, a stepwise accumulation of SNPs in the isolates with gentamicin exposure was observed and could be phylogenomically discerned (Figure 1). The mutations had been acquired stepwise in the following order: the pilE gene, followed by ubiM, lptD, pglD, opa, unknown membrane protein, pgsA, fusA and finally cysB. The mutations fusA G1560A (EF-G M520I) and ubiM D186N resulted in the most substantial increase in gentamicin MICs (Figure 1).

Phylogenomic tree of the gentamicin-susceptible N. gonorrhoeae WHO X35 parent strain and the in vitro gentamicin-exposed WHO X2–100 mutants, including the gentamicin-resistant WHO X30–100 mutants. The coloured boxes show the different MICs (mg/L) of gentamicin and the black boxes represent presence of mutations in the genes. This figure appears in colour in the online version of JAC and in black and white in the print version of JAC.
Figure 1.

Phylogenomic tree of the gentamicin-susceptible N. gonorrhoeae WHO X35 parent strain and the in vitro gentamicin-exposed WHO X2–100 mutants, including the gentamicin-resistant WHO X30–100 mutants. The coloured boxes show the different MICs (mg/L) of gentamicin and the black boxes represent presence of mutations in the genes. This figure appears in colour in the online version of JAC and in black and white in the print version of JAC.

Table 2.

All SNPs (n = 19) selected by gentamicin exposure in the WHO X10–100 mutants

UniProt locus tag/geneMutation (genetic alterations)Protein/functionFunction
pglDG1367A (G456D)Protein-PII uridylyltransferasePilin glycosylation protein PglD
pgsAC211T (D71F)CDP-diacylglycerol–glycerol-3-phosphate 3-phosphatidyltransferasePhospholipid biosynthesis
NGO_0974C1008A (premature stop-codon)Membrane proteinLysophospholipid transporter
OpaG90A, G91A (-25 and -26 opa mutation)Opacity proteinOuter membrane protein
pilEA812C, C813G, T815C, T816A, G817T, A818C, G824A, C825A, G827A (N121T, F122S, D123S, S125K)Fimbrial proteinType IV major pilin protein PilE
cysBC550T (P184S)Transcriptional regulator CysB
ubiMG556A (D186N)5-demethoxyubiquinol-8 5-hydroxylase UbiMUbiquinone synthesis involved in electron transport and oxidative stress defence
lptDC2066T (S689F)LOS-assembly proteinLPS synthesis
fusAG1560A, G1904T (M520I, R635L)Elongation factor GCentral role in protein synthesis
UniProt locus tag/geneMutation (genetic alterations)Protein/functionFunction
pglDG1367A (G456D)Protein-PII uridylyltransferasePilin glycosylation protein PglD
pgsAC211T (D71F)CDP-diacylglycerol–glycerol-3-phosphate 3-phosphatidyltransferasePhospholipid biosynthesis
NGO_0974C1008A (premature stop-codon)Membrane proteinLysophospholipid transporter
OpaG90A, G91A (-25 and -26 opa mutation)Opacity proteinOuter membrane protein
pilEA812C, C813G, T815C, T816A, G817T, A818C, G824A, C825A, G827A (N121T, F122S, D123S, S125K)Fimbrial proteinType IV major pilin protein PilE
cysBC550T (P184S)Transcriptional regulator CysB
ubiMG556A (D186N)5-demethoxyubiquinol-8 5-hydroxylase UbiMUbiquinone synthesis involved in electron transport and oxidative stress defence
lptDC2066T (S689F)LOS-assembly proteinLPS synthesis
fusAG1560A, G1904T (M520I, R635L)Elongation factor GCentral role in protein synthesis
Table 2.

All SNPs (n = 19) selected by gentamicin exposure in the WHO X10–100 mutants

UniProt locus tag/geneMutation (genetic alterations)Protein/functionFunction
pglDG1367A (G456D)Protein-PII uridylyltransferasePilin glycosylation protein PglD
pgsAC211T (D71F)CDP-diacylglycerol–glycerol-3-phosphate 3-phosphatidyltransferasePhospholipid biosynthesis
NGO_0974C1008A (premature stop-codon)Membrane proteinLysophospholipid transporter
OpaG90A, G91A (-25 and -26 opa mutation)Opacity proteinOuter membrane protein
pilEA812C, C813G, T815C, T816A, G817T, A818C, G824A, C825A, G827A (N121T, F122S, D123S, S125K)Fimbrial proteinType IV major pilin protein PilE
cysBC550T (P184S)Transcriptional regulator CysB
ubiMG556A (D186N)5-demethoxyubiquinol-8 5-hydroxylase UbiMUbiquinone synthesis involved in electron transport and oxidative stress defence
lptDC2066T (S689F)LOS-assembly proteinLPS synthesis
fusAG1560A, G1904T (M520I, R635L)Elongation factor GCentral role in protein synthesis
UniProt locus tag/geneMutation (genetic alterations)Protein/functionFunction
pglDG1367A (G456D)Protein-PII uridylyltransferasePilin glycosylation protein PglD
pgsAC211T (D71F)CDP-diacylglycerol–glycerol-3-phosphate 3-phosphatidyltransferasePhospholipid biosynthesis
NGO_0974C1008A (premature stop-codon)Membrane proteinLysophospholipid transporter
OpaG90A, G91A (-25 and -26 opa mutation)Opacity proteinOuter membrane protein
pilEA812C, C813G, T815C, T816A, G817T, A818C, G824A, C825A, G827A (N121T, F122S, D123S, S125K)Fimbrial proteinType IV major pilin protein PilE
cysBC550T (P184S)Transcriptional regulator CysB
ubiMG556A (D186N)5-demethoxyubiquinol-8 5-hydroxylase UbiMUbiquinone synthesis involved in electron transport and oxidative stress defence
lptDC2066T (S689F)LOS-assembly proteinLPS synthesis
fusAG1560A, G1904T (M520I, R635L)Elongation factor GCentral role in protein synthesis

Transformation experiments were focused on fusA because an EF-G A563V substitution was recently shown to cause low-level resistance to gentamicin (MIC = 32 mg/L) in N. gonorrhoeae34 and the two EF-G amino acid substitutions selected in the present study appeared to substantially increase the MICs of gentamicin. Accordingly, the impact of these identified amino acid alterations in EF-G on the MICs of gentamicin and other aminoglycosides was examined, i.e. by transforming R635L alone (fusA G1904T) and R635L (fusA G1904T) in combination with M520I (fusA G1560A) to the WHO M and WHO X reference strains.35 The MICs of gentamicin and the other aminoglycosides in the transformants are summarized in Table 3. Briefly, the R635L transformants, i.e. WHO MR635L and WHO XR635L, showed a 2- to 3-fold increase in the gentamicin MIC (to 16 and 32 mg/L, respectively) as well as in the MICs of kanamycin, tobramycin, amikacin and streptomycin. The M520I + R635L transformants were both resistant to gentamicin (MIC = 32 mg/L) and showed high MICs (≥256 mg/L) of kanamycin, amikacin and streptomycin (Table 3).

Table 3.

MICs (mg/L) of aminoglycosides in N. gonorrhoeae reference strains WHO M35 and WHO X35 and corresponding transformants expressing elongation factor G R635L or M520I + R635L substitutions

Etest (mg/L)
IsolateGentamicinKanamycinAmikacinTobramycinStreptomycin
WHO M41632816
WHO MR635L1612825632128
WHO MM520I + R635L3225625632256
WHO X41632416
WHO XR635L3212812832128
WHO XM520I + R635L32>256>25632>256
Etest (mg/L)
IsolateGentamicinKanamycinAmikacinTobramycinStreptomycin
WHO M41632816
WHO MR635L1612825632128
WHO MM520I + R635L3225625632256
WHO X41632416
WHO XR635L3212812832128
WHO XM520I + R635L32>256>25632>256
Table 3.

MICs (mg/L) of aminoglycosides in N. gonorrhoeae reference strains WHO M35 and WHO X35 and corresponding transformants expressing elongation factor G R635L or M520I + R635L substitutions

Etest (mg/L)
IsolateGentamicinKanamycinAmikacinTobramycinStreptomycin
WHO M41632816
WHO MR635L1612825632128
WHO MM520I + R635L3225625632256
WHO X41632416
WHO XR635L3212812832128
WHO XM520I + R635L32>256>25632>256
Etest (mg/L)
IsolateGentamicinKanamycinAmikacinTobramycinStreptomycin
WHO M41632816
WHO MR635L1612825632128
WHO MM520I + R635L3225625632256
WHO X41632416
WHO XR635L3212812832128
WHO XM520I + R635L32>256>25632>256

The genes (locus_tag) cysB (C7S01_RS08790), fusA (C7S01_RS1176), lptD (C7S01_RS10025), pglD (C7S01_RS00555), pgsA (C7S01_RS01250) and ubiM (C7S01_RS09765) were queried among 33 306 publicly available gonococcal genomes for the in vitro selected mutations associated with gentamicin resistance (Table 2). None of the mutations selected in cysB, lptD, pglD, pgsA and ubiM were found. However, two isolates from Ontario, Canada contained a fusA R635G mutation,48 and these isolates would be of interest to investigate further.

In vitro selected high-level gentamicin-resistant N. gonorrhoeae isolates have impaired biofitness

In general, compared with colonies of the gentamicin-susceptible WHO X35 parent strain the colonies of the high-level gentamicin-resistant WHO X100 mutant were substantially smaller and slow-growing on the GCAGP agar plates, i.e. grew as SCVs. The WHO X100 mutant and WHO X parent strain were subsequently cultured in separate HFIM cartridges as well as cocultured in the same HFIM cartridge for 7 days to evaluate if the in vitro selected high-level gentamicin-resistant WHO X100 mutant showed any impaired bacterial growth and biofitness. The growth of WHO X was maintained at approximately 1010 cfu/mL during the 7 day experiment, while the growth of the WHO X100 mutant was substantially lower from 24 h and to the end of the 7 day experiment (Figure 2a). The CI further verified that the WHO X parent strain outcompeted the high-level gentamicin-resistant WHO X100 mutant. Accordingly, despite the WHO X100 mutant starting to recover in growth after 96 h, which may indicate emergence of compensatory mutations restoring biofitness, it clearly suffered from impaired biofitness compared to the WHO X parent strain (Figure 2b).

Growth curves for WHO X35 parent strain and the in vitro selected high-level gentamicin-resistant WHO X100 mutant and competition experiment using coculture in an HFIM for gonorrhoea.36,43 (a) shows the cfu/mL of WHO X compared with WHO X100 during growth for 7 days and (b) displays the CI36,43 during coculture for 7 days, which illustrates how WHO X is outcompeting WHO X100 (especially during the first 4 days). A negative value indicates that the in vitro selected high-level gentamicin-resistant WHO X100 mutant has a lower biofitness compared with the gentamicin-susceptible WHO X35 parent strain. This figure appears in colour in the online version of JAC and in black and white in the print version of JAC.
Figure 2.

Growth curves for WHO X35 parent strain and the in vitro selected high-level gentamicin-resistant WHO X100 mutant and competition experiment using coculture in an HFIM for gonorrhoea.36,43 (a) shows the cfu/mL of WHO X compared with WHO X100 during growth for 7 days and (b) displays the CI36,43 during coculture for 7 days, which illustrates how WHO X is outcompeting WHO X100 (especially during the first 4 days). A negative value indicates that the in vitro selected high-level gentamicin-resistant WHO X100 mutant has a lower biofitness compared with the gentamicin-susceptible WHO X35 parent strain. This figure appears in colour in the online version of JAC and in black and white in the print version of JAC.

Low- and high-level gentamicin resistance is associated with alterations in domain IV and V in EF-G

The 3D structure of EF-G has five domains (Figure 3). The amino acid substitutions M520I and R635L are found in domains IV and V, respectively. EF-G M520 primarily interacts with V497 and mutagenesis analysis shows that the M520I amino acid substitution in domain IV has a similar interaction with V497 and both methionine (M) and isoleucine (I) are part of the aspartate family of amino acids. Interestingly, the M520I substitution occurs in the ribosomal A-site, which is essential for EF-G as well as for gentamicin. Similarly, R635 in domain V interacts with L631 and N632 via two hydrogen bonds, while the amino acid substitution R635L results in a loss of one hydrogen bond and only interacts with L631 (Figure S1, available as Supplementary data at JAC Online). Both M520I and R635L could potentially cause conformational changes in the protein structure indicated by the number of steric clashes that the mutated amino acids present (Figure S2). Furthermore, a single binding pocket was identified between domain III and domain V where ribosomal protein S12 directly interacts with EF-G, a critical element for the stabilization of EF-G while bound to the ribosome49,50 (Figure S3).

Two views of the same elongated EF-G protein structure predicted using AlphaFold. The domains II–V as well as the GTP binding domain are coloured. The amino acid substitutions causing low-level gentamicin resistance in a previous report34 are pointed out and amino acids found to be altered in the present study are in bold. This figure appears in colour in the online version of JAC and in black and white in the print version of JAC.
Figure 3.

Two views of the same elongated EF-G protein structure predicted using AlphaFold. The domains II–V as well as the GTP binding domain are coloured. The amino acid substitutions causing low-level gentamicin resistance in a previous report34 are pointed out and amino acids found to be altered in the present study are in bold. This figure appears in colour in the online version of JAC and in black and white in the print version of JAC.

Discussion

In the absence of a gonococcal vaccine, management and control of gonorrhoea rely on effective, affordable and accessible antimicrobial treatment. As ceftriaxone- and azithromycin-resistant clones are spreading internationally and the accumulation of AMR determinants in N. gonorrhoeae to these current as well as previous first-line empirical treatments for gonorrhoea continues,2–4,9–25 there is an urgent need for novel antimicrobials and/or repurposed and optimized use of existing antimicrobials. Gentamicin, recommended in combination with azithromycin as an alternative treatment for gonorrhoea,5–8,22 might become much more frequently used and verified clinical N. gonorrhoeae isolates with in vitro resistance to gentamicin (>16 mg/L) have also been mainly lacking. In fact, the in vitro susceptibility to gentamicin has remained remarkably stable in the global gonococcal population with a narrow MIC range (4–8 mg/L) over time.28–30,51 In Malawi, gentamicin (240 mg IM) plus doxycycline (100 mg two times per day for 7 days) have been used as first-line empirical treatment for urethritis and cervicitis for more than 25 years and although a recent study identified several gonorrhoea treatment failures, no gonococcal isolates with substantially elevated MICs of gentamicin were detected.29 Nevertheless, the surveillance of gonococcal susceptibility to gentamicin internationally remains scarce. Accordingly, enhanced, quality-assured surveillance of gentamicin MICs in N. gonorrhoeae is crucial internationally. Furthermore, it would be exceedingly valuable if EUCAST or CLSI recommended evidence-based clinical gentamicin resistance breakpoints for N. gonorrhoeae.

To predict gentamicin-resistance mutations that might evolve in N. gonorrhoeae if gentamicin becomes more widely used internationally, we selected for gentamicin resistance in WHO X35  in vitro using serial culture passages on gentamicin-gradient agar plates. This resulted in mutants of WHO X that expressed up to 5-fold increases in the MICs of gentamicin as well as of other aminoglycosides (Table 1). We subsequently found several mutations in the WHO X mutants with increased MICs of the aminoglycosides (WHO X10–100) and showed that high-level gentamicin resistance in N. gonorrhoeae arises in vitro in a stepwise accumulative manner (Figure 1). The most substantial increases in the gentamicin MIC were caused by mutations in fusA and ubiM encoding EF-G and UbiM involved in the biosynthetic pathway of ubiquinone (also called coenzyme Q),52 respectively. Similarly, Holley et al.34 recently showed that the A563V alteration in domain IV of EF-G conferred low-level gentamicin resistance (MIC = 32 mg/L); however, a mutation in ubiB in the gentamicin-resistant mutant was found as well.

Herein, we found that the ubiM gene mutation D186N had evolved in WHO X mutants with low-level gentamicin resistance (MIC = 32 mg/L) without any fusA mutations. Subsequently, a stepwise accumulation of mutations in several genes were observed before the final fusA (EF-GM520I) mutation caused high-level gentamicin resistance (Figure 1). The coenzyme Q is involved in aerobic respiration53 with a highly conserved pathway for biosynthesis, which involves several ubi genes54 and plays a major role in the regulatory network that allows facultative anaerobic bacteria to sense various respiratory growth conditions. Furthermore, the proton-pump NADH-ubiquinone oxidoreductase (complex I) catalyses electron transfer via proton translocation from the negative inner to the positive outer side of the membrane, therefore generating a proton-motive force (PMF). PMF is essential for the import of antibiotics55–57 and consequently, ubi mutations cause low-level resistance to aminoglycosides by suboptimal transport of the positively charged antimicrobial molecule into the cell. For example, ubi mutations have led to resistance to aminoglycosides in E. coli, specifically gentamicin and streptomycin.58,59 Similar mutations in several ubi genes may be unrecognized AMR determinants also in N. gonorrhoeae and the potential roles of coenzyme Q and ubi gene mutations need to be further investigated.

During protein synthesis, EF-G binds the ribosomal complex and facilitates the movement (translocation) of tRNAs and mRNA by one codon. EF-G domain IV binds to the 30S ribosomal A-site and to ribosomal protein S12 via domains III and V, controlling the conformation of EF-G during translocation.60 Interestingly, gentamicin binds the 30S ribosomal A-site RNA and inhibits translocation61 and most likely competes with the binding of EF-G. We propose that mutations found in EF-G in domain IV (M520I), the binding site to the 30S ribosome, is the cause of high-level resistance to gentamicin by outcompeting the antimicrobial for the A-site. However, the high-level gentamicin resistance caused SCVs of WHO X, as previously reported in gentamicin-resistant S. aureus,62 and competitive coculture in an HFIM for gonorrhoea showed that these SCVs also had impaired biofitness compared with WHO X.

In conclusion, we describe the first high-level gentamicin-resistant N. gonorrhoeae isolate (MIC = 128 mg/L), which was selected in vitro through experimental evolution. The most substantial increases in the MICs of gentamicin and other aminoglycosides such as kanamycin, tobramycin and amikacin were caused by resistance mutations in fusA (G1560A and G1904T encoding EF-G M520I and R635L, respectively) and ubiM (D186N). However, the high-level gentamicin-resistant N. gonorrhoeae mutant showed impaired biofitness. Additional studies are required to confirm the role of all the mutations potentially involved in gentamicin resistance and to examine whether N. gonorrhoeae could restore biofitness through compensatory genetic events during continued exposure to gentamicin.

Funding

This project was supported by the Örebro County Council Research Committee and the Foundation for Medical Research at Örebro University Hospital, Örebro, Sweden. NIH grant AI47609 supports W.M.S. W.M.S. is the recipient of a Senior Research Career Scientist award from the Medical Research Service of the Department of Veterans Affairs.

Transparency declarations

None to declare. The contents of this article are solely the responsibility of the authors and do not necessarily reflect the official views of the National Institutes of Health, the Department of Veterans Affairs or the United States government.

Supplementary data

Figures S1–S3 are available as Supplementary data at JAC Online.

References

1

WHO
.
Global progress report on HIV, viral hepatitis and sexually transmitted infections
,
2021
. https://www.who.int/publications/i/item/9789240027077.

2

Unemo
 
M
,
Shafer
 
WM
.
Antimicrobial resistance in Neisseria gonorrhoeae in the 21st century: past, evolution, and future
.
Clin Microbiol Rev
 
2014
;
27
:
587
613
. https://doi.org/10.1128/CMR.00010-14

3

Golparian
 
D
,
Harris
 
SR
,
Sánchez-Busó
 
L
 et al.  
Genomic evolution of Neisseria gonorrhoeae since the preantibiotic era (1928–2013): antimicrobial use/misuse selects for resistance and drives evolution
.
BMC Genomics
 
2020
;
21
:
116
. https://doi.org/10.1186/s12864-020-6511-6

4

Unemo
 
M
,
Seifert
 
HS
,
Hook
 
EW
 III
 et al.  
Gonorrhoea
.
Nat Rev Dis Primers
 
2019
;
5
:
79
. https://doi.org/10.1038/s41572-019-0128-6

5

Unemo
 
M
,
Ross
 
JDC
,
Serwin
 
AB
 et al.  
2020 European guideline for the diagnosis and treatment of gonorrhoea in adults
.
Int J STD AIDS
 
2020
;
956462420949126
. https://doi.org/10.1177/0956462420949126

6

WHO
.
WHO Guidelines for the Treatment of Neisseria gonorrhoeae
.
2016
. https://apps.who.int/iris/bitstream/handle/10665/246114/9789241549691-eng.pdf.

7

Workowski
 
KA
,
Bachmann
 
LH
,
Chan
 
PA
 et al.  
Sexually transmitted infections treatment guidelines, 2021
.
MMWR Recomm Rep
 
2021
;
70
:
1
187
. https://doi.org/10.15585/mmwr.rr7004a1

8

Fifer
 
H
,
Saunders
 
J
,
Soni
 
S
 et al.  
2018 UK national guideline for the management of infection with Neisseria gonorrhoeae
.
Int J STD AIDS
 
2020
;
31
:
4
15
. https://doi.org/10.1177/0956462419886775

9

Eyre
 
DW
,
Sanderson
 
ND
,
Lord
 
E
 et al.  
Gonorrhoea treatment failure caused by a Neisseria gonorrhoeae strain with combined ceftriaxone and high-level azithromycin resistance, England, February 2018
.
Euro Surveill
 
2018
;
23
:
1800323
. https://doi.org/10.2807/1560-7917.ES.2018.23.27.1800323

10

Whiley
 
DM
,
Jennison
 
A
,
Pearson
 
J
 et al.  
Genetic characterisation of Neisseria gonorrhoeae resistant to both ceftriaxone and azithromycin
.
Lancet Infect Dis
 
2018
;
18
:
717
8
. https://doi.org/10.1016/S1473-3099(18)30340-2

11

Pleininger
 
S
,
Indra
 
A
,
Golparian
 
D
 et al.  
Extensively drug-resistant (XDR) Neisseria gonorrhoeae causing possible gonorrhoea treatment failure with ceftriaxone plus azithromycin in Austria, April 2022
.
Euro Surveill
 
2022
;
27
:
2200455
. https://doi.org/10.2807/1560-7917.ES.2022.27.24.2200455

12

Nakayama
 
S-I
,
Shimuta
 
K
,
Furubayashi
 
K-I
 et al.  
New ceftriaxone- and multidrug-resistant Neisseria gonorrhoeae strain with a novel mosaic penA gene isolated in Japan
.
Antimicrob Agents Chemother
 
2016
;
60
:
4339
41
. https://doi.org/10.1128/AAC.00504-16

13

Lahra
 
MM
,
Martin
 
I
,
Demczuk
 
W
 et al.  
Cooperative recognition of internationally disseminated ceftriaxone-resistant Neisseria gonorrhoeae strain
.
Emerg Infect Dis
 
2018
;
24
:
735
40
. https://doi.org/10.3201/eid2404.171873

14

Eyre
 
DW
,
Town
 
K
,
Street
 
T
 et al.  
Detection in the United Kingdom of the Neisseria gonorrhoeae FC428 clone, with ceftriaxone resistance and intermediate resistance to azithromycin, October to December 2018
.
Euro Surveill
 
2019
;
24
:
1900147
. https://doi.org/10.2807/1560-7917.ES.2019.24.10.1900147

15

Day
 
M
,
Pitt
 
R
,
Mody
 
N
 et al.  
Detection of 10 cases of ceftriaxone-resistant Neisseria gonorrhoeae in the United Kingdom, December 2021 to June 2022
.
Euro Surveill
 
2022
;
27
:
2200803
. https://doi.org/10.2807/1560-7917.ES.2022.27.46.2200803

16

Golparian
 
D
,
Rose
 
L
,
Lynam
 
A
 et al.  
Multidrug-resistant Neisseria gonorrhoeae isolate, belonging to the internationally spreading Japanese FC428 clone, with ceftriaxone resistance and intermediate resistance to azithromycin, Ireland, August 2018
.
Euro Surveill
 
2018
;
23
:
1800617
. https://doi.org/10.2807/1560-7917.ES.2018.23.47.1800617

17

Lee
 
K
,
Nakayama
 
S-I
,
Osawa
 
K
 et al.  
Clonal expansion and spread of the ceftriaxone-resistant Neisseria gonorrhoeae strain FC428, identified in Japan in 2015, and closely related isolates
.
J Antimicrob Chemother
 
2019
;
74
:
1812
9
. https://doi.org/10.1093/jac/dkz129

18

Trinh
 
TM
,
Nguyen
 
TT
,
Le
 
TV
 et al.  
Neisseria gonorrhoeae FC428 subclone, Vietnam, 2019–2020
.
Emerg Infect Dis
 
2022
;
28
:
432
5
. https://doi.org/10.3201/eid2802.211788

19

Lin
 
X
,
Chen
 
W
,
Xie
 
Q
 et al.  
Dissemination and genome analysis of high-level ceftriaxone-resistant penA60.001 Neisseria gonorrhoeae strains from the Guangdong Gonococcal Antibiotics Susceptibility Programme (GD-GASP), 2016-2019
.
Emerg Microbes Infect
 
2022
;
11
:
344
50
. https://doi.org/10.1080/22221751.2021.2011618

20

Terkelsen
 
D
,
Tolstrup
 
J
,
Johnsen
 
CH
 et al.  
Multidrug-resistant Neisseria gonorrhoeae infection with ceftriaxone resistance and intermediate resistance to azithromycin, Denmark, 2017
.
Euro Surveill
 
2017
;
22
:
17-00659
. https://doi.org/10.2807/1560-7917.ES.2017.22.42.17-00659

21

Berenger
 
BM
,
Demczuk
 
W
,
Gratrix
 
J
 et al.  
Genetic characterization and enhanced surveillance of ceftriaxone-resistant Neisseria gonorrhoeae strain, Alberta, Canada, 2018
.
Emerg Infect Dis
 
2019
;
25
:
1660
7
. https://doi.org/10.3201/eid2509.190407

22

Unemo
 
M
,
Lahra
 
MM
,
Escher
 
M
 et al.  
WHO global antimicrobial resistance surveillance for Neisseria gonorrhoeae 2017–18: a retrospective observational study
.
Lancet Microbe
 
2021
;
2
:
e627
36
. https://doi.org/10.1016/S2666-5247(21)00171-3

23

Sánchez-Busó
 
L
,
Cole
 
MJ
,
Spiteri
 
G
 et al.  
Europe-wide expansion and eradication of multidrug-resistant Neisseria gonorrhoeae lineages: a genomic surveillance study
.
Lancet Microbe
 
2022
;
3
:
e452
63
. https://doi.org/10.1016/S2666-5247(22)00044-1

24

Gernert
 
KM
,
Seby
 
S
,
Schmerer
 
MW
 et al.  
Azithromycin susceptibility of Neisseria gonorrhoeae in the USA in 2017: a genomic analysis of surveillance data
.
Lancet Microbe
 
2020
;
1
:
e154
64
. https://doi.org/10.1016/S2666-5247(20)30059-8

25

Lin
 
X
,
Qin
 
X
,
Wu
 
X
 et al.  
Markedly increasing antibiotic resistance and dual treatment of Neisseria gonorrhoeae isolates in Guangdong, China, from 2013 to 2020
.
Antimicrob Agents Chemother
 
2022
;
66
:
e0229421
. https://doi.org/10.1128/aac.02294-21

27

Kanj
 
SS
,
Bassetti
 
M
,
Kiratisin
 
P
 et al.  
Clinical data from studies involving novel antibiotics to treat multidrug-resistant Gram-negative bacterial infections
.
Int J Antimicrob Agents
 
2022
;
60
:
106633
. https://doi.org/10.1016/j.ijantimicag.2022.106633

28

Brown
 
LB
,
Krysiak
 
R
,
Kamanga
 
G
 et al.  
Neisseria gonorrhoeae antimicrobial susceptibility in Lilongwe, Malawi, 2007
.
Sex Transm Dis
 
2010
;
37
:
169
72
. https://doi.org/10.1097/OLQ.0b013e3181bf575c

29

Matoga
 
M
,
Chen
 
JS
,
Krysiak
 
R
 et al.  
Gentamicin susceptibility in Neisseria gonorrhoeae and treatment outcomes for urogenital gonorrhea after 25 years of sustained gentamicin use in Malawi
.
Sex Transm Dis
 
2022
;
49
:
251
6
. https://doi.org/10.1097/OLQ.0000000000001580

30

Chisholm
 
SA
,
Quaye
 
N
,
Cole
 
MJ
 et al.  
An evaluation of gentamicin susceptibility of Neisseria gonorrhoeae isolates in Europe
.
J Antimicrob Chemother
 
2011
;
66
:
592
5
. https://doi.org/10.1093/jac/dkq476

31

Shaw
 
KJ
,
Rather
 
PN
,
Hare
 
RS
 et al.  
Molecular genetics of aminoglycoside resistance genes and familial relationships of the aminoglycoside-modifying enzymes
.
Microbiol Rev
 
1993
;
57
:
138
63
. https://doi.org/10.1128/mr.57.1.138-163.1993

32

Mingeot-Leclercq
 
MP
,
Glupczynski
 
Y
,
Tulkens
 
PM
.
Aminoglycosides: activity and resistance
.
Antimicrob Agents Chemother
 
1999
;
43
:
727
37
. https://doi.org/10.1128/AAC.43.4.727

33

Maness
 
MJ
,
Foster
 
GC
,
Sparling
 
PF
.
Ribosomal resistance to streptomycin and spectinomycin in Neisseria gonorrhoeae
.
J Bacteriol
 
1974
;
120
:
1293
9
. https://doi.org/10.1128/jb.120.3.1293-1299.1974

34

Holley
 
CL
,
Dhulipala
 
V
,
Balthazar
 
JT
 et al.  
A single amino acid substitution in elongation factor G can confer low-level gentamicin resistance in Neisseria gonorrhoeae
.
Antimicrob Agents Chemother
 
2022
;
66
:
e0025122
. https://doi.org/10.1128/aac.00251-22

35

Unemo
 
M
,
Golparian
 
D
,
Sánchez-Busó
 
L
 et al.  
The novel 2016 WHO Neisseria gonorrhoeae reference strains for global quality assurance of laboratory investigations: phenotypic, genetic and reference genome characterization
.
J Antimicrob Chemother
 
2016
;
71
:
3096
108
. https://doi.org/10.1093/jac/dkw288

36

Jacobsson
 
S
,
Golparian
 
D
,
Oxelbark
 
J
 et al.  
Pharmacodynamic evaluation of lefamulin in the treatment of gonorrhea using a hollow fiber infection model simulating Neisseria gonorrhoeae infections
.
Front Pharmacol
 
2022
;
13
:
1035841
. https://doi.org/10.3389/fphar.2022.1035841

37

Berglund
 
T
,
Unemo
 
M
,
Olcén
 
P
 et al.  
One year of Neisseria gonorrhoeae isolates in Sweden: the prevalence study of antibiotic susceptibility shows relation to the geographic area of exposure
.
Int J STD AIDS
 
2002
;
13
:
109
14
. https://doi.org/10.1258/0956462021924730

38

CLSI
.
Performance Standards for Antimicrobial Susceptibility Testing—Thirty-Second Edition: M100
.
2022
.

39

Weinberg
 
ED
.
Gradient agar plates
.
Am Biol Teach
 
1959
;
21
:
347
50
. https://doi.org/10.2307/4439185

40

Tamura
 
K
,
Stecher
 
G
,
Kumar
 
S
.
MEGA11: molecular evolutionary genetics analysis version 11
.
Mol Biol Evol
 
2021
;
38
:
3022
7
. https://doi.org/10.1093/molbev/msab120

41

Argimón
 
S
,
Abudahab
 
K
,
Goater
 
RJE
 et al.  
Microreact: visualizing and sharing data for genomic epidemiology and phylogeography
.
Microb Genom
 
2016
;
2
:
e000093
.

42

Golparian
 
D
,
Jacobsson
 
S
,
Sánchez-Busó
 
L
 et al.  
GyrB in silico mining in 27 151 global gonococcal genomes from 1928–2021 combined with zoliflodacin in vitro testing of 71 international gonococcal isolates with different GyrB, ParC and ParE substitutions confirms high susceptibility
.
J Antimicrob Chemother
 
2022
;
78
:
150
4
. https://doi.org/10.1093/jac/dkac366

43

Jacobsson
 
S
,
Golparian
 
D
,
Oxelbark
 
J
 et al.  
Pharmacodynamic evaluation of zoliflodacin treatment of Neisseria gonorrhoeae strains with amino acid substitutions in the zoliflodacin target GyrB using a dynamic hollow fiber infection model
.
Front Pharmacol
 
2022
;
13
:
874176
. https://doi.org/10.3389/fphar.2022.874176

44

Vincent
 
LR
,
Kerr
 
SR
,
Tan
 
Y
 et al.  
In vivo-selected compensatory mutations restore the fitness cost of mosaic penA alleles that confer ceftriaxone resistance in Neisseria gonorrhoeae
.
mBio
 
2018
;
9
:
e01905-17
. https://doi.org/10.1128/mBio.01905-17

45

Jumper
 
J
,
Evans
 
R
,
Pritzel
 
A
 et al.  
Highly accurate protein structure prediction with AlphaFold
.
Nature
 
2021
;
596
:
583
9
. https://doi.org/10.1038/s41586-021-03819-2

46

Mistry
 
J
,
Chuguransky
 
S
,
Williams
 
L
 et al.  
Pfam: the protein families database in 2021
.
Nucleic Acids Res
 
2021
;
49
:
D412
9
. https://doi.org/10.1093/nar/gkaa913

47

Tian
 
W
,
Chen
 
C
,
Lei
 
X
 et al.  
CASTp 3.0: computed atlas of surface topography of proteins
.
Nucleic Acids Res
 
2018
;
46
:
W363
7
. https://doi.org/10.1093/nar/gky473

48

Demczuk
 
W
,
Lynch
 
T
,
Martin
 
I
 et al.  
Whole-genome phylogenomic heterogeneity of Neisseria gonorrhoeae isolates with decreased cephalosporin susceptibility collected in Canada between 1989 and 2013
.
J Clin Microbiol
 
2015
;
53
:
191
200
. https://doi.org/10.1128/JCM.02589-14

49

Panecka
 
J
,
Mura
 
C
,
Trylska
 
J
.
Interplay of the bacterial ribosomal A-site, S12 protein mutations and paromomycin binding: a molecular dynamics study
.
PLoS One
 
2014
;
9
:
e111811
. https://doi.org/10.1371/journal.pone.0111811

50

Macé
 
K
,
Giudice
 
E
,
Chat
 
S
 et al.  
The structure of an elongation factor G-ribosome complex captured in the absence of inhibitors
.
Nucleic Acids Res
 
2018
;
46
:
3211
7
. https://doi.org/10.1093/nar/gky081

51

Mann
 
LM
,
Kirkcaldy
 
RD
,
Papp
 
JR
 et al.  
Susceptibility of Neisseria gonorrhoeae to gentamicin—gonococcal isolate surveillance project, 2015–2016
.
Sex Transm Dis
 
2018
;
45
:
96
8
. https://doi.org/10.1097/OLQ.0000000000000693

52

Wong
 
K
,
Gehring
 
K
.
Deciphering the catalytic mechanism of bacterial ubiquitination
.
Nature
 
2018
;
557
:
644
5
. https://doi.org/10.1038/d41586-018-05250-6

53

Schoepp-Cothenet
 
B
,
Lieutaud
 
C
,
Baymann
 
F
 et al.  
Menaquinone as pool quinone in a purple bacterium
.
Proc Nat Acad Sci U S A
 
2009
;
106
:
8549
54
. https://doi.org/10.1073/pnas.0813173106

54

Meganathan
 
R
.
Ubiquinone biosynthesis in microorganisms
.
FEMS Microbiol Lett
 
2001
;
203
:
131
9
. https://doi.org/10.1111/j.1574-6968.2001.tb10831.x

55

Taber
 
HW
,
Mueller
 
JP
,
Miller
 
PF
 et al.  
Bacterial uptake of aminoglycoside antibiotics
.
Microbiol Rev
 
1987
;
51
:
439
57
. https://doi.org/10.1128/mr.51.4.439-457.1987

56

Davis
 
BD
,
Chen
 
LL
,
Tai
 
PC
.
Misread protein creates membrane channels: an essential step in the bactericidal action of aminoglycosides
.
Proc Nat Acad Sci U S A
 
1986
;
83
:
6164
8
. https://doi.org/10.1073/pnas.83.16.6164

57

Ezraty
 
B
,
Vergnes
 
A
,
Banzhaf
 
M
 et al.  
Fe–S cluster biosynthesis controls uptake of aminoglycosides in a ROS-less death pathway
.
Science
 
2013
;
340
:
1583
7
. https://doi.org/10.1126/science.1238328

58

Collis
 
CM
,
Grigg
 
GW
.
An Escherichia coli mutant resistant to phleomycin, bleomycin, and heat inactivation is defective in ubiquinone synthesis
.
J Bacteriol
 
1989
;
171
:
4792
8
. https://doi.org/10.1128/jb.171.9.4792-4798.1989

59

Bryan
 
LE
,
Van Den Elzen
 
HM
.
Effects of membrane-energy mutations and cations on streptomycin and gentamicin accumulation by bacteria: a model for entry of streptomycin and gentamicin in susceptible and resistant bacteria
.
Antimicrob Agents Chemother
 
1977
;
12
:
163
77
. https://doi.org/10.1128/AAC.12.2.163

60

Lin
 
J
,
Gagnon
 
MG
,
Bulkley
 
D
 et al.  
Conformational changes of elongation factor G on the ribosome during tRNA translocation
.
Cell
 
2015
;
160
:
219
27
. https://doi.org/10.1016/j.cell.2014.11.049

61

Fourmy
 
D
,
Recht
 
MI
,
Blanchard
 
SC
 et al.  
Structure of the A site of Escherichia coli 16S ribosomal RNA complexed with an aminoglycoside antibiotic
.
Science
 
1996
;
274
:
1367
71
. https://doi.org/10.1126/science.274.5291.1367

62

Lannergård
 
J
,
von Eiff
 
C
,
Sander
 
G
 et al.  
Identification of the genetic basis for clinical menadione-auxotrophic small-colony variant isolates of Staphylococcus aureus
.
Antimicrob Agents Chemother
 
2008
;
52
:
4017
22
. https://doi.org/10.1128/AAC.00668-08

This article is published and distributed under the terms of the Oxford University Press, Standard Journals Publication Model (https://dbpia.nl.go.kr/pages/standard-publication-reuse-rights)

Supplementary data