Neisseria gonorrhoeae is a global threat to public health due to the accumulation of antimicrobial resistance mechanisms. ST-1901 is an internationally important sequence type (ST) because of its high incidence and the usual occurrence of chromosomally determined resistance. In this study, we describe the evolution of the ST-1901 and its single locus variants in Rio de Janeiro from 2006 to 2022. We analyzed 82 N. gonorrhoeae isolates according to antimicrobial susceptibility profile, resistance mechanisms, molecular typing, and phylogenetics. Six different single locus variants were detected. Phylogenetic analysis identified five clades, which share similar characteristics. Resistance rates for penicillin and tetracycline decreased due to the lower occurrence of resistance plasmids, but intermediary resistance to penicillin rose. Resistance to ciprofloxacin remained high throughout all clades and the years of the study. Regarding resistance to azithromycin, alterations in mtrR promoter and gene, and 23S rRNA encoding gene rrl were detected, with a notable rise in the incidence of C2611T mutations in more recent years occurring in 4 out of 5 clades. In contrast, beta-lactam resistance associated penA 34 mosaic was found only in one persisting clade (Clade D), as well as unique G45D and A39T mutations in mtrR gene and its promoter (Nm-Like) were found in only Clade B. Taken together, these data suggest that ST-1901, a persistently circulating lineage of N. gonorrhoeae in Rio de Janeiro, has undergone changes over the years and may evolve to develop resistance to the current recommended dual therapy adopted in Brazil, ceftriaxone and azithromycin.

Neisseria gonorrhoeae multi-locus sequence type (ST) 9363 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 genogroup genomes to shed light on the emergence and dissemination of this strain. The global population structure of ST-9363 genogroup falls into three lineages: Basal, European, and North American; with 32 clades within all lineages. Although, ST-9363 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 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 genogroup, it remains critically important to continue monitoring it, as it could acquire additional resistance markers.Competing Interest StatementThe authors have declared no competing interest.

Recent reports suggest that mosaic-like sequences within the mtr (multiple transferable resistance) efflux pump locus of Neisseria gonorrhoeae likely originating from commensal Neisseria sp. by transformation can increase the ability of gonococci to resist structurally diverse antimicrobials. Thus, acquisition of numerous nucleotide changes within the mtrR gene encoding the transcriptional repressor (MtrR) of the mtrCDE efflux pump-encoding operon or overlapping promoter region for both along with those that cause amino acid changes in the MtrD transporter protein were recently reported to decrease gonococcal susceptibility to numerous antimicrobials, including azithromycin (Azi) (Wadsworth et al. 2018. MBio. doi.org/10.1128/mBio.01419-18). We performed detailed genetic and molecular studies to define the mechanistic basis for why such strains can exhibit decreased susceptibility to MtrCDE antimicrobial substrates including Azi. We report that a strong cis-acting transcriptional impact of a single nucleotide change within the -35 hexamer of the mtrCDE promoter as well gain-of-function amino acid changes at the C-terminal region of MtrD can mechanistically account for the decreased antimicrobial susceptibility of gonococci with a mosaic-like mtr locus.IMPORTANCE Historically, after introduction of an antibiotic for treatment of gonorrhea, strains of N. gonorrhoeae emerge that display clinical resistance due to spontaneous mutation or acquisition of resistance genes. Genetic exchange between members of the Neisseria genus occurring by transformation can cause significant changes in gonococci that impact the structure of an antibiotic target or expression of genes involved in resistance. The results presented herein provide a framework for understanding how mosaic-like DNA sequences from commensal Neisseria that recombine within the gonococcal mtr efflux pump locus function to decrease bacterial susceptibility to antimicrobials including antibiotics used in therapy of gonorrhea.

