Neisseria gonorrhoeae Population Genomics: Use of the Gonococcal Core Genome to Improve Surveillance of Antimicrobial Resistance

Development and implementation of a core genome multilocus sequence typing scheme for Haemophilus influenzae

Haemophilus influenzae is part of the human nasopharyngeal microbiota and a pathogen causing invasive disease. The extensive genetic diversity observed in H. influenzae necessitates discriminatory analytical approaches to evaluate its population structure. This study developed a core genome multilocus sequence typing (cgMLST) scheme for H. influenzae using pangenome analysis tools and validated the cgMLST scheme using datasets consisting of complete reference genomes (N = 14) and high-quality draft H. influenzae genomes (N = 2297). The draft genome dataset was divided into a development dataset (N = 921) and a validation dataset (N = 1376). The development dataset was used to identify potential core genes, and the validation dataset was used to refine the final core gene list to ensure the reliability of the proposed cgMLST scheme. Functional classifications were made for all the resulting core genes. Phylogenetic analyses were performed using both allelic profiles and nucleotide sequence alignments of the core genome to test congruence, as assessed by Spearman's correlation and ordinary least square linear regression tests. Preliminary analyses using the development dataset identified 1067 core genes, which were refined to 1037 with the validation dataset. More than 70% of core genes were predicted to encode proteins essential for metabolism or genetic information processing. Phylogenetic and statistical analyses indicated that the core genome allelic profile accurately represented phylogenetic relatedness among the isolates (R 2 = 0.945). We used this cgMLST scheme to define a high-resolution population structure for H. influenzae, which enhances the genomic analysis of this clinically relevant human pathogen.

Genotypic cluster analysis of Neisseria gonorrhoeae reveals a spectrum of sexual mixing including among HIV-serodiscordant men who have sex with men

BACKGROUND

Men who have sex with men (MSM) are more vulnerable to acquiring sexually transmitted infections (STIs). In 2019, for instance, 74% of European Neisseria gonorrhoeae (Ng) cases among males affected MSM. A recent report by the World Health Organization showed that most of the 2020' interim targets to end STIs by 2030 had not been met. A broadened understanding of STI transmission networks could guide future elimination strategies and reduce the STI burden. Therefore, we used whole-genome sequencing (WGS) to determine Ng-clusters and assess sexual mixing.

METHODS

WGS was performed on Ng-isolates collected at the Medical University of Vienna, Austria and was used for core genome multi-locus sequencing typing cluster analysis. Epidemiologic and infection-specific details were extracted from medical records.

RESULTS

Genomic analysis and demographic data were available for 415 isolates, and 43.9% (182/415) were allocated to 31 Ng-clusters. Nine clusters comprised samples from heterosexual individuals only (women N = 4, human immunodeficiency virus (HIV)-negative men N = 49, HIV-positive man N = 1), nine clusters included MSM only (HIV-negative N = 22, HIV-positive N = 13) and 13 clusters included both heterosexuals and MSM (HIV-negative N = 75, HIV-positive N = 18). Current use of HIV pre-exposure prophylaxis (PrEP) was reported by 22.8% of MSM. In multivariate analysis, only 'MSM' predicted clustering with isolates from HIV-positive individuals (adjusted odds ratio 10.24 (95% CI 5.02-20.90)).

CONCLUSIONS

Sexual mixing of HIV-positive, HIV-negative MSM and non-MSM was frequently observed. Furthermore, HIV-serodiscordant clustering highlights the importance of PrEP rollout to avert HIV transmission. Our findings can inform future STI prevention strategies and continuous surveillance efforts are required to keep up with transmission dynamics.

Antimicrobial resistance and heterogeneity of Neisseria gonorrhoeae isolated from patients attending sexually transmitted infection clinics in Lusaka, Zambia

BACKGROUND

