Molecular Characteristics and the Effect of Mutations in Different Sites of the rplD Gene Among Clinical Isolates of Azithromycin Resistance Neisseria gonorrhoeae in Eastern China

authors:

avatar WangQiang Hu ORCID 1 , avatar Han Chen 1 , avatar Dan Zhang ORCID 1 , avatar Fangyou Yu 2 , avatar Mingpeng Hu 2 , 3 , * , avatar Zhou Zheng ORCID 1 , **

Department of Clinical Laboratory, Key Laboratory of Clinical Laboratory Diagnosis and Translational Research of Zhejiang Province, the First Affiliated Hospital of Wenzhou Medical University, Wenzhou, China
Department of Laboratory Medicine, Shanghai Pulmonary Hospital, Tongji University School of Medicine, Shanghai, China
Ruian Hospital of Traditional Chinese Medicine, Wenzhou, China
Corresponding Authors:

How To Cite Hu W, Chen H, Zhang D, Yu F, Hu M, et al. Molecular Characteristics and the Effect of Mutations in Different Sites of the rplD Gene Among Clinical Isolates of Azithromycin Resistance Neisseria gonorrhoeae in Eastern China. Jundishapur J Microbiol. 2023;16(9):e140659. https://doi.org/10.5812/jjm-140659.

Abstract

Background:

The emergence of resistance to azithromycin complicates the treatment of Neisseria gonorrhoeae, the causative agent of gonorrhea.

Objectives:

The main objective of this study was to investigate the impact of mutations in different sites of the rplD gene on azithromycin resistance (AZM-R) and the molecular characteristics of N. gonorrhoeae. This study aimed to provide evidence for drug resistance and transmission.

Methods:

A total of 37 isolates of N. gonorrhoeae were collected within January 2018 to December 2020. These isolates were obtained from urine, semen, or vaginal secretions of different patients. Azithromycin resistance was assessed, and genes associated with AZM-R, namely rplD, and rplV, were analyzed using polymerase chain reaction (PCR) and deoxyribonucleic acid (DNA) sequencing. All clinical isolates were characterized through multi-locus sequence typing (MLST).

Results:

The study comprised 21 azithromycin-resistant N. gonorrhoeae isolates, with two of them demonstrating high resistance, indicated by a minimum inhibitory concentration (MIC) greater than 256 µg/mL. Additionally, 16 sensitive isolates were included in the study. Ten isolates were observed to have rplD point mutations, including mutations such as G70D, G70S, G68D, and A43T. No mutations were detected in rplV. The rate of point mutations in rplD was significantly different between the azithromycin-sensitive (AZM-S) group and the AZM-R group (P < 0.05). Among the 37 isolates studied, 12 distinct MLST types were identified and further grouped into four different MLST goeBURST groups. The two isolates with high-level AZM-R were ST1901 and ST1588, respectively.

Conclusions:

The clinical isolates of N. gonorrhoeae from Wenzhou, Eastern China, exhibit significant genetic diversity and a relatively high prevalence of AZM-R. Mutations in the rplD gene were identified, which reduced susceptibility to macrolides and were significantly associated with increased AZM-R.

1. Background

The increasing antibiotic resistance of Neisseria gonorrhoeae, the pathogen responsible for gonorrhea, is a growing global concern. According to the US Centers for Disease Control and Prevention (CDC), the incidence rate of gonorrhea increased by 111% from 2009 to 2020 (1). Gonorrhea is one of the most common bacterial sexually transmitted infections (STIs), with an estimated 82,400,000 cases worldwide in 2020, as reported by the World Health Organization (WHO) (2). The current empirical treatment for gonorrhea in the United States consists of a combination of ceftriaxone and azithromycin (3). However, the emergence of N. gonorrhoeae isolates with high-level azithromycin resistance (AZM-R) poses a threat to the sustainability of this treatment (4, 5). The rapid emergence of AZM-R has led some countries, such as the United Kingdom, to consider this treatment endangered and recommend ceftriaxone monotherapy instead (3).

Azithromycin is a 15-membered macrolide that inhibits protein synthesis by binding to the bacterial 50S large ribosome. Azithromycin functions by inhibiting the peptide transfer/translocation step, which blocks the 50S peptide exit tunnel, leading to the release of incomplete peptides from the ribosomes (Figure 1) (6-9). Ribosomal proteins L4 and L22, products of the rplD and rplV genes, respectively, bind to domain I of 23S rRNA to form a channel and facilitate the entry of macrolide antibiotics (10). Mutated L4 and L22 proteins can alter domains II, III, and V, affecting the susceptibility of microorganisms to macrolides (11).

