Association of Neisseria gonorrhoeae Plasmids With Distinct Lineages and The Economic Status of Their Country of Origin

Integrated molecular, phenotypic and epidemiological surveillance of antimicrobial resistance in Neisseria gonorrhoeae in Germany

Numbers of infections with Neisseria gonorrhoeae are among the top three sexually transmitted infections (STI) worldwide. In addition, the emergence and spread of antimicrobial resistance (AMR) in Neisseria gonorrhoeae pose an important public-health issue. The integration of genomic, phenotypic and epidemiological data to monitor Neisseria gonorrhoeae fosters our understanding of the emergence and spread of AMR in Neisseria gonorrhoeae and helps to inform therapy guidelines and intervention strategies. Thus, the Gonococcal resistance surveillance (Go-Surv-AMR) was implemented at the Robert Koch Institute in Germany in 2021 to obtain molecular, phenotypic and epidemiological data on Neisseria gonorrhoeae isolated in Germany. Here, we describe the structure and aims of Go-Surv-AMR. Furthermore, we point out future directions of Go-Surv-AMR to improve the integrated genomic surveillance of Neisseria gonorrhoeae. In this context we discuss current and prospective sequencing approaches and the information derived from their application. Moreover, we highlight the importance of combining phenotypic and WGS data to monitor the evolution of AMR in Neisseria gonorrhoeae in Germany. The implementation and constant development of techniques and tools to improve the genomic surveillance of Neisseria gonorrhoeae will be important in coming years.

Antimicrobial susceptibility profile of Neisseria gonorrhoeae from patients attending a medical laboratory, Institut Pasteur de Madagascar between 2014 and 2020: phenotypical and genomic characterisation in a subset of Neisseria gonorrhoeae isolates

OBJECTIVES

Antimicrobial-resistant Neisseria gonorrhoeae (NG) is a concern. Little is known about antimicrobial susceptibility profiles and associated genetic resistance mechanisms of NG in Madagascar. We report susceptibility data of NG isolates obtained by the medical laboratory (CBC) of the Institut Pasteur de Madagascar, Antananarivo, Madagascar, during 2014-2020. We present antimicrobial resistance mechanisms data and phenotype profiles of a subset of isolates.

METHODS

We retrieved retrospective data (N=395) from patients with NG isolated during 2014-2020 by the CBC. We retested 46 viable isolates including 6 found ceftriaxone and 2 azithromycin resistant, as well as 33 isolated from 2020. We determined minimal inhibitory concentrations for ceftriaxone, ciprofloxacin, azithromycin, penicillin, tetracycline and spectinomycin using Etest. We obtained whole-genome sequences and identified the gene determinants associated with antimicrobial resistance and the sequence types (STs).

RESULTS

Over the study period, ceftriaxone-resistant isolates exceeded the threshold of 5% in 2017 (7.4% (4 of 54)) and 2020 (7.1% (3 of 42)). All retested isolates were found susceptible to ceftriaxone, azithromycin and spectinomycin, and resistant to ciprofloxacin. The majority were resistant to penicillin (83% (38 of 46)) and tetracycline (87% (40 of 46)). We detected chromosomal mutations associated with antibiotic resistance in gyrA, parC, penA, ponA, porB and mtrR genes. None of the retested isolates carried the mosaic penA gene. The high rate of resistance to penicillin and tetracycline is explained by the presence of bla TEM (94.7% (36 of 38)) and tetM (97.5% (39 of 40)). We found a high number of circulating multilocus STs. Almost half of them were new types, and one new type was among the four most predominant.

CONCLUSIONS

Our report provides a detailed dataset obtained through phenotypical and genotypical methods which will serve as a baseline for future surveillance of NG. We could not confirm the occurrence of ceftriaxone-resistant isolates. Our results highlight the importance of implementing quality-assured gonococcal antimicrobial resistance surveillance in Madagascar.

Evolution and exchange of plasmids in pathogenic Neisseria

IMPORTANCE