Antimicrobial resistance (AMR) of Neisseria gonorrhoeae is a threat to public health as strains have developed resistance to antimicrobials available for the treatment of gonorrhea. Whole genome sequencing (WGS) can detect and predict antimicrobial resistance to enhance the control and prevention of gonorrhea. Data on the molecular epidemiology of N. gonorrhoeae is sparse in Zambia. This study aimed to determine the genetic diversity of N. gonorrhoeae isolated from patients attending sexually transmitted infection (STI) clinics in Lusaka, Zambia.

METHODS

A cross-sectional study that sequenced 38 N. gonorrhoeae isolated from 122 patients with gonorrhea from 2019 to 2020 was conducted. The AMR profiles were determined by the E-test, and the DNA was extracted using the NucliSens easyMaG magnetic device. Whole genome sequencing was performed on the Illumina NextSeq550 platform. The Bacterial analysis pipeline (BAP) that is readily available at: https://cge.cbs.dtu.dk/services/CGEpipeline-1.1 was used for the identification of the species, assembling the genome, multi-locus sequence typing (MLST), detection of plasmids and AMR genes. Phylogeny by single nucleotide polymorphisms (SNPs) was determined with the CCphylo dataset.

RESULTS

The most frequent STs with 18.4% of isolates each were ST7363, ST1921 and ST1582, followed by ST1583 (13%), novel ST17026 (7.9%), ST1588 (7.9%), ST1596 (5.3%), ST11181 (5.3%), ST11750 (2.6/%) and ST11241 (2.6%) among the 38 genotyped isolates. The blaTeM-1B and tetM (55%) was the most prevalent combination of AMR genes, followed by blaTeM-1B (18.4%), tetM (15.8%), and the combination of blaTeM-1B, ermT, and tetL was 2.6% of the isolates. The AMR phenotypes were predicted in ciprofloxacin, penicillin, tetracycline, azithromycin, and cefixime. The combination of mutations 23.7% was gryA (S91F), parC (E91G), ponA (L421) and rpsJ (V57M), followed by 18.4% in gyrA (S91F), ponA (L421P), rpsJ (V57M), and 18.4% in gyrA (D95G, S91F), ponA (L421P), and rpsJ (V57M). The combinations in gyrA (D95G, S91F) and rpsJ (V57M), and gyrA (D95G, S91F), parC (E91F), ponA (L421P) and rpsJ (V57M) were 13.2% each of the isolates. Plasmid TEM-1 (84.2%), tetM (15.8%), and gonococcal genetic island (GGI) was detected in all isolates.

CONCLUSION

This study revealed remarkable heterogeneity of N. gonorrhoeae with blaTEM-1, tetM, ponA, gyrA, and parC genes associated with high resistance to penicillin, tetracycline, and ciprofloxacin demanding revision of the standard treatment guidelines and improved antimicrobial stewardship in Zambia.

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.

Characterizing the diversity and commensal origins of penA mosaicism in the genus Neisseria

Mosaic penA alleles formed through horizontal gene transfer (HGT) have been instrumental to the rising incidence of ceftriaxone-resistant gonococcal infections. Although interspecies HGT of regions of the penA gene between Neisseria gonorrhoeae and commensal Neisseria species has been described, knowledge concerning which species are the most common contributors to mosaic penA alleles is limited, with most studies examining only a small number of alleles. Here, we investigated the origins of recombinant penA alleles through in silico analyses that incorporated 1700 penA alleles from 35 513 Neisseria isolates, comprising 15 different Neisseria species. We identified Neisseria subflava and Neisseria cinerea as the most common source of recombinant sequences in N. gonorrhoeae penA. This contrasted with Neisseria meningitidis penA, for which the primary source of recombinant DNA was other meningococci, followed by Neisseria lactamica. Additionally, we described the distribution of polymorphisms implicated in antimicrobial resistance in penA, and found that these are present across the genus. These results provide insight into resistance-related changes in the penA gene across human-associated Neisseria species, illustrating the importance of genomic surveillance of not only the pathogenic Neisseria, but also of the oral niche-associated commensals from which these pathogens are sourcing key genetic variation.

