Using crude whole-genome assemblies of Neisseria gonorrhoeae as a platform for strain analysis: clonal spread of gonorrhea infection in Saskatchewan, Canada

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.

Genomic Analysis Reveals Antibiotic-Susceptible Clones and Emerging Resistance in Neisseria gonorrhoeae in Saskatchewan, Canada

Whole-genome sequencing was used to identify mutations in antibiotic resistance-conferring genes to compare susceptibility predictions with MICs and to ascertain strain types in 99 isolates of Neisseria gonorrhoeae Genotypes associated with susceptibility, as well as MIC creep or emerging resistance, were noted. Phylogenomic analysis revealed three distinctive clades and putative gonococcal transmission linkages involving a tetracycline-resistant N. gonorrhoeae outbreak and the clonal spread of susceptible isolates in men.

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.

Gen2Epi: an automated whole-genome sequencing pipeline for linking full genomes to antimicrobial susceptibility and molecular epidemiological data in Neisseria gonorrhoeae

BACKGROUND

Recent adva1nces in whole genome sequencing (WGS) based technologies have facilitated multi-step applications for predicting antimicrobial resistance (AMR) and investigating the molecular epidemiology of Neisseria gonorrhoeae. However, generating full scaffolds of N. gonorrhoeae genomes from short reads, and the assignment of molecular epidemiological information (NG-MLST, NG-MAST, and NG-STAR) to multiple assembled samples, is challenging due to required manual tasks such as annotating antimicrobial resistance determinants with standard nomenclature for a large number of genomes.

RESULTS

We present Gen2Epi, a pipeline that assembles short reads into full scaffolds and automatically assigns molecular epidemiological and AMR information to the assembled genomes. Gen2Epi is a command-line tool integrating third-party software and tailored specifically for N. gonorrhoeae. For its evaluation, the Gen2Epi pipeline successfully assembled the WGS short reads from 1484 N. gonorrhoeae samples into full-length genomes for both chromosomes and plasmids and was able to assign in silico molecular determinant information to each dataset automatically. The assemblies were generated using raw as well as trimmed short reads. The median genome coverage of full-length scaffolds and "N" statistics (N50, NG50, and NGA50) were higher than, or comparable to, previously published results and the scaffolding process improved the quality of the draft genome assemblies. Molecular antimicrobial resistant (AMR) determinants identified by Gen2Epi reproduced information for the 1484 samples as previously reported, including NG-MLST, NG-MAST, and NG-STAR molecular sequence types.

CONCLUSIONS

Gen2Epi can be used to assemble short reads into full-length genomes and assign accurate molecular marker and AMR information automatically from NG-STAR, NG-MAST, and NG-MLST. Gen2Epi is publicly available under "CC BY-NC 2.0 CA" Creative Commons licensing as a VirtualBox image containing the constituent software components running on the LINUX operating system (CentOS 7). The image and associated documentation are available via anonymous FTP at ftp://www.cs.usask.ca/pub/combi or ftp://ftp.cs.usask.ca/pub/combi.

Evaluation of a Hydrogel-Based Diagnostic Approach for the Point-of-Care Based Detection of Neisseria gonorrhoeae

Eleven primer pairs were developed for the identification of Neisseria gonorrhoeae. The sensitivity and specificity of these primers were evaluated by Real Time (RT)-PCR melt curve analyses with DNA from 145 N. gonorrhoeae isolates and 40 other Neisseria or non-Neisseria species. Three primer pairs were further evaluated in a hydrogel-based RT-PCR detection platform, using DNA extracted from 50 N. gonorrhoeae cultures. We observed 100% sensitivity and specificity in the hydrogel assay, confirming its potential as a point-of-care test (POCT) for N. gonorrhoeae diagnosis.

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.

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.

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.

Multiplex Real-Time PCR Assay for Simultaneous Identification of Neisseria gonorrhoeae and Its Ciprofloxacin Susceptibility Status