Horizontal gene transfer (HGT) is a major influence in driving the spread of antimicrobial resistance (AMR) in many bacteria. A conjugative plasmid which is widespread in Neisseria gonorrhoeae, pConj, prevented the use of tetracycline/doxycycline for treating gonococcal infection. Here, we show that pConj evolved in the related pathogen, Neisseria meningitidis, and has been acquired by the gonococcus from the meningococcus on multiple occasions. Following its initial acquisition, pConj spread to different gonococcal lineages; changes in the plasmid's conjugation machinery associated with another transfer event limit spread in the gonococcal populations. Our findings have important implications for the use of doxycycline to prevent bacterial sexually transmitted disease which is likely to exacerbate the spread of AMR through HGT in pathogenic bacteria.

Phenotypic and genotypic characterization of Neisseria gonorrhoeae isolates from Yaoundé, Cameroon, 2019 to 2020

This study investigated antimicrobial resistance (AMR) phenotypes and genotypes exhibited by Neisseria gonorrhoeae from Yaoundé, Cameroon. AMR to tetracycline, penicillin and ciprofloxacin was observed although none of the isolates had reduced susceptibility to azithromycin, cefixime or ceftriaxone. Whole genome sequence (WGS) data were obtained and, using a threshold of 300 or fewer locus differences in the N. gonorrhoeae core gene multilocus sequence typing (cgMLST) scheme, four distinct core genome lineages were identified. Publicly available WGS data from 1355 gonococci belonging to these four lineages were retrieved from the PubMLST database, allowing the Cameroonian isolates to be examined in the context of existing data and compared with related gonococci. Examination of AMR genotypes in this dataset found an association between the core genome and AMR with, for example, isolates belonging to the core genome group, Ng_cgc_300 : 21, possessing GyrA and ParC alleles with amino acid substitutions conferring high-level resistance to ciprofloxacin while lineages Ng_cgc_300 : 41 and Ng_cgc_300 : 243 were predicted to be susceptible to several antimicrobials. A core genome lineage, Ng_cgc_300 : 498, was observed which largely consisted of gonococci originating from Africa. Analyses from this study demonstrate the advantages of using the N. gonorrhoeae cgMLST scheme to find related gonococci to carry out genomic analyses that enhance our understanding of the population biology of this important pathogen.

There are three major Neisseria gonorrhoeae β-lactamase plasmid variants which are associated with specific lineages and carry distinct TEM alleles

Neisseria gonorrhoeae is a significant threat to global health with an estimated incidence of over 80 million cases each year and high levels of antimicrobial resistance. The gonococcal β-lactamase plasmid, pbla, carries the TEM β-lactamase, which requires only one or two amino acid changes to become an extended-spectrum β-lactamase (ESBL); this would render last resort treatments for gonorrhoea ineffective. Although pbla is not mobile, it can be transferred by the conjugative plasmid, pConj, found in N. gonorrhoeae. Seven variants of pbla have been described previously, but little is known about their frequency or distribution in the gonococcal population. We characterised sequences of pbla variants and devised a typing scheme, Ng_pblaST that allows their identification from whole genome short-read sequences. We implemented Ng_pblaST to assess the distribution of pbla variants in 15 532 gonococcal isolates. This demonstrated that only three pbla variants commonly circulate in gonococci, which together account for >99 % of sequences. The pbla variants carry different TEM alleles and are prevalent in distinct gonococcal lineages. Analysis of 2758 pbla-containing isolates revealed the co-occurrence of pbla with certain pConj types, indicating co-operativity between pbla and pConj variants in the spread of plasmid-mediated AMR in N. gonorrhoeae. Understanding the variation and distribution of pbla is essential for monitoring and predicting the spread of plasmid-mediated β-lactam resistance in N. gonorrhoeae.

Multisite transformation in Neisseria gonorrhoeae: insights on transformations mechanisms and new genetic modification protocols