Distributed genotyping and clustering of Neisseria strains reveal continual emergence of epidemic meningococcus over a century

Core genome multilocus sequence typing (cgMLST) is commonly used to classify bacterial strains into different types, for taxonomical and epidemiological applications. However, cgMLST schemes require central databases for the nomenclature of new alleles and sequence types, which must be synchronized worldwide and involve increasingly intensive calculation and storage demands. Here, we describe a distributed cgMLST (dcgMLST) scheme that does not require a central database of allelic sequences and apply it to study evolutionary patterns of epidemic and endemic strains of the genus Neisseria. We classify 69,994 worldwide Neisseria strains into multi-level clusters that assign species, lineages, and local disease outbreaks. We divide Neisseria meningitidis into 168 endemic lineages and three epidemic lineages responsible for at least 9 epidemics in the past century. According to our analyses, the epidemic and endemic lineages experienced very different population dynamics in the past 100 years. Epidemic lineages repetitively emerged from endemic lineages, disseminated worldwide, and apparently disappeared rapidly afterward. We propose a stepwise model for the evolutionary trajectory of epidemic lineages in Neisseria, and expect that the development of similar dcgMLST schemes will facilitate epidemiological studies of other bacterial pathogens.

Sustained Transmission of Neisseria gonorrhoeae Strains with High-Level Azithromycin Resistance (MIC ≥ 256 μg/mL) in Argentina, 2018 to 2022

Azithromycin combined with ceftriaxone is the recommended dual therapy for uncomplicated gonorrhea in many countries. Nevertheless, the increasing prevalence of azithromycin resistance compromises the effectiveness of this treatment strategy. From 2018 to 2022, we collected 13 gonococcal isolates with high-level azithromycin resistance (MIC ≥ 256 μg/mL) across Argentina. Whole-genome sequencing revealed that these isolates were mainly represented by the internationally spreading Neisseria gonorrhoeae multi-antigen sequence typing (NG-MAST) genogroup G12302, containing the 23S rRNA A2059G mutation (in all four alleles) together with mosaic mtrD and mtrR promoter 2 loci. This information is important to develop targeted public health policies to control the spread of azithromycin-resistant N. gonorrhoeae in Argentina and internationally. IMPORTANCE Azithromycin resistance in Neisseria gonorrhoeae has been increasing in numerous populations worldwide, which is of concern, as azithromycin is part of the recommended dual treatment in many countries. Here, we report 13 N. gonorrhoeae isolates with high-level azithromycin resistance (MIC ≥ 256 μg/mL). This study observed that high-level azithromycin-resistant gonococcal strains have shown sustained transmission in Argentina and are related to the successful international clone NG-MAST G12302. Genomic surveillance together with real-time tracing and data-sharing networks will be crucial in controlling the spread of azithromycin resistance in gonococcus.

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.

Biochemical Atypia in Russian Neisseria gonorrhoeae Clinical Isolates Belonging to the G807 NG-MAST Genogroup/ST1594 MLST