The ribosome nascent chain tunnel environment. Staphylococcus aureus ribosome, where the large subunit is shown in light gray, and the small subunit is shown in dark gray. The A-site, P-site, and E-site docked tRNA molecules are shown in blue, green, and magenta, respectively. The surface of a nascent chain within the peptide exit tunnel is shown in yellow. Bound macrolide is shown in red, and uL4, uL22, uL23, and uL24 are shown in brown, orange, teal, and khaki, respectively. The figure was reproduced with permission from Halfon et al. (9)
The ribosome nascent chain tunnel environment. Staphylococcus aureus ribosome, where the large subunit is shown in light gray, and the small subunit is shown in dark gray. The A-site, P-site, and E-site docked tRNA molecules are shown in blue, green, and magenta, respectively. The surface of a nascent chain within the peptide exit tunnel is shown in yellow. Bound macrolide is shown in red, and uL4, uL22, uL23, and uL24 are shown in brown, orange, teal, and khaki, respectively. The figure was reproduced with permission from Halfon et al. (9)

Mutations in L4 and L22 have been confirmed to confer resistance to macrolides in various bacteria, such as Brucella, Escherichia coli, and Streptococcus pneumoniae (12-14). Although no mutation related to the 50S ribosome protein L22 has been reported in N. gonorrhoeae, point mutations of the 50S ribosome protein L4, specifically the G70D mutation, have been implicated in macrolide resistance (15, 16). Additional studies have identified other L4 mutations at amino acid positions 68 (G68D and G68C), 69 (T69I), and 70 (G70R, G70S, and G70A), which are associated with decreased susceptibility to azithromycin (16). These mutations are located at the end of the L4 ring near the macrolide binding site and contribute significantly to the development of AZM-R (17).

Multi-locus sequence typing (MLST), a method that characterizes isolates by amplifying sequences from seven housekeeping genes, has been widely used to study the dynamics of N. gonorrhoeae infection and drug-resistant isolates. Multi-locus sequence typing is frequently employed to trace the spread of N. gonorrhoeae isolates, identify lineages, and detect short-term transmission chains (18, 19). It is crucial to deepen our understanding of the emergence and dynamics of drug-resistant N. gonorrhoeae isolates at the national and global levels.

2. Objectives

The objective of this study was to examine the mutation rate and diversity of the rplD gene in N. gonorrhoeae within the Wenzhou area and to assess the extent to which rplD mutations contribute to the consistently elevated azithromycin minimum inhibitory concentration (MIC) levels. Additionally, this study aimed to understand the genetic characteristics of N. gonorrhoeae and analyze the variations in AZM-R among isolates with different genetic backgrounds. Ultimately, the results of this study will serve as a valuable foundation for the prevention and treatment of gonorrhea.

3. Methods

3.1. Bacterial Isolates

This study was conducted within January 2018 and December 2020, during which 37 non-repetitive N. gonorrhoeae isolates were randomly selected from the First Affiliated Hospital of Wenzhou Medical University, Wenzhou, Zhejiang, China. The isolates originated from 5 females and 32 males, aged 18 - 66 years (median age: 31 years), with urine and secretions being the main types of specimens. All experimental isolates were confirmed as N. gonorrhoeae through the fully automated rapid microbial mass spectrometry detection system (VITEK MS, from bioMérieux, France). These isolates included 21 AZM-R isolates and 16 azithromycin-sensitive (AZM-S) isolates.

3.2. Azithromycin Susceptibility Testing

The susceptibility of N. gonorrhoeae to azithromycin was evaluated through MIC determination (Wenzhou Kangtai Biotechnology Co., Ltd., China). All test results were interpreted following the guidelines of the Clinical and Laboratory Standards Institute (CLSI M100ed30) for all antibiotics (20), where MIC ≤ 1 mg/L was deemed sensitive, MIC > 1 mg/L was considered resistant, and MIC ≥ 256 mg/L indicated high-level drug resistance. Individual N. gonorrhoeae isolates were isolated and purified on Thayer-Martin (T-M) selective media. Azithromycin E-test experiments were conducted using GC plates supplemented with 1% IsoVitaleX, vancomycin, colistin, nystatin, and trimethoprim selective supplements in a 5% CO2 incubator at 35°C. This helped determine either the area diameter or MIC. The quality control strain was N. gonorrhoeae ATCC 49226 from the clinical laboratory of the Ministry of Health of the People’s Republic of China. All isolates were stored in glycerol broth at -80°C.