Natural transformation, or the uptake of naked DNA from the external milieu by bacteria, holds a unique place in the history of biology. This is both the beginning of the realization of the correct chemical nature of genes and the first technical step to the molecular biology revolution that sees us today able to modify genomes almost at will. Yet the mechanistic understanding of bacterial transformation still presents many blind spots and many bacterial systems lag behind power horse model systems like Escherichia coli in terms of ease of genetic modification. Using Neisseria gonorrhoeae as a model system and using transformation with multiple DNA molecules, we tackle in this paper both some aspects of the mechanistic nature of bacterial transformation and the presentation of new molecular biology techniques for this organism. We show that similarly to what has been demonstrated in other naturally competent bacteria, Neisseria gonorrhoeae can incorporate, at the same time, different DNA molecules modifying DNA at different loci within its genome. In particular, co-transformation of a DNA molecule bearing an antibiotic selection cassette and another non-selected DNA piece can lead to the integration of both molecules in the genome while selecting only through the selective cassette at percentages above 70%. We also show that successive selections with two selection markers at the same genetic locus can drastically reduce the number of genetic markers needed to do multisite genetic modifications in Neisseria gonorrhoeae. Despite public health interest heightened with the recent rise in antibiotic resistance, the causative agent of gonorrhea still does not possess a plethora of molecular techniques. This paper will extend the techniques available to the Neisseria community while providing some insights into the mechanisms behind bacterial transformation in Neisseria gonorrhoeae. We are providing a suite of new techniques to quickly obtain modifications of genes and genomes in the Neisserial naturally competent bacteria.

Evolution, persistence, and host adaption of a gonococcal AMR plasmid that emerged in the pre-antibiotic era

Plasmids are diverse extrachromosomal elements significantly contributing to interspecies dissemination of antimicrobial resistance (AMR) genes. However, within clinically important bacteria, plasmids can exhibit unexpected narrow host ranges, a phenomenon that has scarcely been examined. Here we show that pConj is largely restricted to the human-specific pathogen, Neisseria gonorrhoeae. pConj can confer tetracycline resistance and is central to the dissemination of other AMR plasmids. We tracked pConj evolution from the pre-antibiotic era 80 years ago to the modern day and demonstrate that, aside from limited gene acquisition and loss events, pConj is remarkably conserved. Notably, pConj has remained prevalent in gonococcal populations despite cessation of tetracycline use, thereby demonstrating pConj adaptation to its host. Equally, pConj imposes no measurable fitness costs and is stably inherited by the gonococcus. Its maintenance depends on the co-operative activity of plasmid-encoded Toxin:Antitoxin (TA) and partitioning systems rather than host factors. An orphan VapD toxin encoded on pConj forms a split TA with antitoxins expressed from an ancestral co-resident plasmid or a horizontally-acquired chromosomal island, potentially explaining pConj's limited distribution. Finally, ciprofloxacin can induce loss of this highly stable plasmid, reflecting epidemiological evidence of transient local falls in pConj prevalence when fluoroquinolones were introduced to treat gonorrhoea.

Whole-Genome Sequencing of Neisseria gonorrhoeae Isolates with Antimicrobial Resistance Plasmids from a Homeless Community in Colombia

Neisseria gonorrhoeae is a pathogenic bacterium causing sexually transmitted infections, and it is associated with high antibiotic resistance rates. Here, we describe the genome sequences of four isolates from a homeless community in Colombia.

Whole-genome sequence analysis of high-level penicillin-resistant strains and antimicrobial susceptibility of Neisseria gonorrhoeae clinical isolates from Thailand

Background

The increasing rate of antimicrobial-resistant Neisseria gonorrhoeae poses a considerable public health threat due to the difficulty in treating gonococcal infections. This study examined antimicrobial resistance (AMR) to drugs recommended for gonorrhea treatment between 2015 and 2017, and the AMR determinants and genetic compositions of plasmids in 3 gonococcal strains with high-level penicillin resistance.

Methods

We collected 117 N. gonorrhoeae isolates from patients with gonococcal infections who attended Siriraj Hospital, Bangkok, Thailand, between 2015 and 2017. Minimum inhibitory concentrations (MICs) of penicillin, tetracycline, ciprofloxacin, azithromycin, spectinomycin, cefixime, and ceftriaxone were determined by the agar dilution method. PCR amplification and sequencing of 23S rRNA and mtrR (a negative regulator of MtrCDE efflux pump) were performed. Whole genomes of 3 PPNG strains with high-level penicillin resistance (MIC ≥ 128 μg/ml) were sequenced using Illumina and Nanopore sequencing platforms.