Many current gonococcal clinical isolates in Russia show atypical taxonomically significant biochemical activity, which leads to species misidentification. Molecular typing of such cultures according Neisseria gonorrhoeae multiantigen sequence typing (NG-MAST) and multilocus sequence typing (MLST) protocols assigned them to the G807 NG-MAST GENOGROUP/ST1594 MLST that has been predominant in Russia in recent years. The goal of the study was to analyze the molecular mechanisms of biochemical atypia in N. gonorrhoeae clinical isolates characterized as the members of G807 NG-MAST GENOGROUP/ST1594 MLST. Sixteen isolates of this genogroup were included in the study, eight showed defective amino acid metabolism or loss of D-glucose fermentation. Comparative bioinformatic analysis based on WGS data divided these isolates into two clusters strictly associated with typical or atypical biochemical activity. Cultures with defective amino acid metabolism had a 5-nucleotide insertion in the pip-gene that caused a stop codon and led to synthesis of the non-functional enzyme. Comparison of the sequenced genomes with publicly available N. gonorrhoeae genomes showed the rarity of this insertion. In the global N. gonorrhoeae phylogenetic tree the G807 NG-MAST GENOGROUP/ST1594 MLST forms a distinct branch characterized by 170 SNPs, most of which are non-synonymous. We hypothesized a unique strategy for G807 NG-MAST GENOGROUP/ST1594 MLST clone persistence in the global N. gonorrhoeae population via escape of antimicrobial therapy due to diagnostic misidentification.

Molecular characterisation of Neisseria gonorrhoeae associated with disseminated gonococcal infections in Queensland, Australia: a retrospective surveillance study

OBJECTIVES

Gonorrhoea caused by Neisseria gonorrhoeae is the second most notified sexually transmitted infection (STI) in Australia and the case numbers for this STI have been increasing globally. Progressive gonococcal infection may lead to disseminated gonococcal infection (DGI), which causes significant morbidity among patients. This study aims to examine the genetic diversity of N. gonorrhoeae isolates collected in Queensland from January 2010 to August 2015 and to determine factors associated with DGI in Queensland.

DESIGN

Retrospective surveillance study for epidemiological purposes.

SETTING

All gonorrhoeae isolates referred by private and public pathology laboratories to the state of Queensland, Australia Neisseria reference laboratory.

METHODS

Between January 2010 and August 2015, 3953 N. gonorrhoeae isolates from both metropolitan and regional Queensland infections were typed with NG-MAST (N. gonorrhoeae multiantigen sequence typing) to assess the genetic diversity between strains. Whole-genome sequencing (WGS) was used to investigate strain-related factors associated with DGI.

RESULTS

ST6876 was the most common NG-MAST type, detected in 7.6% of the isolates. DGI was significantly more likely in females <30 years (OR 13.02, p<0.0001) and in older males >30 years (OR 6.04, p<0.0001), with most cases originating from North Queensland (OR 8.5, p<0.0001). Strains harbouring PIA class of porB type were associated with DGI (OR 33.23, p<0.0001).

CONCLUSION

Genotyping techniques, such as NG-MAST and WGS, are proving instrumental in providing an insight into the population structure of N. gonorrhoeae, and genetic mechanisms of pathogenesis, such as for DGI.

Association of Phylogenomic Relatedness among Neisseria gonorrhoeae Strains with Antimicrobial Resistance, Austria, 2016-2020

We investigated genomic determinants of antimicrobial resistance in 1,318 Neisseria gonorrhoeae strains isolated in Austria during 2016–2020. Sequence type (ST) 9363 and ST11422 isolates had high rates of azithromycin resistance, and ST7363 isolates correlated with cephalosporin resistance. These results underline the benefit of genomic surveillance for antimicrobial resistance monitoring.

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.

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.

Global emergence and dissemination of Neisseria gonorrhoeae ST-9363 isolates with reduced susceptibility to azithromycin