3.3. Extraction of Bacterial Genomic Deoxyribonucleic Acid

The genomic deoxyribonucleic acid (gDNA) of N. gonorrhoeae was extracted through boiling. Moreover, 3 - 5 purified N. gonorrhoeae colonies were suspended in 200 μL of sterile distilled water, boiled in water at 100°C for 10 minutes, and then centrifuged at 13800 g at 4°C for 15 minutes (10). The resulting supernatant contained gDNA, which was stored at -20°C for backup purposes.

3.4. Amplification of the rplD and rplV Genes

The primers for rplD and rplV have been previously published and are listed in Table 1 (10). The polymerase chain reaction (PCR) conditions for rplD involved a pre-denaturation step at 94°C for 4 minutes, followed by denaturation at 94°C for 30 seconds, annealing at 60°C for 30 seconds, and extension at 72°C for 45 seconds, with a total of 30 cycles. For rplV amplification, the annealing temperature was adjusted to 55°C, and the extension time was increased to 1 minute; nevertheless, the other conditions remained unchanged (21). The amplified PCR products were subjected to electrophoresis at 110 V for 30 minutes, photographed with a gel imager, and confirmed as single bands. Subsequently, the PCR products were purified using magnetic beads by the sequencing company (Beijing Qingke New Industry Biotechnology Co., Ltd. Hangzhou Branch, China). First-generation sequencing technology was employed for sequencing. The sequencing results were compared to DNA sequences using BLAST to identify genetic mutations.

Table 1.

The rplD and rplV Gene Primer Sequences and Target Fragment Sizes

DesignationPCR Primer Sequence (5ˊ→3ˊ)Product Size (bp)Annealing Temperature (°C)
rplDFCAGCGATGTTGTAGTTCGT48260
RACTCAAGTAATCTTGGCGC
rplVFTCAGCGACAATATGGTTGGT46855
RAGCCCAGTCTTTAGTTACC

3.5. Neisseria gonorrhoeae MLST

Neisseria gonorrhoeae isolates in the present study were genotyped using MLST. Seven housekeeping genes, namely adk, abcZ, aroE, fumC, pgm, gdh, and pdhC of N. gonorrhoeae, were amplified and sequenced. Sequence types (STs) were identified using the online database at the Pasteur Institute MLST typing website for N. gonorrhoeae. The evolutionary relationship between the isolates was analyzed using the platform-independent Java software PHYLOViZ at the single-locus variant (SLV) level with the goeBURST algorithm (22).

3.6. Statistical Analysis

SPSS software (version 22.0) was used to statistically analyze the difference in rplD point mutations between the AZM-S and AZM-R groups of N. gonorrhoeae, with a P-value less than 0.05 indicating significance.

4. Results

4.1. Antimicrobial Susceptibility

Among the 37 isolates of N. gonorrhoeae, the MIC values of azithromycin ranged from 0.125 to 256 μg/mL. Sixteen isolates were classified as sensitive, accounting for 43.24% of the total; however, 21 isolates were resistant, making up 56.76%. Two isolates exhibited extremely high MIC values exceeding 256 μg/mL, categorizing them as high-level AZM-R N. gonorrhoeae.

4.2. Comparison of rplD and rplV Genes

The current study analyzed the DNA of N. gonorrhoeae through PCR amplification and identified 10 isolates with rplD point mutations, resulting in a mutation rate of 27.03%. These mutations were found at four specific sites: G70D, G70S, G68D, and A43T (Table 2). It is worth noting that the A43T mutation was detected for the first time in this study. Except for A43T, the remaining mutations were confirmed to be non-synonymous and were located at the end of the L4 loop, which is near the macrolide-binding site. All these mutations were associated with reduced susceptibility to azithromycin, as suggested by previous studies (5, 14, 16).

Table 2.

Characterization and Molecular Typing Results of Clinical Isolates of Neisseria gonorrhoeae

