Genomic analyses of Neisseria gonorrhoeae reveal an association of the gonococcal genetic island with antimicrobial resistance

penA profile of Neisseria gonorrhoeae in Guangdong, China: Novel penA alleles are related to decreased susceptibility to ceftriaxone or cefixime

BACKGROUND

Resistance to extended-spectrum cephalosporins (ESCs) has become a public health concern with the spread of Neisseria gonorrhoeae and increasing antimicrobial resistance. Mutation of penA, encoding penicillin-binding protein 2, represents a mechanism of ESC resistance. This study sought to assess penA alleles and mutations associated with decreased susceptibility (DS) to ESCs in N. gonorrhoeae.

MATERIALS AND METHODS

In 2021, 347 gonococci were collected in Guangdong, China. Minimum inhibitory concentations (MICs) of ceftriaxone and cefixime were determined, and whole-genome sequencing and phylogenetic analysis were performed. Multi-locus sequence typing (MLST) and conventional resistance determinants such as penA, mtrR, PonA and PorB were analysed. penA was genotyped and sequence-aligned using PubMLST.

RESULTS

Genome-wide phylogenetic analysis revealed that the prevalence of DS to ESCs was highest in Clade 11.1 (100.0%), Clade 2 (66.7%) and Clade 0 (55.7%), and the leading cause was strains with penA-60.001 or new penA alleles in clades. The penA phylogenetic tree is divided into two branches: non-mosaic penA and mosaic penA. The latter contained penA-60.001, penA-10 and penA-34. penA profile analysis indicated that A311V and T483S are closely related to DS to ESCs in mosaic penA. The new alleles NEIS1753_2840 and NEIS1753_2837 are closely related to penA-60.001, with DS to ceftriaxone and cefixime of 100%. NEIS1753_2660, a derivative of penA-10 (A486V), has increased DS to ceftriaxone. NEIS1753_2846, a derivative of penA-34.007 (G546S), has increased DS to cefixime.

CONCLUSION

This study identified critical penA alleles related to elevated MICs, and trends of gonococcus-evolved mutated penA associated with DS to ESCs in Guangdong.

Dinner date: Neisseria gonorrhoeae central carbon metabolism and pathogenesis

Neisseria gonorrhoeae, the causative agent of the sexually transmitted infection gonorrhea, is a human-adapted pathogen that does not productively infect other organisms. The ongoing relationship between N. gonorrhoeae and the human host is facilitated by the exchange of nutrient resources that allow for N. gonorrhoeae growth in the human genital tract. What N. gonorrhoeae 'eats' and the pathways used to consume these nutrients have been a topic of investigation over the last 50 years. More recent investigations are uncovering the impact of N. gonorrhoeae metabolism on infection and inflammatory responses, the environmental influences driving N. gonorrhoeae metabolism, and the metabolic adaptations enabling antimicrobial resistance. This mini-review is an introduction to the field of N. gonorrhoeae central carbon metabolism in the context of pathogenesis. It summarizes the foundational work used to characterize N. gonorrhoeae central metabolic pathways and the effects of these pathways on disease outcomes, and highlights some of the most recent advances and themes under current investigation. This review ends with a brief description of the current outlook and technologies under development to increase understanding of how the pathogenic potential of N. gonorrhoeae is enabled by metabolic adaptation.

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.

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

Introduction

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

Methods

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

Results and Discussion

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

Gonococcal Genetic Island in the Global Neisseria gonorrhoeae Population: A Model of Genetic Diversity and Association with Resistance to Antimicrobials

The aim of this work was to study the genetic diversity of the gonococcal genetic island (GGI) responsible for the type IV secretion system (T4SS) and the association of a functionally active GGI with antimicrobial resistance. An analysis of the GGI in a sample of 14,763 genomes of N. gonorrhoeae isolates from the Pathogenwatch database collected in 1996-2019 from 68 countries was performed. A model of GGI's genetic diversity that divides the global gonococcal population into fifty-one GGI clusters and three GGI superclusters based on the allele type of the traG gene and substitutions of the atlA and ych genes for eppA and ych1 has been proposed, reflecting differences among isolates in the T4SS functionality. The NG-MAST and MLST typing schemes (with accuracies of 91% and 83%, respectively) allowed the determination of both the presence of a GGI and the GGI cluster and, correspondingly, the structure of the GGI and the ability to secrete DNA. A statistically significant difference in the proportion of N. gonorrhoeae isolates resistant to ciprofloxacin, cefixime, tetracycline, and penicillin was found when comparing populations with a functional and a non-functional GGI. The presence of a functional GGI did not affect the proportion of azithromycin-resistant isolates.

Molecular Mechanisms of Drug Resistance and Epidemiology of Multidrug-Resistant Variants of Neisseria gonorrhoeae

The paper presents various issues related to the increasing drug resistance of Neisseria gonorrhoeae and the occurrence and spread of multidrug-resistant clones. One of the most important is the incidence and evolution of resistance mechanisms of N. gonorrhoeae to beta-lactam antibiotics. Chromosomal resistance to penicillins and oxyimino-cephalosporins and plasmid resistance to penicillins are discussed. Chromosomal resistance is associated with the presence of mutations in the PBP2 protein, containing mosaic variants and nonmosaic amino acid substitutions in the transpeptidase domain, and their correlation with mutations in the mtrR gene and its promoter regions (the MtrCDE membrane pump repressor) and in several other genes, which together determine reduced sensitivity or resistance to ceftriaxone and cefixime. Plasmid resistance to penicillins results from the production of beta-lactamases. There are different types of beta-lactamases as well as penicillinase plasmids. In addition to resistance to beta-lactam antibiotics, the paper covers the mechanisms and occurrence of resistance to macrolides (azithromycin), fluoroquinolones and some other antibiotics. Moreover, the most important epidemiological types of multidrug-resistant N. gonorrhoeae, prevalent in specific years and regions, are discussed. Epidemiological types are defined as sequence types, clonal complexes and genogroups obtained by various typing systems such as NG-STAR, NG-MAST and MLST. New perspectives on the treatment of N. gonorrhoeae infections are also presented, including new drugs active against multidrug-resistant strains.

Genomic sequence of the non-pathogen Neisseria sp. strain MA1-1 with antibiotic resistance and virulence factors isolated from a head and neck cancer patient

Recent research has claimed virulence factors or antimicrobial resistance in commensal or non-pathogenic Neisseria spp. This study aimed to isolate and analyze commensal microorganisms related to the genus Neisseria from the oral cavity of a patient with head and neck cancer. We successfully isolated strain MA1-1 and identified its functional gene contents. Although strain MA1-1 was related to Neisseria flava based on 16S rRNA gene sequence similarity, genomic relatedness analysis revealed that strain MA1-1 was closely related to Neisseria mucosa, reported as a commensal Neisseria species. The strain MA1-1 genome harbored genes for microaerobic respiration and the complete core metabolic pathway with few transporters for nutrients. A number of genes have been associated with virulence factors and resistance to various antibiotics. In addition, the comparative genomic analysis showed that most genes identified in the strain MA1-1 were shared with other Neisseria spp. including two well-known pathogens, Neisseria gonorrhoeae and Neisseria meningitidis. This indicates that the gene content of intra-members of the genus Neisseria has been evolutionarily conserved and is stable, with no gene recombination with other microbes in the host. Finally, this study provides more fundamental interpretations for the complete gene sequence of commensal Neisseria spp. and will contribute to advancing public health knowledge.

