A Case-Control Study of Molecular Epidemiology in Relation to Azithromycin Resistance in Neisseria gonorrhoeae Isolates Collected in Amsterdam, the Netherlands, between 2008 and 2015

Molecular epidemiology of penicillinase-producing <i>Neisseria gonorrhoeae</i> isolates and their <i>bla</i>_TEM-135 gene variant in Bangkok, Thailand, 2015–2017

Penicillinase-producing Neisseria gonorrhoeae (PPNG) possessing bla_TEM-135 is a serious public health threat. With only a single change in the amino acid sequence, bla_TEM-135 could evolve into a TEM-type extended-spectrum beta-lactamase (ESBL), which hydrolyzes extended-spectrum cephalosporins, including ceftriaxone and cefixime. We investigated the molecular epidemiological characteristics, types of plasmids in PPNG isolates, and prevalence of PPNG clinical isolates producing TEM-135 beta-lactamases. N. gonorrhoeae multi-antigen sequence typing (NG-MAST) was used to determine the molecular epidemiological characteristics of 99 PPNG isolates collected from 2015 to 2017. A mismatch amplification mutation assay was used to examine the bla_TEM-135 gene prevalence. Of the 89 identified NG-MAST sequence types, 65 (73.0%) were novel. Only 17.7% (43/243) of PPNG isolates belonged to 16 genogroups. The most frequent plasmid was African, followed by Rio/Toronto, and Asian. The bla_TEM-135 allele was found in Rio/Toronto plasmids. The bla_TEM-135 allele was present in 23.2% (23/99) of the PPNG isolates. PPNG isolates expressing TEM-135 beta-lactamase exhibited significantly higher penicillin MIC (minimum inhibitory concentration) values than TEM-1 PPNG isolates. The PPNG isolates showed high genetic diversity and a high proportion of bla_TEM-135 alleles. Mutation of the bla_TEM-135 allele is worrisome as only one mutation could cause TEM-1 to evolve into an ESBL variant that degrades ceftriaxone. Ongoing surveillance of bla_TEM-135 and new PPNG isolates is imperative.

Trends and regional variations of gonococcal antimicrobial resistance in the Netherlands, 2013 to 2019

Background

Gonococcal antimicrobial resistance is emerging worldwide and is monitored in the Netherlands in 18 of 24 Sexual Health Centres (SHC).

Aim

To report trends, predictors and regional variation of gonococcal azithromycin resistance (AZI-R, minimum inhibitory concentration (MIC) > 1 mg/L) and ceftriaxone decreased susceptibility (CEF-DS, MIC > 0.032 mg/L) in 2013–2019.

Methods

SHC reported data on individual characteristics, sexually transmitted infection diagnoses, and susceptibility testing (MIC, measured by Etest). We used multilevel logistic regression analysis to identify AZI-R/CEF-DS predictors, correcting for SHC region. Population differences’ effect on regional variance of AZI-R and CEF-DS was assessed with a separate multilevel model.

Results

The study included 13,172 isolates, predominantly (n = 9,751; 74%) from men who have sex with men (MSM). Between 2013 and 2019, annual proportions of AZI-R isolates appeared to increase from 2.8% (37/1,304) to 9.3% (210/2,264), while those of CEF-DS seemed to decrease from 7.0% (91/1,306) to 2.9% (65/2,276). Among SHC regions, 0.0‒16.9% isolates were AZI-R and 0.0−7.0% CEF-DS; population characteristics could not explain regional variance. Pharyngeal strain origin and consultation year were significantly associated with AZI-R and CEF-DS for MSM, women, and heterosexual men. Among women and heterosexual men ≥ 4 partners was associated with CEF-DS, and ≥ 10 with AZI-R.

Conclusions

No resistance or decreasing susceptibility was found for CEF, the first line gonorrhoea treatment in the Netherlands. Similar to trends worldwide, AZI-R appeared to increase. Regional differences between SHC support nationwide surveillance with regional-level reporting. The increased risk of resistance/decreased susceptibility in pharyngeal strains underlines the importance of including extragenital infections in gonococcal resistance surveillance.

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.

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.

Emergence of a Neisseria gonorrhoeae clone with reduced cephalosporin susceptibility between 2014 and 2019 in Amsterdam, The Netherlands, revealed by genomic population analysis

BACKGROUND

Emerging resistance to cephalosporins in Neisseria gonorrhoeae (Ng) is a major public health threat, since these are considered antibiotics of last resort. Continuous surveillance is needed to monitor the circulation of resistant strains and those with reduced susceptibility.

OBJECTIVES

For the purpose of epidemiological surveillance, genomic population analysis was performed on Ng isolates from Amsterdam with a focus on isolates with reduced susceptibility to ceftriaxone.