GroupsStrain No.AZM MIC (µg/mL)rplD Mutation (L4 Protein)rplV Mutation (L22 Protein)MLST ST
AZM-R (n = 21)24G70DNo mutation9899
38No mutationNo mutation8123
52No mutationNo mutation7363
64No mutationNo mutation1901
91.25No mutationNo mutation1928
1316G70SNo mutation7367
1616G70SNo mutation7367
211.25No mutationNo mutation7367
25> 256No mutationNo mutation1901
272G70SNo mutation7367
294G70DNo mutation14421
311.25No mutationNo mutation7822
322G68DNo mutation7822
341.25No mutationNo mutation1901
37> 256A43TNo mutation1588
382No mutationNo mutation1928
391.25No mutationNo mutation14283
414G70DNo mutation15251
421.25No mutationNo mutation8123
432G68DNo mutation7822
448No mutationNo mutation8123
AZM-S (n = 16)10.25No mutationNo mutation7822
81No mutationNo mutation8123
120.5No mutationNo mutation7822
170.125No mutationNo mutation1928
181No mutationNo mutation7363
190.5No mutationNo mutation1928
201No mutationNo mutation1901
220.125No mutationNo mutation1588
230.032No mutationNo mutation1600
240.125No mutationNo mutation1588
260.5No mutationNo mutation8123
281No mutationNo mutation1928
300.125No mutationNo mutation8123
330.25G70SNo mutation7367
360.25No mutationNo mutation1600
400.5No mutationNo mutation8123
P0.015

When comparing the AZM-R and AZM-S groups based on whether the MIC value of azithromycin exceeded 1 μg/mL, it was observed that the point mutation rate of rplD in the AZM-S group was significantly lower than in the AZM-R group (P < 0.05). This indicates that the presence of point mutations in rplD was more common in isolates resistant to azithromycin than those sensitive to azithromycin. Another noteworthy discovery was that one of the two high-level AZM-R isolates exhibited an rplD A43T mutation; nonetheless, no mutations were detected in the other. Lastly, the study did not identify any rplV mutations in the analyzed isolates.

4.3. Molecular Epidemiologic Typing

Among the 37 clinical isolates of N. gonorrhoeae, MLST identified 12 different STs. The distribution of these STs was as follows: ST8123 accounted for 18.91% (7/37), ST7822 for 13.51% (5/37), ST7367 for 13.51% (5/37), and ST1928 for 13.51% (5/37). The STs of two high-level AZM-R isolates were identified as ST1901 and ST1588. Using goeBURST minimum spanning tree analysis on the MLST allelic spectrum of the seven MLST housekeeping genes of N. gonorrhoeae, the genetic relatedness of the MLSTs in this experiment was compared to all the MLSTs of N. gonorrhoeae (Figure 2). Based on these genetic relationships, the MLSTs identified in this study were categorized into four large clusters, as shown in Figure 3. Cluster A isolates consisted of 6 different ST types centered around ST7363, with ST14283, ST1600, and ST1823 representing single locus variants of ST7363. Cluster B isolates comprised four different ST types centered around ST1901, with ST7822 being a single locus variant of ST1901. Cluster C and D isolates each contained only one type.

Minimum spanning tree for all Neisseria gonorrhoeae multi-locus sequence typing (MLST) sequence types (STs). The graph was analyzed by PHYLOVIZ using the goeBURST algorithm at the single-locus variant (SLV) level. The types of this study are marked in red dots.
Minimum spanning tree for all Neisseria gonorrhoeae multi-locus sequence typing (MLST) sequence types (STs). The graph was analyzed by PHYLOVIZ using the goeBURST algorithm at the single-locus variant (SLV) level. The types of this study are marked in red dots.
Genetic relatedness of multi-locus sequence typing (MLST) determined by goeBURST minimum spanning tree analysis using the Neisseria gonorrhoeae MLST allelic profiles of seven MLST housekeeping genes
Genetic relatedness of multi-locus sequence typing (MLST) determined by goeBURST minimum spanning tree analysis using the Neisseria gonorrhoeae MLST allelic profiles of seven MLST housekeeping genes

The two isolates highly resistant to azithromycin, ST1588 and ST1901, differed at the four loci adk, aroE, fumC, and gdh. They belonged to cluster A and cluster B, respectively, indicating that they originated from distinct genetic backgrounds, as shown in Figure 2. This suggests that high-level AZM-R N. gonorrhoeae can emerge from different genetic backgrounds, potentially necessitating different treatment strategies.

5. Discussion