This study characterized high-quality whole-genome sequences of a sentinel, surveillance-based collection of 1710 Neisseria gonorrhoeae (GC) isolates from 2019 collected in the USA as part of the Gonococcal Isolate Surveillance Project (GISP). It aims to provide a detailed report of strain diversity, phylogenetic relationships and resistance determinant profiles associated with reduced susceptibilities to antibiotics of concern. The 1710 isolates represented 164 multilocus sequence types and 21 predominant phylogenetic clades. Common genomic determinants defined most strains' phenotypic, reduced susceptibility to current and historic antibiotics (e.g. bla (TEM) plasmid for penicillin, tetM plasmid for tetracycline, gyrA for ciprofloxacin, 23S rRNA and/or mosaic mtr operon for azithromycin, and mosaic penA for cefixime and ceftriaxone). The most predominant phylogenetic clade accounted for 21 % of the isolates, included a majority of the isolates with low-level elevated MICs to azithromycin (2.0 µg ml(-1)), carried a mosaic mtr operon and variants in PorB, and showed expansion with respect to data previously reported from 2018. The second largest clade predominantly carried the GyrA S91F variant, was largely ciprofloxacin resistant (MIC ≥1.0 µg ml(-1)), and showed significant expansion with respect to 2018. Overall, a low proportion of isolates had medium- to high-level elevated MIC to azithromycin ((≥4.0 µg ml(-1)), based on C2611T or A2059G 23S rRNA variants). One isolate carried the penA 60.001 allele resulting in elevated MICs to cefixime and ceftriaxone of 1.0 µg ml(-1). This high-resolution snapshot of genetic profiles of 1710 GC sequences, through a comparison with 2018 data (1479 GC sequences) within the sentinel system, highlights change in proportions and expansion of select GC strains and the associated genetic mechanisms of resistance. The knowledge gained through molecular surveillance may support rapid identification of outbreaks of concern. Continued monitoring may inform public health responses to limit the development and spread of antibiotic-resistant gonorrhoea.

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.

BACKGROUND: Sexual networks are difficult to construct due to incomplete sexual partner data. The proximity of people within a network may be inferred from genetically similar infections. We explored genomic data combined with partner services investigation (PSI) data to extend our understanding of sexual networks affected by Neisseria gonorrhoeae (NG). METHODS: We used 2017-2019 PSI and whole-genome sequencing (WGS) data from eight jurisdictions participating in CDC's Strengthening the United States Response to Resistant Gonorrhea (SURRG) project. Clusters were identified from sexual contacts and through genetically similar NG isolates. Sexual mixing patterns were characterized by describing the clusters by the individual's gender and gender of their sex partners. RESULTS: Our study included 4,627 diagnoses of NG infection (81% sequenced), 2,455 people received a PSI, 393 people were negative contacts of cases, and 495 contacts with unknown NG status. We identified 823 distinct clusters using PSI data combined with WGS data. Of cases that were not linked to any other case using PSI data, 37% were linked when using WGS data. Overall, 40% of PSI cases were allocated to a larger cluster when PSI and WGS data were combined compared with PSI data alone. Mixed clusters containing women, men who report sex with women, and men who report sex with men were common when using the WGS data either alone or in combination with the PSI data. CONCLUSIONS: Combining PSI and WGS data improves our understanding of sexual network connectivity.

BACKGROUND: The prevalence of Neisseria gonorrhoeae (GC) isolates with elevated minimum inhibitory concentrations (MICs) to various antibiotics continues to rise in the U.S. and globally. Genomic analysis provides a powerful tool for surveillance of circulating strains, antimicrobial resistance determinants, and understanding of transmission through a population. METHODS: GC isolates collected from the U.S. Gonococcal Isolate Surveillance Project (GISP) in 2018 (n=1479) were sequenced and characterized. Whole genome sequencing was used to identify sequence types, antimicrobial resistance profiles, and phylogenetic relationships across demographic and geographic populations. RESULTS: Genetic characterization identified that (1) 80% of the GC isolates were represented in 33 multilocus sequence types, (2) isolates clustered in 23 major phylogenetic clusters with select phenotypic and demographic prevalence, and (3) common antimicrobial resistance determinants associated with low-level or high-level decreased susceptibility or resistance to relevant antibiotics. CONCLUSIONS: Characterization of this 2018 GISP genomic dataset, which is the largest U.S. whole genome sequence data set to date, sets the basis for future prospective studies, and establishes a genomic baseline of GC populations for local and national monitoring.