Antibiotic Resistance in Neisseria gonorrhoeae: Challenges in Research and Treatment

Gonococcal infection caused by the Gram-negative bacteria Neisseria gonorrhoeae is one of the most common sexually transmitted infections (STIs) worldwide [...]

Comparative Whole-Genome Analysis of Neisseria gonorrhoeae Isolates Revealed Changes in the Gonococcal Genetic Island and Specific Genes as a Link to Antimicrobial Resistance

Comparative whole-genome analysis was performed for Neisseria gonorrhoeae isolates belonging to the Neisseria gonorrhoeae multiantigen sequence typing (NG-MAST) types predominant worldwide — 225, 1407, 2400, 2992, and 4186 — and to genogroup 807, the most common genogroup in the Russian Federation. Here, for the first time, the complete genomes of 25 N. gonorrhoeae isolates from genogroup 807 were obtained. For NG-MAST types 225, 1407, 2400, 2992, and 4186, genomes from the Pathogenwatch database were used. The phylogenetic network constructed for 150 genomes showed that the clustering according to NG-MAST type corresponded to the clustering according to genome. Comparisons of genomes of the six sequence types revealed 8-20 genes specific to each sequence type, including the loci for phase variations and genetic components of the gonococcal genetic island (GGI). NG-MAST type 2992 and 4186 isolates either lacked the GGI or carried critical mutations in genes essential for DNA secretion. In all analyzed genogroup 807 isolates, substitution of the essential atlA gene with the eppA gene was found, accompanied by a change in the traG allele, replacement of the ych gene with ych1, and the absence of the exp1 gene, which is likely to result in loss of GGI functionality. For the NG-MAST type 225, 1407 and 2400 isolates, no premature stop codons or reading frameshifts were found in the genes essential for GGI function. A relationship between isolate susceptibility to ciprofloxacin, penicillin, tetracycline and the presence of lesions in GGI genes necessary for DNA secretion was established. The N. gonorrhoeae evolutionary pathways, which allow a particular sequence type to maintain long-term predominance in the population, may include changes in genes responsible for adhesion and virulence, changes in the GGI structure, preservation of genes carrying drug resistance determinants, and changes in genes associated with host adaptation or encoding enzymes of biochemical pathways.

Transcriptional and Translational Responsiveness of the Neisseria gonorrhoeae Type IV Secretion System to Conditions of Host Infections

The type IV secretion system of Neisseria gonorrhoeae translocates single-stranded DNA into the extracellular space, facilitating horizontal gene transfer and initiating biofilm formation. Expression of this system has been observed to be low under laboratory conditions, and multiple levels of regulation have been identified. We used a translational fusion of lacZ to traD, the gene for the type IV secretion system coupling protein, to screen for increased type IV secretion system expression. We identified several physiologically relevant conditions, including surface adherence, decreased manganese or iron, and increased zinc or copper, which increase gonococcal type IV secretion system protein levels through transcriptional and/or translational mechanisms. These metal treatments are reminiscent of the conditions in the macrophage phagosome. The ferric uptake regulator, Fur, was found to repress traD transcript levels but to also have a second role, acting to allow TraD protein levels to increase only in the absence of iron. To better understand type IV secretion system regulation during infection, we examined transcriptomic data from active urethral infection samples from five men. The data demonstrated differential expression of 20 of 21 type IV secretion system genes during infection, indicating upregulation of genes necessary for DNA secretion during host infection.

Hydrogel Droplet Microarray for Genotyping Antimicrobial Resistance Determinants in Neisseria gonorrhoeae Isolates

A multiplex assay based on a low-density hydrogel microarray was developed to identify genomic substitutions in N. gonorrhoeae that determine resistance to the currently recommended treatment agents ceftriaxone and azithromycin and the previously used drugs penicillin, tetracycline, and ciprofloxacin. The microarray identifies 74 drug resistance determinants in the N. gonorrhoeae penA, ponA, porB, gyrA, parC, rpsJ, mtrR, bla_TEM, tetM, and 23S rRNA genes. The hydrogel elements were formed by automated dispensing of nanoliter-volume droplets followed by UV-induced copolymerization of NH₂-containing oligonucleotides with gel-forming monomers. Polybutylene terephthalate plates without special modifications were used as microarray substrates. Sequences and concentrations of immobilized oligonucleotides, gel composition, and hybridization conditions were carefully selected, and the median discrimination ratio ranged from 2.8 to 29.4, allowing unambiguous identification of single-nucleotide substitutions. The mutation identification results in a control sample of 180 N. gonorrhoeae isolates were completely consistent with the Sanger sequencing results. In total, 648 clinical N. gonorrhoeae isolates obtained in Russia during the last 5 years were analyzed and genotyped using these microarrays. The results allowed us to draw conclusions about the present situation with antimicrobial susceptibility of N. gonorrhoeae in Russia and demonstrated the possibility of using hydrogel microarrays to control the spread of antibiotic resistance.

Antimicrobial Resistance Mechanisms, Multilocus Sequence Typing, and NG-STAR Sequence Types of Diverse Neisseria gonorrhoeae Isolates in KwaZulu-Natal, South Africa

Antimicrobial resistance (AMR) is a major challenge to managing infectious diseases. Africa has the highest incidence of gonorrhoea, but there is a lack of comprehensive data from sparse surveillance programs. This study investigated the molecular epidemiology and AMR profiles of Neisseria gonorrhoeae isolates in KwaZulu-Natal province (KZN), South Africa. Repository isolates from patients attending public health care clinics for sexually transmitted infection (STI) care were used for phenotypic and genotypic analysis. An Etest was performed to determine antimicrobial susceptibility. Whole-genome sequencing (WGS) was used to determine epidemiology and to predict susceptibility by detecting resistance-associated genes and mutations. Among the 61 isolates, multiple sequence types were identified. Six isolates were novel, as determined by multilocus sequence typing. N. gonorrhoeae sequence typing for antimicrobial resistance (NG-STAR) determined 48 sequence types, of which 35 isolates had novel antimicrobial profiles. Two novel penA alleles and eight novel mtrR alleles were identified. Point mutations were detected in gyrA, parC, mtrR, penA, ponA, and porB1. This study revealed a high prevalence of AMR (penicillin 67%, tetracycline 89%, and ciprofloxacin 52%). However, spectinomycin, cefixime, ceftriaxone, and azithromycin remained 100% effective. This study is one of the first to comprehensively describe the epidemiology and AMR of N. gonorrhoeae in KZN, South Africa and Africa, using WGS. KZN has a wide strain diversity and most of these sequence types have been detected in multiple countries; however, more than half of our isolates have novel antimicrobial profiles. Continued surveillance is crucial to monitor the emergence of resistance to cefixime, ceftriaxone, and azithromycin.

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.