Neisseria gonorrhoeae multi-locus sequence type (ST) 9363 core-genogroup isolates have been associated with reduced azithromycin susceptibility (AZMrs) and show evidence of clonal expansion in the U.S. Here we analyze a global collection of ST-9363 core-genogroup genomes to shed light on the emergence and dissemination of this strain. The global population structure of ST-9363 core-genogroup falls into three lineages: Basal, European, and North American; with 32 clades within all lineages. Although, ST-9363 core-genogroup is inferred to have originated from Asia in the mid-19th century; we estimate the three modern lineages emerged from Europe in the late 1970s to early 1980s. The European lineage appears to have emerged and expanded from around 1986 to 1998, spreading into North America and Oceania in the mid-2000s with multiple introductions, along with multiple secondary reintroductions into Europe. Our results suggest two separate acquisition events of mosaic mtrR and mtrR promoter alleles: first during 2009-2011 and again during the 2012-2013 time, facilitating the clonal expansion of this core-genogroup with AZMrs in the U.S. By tracking phylodynamic evolutionary trajectories of clades that share distinct demography as well as population-based genomic statistics, we demonstrate how recombination and selective pressures in the mtrCDE efflux operon granted a fitness advantage to establish ST-9363 as a successful gonococcal lineage in the U.S. and elsewhere. Although it is difficult to pinpoint the exact timing and emergence of this young core-genogroup, it remains critically important to continue monitoring it, as it could acquire additional resistance markers.

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.

Phylogenomic analysis reveals persistence of gonococcal strains with reduced-susceptibility to extended-spectrum cephalosporins and mosaic penA-34

The recent emergence of strains of Neisseria gonorrhoeae associated with treatment failures to ceftriaxone, the foundation of current treatment options, has raised concerns over a future of untreatable gonorrhea. Current global data on gonococcal strains suggest that several lineages, predominately characterized by mosaic penA alleles, are associated with elevated minimum inhibitory concentrations (MICs) to extended spectrum cephalosporins (ESCs). Here we report on whole genome sequences of 813 N. gonorrhoeae isolates collected through the Gonococcal Isolate Surveillance Project in the United States. Phylogenomic analysis revealed that one persisting lineage (Clade A, multi-locus sequence type [MLST] ST1901) with mosaic penA-34 alleles, contained the majority of isolates with elevated MICs to ESCs. We provide evidence that an ancestor to the globally circulating MLST ST1901 clones potentially emerged around the early to mid-20th century (1944, credibility intervals [CI]: 1935-1953), predating the introduction of cephalosporins, but coinciding with the use of penicillin. Such results indicate that drugs with novel mechanisms of action are needed as these strains continue to persist and disseminate globally.

Using Neisseria meningitidis genomic diversity to inform outbreak strain identification

Meningococcal disease is a life-threatening illness caused by the human-restricted bacterium Neisseria meningitidis. Outbreaks in the USA involve at least two cases in an organization or community caused by the same serogroup within three months. Genome comparisons, including phylogenetic analysis and quantification of genome distances can provide confirmatory evidence of pathogen transmission during an outbreak. Interpreting genome distances depends on understanding their distribution both among isolates from outbreaks and among those not from outbreaks. Here, we identify outbreak strains based on phylogenetic relationships among 141 N. meningitidis isolates collected from 28 outbreaks in the USA during 2010-2017 and 1516 non-outbreak isolates collected through contemporaneous meningococcal surveillance. We show that genome distance thresholds based on the maximum SNPs and allele distances among isolates in the phylogenetically defined outbreak strains are sufficient to separate most pairs of non-outbreak isolates into separate strains. Non-outbreak isolate pairs that could not be distinguished from each other based on genetic distances were concentrated in the clonal complexes CC11, CC103, and CC32. Within each of these clonal complexes, phylodynamic analysis identified a group of isolates with extremely low diversity, collected over several years and multiple states. Clusters of isolates with low genetic diversity could indicate increased pathogen transmission, potentially resulting in local outbreaks or nationwide clonal expansions.

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.

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.

Neisseria gonorrhoeae clustering to reveal major European whole-genome-sequencing-based genogroups in association with antimicrobial resistance

