- Research article
- Open Access
Genomic epidemiology and population structure of Neisseria gonorrhoeae from remote highly endemic Western Australian populations
BMC Genomics volume 19, Article number: 165 (2018)
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.
Population structure analysis of the 59 isolates together with 72 isolates from an international collection, revealed six population groups suggesting that N. gonorrhoeae is a weakly clonal species. Two distinct population groups, Aus1 and Aus2, represented 63% of WA isolates and were mostly composed of the remote community isolates that carried no chromosomal AMR genotypes. In contrast, the Western Australian metropolitan isolates were frequently multi-drug resistant and belonged to population groups found in the international database, suggesting international transmission of the isolates.
Our study suggests that the population structure of N. gonorrhoeae is distinct between the communities in remote and metropolitan WA. Given the high rate of AMR in metropolitan regions, ongoing surveillance is essential to ensure the enduring efficacy of the empiric gonorrhoea treatment in remote WA.
Neisseria gonorrhoeae (gonococcus), the causative agent of gonorrhoea, is one of the most common causes of bacterial-associated sexually transmitted infection (STI). Estimated to cause 78 million new infections worldwide, 35 million cases of N. gonorrhoeae occur annually in the Western Pacific Region which includes Australia . Australian states and territories are divided into remote and metropolitan areas, based on the distance to the closest service centers and the population density. Communities living in remote regions of the NT and WA have the highest notification rate of gonorrhoea in Australia [2, 3]. In 2015, gonorrhoea notifications were reported to be 959 per 100,000 population of people living in the remote WA Kimberley region, compared to 85 per 100,000 population for Western Australia .
Gonococcal vaccines are currently not available and antimicrobial therapy is the only effective option for treating gonorrhoea. However, the World Health Organization (WHO) Global Gonococcal Antimicrobial Surveillance Programme (GASP) has reported the efficacy of many antibiotic treatments such as penicillin, tetracycline, and ciprofloxacin is lower than 95% indicating these antibiotics can no longer be used for empirical treatment in almost all GASP countries . This has resulted in a shift towards dual antimicrobial therapy, mainly ceftriaxone and azithromycin. Among WHO GASP countries, isolates with decreased susceptibility or resistance to the extended-spectrum cephalosporins (ESCs) and isolates with resistance to azithromycin have been reported in 51% and 75% of these countries respectively . To understand the epidemiology underlying the appearance and spread of antimicrobial resistance (AMR) in gonococci, several molecular approaches such as whole genome sequencing (WGS), multilocus sequence typing (MLST), N. gonorrhoeae multi-antigen sequence typing (NG-MAST) and most recently N. gonorrhoeae sequence typing for antimicrobial resistance (NG-STAR) have been developed [7,8,9]. Resistance to ESCs has spread primarily through clonal expansion and is highly correlated with the presence of the penA mosaic allele in NG-MAST sequence type (ST)1407, MLST STs 1901 and 7363 and NG-STAR ST90 clusters [8, 10, 11]. Similarly, emergent quinolone resistant strains have spread clonally in the NG-MAST ST225 lineage which is associated with gyrA and parC chromosomal mutations [11, 12]. Azithromycin resistance, which initially appeared sporadically, has rapidly spread through local sexual networks in geographically distinct regions . Multi-drug resistance (MDR) has been increasing clonally through a limited number of STs and phylogenetic clusters in the last 4 years with resistance to cephalosporins and azithromycin emerging a limited number of times [12, 14, 15].
The prevalence of AMR amongst gonococcal isolates in metropolitan areas of WA is similar to global levels, with 20% of isolates being penicillin and ciprofloxacin resistant; while ESC decreased susceptibility and azithromycin resistant isolates are sporadically reported . In contrast, the prevalence of AMR amongst isolates from the remote areas of WA is lower at 2–3% for penicillin and ciprofloxacin resistance while no resistance to ESCs or azithromycin has been detected [16, 17]. The lack of AMR in isolates from remote regions of Australia may be explained by the limited international contacts and by recent studies utilizing NG-MAST and MLST that have revealed gonococcal isolates in the remote areas are genetically distinct from metropolitan isolates which may be due to the different risk groups in each area: men who have sex with men (MSM) networks from the metropolitan areas versus the heterosexual indigenous communities living in remote areas [18, 19].
Our study aims to elucidate the genomic epidemiology and population structure of N. gonorrhoeae in WA using isolates collected from 2011 to 2013. Moreover, the study aims to determine the association between the N. gonorrhoeae types found in WA according to the molecular typing schemes and their antimicrobial susceptibility.
In WA, gonorrhoea is a Department of Health notifiable infectious disease and all isolates are referred to the Western Australian Gonococcal Surveillance Programme (WAGSP) laboratory for antimicrobial susceptibility testing. Antimicrobial resistance profiles of all isolates to penicillin, spectinomycin, azithromycin, ciprofloxacin, ceftriaxone and high-level resistance to tetracycline (tetracycline HLR) were assessed by the agar dilution method and interpreted using the Calibrated Dichotomous Sensitivity guidelines . Decreased susceptibility to ceftriaxone (0.06–0.126 mg/L) is confirmed by the E-test (bioMérieux, France). Penicillinase production is detected using the nitrocefin test (Oxoid, Australia). Anatomical isolation site, geographical location and postcode are available from the Communicable Disease Control Directorate.
Fifty nine N. gonorrhoeae isolates from patients living in the remote (n = 33) and metropolitan population centers (n = 26) of WA from 2011 and 2013 were obtained from the WAGSP laboratory . Isolates were stored in GC broth with 20% glycerol at − 80 °C and were passaged fewer than five times and were cultured under aerobic conditions with 5% CO2 at 37 °C on GC agar (Oxoid, Australia) supplemented with 0.4% glucose, 0.01% glutamine, 0.2 mg/L of cocarboxylase and 5 mg/L of iron (III) nitrate.
DNA extraction, genome sequencing and assembly
Genomic DNA extraction was performed using the DNeasy Blood and Tissue Kit (Qiagen, Germany) and the extract was stored as per the manufacturer’s instructions for DNA extraction from Gram negative bacteria. Genome sequencing of the 59 isolates was performed using the Illumina MiSeq platform (Illumina, USA) with 2 × 300 base pair read lengths. The targeted sequencing depth was 120 with a minimum Phred quality score of 30. Reads were de novo assembled using SPAdes genome assembler version 9.0 . The quality of the assembled genomes was assessed using the Quast genome assembly evaluation tool . Sequencing and assembly quality statistics of the 59 WA N. gonorrhoeae isolates are shown in Additional file 1. Bacterial Isolates Genome Sequence database (BIGSdb) genomics platform tools – hosted on www.pubmlst.org/neisseria – were used for annotation and genome wide analysis of the assembled isolates . The core genome consisted of 1427 genes that were present in all of the 59 N. gonorrhoeae isolates included in this analysis. The core genome was determined using the genome comparator tool at BIGSdb using the “all loci scheme” with 100% core threshold and excluding incomplete loci. Once the analysis was completed a list of complete variable genes that were present in all isolates was used as the core genome. To determine if the N. gonorrhoeae isolates circulating in remote WA are genetically different from those circulating in the WA metropolitan areas, 72 N. gonorrhoeae genomes from diverse international locations (available at PubMLST Neisseria spp. isolates database) from a recent epidemiological study by Ezewudo et al.  were included. A MAFFT alignment of the concatenated core genome of 1005 genes (excluding the invariant loci of the core genome) was constructed using the genome comparator tool at BIGSdb. Invariant genes were removed to reduce computation load as they will not provide cladistic information . RAxML version 8.0 software was used to create a maximum likelihood core genome phylogenetic tree using the GTRGAMMA model, which combines a GTR (Generalized Time Reversible) model for the rate of substitutions between nucleotides at a site, with a Gamma distribution model for substitution-rate heterogeneity between sites). In this model four discrete rate categories are used providing an acceptable balance between speed and accuracy of analysis. The majority-rule consensus tree was generated from 200 bootstrapped replicates of the model . FigTree software v1.4.2 was used to visualise the generated Newick trees . The TempEst program was used to compute the most likely root for the tree . In a larger scale analysis, the core genome MAFFT alignment (90% core threshold, 2182 genes) of 1053 N. gonorrhoeae isolates including the 59 WA and the 72 international N. gonorrhoeae isolates was used to construct a neighbor joining phylogeny with 100 bootstrapped replicates of the model using MEGA7: Molecular Evolutionary Genetics Analysis version 7.0 for bigger datasets . These 1053 N. gonorrhoeae isolates were selected using maximum variation sampling covering all available MLST STs, locations and years of isolates available at PubMLST database.
Sequence typing, antimicrobial genetic markers and genomic islands
Using the Short Read Sequence Typing for Bacterial Pathogens (SRST2) software [30, 31]. MLST was performed by comparing the assembled sequences of the seven housekeeping loci (abcZ, adk, aroE, fumC, gdh, pdhC and pgm) to the reference MLST profiles on the PubMLST database (http://pubmlst.org/neisseria/) . Pileup and scores files generated by SRST2 were used for manual curation. All novel allelic combinations were referred to the PubMLST website curator to be assigned a ST.
NG-MAST was performed on the sequences produced by whole genome sequencing using the genome comparator tool at BIGSdb ((http://pubmlst.org/neisseria/) . The sequences were trimmed, as described in the NG-MAST website (http://www.ng-mast.net/) and submitted to the NG-MAST database for ST determination. An integer was assigned to porB (encoding major outer membrane protein porin) and tbpB (encoding transferrin binding protein B). All novel alleles were referred to the NG-MAST website curator to be assigned an allelic number and ST.
NG-STAR, a novel molecular antimicrobial resistance typing scheme based on the sequence of seven genes associated with antimicrobial resistance in N. gonorrhoeae (penA, mtrR, porB, ponA, gyrA, parC and 23S rRNA) , was performed by comparing each allele to the publicly accessible database at https://ngstar.canada.ca by SRST2. All novel alleles and new allele combinations were referred to the NG-STAR curator to be assigned an allelic number and ST.
Plasmids were assembled using SPAdes as single contigs then compared to their corresponding references (β-lactamase producing pJD4 and the gonococcal cryptic pJD1). The variable regions between β-lactamase plasmids were aligned and were assigned based upon the method of Trembizki et al. . The three conjugative plasmids found in N. gonorrhoeae, a 39 kb plasmid (pEP5233) (GenBank accession number GU479465.1) and two 42 kb tetM-positive, the Dutch (pEP5289) and the American (pEP5050) plasmids (GenBank accession numbers GU479466.1 and GU479464.1)  were also identified as single contigs. Prior to plasmid alignment the Cyclic DNA Sequence Aligner software was used to find the optimal rotation for the circular plasmid sequences .
The genome comparator tool at PubMLST was used to detect the existence of defined loci of conjugative plasmids and gonococcal genetic islands (GGI) in http://pubmlst.org/neisseria/. The AMR genes 23 S rRNA, ponA, penA, porB, gyrA, parC, mtrR, promtrR, bla-TEM, tetM and ftsX were detected using the genome comparator tool and confirmed using SRST2. Novel alleles were aligned and notionally translated to amino acid sequences to enable detection of amino acid substitutions.
Population structure analysis
Hierarchical clustering of sequences was inferred by using the hierBAPS tandem command line program implemented in BAPS v6.0 to estimate the population structure . Concatenated core genome MAFFT alignments for all 131 N. gonorrhoeae isolates in the dataset indicating variable loci that were present in 100% of the isolates in the dataset in FASTA format were used for hierarchical clustering. Clustering was performed with two levels in the hierarchy using k = 40 as the prior upper bound for the number of clusters. Deeper levels of clustering was performed based on the first result using k = 3 to k = 7 K value and the best k value in both analyses was six. For the detection and representation of recombination between populations, the hierBAPS output file was converted to BAPS format and the ‘admixture with pre-defined populations’ approach was used in the BAPS software .
Categorical variables were examined using the Fisher’s Exact test. GraphPad Prism 7 (GraphPad Software Inc., California) was used to perform the analyses. A 5% level of confidence was used and statistical significance was determined with a p value of < 0.05. Wallace coefficients measure the extent of congruity between different metrics covering the same data. Adjusted Wallace coefficients (AW)  were used to determine the congruence of the three typing methods using the online tool available at http://www.comparingpartitions.info/index.php?link=Tool.
Molecular epidemiological typing indicates novel sequence types in remote areas
The core genome phylogeny revealed three persistent and stable genetic clusters, clusters A, B and C, which included 37 of the 59 isolates (Fig. 1). Clusters A and B were phylogenetically related and contained 28 isolates of which 78% were from remote WA (n = 22, Fig. 1). Cluster C contained 9 isolates of which 67% were from remote WA (n = 6). Approximately three quarters (17/22) of the non-A, B, C cluster isolates were collected from metropolitan areas. A significant association (p = 0.0001) of the remote isolates with the three clusters A, B and C was observed using Fisher’s exact test.
MLST identified a total of 23 STs, of which eight were new to the international PubMLST database (ST12039–12046), and these eight STs accounted for 51% of isolates (n = 30/59) of the collection. The MLST STs found in clusters A (ST12045 and ST12046) and B (ST7363, ST12040 and ST12044) were closely related with isolates sharing five of the seven alleles found in ST7363. Cluster C contained two MLST STs (ST12042 and ST12043) which shared six of the seven alleles. The 22 isolates that did not group into clusters A, B or C shared 17 different MLST STs of which two had not been described before (ST12039, ST12041) [Additional file 2].
NG-MAST typing identified 32 STs. Three quarters of the NG-MAST STs (n = 21/32) were singletons. NG-MAST STs with the most isolates were ST758 and ST9716. Both STs were cluster specific, with ST758 making up eight of 13 isolates identified in cluster A and ST9716 making up eight of the nine isolates identified in cluster C [Additional file 2]. The majority of single NG-MAST STs (17/21) were not present in the three clusters.
NG-STAR identified 26 STs. Only seven of the NG-STAR STs have previously been described and these were ST73, ST85, ST90, ST139, ST177, ST178, and ST231. Of the nineteen NG-STAR STs, four had new combinations (NCs) of previously reported alleles and 15 were novel having at least one novel locus (NVs). Table 1 defines all the NC and NV NG-STAR STs identified in this study. NG-STAR NCs STs were not found in clusters A, B and C. The 15 novel NG-STAR STs possessed novel mtrR alleles. Three novel profiles have been designated NG-STAR ST715, ST754, and ST755 which have not been described elsewhere. NG-STAR ST755 was only found in cluster A and accounted for 61% (n = 8/13) of isolates in this cluster. NG-STAR ST754 characterised three metropolitan isolates in cluster B. In addition, NVs ST11 and ST12 had a novel gyrA allele, ST1 had a novel penA allele, and ST15 possessed a novel porB allele. All isolates in clusters A, B and C had novel NG-STAR STs [Additional file 2].
Population structure analysis identifies two unique cluster groups in Western Australia
Since there were many novel STs in the West Australian dataset, core phylogeny of these isolates was compared to 72 isolates in an international dataset of strains characterized by Ezewudo et al. . The total dataset of 131 N. gonorrhoeae isolates formed six genetic groups after structure grouping by heirBAPs (Fig. 2). International isolates appeared in four genetic groups Int1-red, Int2-green, Int3-yellow and Int4-gray. While Ezewudo et al.  [Additional file 3] originally identified five structure groups, Groups 3 and 5 from that study which were noted to be highly similar collapsed into one group, Int-1, in this study. Two new groups, Australian Group 1 (Aus1-blue) and Australian Group 2 (Aus2-cyan), were formed of isolates only from the WA dataset.
Int1 was the largest structure group containing 53 isolates representing 35 different MLST STs, distributed among other structure groups in the phylogeny. In this structure group, 16 isolates were from WA of which 68% (n = 11) were metropolitan isolates and 62% (n = 10) were antibiotic resistant isolates. Int2 contained 15 isolates from the international dataset and 3 isolates from WA dataset with the majority of these being representatives of the MLST ST1901 (n = 16). In addition all the WA isolates in Int2 were also NG-STAR ST90 which was the most common ST in the NG-STAR study of Demczuk et al. . Int3 contained 12 isolates in total of which 3 isolates were from metropolitan Perth. The WA isolates in Int3 belonged to NG-STAR ST178 but were closely related to MLST STs that were not associated with AMR. All WA isolates in Int2 and Int3 were isolated from the metropolitan area and contained no novel MLST STs, NG-MAST or NG-STAR STs. Aus1 contained the 28 isolates of clusters A and B and Aus2 contained the nine cluster C isolates. The three clusters also contained most of the isolates with novel MLST, NG-MAST and NG-STAR STs.
The BAPS admixture analysis (Fig. 2b) shows Int1 is characterized by significant recombination and appears to be a nexus for gene exchange with isolates scattered across the phylogenetic tree . Three WA isolates (ExNg242, ExNg248 and ExNg316) from Int1 showed no recombination with other structure groups and were identical by MLST (ST7359), NG-MAST (ST4186) and NG-STAR (ST231) suggesting clonal expansion within the recombinant group. In contrast, isolates from the other structure groups formed tight clusters. Clonal expansion was observed in five hierBAPS groups, most significantly in Aus2, as illustrated by the core genome phylogeny and BAPS admixture analysis. Based on adjusted Wallace coefficients (AW), the isolates grouped by MLST have a 97% chance of grouping by hierBAPS, 53% by NG-MAST and 42% by NG-STAR.
To further validate the clusters identified by heirBAPs using the 131 strain dataset, a Neighbor joining tree was constructed using the core genome of 1053 N. gonorrhoeae isolates sourced from PubMLST. Australian structure groups (Aus1 and Aus2) still formed two unique clusters with Aus1 being supplemented with 11 isolates from Queensland (Australia) (Fig. 3).
Australian structure groups are characterized by an absence of antibiotic resistance and the presence of the gonococcal genetic island
The association between the genetic groups, AMR and selected genomic features of the 59 WA isolates were analyzed. The gonococcal genetic island (GGI) was identified in 86% (50/59) of isolates (Fig. 4). Conjugative plasmids, which often carry the tetM determinant, were present in 42% (24/59) of the isolates. The markerless 39 kb plasmid (pEP5233) was found in two Int1, three Int3 and three Aus1 isolates. pEP5289 (Dutch), which carries tetM, was found in five Int1 tetracycline HLR isolates. pEP5050 (American) was found in all Aus2 isolates and two Int1 isolates, and apart from ExNg314 all isolates were tetracycline HLR. Two previously identified beta-lactamase plasmid types, African (pJD5) and Rio/Torino plasmid (pJD7) were found in seven Int1 and two Int3 penicillinase-producing N. gonorrhoeae (PPNG) isolates [Additional file 2].
Other than tetracycline resistance in Aus2, all Aus1 and Aus2 isolates did not show resistance to beta-lactams, quinolones or macrolides (Fig. 4) and lacked previously reported mutations associated with these antibiotics (Additional file 2). However, the Aus1 and Aus2 isolates did possess an A39T mutation in the coding sequence of the repressor gene (mtrR, NEIS1635) which has been shown to result in the over-expression of the MtrCDE efflux pump system but not with wild type mtrR promoter .
Twenty-two WA isolates clustered in Int1, Int2 and Int3 structure groups. Approximately two thirds (62%, 16/22) were resistant to one or more antibiotics. A penA type 63 mosaic allele (NEIS1753 allele 1147) was identified in one susceptible strain from Int1. Int2 isolates contained NG-STAR ST90 that has been reported associated with ESCs and ciprofloxacin resistance. Decreased susceptibility to the ESCs due to the penA (NEIS1753) motif XXXIV was restricted to Int2. Int2 isolates also contained GGI and (pEP5233) conjugative plasmids. Int3 did not contain any chromosomal mutations affecting beta-lactamase, quinolone and macrolide susceptibility consistent with the antibiotic sensitive phenotype. mtrR with an internal stop codon but no A39T or G45D mutations was found in all Int3 isolates [Additional file 2].
All isolates with high-level resistance to ciprofloxacin (QRNG HLR) had mutations in gyrA (NEIS1320) and parC, (NEIS1525) conferring resistance to quinolones. High-level resistance to azithromycin (AziRNG HLR) mediated by the A2059G mutation in genes encoding for 23S rRNA was identified in AziRNG HLR isolate ExNg321. A recently reported R251H substitution in the ftsX gene (NEIS2146 allele 18) which confers ESCs resistance [40, 41], was found in five extra-genital N. gonorrhoeae isolates however only two of these isolates demonstrated decreased susceptibility to ceftriaxone.
The high rate of recombination within the N. gonorrhoeae population has led to the presumption that this species has a non-clonal population structure . However, clonal spread and persistence of N. gonorrhoeae strains has been described in many whole genome sequence-based epidemiological studies from Europe, France, USA, Canada, and recently from a global collection of N. gonorrhoeae [11, 12, 14, 15, 43]. Population structure analysis of N. gonorrhoeae isolates in the study dataset has provided some evidence for both recombination and clonality in N. gonorrhoeae, suggesting N. gonorrhoeae is a weakly clonal species similar to H. pylori and Bacillus sphaericus [12, 44, 45]. However, an analysis of an expanded dataset over a prolonged period of time would provide a much clearer understanding of the N. gonorrhoeae population structure.
In the remote areas of WA, when compared to the global population, gonococcal isolates formed two distinctive structure groups (Aus1 and Aus2). Aus1 was characterized by broad antibiotic-susceptibility and was mostly composed of isolates related to novel MLST, NG-MAST and NG-STAR STs. The novelty of the MLST, NG-MAST ST and NG-STAR designations is significant as the isolates were compared to over 3700 isolates found in the PubMLST database, which represents collections from around the world, particularly Europe and the US. All isolates in Aus1 belonged to two core genome clusters, A and B. Isolates belonging to MLST ST7367, which have been reported from global sources, are associated with decreased susceptibility to ceftriaxone and have been shown to cluster with Int1 by Ezewudo et al. . However, the MLST ST7363 isolates of Aus1 (cluster A) formed a discrete cluster group apart from Int1, suggesting they are unique. In cluster A of Aus1, MLST ST12045 isolates were associated with a previously reported NG-MAST ST758 reported from Russia which is associated with decreased susceptibility to ceftriaxone . However, the Aus1 MLST ST12045 isolates were susceptible to ceftriaxone and possessed novel NG-STAR profiles. Collectively this suggests Aus1 is a distinctive population of isolates circulating in remote Australian regions. Aus2 represents a local dissemination of a highly clonal tetracycline HLR population which was characterized by novel MLST, NG-MAST and NG-STAR assignments. The presence of pEP5050 (American) conjugative plasmid in Aus2 suggests a global introduction of the initial strain to the region followed by clonal persistence of this group. One Aus2 isolate that was not tetracycline HLR, ExNg314, possessed an intact conjugative plasmid carrying the tetM determinant but did not have the V57 M in rspJ gene, unlike the other tetracycline HLR isolates in this group . This suggests the V57 M mutation in combination with the plasmid mediated resistance determinant is required to achieve high-level tetracycline resistance. The absence of other mutations in the AMR genetic markers in the Aus groups suggests no independent emergence of resistance (de novo) since divergence from its ancestor. Lastly Aus1 and Aus2 were characterized by possession of porB1a (68%, 27/40) which is associated with asymptomatic urethral infections and disseminated gonococcal infection (DGI) [48, 49]. Though porB1a is not associated with MDR strains, monitoring these isolates is crucial as DGI is considered a common cause of complicated gonococcal infection and infertility [49, 50].
Apart from Aus1 and Aus2, 22/59 isolates from WA in this period clustered with the Int1, Int2 and Int3 structure groups. WA isolates from Int1 were retrieved primarily from metropolitan areas. The strains are genetically unrelated and highly recombinant, and suggest the sporadic emergence of various N. gonorrhoeae clusters or represent sporadic incursions of international strains into WA. The ciprofloxacin resistant isolate ExNg232 which had a NG-MAST 225 and NG-STAR177 STs, was previously reported in Europe and the US [11, 12], suggesting an introduction event to the metropolitan areas of WA. In another instance, a high-level azithromycin resistant isolate (ExNg321) had a NG-MAST ST4850 which differed by one SNP from a ST1866 AziRNG HLR isolate reported in China and a MLST ST12039 that was a single locus variant of a ST10899 AziRNG HLR isolate reported in Canada [24, 51]. This would suggest that ExNg321 is a foreign introduction into WA and may be one of the earliest isolates to cause a recent outbreak of azithromycin resistant strains in South Australia that led to an updated treatment recommendation to avoid single agent azithromycin treatment . The antibiotic sensitive isolates of Int1 from MLST7359 have recently been linked to increasing rates of gonococcal infections among women in heterosexual networks in metropolitan areas of New South Wales in Australia .
Three isolates that clustered with Int2 represented the global AMR cluster MLST ST1901 and the two related STs (matching at six MLST loci) that are associated with decreased susceptibility and resistance to ceftriaxone [11, 12, 41]. One of the isolates (ExNg304) belonged to NG-STAR ST90, which has been reported to be associated with decreased susceptibility and resistance to ceftriaxone and ciprofloxacin. However, although ExNg304 had all the requisite mutations for this phenotype, the isolate was phenotypically susceptible. Lastly, all the WA isolates in Int3 were NG-STAR ST178, a cluster that is antibiotic susceptible and does not contain any MLST or NG-MAST STs associated with resistance . This appears to represent a potential persistent susceptible cluster that has been circulating globally since 1998 .
The gonococcal genetic island (GGI) reported by Harrison et al.  has been significantly associated with different MDR core genome clusters such as clusters of MLST ST1901 and ST1508. In the present study, GGI was identified in most isolates (50/59, 84%) in WA N. gonorrhoeae populations, being found in 68% (13/19) of isolates in MDR groups Int1 and Int2 and 95% (38/40) of isolates in the non-MDR groups Aus1, Aus2 and Int3. This suggests that the presence of GGI in a given structure group may provide a fitness advantage and contribute to the stability of that cluster, which occurred before the acquisition of the MDR loci by natural transformation or mutation.
The population structure of a small set of gonococcal isolates from WA has revealed the presence of two unique clusters: Aus1 and Aus2. Due to the small sample size and the lack of representative isolates from Australia and SouthEast Asia in the PubMLST database, it remains to be determined if these two clusters are more widely dispersed in this geographic region. However, the clustering of 11 gonococcal isolates from Queensland into Aus1 is suggestive of a broader distribution across Australia which will be confirmed as much larger studies in these regions are underway. The persistence of the antibiotic susceptible cluster Aus1 in WA could be due to azithromycin and amoxicillin dual therapy in remote WA, thereby reducing the selective pressure for AMR or could be the result of the remoteness of this region which has impeded incursion of AMR strains . The NG-STAR typing scheme correlates well with core genome phylogeny of N. gonorrhoeae and is sufficient for high throughput surveillance of AMR. However with the recent report that Bexsero® vaccination may reduce the rate of gonorrhoea , WGS is needed for antigenic profiling and monitoring the persistence of antigenic combinations over short timeframes in response to vaccination in clinical trials and given its high level of discrimination is the most accurate way of defining sexual networks.
Australian group 2
Azithromycin-resistant N. gonorrhoeae
Bacterial Isolates Genome Sequence database
Disseminated gonococcal infection
Gonococcal Antimicrobial Surveillance Programme
Gonococcal genetic island
Minimum inhibitory concentration
Multilocus sequence typing
Men who have sex with men
N. gonorrhoeae multi antigen sequence typing
N. gonorrhoeae sequence typing for antimicrobial resistance
Penicillinase-producing N. gonorrhoeae
Quinolone-resistant N. gonorrhoeae
Western Australian Gonococcal Surveillance Programme
Whole genome sequencing
Unemo M, del Rio C, Shafer WM. Antimicrobial resistance expressed by Neisseria gonorrhoeae: a major global public health problem in the 21st century. Microbiol Spec. 2016;4(3):1-32.
Roberts-Witteveen A, Pennington K, Higgins N, Lang C, Lahra M, Waddell R, Kaldor J. Epidemiology of gonorrhoea notifications in Australia, 2007-12. Sex Health. 2014;11(4):324–31.
Whiley DM, Trembizki E, Buckley C, Freeman K, Baird RW, Beaman M, Chen M, Donovan B, Kundu RL, Fairley CK et al. Molecular antimicrobial resistance surveillance for Neisseria gonorrhoeae, Northern Territory, Australia. Emerg Infect Dis. 2017;23(9):1478-85.
Mitchell K, Giele C, Byron M. The epidemiology of notifiable sexually transmitted infections and blood-borne viruses in Western Australia 2015. Western Australia: Department of Health; 2016.
WHO guidelines for the treatment of Neisseria gonorrhoeae [http://www.who.int/reproductivehealth/publications/rtis/gonorrhoea-treatment-guidelines/en/]. Accessed 15 Jun 2017.
Wi T, Lahra MM, Ndowa F, Bala M, Dillon J-AR, Ramon-Pardo P, Eremin SR, Bolan G, Unemo M. Antimicrobial resistance in Neisseria gonorrhoeae: global surveillance and a call for international collaborative action. PLoS Med. 2017;14(7):e1002344.
Goire N, Lahra MM, Chen M, Donovan B, Fairley CK, Guy R, Kaldor J, Regan D, Ward J, Nissen MD et al. Molecular approaches to enhance surveillance of gonococcal antimicrobial resistance. Nat Rev Microbiol. 2014;12(3):223–9.
Demczuk W, Sidhu S, Unemo M, Whiley DM, Allen VG, Dillon JR, Cole M, Seah C, Trembizki E, Trees DL, et al. Neisseria gonorrhoeae sequence typing for antimicrobial resistance (NG-STAR): a novel antimicrobial resistance multilocus typing scheme for tracking the global dissemination of N. gonorrhoeae strains. J Clin Microbiol. 2017;55(5):1454-68.
Whiley DM, Goire N, Rahimi F, Lahra MM, Limnios AE, Nissen MD, Sloots TP. Real-time PCR genotyping of Neisseria gonorrhoeae isolates using 14 informative single nucleotide polymorphisms on gonococcal housekeeping genes. J Antimicrob Chemother. 2013;68(2):322–8.
Deguchi T, Yasuda M, Hatazaki K, Kameyama K, Horie K, Kato T, Mizutani K, Seike K, Tsuchiya T, Yokoi S, et al. New clinical strain of Neisseria gonorrhoeae with decreased susceptibility to ceftriaxone, Japan. Emerg Infect Dis. 2016;22(1):142.
Chisholm S, Unemo M, Quaye N, Johansson E, Cole M, Ison C, Van de Laar MJ. Molecular epidemiological typing within the European Gonococcal Antimicrobial Resistance Surveillance Programme reveals predominance of a multidrug-resistant clone. Euro Surveill. 2013;18(3):1-10.
Grad YH, Harris SR, Kirkcaldy RD, Green AG, Marks DS, Bentley SD, Trees D, Lipsitch M. Genomic epidemiology of gonococcal resistance to extended-spectrum cephalosporins, macrolides, and fluoroquinolones in the United States, 2000–2013. J Infect Dis. 2016;214(10):1579–87.
Papp JR, Abrams AJ, Nash E, Katz AR, Kirkcaldy RD, O’Connor NP, O’Brien PS, Harauchi DH, Maningas EV, Soge OO. Azithromycin resistance and decreased ceftriaxone susceptibility in Neisseria gonorrhoeae, Hawaii, USA. Emerg Infect Dis. 2017;23(5):830-32.
Demczuk W, Martin I, Peterson S, Bharat A, Van Domselaar G, Graham M, Lefebvre B, Allen V, Hoang L, Tyrrell G et al. Genomic epidemiology and molecular resistance mechanisms of azithromycin resistant Neisseria gonorrhoeae in Canada from 1997 to 2014. J Clin Microbiol. 2016;54(5):1304-13.
Belkacem A, Jacquier H, Goubard A, Mougari F, La Ruche G, Patey O, Micaëlo M, Semaille C, Cambau E, Bercot B. Molecular epidemiology and mechanisms of resistance of azithromycin-resistant Neisseria gonorrhoeae isolated in France during 2013–14. J Antimicrob Chemother. 2016;71(9):2471-8.
Lahra MM and Enriquez RP. Australian Gonococcal Surveillance Programme annual report, 2015. Commun Dis Intell. 2017;41(1):E60–7.
Lahra MM. Australian Gonococcal Surveillance Programme annual report, 2014. Commun Dis Intell Q Rep. 2015;39(3):E347–54.
Trembizki E, Wand H, Donovan B, Chen M, Fairley CK, Freeman K, Guy R, Kaldor JM, Lahra MM, Lawrence A, et al. The molecular epidemiology and antimicrobial resistance of Neisseria gonorrhoeae in Australia: a Nationwide cross-sectional study, 2012. Clin Infect Dis. 2016;63(12):1591–8.
O’Reilly LC, Goire N, Fisk RE, Speers DJ. Molecular epidemiology of Neisseria gonorrhoeae using multi-antigen sequence typing and pulse-field gel electrophoresis in highly endemic Western Australian populations. BMC Infect Dis. 2015;15(1):272.
Bell S, Pham J, Newton P, Nguyen T. Antibiotic susceptibility testing by the CDS method. A manual for Medical and Veterinary Laboratories 2013. 2013, Seventh Edition.
Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, Lesin VM, Nikolenko SI, Pham S, Prjibelski AD, et al. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol. 2012;19(5):455–77.
Gurevich A, Saveliev V, Vyahhi N, Tesler G. QUAST: quality assessment tool for genome assemblies. Bioinformatics. 2013;29(8):1072–5.
Jolley KA, Maiden MCJ. BIGSdb: scalable analysis of bacterial genome variation at the population level. BMC Bioinformatics. 2010;11(1):595.
Ezewudo MN, Joseph SJ, Castillo-Ramirez S, Dean D, Del Rio C, Didelot X, Dillon J-A, Selden RF, Shafer WM, Turingan RS, et al. Population structure of Neisseria gonorrhoeae based on whole genome data and its relationship with antibiotic resistance. PeerJ. 2015;3:e806.
Wise MJ. dCITE: Measuring necessary cladistic information can help you reduce polytomy artefacts in trees. PLoS One. 2016;11(11):e0166991.
Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014;30(9):1312–3.
FigTree - Produce high-quality figures of phylogenetic trees [http://tree.bio.ed.ac.uk/software/figtree/]. Accessed 23 Oct 2016.
Rambaut A, Lam TT, Max Carvalho L, Pybus OG. Exploring the temporal structure of heterochronous sequences using TempEst (formerly path-O-gen). Virus Evol. 2016;2(1):vew007.
Kumar S, Stecher G, Tamura K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33(7):1870–4.
Inouye M, Dashnow H, Raven L-A, Schultz MB, Pope BJ, Tomita T, Zobel J, Holt KE. SRST2: rapid genomic surveillance for public health and hospital microbiology labs. Genome Medicine. 2014;6(11):90.
Inouye M, Conway TC, Zobel J, Holt KE. Short read sequence typing (SRST): multi-locus sequence types from short reads. BMC Genomics. 2012;13:338.
Martin IM, Ison CA, Aanensen DM, Fenton KA, Spratt BG. Rapid sequence-based identification of gonococcal transmission clusters in a large metropolitan area. J Infect Dis. 2004;189(8):1497–505.
Trembizki E, Buckley C, Lawrence A, Lahra M, Whiley D. Characterization of a novel Neisseria gonorrhoeae penicillinase-producing plasmid isolated in Australia in 2012. Antimicrob Agents Chemother. 2014;58(8):4984–5.
Pachulec E, van der Does C. Conjugative plasmids of Neisseria gonorrhoeae. PLoS One. 2010;5(4):e9962.
Fernandes F, Pereira L, Freitas AT. CSA: an efficient algorithm to improve circular DNA multiple alignment. BMC Bioinformatics. 2009;10(1):230.
Cheng L, Connor TR, Sirén J, Aanensen DM, Corander J. Hierarchical and spatially explicit clustering of DNA sequences with BAPS software. Mol Biol Evol. 2013;30(5):1224–8.
Corander J, Marttinen P, Sirén J, Tang J. Enhanced Bayesian modelling in BAPS software for learning genetic structures of populations. BMC Bioinformatics. 2008;9(1):539.
Severiano A, Pinto FR, Ramirez M, Carriço JA. Adjusted Wallace coefficient as a measure of congruence between typing methods. J Clin Microbiol. 2011;49(11):3997–4000.
Thakur S, Starnino S, Horsman G, Levett P, Dillon J. Unique combined penA/mtrR/porB mutations and NG-MAST strain types associated with ceftriaxone and cefixime MIC increases in a ‘susceptible’Neisseria gonorrhoeae population. J Antimicrob Chemother. 2014;69(6):1510–6.
Gong Z, Lai W, Liu M, Hua Z, Sun Y, Xu Q, Xia Y, Zhao Y, Xie X. Novel genes related to ceftriaxone resistance found among ceftriaxone-resistant Neisseria gonorrhoeae strains selected in vitro. Antimicrob Agents Chemother. 2016;60(4):2043–51.
de Curraize C, Kumanski S, Micaëlo M, Fournet N, La Ruche G, Meunier F, Amarsy-Guerle R, Jacquier H, Cambau E, Goubard A. Ceftriaxone-resistant Neisseria gonorrhoeae isolates (2010–2014) in France characterized using whole genome sequencing. Antimicrob Agents Chemother. 2016;60(11):6962-4.
O'Rourke M, Stevens E. Genetic structure of Neisseria gonorrhoeae populations: a non-clonal pathogen. Microbiology. 1993;139(11):2603–11.
Vidovic S, Caron C, Taheri A, Thakur SD, Read TD, Kusalik A, Dillon J-AR. Using crude whole-genome assemblies of Neisseria gonorrhoeae as a platform for strain analysis: clonal spread of gonorrhea infection in Saskatchewan, Canada. J Clin Microbiol. 2014;52(10):3772–6.
Achtman M, Azuma T, Berg DE, Ito Y, Morelli G, Pan ZJ, Suerbaum S, Thompson SA, Van Der Ende A, Van Doorn LJ. Recombination and clonal groupings within Helicobacter pylori from different geographical regions. Mol Microbiol. 1999;32(3):459–70.
Ge Y, Hu X, Zheng D, Wu Y, Yuan Z. Allelic diversity and population structure of Bacillus sphaericus revealed by multilocus sequence typing. Appl Environ Microbiol. 2011;77(15):5553-6.
Unemo M, Vorobieva V, Firsova N, Ababkova T, Leniv I, Haldorsen B, Fredlund H, Skogen V. Neisseria gonorrhoeae population in Arkhangelsk, Russia: phenotypic and genotypic heterogeneity. Clin Microbiol Infect. 2007;13(9):873–8.
Hu M, Nandi S, Davies C, Nicholas RA. High-level chromosomally mediated tetracycline resistance in Neisseria gonorrhoeae results from a point mutation in the rpsJ gene encoding ribosomal protein S10 in combination with the mtrR and penB resistance determinants. Antimicrob Agents Chemother. 2005;49(10):4327–34.
Rechner C, Kühlewein C, Müller A, Schild H, Rudel T. Host glycoprotein Gp96 and scavenger receptor SREC interact with PorB of disseminating Neisseria gonorrhoeae in an epithelial invasion pathway. Cell Host Microbe. 2007;2(6):393–403.
O'Brien JP, Goldenberg DL, Rice PA. Disseminated gonococcal infection: a prospective analysis of 49 patients and a review of pathophysiology and immune mechanisms. Medicine. 1983;62(6):395–406.
Levens E. Disseminated gonococcal infection. Primary care update for OB/GYNS. 2003;10(5):217–9.
Yuan L-F, Yin Y-P, Dai X-Q, Pearline RV, Xiang Z, Unemo M, Chen X-S. Resistance to azithromycin of Neisseria gonorrhoeae isolates from 2 cities in China. Sex Transm Dis. 2011;38(8):764–8.
Lahra MM, Ward A, Trembizki E, Hermanson J, Clements E, Lawrence A, Whiley D. Treatment guidelines after an outbreak of azithromycin-resistant Neisseria gonorrhoeae in South Australia. Lancet Infect Dis. 2017;17(2):133–4.
Buckley C, Forde BM, Trembizki E, Lahra MM, Beatson SA, Whiley DM. Use of whole genome sequencing to investigate an increase in Neisseria gonorrhoeae infection among women in urban areas of Australia. Sci Rep. 2018;8(1):1503.
Harrison OB, Clemence M, Dillard JP, Tang CM, Trees D, Grad YH, Maiden MCJ. Genomic analyses of Neisseria gonorrhoeae reveal an association of the gonococcal genetic island with antimicrobial resistance. J Infect. 2016;73(6):578-87.
Petousis-Harris H, Paynter J, Morgan J, Saxton P, McArdle B, Goodyear-Smith F, Black S. Effectiveness of a group B outer membrane vesicle meningococcal vaccine against gonorrhoea in New Zealand: a retrospective case-control study. Lancet. 2017;390(10102):1603–10.
Collection of gonococcal isolates is facilitated by the National Neisseria Network in Australia. We thank Dr. Alfred Tay from the Marshall Center for Infectious Diseases at University of Western Australia for laboratory technical support and guidance. We also thank Dr. Jukka Corander from Department of Biostatistics at University of Oslo for his genomic population structure analysis expertise. We thank Mr. Shakeel Mowlaboccus from the Marshall Center for Infectious Diseases at University of Western Australia for the statistical analysis.
BA is supported by the King Abdullah Scholarship program (Saudi Arabia). CMK is supported by National Medical Research Council (APP1078642). The funding body had no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Availability of data and materials
All sequence types and alleles generated and analysed during the current study are available in the publically accessible Bacterial Isolates Genome Sequence database (BIGSdb) hosted on www.pubmlst.org/neisseria. All of the alleles are located in Additional file 2 and can be used to identify the relevant records at BIGSdb.
Ethics approval and consent to participate
In Australia, gonorrhoea is a notifiable infectious disease. De-identified patient samples are sent to the state-based reference laboratory, PathWest, for analysis during routine clinical care. Routine surveillance data is collected by the West Australian Department of Health from patient records and includes site of isolation, date of sample collection and postcode. The linkage of routine surveillance data with the bacterial isolate does not require patient permission since it does not result in the disclosure of data that could identify the patient. The research project was assessed by the UWA Human Research Ethics Committee (HREC) and has the approval number RA/4/20/4010.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Sequencing and assembly quality statistics of the 59 Western Australian N. gonorrhoeae isolates. (PDF 63 kb)
Table S2. WA N. gonorrhoeae isolates metadata, antimicrobial resistance profiles, sequence types and genetic characteristics of antimicrobial resistance markers. (PDF 95 kb)
About this article
Cite this article
Al Suwayyid, B.A., Coombs, G.W., Speers, D.J. et al. Genomic epidemiology and population structure of Neisseria gonorrhoeae from remote highly endemic Western Australian populations. BMC Genomics 19, 165 (2018). https://doi.org/10.1186/s12864-018-4557-5
- Neisseria gonorrhoeae
- Whole genome sequencing
- Population structure
- Western Australia