A real-time PCR (RT-PCR) assay was designed for the simultaneous identification of Neisseria gonorrhoeae and its ciprofloxacin susceptibility status. A SYBR green-based multiplex RT-PCR format was used; it comprised two different forward primers and a common reverse primer to detect single nucleotide polymorphisms (SNPs) in gyrA of N. gonorrhoeae The primer pairs were evaluated for their sensitivity and specificity using genomic DNA from 254 N. gonorrhoeae isolates (82 were ciprofloxacin susceptible and 172 were ciprofloxacin resistant) and 23 non-N. gonorrhoeae species isolates. The performance of the primers was validated using genomic DNA from 100 different N. gonorrhoeae isolates (46 were ciprofloxacin susceptible and 54 were ciprofloxacin resistant) and 52 non-N. gonorrhoeae isolates. The latter panel was revalidated by testing 99 (46 isolates were ciprofloxacin susceptible and 53 isolates were ciprofloxacin resistant) of the N. gonorrhoeae isolates and 23 non-N. gonorrhoeae isolates. These primers detected N. gonorrhoeae and its ciprofloxacin susceptibility status with over 99% sensitivity and specificity for all panels tested. This assay has the potential to be an inexpensive and rapid test for the simultaneous identification of N. gonorrhoeae and its ciprofloxacin susceptibility status.

Antimicrobial resistance in Neisseria gonorrhoeae: history, molecular mechanisms and epidemiological aspects of an emerging global threat

Neisseria gonorrhoeae is the agent of gonorrhea, a sexually transmitted infection with an estimate from The World Health Organization of 78 million new cases in people aged 15-49 worldwide during 2012. If left untreated, complications may include pelvic inflammatory disease and infertility. Antimicrobial treatment is usually effective; however, resistance has emerged successively through various molecular mechanisms for all the regularly used therapeutic agents throughout decades. Detection of antimicrobial susceptibility is currently the most critical aspect for N. gonorrhoeae surveillance, however poorly structured health systems pose difficulties. In this review, we compiled data from worldwide reports regarding epidemiology and antimicrobial resistance in N. gonorrhoeae, and highlight the relevance of the implementation of surveillance networks to establish policies for gonorrhea treatment.

Multidrug-resistant gonorrhea: A research and development roadmap to discover new medicines

Emilie Alirol and colleagues discuss the development of new treatments for gonorrhea.

Antimicrobial resistance in Neisseria gonorrhoeae: Global surveillance and a call for international collaborative action

In a Policy Forum, Teodora Wi and colleagues discuss the challenges of antimicrobial resistance in gonococci.

Genomic sequencing of Neisseria gonorrhoeae to respond to the urgent threat of antimicrobial-resistant gonorrhea

The development of resistance of Neisseria gonorrhoeae to available first-line antibiotics, including penicillins, tetracyclines, fluoroquinolones and cephalosporins, has led to the circulation of multidrug-resistant gonorrhea at a global scale. Advancements in high-throughput whole-genome sequencing (WGS) provide useful tools that can be used to enhance gonococcal detection, treatment and management capabilities, which will ultimately aid in the control of antimicrobial resistant gonorrhea worldwide. In this minireview, we discuss the application of WGS of N. gonorrhoeae to strain typing, phylogenomic, molecular surveillance and transmission studies. We also examine the application of WGS analyses to the public health sector as well as the potential usage of WGS-based transcriptomic and epigenetic methods to identify novel gonococcal resistance mechanisms.

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

OBJECTIVES

Antimicrobial resistance (AMR) threatens our ability to treat the sexually transmitted bacterial infection gonorrhoea. The increasing availability of whole genome sequence (WGS) data from Neisseria gonorrhoeae isolates, however, provides us with an opportunity in which WGS can be mined for AMR determinants.

METHODS