Antimicrobial Resistance Profiling and Phylogenetic Analysis of Neisseria gonorrhoeae Clinical Isolates From Kenya in a Resource-Limited Setting

Background

Africa has one of the highest incidences of gonorrhea. Neisseria gonorrhoeae is gaining resistance to most of the available antibiotics, compromising treatment across the world. Whole-genome sequencing (WGS) is an efficient way of predicting AMR determinants and their spread in the population. Recent advances in next-generation sequencing technologies like Oxford Nanopore Technology (ONT) have helped in the generation of longer reads of DNA in a shorter duration with lower cost. Increasing accuracy of base-calling algorithms, high throughput, error-correction strategies, and ease of using the mobile sequencer MinION in remote areas lead to its adoption for routine microbial genome sequencing. To investigate whether MinION-only sequencing is sufficient for WGS and downstream analysis in resource-limited settings, we sequenced the genomes of 14 suspected N. gonorrhoeae isolates from Nairobi, Kenya.

Methods

Using WGS, the isolates were confirmed to be cases of N. gonorrhoeae (n = 9), and there were three co-occurrences of N. gonorrhoeae with Moraxella osloensis and N. meningitidis (n = 2). N. meningitidis has been implicated in sexually transmitted infections in recent years. The near-complete N. gonorrhoeae genomes (n = 10) were analyzed further for mutations/factors causing AMR using an in-house database of mutations curated from the literature.

Results

We observe that ciprofloxacin resistance is associated with multiple mutations in both gyrA and parC. Mutations conferring tetracycline (rpsJ) and sulfonamide (folP) resistance and plasmids encoding beta-lactamase were seen in all the strains, and tet(M)-containing plasmids were identified in nine strains. Phylogenetic analysis clustered the 10 isolates into clades containing previously sequenced genomes from Kenya and countries across the world. Based on homology modeling of AMR targets, we see that the mutations in GyrA and ParC disrupt the hydrogen bonding with quinolone drugs and mutations in FolP may affect interaction with the antibiotic.

Conclusion

Here, we demonstrate the utility of mobile DNA sequencing technology in producing a consensus genome for sequence typing and detection of genetic determinants of AMR. The workflow followed in the study, including AMR mutation dataset creation and the genome identification, assembly, and analysis, can be used for any clinical isolate. Further studies are required to determine the utility of real-time sequencing in outbreak investigations, diagnosis, and management of infections, especially in resource-limited settings.

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.

Emergence and evolution of antimicrobial resistance genes and mutations in Neisseria gonorrhoeae

BACKGROUND

Antimicrobial resistance in Neisseria gonorrhoeae is a global health concern. Strains from two internationally circulating sequence types, ST-7363 and ST-1901, have acquired resistance to third-generation cephalosporins, mainly due to mosaic penA alleles. These two STs were first detected in Japan; however, the timeline, mechanism, and process of emergence and spread of these mosaic penA alleles to other countries remain unknown.

METHODS

We studied the evolution of penA alleles by obtaining the complete genomes from three Japanese ST-1901 clinical isolates harboring mosaic penA allele 34 (penA-34) dating from 2005 and generating a phylogenetic representation of 1075 strains sampled from 35 countries. We also sequenced the genomes of 103 Japanese ST-7363 N. gonorrhoeae isolates from 1996 to 2005 and reconstructed a phylogeny including 88 previously sequenced genomes.

RESULTS

Based on an estimate of the time-of-emergence of ST-1901 (harboring mosaic penA-34) and ST-7363 (harboring mosaic penA-10), and > 300 additional genome sequences of Japanese strains representing multiple STs isolated in 1996-2015, we suggest that penA-34 in ST-1901 was generated from penA-10 via recombination with another Neisseria species, followed by recombination with a gonococcal strain harboring wildtype penA-1. Following the acquisition of penA-10 in ST-7363, a dominant sub-lineage rapidly acquired fluoroquinolone resistance mutations at GyrA 95 and ParC 87-88, by independent mutations rather than horizontal gene transfer. Data in the literature suggest that the emergence of these resistance determinants may reflect selection from the standard treatment regimens in Japan at that time.

CONCLUSIONS

Our findings highlight how antibiotic use and recombination across and within Neisseria species intersect in driving the emergence and spread of drug-resistant gonorrhea.

Meningococcal Deduced Vaccine Antigen Reactivity (MenDeVAR) Index: a Rapid and Accessible Tool That Exploits Genomic Data in Public Health and Clinical Microbiology Applications

As microbial genomics makes increasingly important contributions to clinical and public health microbiology, the interpretation of whole-genome sequence data by nonspecialists becomes essential. In the absence of capsule-based vaccines, two protein-based vaccines have been used for the prevention of invasive serogroup B meningococcal disease (IMD) since their licensure in 2013 and 2014. These vaccines have different components and different levels of coverage of meningococcal variants. Hence, decisions regarding which vaccine to use in managing serogroup B IMD outbreaks require information about the index case isolate, including (i) the presence of particular vaccine antigen variants, (ii) the expression of vaccine antigens, and (iii) the likely susceptibility of its antigen variants to antibody-dependent bactericidal killing.