Over the years, gonorrhea has developed resistance to all classes of antimicrobials targeting it, especially with the emergence of isolates with reduced susceptibility to azithromycin and ceftriaxone, posing a public health concern due to the rise of incurable gonorrhea infections (23). To prevent further drug resistance, a combination of azithromycin and ceftriaxone has been used for the treatment of gonorrhea (23, 24). However, despite combination therapy, azithromycin resistance has continued to increase, and the emergence of N. gonorrhoeae resistant to both azithromycin and ceftriaxone since 2018 is of even greater concern (25). In 2001, N. gonorrhoeae with high-level resistance to azithromycin (MIC ≥ 256 mg/L) was first isolated in Argentina, and it has been constantly reported in various countries, including England, Italy, the United States, Canada, and China (26). The G70D mutation was the most frequently found point mutation in rplD, which is significantly related to increased AZM-R (16). A recent bacterial genome-wide association study confirmed the role of the rplD G70D mutation in mediating macrolide resistance, showing an average six-fold increase in clarithromycin MICs, with the MICs of azithromycin and erythromycin being 4 times higher than that of the wild-type isolates with the same gene (16).

In addition to detecting the rplD G70D mutation, G70S, G68D, and A43T mutations were also detected in this study. Out of the isolates with rplD mutations, only one N. gonorrhoeae strain with the G70S mutation was sensitive to azithromycin; nevertheless, the other isolates with the G70D, G70S, or G68D mutation were resistant to azithromycin. When comparing the AZM-S group to the AZM-R group, it was observed that the point mutation rate of rplD in the AZM-S group was significantly lower than in the AZM-R group (P < 0.05). This was consistent with previous research results. Although the rplD mutation did not significantly increase the MIC value, it played a role in mediating treatment failure caused by the increase of drug resistance sub-breakpoint (27). Various mutations, such as G68D, G68C, T69I, G70D, G70S, G70A, and G70R, reduced the sensitivity of macrolides by narrowing the peptide exit tunnel and interrupting the binding of macrolide antibiotics to the 50S subunit of ribosome (Figure 1) (16, 28, 29). Therefore, the detection and monitoring of these mutations are crucial in the fight against antibiotic-resistant N. gonorrhoeae.

It was generally believed that mutations in the rplD gene were associated with low-level resistance to macrolides (MIC: 0.5 mg/L); nonetheless, mutations in the rplV gene were associated with low to high-level resistance (MIC: 1.5 to 256 mg/L) (5). No mutations of the rplV gene were observed in this study, which is similar to previous research, possibly due to the relatively few N. gonorrhoeae isolates highly resistant to azithromycin. Of particular interest is the detection of the A43T mutation in the rplD gene in a high-level AZM-R N. gonorrhoeae strain, which, to the best of our knowledge, has not been reported before. The presence of this previously undescribed mutation suggests a potential mechanism for the high resistance of N. gonorrhoeae to azithromycin, although further verification is needed. Multiple previously undescribed mutations in rplD were observed to be associated with higher azithromycin MICs than the G70D mutation (16).

Laumen et al. demonstrated the G70D mutation in one strain with a MIC value of 0.5 mg/L, while the R71C mutation in 2 isolates with MIC values of 24 and 32 mg/L, respectively (5). It is suggested that mutations in rplD/rplV serve as stepping stones to higher-level resistance, with mutations first occurring in the ribosomal genes followed by the MtrCDE efflux pump and 23S rRNA (5). Although AZM-R in N. gonorrhoeae is mainly attributed to 23S rRNA point mutations and overexpression of the MtrCDE efflux pump, the A43T mutation found in this study might represent another mechanism for high resistance. Further research is needed to confirm the potential role of this mutation in antibiotic resistance (8, 30, 31). In this study, MLST typing was used to understand the spread of gonorrhea.

Among all identified MLST types, up to 33.33% (4/12) of STs were represented by only a single isolate, indicating that N. gonorrhoeae, clinically isolated in Wenzhou, China, showed considerable genetic diversity. Shimuta et al. reported that MLST ST7363 and MLST ST1901 N. gonorrhoeae showed their capacity to develop high-level in vitro resistance to ceftriaxone and called them “superbugs” (32). In this study, 5 isolates of ST7363 and 4 isolates of ST1901 were sequenced, accounting for 24.32% (9/37), indicating that the proportion of “superbug” in the Wenzhou area was quite considerable, and there was a risk of multidrug-resistant (MDR) N. gonorrhoeae, which needed to be paid great attention to.

The most common MLST type of AZM-R N. gonorrhoeae was observed to be ST1901, consistent with previous reports from Taiwan and China (16, 32). This type, along with ST7365 and ST1927, was observed to be prevalent in Taiwan, with higher MIC values for ceftriaxone and azithromycin reported for ST1901 isolates (16). The 2 isolates of highly AZM-R N. gonorrhoeae isolated in this experiment were ST1901 and ST1588, belonging to cluster B and cluster A in the goeBURST minimum spanning tree, respectively. These types differed at four loci (adk, aroE, fumC, and gdh), indicating significant genomic differences, which are similar to the previous report (10, 26).

