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 meningitidis (Nm) and Neisseria gonorrhoeae (Ng) are human pathogens that sometimes occupy the same anatomical niche. Ng, the causative agent of gonorrhea, infects 87 million individuals annually worldwide and is an urgent threat due to increasing drug resistance. Ng is a pathogen of the urogenital tract and may infect the oropharyngeal or rectal site, often asymptomatically. Conversely, Nm is an opportunistic pathogen. While often a commensal in the oropharyngeal tract, it is also the leading cause of bacterial meningitis with 1.2 million cases globally, causing significant morbidity and mortality. Horizontal gene transfer (HGT) is likely to occur between Ng and Nm due to their shared anatomical niches and genetic similarity, which poses challenges for accurate detection and treatment. Routine surveillance through the Gonococcal Isolate Surveillance Project and Strengthening the U.S. Response to Resistant Gonorrhea detected six concerning urogenital Neisseria isolates with contradicting species identification in Milwaukee (MIL). While all six isolates were positive for Ng using nucleic acid amplification testing (NAAT) and matrix-assisted laser desorption/ionization time of flight identified the isolates as Ng, two biochemical tests, Gonochek-II and API NH, classified them as Nm. To address this discrepancy, we performed whole-genome sequencing (WGS) using Illumina MiSeq on all isolates and employed various bioinformatics tools. Species detection analysis using BMScan, which uses WGS data, identified all isolates as Ng. Furthermore, Kraken revealed over 98% of WGS reads mapped to the Ng genome and <1% to Nm. Recombination analysis identified putative HGT in all MIL isolates within the γ-glutamyl transpeptidase (ggt) gene, a key component in the biochemical tests used to differentiate between Nm and Ng. Further analysis identified Nm as the source of HGT event. Specifically, the active Nm ggt gene replaced the Ng pseudogenes, ggt1 and ggt2. Together, this study demonstrates that closely related Neisseria species sharing a niche underwent HGT, which led to the misidentification of species following biochemical testing. Importantly, NAAT accurately detected Ng. The misidentification highlights the importance of using WGS to continually evaluate diagnostic or bacterial identification tests.

Unsuccessful treatment of gonorrhoea has not yet occurred in the USA, and cases of gonorrhoea that are non-susceptible to cephalosporins have been rare. In 2019, non-susceptibility to ceftriaxone conferred by the mosaic penA 60.001 allele was found in a Neisseria gonorrhoeae multilocus sequence type (MLST) 1901 isolate from Nevada.1 In this Correspondence, we present two additional US cases of the penA 60.001 allele identified in MLST 8123, an emerging international multidrug non-susceptible N gonorrhoeae lineage. Although these cases responded to ceftriaxone treatment, N gonorrhoeae isolates from the first known patient (case 1) demonstrated in-vitro non-susceptibility to ceftriaxone as well as non-susceptibility or resistance to drugs previously recommended for front-line treatment. | | In August, 2022, N gonorrhoeae grown from urine culture from a patient with urethritis in primary care in Massachusetts displayed non-susceptibility to cephalosporins (the minimum inhibitory concentrations were 1·0 μg/mL for ceftriaxone and >1·0 μg/mL for cefixime by agar dilution; the minimum inhibitory concentration for cefixime was 1·5 μg/mL by gradient strip) and azithromycin and resistance to ciprofloxacin, penicillin, and tetracycline (appendix pp 6–7). Antimicrobial susceptibility testing was done with gradient strips at the state public health laboratory Massachusetts and then confirmed via agar dilution at the US Centers for Disease Control and Prevention (CDC). The patient (case 1) had already been successfully diagnosed on nucleic acid amplification test (NAAT) with gonorrhoea and was given 500 mg ceftriaxone intramuscularly and asked to return to primary care where, 9 days after treatment, he was asymptomatic, had normal results during examination, and tested negative by urine culture and pharyngeal and rectal NAAT recommended by the Massachusetts sexually transmitted diseases programme to document N gonorrhoeae clearance from any site of infection. The patient reported that he had not travelled outside USA in the 60 days before onset of symptoms. He disclosed female sex worker contacts, but insufficient information was provided to trace the contacts.

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.