Results

The proportions of N. gonorrhoeae isolates with resistance were 84.6% for penicillin, 91.5% for tetracycline, and 96.6% for ciprofloxacin. All isolates were susceptible to spectinomycin, azithromycin, cefixime, and ceftriaxone. An adenine deletion within a 13 bp inverted repeat sequence in the mtrR promoter and an H105Y mutation in the mtrR coding region were found in the N. gonorrhoeae isolate with the highest azithromycin MIC value (1 μg/ml). Three high-level penicillin-resistant isolates contained nonmosaic type II penA and had mutations in penB and the mtrR coding region. All isolates with high-level penicillin resistance carried the conjugative plasmids with or without the Dutch type tetM determinant, the beta-lactamase plasmid (Rio/Toronto), and the cryptic plasmid.

Conclusions

The gonococcal population in Thailand showed high susceptibility to ceftriaxone and azithromycin, current dual therapy recommended for gonorrhea treatment. As elevated MIC of azithromycin has been observed in 1 strain of N. gonorrhoeae, expanded and enhanced surveillance of antimicrobial susceptibility and study of genetic resistance determinants are essential to improve treatment guidelines.

The Neisseria gonorrhoeae Accessory Genome and Its Association with the Core Genome and Antimicrobial Resistance

The bacterial accessory genome provides the genetic flexibility needed to facilitate environment and host adaptation. In Neisseria gonorrhoeae, known accessory elements include plasmids which can transfer and mediate antimicrobial resistance (AMR); however, chromosomal accessory genes could also play a role in AMR. Here, the gonococcal accessory genome was characterized using gene-by-gene approaches and its association with the core genome and AMR were assessed. The gonococcal accessory gene pool consisted of 247 genes, which were mainly genes located on large mobile genetic elements, phage associated genes, or genes encoding putative secretion systems. Accessory elements showed similar synteny across genomes, indicating either a predisposition for particular genomic locations or ancestral inheritance that are conserved during strain expansion. Significant associations were found between the prevalence of accessory elements and core genome multi-locus sequence types (cgMLST), consistent with a structured gonococcal population despite frequent horizontal gene transfer (HGT). Increased prevalence of putative DNA exchange regulators was significantly associated with AMR, which included a putative secretion system, methyltransferases and a toxin-antitoxin system. Although frequent HGT results in high genetic diversity in the gonococcus, we found that this is mediated by a small gene pool. In fact, a highly organized genome composition was identified with a strong association between the accessory and core genome. Increased prevalence of DNA exchange regulators in antimicrobial resistant isolates suggests that genetic material exchange plays a role in the development or maintenance of AMR. These findings enhance our understanding of gonococcal genome architecture and have important implications for gonococcal population biology. IMPORTANCE The emergence of antimicrobial resistance (AMR) against third generation cephalosporins in Neisseria gonorrhoeae is a major public health concern, as these are antibiotics of last resort for the effective treatment of gonorrhea. Although the resistance mechanisms against this class of antibiotics have not been entirely resolved, resistance against other classes of antibiotics, such as tetracyclines, is known to be mediated through plasmids, which are known gonococcal extra-chromosomal accessory elements. A complete assessment of the chromosomal accessory genome content and its role in AMR has not yet been undertaken. Here, we comprehensively characterize the gonococcal accessory genome to better understand genome architecture as well as the evolution and mechanisms of AMR in this species.

High prevalence of blaTEM-135 and genetic epidemiology of blaTEM-135-carrying Neisseria gonorrhoeae isolates in Shandong, China, 2017-19

BACKGROUND

Only one additional specific SNP in blaTEM-135 in penicillinase-producing Neisseria gonorrhoeae (PPNG) isolates is required for TEM-135 to evolve into an ESBL that is capable of degrading all extended-spectrum cephalosporins (ESCs).

METHODS

