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 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.

Next-generation sequencing (NGS) of amplicons is used in a wide variety of contexts. Most NGS amplicon sequencing remains overly expensive and inflexible, with library preparation strategies relying upon the fusion of locus-specific primers to full-length adapter sequences with a single identifying sequence or ligating adapters onto PCR products. In Adapterama I, we presented universal stubs and primers to produce thousands of unique index combinations and a modifiable system for incorporating them into Illumina libraries. Here, we describe multiple ways to use the Adapterama system and other approaches for amplicon sequencing on Illumina instruments. In the variant we use most frequently for large-scale projects, we fuse partial adapter sequences (TruSeq or Nextera) onto the 5’ end of locus-specific PCR primers with variable-length tag sequences between the adapter and locus-specific sequences. These fusion primers can be used combinatorially to amplify samples within a 96-well plate (eight forward primers + 12 reverse primers yield 8 × 12 = 96 combinations), and the resulting amplicons can be pooled. The initial PCR products then serve as template for a second round of PCR with dual-indexed iTru or iNext primers (also used combinatorially) to make full-length libraries. The resulting quadruple-indexed amplicons have diversity at most base positions and can be pooled with any standard Illumina library for sequencing. The number of sequencing reads from the amplicon pools can be adjusted, facilitating deep sequencing when required or reducing sequencing costs per sample to an economically trivial amount when deep coverage is not needed. We demonstrate the utility and versatility of our approaches with results from six projects using different implementations of our protocols. Thus, we show that these methods facilitate amplicon library construction for Illumina instruments at reduced cost with increased flexibility. A simple web page to design fusion primers compatible with iTru primers is available at: http://baddna.uga.edu/tools-taggi.html. A fast and easy to use program to demultiplex amplicon pools with internal indexes is available at: https://github.com/lefeverde/Mr_Demuxy.

Next-generation sequencing (NGS) of amplicons is used in a wide variety of contexts. In many cases, NGS amplicon sequencing remains overly expensive and inflexible, with library preparation strategies relying upon the fusion of locus-specific primers to full-length adapter sequences with a single identifying sequence or ligating adapters onto PCR products. In Adapterama I, we presented universal stubs and primers to produce thousands of unique index combinations and a modifiable system for incorporating them into Illumina libraries. Here, we describe multiple ways to use the Adapterama system and other approaches for amplicon sequencing on Illumina instruments. In the variant we use most frequently for large-scale projects, we fuse partial adapter sequences (TruSeq or Nextera) onto the 5' end of locus-specific PCR primers with variable-length tag sequences between the adapter and locus-specific sequences. These fusion primers can be used combinatorially to amplify samples within a 96-well plate (8 forward primers + 12 reverse primers yield 8 × 12 = 96 combinations), and the resulting amplicons can be pooled. The initial PCR products then serve as template for a second round of PCR with dual-indexed iTru or iNext primers (also used combinatorially) to make full-length libraries. The resulting quadruple-indexed amplicons have diversity at most base positions and can be pooled with any standard Illumina library for sequencing. The number of sequencing reads from the amplicon pools can be adjusted, facilitating deep sequencing when required or reducing sequencing costs per sample to an economically trivial amount when deep coverage is not needed. We demonstrate the utility and versatility of our approaches with results from six projects using different implementations of our protocols. Thus, we show that these methods facilitate amplicon library construction for Illumina instruments at reduced cost with increased flexibility. A simple web page to design fusion primers compatible with iTru primers is available at: http://baddna.uga.edu/tools-taggi.html. A fast and easy to use program to demultiplex amplicon pools with internal indexes is available at: https://github.com/lefeverde/Mr_Demuxy.

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.