Chlamydia trachomatis is the most prevalent cause of bacterial sexually transmitted infections (STIs) worldwide. U.S. cases have been steadily increasing for more than a decade in both the urogenital tract and rectum. C. trachomatis is an obligate intracellular bacterium that is not easily cultured, limiting the capacity for genome studies to understand strain diversity and emergence among various patient populations globally. While Agilent SureSelectXT target-enrichment RNA bait libraries have been developed for whole-genome enrichment and sequencing of C. trachomatis directly from clinical urine, vaginal, conjunctival and rectal samples, efficiencies are only 60-80% for ≥95-100% genome coverage. We therefore re-designed and expanded the RNA bait library to augment enrichment of the organism from clinical samples to improve efficiency. We describe the expanded library, the limit of detection for C. trachomatis genome copy input, and the 100% efficiency and high-resolution of generated genomes where genomic recombination among paired vaginal and rectal specimens from four patients was identified. This workflow provides a robust approach for discerning genomic diversity and advancing our understanding of the molecular epidemiology of contemporary C. trachomatis STIs across sample types, among geographic populations, sexual networks, and outbreaks associated with proctitis/proctocolitis among women and men who have sex with men.Importance Chlamydia trachomatis is an obligate intracellular bacterium that is not easily cultured, and there is limited information on rectal C. trachomatis transmission and its impact on morbidity. To improve efficiency of previous studies involving whole genome target enrichment and sequencing of C. trachomatis directly from clinical urine, vaginal, conjunctival, and rectal specimens, we expanded the RNA bait library to augment enrichment of the organism from clinical samples. We demonstrate an increased efficiency in the percentage of reads mapping to C. trachomatis. We show the new system is sensitive for near identical genomes of C. trachomatis from two body sites in four women. Further, we provide a robust genomic epidemiologic approach to advance our understanding of C. trachomatis strains causing ocular, urogenital and rectal infections, and to explore geo-sexual networks, outbreaks of colorectal infections among women and men who have sex with men, and the role of these strains in morbidity.Competing Interest StatementThe authors have declared no competing interest.

Disseminated gonococcal infection (DGI) is a rare complication caused by the systemic dissemination of Neisseria gonorrhoeae to normally sterile anatomical sites. Little is known about the genetic diversity of DGI gonococcal strains and how they relate to other gonococcal strains causing uncomplicated mucosal infections. We used whole genome sequencing to characterize DGI isolates (n = 30) collected from a surveillance system in Georgia, United States, during 2017-2020 to understand phylogenetic clustering among DGI as well as uncomplicated uro- and extragenital gonococcal infection (UGI) isolates (n = 110) collected in Fulton County, Georgia, during 2017-2019. We also investigated the presence or absence of genetic markers related to antimicrobial resistance (AMR) as well as surveyed the genomes for putative virulence genetic factors associated with normal human-serum (NHS) resistance that might facilitate DGI. We found that DGI strains demonstrated significant genetic variability similar to the population structure of isolates causing UGI, with sporadic incidences of geographically clustered DGI strains. DGI isolates contained various AMR markers and genetic mechanisms associated with NHS resistance. DGI isolates had a higher frequency of the porB1A allele compared with UGI (67% vs 9%, P < .0001); however, no single NHS resistance marker was found in all DGI isolates. Continued DGI surveillance with genome-based characterization of DGI isolates is necessary to better understand specific factors that promote systemic dissemination.

High-throughput sequencing has become commonplace in evolutionary studies. Large, rapidly collected genomic datasets are used to capture biodiversity and for monitoring global and national scale disease transmission patterns, among many other applications. Updating homologous sequence datasets with new samples is cumbersome, requiring excessive program runtimes and data processing. We describe Extensiphy, a bioinformatics tool to efficiently update multiple sequence alignments with whole-genome short-read data. Extensiphy performs reference based sequence assembly and alignment in one process while maintaining the alignment length of the original alignment. Input data-types for Extensiphy are any multiple sequence alignment in fasta format and whole-genome, short-read fastq sequences. To validate Extensiphy, we compared its results to those produced by two other methods that construct whole-genome scale multiple sequence alignments. We measured our comparisons by analysing program runtimes, base-call accuracy, dataset retention in the presence of missing data and phylogenetic accuracy. We found that Extensiphy rapidly produces high-quality updated sequence alignments while preventing alignment shrinkage due to missing data. Phylogenies estimated from alignments produced by Extensiphy show similar accuracy to other commonly used alignment construction methods. Extensiphy is suitable for updating large sequence alignments and is ideal for studies of biodiversity, ecology and epidemiological monitoring efforts. © 2021 The Authors. Methods in Ecology and Evolution published by John Wiley & Sons Ltd on behalf of British Ecological Society

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.

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.