METHODS

WGS data were obtained from 318 isolates from Amsterdam, the Netherlands between 2014 and 2019. Isolates were typed according to MLST, Ng Multi-Antigen Sequence Typing (NG-MAST) and Ng Sequence Typing for Antimicrobial Resistance (NG-STAR) schemes and additional resistance markers were identified. Phylogenetic trees were created to identify genetic clusters and to compare Dutch and non-Dutch MLST7827 isolates.

RESULTS

MLST7363 and MLST1901 were the predominant strains having reduced susceptibility to ceftriaxone during 2014-16; MLST7827 emerged and dominated during 2017-19. NG-STAR38 and NG-MAST2318/10386 were predominant among MLST7827 isolates. MLST7827 reduced susceptibility isolates carried a non-mosaic 13.001 penA allele with an A501V mutation and porB1b G120K/A121D mutations, which were lacking in susceptible MLST7827 isolates. Phylogenetic analysis of all publicly available MLST7827 isolates showed strong genetic clustering of Dutch and other European MLST7827 isolates.

CONCLUSIONS

MLST7827 isolates with reduced ceftriaxone susceptibility have emerged during recent years in Amsterdam. Co-occurrence of penA A501V and porB1b G120K/A121D mutations was strongly associated with reduced susceptibility to ceftriaxone. Genetic clustering of Dutch and other European MLST7827 isolates indicates extensive circulation of this strain in Europe. Close monitoring of the spread of this strain having an alarming susceptibility profile is needed.

Association of Neisseria gonorrhoeae genogroups and specific PBP2/MtrR/PorB mutation patterns with susceptibility to penicillin in a susceptible gonococcal population

Objectives

To ascertain whether the antimicrobial susceptibility of Neisseria gonorrhoeae isolates with differing susceptibilities to penicillin is associated with genogroups (GGs) and combined mutation patterns in PBP2 (penA), the multiple transfer resistance repressor (MtrR; mtrR) and porin B (PorB; porB).

Methods

The susceptibility of 146 clinical N. gonorrhoeae isolates to penicillin was determined using the agar dilution method and the interpretation criteria of CLSI. The DNA sequences of penA, mtrR and porB in isolates were compared with WT sequences and mutation patterns were determined. Isolates were typed by N. gonorrhoeae multi-antigen sequence typing (NG-MAST) and STs were grouped into specific GGs.

Results

The isolates tested carried 9 mutation patterns in PBP2 and 12 mutation patterns in each of MtrR and PorB. Of the 146 isolates, 121 (82.9%) were grouped into 13 different GGs. Isolates with penicillin MICs of 0.03-0.06 mg/L were significantly associated with GG25 (P < 0.05) and PBP2/MtrR/PorB mutation pattern I/WT/WT (P < 0.01). Isolates with a penicillin MIC of 1.0 mg/L were associated (P < 0.05) with: (i) GG3655 and mutation pattern XXII/A-;G45D/G120K;A121N; (ii) GG921 and mutation pattern IX/G45D/G120D;A121N; and (iii) GG1109 and mutation pattern IX/G45D/WT. Sixty percent (9/15) of penicillin-resistant isolates (MIC ≥2 mg/L) were GG3654 (P < 0.0001) and carried mutation pattern IX/G45D/G120K;A121D or IX/G45D/G120D;A121D (P < 0.05).

Conclusions

Specific mutation patterns in PBP2/MtrR/PorB were associated with specific GGs and penicillin susceptibility. This approach of typing strains and resistance patterns is ideal for predicting antimicrobial resistance and should be used in instances in which gonococcal culture is not available but DNA can be obtained from clinical specimens.

Emergence and Spread of Neisseria gonorrhoeae Strains with High-Level Resistance to Azithromycin in Taiwan from 2001 to 2018

A total of 598 Neisseria gonorrhoeae isolates obtained from patients in Taiwan from 2001 to 2018 were evaluated. The MICs of ceftriaxone (CRO) and azithromycin (AZM) against the isolates were determined by the agar dilution method. N. gonorrhoeae isolates with AZM MICs of ≥1 μg/ml were identified and characterized by the presence of AZM resistance determinants. For high-level AZM-resistant (AZM-HLR) isolates (MIC ≥ 256 μg/ml), genotyping was performed using multilocus sequence typing (MLST) and N. gonorrhoeae multiantigen sequence typing (NG-MAST). Among the N. gonorrhoeae isolates studied, 8.7% (52/598) exhibited AZM MICs of ≥1 μg/ml. Thirteen of the 52 isolates contained A2059G (23S rRNA NG-STAR type 1) or C2611T (23S rRNA NG-STAR type 2) mutations. The prevalence of the A2059G mutation was higher in AZM-HLR isolates (P < 0.001). The -35A deletion in the promoter region of the mtrR gene did not differ between AZM-HLR isolates (100%, 10/10) and the isolates with AZM MICs of 1 μg/ml to 64 μg/ml (95.2%, 40/42) (P = 1.000). The presence of mutations in the mtrR coding region was significantly different between these two groups at 90% (9/10) and 26.2% (11/42), respectively (P < 0.001). The AZM-HLR isolates, all carrying four mutated A2059G alleles, a -35A deletion, and G45D, were classified as MLST 12039/10899 and NG-MAST 1866/16497. In conclusion, Taiwan is among the countries reporting gonococci with high-level resistance to AZM so that a single dose of 1 g ceftriaxone intramuscularly as the first choice for management of N. gonorrhoeae infection should be evaluated.