As microbial genomics makes increasingly important contributions to clinical and public health microbiology, the interpretation of whole-genome sequence data by nonspecialists becomes essential. In the absence of capsule-based vaccines, two protein-based vaccines have been used for the prevention of invasive serogroup B meningococcal disease (IMD) since their licensure in 2013 and 2014. These vaccines have different components and different levels of coverage of meningococcal variants. Hence, decisions regarding which vaccine to use in managing serogroup B IMD outbreaks require information about the index case isolate, including (i) the presence of particular vaccine antigen variants, (ii) the expression of vaccine antigens, and (iii) the likely susceptibility of its antigen variants to antibody-dependent bactericidal killing. To obtain this information requires a multitude of laboratory assays, impractical in real-time clinical settings, where the information is most urgently needed. To facilitate assessment for public health and clinical purposes, we synthesized genomic and experimental data from published sources to develop and implement the Meningococcal Deduced Vaccine Antigen Reactivity (MenDeVAR) Index, which is publicly available on PubMLST (https://pubmlst.org). Using whole-genome sequences or individual gene sequences obtained from IMD isolates or clinical specimens, the MenDeVAR Index provides rapid evidence-based information on the presence and possible immunological cross-reactivity of different meningococcal vaccine antigen variants. The MenDeVAR Index enables practitioners who are not genomics specialists to assess the likely reactivity of vaccines for individual cases, outbreak management, or the assessment of public health vaccine programs. The MenDeVAR Index has been developed in consultation with, but independently of, both the 4CMenB (Bexsero; GSK) and rLP2086 (Trumenba; Pfizer, Inc.) vaccine manufacturers.

Increased power from conditional bacterial genome-wide association identifies macrolide resistance mutations in Neisseria gonorrhoeae

The emergence of resistance to azithromycin complicates treatment of Neisseria gonorrhoeae, the etiologic agent of gonorrhea. Substantial azithromycin resistance remains unexplained after accounting for known resistance mutations. Bacterial genome-wide association studies (GWAS) can identify novel resistance genes but must control for genetic confounders while maintaining power. Here, we show that compared to single-locus GWAS, conducting GWAS conditioned on known resistance mutations reduces the number of false positives and identifies a G70D mutation in the RplD 50S ribosomal protein L4 as significantly associated with increased azithromycin resistance (p-value = 1.08 × 10-11). We experimentally confirm our GWAS results and demonstrate that RplD G70D and other macrolide binding site mutations are prevalent (present in 5.42% of 4850 isolates) and widespread (identified in 21/65 countries across two decades). Overall, our findings demonstrate the utility of conditional associations for improving the performance of microbial GWAS and advance our understanding of the genetic basis of macrolide resistance.

Antimicrobial Resistance of Neisseria gonorrhoeae Isolates from High-Risk Men in Johannesburg, South Africa

Neisseria gonorrhoeae antimicrobial drug resistance has emerged worldwide; however, the situation in sub-Saharan Africa is not well documented. We investigated the molecular epidemiology and occurrence of antimicrobial resistance in Neisseria gonorrhoeae infections in two core transmission groups of men in Johannesburg, South Africa. We recruited men who have sex with men (MSM) presenting with urethral discharge and men with recurrent episodes of urethral discharge. Molecular testing and culture for N. gonorrhoeae were performed, followed by antimicrobial susceptibility testing. Whole-genome sequencing (WGS) was used to identify resistance-conferring mutations and to determine the genetic relatedness of the isolates. In all, 51 men were recruited; 42 (82%) had N. gonorrhoeae infections. Most gonococcal isolates were resistant to ciprofloxacin (78%) and tetracycline (74%); 33% were penicillin resistant. All gonococcal isolates were susceptible to cephalosporins and spectinomycin. Azithromycin resistance was observed in 4 (15%) isolates (epidemiological cutoff), all with mutations in the mtrR promoter region. Most of the isolates (19/27) harbored the gonococcal genetic island, which is associated with antimicrobial resistance. WGS revealed a diverse epidemic with mostly novel NG-STAR (70%) and NG-MAST (70%) sequence types. Thus, we demonstrate a high prevalence of antimicrobial resistance in Neisseria gonorrhoeae strains obtained from high-risk men in South Africa. The introduction of diagnostics and scale-up of surveillance are warranted to prevent the emergence of multidrug-resistant infections.

Identification of integrative and conjugative elements in pathogenic and commensal Neisseriaceae species via genomic distributions of DNA uptake sequence dialects

Mobile genetic elements (MGEs) are key factors responsible for dissemination of virulence determinants and antimicrobial-resistance genes amongst pathogenic bacteria. Conjugative MGEs are notable for their high gene loads donated per transfer event, broad host ranges and phylogenetic ubiquity amongst prokaryotes, with the subclass of chromosomally inserted integrative and conjugative elements (ICEs) being particularly abundant. The focus on a small number of model systems has biased the study of ICEs towards those conferring readily selectable phenotypes to host cells, whereas the identification and characterization of integrated cryptic elements remains challenging. Even though antimicrobial resistance and horizontally acquired virulence genes are major factors aggravating neisserial infection, conjugative MGEs of Neisseria gonorrhoeae and Neisseria meningitidis remain poorly characterized. Using a phenotype-independent approach based on atypical distributions of DNA uptake sequences (DUSs) in MGEs relative to the chromosomal background, we have identified two groups of chromosomally integrated conjugative elements in Neisseria: one found almost exclusively in pathogenic species possibly deriving from the genus Kingella, the other belonging to a group of Neisseria mucosa-like commensals. The former element appears to enable transfer of traditionally gonococcal-specific loci such as the virulence-associated toxin-antitoxin system fitAB to N. meningitidis chromosomes, whilst the circular form of the latter possesses a unique attachment site (attP) sequence seemingly adapted to exploit DUS motifs as chromosomal integration sites. In addition to validating the use of DUS distributions in Neisseriaceae MGE identification, the >170 identified ICE sequences provide a valuable resource for future studies of ICE evolution and host adaptation.

Evidence of Recent Genomic Evolution in Gonococcal Strains With Decreased Susceptibility to Cephalosporins or Azithromycin in the United States, 2014-2016

BACKGROUND

Given the lack of new antimicrobials or a vaccine, understanding the evolutionary dynamics of Neisseria gonorrhoeae is a significant public and global health priority. We investigated the emergence and spread of gonococcal strains with decreased susceptibility to cephalosporins and azithromycin using detailed genomic analyses of gonococcal isolates collected in the United States, 2014-2016.

METHODS

We sequenced genomes of 649 isolates collected through the Gonococcal Isolate Surveillance Project. We examined the genetic relatedness of isolates and assessed associations between clades and various genotypic and phenotypic combinations.

RESULTS

We identified a large and clonal lineage of strains (MLST ST9363) associated with elevated azithromycin minimum inhibitory concentration (AZIem), characterized by a mosaic mtr locus (C substitution in the mtrR promoter, mosaic mtrR and mtrD). Mutations in 23S rRNA were sporadically distributed among AZIem strains. Another clonal group (MLST ST1901) possessed 7 unique PBP2 patterns, and it shared common mutations in other genes associated with cephalosporin resistance.

CONCLUSIONS

Whole-genome sequencing methods can enhance monitoring of antimicrobial resistant gonococcal strains by identifying gonococcal populations containing mutations of concern. These methods could inform the development of point-of-care diagnostic tests designed to determine the specific antibiotic susceptibility profile of a gonococcal infection in a patient.

Genomic evolution of Neisseria gonorrhoeae since the preantibiotic era (1928-2013): antimicrobial use/misuse selects for resistance and drives evolution

BACKGROUND

Multidrug-resistant Neisseria gonorrhoeae strains are prevalent, threatening gonorrhoea treatment globally, and understanding of emergence, evolution, and spread of antimicrobial resistance (AMR) in gonococci remains limited. We describe the genomic evolution of gonococci and their AMR, related to the introduction of antimicrobial therapies, examining isolates from 1928 (preantibiotic era) to 2013 in Denmark. This is, to our knowledge, the oldest gonococcal collection globally.

METHODS

Lyophilised isolates were revived and examined using Etest (18 antimicrobials) and whole-genome sequencing (WGS). Quality-assured genome sequences were obtained for 191 viable and 40 non-viable isolates and analysed with multiple phylogenomic approaches.

RESULTS

Gonococcal AMR, including an accumulation of multiple AMR determinants, started to emerge particularly in the 1950s-1970s. By the twenty-first century, resistance to most antimicrobials was common. Despite that some AMR determinants affect many physiological functions and fitness, AMR determinants were mainly selected by the use/misuse of gonorrhoea therapeutic antimicrobials. Most AMR developed in strains belonging to one multidrug-resistant (MDR) clade with close to three times higher genomic mutation rate. Modern N. gonorrhoeae was inferred to have emerged in the late-1500s and its genome became increasingly conserved over time.

CONCLUSIONS

WGS of gonococci from 1928 to 2013 showed that no AMR determinants, except penB, were in detectable frequency before the introduction of gonorrhoea therapeutic antimicrobials. The modern gonococcus is substantially younger than previously hypothesized and has been evolving into a more clonal species, driven by the use/misuse of antimicrobials. The MDR gonococcal clade should be further investigated for early detection of strains with predispositions to develop and maintain MDR and for initiation of public health interventions.

Mosaic structure of the penA gene in the oropharynx of men who have sex with men negative for gonorrhoea

The oropharynx represents a crucial site for the emergence of multi-drug resistance in Neisseria gonorrhoeae. The mosaic penA alleles, associated with decreased susceptibility to cephalosporins, have emerged by DNA recombination with partial penA genes, particularly those from commensal pharyngeal Neisseria species. Here, we investigated the prevalence of the mosaic structure of the penA gene in the oropharynx of men who have sex with men testing negative for pharyngeal gonorrhoea. From January 2016 to June 2018, 351 gonorrhoea-negative men who have sex with men attending a sexually transmitted infection clinic in Italy were enrolled. Pharyngeal swabs underwent a real-time polymerase chain reaction (PCR) for the detection of the mosaic penA gene. In case of positivity, PCR products were sequenced and searched against several sequences of Neisseria strains. Overall, 31 patients (8.8%) were found positive for the presence of the mosaic penA gene. The positivity was significantly associated with previous cases of pharyngeal gonorrhoea (relative risk [RR]: 3.56, 95% confidence interval 1.44–8.80) and with recent exposure to beta-lactams (RR: 4.29, 95% confidence interval 2.20–8.38). All penA-positive samples showed a high relatedness (90–99%) with mosaic-positive Neisseria strains. Our data underline that commensal Neisseria species of the oropharynx may be a significant reservoir for genetic material conferring antimicrobial resistance in N. gonorrhoeae.

Characterization of antimicrobial resistance genes from Neisseria gonorrhoeae positive remnant Aptima urine specimens

Aim: To ascertain the antimicrobial resistance and strain types (STs) of Neisseria gonorrhoeae from 50 remnant Aptima urine specimens using molecular methods. Methods: Mutations predictive of resistance to six antibiotics were identified in eight genes. STs were determined using NG-MAST and NG-STAR. Results: All eight antimicrobial resistance genes could be characterized in 36 specimens. A total of 17 specimens were predicted to be susceptible to all antibiotics, including ceftriaxone. Decreased susceptibility to cefixime and ciprofloxacin resistance was predicted in 11 specimens (PBP2 type 34.001). Overall, 38/50 specimens were predicted to be ciprofloxacin susceptible; three were azithromycin resistant. Nineteen NG-MAST and 21 NG-STAR STs were noted. Conclusion: Molecular analysis of remnant Aptima specimens enabled the prediction of emerging gonococcal cefixime and azithromycin resistance which would otherwise have been undetected.

Genomic Characterization of Gonococci from Different Anatomic Sites, Italy, 2007-2014

In recent decades, Neisseria gonorrhoeae has developed resistance to several antimicrobial classes. Molecular epidemiology approaches are useful for detecting emerging, often resistant, gonococcal clones. In this study, 67 N. gonorrhoeae isolates from different anatomic sites, collected over 8 years in Italy, were analyzed by whole genome sequencing (WGS). WGS was performed using the Illumina NextSeq 500 platform. Phylogenetic analysis was based on core single nucleotide polymorphism (SNP) and core genome multilocus sequence typing (cgMLST). N. gonorrhoeae multi-antigen sequence typing (NG-MAST), MLST, and N. gonorrhoeae sequence typing for antimicrobial resistance (NG-STAR) were carried out in silico using WGS data. Antimicrobial susceptibility against a four-drug panel was evaluated using a gradient diffusion method. Overall, gonococci clustered in accordance with NG-MAST, MLST, NG-STAR, and antimicrobials susceptibility profiles, but not with the site of isolation, HIV status, and patient sexual orientation. Phylogenetic analysis identified nine clades: two of them were the predominant and including gonococci of G1407 and G2400 genogroups.

Secretion of Chromosomal DNA by the Neisseria gonorrhoeae Type IV Secretion System

Approximately 80% of Neisseria gonorrhoeae and 17.5% of Neisseria meningitidis clinical isolates carry a ~59 kb genomic island known as the gonococcal genetic island (GGI). About half of the GGI consists of genes encoding a type IV secretion system (T4SS), and most of these genes are clustered in a ~28 kb region at one end of the GGI. Two additional genes (parA and parB) are found at the other end of the island. The remainder of the GGI consists mostly of hypothetical proteins, with several being identified as DNA-binding or DNA-processing proteins. The T4SS genes show similarity to those of the F-plasmid family of conjugation systems, with similarity in gene order and a low but significant level of sequence identity for the encoded proteins. However, several GGI-encoded proteins are unique from the F-plasmid system, such as AtlA, Yag, and TraA. Interestingly, the gonococcal T4SS does not act as a conjugation system. Instead, this T4SS secretes ssDNA into the extracellular milieu, where it can serve to transform highly competent Neisseria species, thereby increasing the transfer of genetic information. Although many of the T4SS proteins are expressed at low levels, this system has been implicated in several cellular processes. The secreted ssDNA is involved in the initial stages of biofilm formation, and the presence of the T4SS enables TonB-independent intracellular survival of N. gonorrhoeae strains during infection of cervical cells. Other GGI-like T4SSs have been identified in several other α-, β-, and γ-proteobacteria, but the function of these GGI-like T4SSs is unknown. Remarkably, the presence of the GGI is related to resistance to several antibiotics. Here, we describe our current knowledge about the GGI and its unique ssDNA-secreting T4SS.

Closed Genome Sequences of Clinical Neisseria gonorrhoeae Strains Obtained from Combined Oxford Nanopore and Illumina Sequencing

Neisseria gonorrhoeae is the etiological agent of gonorrhea, the second most common notifiable disease in the United States. Here, we used a hybrid approach combining Oxford Nanopore Technologies MinION and Illumina MiSeq sequencing data to obtain closed genome sequences of nine clinical N. gonorrhoeae isolates.

Neisseria gonorrhoeae is the etiological agent of gonorrhea, the second most common notifiable disease in the United States. Here, we used a hybrid approach combining Oxford Nanopore Technologies MinION and Illumina MiSeq sequencing data to obtain closed genome sequences of nine clinical N. gonorrhoeae isolates.

Genomic epidemiology and antimicrobial resistance of Neisseria gonorrhoeae in New Zealand

Background

Antimicrobial-resistant Neisseria gonorrhoeae is a major threat to public health. No studies to date have examined the genomic epidemiology of gonorrhoea in the Western Pacific Region, where the incidence of gonorrhoea is particularly high.

Methods

A population-level study of N. gonorrhoeae in New Zealand (October 2014 to May 2015). Comprehensive susceptibility testing and WGS data were obtained for 398 isolates. Relatedness was inferred using phylogenetic trees, and pairwise core SNPs. Mutations and genes known to be associated with resistance were identified, and correlated with phenotype.

Results

Eleven clusters were identified. In six of these clusters, >25% of isolates were from females, while in eight of them, >15% of isolates were from females. Drug resistance was common; 98%, 32% and 68% of isolates were non-susceptible to penicillin, ciprofloxacin and tetracycline, respectively. Elevated MICs to extended-spectrum cephalosporins (ESCs) were observed in 3.5% of isolates (cefixime MICs ≥ 0.12 mg/L, ceftriaxone MICs ≥ 0.06 mg/L). Only nine isolates had penA XXXIV genotypes, three of which had decreased susceptibility to ESCs (MIC = 0.12 mg/L). Azithromycin non-susceptibility was identified in 43 isolates (10.8%); two of these isolates had 23S mutations (C2611T, 4/4 alleles), while all had mutations in mtrR or its promoter.

Conclusions

The high proportion of females in clusters suggests transmission is not exclusively among MSM in New Zealand; re-assessment of risk factors for transmission may be warranted in this context. As elevated MICs of ESCs and/or azithromycin were found in closely related strains, targeted public health interventions to halt transmission are urgently needed.

Using Machine Learning To Predict Antimicrobial MICs and Associated Genomic Features for Nontyphoidal Salmonella

Nontyphoidal Salmonella species are the leading bacterial cause of foodborne disease in the United States. Whole-genome sequences and paired antimicrobial susceptibility data are available for Salmonella strains because of surveillance efforts from public health agencies. In this study, a collection of 5,278 nontyphoidal Salmonella genomes, collected over 15 years in the United States, was used to generate extreme gradient boosting (XGBoost)-based machine learning models for predicting MICs for 15 antibiotics. The MIC prediction models had an overall average accuracy of 95% within ±1 2-fold dilution step (confidence interval, 95% to 95%), an average very major error rate of 2.7% (confidence interval, 2.4% to 3.0%), and an average major error rate of 0.1% (confidence interval, 0.1% to 0.2%). The model predicted MICs with no a priori information about the underlying gene content or resistance phenotypes of the strains. By selecting diverse genomes for the training sets, we show that highly accurate MIC prediction models can be generated with less than 500 genomes. We also show that our approach for predicting MICs is stable over time, despite annual fluctuations in antimicrobial resistance gene content in the sampled genomes. Finally, using feature selection, we explore the important genomic regions identified by the models for predicting MICs. To date, this is one of the largest MIC modeling studies to be published. Our strategy for developing whole-genome sequence-based models for surveillance and clinical diagnostics can be readily applied to other important human pathogens.

Genome-wide discovery of epistatic loci affecting antibiotic resistance in Neisseria gonorrhoeae using evolutionary couplings

Genome analysis should allow the discovery of interdependent loci that together cause antibiotic resistance. In practice, however, the vast number of possible epistatic interactions erodes statistical power. Here, we extend an approach that has been successfully used to identify epistatic residues in proteins to infer genomic loci that are strongly coupled. This approach reduces the number of tests required for an epistatic genome-wide association study of antibiotic resistance and increases the likelihood of identifying causal epistasis. We discovered 38 loci and 240 epistatic pairs that influence the minimum inhibitory concentrations of 5 different antibiotics in 1,102 isolates of Neisseria gonorrhoeae that were confirmed in a second dataset of 495 isolates. Many known resistance-affecting loci were recovered; however, the majority of associations occurred in unreported genes, such as murE. About half of the discovered epistasis involved at least one locus previously associated with antibiotic resistance, including interactions between gyrA and parC. Still, many combinations involved unreported loci and genes. While most variation in minimum inhibitory concentrations could be explained by identified loci, epistasis substantially increased explained phenotypic variance. Our work provides a systematic identification of epistasis affecting antibiotic resistance in N. gonorrhoeae and a generalizable approach for epistatic genome-wide association studies.

Open-access bacterial population genomics: BIGSdb software, the PubMLST.org website and their applications

The PubMLST.org website hosts a collection of open-access, curated databases that integrate population sequence data with provenance and phenotype information for over 100 different microbial species and genera.  Although the PubMLST website was conceived as part of the development of the first multi-locus sequence typing (MLST) scheme in 1998 the software it uses, the Bacterial Isolate Genome Sequence database (BIGSdb, published in 2010), enables PubMLST to include all levels of sequence data, from single gene sequences up to and including complete, finished genomes.  Here we describe developments in the BIGSdb software made from publication to June 2018 and show how the platform realises microbial population genomics for a wide range of applications.  The system is based on the gene-by-gene analysis of microbial genomes, with each deposited sequence annotated and curated to identify the genes present and systematically catalogue their variation.  Originally intended as a means of characterising isolates with typing schemes, the synthesis of sequences and records of genetic variation with provenance and phenotype data permits highly scalable (whole genome sequence data for tens of thousands of isolates) means of addressing a wide range of functional questions, including: the prediction of antimicrobial resistance; likely cross-reactivity with vaccine antigens; and the functional activities of different variants that lead to key phenotypes.  There are no limitations to the number of sequences, genetic loci, allelic variants or schemes (combinations of loci) that can be included, enabling each database to represent an expanding catalogue of the genetic variation of the population in question.  In addition to providing web-accessible analyses and links to third-party analysis and visualisation tools, the BIGSdb software includes a RESTful application programming interface (API) that enables access to all the underlying data for third-party applications and data analysis pipelines.

Genomic surveillance of Neisseria gonorrhoeae to investigate the distribution and evolution of antimicrobial-resistance determinants and lineages

The first extensively drug resistant (XDR) Neisseria gonorrhoeae strain with high resistance to the extended-spectrum cephalosporin ceftriaxone was identified in 2009 in Japan, but no other strain with this antimicrobial-resistance profile has been reported since. However, surveillance to date has been based on phenotypic methods and sequence typing, not genome sequencing. Therefore, little is known about the local population structure at the genomic level, and how resistance determinants and lineages are distributed and evolve. We analysed the whole-genome sequence data and the antimicrobial-susceptibility testing results of 204 strains sampled in a region where the first XDR ceftriaxone-resistant N. gonorrhoeae was isolated, complemented with 67 additional genomes from other time frames and locations within Japan. Strains resistant to ceftriaxone were not found, but we discovered a sequence type (ST)7363 sub-lineage susceptible to ceftriaxone and cefixime in which the mosaic penA allele responsible for reduced susceptibility had reverted to a susceptible allele by recombination. Approximately 85 % of isolates showed resistance to fluoroquinolones (ciprofloxacin) explained by linked amino acid substitutions at positions 91 and 95 of GyrA with 99 % sensitivity and 100 % specificity. Approximately 10 % showed resistance to macrolides (azithromycin), for which genetic determinants are less clear. Furthermore, we revealed different evolutionary paths of the two major lineages: single acquisition of penA X in the ST7363-associated lineage, followed by multiple independent acquisitions of the penA X and XXXIV in the ST1901-associated lineage. Our study provides a detailed picture of the distribution of resistance determinants and disentangles the evolution of the two major lineages spreading worldwide.

Identification of Novel Neisseria gonorrhoeae Lineages Harboring Resistance Plasmids in Coastal Kenya

Background

Africa has the highest incidence of gonorrhea in the world. However, little is known about gonococcal populations in this continent or mechanisms of antimicrobial resistance (AMR).

Methods

Whole-genome sequence data were analyzed from 103 Neisseria gonorrhoeae isolates from 73 patients, mainly men who have sex with men, from coastal Kenya. We annotated loci, defined the core genome, defined mechanisms of AMR, and performed phylogenetic analysis. For patients with multiple episodes of gonorrhea, we determined whether infections occurred with related strains.

Results

We identified 3 clusters of isolates that are phylogenetically distinct from isolates found elsewhere. Plasmids were virtually ubiquitous: pTetM and pblaTEM were found in 97%, and 55% of isolates, respectively. This was associated with high doxycycline use for undiagnosed sexually transmitted infections. Twenty-three percent of multiple episodes of gonorrhea in the same individual were caused by a related strain, suggesting inadequate treatment or reinfection.

Conclusions

The prevalence of plasmid-mediated AMR in Kenyan gonococci contrasts with that in wealthy countries, where AMR is largely chromosomally mediated. Antimicrobials have a profound effect on the maintenance of lineages harboring plasmids. Doxycycline can select for tetracycline and penicillin resistance, through plasmid cooperation. Understanding the mechanisms of AMR in high-risk groups is required to inform treatment strategies.

Transcriptome Analysis of Neisseria gonorrhoeae during Natural Infection Reveals Differential Expression of Antibiotic Resistance Determinants between Men and Women

Recent emergence of antimicrobial resistance of Neisseria gonorrhoeae worldwide has resulted in limited therapeutic choices for treatment of infections caused by this organism. We performed global transcriptomic analysis of N. gonorrhoeae in subjects with gonorrhea who attended a Nanjing, China, sexually transmitted infection (STI) clinic, where antimicrobial resistance of N. gonorrhoeae is high and increasing. We found that N. gonorrhoeae transcriptional responses to infection differed in genital specimens taken from men and women, particularly antibiotic resistance gene expression, which was increased in men. These sex-specific findings may provide a new approach to guide therapeutic interventions and preventive measures that are also sex specific while providing additional insight to address antimicrobial resistance of N. gonorrhoeae.

Neisseria gonorrhoeae is a bacterial pathogen responsible for the sexually transmitted infection gonorrhea. Emergence of antimicrobial resistance (AMR) of N. gonorrhoeae worldwide has resulted in limited therapeutic choices for this infection. Men who seek treatment often have symptomatic urethritis; in contrast, gonococcal cervicitis in women is usually minimally symptomatic, but may progress to pelvic inflammatory disease. Previously, we reported the first analysis of gonococcal transcriptome expression determined in secretions from women with cervical infection. Here, we defined gonococcal global transcriptional responses in urethral specimens from men with symptomatic urethritis and compared these with transcriptional responses in specimens obtained from women with cervical infections and in vitro-grown N. gonorrhoeae isolates. This is the first comprehensive comparison of gonococcal gene expression in infected men and women. RNA sequencing analysis revealed that 9.4% of gonococcal genes showed increased expression exclusively in men and included genes involved in host immune cell interactions, while 4.3% showed increased expression exclusively in women and included phage-associated genes. Infected men and women displayed comparable antibiotic-resistant genotypes and in vitro phenotypes, but a 4-fold higher expression of the Mtr efflux pump-related genes was observed in men. These results suggest that expression of AMR genes is programed genotypically and also driven by sex-specific environments. Collectively, our results indicate that distinct N. gonorrhoeae gene expression signatures are detected during genital infection in men and women. We propose that therapeutic strategies could target sex-specific differences in expression of antibiotic resistance genes.

IMPORTANCE Recent emergence of antimicrobial resistance of Neisseria gonorrhoeae worldwide has resulted in limited therapeutic choices for treatment of infections caused by this organism. We performed global transcriptomic analysis of N. gonorrhoeae in subjects with gonorrhea who attended a Nanjing, China, sexually transmitted infection (STI) clinic, where antimicrobial resistance of N. gonorrhoeae is high and increasing. We found that N. gonorrhoeae transcriptional responses to infection differed in genital specimens taken from men and women, particularly antibiotic resistance gene expression, which was increased in men. These sex-specific findings may provide a new approach to guide therapeutic interventions and preventive measures that are also sex specific while providing additional insight to address antimicrobial resistance of N. gonorrhoeae.

Applications of genomics to slow the spread of multidrug-resistant Neisseria gonorrhoeae

Infections with Neisseria gonorrhoeae, a sexually transmitted pathogen that causes urethritis, cervicitis, and more severe complications, are increasing. Gonorrhea is typically treated with antibiotics; however, N. gonorrhoeae has rapidly acquired resistance to many antibiotic classes, and lineages with reduced susceptibility to the currently recommended therapies are emerging worldwide. In this review, we discuss the contributions of whole genome sequencing (WGS) to our understanding of resistant N. gonorrhoeae. Genomics has illuminated the evolutionary origins and population structure of N. gonorrhoeae and the magnitude of horizontal gene transfer within and between Neisseria species. WGS can be used to predict the susceptibility of N. gonorrhoeae based on known resistance determinants, track the spread of these determinants throughout the N. gonorrhoeae population, and identify novel loci contributing to resistance. WGS has also allowed more detailed epidemiological analysis of transmission of N. gonorrhoeae between individuals and populations than previously used typing methods. Ongoing N. gonorrhoeae genomics will complement other laboratory techniques to understand the biology and evolution of the pathogen, improve diagnostics and treatment in the clinic, and inform public health policies to limit the impact of antibiotic resistance.

Neisseria genomics: current status and future perspectives

High-throughput whole genome sequencing has unlocked a multitude of possibilities enabling members of the Neisseria genus to be examined with unprecedented detail, including the human pathogens Neisseria meningitidis and Neisseria gonorrhoeae. To maximise the potential benefit of this for public health, it is becoming increasingly important to ensure that this plethora of data are adequately stored, disseminated and made readily accessible. Investigations facilitating cross-species comparisons as well as the analysis of global datasets will allow differences among and within species and across geographic locations and different times to be identified, improving our understanding of the distinct phenotypes observed. Recent advances in high-throughput platforms that measure the transcriptome, proteome and/or epigenome are also becoming increasingly employed to explore the complexities of Neisseria biology. An integrated approach to the analysis of these is essential to fully understand the impact these may have in the Neisseria genus. This article reviews the current status of some of the tools available for next generation sequence analysis at the dawn of the ‘post-genomic’ era.

Mobile genetic elements in Neisseria gonorrhoeae: movement for change

Neisseria gonorrhoeae, the causative agent of the sexually transmitted disease gonorrhoeae, possesses several mobile genetic elements (MGEs). The MGEs such as transposable elements mediate intrachromosomal rearrangements, while plasmids and the gonococcal genetic island are involved in interchromosomal gene transfer. Additionally, gonococcal MGEs serve as hotspots for recombination and integration of other genetic elements such as bacteriophages, contribute to gene regulation or spread genes through gonococcal populations by horizontal gene transfer. In this review, we summarise the literature on the structure and biology of MGEs and discuss how these genetic elements may play a role in the pathogenesis and spread of antimicrobial resistance in N. gonorrhoeae. Although an abundance of information about gonococcal MGEs exists (mainly from whole genome sequencing and bioinformatic analysis), there are still many open questions on how MGEs influence the biology of N. gonorrhoeae.

Genomic epidemiology and population structure of Neisseria gonorrhoeae from remote highly endemic Western Australian populations

Background

Neisseria gonorrhoeae causes gonorrhoea, the second most commonly notified sexually transmitted infection in Australia. One of the highest notification rates of gonorrhoea is found in the remote regions of Western Australia (WA). Unlike isolates from the major Australian population centres, the remote community isolates have low rates of antimicrobial resistance (AMR).

Population structure and whole-genome comparison of 59 isolates from the Western Australian N. gonorrhoeae collection were used to investigate relatedness of isolates cultured in the metropolitan and remote areas. Core genome phylogeny, multilocus sequencing typing (MLST), N. gonorrhoeae multi-antigen sequence typing (NG-MAST) and N. gonorrhoeae sequence typing for antimicrobial resistance (NG-STAR) in addition to hierarchical clustering of sequences were used to characterize the isolates.

Results

Population structure analysis of the 59 isolates together with 72 isolates from an international collection, revealed six population groups suggesting that N. gonorrhoeae is a weakly clonal species. Two distinct population groups, Aus1 and Aus2, represented 63% of WA isolates and were mostly composed of the remote community isolates that carried no chromosomal AMR genotypes. In contrast, the Western Australian metropolitan isolates were frequently multi-drug resistant and belonged to population groups found in the international database, suggesting international transmission of the isolates.

Conclusions

Our study suggests that the population structure of N. gonorrhoeae is distinct between the communities in remote and metropolitan WA. Given the high rate of AMR in metropolitan regions, ongoing surveillance is essential to ensure the enduring efficacy of the empiric gonorrhoea treatment in remote WA.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4557-5) contains supplementary material, which is available to authorized users.

