Antibiotic Treatment Regimes as a Driver of the Global Population Dynamics of a Major Gonorrhea Lineage

Cyr SB, Schlanger K, Raphael BH, Kersh EN, Joseph SJ. Emergence and evolution of mosaic penA-60 and penA-237 alleles in a Neisseria gonorrhoeae core genogroup that was historically susceptible to extended spectrum cephalosporins

Introduction

Neisseria gonorrhoeae (Ng) has successively developed resistance to all previously recommended antimicrobial therapies, with ceftriaxone being the last option for monotherapy of gonorrhea. Global emergence and international spread of the FC428 clone derived mosaic penA-60 allele, associated with highlevel ceftriaxone minimum inhibitory concentrations (MICs) in non FC428 clone Ng lineages, has become an increasing concern. The penA-60 allele carrying Ng was first identified in the U.S. in Las Vegas, Nevada (2019; GCWGS-102723), with a multi-locus sequence type (MLST)-1901 strain, in a non FC428 clone Ng lineage, which is associated with a historically ceftriaxone susceptible core genogroup. Later in 2022, an allele genetically similar to penA-60, mosaic penA-237, was identified in the UK (H22-722) and France (F92) with high-level ceftriaxone MICs and both belonged to MLST-1901.

Methods

In this study, we assessed phylogenomic relatedness and antimicrobial resistance (AMR) determinant profiles of these three isolates with high-level ceftriaxone MICs among a global collection of 2,104 genomes belonging to the MLST-1901 core genome cluster group 31, which includes strains separated by a locus threshold of 200 or fewer differences (Ng_cgc_200). Recombination events in and around the penA coding region were catalogued and potential sources of inter species recombinant DNA were also inferred.

Results

The global population structure of MLST-1901 core genogroup falls into 4 major lineages. Isolates GCWGS-10723, F92, and H22-722 clustered within Lineage 1, which was dominated by non-mosaic penA-5 alleles. These three isolates formed a clade within Lineage 1 that consisted of isolates from North America and southeast Asia. Neisseria subflava and Neisseria sicca were identified as likely progenitors of two independent recombination events that may have led to the generation of mosaic penA-60 and penA-237, within a possible non-mosaic penA-5 background.

Discussions

Our study suggests that there are multiple evolutionary pathways that could generate concerning mosaic penA alleles via homologous recombination of historically susceptible Ng lineages with Neisseria commensals. Enhanced surveillance of gonococcal strains and Neisseria commensals is crucial for understanding of the evolution of AMR, particularly in less-studied regions (e.g., Asia), where high-level ceftriaxone MICs and multi-drug resistance are more prevalent.

Neisseria gonorrhoeae ST-1901 in Rio de Janeiro from 2006 to 2022: Phylogeny and antimicrobial resistance evolution of a well-succeeded pathogen

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 four of five clades. In contrast, β-lactam resistance associated penA 34 mosaic was found only in one persisting clade (Clade D), and unique G45D and A39T mutations in mtrR gene and its promoter (Nm-Like) were found only in 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, namely, ceftriaxone and azithromycin.

Horizontal Gene Transfer of Fluoroquinolone Resistance-Conferring Genes From Commensal Neisseria to Neisseria gonorrhoeae: A Global Phylogenetic Analysis of 20,047 Isolates

Antimicrobial resistance in Neisseria gonorrhoeae is an important global health concern. The genetically related commensal Neisseria act as a reservoir of resistance genes, and horizontal gene transfer (HGT) has been shown to play an important role in the genesis of resistance to cephalosporins and macrolides in N. gonorrhoeae. In this study, we evaluated if there was evidence of HGT in the genes gyrA/gyrB and parC/parE responsible for fluoroquinolone resistance. Even though the role of gyrB and parE in quinolone resistance is unclear, the subunits gyrB and parE were included as zoliflodacin, a promising new drug to treat N. gonorrhoeae targets the gyrB subunit. We analyzed a collection of 20,047 isolates; 18,800 N. gonorrhoeae, 1,238 commensal Neisseria spp., and nine Neisseria meningitidis. Comparative genomic analyses identified HGT events in genes, gyrA, gyrB, parC, and parE. Recombination events were predicted in N. gonorrhoeae and Neisseria commensals. Neisseria lactamica, Neisseria macacae, and Neisseria mucosa were identified as likely progenitors of the HGT events in gyrA, gyrB, and parE, respectively.

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.

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.

Donnan Potential across the Outer Membrane of Gram-Negative Bacteria and Its Effect on the Permeability of Antibiotics

The cell envelope structure of Gram-negative bacteria is unique, composed of two lipid bilayer membranes and an aqueous periplasmic space sandwiched in between. The outer membrane constitutes an extra barrier to limit the exchange of molecules between the cells and the exterior environment. Donnan potential is a membrane potential across the outer membrane, resulted from the selective permeability of the membrane, which plays a pivotal role in the permeability of many antibiotics. In this review, we discussed factors that affect the intensity of the Donnan potential, including the osmotic strength and pH of the external media, the osmoregulated periplasmic glucans trapped in the periplasmic space, and the displacement of cell surface charges. The focus of our discussion is the impact of Donnan potential on the cellular permeability of selected antibiotics including fluoroquinolones, tetracyclines, β-lactams, and trimethoprim.

Evidence of Horizontal Gene Transfer of 50S Ribosomal Genes rplB, rplD, and rplY in Neisseria gonorrhoeae

Horizontal gene transfer (HGT) in the penA and multidrug efflux pump genes has been shown to play a key role in the genesis of antimicrobial resistance in Neisseria gonorrhoeae. In this study, we evaluated if there was evidence of HGT in the genes coding for the ribosomal proteins in the Neisseria genus. We did this in a collection of 11,659 isolates of Neisseria, including N. gonorrhoeae and commensal Neisseria species (N. cinerea, N. elongata, N. flavescens, N. mucosa, N. polysaccharea, and N. subflava). Comparative genomic analyses identified HGT events in three genes: rplB, rplD, and rplY coding for ribosomal proteins L2, L4 and L25, respectively. Recombination events were predicted in N. gonorrhoeae and N. cinerea, N. subflava, and N. lactamica were identified as likely progenitors. In total, 2,337, 2,355, and 1,127 isolates possessed L2, L4, and L25 HGT events. Strong associations were found between HGT in L2/L4 and the C2597T 23S rRNA mutation that confers reduced susceptibility to macrolides. Whilst previous studies have found evidence of HGT of entire genes coding for ribosomal proteins in other bacterial species, this is the first study to find evidence of HGT-mediated chimerization of ribosomal proteins.

Choosing New Therapies for Gonorrhoea: We Need to Consider the Impact on the Pan-Neisseria Genome. A Viewpoint

The development of new gonorrhoea treatment guidelines typically considers the resistance-inducing effect of the treatment only on Neisseria gonorrhoeae. Antimicrobial resistance in N. gonorrhoeae has, however, frequently first emerged in commensal Neisseria species and then been passed on to N. gonorrhoeae via transformation. This creates the rationale for considering the effect of gonococcal therapies on resistance in commensal Neisseria. We illustrate the benefits of this pan-Neisseria strategy by evaluating three contemporary treatment options for N. gonorrhoeae-ceftriaxone plus azithromycin, monotherapy with ceftriaxone and zoliflodacin.