Epidemiological surveillance study of gonococcal infection in Northern Spain

BACKGROUND

Treatment of gonorrhoea is threatened by antimicrobial resistance, and decreased susceptibility to recommended therapies is emerging. Thus, gonococcal infection (GI) is becoming a public health problem. The objectives of the present study were to monitor the antimicrobial sensitivity in Neisseria gonorrhoeae (NG) during 2011-2015 and to study their genogroups.

METHODS

Antimicrobial susceptibility was studied by disc diffusion, in addition to the agar dilution method for cefixime and ceftriaxone and the Etest® for azithromycin. Genotyping was performed by the NG multi-antigen sequence typing (NG-MAST) method. Genogroups of closely related sequence types (STs) were defined.

RESULTS

All the strains were susceptible to cefixime, ceftriaxone and gentamicin and 1.8% of the strains were resistant to azithromycin. A total of 531 STs and 6 genotypes (Gs) were identified during 2012-2015 period. G2992 was the largest and was associated with resistance to azithromycin, and with men who have sex with men (MSM), alongside G2400. G1407 and G2400 strains were related to high minimum inhibitory concentration (MICs) to cefixime and G1407 also to ceftriaxone. For the first time, G1861 and G2018 were described and associated with ciprofloxacin resistance and G2018 also with high MICs to ceftriaxone.

CONCLUSION

Molecular typing is a useful tool to predict antimicrobial resistance. These results show the need to develop novel antimicrobials or to design new antimicrobial therapies based on drugs that show their efficacy against GI. This also highlights the importance of developing sexually transmitted infection (STI) surveillance in homosexual populations.

Azithromycin-resistant Neisseria gonorrhoeae spreading amongst men who have sex with men (MSM) and heterosexuals in New South Wales, Australia, 2017

Objectives

To identify the genetic basis of resistance as well as to better understand the epidemiology of a recent surge in azithromycin-resistant Neisseria gonorrhoeae in New South Wales, Australia.

Methods

Azithromycin-resistant N. gonorrhoeae isolates (n = 118) collected from 107 males, 10 females and 1 transsexual between January and July 2017 were genotyped using a previously described iPLEX method. The results were compared with phenotypic resistance profiles and available patient data.

Results

The iPLEX results revealed 10 different N. gonorrhoeae genotypes (designated AZI-G1 to AZI-G10) of which three were responsible for the majority of infections; AZI-G10 (74.6%, 88 isolates; 87 males and 1 transsexual), AZI-G4 (11.0%, 13 isolates; 7 males and 6 females) and AZI-G7 (6.8%, 8 isolates; 7 males and 1 female). The observed resistance was attributable to one of two different azithromycin resistance mechanisms; the 23S rRNA C2611T mutation was identified in 24% of isolates, whereas the majority of resistance (76%) was associated with a meningococcal-type mtrR variant. Additionally, one isolate was found to harbour both the 23S rRNA C2611T mutation and a type XXXIV mosaic penA sequence associated with cephalosporin resistance.

Conclusions

These data indicate outbreaks of azithromycin-resistant gonococci amongst networks of MSM and heterosexuals in New South Wales. The results also provide further evidence that azithromycin may soon be an ineffective treatment option for gonococcal infection and highlight the urgent need to explore alternative therapies.

Emergence of Neisseria gonorrhoeae Strains Harboring a Novel Combination of Azithromycin-Attenuating Mutations