Whole-genome sequencing reveals transmission of gonococcal antibiotic resistance among men who have sex with men: an observational study

Objectives

Drug-resistant Neisseria gonorrhoeae are now a global public health threat. Direct transmission of antibiotic-resistant gonococci between individuals has been proposed as a driver for the increased transmission of resistance, but direct evidence of such transmission is limited. Whole-genome sequencing (WGS) has superior resolution to investigate outbreaks and disease transmission compared with traditional molecular typing methods such as multilocus sequence typing (MLST) and N. gonorrhoeae multiantigen sequence (NG-MAST). We therefore aimed to systematically investigate the transmission of N. gonorrhoeae between men in sexual partnerships using WGS to compare isolates and their resistance to antibiotics at a genome level.

Methods

458 couples from a large prospective cohort of men who have sex with men (MSM) tested for gonorrhoea together between 2005 and 2014 were included, and WGS was conducted on all isolates from couples where both men were culture-positive for N. gonorrhoeae. Resistance-determining sequences were identified from genome assemblies, and comparison of isolates between and within individuals was performed by pairwise single nucleotide polymorphism and pangenome comparisons, and in silico predictions of NG-MAST and MLST.

Results

For 33 of 34 (97%; 95% CI 85% to 100%) couples where both partners were positive for gonorrhoea, the resistance-determining genes and mutations were identical in isolates from each partner (94 isolates in total). Resistance determinants in isolates from 23 of 23 (100%; 95% CI 86% to 100%) men with multisite infections were also identical within an individual. These partner and within-host isolates were indistinguishable by NG-MAST, MLST and whole genomic comparisons.