Identification of the blaTEM-135 gene was achieved by performing a mismatch amplification mutation assay PCR. PPNG isolates were then genotyped using MLST and N. gonorrhoeae multiantigen sequence typing (NG-MAST) techniques. Finally, phylogenetic analyses based on MLST and NG-MAST were performed on all of the PPNG isolates, and genogroups were classified according to the phylogenetic relatedness on the basis of each MLST and NG-MAST gene typing and the blaTEM allele varieties.

RESULTS

Among the 197 PPNG isolates, one (0.5%) isolate possessed the blaTEM-220 allele, 90 (45.7%) possessed the blaTEM-135 allele, and 106 (53.8%) possessed the blaTEM-1 allele. Among the blaTEM-135-carrying PPNG isolates, MLST ST8109 (n = 30) and NG-MAST ST12199 (n = 5), ST12754 (n = 5) and ST15073 (n = 5) were the most prevalent in the Shandong province of China. The phylogenetic analyses revealed that MLST ST8109, ST8140 and ST14417, mainly carrying the blaTEM-135 allele, belonged to the same genogroup, G1, while NG-MAST ST12199 and ST15073, mainly carrying the blaTEM-135 allele, belonged to the same genogroup, Ga.

CONCLUSIONS

This study has shown that the molecular epidemiology of PPNG isolates carrying the blaTEM-135 allele in Shandong has a dynamically changing status. Therefore, it is very pressing to continuously monitor the prevalence and mutation of the blaTEM-135 allele and the genetic epidemiology of PPNG isolates carrying the blaTEM-135 allele in this district.

Markerless gene editing in Neisseria gonorrhoeae

Neisseria gonorrhoeae, the gonococcus, is a pathogen of major public health concern, but sophisticated approaches to gene manipulation are limited for this species. For example, there are few methods for generating markerless mutations, which allow the generation of precise point mutations and deletions without introducing additional DNA sequence. Markerless mutations are central to studying pathogenesis, the spread of antimicrobial resistance (AMR) and for vaccine development. Here we describe the use of galK as a counter-selectable marker that can be used for markerless mutagenesis in N. gonorrhoeae. galK encodes galactokinase, an enzyme that metabolizes galactose in bacteria that can utilize it as a sole carbon source. GalK can also phosphorylate a galactose analogue, 2-deoxy-galactose (2-DOG), into a toxic, non-metabolisable intermediate, 2-deoxy-galactose-1-phosphate. We utilized this property of GalK to develop a markerless approach for mutagenesis in N. gonorrhoeae. We successfully deleted both chromosomally and plasmid-encoded genes, that are important for gonococcal vaccine development and studies of AMR spread. We designed a positive-negative selection cassette, based on an antibiotic resistance marker and galK, that efficiently rendered N. gonorrhoeae susceptible to growth on 2-DOG. We then adapted the galK-based counter-selection and the use of 2-DOG for markerless mutagenesis, and applied biochemical and phenotypic analyses to confirm the absence of target genes. We show that our markerless mutagenesis method for N. gonorrhoeae has a high success rate, and should be a valuable gene editing tool in the future.

Mobility of extracellular DNA within gonococcal colonies

Transformation enables bacteria to acquire genetic information from extracellular DNA (eDNA). Close proximity between bacteria in colonies and biofilms may inhibit escape of eDNA from the colony but it also hinders its diffusion between donor and recipient. In this study, we investigate the mobility of DNA within colonies formed by Neisseria gonorrhoeae, and relate it to transformation efficiency. We characterize the penetration dynamics of fluorescent DNA into the colony at a time scale of hours and find that 300 bp fragments diffuse through the colony without hindrance. For DNA length exceeding 3 kbp, a concentration gradient between the edge and the center of the colony develops, indicating hindered diffusion. Accumulation of DNA within the colony increases with increasing DNA length. The presence of the gonococcal DNA uptake sequence (DUS), which mediates specific binding to type 4 pili (T4P) and uptake into the cell, steepens the radial concentration gradient within the colony, suggesting that the DUS reduces DNA mobility. In particular, DNA of N. gonorrhoeae containing multiple DUS is trapped at the periphery. Under conditions, where DUS containing DNA fragments readily enter the colony center, we investigate the efficiency of transformation. We show that despite rapid diffusion of DNA, the transformation is limited to the edge of young colonies. We conclude that DNA mobility depends on DNA length and specific binding mediated by the DUS, resulting in restricted mobility of gonococcal DNA. Yet gonococcal colonies accumulate DNA, and may therefore act as a reservoir for eDNA.