Chlamydia trachomatis, an obligately intracellular bacterium, is the most prevalent cause of bacterial sexually transmitted infections (STIs) worldwide. Numbers of U.S. infections of the urogenital tract and rectum have increased annually. Because C. trachomatis is not easily cultured, comparative genomic studies are limited, restricting our understanding of strain diversity and emergence among populations globally. While Agilent SureSelect(XT) target enrichment RNA bait libraries have been developed for whole-genome enrichment and sequencing of C. trachomatis directly from clinical urine, vaginal, conjunctival, and rectal samples, public access to these libraries is not available. We therefore designed an RNA bait library (34,795 120-mer probes based on 85 genomes, versus 33,619 probes using 74 genomes in a previous one) to augment organism sequencing from clinical samples that can be shared with the scientific community, enabling comparison studies. We describe the library and limit of detection for genome copy input, and we present results of 100% efficiency and high-resolution determination of recombination and identical genomes within vaginal-rectal specimen pairs in women. This workflow provides a robust approach for discerning genomic diversity and advancing our understanding of the molecular epidemiology of contemporary C. trachomatis STIs across sample types, geographic populations, sexual networks, and outbreaks associated with proctitis/proctocolitis among women and men who have sex with men.IMPORTANCE Chlamydia trachomatis is an obligate intracellular bacterium that is not easily cultured, which limits our understanding of urogenital and rectal C. trachomatis transmission and impact on morbidity. To provide a publicly available workflow for whole-genome target enrichment and sequencing of C. trachomatis directly from clinical urine, vaginal, conjunctival, and rectal specimens, we developed and report on an RNA bait library to enrich the organism from clinical samples for sequencing. We demonstrate an increased efficiency in the percentage of reads mapping to C. trachomatis and identified recombinant and identical C. trachomatis genomes in paired vaginal-rectal samples from women. Our workflow provides a robust genomic epidemiologic approach to advance our understanding of C. trachomatis strains causing ocular, urogenital, and rectal infections and to explore geo-sexual networks, outbreaks of colorectal infections among women and men who have sex with men, and the role of these strains in morbidity.

Contaminants such as heavy metals may contribute to the dissemination of antimicrobial resistance (AMR) by enriching resistance gene determinants via co-selection mechanisms. In the present study, a survey was performed on soils collected from four areas at the Savannah River Site (SRS), South Carolina, USA, with varying contaminant profiles: relatively pristine (Upper Three Runs), heavy metals (Ash Basins), radionuclides (Pond B) and heavy metal and radionuclides (Tim's Branch). Using 16S rRNA gene amplicon sequencing, we explored the structure and diversity of soil bacterial communities. Sites with legacies of metal and/or radionuclide contamination displayed significantly lower bacterial diversity compared to the reference site. Metagenomic analysis indicated that multidrug and vancomycin antibiotic resistance genes (ARGs) and metal resistance genes (MRGs) including those associated with copper, arsenic, iron, nickel and zinc were prominent in all soils including the reference site. However, significant differences were found in the relative abundance and diversity of certain ARGs and MRGs in soils with metal/radionuclide contaminated soils compared to the reference site. Co-occurrence patterns revealed significant ARG/MRG subtypes in predominant soil taxa including Acidobacteriaceae, Bradyrhizobium, Mycobacterium, Streptomyces, Verrumicrobium, Actinomadura and Solirubacterales. Overall, the study emphasizes the potential risk of human activities on the dissemination of AMR in the environment.

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.

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.