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.

Neisseria spp. possess four genogroups of filamentous prophages, termed Nf1 to 4. A filamentous bacteriophage from the Nf1 genogroup termed meningococcal disease-associated phage (MDA φ) is associated with clonal complexes of Neisseria meningitidis that cause invasive meningococcal disease. Recently, we recovered an isolate of Neisseria gonorrhoeae (ExNg63) from a rare case of gonococcal meningitis, and found that it possessed a region with 90% similarity to Nf1 prophages, specifically, the meningococcal MDA φ. This led to the hypothesis that the Nf1 prophage may be more widely distributed amongst the genus Neisseria. An analysis of 92 reference genomes revealed the presence of intact Nf1 prophages in the commensal species, Neisseria lactamica and Neisseria cinerea in addition to the pathogen N. gonorrhoeae. In N. gonorrhoeae, Nf1 prophages had a restricted distribution but were present in all representatives of MLST ST1918. Of the 160 phage integration sites identified, only one common insertion site was found between one isolate of N. gonorrhoeae and N. meningitidis. There was an absence of any obvious conservation of the receptor for prophage entry, PilE, suggesting that the phage may have been obtained by natural transformation. An examination of the restriction modification systems and mutated mismatch repair systems with prophage presence suggested that there was no obvious preference for these hosts. A timed phylogeny inferred that N. meningitidis was the donor of the Nf1 prophages in N. lactamica and N. gonorrhoeae. Further work is required to determine whether Nf1 prophages are active and can act as accessory colonization factors in these species.

In Australia, gonococcal isolates are monitored for antimicrobial susceptibilities. In Western Australia (WA), gonorrhoea notification rates increased by 63 % between 2013 and 2016, with the steepest increase occurring between 2015 and 2016, before stabilizing at this higher baseline between 2017 and 2020. This increased prevalence was associated with antimicrobial-susceptible (AMS) lineages. To understand the provenance of these isolates causing gonorrhoea in WA, whether they were introduced or expanded from endogenous lineages, 741 isolates were collected in 2017 and characterized by both iPLEX typing and whole genome sequencing (WGS). Antibiograms and genocoding of the isolates revealed that AMS isolates were most prevalent in the remote regions, while the urban/rural regions were characterized by antimicrobial-resistant (AMR) isolates. iPLEX typing identified 78 iPLEX genotypes (WA-1 to WA-78) of which 20 accounted for over 88 % of isolates. WA-10 was the most frequently identified genotype in the urban/rural regions whilst WA-29 was the most frequently identified genotype in the remote regions. Genotypes WA-38, WA-52 and WA-13 accounted for 81 % (n=36/44) of the azithromycin-resistant N. gonorrhoeae (AziR) isolates. A representative isolate of each iPLEX genotype and AMR biotype was whole genome sequenced and analysed using MLST, NG-MAST and NG-STAR, and the novel core genome clustering Ng_cgc_400 typing scheme. Five predominant Bayesian population groups (termed BPG-1 to 5) were identified in the study collection. BPG-1 and BPG-2 were associated with AMS isolates from the remote regions. BPG-1 and BPG-2 were shown to be unique to the remote regions based on a minimum spanning tree against 4000 international isolates. AMS isolates in urban/rural regions were dominated by international lineages. AziR and Cef DS (decreased susceptibility to ceftriaxone) was concentrated in three urban/rural genomic groups (BPG-3, 4 and 5). Azithromycin minimum inhibitory concentrations (0.5–16 mg l−1) correlated with the accumulation of mtrR mutations or/and the fraction of 23S rRNA C2611T mutated copies. The majority of isolates in BPG-3, 4 and 5 could be correlated with known AMR lineages circulating globally and nationally. In conclusion, the surge in AMS isolates in WA in 2017 was due to importation of international AMS lineages into urban/rural regions, whilst the local AMS lineages persisted largely in the remote regions. Bridging between the urban/rural and remote regions was relatively rare, but continued surveillance is required to prevent ingress of AMR strains/lineages into the remote regions of WA.