We would like to extend our gratitude to Drs. Deng and Klausner for their interest in our article, in which we phylogenetically characterized gonococcal isolates collected through national sentinel surveillance and circulating in the United States, between 2014 and 2016 [1]. We identified 2 major subpopulations of strains associated with reduced susceptibility to either azithromycin (multilocus sequence typing [MLST] ST9363) or cephalosporins (MLST ST1901), and we detailed the evolution of several strains that possessed mutations that were not observed in the United States in 2000–2013. | | Antimicrobial resistance (AMR) in the gonococcus represents a major public health threat, chiefly due to its rapidly evolving nature and the increasingly limited number of available treatment options. Although routine surveillance of AMR in circulating strains is essential in monitoring this threat, an equally important practice involves the promotion of antibiotic stewardship in limiting the spread of resistance. It is notable that the Centers for Disease Control and Prevention’s Gonococcal Isolate Surveillance Project, a historic sentinel surveillance program, has tracked AMR trends for over 30 years and, more recently, included whole-genome sequencing (WGS) of a subset of isolates [2]. As a result, the implementation of WGS has generated massive quantities of genomic data that can be used in tandem with data from traditional antimicrobial susceptibility testing to examine AMR trends. In addition, the release of these data to the public repositories has facilitated the ability of other investigators to conduct a myriad of follow-up studies.

BACKGROUND: The number of cases of gonorrhoea in the USA and worldwide caused by Neisseria gonorrhoeae is increasing (555 608 reported US cases in 2017, and 87 million cases worldwide in 2016). Many countries report declining in vitro susceptibility of azithromycin, which is a concern because azithromycin and ceftriaxone are the recommended dual treatment in many countries. We aimed to identify strain types associated with decreased susceptibility to azithromycin. METHODS: We did a genomic analysis of N gonorrhoeae isolates obtained by the US Gonococcal Isolate Surveillance Project. Isolates were whole-genome sequenced based on decreased susceptibility to azithromycin (minimal inhibitory concentration [MIC] ≥2 μg/mL, using agar dilution antibiotic susceptibility testing) and geographical representation. Bioinformatic analyses established genomic diversity, strain population dynamics, and antimicrobial resistance profiles. FINDINGS: 410 isolates were sorted into more than 20 unique phylogenetic clades. One predominant persistent clade (consisting of 97 isolates) included the most isolates with azithromycin MICs of 2 μg/mL or higher (61 of 97 [63%] vs 59 of 311 [19%]; p<0·0001) and carried a mosaic mtr (multiple transferable resistance) locus (68 of 97 [70%] vs two of 313 [1%]; p<0·0001). Of the remaining 313 isolates, 57 (18%) had decreased susceptibility to azithromycin (MIC ≥4 μg/mL), which was attributed to 23S rRNA variants (56 of 57 [98%]) and formed phylogenetically diverse clades, showing various levels of clonal expansion. INTERPRETATION: Reduced azithromycin susceptibility was associated with expanding and persistent clades harbouring two well described resistance mechanisms, mosaic mtr locus and 23S rRNA variants. Understanding the role of recombination, particularly within the mtr locus, on the fitness and expansion of strains with decreased susceptibility has important implications for the public health response to minimise gonorrhoea transmission. FUNDING: US Centers for Disease Control and Prevention (CDC), CDC Combating Antibiotic Resistant Bacteria initiative, Oak Ridge Institute for Science Education, US Department of Energy/CDC/Emory University, National Institutes of Health, and Biomedical Laboratory Research and Development Service of the US Department of Veterans Affairs.