Chromosomal and plasmid genes implicated in AMR were catalogued on the PubMLST Neisseria database (http://pubmlst.org/neisseria). AMR genotypes were identified in WGS from 289 gonococci for which MICs against several antimicrobial compounds had been determined. Whole genome comparisons were undertaken using whole genome MLST (wgMLST).

RESULTS

Clusters of isolates with distinct AMR genotypes were apparent following wgMLST analysis consistent with the occurrence of genome wide genetic variation. This included the presence of the gonococcal genetic island (GGI), a type 4 secretion system shown to increase recombination and for which possession was significantly associated with AMR to multiple antimicrobials.

CONCLUSIONS

Evolution of the gonococcal genome occurs in response to antimicrobial selective pressure resulting in the formation of distinct N. gonorrhoeae populations evidenced by the wgMLST clusters seen here. Genomic islands offer selective advantages to host bacteria and possession of the GGI may, not only facilitate the spread of AMR in gonococcal populations, but may also confer fitness advantages.

Whole-genome sequencing targets drug-resistant bacterial infections

During the past two decades, the technological progress of whole-genome sequencing (WGS) had changed the fields of Environmental Microbiology and Biotechnology, and, currently, is changing the underlying principles, approaches, and fundamentals of Public Health, Epidemiology, Health Economics, and national productivity. Today’s WGS technologies are able to compete with conventional techniques in cost, speed, accuracy, and resolution for day-to-day control of infectious diseases and outbreaks in clinical laboratories and in long-term epidemiological investigations. WGS gives rise to an exciting future direction for personalized Genomic Epidemiology. One of the most vital and growing public health problems is the emerging and re-emerging of multidrug-resistant (MDR) bacterial infections in the communities and healthcare settings, reinforced by a decline in antimicrobial drug discovery. In recent years, retrospective analysis provided by WGS has had a great impact on the identification and tracking of MDR microorganisms in hospitals and communities. The obtained genomic data are also important for developing novel easy-to-use diagnostic assays for clinics, as well as for antibiotic and therapeutic development at both the personal and population levels. At present, this technology has been successfully applied as an addendum to the real-time diagnostic methods currently used in clinical laboratories. However, the significance of WGS for public health may increase if: (a) unified and user-friendly bioinformatics toolsets for easy data interpretation and management are established, and (b) standards for data validation and verification are developed. Herein, we review the current and future impact of this technology on diagnosis, prevention, treatment, and control of MDR infectious bacteria in clinics and on the global scale.

Population structure of Neisseria gonorrhoeae based on whole genome data and its relationship with antibiotic resistance

Neisseria gonorrhoeae is the causative agent of gonorrhea, a sexually transmitted infection (STI) of major importance. As a result of antibiotic resistance, there are now limited options for treating patients. We collected draft genome sequence data and associated metadata data on 76 N. gonorrhoeae strains from around the globe and searched for known determinants of antibiotics resistance within the strains. The population structure and evolutionary forces within the pathogen population were analyzed. Our results indicated a cosmopolitan gonoccocal population mainly made up of five subgroups. The estimated ratio of recombination to mutation (r/m = 2.2) from our data set indicates an appreciable level of recombination occurring in the population. Strains with resistance phenotypes to more recent antibiotics (azithromycin and cefixime) were mostly found in two of the five population subgroups.

Whole-genome phylogenomic heterogeneity of Neisseria gonorrhoeae isolates with decreased cephalosporin susceptibility collected in Canada between 1989 and 2013

A large-scale, whole-genome comparison of Canadian Neisseria gonorrhoeae isolates with high-level cephalosporin MICs was used to demonstrate a genomic epidemiology approach to investigate strain relatedness and dynamics. Although current typing methods have been very successful in tracing short-chain transmission of gonorrheal disease, investigating the temporal evolutionary relationships and geographical dissemination of highly clonal lineages requires enhanced resolution only available through whole-genome sequencing (WGS). Phylogenomic cluster analysis grouped 169 Canadian strains into 12 distinct clades. While some N. gonorrhoeae multiantigen sequence types (NG-MAST) agreed with specific phylogenomic clades or subclades, other sequence types (ST) and closely related groups of ST were widely distributed among clades. Decreased susceptibility to extended-spectrum cephalosporins (ESC-DS) emerged among a group of diverse strains in Canada during the 1990s with a variety of nonmosaic penA alleles, followed in 2000/2001 with the penA mosaic X allele and then in 2007 with ST1407 strains with the penA mosaic XXXIV allele. Five genetically distinct ESC-DS lineages were associated with penA mosaic X, XXXV, and XXXIV alleles and nonmosaic XII and XIII alleles. ESC-DS with coresistance to azithromycin was observed in 5 strains with 23S rRNA C2599T or A2143G mutations. As the costs associated with WGS decline and analysis tools are streamlined, WGS can provide a more thorough understanding of strain dynamics, facilitate epidemiological studies to better resolve social networks, and improve surveillance to optimize treatment for gonorrheal infections.