This article has Open Peer Review reports available.
Molecular epidemiology of Neisseria gonorrhoeae using multi-antigen sequence typing and pulse-field gel electrophoresis in highly endemic Western Australian populations
© O'Reilly et al. 2015
Received: 6 March 2015
Accepted: 17 June 2015
Published: 15 July 2015
The remote and indigenous populations of Western Australia (WA) have one of the highest notification rates of gonorrhoea in the world. Despite this, the low rate of antimicrobial resistance in Neisseria gonorrhoeae from these regions permits the use of amoxycillin as empirical therapy. We describe the first molecular epidemiological study of gonococci isolated from this population using two different typing platforms.
Pulse-field gel electrophoresis (PFGE), Neisseria gonorrhoeae multi-antigen sequence typing (NG-MAST) and antimicrobial susceptibility tests were performed on 128 consecutive N. gonorrhoeae isolates cultured between January 2011 and December 2013. To highlight clusters isolates were evaluated based on their tbpB sequence types.
No predominant NG-MAST or PFGE types were found. A total of 67 distinct PFGE pulsotypes were identified amongst the 128 isolates in this study with 20 PFGE pulsotypes representing 78 isolates. A total of 59 NG-MAST sequence types were found, represented by 45 porB alleles and 28 tbpB alleles with 13 tbpB genomogroups from 45 NG-MAST sequence types. TbpB genomogroup 29, represented by 45 isolates, was by far the most common genomogroup overall.
Results from this study suggest that gonococcal epidemiology in WA is quite different between remote regions and major population centres and, in some cases, geographically restricted. It is likely that isolates originating from endemic regions of WA mostly represent independent, small sexual networks with an infrequent interchange between other communities and regions. Given the high rate of antimicrobial resistance elsewhere in Australia, ongoing surveillance is essential to ensure the enduring efficacy of amoxycillin empiric use in the remote regions of WA.
Neisseria gonorrhoeae, the etiological agent of the sexually transmitted disease gonorrhoea, has emerged as a challenging pathogen due to its propensity to acquire resistance to a wide range of antimicrobial agents. With the world-wide spread of antimicrobial resistant strains, epidemiological surveillance has acquired new importance for the control of gonorrhoea . The remote and indigenous populations of Western Australia (WA), Australia’s largest state by area, have one of the highest notification rates of gonorrhoea in the world yet the circulating gonococci in these populations harbour very low levels of penicillin resistance . As a result, amoxycillin combined with probenecid remains successful for the empiric management of gonorrhoeae acquired locally in these remote regions . Moreover, to ensure the ongoing success of this treatment strategy surveillance for strains associated with antimicrobial resistance is required. This is especially so given the ongoing risk of resistant strain incursions from the more heavily populated metropolitan areas of WA and from the influx of inter-state and overseas workers for resource-mining related activities in these regions.
Most bacterial genotyping methods such as multi-locus sequence typing and pulse field-gel electrophoresis (PFGE) require live cultures. However, much of the gonococcal disease notification from remote region WA is based on molecular detection due to the difficulties in obtaining viable gonococci from remote regions. In this study N. gonorrhoeae multi-antigen sequence typing (NG-MAST) was compared to PFGE as NG-MAST has been successfully implemented on non-cultured clinical specimens . Both methods have been used to identify transmission clusters to aid medico-legal investigations, but there is a lack of data comparing the performance of these methods for gonococcal typing [5–11]. In addition, tbpB-sequence based analysis was also performed as recent studies suggest that of the two loci used in NG-MAST study, tbpB exhibits lower variability and therefore can show distribution patterns of isolates. We examined the utility of tbpB sequence-based classification in determining disease transmission networks in these endemic populations [12, 13]. By employing PFGE, NG-MAST and tbpB-sequence based analysis, this study aimed to ascertain the molecular epidemiology of circulating gonococci from remote regional WA over the 3-year period from 2011 to 2013.
All N. gonorrhoeae isolates were cultured using standard microbiology laboratory protocols. The identification of bacterial colonies resembling N. gonorrhoeae was performed using at least two of the following identification platforms, Vitek2 NHI cards (Biomerieux), Phadebact gonococcal agglutination tests and Matrix-assisted laser desorption-ionisation test-of-flight (Bruker Deltonics, Bremen, Germany) following manufacturers’ instructions. The identity of each isolate was confirmed using a previously described porA gene specific N. gonorrhoeae real time PCR assay . Antimicrobial susceptibility testing was performed on all isolates using the Etest method for penicillin and ceftriaxone following Central Laboratory Standards Institute guidelines and beta-lactamase activity was assayed using nitrocefin discs (Oxoid, Australia).
NG-MAST was performed as previously described  with modification. Briefly, each reaction mix consisted of 12 μL of PCR mix containing 0.2 mM of each forward and reverse primer (porB-Forward- CAAGAAGACCTCGGCAA and porB-Reverse- CCGACAACCACTTGGT for porB reaction and tbpB-Forward- CGTTGTCGGCAGCGCGAAAAC and tbpB-Reverse- TTCATCGGTGCGCTCGCCTTG for tbpB reaction), 0.5 units of DNA Polymerase (Applied Biosystems, USA), 2 mM Magnesium chloride (Invitrogen, USA), 1 unit of 10x PCR buffer (Applied Biosystems, USA), 0.01 % Bovine Serum Albumin (Sigma Aldrich, USA), 0.2 μM dNTP pool (Fisher Biotec, Australia), and 1 μM of 5-Carboxy-X-rhodamine succinimidyl ester (ROX) SYTO9 dye (Thermofisher Scientific, Australia). DNA lysate prepared by heat-denaturing gonococcal colonies (heating to 100 °C for 15 min) with 8 μL added to the PCR mix for a final reaction volume of 20 μL. Thermocycling was performed on a RotorGene 6000 real-time PCR instrument (Corbett Life Science, Australia) with the following cycling conditions; 95 °C for 15 min, followed by 45 cycles of 95 °C for 30 s, 60 °C for 30 s, and 72 °C for 1 min. PCR product formation was indicated by the presence of real-time PCR amplification curves. The amplified porB and tbpB sequences were then submitted for DNA sequencing using the ABI Prism BigDye Terminator v3.1 system (Applied Biosystems, Foster City, CA) according to the manufacturer’s instructions and the sequence data analysed using an ABI Prism 3130XL 16-channel Genetic Analyzer. The porB and tbpB sequences were trimmed to 490 bp and 390 bp respectively starting at particular sites of sequence conservation (TTGAA at the pre-loop 3 to loop 6 for porB and CGTCTGAA for tbpB) and then submitted to the NG-MAST website (http://www.ng-mast.net/). For all isolates returning previously unidentified (new) NG-MAST sequence types, the respective porB and tbpB sequences were submitted for addition to the NG-MAST database and assignment of a new sequence type (ST).
Samples were prepared using the standard PFGE method (PulseNet USA protocol, Centers for Disease Control and Prevention, Atlanta, Georgia. USA. 2004). Digestion with Nhel (Promega, Madison, USA) was performed overnight at 37 °C in 300 μL of fresh buffer containing 30 units of restriction endonuclease. The digested plugs were sealed into slots in 1 % agarose gel (Seakem Gold Agarose, Lonza, Rockland USA) and subjected to electrophoresis in a contour-clamped homogeneous electric field apparatus (BioRad, Hercules, California, USA). The electrophoresis was performed with pulse times ramping from 2.2 to 35 s for 18.3 h. Gels were stained with ethidium bromide. The Nhel enzyme produced approximately 15 to 20 distinct DNA fragments. Strains were considered identical if no fragment differences occurred. Salmonella enteritidis ser Braenderup H9812 was used as the reference ladder for normalisation. The results were analyzed using the BioNumerics software package version 6.6 (AppliedMaths, Kortrijk, Belgium) and the dendrograms were calculated using the unweighted-pair group method using the Dice coefficient. Pulsotypes (PT) were assigned numbers consecutively based on differences of more than one band in pulse-field patterns upon visual inspection.
Diversity of both STs and PT were calculated using Simpsons’ index of diversity with an online calculator (http://darwin.phyloviz.net/ComparingPartitions/index.php?link=Home) as per a previously published method . Discriminatory indices for PFGE and NGMAST were determined at 0.97 and 92.5 respectively.
No human material and only bacterial isolates originating from clinical specimens were used in this study. The use of human data in this study did not require an approval from the ethics committee of our organisation (Western Australian Health Department and the Queen Elizabeth II Medical Centre, Perth) as only subject data already obtained for diagnostic testing were used, and all data were de-identified.
Bacterial identification and susceptibility tests
All isolates were confirmed as N. gonorrhoeae by both phenotypic methods and PCR. Reduced susceptibility or resistance to ceftriaxone was not identified. A total of 17 isolates (13 %) were resistant to penicillin, with 13 identified as penicillinase-producers and 4 exhibiting chromosomally mediated resistance. As expected, penicillin resistance was identified more often in metropolitan isolates (23 %) compared to remote region isolates (8 %).
Distribution of NGMAST types, tbpB groups and their associated pulsotypes by year
Number of Isolates/year
Associated NGMAST types (n)1
tbpB Group (n)
Associated pulsotypes (n)1
10111, 10115, 10118
Group 21 (3)
5, 12, 23
Group 25 (2)
758 (25), 1498, 3042, 4822, 5961, 8063 (4) 7206 (3), 7353, 7707 (2), 10102, 10107, 10115 (2), 10123, 10124
Group 29 (45)
1, 2(2), 3(2), 5(3), 8(4), 12 (10), 13, 15 (3), 13, 16, 17, 19, 20, 22, 24, 34(2), 42, 49, 54, 65, 66, 68, 73(2), 76, 77
21, 5267, 8842, 5785, 10108
Group 33 (5)
29, 30, 46, 51, 73
Group 110 (2)
Group 197 (2)
4168 (3), 10121
Group 241 (4)
28(2), 47, 58
7126 (12), 7803 (2)
Group 349 (14)
4(6), 10(2), 21, 39(3), 50, 73
10080(4), 10103 (3), 10104, 10110, 10122, 10126
Group 1330 (11)
5(2), 12, 13, 14, 18, 73(5)
Group 1344 (2)
10116 (2), 10117 (2)
Group 1408 (4)
25, 48, 72(2)
Group 1481 (3)
9716 (15), 10113
Group 1749 (16)
8, 9(5), 12, 43, 44, 60, 71(5), 74
Classification based on tbpB sequences
There were 14 porB and 13 tbpB alleles that were associated with more than one isolate, accounting for 98 (76 %) and 113 (88 %) of all isolates respectively. The 13 tbpB-groups (Table 1) represented 45 of the STs. TbpB-group 29, represented by 45 isolates, was by far the most common genomogroup overall, being the most common group in the Metropolitan (n = 8), Kimberley (n = 22) and Mid-West (n = 12) regions (Fig. 1). In contrast, tbpB-groups 349 (n = 9) and 1749 (n = 10) were the most common in the Pilbara with only 3 tbpB-group 29 isolates from this region.
We have described the first molecular epidemiological study of N. gonorrhoeae isolates from WA encompassing regions with some of the highest disease notification rates of gonorrhoea in the world. This study also represents the largest collection of N. gonorrhoeae isolates studied using the NG-MAST and PFGE methods in parallel.
The great diversity observed amongst the NG-MAST types in this study is most likely due to the hyper-variable nature of the target loci, especially the porB locus, evolving under selection pressure over years. Therefore we classified our isolates based on the relatively more conserved tbpB gene as this was more likely to demonstrate linkages over time . Although NG-MAST still remains crucial for international comparison and tracking of strains of epidemiological importance, utility of single locus based typing for identifying closely-related clusters over a short period of time has been highlighted in studies from Canada and Germany [12, 13]. Of further interest was the gradual emergence of tbpB-group 1749 in this study, which was represented by 6 different PT whereas tbpB group 39 maintained relative stability in number over this three year period (Table 1).
Transmission of sexually transmitted disease like gonorrhoea is driven by a ‘core group’ of high risk individuals with high rates of infection and frequent partner change. A second ‘bridging group’, comprised of a population with lower risks, acts as the link between the general population and the core group [13, 16]. Identifying ‘core-group’ clusters is therefore very important in gonococcal disease transmission studies and is of particular relevance in remote communities where close-knit sexual networks are likely to exist . The tbpB-group based classification identified the presence of four large clusters- namely tbpB-groups 29, 349, 1330 and 1749 which accounted for 67 % percent of all isolates in the study (Table 1). TbpB-group 29, the most numerous genomogroup in WA, could represent such a cluster in the Kimberley region where 11 out of 22 tbpB-group 29 isolates from this region originated from the same local precinct. Similarly, 8 out of the 12 tbpB-group 29 isolates from the Mid-West region were from two neighbouring precincts. In the Pilbara region, tbpB-groups 349 and 1749 together represented 19 isolates, 16 of which could be traced back to just two precincts. Of further interest, tbpB-group 29 was infrequently found amongst isolates from the Pilbara region (n = 3) despite its predominance in the surrounding Kimberley and Mid-West regions. Based on this, NGMAST has the potential to provide highly useful epidemiological information in our population, especially if applied on non-cultured clinical specimens which comprise the majority of notifications in WA.
Surprisingly, there was a major discordance between the PFGE and NG-MAST profiles of the isolates which remained unresolved despite reclassification into tbpB-groups. The constantly evolving NG-MAST loci are known to undergo changes to their porB and tbpB alleles over relatively short periods of time which would explain their distribution over a wide range of PT. Although PFGE exhibited higher discriminatory power compared to NG-MAST, based on this study, use of both methods is likely to yield more definitive results in medico-legal investigations [7, 8]. Further comparisons of the performance of these two methods over a larger population would possibly yield more conclusive results.
The results from this study are most consistent with the hypothesis that the majority of the isolates originated from independent, small sexual networks with an infrequent interchange between other communities and regions. In terms of the N. gonorrhoeae genomogroups, it appears that some genomogroups were maintained in high transmission cycles, some were disappearing whereas others were establishing efficient transmission following recent introduction. No association could be made between antimicrobial susceptibility patterns and the ST, tbpB-group or PT.
A limitation of this study was that the cultured isolates included in this study comprise only a proportion of total gonorrhoea notifications in WA since most notifications in WA are based on molecular laboratory confirmation only. According to the Australian Gonococcal Surveillance Program (AGSP) the total number of notifications of gonorrhoeae in WA in the year 2012 was 548, thus the 38 isolates from 2012 in this study represented approximately 7 % . Despite this, our study was able to show the diversity of N. gonorrhoeae strains in WA and suggest that the epidemiology of the various stains is quite different and, in some cases, geographically restricted. The application of NG-MAST or tbpB allele analysis to urogenital specimens uncultured but with detectable N. gonorrhoea DNA is worthy of further study as this would provide a more reliable picture of the prevalence and epidemiology of WA N. gonorrhoeae strains.
The low rates of antimicrobial resistance amongst the gonococci from remote and Aboriginal communities in WA has been attributed to both empirical use of dual amoxicillin/probenecid and azithromycin therapy as well as the social and geographical isolation from the mainstream urban populations [3, 18]. Incursion of strains harbouring antimicrobial resistance markers would therefore pose a serious public health challenge in this population. We did not find STs associated with international spread of multi-drug resistance phenotypes, such as ST 1407 or ST 225, in our remote populations although these strains have been described from the Eastern States of Australia [19, 20]. Therefore, given the history of spread of gonococcal antimicrobial resistance, a heightened vigil supported by a robust surveillance system is essential to ensure first line agents such as amoxycillin remain effective for our population.
- Goire N, Lahra MM, Chen M, Donovan B, Fairley CK, Guy R, et al. Molecular approaches to enhance surveillance of gonococcal antimicrobial resistance. Nat Rev Microbiol. 2014;12(3):223–9. doi:10.1038/nrmicro3217.View ArticlePubMedGoogle Scholar
- Lahra MM. Australian Gonococcal Surveillance Programme annual report, 2012. Commun Dis Intell Q Rep. 2013;37(3):E233–9.PubMedGoogle Scholar
- Speers DJ, Fisk RE, Goire N, Mak DB. Non-culture Neisseria gonorrhoeae molecular penicillinase production surveillance demonstrates the long-term success of empirical dual therapy and informs gonorrhoea management guidelines in a highly endemic setting. J Antimicrob Chemother. 2014;69(5):1243–7. doi:10.1093/jac/dkt501.View ArticlePubMedGoogle Scholar
- Whiley DM, Goire N, Ray ES, Limnios A, Lambert SB, Nissen MD, et al. Neisseria gonorrhoeae multi-antigen sequence typing using non-cultured clinical specimens. Sex Transm Infect. 2010;86(1):51–5. doi:10.1136/sti.2009.037689.View ArticlePubMedGoogle Scholar
- Carannante A, Prignano G, Cusini M, Matteelli A, Dal Conte I, Ghisetti V, et al. Cefixime and ceftriaxone susceptibility of Neisseria gonorrhoeae in Italy from 2006 to 2010. Clin Microbiol Infect. 2012;18(6):558–64. doi:10.1111/j.1469-0691.2011.03619.x.View ArticlePubMedGoogle Scholar
- Lundback D, Fredlund H, Berglund T, Wretlind B, Unemo M. Molecular epidemiology of Neisseria gonorrhoeae- identification of the first presumed Swedish transmission chain of an azithromycin-resistant strain. APMIS. 2006;114(1):67–71. doi:10.1111/j.1600-0463.2006.apm_332.x.View ArticlePubMedGoogle Scholar
- Martin IM, Foreman E, Hall V, Nesbitt A, Forster G, Ison CA. Non-cultural detection and molecular genotyping of Neisseria gonorrhoeae from a piece of clothing. J Med Microbiol. 2007;56(Pt 4):487–90. doi:10.1099/jmm.0.46956-0.View ArticlePubMedGoogle Scholar
- Sathirareuangchai S, Phuangphung P, Leelaporn A, Boon-Yasidhi V. The usefulness of Neisseria gonorrhoeae strain typing by Pulse-Field Gel Electrophoresis (PFGE) and DNA detection as the forensic evidence in child sexual abuse cases: a case series. Int J Leg Med. 2014. doi:10.1007/s00414-014-1007-z.Google Scholar
- Starnino S, Suligoi B, Regine V, Bilek N, Stefanelli P, Dal Conte I, et al. Phenotypic and genotypic characterization of Neisseria gonorrhoeae in parts of Italy: detection of a multiresistant cluster circulating in a heterosexual network. Clin Microbiol Infect. 2008;14(10):949–54. doi:10.1111/j.1469-0691.2008.02071.x.View ArticlePubMedGoogle Scholar
- Unemo M, Berglund T, OlcEn P, Fredlund H. Pulsed-field gel electrophoresis as an epidemiologic tool for Neisseria gonorrhoeae: identification of clusters within serovars. Sex Transm Dis. 2002;29(1):25–31.View ArticlePubMedGoogle Scholar
- Unemo M, Dillon JA. Review and international recommendation of methods for typing neisseria gonorrhoeae isolates and their implications for improved knowledge of gonococcal epidemiology, treatment, and biology. Clin Microbiol Rev. 2011;24(3):447–58. doi:10.1128/cmr.00040-10.View ArticlePubMedPubMed CentralGoogle Scholar
- Horn NN, Kresken M, Korber-Irrgang B, Gottig S, Wichelhaus C, Wichelhaus TA. Antimicrobial susceptibility and molecular epidemiology of Neisseria gonorrhoeae in Germany. Int J Med Microbiol. 2014;304(5-6):586–91. doi:10.1016/j.ijmm.2014.04.001.View ArticlePubMedGoogle Scholar
- Thakur SD, Levett PN, Horsman GB, Dillon JA. Molecular epidemiology of Neisseria gonorrhoeae isolates from Saskatchewan, Canada: utility of NG-MAST in predicting antimicrobial susceptibility regionally. Sex Transm Infect. 2014;90(4):297–302. doi:10.1136/sextrans-2013-051229.View ArticlePubMedGoogle Scholar
- Whiley DM, Anderson TP, Barratt K, Beaman MH, Buda PJ, Carter M, et al. Evidence that the gonococcal porA pseudogene is present in a broad range of Neisseria gonorrhoeae strains; suitability as a diagnostic target. Pathology. 2006;38(5):445–8. doi:10.1080/00313020600928253.View ArticlePubMedGoogle Scholar
- Hunter PR, Gaston MA. Numerical index of the discriminatory ability of typing systems: an application of Simpson’s index of diversity. J Clin Microbiol. 1988;26(11):2465–6.PubMedPubMed CentralGoogle Scholar
- World Health Organization DoRHaR. Global strategy for the prevention and control of sexually transmitted infections: 2006 - 2015. Breaking the chain of transmission. In: Organization WH, editor. Geneva, Switzerland; 2007.Google Scholar
- Bowden FJ. Controlling HIV, in Indigenous Australians. Med J Aust. 2005;183(3):116–7.PubMedGoogle Scholar
- Whiley DM, Goire N, Lahra MM, Donovan B, Limnios AE, Nissen MD, et al. The ticking time bomb: escalating antibiotic resistance in Neisseria gonorrhoeae is a public health disaster in waiting. J Antimicrob Chemother. 2012;67(9):2059–61. doi:10.1093/jac/dks188.View ArticlePubMedGoogle Scholar
- Goire N, Freeman K, Tapsall JW, Lambert SB, Nissen M, Sloots T, et al. Enhancing gonococcal antimicrobial resistance surveillance: a real-time PCR assay for detection of penicillinase-producing Neisseria gonorrhoeae by use of noncultured clinical samples. J Clin Microbiol. 2011;2:513–8.View ArticleGoogle Scholar
- Goire N, Freeman K, Lambert SB, Nimmo GR, Limnios AE, Lahra MM, et al. The influence of target population on nonculture-based detection of markers of Neisseria gonorrhoeae antimicrobial resistance. Sex Health. 2012;9(5):422–9. doi:10.1071/sh12026.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.