In 2016, the proportion of Neisseria gonorrhoeae isolates with reduced susceptibility to azithromycin rose to 3.6%. A phylogenetic analysis of 334 N. gonorrhoeae isolates collected in 2016 revealed a single, geographically diverse lineage of isolates with MICs of 2-16 mug/mL that carried a mosaic-like mtr locus, whereas the majority of isolates with MICs >/= 16 mug/mL appeared sporadically and carried 23S rRNA mutations. Continued molecular surveillance of N. gonorrheae will identify new resistance mechanisms.

US gonorrhea rates are rising, and antibiotic-resistant Neisseria gonorrhoeae (AR-Ng) is an urgent public health threat. Since implementation of nucleic acid amplification tests for Ng identification, capacity for culturing Ng in the US has declined, along with the ability to perform culture-based antimicrobial susceptibility testing (AST). Yet, AST is critical for detecting and monitoring AR-Ng. In 2016, CDC established the Antibiotic Resistance Laboratory Network (AR Lab Network) to shore up national capacity for detecting several resistance threats including Ng. AR-Ng testing, a sub-activity of CDC's AR Lab Network, is performed in a tiered network of approximately 35 local laboratories, four regional laboratories (state public health laboratories in MD, TN, TX, WA), and CDC's national reference laboratory. Local laboratories receive specimens from approximately 60 clinics associated with the Gonococcal Isolate Surveillance Project (GISP), enhanced GISP (eGISP), and Strengthening the U.S. Response to Resistant Gonorrhea (SURRG). They isolate and ship up to 20,000 isolates to regional laboratories for culture-based agar dilution AST with seven antibiotics and for whole genome sequencing of up to 5,000 isolates. The CDC further examines concerning isolates and monitors genetic AR markers. During 2017 and 2018, the network tested 8,214 and 8,628 Ng isolates, and CDC received 531 and 646 concerning isolates, and 605 and 3,159 sequences, respectively. In summary, the AR Lab Network supported laboratory capacity for Ng-AST and associated genetic marker detection, expanding pre-existing notification and analysis systems for resistance detection. Continued, robust AST and genomic capacity can help inform national public health monitoring and intervention.