Although ST1588 has not been previously reported as a highly AZM-R strain, it was closely related to the “superbug” ST7363, belonging to the same cluster A and differing by only 2 loci (fumC and gdh). Therefore, increased surveillance of this type is also recommended. The international emergence and spread of antibiotic-resistant N. gonorrhoeae emphasize the need to develop epidemiologically relevant surveillance systems, which not only closely track the spread of known resistant isolates but also rapidly detect novel resistance mechanisms. Furthermore, it should develop rapid and reliable genetic methods to predict antibiotic resistance. In a clinical setting, it is of paramount importance to apply appropriate therapies to limit clonal expansion, thereby having the opportunity to counteract the evolution and spread of gonorrhea.

5.1. Conclusions

Although this study has limitations, such as a relatively small overall number of N. gonorrhoeae isolates and a very small number of AZM-R isolates, especially highly resistant isolates, the mutations discovered in the rplD gene, particularly the newly identified A43T mutation, hold great significance in the research on N. gonorrhoeae resistance to azithromycin in East China. It is necessary to be vigilant about the spread of MLST ST1901, ST7363, and ST1588, especially ST1901 clones.

References

  • 1.

    Mlynarczyk-Bonikowska B, Kowalewski C, Krolak-Ulinska A, Marusza W. Molecular mechanisms of drug resistance and epidemiology of multidrug-resistant variants of Neisseria gonorrhoeae. Int J Mol Sci. 2022;23(18):10499. [PubMed ID: 36142410]. [PubMed Central ID: PMC9505821]. https://doi.org/10.3390/ijms231810499.

  • 2.

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

  • 3.

    Fifer H, Saunders J, Soni S, Sadiq ST, FitzGerald M. British Association for Sexual Health and HIV national guideline for the management of infection with Neisseria gonorrhoeae (2019). 2019. Available from: https://www.bashhguidelines.org/media/1208/gc-2019.pdf.

  • 4.

    Salmeron P, Moreno-Mingorance A, Trejo J, Amado R, Vinado B, Cornejo-Sanchez T, et al. Emergence and dissemination of three mild outbreaks of Neisseria gonorrhoeae with high-level resistance to azithromycin in Barcelona, 2016-18. J Antimicrob Chemother. 2021;76(4):930-5. [PubMed ID: 33367806]. https://doi.org/10.1093/jac/dkaa536.

  • 5.

    Laumen JGE, Manoharan-Basil SS, Verhoeven E, Abdellati S, De Baetselier I, Crucitti T, et al. Molecular pathways to high-level azithromycin resistance in Neisseria gonorrhoeae. J Antimicrob Chemother. 2021;76(7):1752-8. [PubMed ID: 33760080]. https://doi.org/10.1093/jac/dkab084.

  • 6.

    Wilson DN. Ribosome-targeting antibiotics and mechanisms of bacterial resistance. Nat Rev Microbiol. 2014;12(1):35-48. [PubMed ID: 24336183]. https://doi.org/10.1038/nrmicro3155.

  • 7.

    Unemo M, Del Rio C, Shafer WM. Antimicrobial resistance expressed by Neisseria gonorrhoeae: A major global public health problem in the 21st century. Microbiol Spectr. 2016;4(3). [PubMed ID: 27337478]. [PubMed Central ID: PMC4920088]. https://doi.org/10.1128/microbiolspec.EI10-0009-2015.

  • 8.

    Unemo M, Shafer WM. Antimicrobial resistance in Neisseria gonorrhoeae in the 21st century: Past, evolution, and future. Clin Microbiol Rev. 2014;27(3):587-613. [PubMed ID: 24982323]. [PubMed Central ID: PMC4135894]. https://doi.org/10.1128/CMR.00010-14.

  • 9.

    Halfon Y, Matzov D, Eyal Z, Bashan A, Zimmerman E, Kjeldgaard J, et al. Exit tunnel modulation as resistance mechanism of S. aureus erythromycin resistant mutant. Sci Rep. 2019;9(1):11460. [PubMed ID: 31391518]. [PubMed Central ID: PMC6685948]. https://doi.org/10.1038/s41598-019-48019-1.

  • 10.

    Zheng Z, Liu L, Shen X, Yu J, Chen L, Zhan L, et al. Antimicrobial resistance and molecular characteristics among Neisseria gonorrhoeae clinical isolates in a Chinese tertiary hospital. Infect Drug Resist. 2019;12:3301-9. [PubMed ID: 31695449]. [PubMed Central ID: PMC6815782]. https://doi.org/10.2147/IDR.S221109.

  • 11.

    Kubanov AA, Runina AV, Chestkov AV, Kudryavtseva AV, Pekov YA, Korvigo IO, et al. Whole-genome sequencing of Russian Neisseria gonorrhoeae isolates related to ST 1407 genogroup. Acta Naturae. 2018;10(3):68-76. [PubMed ID: 30397529]. [PubMed Central ID: PMC6209400].

  • 12.

    Halling SM, Jensen AE. Intrinsic and selected resistance to antibiotics binding the ribosome: analyses of Brucella 23S rrn, L4, L22, EF-Tu1, EF-Tu2, efflux and phylogenetic implications. BMC Microbiol. 2006;6:84. [PubMed ID: 17014718]. [PubMed Central ID: PMC1617103]. https://doi.org/10.1186/1471-2180-6-84.

  • 13.

    Gregory ST, Dahlberg AE. Erythromycin resistance mutations in ribosomal proteins L22 and L4 perturb the higher order structure of 23 S ribosomal RNA. J Mol Biol. 1999;289(4):827-34. [PubMed ID: 10369764]. https://doi.org/10.1006/jmbi.1999.2839.

  • 14.

    Tait-Kamradt A, Davies T, Cronan M, Jacobs MR, Appelbaum PC, Sutcliffe J. Mutations in 23S rRNA and ribosomal protein L4 account for resistance in pneumococcal strains selected in vitro by macrolide passage. Antimicrob Agents Chemother. 2000;44(8):2118-25. [PubMed ID: 10898684]. [PubMed Central ID: PMC90022]. https://doi.org/10.1128/AAC.44.8.2118-2125.2000.

  • 15.

    Manoharan-Basil SS, Laumen JGE, Van Dijck C, De Block T, De Baetselier I, Kenyon C. Evidence of horizontal gene transfer of 50S ribosomal genes rplB, rplD, and rplY in Neisseria gonorrhoeae. Front Microbiol. 2021;12:683901. [PubMed ID: 34177869]. [PubMed Central ID: PMC8222677]. https://doi.org/10.3389/fmicb.2021.683901.

  • 16.

    Ma KC, Mortimer TD, Duckett MA, Hicks AL, Wheeler NE, Sanchez-Buso L, et al. Increased power from conditional bacterial genome-wide association identifies macrolide resistance mutations in Neisseria gonorrhoeae. Nat Commun. 2020;11(1):5374. [PubMed ID: 33097713]. [PubMed Central ID: PMC7584619]. https://doi.org/10.1038/s41467-020-19250-6.

  • 17.

    Diner EJ, Hayes CS. Recombineering reveals a diverse collection of ribosomal proteins L4 and L22 that confer resistance to macrolide antibiotics. J Mol Biol. 2009;386(2):300-15. [PubMed ID: 19150357]. [PubMed Central ID: PMC2644216]. https://doi.org/10.1016/j.jmb.2008.12.064.

  • 18.

    Grad YH, Kirkcaldy RD, Trees D, Dordel J, Harris SR, Goldstein E, et al. Genomic epidemiology of Neisseria gonorrhoeae with reduced susceptibility to cefixime in the USA: A retrospective observational study. Lancet Infect Dis. 2014;14(3):220-6. [PubMed ID: 24462211]. [PubMed Central ID: PMC4030102]. https://doi.org/10.1016/S1473-3099(13)70693-5.

  • 19.

    Demczuk W, Lynch T, Martin I, Van Domselaar G, Graham M, Bharat A, 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(1):191-200. [PubMed ID: 25378573]. [PubMed Central ID: PMC4290921]. https://doi.org/10.1128/JCM.02589-14.

  • 20.

    Clinical and Laboratory Standards Institute. M100: Performance standards for antimicrobial susceptibility testing, 30th ed. 2020. Available from: https://clsi.org/media/3481/m100ed30_sample.pdf.

  • 21.

    Allen VG, Farrell DJ, Rebbapragada A, Tan J, Tijet N, Perusini SJ, et al. Molecular analysis of antimicrobial resistance mechanisms in Neisseria gonorrhoeae isolates from Ontario, Canada. Antimicrob Agents Chemother. 2011;55(2):703-12. [PubMed ID: 21098249]. [PubMed Central ID: PMC3028768]. https://doi.org/10.1128/AAC.00788-10.

  • 22.

    Demczuk W, Martin I, Peterson S, Bharat A, Van Domselaar G, Graham M, et al. Genomic epidemiology and molecular resistance mechanisms of azithromycin-resistant Neisseria gonorrhoeae in Canada from 1997 to 2014. J Clin Microbiol. 2016;54(5):1304-13. [PubMed ID: 26935729]. [PubMed Central ID: PMC4844716]. https://doi.org/10.1128/JCM.03195-15.

  • 23.

    Wi T, Lahra MM, Ndowa F, Bala M, Dillon JR, Ramon-Pardo P, et al. Antimicrobial resistance in Neisseria gonorrhoeae: Global surveillance and a call for international collaborative action. PLoS Med. 2017;14(7):e1002344. [PubMed ID: 28686231]. [PubMed Central ID: PMC5501266]. https://doi.org/10.1371/journal.pmed.1002344.

  • 24.

    Unemo M; European Sti Guidelines Editorial Board. The ‘2012 European guideline on the diagnosis and treatment of gonorrhoea in adults’ recommends dual antimicrobial therapy. Eurosurveillance. 2012;17(47):pii=20323. https://doi.org/10.2807/ese.17.47.20323-en.

  • 25.

    European Centre for Disease Prevention and Control. Extensively drug-resistant (XDR) Neisseria gonorrhoeae in the United Kingdom and Australia. 2018. Available from: https://www.ecdc.europa.eu/sites/default/files/documents/RRA-Gonorrhoea%2C%20Antimicrobial%20resistance-United%20Kingdom%2C%20Australia.pdf.

  • 26.

    Jacobsson S, Golparian D, Cole M, Spiteri G, Martin I, Bergheim T, et al. WGS analysis and molecular resistance mechanisms of azithromycin-resistant (MIC > 2 mg/L) Neisseria gonorrhoeae isolates in Europe from 2009 to 2014. J Antimicrob Chemother. 2016;71(11):3109-16. [PubMed ID: 27432597]. https://doi.org/10.1093/jac/dkw279.

  • 27.

    Colangeli R, Jedrey H, Kim S, Connell R, Ma S, Chippada Venkata UD, et al. Bacterial factors that predict relapse after tuberculosis therapy. N Engl J Med. 2018;379(9):823-33. [PubMed ID: 30157391]. [PubMed Central ID: PMC6317071]. https://doi.org/10.1056/NEJMoa1715849.

  • 28.

    Grad YH, Harris SR, Kirkcaldy RD, Green AG, Marks DS, Bentley SD, et al. Genomic epidemiology of gonococcal resistance to extended-spectrum cephalosporins, macrolides, and fluoroquinolones in the United States, 2000-2013. J Infect Dis. 2016;214(10):1579-87. [PubMed ID: 27638945]. [PubMed Central ID: PMC5091375]. https://doi.org/10.1093/infdis/jiw420.

  • 29.

    Belkacem A, Jacquier H, Goubard A, Mougari F, La Ruche G, Patey O, et al. Molecular epidemiology and mechanisms of resistance of azithromycin-resistant Neisseria gonorrhoeae isolated in France during 2013-14. J Antimicrob Chemother. 2016;71(9):2471-8. [PubMed ID: 27301565]. https://doi.org/10.1093/jac/dkw182.

  • 30.

    Chisholm SA, Dave J, Ison CA. High-level azithromycin resistance occurs in Neisseria gonorrhoeae as a result of a single point mutation in the 23S rRNA genes. Antimicrob Agents Chemother. 2010;54(9):3812-6. [PubMed ID: 20585125]. [PubMed Central ID: PMC2935028]. https://doi.org/10.1128/AAC.00309-10.

  • 31.

    Cousin Jr SL, Whittington WL, Roberts MC. Acquired macrolide resistance genes and the 1 bp deletion in the mtrR promoter in Neisseria gonorrhoeae. J Antimicrob Chemother. 2003;51(1):131-3. [PubMed ID: 12493797]. https://doi.org/10.1093/jac/dkg040.

  • 32.

    Shimuta K, Watanabe Y, Nakayama S, Morita-Ishihara T, Kuroki T, Unemo M, et al. Emergence and evolution of internationally disseminated cephalosporin-resistant Neisseria gonorrhoeae clones from 1995 to 2005 in Japan. BMC Infect Dis. 2015;15:378. [PubMed ID: 26381611]. [PubMed Central ID: PMC4574456]. https://doi.org/10.1186/s12879-015-1110-x.