The nimbleness of Neisseria gonorrhoeae to evade the effect of antibiotics has perpetuated the fight against antibiotic-resistant gonorrhea for more than 80 years. The ability to develop resistance to antibiotics is attributable to its indiscriminate nature in accepting and integrating exogenous DNA into its genome. Here, we provide data demonstrating a novel combination of the 23S rRNA A2059G mutation with a mosaic-multiple transferable resistance (mosaic-mtr) locus haplotype in 14 N. gonorrhoeae isolates with high-level azithromycin MICs (≥256 μg/ml), a combination that may confer more fitness than in previously identified isolates with high-level azithromycin resistance. To our knowledge, this is the first description of N. gonorrhoeae strains harboring this novel combination of resistance determinants. These strains were isolated at two independent jurisdictions participating in the Gonococcal Isolate Surveillance Project (GISP) and in the Strengthening the U.S. Response to Resistant Gonorrhea (SURRG) project. The data suggest that the genome of N. gonorrhoeae continues to shuffle its genetic material. These findings further illuminate the genomic plasticity of N. gonorrhoeae, which allows this pathogen to develop mutations to escape the inhibitory effects of antibiotics.

Culture-free genotyping of Neisseria gonorrhoeae revealed distinct strains at different anatomical sites in a quarter of patients, the Netherlands, 2012 to 2016

BackgroundGenotyping of Neisseria gonorrhoeae (NG) is essential for surveillance to monitor NG transmission and dissemination of resistant strains. Current genotyping methods rely on bacterial culture which frequently fails.AimOur aim was to develop a culture-free genotyping method that is compatible with the widely used N. gonorrhoeae multi-antigen sequence typing (NG-MAST) database, which facilitates genotyping of NG detected at separate anatomical sites in individual patients.MethodsSpecific primers for both PCR targets porB and tbpB were designed and technically validated by assessing the analytical sensitivity, cross-reactivity with 32 non-gonoccocal Neisseria species, and concordance with NG-MAST. Clinical application was assessed on 205 paired samples from concurrent NG infections at different anatomical sites of 98 patients (81 men who have sex with men and 17 women) visiting our sexually transmitted infections clinic.ResultsTyping could be consistently performed on samples with a PCR quantification cycle (Cq) value <35. Furthermore, the method showed no cross-reactivity and was concordant with NG-MAST. Culture-free NG-MAST improved the typing rate from 62% (59/95) for cultured samples to 94% (89/95) compared with culture-dependent NG-MAST. Paired samples of 80 of 98 patients were genotyped, revealing distinct NG strains in separate anatomical sites in 25% (20/80) of the patients.ConclusionsThis NG-specific genotyping method can improve NG surveillance as it facilitates genotyping of non-culturable and extra-genital samples. Furthermore, 25% of patients were infected with multiple NG strains, which is missed in current culture-dependent surveillance. Including non-culturable and concurrent NG infections in surveillance informs actions on dissemination of multidrug-resistant NG strains.

Molecular epidemiology of Neisseria gonorrhoeae strains circulating in Indonesia using multi-locus variable number tandem repeat analysis (MLVA) and Neisseria gonorrhoeae multi-antigen sequence typing (NG-MAST) techniques

BACKGROUND

Control of gonorrhea in resource-limited countries, such as Indonesia, is mostly unsuccessful. Examining Neisseria gonorrhoeae (Ng) transmission networks using strain typing might help prioritizing public health interventions.

METHODS

In 2014, urogenital Ng strains were isolated from clients of sexually transmitted infection clinics in three Indonesian cities. Strains were typed using Multiple-Locus Variable Number Tandem Repeat (VNTR) Analysis (MLVA) and Ng Multi-Antigen Sequence Typing (NG-MAST) at the Public Health Service, Amsterdam, the Netherlands, and compared to Dutch strains collected from 2012 to 2015. Minimum spanning trees (MSTs) were constructed using MLVA profiles incorporating demographics and NG-MAST genogroups. A cluster was defined as ≥5 strains differing in ≤1 VNTR locus. The concordance between MLVA and NG-MAST was examined with the adjusted Wallace coefficients (AW).

RESULTS

We collected a total of 78 Indonesian strains from Yogyakarta (n = 44), Jakarta (n = 25), and Denpasar (n = 9). Seven MLVA clusters and 16 non-clustered strains were identified. No cluster was specific for any geographic area, risk group or age group. Combined with 119 contemporary Dutch strains, 8 MLVA clusters were identified, being four clusters of Indonesian strains, two of Dutch strains, and two of both Indonesian and Dutch strains. Indonesian strains (79.5%) were more often clustered compared to Dutch strains (24.3%). The most prevalent NG-MAST genogroups among Indonesian strains was G1407 (51.3%) and among Dutch strains was G2992 (19.3%). In Indonesian strains, the AW [95% confidence interval] for MLVA to NG-MAST was 0.07 [0.00-0.27] and for NG-MAST to MLVA was 0.03 [0.00-0.12].

CONCLUSION

Indonesian Ng strains are more often clustered than Dutch strains, but show no relation with geographical area, risk group, or age group, suggesting a more clonal Ng epidemic in Indonesia. Some similar Ng strains circulate in both Indonesia and the Netherlands.