•

DNA fragments encompassing the length of a typical operon efficiently penetrate bacterial colonies.

•

Bacterial colonies accumulate eDNA with an efficiency that depends on the length and the DNA uptake sequence.

•

Genomic DNA from a distinct species spreads efficiently through gonococcal colonies, while gonococcal DNA and DNA from a closely related species are trapped.

•

Transformation is most efficient at the periphery of freshly assembled gonococcal colonies.

Recent advances in understanding and combatting Neisseria gonorrhoeae: a genomic perspective

The sexually transmitted infection (STI) gonorrhoea remains a major global public health concern. The World Health Organization (WHO) estimates that 87 million new cases in individuals who were 15 to 49 years of age occurred in 2016. The growing number of gonorrhoea cases is concerning given the rise in gonococci developing antimicrobial resistance (AMR). Therefore, a global action plan is needed to facilitate surveillance. Indeed, the WHO has made surveillance leading to the elimination of STIs (including gonorrhoea) a global health priority. The availability of whole genome sequence data offers new opportunities to combat gonorrhoea. This can be through (i) enhanced surveillance of the global prevalence of AMR, (ii) improved understanding of the population biology of the gonococcus, and (iii) opportunities to mine sequence data in the search for vaccine candidates. Here, we review the current status in Neisseria gonorrhoeae genomics. In particular, we explore how genomics continues to advance our understanding of this complex pathogen.

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

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

A community-driven resource for genomic epidemiology and antimicrobial resistance prediction of Neisseria gonorrhoeae at Pathogenwatch

BACKGROUND

Antimicrobial-resistant (AMR) Neisseria gonorrhoeae is an urgent threat to public health, as strains resistant to at least one of the two last-line antibiotics used in empiric therapy of gonorrhoea, ceftriaxone and azithromycin, have spread internationally. Whole genome sequencing (WGS) data can be used to identify new AMR clones and transmission networks and inform the development of point-of-care tests for antimicrobial susceptibility, novel antimicrobials and vaccines. Community-driven tools that provide an easy access to and analysis of genomic and epidemiological data is the way forward for public health surveillance.

METHODS

Here we present a public health-focussed scheme for genomic epidemiology of N. gonorrhoeae at Pathogenwatch ( https://pathogen.watch/ngonorrhoeae ). An international advisory group of experts in epidemiology, public health, genetics and genomics of N. gonorrhoeae was convened to inform on the utility of current and future analytics in the platform. We implement backwards compatibility with MLST, NG-MAST and NG-STAR typing schemes as well as an exhaustive library of genetic AMR determinants linked to a genotypic prediction of resistance to eight antibiotics. A collection of over 12,000 N. gonorrhoeae genome sequences from public archives has been quality-checked, assembled and made public together with available metadata for contextualization.

RESULTS

AMR prediction from genome data revealed specificity values over 99% for azithromycin, ciprofloxacin and ceftriaxone and sensitivity values around 99% for benzylpenicillin and tetracycline. A case study using the Pathogenwatch collection of N. gonorrhoeae public genomes showed the global expansion of an azithromycin-resistant lineage carrying a mosaic mtr over at least the last 10 years, emphasising the power of Pathogenwatch to explore and evaluate genomic epidemiology questions of public health concern.

CONCLUSIONS

The N. gonorrhoeae scheme in Pathogenwatch provides customised bioinformatic pipelines guided by expert opinion that can be adapted to public health agencies and departments with little expertise in bioinformatics and lower-resourced settings with internet connection but limited computational infrastructure. The advisory group will assess and identify ongoing public health needs in the field of gonorrhoea, particularly regarding gonococcal AMR, in order to further enhance utility with modified or new analytic methods.