The human pathogens Neisseria gonorrhoeae and Neisseria meningitidis share high genome identity. Retrospective analysis of surveillance data from New Zealand indicates the potential cross-protective effect of outer membrane vesicle (OMV) meningococcal serogroup B vaccine (MeNZB) against N. gonorrhoeae A licensed OMV-based MenB vaccine, MenB-4C, consists of a recombinant FHbp, NhbA, NadA, and the MeNZB OMV. Previous work has identified several abundantly expressed outer membrane proteins (OMPs) as major components of the MenB-4C OMV with high sequence similarity between N. gonorrhoeae and N. meningitidis, suggesting a mechanism for cross-protection. To build off these findings, we performed comparative genomic analysis on 970 recent N. gonorrhoeae isolates collected through a U.S surveillance system against N. meningitidis serogroup B (NmB) reference sequences. We identified 1,525 proteins that were common to both Neisseria species, of which 57 proteins were predicted to be OMPs using in silico methods. Among the MenB-4C antigens, NhbA showed moderate sequence identity (73%) to the respective gonococcal homolog, was highly conserved within N. gonorrhoeae, and was predicted to be surface expressed. In contrast, the gonococcal FHbp was predicted not to be surface expressed, while NadA was absent in all N. gonorrhoeae isolates. Our work confirmed recent observations (E. A. Semchenko, A. Tan, R. Borrow, and K. L. Seib, Clin Infect Dis, 2018, https://doi.org/10.1093/cid/ciy1061) and describes homologous OMPs from a large panel of epidemiologically relevant N. gonorrhoeae strains in the United States against NmB reference strains. Based on our results, we report a set of OMPs that may contribute to the previously observed cross-protection and provide potential antigen targets to guide the next steps in gonorrhea vaccine development.IMPORTANCE Gonorrhea, a sexually transmitted disease, causes substantial global morbidity and economic burden. New prevention and control measures for this disease are urgently needed, as strains resistant to almost all classes of antibiotics available for treatment have emerged. Previous reports demonstrate that cross-protection from gonococcal infections may be conferred by meningococcal serogroup B (MenB) outer membrane vesicle (OMV)-based vaccines. Among 1,525 common proteins shared across the genomes of both N. gonorrhoeae and N. meningitidis, 57 proteins were predicted to be surface expressed (outer membrane proteins [OMPs]) and thus preferred targets for vaccine development. The majority of these OMPs showed high sequence identity between the 2 bacterial species. Our results provide valuable insight into the meningococcal antigens present in the current OMV-containing MenB-4C vaccine that may contribute to cross-protection against gonorrhea and may inform next steps in gonorrhea vaccine development.

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 from 2014 to 2016. METHODS: We sequenced the genomes of 649 isolates collected through the Gonococcal Isolate Surveillance Project (GISP). 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 MICs (AZI em), characterized by a mosaic mtr locus (C-substitution in the mtrR promoter, mosaic mtrR and mtrD). Mutations in 23S rRNA were sporadically distributed among AZI em strains. Another clonal group (MLST ST1901) possessed seven 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 within a patient.

Recent reports suggest that mosaic-like sequences within the mtr (multiple transferable resistance) efflux pump locus of Neisseria gonorrhoeae, likely originating from commensal Neisseria sp. by transformation, can increase the ability of gonococci to resist structurally diverse antimicrobials. Thus, acquisition of numerous nucleotide changes within the mtrR gene encoding the transcriptional repressor (MtrR) of the mtrCDE efflux pump-encoding operon or overlapping promoter region for both along with those that cause amino acid changes in the MtrD transporter protein were recently reported to decrease gonococcal susceptibility to numerous antimicrobials, including azithromycin (Azi) (C. B. Wadsworth, B. J. Arnold, M. R. A. Satar, and Y. H. Grad, mBio 9:e01419-18, 2018, https://doi.org/10.1128/mBio.01419-18). We performed detailed genetic and molecular studies to define the mechanistic basis for why such strains can exhibit decreased susceptibility to MtrCDE antimicrobial substrates, including Azi. We report that a strong cis-acting transcriptional impact of a single nucleotide change within the -35 hexamer of the mtrCDE promoter as well gain-of-function amino acid changes at the C-terminal region of MtrD can mechanistically account for the decreased antimicrobial susceptibility of gonococci with a mosaic-like mtr locus.IMPORTANCE Historically, after introduction of an antibiotic for treatment of gonorrhea, strains of N. gonorrhoeae emerge that display clinical resistance due to spontaneous mutation or acquisition of resistance genes. Genetic exchange between members of the Neisseria genus occurring by transformation can cause significant changes in gonococci that impact the structure of an antibiotic target or expression of genes involved in resistance. The results presented here provide a framework for understanding how mosaic-like DNA sequences from commensal Neisseria that recombine within the gonococcal mtr efflux pump locus function to decrease bacterial susceptibility to antimicrobials, including antibiotics used in therapy of gonorrhea.

Increasingly applied to identify mutations conferring antimicrobial resistance (AMR), disease outbreaks, and pathways of disease spread, whole genome sequencing (WGS)—the process of determining the complete DNA sequence of an organism’s genome at a single time—has emerged as a powerful tool for public health. Genomic analyses played central roles in recent outbreak investigations, such as of a high-profile outbreak of carbapenem-resistant Klebsiella pneumoniae at the US National Institutes of Health Clinical Center, the 2010 outbreak of cholera in Haiti, the 2014–2015 HIV outbreak in Indiana, the epidemic of Zika virus in the Americas, and large outbreaks of foodborne and waterborne illness.1–7 Whole genome sequencing findings have informed development of novel molecular diagnostics and explorations of human microbiomes.8,9 Whereas DNA sequencing methods were painstakingly performed manually decades ago, the development of automated methods in the 1990s, followed by rapidly accelerating speed of sequencing, plummeting cost, increasing computational capacity, growing number of sequences in publically available repositories (e.g., GenBank), and increasing availability of bioinformatics tools in the past decade, have supported a dramatic expansion of WGS.