Conclusions

These data support the transmission of antibiotic-resistant strains between sexual partners as a key driver of resistance rates in gonorrhoea among MSM. This improved understanding of the transmission dynamics of N. gonorrhoeae between sexual partners will inform treatment and prevention guidelines.

Clonal Expansion of New Penicillin-Resistant Clade of Neisseria meningitidis Serogroup W Clonal Complex 11, Australia

In Western Australia, Neisseria meningitidis serogroup W clonal complex 11 became the predominant cause of invasive meningococcal disease in 2016. We used core-genome analysis to show emergence of a penicillin-resistant clade that had the penA_253 allele. This new penicillin-resistant clade might affect treatment regimens for this disease.

Genomic analysis of urogenital and rectal Neisseria meningitidis isolates reveals encapsulated hyperinvasive meningococci and coincident multidrug-resistant gonococci

Objective

Invasive meningococcal disease (IMD) outbreaks in men who have sex with men (MSM) have been associated with meningococcal colonisation of the urethra and rectum, but little is known about this colonisation or co-colonisation with the closely related gonococcus. Whole genome sequencing (WGS) was employed to explore these phenomena.

Methods

Meningococci isolated from the urogenital tract and rectum (n=23) and coincident gonococci (n=14) were analysed by WGS along with contemporary meningococci from IMD (n=11). All isolates were obtained from hospital admissions in Brighton, UK, 2011–2013. Assembled WGS were deposited in the PubMLST/neisseria database (http://pubmlst.org/neisseria) and compared at genomic loci common to gonococci or meningococci.

Results

As expected, most meningococci from IMD were encapsulated and belonged to hyperinvasive lineages. So too were meningococci found in the urogenital tract and rectum, contrasting to those asymptomatically carried in the nasopharynx where such meningococci are rare. Five hyperinvasive meningococcal lineages and four distinct gonococcal genotypes were recovered, including multiresistant ST-1901 (NG MAST-1407) gonococci.

Conclusions

These data were consistent with a predisposition for potentially virulent encapsulated hyperinvasive meningococci to colonise the urethra and rectum, which suggests their involvement in MSM IMD outbreaks. The coincidence of multiresistant gonococci raises wider public health concerns.