Neisseria gonorrhoeae, the bacterium responsible for the sexually transmitted disease gonorrhoea, has shown an extraordinary ability to develop antimicrobial resistance (AMR) to multiple classes of antimicrobials. With no available vaccine, managing N. gonorrhoeae infections demands effective preventive measures, antibiotic treatment and epidemiological surveillance. The latter two are progressively being supported by the generation of whole-genome sequencing (WGS) data on behalf of national and international surveillance programmes. In this context, this study aims to perform N. gonorrhoeae clustering into genogroups based on WGS data, for enhanced prospective laboratory surveillance. Particularly, it aims to identify the major circulating WGS-genogroups in Europe and to establish a relationship between these and AMR. Ultimately, it enriches public databases by contributing with WGS data from Portuguese isolates spanning 15 years of surveillance. A total of 3791 carefully inspected N. gonorrhoeae genomes from isolates collected across Europe were analysed using a gene-by-gene approach (i.e. using cgMLST). Analysis of cluster composition and stability allowed the classification of isolates into a two-step hierarchical genogroup level determined by two allelic distance thresholds revealing cluster stability. Genogroup clustering in general agreed with available N. gonorrhoeae typing methods [i.e. MLST (multilocus sequence typing), NG-MAST (N. gonorrhoeae multi-antigen sequence typing) and PubMLST core-genome groups], highlighting the predominant genogroups circulating in Europe, and revealed that the vast majority of the genogroups present a dominant AMR profile. Additionally, a non-static gene-by-gene approach combined with a more discriminatory threshold for potential epidemiological linkage enabled us to match data with previous reports on outbreaks or transmission chains. In conclusion, this genogroup assignment allows a comprehensive analysis of N. gonorrhoeae genetic diversity and the identification of the WGS-based genogroups circulating in Europe, while facilitating the assessment (and continuous monitoring) of their frequency, geographical dispersion and potential association with specific AMR signatures. This strategy may benefit public-health actions through the prioritization of genogroups to be controlled, the identification of emerging resistance carriage, and the potential facilitation of data sharing and communication.

Sudden emergence of a Neisseria gonorrhoeae clade with reduced susceptibility to extended-spectrum cephalosporins, Norway

Neisseria gonorrhoeae multilocus sequence type (ST)-7827 emerged in a dramatic fashion in Norway in the period 2016-2018. Here, we aim to shed light on the provenance and expansion of this ST. ST-7827 was found to be polyphyletic, but the majority of members belonged to a monophyletic clade we termed PopPUNK cluster 7827 (PC-7827). In Norway, both PC-7827 and ST-7827 isolates were almost exclusively isolated from men. Phylogeographical analyses demonstrated an Asian origin of the genogroup, with multiple inferred exports to Europe and the USA. The genogroup was uniformly resistant to fluoroquinolones, and associated with reduced susceptibility to both azithromycin and the extended-spectrum cephalosporins (ESCs) cefixime and ceftriaxone. From a genetic background including the penA allele 13.001, associated with reduced ESC susceptibility, we identified repeated events of acquisition of porB alleles associated with further reduction in ceftriaxone susceptibility. Transmission of the strain was significantly reduced in Norway in 2019, but our results indicate the existence of a recently established global reservoir. The worrisome drug-resistance profile and rapid emergence of PC-7827 calls for close monitoring of the situation.

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

Plasmids are vehicles for horizontal gene transfer between bacteria, and in Neisseria gonorrhoeae plasmids can mediate high-level antimicrobial resistance (AMR). Using genomic and phylogenetic analyses, we show that plasmids are widespread in a collection of 3724 gonococcal isolates from 56 countries, and characterized the conjugative, β-lactamase and cryptic plasmids. We found that variants of the conjugative plasmid (which can mediate tetracycline resistance) and the β-lactamase plasmid expressing TEM-135 are associated with distinct gonococcal lineages. Furthermore, AMR plasmids are significantly more prevalent in gonococci from less wealthy countries, highlighting the need for further studies. More than 94% of gonococci possess the cryptic plasmid, with its absence correlated with the presence of a novel chromosomal type IV secretion system. Our results reveal the extent of plasmid-mediated AMR in the gonococcus, particularly in less wealthy countries, where diagnostic and therapeutic options can be limited, and highlight the risk of their global spread.