https://ojs.wpro.who.int/ 1WPSAR Vol 12, No 1, 2021 | doi: 10.5365/wpsar.2020.11.1.005 Original Research a Antimicrobial Resistance Surveillance Reference Laboratory, Research Institute for Tropical Medicine, Muntinlupa, Philippines. b Centre for Genomic Pathogen Surveillance, Wellcome Genome Campus, Hinxton, England, United Kingdom. c Brigham and Women’s Hospital, Boston (MA), USA. d University of St Andrews School of Medicine, St Andrews, Scotland, United Kingdom. e Centre for Genomic Pathogen Surveillance, Big Data Institute, University of Oxford, Oxford, England, United Kingdom. † These authors contributed equally to this work. * These authors contributed equally to this work. Published: 26 February 2021 doi: 10.5365/wpsar.2020.11.1.005 N eisseria gonorrhoeae is a leading cause of sexually transmitted infections, with an estimated 78 million cases of gonorrhoea each year worldwide, including 35.2 million in the WHO Western Pacific Region.1 In the Philippines, the prevalence of gonorrhoea in 2002 was reported to be <2% for both men and women,2 with higher rates of 7.7% and 10.8% among men who have sex with men at two different sites in 2005.3 N. gonorrhoeae has developed resistance to first-line antibiotics such as sulfonamides, penicillins, tetracyclines, macrolides, fluoroquinolones and early cephalosporins. Currently recommended monotherapy for gonorrhoea is limited to one last effective class of an- timicrobials, the extended-spectrum cephalosporins (e.g. cefixime and ceftriaxone); however, because of the recent emergence of resistance to these drugs, dual therapy with the injectable ceftriaxone plus oral azithromycin is the recommended treatment in many countries.4 While resistance to azithromycin has also increased globally,1 resistance to the dual therapy remains low.5 The increase in N. gonorrhoeae infections resistant to front-line antibiotics triggered a global action plan from WHO to control the spread and impact of gonococcal resistance and a call for international collaborative ac- tion, especially in the Western Pacific Region.1 The WHO Antimicrobial-resistant Neisseria gonorrhoeae is a major threat to public health and is of particular concern in the Western Pacific Region, where the incidence of gonorrhoea is high. The Antimicrobial Resistance Surveillance Program (ARSP) has been capturing information on resistant gonorrhoea since 1996, but genomic epidemiology studies on this pathogen are lacking in the Philippines. We sequenced the whole genomes of 21 N. gonorrhoeae isolates collected in 2013–2014 by ARSP. The multilocus sequence type, multiantigen sequence type, presence of determinants of antimicrobial resistance and relatedness among the isolates were all derived from the sequence data. The concordance between phenotypic and genotypic resistance was also determined. Ten of 21 isolates were resistant to penicillin, ciprofloxacin and tetracycline, due mainly to the presence of the blaTEM gene, the S91F mutation in the gyrA gene and the tetM gene, respectively. None of the isolates was resistant to ceftriaxone or cefixime. The concordance between phenotypic and genotypic resistance was 92.38% overall for five antibiotics in four classes. Despite the small number of isolates studied, they were genetically diverse, as shown by the sequence types, the N. gonorrhoeae multiantigen sequence typing types and the tree. Comparison with global genomes placed the Philippine genomes within global lineage A and led to the identification of an international transmission route. This first genomic survey of N. gonorrhoeae isolates collected by ARSP will be used to contextualize prospective surveillance. It highlights the importance of genomic surveillance in the Western Pacific and other endemic regions for understanding the spread of drug-resistant gonorrhoea worldwide. Genomic surveillance of Neisseria gonorrhoeae in the Philippines, 2013–2014 Manuel C. Jamoralin Jr,a† Silvia Argimón,b† Marietta L. Lagrada,a Alfred S. Villamin,a Melissa L. Masim,a June M. Gayeta,a Karis D. Boehme,a Agnettah M. Olorosa,a Sonia B. Sia,a Charmian M. Hufano,a Victoria Cohen,b Lara T. Hernandez,a Benjamin Jeffrey,b Khalil Abudahab,b John Stelling,c Matthew T.G. Holden,d David M. Aanensenb,e,* and Celia C. Carlosa,* on behalf of the Philippines Antimicrobial Resistance Surveillance Program Correspondence to Celia Carlos (ccarlosphl@gmail.com) and David M. Aanensen (David.aanensen@bdi.ox.ac.uk) WPSAR Vol 12, No 1, 2021 | doi: 10.5365/wpsar.2020.11.1.005 https://ojs.wpro.who.int/2 Villamin et alN. gonorrhoeae surveillance in the Philippines antimicrobials representing four different classes, namely penicillin (PEN), ciprofloxacin (CIP), tetracycline (TCY), ceftriaxone (CRO) and cefixime (CFX), by disc diffusion and gradient diffusion (Etest, BioMerieux). The minimum inhibitory concentrations were interpreted as resistant, intermediate or susceptible according to the interpretative criteria in the Performance Standards for Antimicrobial Susceptibility Testing (26th edition) of the Clinical and Laboratory Standards Institute (CLSI).12 All isolates were screened for β-lactamase production on cefinase paper discs (BD BBL). Azithromycin was not included in the panel of antibiotics, because, at the time isolates were collected, it was not in the treatment guidelines of the Philippines and no CLSI breakpoint was available. DNA extraction and whole-genome sequencing DNA was extracted from a single colony of each isolate with a Wizard® Genomic DNA Purification Kit (Pro- mega), and the quantity and quality were determined with a Quantus Fluorometer (Promega) with picogreen and a NanoDrop 2000c spectrophotometer (Thermo Fisher Scientific). The DNA extracts were then shipped to Wellcome Trust Sanger Institute for sequencing on the Illumina HiSeq platform with 100-bp paired-end reads. Raw sequence data were deposited in the European Nucleotide Archive under the study accession number PRJEB17615. Run accessions are provided on the Mi- croreact projects. Bioinformatics analysis Genome quality was evaluated according to metrics gen- erated from assemblies, annotation files and the align- ment of the isolates to the reference genome of N. gonor- rhoeae strain TCDC-NG08107 (accession CP002441.1), as previously described.13 Annotated assemblies were produced with the pipeline previously described.14 We included 21 high-quality N. gonorrhoeae genomes in this study. The MLST sequence types (STs) and NG-MAST types, as well as the presence of AMR determinants (known genes or mutations) and clustering of the isolates according to genetic similarity (tree), were predicted in silico from genome assemblies with Pathogenwatch (https://www.sanger.ac.uk/tool/pathogenwatch/).15 In parallel, the evolutionary relations among the isolates were inferred from single nucleotide polymorphisms Gonococcal Antimicrobial Surveillance Programme has operated in the Western Pacific and South-East Asian regions since 1992, but surveillance of gonococcal anti- microbial resistance (AMR) remains limited in the Asia– Pacific region.6 In a recent report, 18 of 21 countries in the Asia and the Pacific reported isolates with decreased susceptibility to ceftriaxone and/or isolates resistant to azithromycin between 2011 and 2016.6 The Antimi- crobial Resistance Surveillance Program (ARSP) of the Philippines Department of Health has been contributing AMR surveillance data to the Western Pacific Gonococcal Antimicrobial Surveillance Programme since 1996 and did not confirm isolates with decreased susceptibility or resistance to these antibiotics during 2011–2016,6 while high gonococcal resistance rates against other first-line antibiotics have long been reported (Fig. 1). Continuous surveillance is thus key to detecting potential emergence or introduction of resistance to current treatment options. Molecular methods for defining the epidemiology of gonococci include both N. gonorrhoeae multiantigen sequence typing (NG-MAST)7 and multilocus sequence typing (MLST),8 although NG-MAST is more widely used to investigate specific gonococcal AMR phenotypes.8 Whole-genome sequencing (WGS) was recently shown to provide better resolution and accuracy than NG-MAST or MLST.9 Good understanding of the population structure and the mechanisms of resistance of N. gonorrhoeae in the Philippines would allow detection of high-risk clones associated with high-risk groups and contribute to the clinical management of gonococcal-related diseases and the creation of policies to prevent the spread of drug resistance.10,11 Here, we describe the results of the first genomic survey of gonococcal isolates in the Philippines. METHODS Bacterial isolates A total of 51 N. gonorrhoeae isolates were collected at ARSP sentinel sites in 2013 and 2014 (Table 1). Of the 36 isolates referred to the ARSP reference laboratory for confirmation, 22 isolates from seven sentinel sites were resuscitated and submitted for WGS. Antimicrobial susceptibility testing All N. gonorrhoeae isolates in this study were tested at the ARS reference laboratory for susceptibility to five WPSAR Vol 12, No 1, 2021 | doi: 10.5365/wpsar.2020.11.1.005https://ojs.wpro.who.int/ 3 N. gonorrhoeae surveillance in the PhilippinesVillamin et al (SNPs) by mapping the paired-end reads to the reference genome of N. gonorrhoeae strain FA 1090 (accession AE004969.1), as described in detail previously.13 Mobile genetic elements described in the FA 1090 genome16 were masked in the alignment of pseudogenomes with a script available at https://github.com/sanger-pathogens/ remove_blocks_from_aln. SNPs were extracted with snp-sites,17 and a maximum likelihood phylogenetic tree was generated from 7518 SNP positions with RAxML,18 the generalized time-reversible model and the GAMMA method of correction for among-site rate variation, with 500 bootstrap replicates. To complement the Pathogenwatch AMR results, known AMR determinants were identified from raw sequence reads with ARIBA19 and a curated database of known resistance genes and mutations available at https://github.com/martinghunt/ariba-publication/tree/ master/N_gonorrhoeae/Ref. The combined genotypic predictions of AMR (test) were compared with the phe- notypic results (reference), and the concordance be- tween the two methods was computed for each of five antibiotics (105 total comparisons). Isolates found to be resistant or to have reduced susceptibility (intermedi- ate) were pooled as non-susceptible for comparison purposes. An isolate with the same outcome for both the test and reference (i.e. both susceptible or both non-susceptible) was counted as a concordant isolate. Concordance was the number of concordant isolates Fig. 1. Annual resistance rates of N. gonorrhoeae between 2000 and 2014 for penicillin (PEN), ciprofloxacin (CIP) and tetracycline (TCY) Table 1. Number of N. gonorrhoeae isolates analysed by ARSP and referred to the reference labo- ratory during 2013 and 2014, isolates sub- mitted for WGS and high-quality N. gonor- rhoeae genomes obtained, by sentinel site and AMR profile CVM: Cagayan Valley Medical Center, DMC: Southern Philippines Medical Center, EVR: Eastern Visayas Regional Medical Center, MMH: Corazon Locsin Montelibano Memorial Regional Hospital, NMC: Northern Mindanao Medical Center, VSM: Vicente Sotto Memorial Medical Center, ZPH: Zamboanga Del Norte Medical Center PEN: penicillin, CIP: ciprofloxacin, TCY: tetracycline Number of isolates 2013 2014 Total Total ARSP 24 27 51 Referred to reference laboratory 16 20 36 Submitted for WGS 8 14 22 High-quality genomes 8 13 21 By sentinel site CVM 0 1 1 DMC 1 0 1 EVR 1 0 1 MMH 1 1 2 NMC 0 1 1 VSM 5 6 11 ZPH 0 4 4 By AMR profile PEN CIP TCY 2 8 10 PEN CIP 4 4 8 PEN 2 0 2 CIP TCY 0 1 1 WPSAR Vol 12, No 1, 2021 | doi: 10.5365/wpsar.2020.11.1.005 https://ojs.wpro.who.int/4 Villamin et alN. gonorrhoeae surveillance in the Philippines mtrR promoter and the non-synonymous substitutions in ponA (L421P) and porB (A121D), which may contribute to high-level penicillin resistance.22 The partially assembled blaTEM gene was also detected in this genome but was not considered in the concordance analysis. We also identified the non-synonymous mutation in the folP gene (R228S), which confers resistance to sulfonamides in 20 of 21 genomes; however, isolates are not routinely tested for resistance to this antibiotic class. The concordance between antimicrobial susceptibil- ity testing results and genotypic predictions (Table 3) was >95% for all antibiotics except tetracycline (66.67%), resulting in an overall concordance of 92.38% (Table 3). The discrepancies were attributed to seven false-positive results (isolates with a susceptible phenotype but with known AMR determinants in their genomes), all of which contained the mutation V57M in the rpsJ gene; two iso- lates also carried the tetM gene.23 over the total number of isolates assessed (expressed as a percentage). The maximum likelihood tree, genotyping results and AMR predictions and the metadata collected from the sentinel sites were visualized with Microreact.20 To contextualize the genomes from this study with publicly available global genomes, we combined them with two surveys available on Pathogenwatch, a European survey of 1054 genomes10 and a global survey of 395 genomes.21 Ethics statement Ethical approval is not applicable, as we used archived bacterial samples processed by ARSP. No identifiable data were used in this study. RESULTS Demographic and clinical characteristics of the N. gonorrhoeae isolates The 21 genomes included in this study represented seven sentinel sites, with Vicente Sotto Memorial Medical Center (VSM) contributing the most isolates (n = 11). The high- est incidence was in the age group 15–24 years (47.6%, n = 10), followed by the age groups 5–14 years (28.6%, n = 6), 25–34 years (14.3%, n = 3) and 1–4 years (9.5%, n = 2, Table 2). The numbers of isolates from females and males were almost equal (n = 11 and n = 10, respectively). The most frequent specimen source was the vagina (n = 8) for female patients and penile dis- charge (n = 5) for males. All the patients were outpatients (n = 21). Concordance between phenotypic and geno- typic AMR Isolates were tested for susceptibility to five antibiotics representing four classes. In line with the resistance trends shown in Fig. 1, the most prevalent resistance profile was PEN CIP TCY, identified in 10 isolates from five sentinel sites and linked mainly to the presence of the blaTEM gene, the S91F mutation in GyrA and the tetM gene, respectively (Ta- ble 3). One penicillin-resistant isolate (13ARS_DMC0024) harboured three mutations, the –57delA mutation in the Table 2. Demographic and clinical characteristics of 21 N. gonorrhoeae isolates Characteristic No. of isolates Sex Male 10 Female 11 Age (years) 1–4 2 5–14 6 15–24 10 25–34 3 Patient type Inpatient 0 Outpatient 21 Location Community 21 Hospital 0 Specimen origin Vagina 8 Penile discharge 5 Cervix 3 Genital discharge, male 2 Decubitus ulcer 1 Urine 1 Other 1 WPSAR Vol 12, No 1, 2021 | doi: 10.5365/wpsar.2020.11.1.005https://ojs.wpro.who.int/ 5 N. gonorrhoeae surveillance in the PhilippinesVillamin et al Table 3. Comparison of genomic predictions of antibiotic resistance with susceptibility testing at the ARS refer- ence laboratory Antibiotic class Antibiotic Resistant isolates False positive False negative Concordance (%) Resistance genes/ SNPs Cephalosporin Cefixime 0 0 0 100 Cephalosporin Ceftriaxone 0 0 0 100 Penicillin Penicillin 18 0 1 95.24 blaTEM, mtrR_ promoter_–57de- lA, mtrR_A39T, mtrR_disrupted, penA_ins346D, ponA1_L421P, porB1b_A121D, porB1b_G120K Tetracycline Tetracycline 14 7 0 66.67 tetM, rpsJ-V57M, mtrR_promoter_ –57delA, mtrR_ A39T, mtrR_dis- rupted Fluoroquinolone Ciprofloxacin 19 0 0 100 gyrA_D95G, gyrA_D95A, gyrA_S91F, parC_D86N, parC_S87N, parC_E91K contained isolates resistant to ciprofloxacin harbouring the gyrA_S91F mutation, clade II was characterized by the presence of gyrA_D95G alone or in combination with the one mutation in parC, while clade III was characterized by the presence of gyrA_D95A with one or two mutations in parC in all but one isolate. Clades II and III showed differ- ent geographical distributions, although both were present at the VSM and ZPH sentinel sites, which submitted the most isolates. The genome from sentinel site Southern Philippines Medical Center (DMC) with a deletion in the mtrR promoter (–57delA) was found on a separate branch (I) in the tree and also carried a different complement of AMR determinants, indicating that it is genetically distinct from the others (Fig. 2). N. gonorrhoeae from the Philippines in the global context The Philippine genomes were contextualized with two recently published collections10,21 available in Pathogen- watch. A recent global collection of 395 genomes from 58 countries (including the Philippines) showed two major lineages with different evolutionary strategies. Most of the genomes in this study were found within a subclade of lineage A, a multidrug-resistant lineage associated with infection in high-risk sexual networks (Fig. 3A), and mixed with genomes from Europe, Pakistan and South-East Asia, including genomes previously isolated in the Philippines Genotypic findings In silico genotyping A total of 15 different STs were identified, with only four (9364, 10316, 1582 and 8133) represented by more than one isolate. Nine genomes were assigned to eight known NG-MAST types and 12 genomes to nine novel types. Only three sentinel sites, namely Corazon Locsin Montelibano Memorial Regional Hospital (MMH), VSM and Zamboanga Del Norte Medical Center (ZPH), submitted more than one isolate, and all were represented by almost as many STs as isolates submitted. Details of the numbers and the most common STs and NG-MAST types found at each sentinel site are shown in Table 4. Population structure of N. gonorrhoeae in the Philip- pines The diverse gonococcal population was represented by a tree with three deep-branching clades and no clear geographical signal (Fig. 2). Clades II and III were characterized by a different repertoire of AMR genes and mutations. Clade II contained mostly isolates susceptible to or with reduced susceptibility to tetracycline and with the V57M mutation in the rpsJ gene alone, while clade III was composed of isolates resistant to tetracycline and also containing the tetM gene. Similarly, while both clades WPSAR Vol 12, No 1, 2021 | doi: 10.5365/wpsar.2020.11.1.005 https://ojs.wpro.who.int/6 Villamin et alN. gonorrhoeae surveillance in the Philippines Table 4. Distribution of isolates, STs, NG-MAST types, resistance profiles and resistance genes and mutations at the seven sentinel sites. The genetic resistance mechanisms for all the isolates from each site are listed. Fig. 2. Genomic surveillance of N. gonorrhoeae from the Philippines, 2013–2014 CVM: Cagayan Valley Medical Center; DMC: Southern Philippines Medical Center; EVR: Eastern Visayas Regional Medical Center; MMH: Corazon Locsin Monteli- bano Memorial Regional Hospital; NMC: Northern Mindanao Medical Center; VSM: Vicente Sotto Memorial Medical Center; ZPH: Zamboanga Del Norte Medical Center PEN: penicillin, CIP: ciprofloxacin, TCY: tetracycline Maximum likelihood phylogenetic tree of 21 isolates from the Philippines, inferred from 7518 SNP sites. The tree leaves are coloured by sentinel site (map). The tree is annotated with the MLST ST as per the legend, the resistance phenotype for five antibiotics (red: resistant, yellow: intermediate, green: susceptible), and the distribution of known AMR mechanisms (red: present, grey: not found). The tree branches are annotated with the bootstrap values and clade designations (I, II and III). The scale bar represents the number of SNPs per variable site. The data are available at https://microreact.org/project/ARSP_NGO_2013-14 and https://patho- gen.watch/collection/flnkisnu6giu-arspngo2013-2014. Laboratory No. of isolates STs (no. of isolates) NG-MAST types (no. of isolates) Resistance profile± (no. of isolates) Resistance genes and SNPs CVM 1 8780 Novel PEN CIP TCY blaTEM, penA_ins346D, ponA1_L421P, mtrR_ A39T, gyrA_S91F, gyrA_D95A, parC_S87N, tetM, rpsJ-V57M DMC 1 1903 Novel PEN CIP TCY penA_ins346D, ponA1_L421P, porB1b_ A121D, mtrR_promoter_ –57delA, gyrA_D95G, gyrA_S91F, parC_S87N, parC_E91K, tetM, rpsJ-V57M EVR 1 1582 13 796 PEN CIP blaTEM, penA_ins346D, mtrR_A39T, gyrA_ D95G, gyrA_S91F, rpsJ-V57M MMH 2 11 956, 11 208 1631, novel PEN CIP blaTEM, penA_ins346D, mtrR_A39T, porB1b_ A121D, porB1b_G120K, gyrA_D95G, gyrA_ S91F, parC_D86N, rpsJ_V57M NMC 1 1588 Novel PEN CIP TCY blaTEM, penA_ins346D, mtrR_A39T, ponA1_ L421P, gyrA_S91F, gyrA_D95A tetM, rpsJ- V57M VSM 11 8133 (2), 9364 (2), 10 316 (1), 1587 (1), 11 431 (1), 11 963 (1), 1582 (1), 11 956 (1), 15 234 (1) 2187(2), 1691 (1), 1498 (1),1797 (1), novel (6) PEN CIP (4) PEN CIP TCY (4) CIP TCY (1) PEN (2) blaTEM, penA_ins346D, mtrR_disrupted, mtrR_A39T, ponA1_L421P, porB1b_A121D, porB1b_G120K, gyrA_D95G, gyrA_D95A, gyrA_S91F, parC_D86N, parC_S87N, parC_ E91K, tetM, rpsJ-V57M ZPH 4 10 316 (2), 9364 (1), 8130 (1) 2080 (1), 11 821 (1), novel (2) PEN CIP TCY (3) PEN CIP (1) blaTEM, mtrR_disrupted, mtrR_A39T, porB1b_ A121D, gyrA_D95G, gyrA_D95A, gyrA_S91F, parC_D86N, parC_S87N, parC_E91K, tetM, rpsJ-V57M ST legend S T P E N C IP TC Y C FM C R O b la T E M m tr R _p ro m ot er _a –5 7 de l m tr R _A 3 9 T m tr R _d is ru p te d p en A _i ns 3 4 6 D p on A 1 _L 4 2 1 P p or B 1 b _A 1 2 1 D p or B 1 b _G 1 2 0 K gy rA _d 9 5 G gy rA _D 9 5 A gy rA _S 9 1 F p ar C _D 8 6 N p ar C _S 8 7 N p ar C _E 9 1 K te tM rp sJ _V 5 7 M 100 100 II I III 100 100 100 100 100 100 95 WPSAR Vol 12, No 1, 2021 | doi: 10.5365/wpsar.2020.11.1.005https://ojs.wpro.who.int/ 7 N. gonorrhoeae surveillance in the PhilippinesVillamin et al Fig. 3. Global context of N. gonorrhoeae from the Philippines Fig. 3A. Tree of 416 genomes from this study (n = 21) and a recent global collection (n = 395), inferred with Pathogenwatch from 22 558 variant sites in 1542 core genes. Purple nodes indicate genomes from the Philippines from this study (circles) or the global collection (squares). The scale bar represents the number of SNPs. The data are available at https://pathogen.watch/collection/wrdonfwhju6f-arsp-ngo-2013-2014-global-context. Fig. 3B. Detail of the subtree of closely related genomes from the Philippines (orange nodes) and Norway (blue node), inferred with Pathogenwatch. The tree branches are annotated with the number of pairwise SNP differences between isolates. The metadata blocks indicate the ST, NG-MAST type and the presence (red blocks) of seven AMR determinants. The full collection is available at https://pathogen.watch/collection/xtusqgwqhxcy-arsp-ngo-2013-2014-european-context. two genomes is consistent with isolates from the same location,10 suggesting an epidemiological link and a route of international transmission. The gonococcal reference laboratory at the Norwegian Institute of Public Health confirmed that the isolate was also resistant to penicillin and ciprofloxacin and susceptible to extended-spectrum cephalosporins (tetracycline not tested) and that the Norwegian male had visited the Philippines and claimed to have contracted gonorrhoea during his stay (personal communication). DISCUSSION WGS showed that the N. gonorrhoeae genomes from the Philippines are genetically diverse and carry a variety of AMR determinants, such as chromosomal mutations and acquired genes. The concordance between phenotypic (1998 and 200821). The genome from DMC was found within a separate subclade with genomes from Europe, India and Pakistan, which also carried the mtrR_–57delA promoter deletion. We further contextualized our isolates with 1054 ge- nomes from 20 countries collected in 2013 in a European survey.10 Notably, the genome of strain 14ARS_VSM0347 isolated from a female on 25/04/2014 in the Philippines was highly similar to that of ECDC_GC_088 isolated from a male on 31/12/2013 in Norway (Fig. 3B), with no SNP differences found in the Pathogenwatch core genome (1542 genes). One SNP was identified in the reference- based alignment of the two pseudogenomes, confirming that these two genomes are highly similar. The two strains also shared the same complement of AMR genes and mutations (Fig. 3B). The genetic distance between the Lineage B Lineage AA B 723 0 MLST ST NG-MAST type 11956 11963 1498 1631 WPSAR Vol 12, No 1, 2021 | doi: 10.5365/wpsar.2020.11.1.005 https://ojs.wpro.who.int/8 Villamin et alN. gonorrhoeae surveillance in the Philippines portance of genomic surveillance in the Western Pacific and other endemic regions to understand the spread of drug-resistant gonorrhoea worldwide. Acknowledgements We are grateful to Martin Steinbakk at the Norwegian Institute of Public Health and to Karianne Wiger Gam- melsrund and the microbiology laboratory at the University Hospital of North Norway for retrieving the information on the Norwegian isolate. We also thank Dr Leonor Sánchez-Busó and Dr Magnus Unemo for critical reading of the manuscript. Funding This work was supported by a Newton Fund award from the Medical Research Council (United Kingdom) MR/ N019296/1 and the Philippine Council for Health Re- search and Development. Additional support was provid- ed by the National Institute for Health Research (United Kingdom), the Global Health Research Unit on genomic surveillance of AMR (16/136/111) and by research grant U01CA207167 from the US National Institutes of Health. Conflicts of interest The authors declare no conflicts of interest. References 1. Wi T, Lahra MM, Ndowa F, Bala M, Dillon JR, Ramon-Pardo P, et al. Antimicrobial resistance in Neisseria gonorrhoeae: Global surveillance and a call for international collaborative action. PLoS Med. 2017;14(7):e1002344. 2. Wi TEC, Saniel OP, Ramos ER, MacNeil J, Neilsen G. RTI/STI prevalence in selected sites in the Philippines. Durham (NC): FHI Foundation; 2002. 3. External review of the national health sector response to HIV and sexually transmitted infections 2013: Republic of the Phil- ippines. Manila: WHO Regional Office foe the Western Pacific; 2015. 4. Unemo M, Shafer WM. Antimicrobial resistance in Neisseria gon- orrhoeae in the 21st century: past, evolution, and future. Clin Mi- crobiol Rev. 2014;27(3):587–613. 5. Fifer H, Natarajan U, Jones L, Alexander S, Hughes G, Golparian D, et al. Failure of dual antimicrobial therapy in treatment of gonor- rhea. N Engl J Med. 2016;374(25):2504–6. 6. George CRR, Enriquez RP, Gatus BJ, Whiley DM, Lo YR, Ishikawa N, et al. Systematic review and survey of Neisseria gonorrhoeae ceftriaxone and azithromycin susceptibility data in the Asia Pa- cific, 2011 to 2016. PLoS One. 2019;14(4):e0213312. 7. 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. and genotypic resistance was high (>95%) for most anti- biotics but only 66% for tetracycline. Susceptible isolates carrying only the rpsJ_V57M mutation, which confers low-level tetracycline resistance, have been reported pre- viously.24 The two susceptible isolates with a full-length tetM gene reported by Pathogenwatch were re-tested by disc diffusion and their susceptibility confirmed. The discrepancy, which could be explained for example by the lack of expression of the gene, errors in susceptibility testing or low-level DNA contamination, will be further investigated. Although N. gonorrhoeae remains largely susceptible to extended-spectrum cephalosporins and azithromycin in the Philippines, we identified one isolate (13ARS_DMC0024) with a –57delA mutation in the mtrR promoter, which results in overexpression of the MtrCDE efflux pump and, in combination with other mutations, can increase the minimum inhibitory concentrations of azithromycin (as well as penicillin and tetracycline).4,22 Prospective surveillance with WGS can detect in real-time the acquisition of additional mutations that could result in decreased susceptibility. A recent report suggested rapid recent intercontinen- tal transmission of gonorrhoea, with common introduc- tions from Asia to the rest of the world.21 In support of this finding, the Philippine genomes, most of which were within a subclade of global lineage A, were interspersed with those from other countries (Fig. 3A). In addition, we found evidence of an introduction event from the Philip- pines to Norway associated with travel of a Norwegian male to the Philippines (Fig. 3B). Global lineage A is associated with infection in high- risk sexual networks.21 The combination of WGS with epidemiological information can reveal transmission routes and risk factors, which can be used to design better control measures.25 The small size of our retrospective data set (n = 21) and the linked epidemiological data did not permit any inferences about sexual networks or risk factors, which is a limitation of this study. The number of reported isolates by ARSP has, however, since increased to more than 100 per year, and data on risk factors are also collected, which will allow a more comprehensive analysis of the population diversity and of risk factors in future reports. Our results represent the first genomic survey of N. gonorrhoeae isolates collected by ARSP and will constitute the background for contextualizing continuous prospective surveillance. In addition, it highlights the im- WPSAR Vol 12, No 1, 2021 | doi: 10.5365/wpsar.2020.11.1.005https://ojs.wpro.who.int/ 9 N. gonorrhoeae surveillance in the PhilippinesVillamin et al 17. Page AJ, Taylor B, Delaney AJ, Soares J, Seemann T, Keane JA, et al. SNP-sites: rapid efficient extraction of SNPs from multi- FASTA alignments. Microb Genom. 2016;2(4):e000056. 18. Stamatakis A. RAxML version 8: a tool for phylogenetic anal- ysis and post-analysis of large phylogenies. Bioinformatics. 2014;30(9):1312–3. 19. Hunt M, Mather AE, Sanchez-Buso L, Page AJ, Parkh- ill J, Keane JA, et al. ARIBA: rapid antimicrobial resistance genotyping directly from sequencing reads. Microb Genom. 2017;3(10):e000131. 20. Argimon S, Abudahab K, Goater RJ, Fedosejev A, Bhai J, Glasner C, et al. Microreact: visualizing and sharing data for genomic epidemiology and phylogeography. Microb Genom. 2016;2(11):e000093. 21. Sánchez-Busó L, Golparian D, Corander J, Grad YH, Ohnishi M, Flemming R, et al. The impact of antimicrobials on gonococcal evolution. Nat Microbiol. 2019;4(11):1941–50. 22. Mortimer TD, Grad YH. Applications of genomics to slow the spread of multidrug-resistant Neisseria gonorrhoeae. Ann N Y Acad Sci. 2019;1435(1):93–109. 23. Eyre DW, De Silva D, Cole K, Peters J, Cole MJ, Grad YH, et al. WGS to predict antibiotic MICs for Neisseria gonorrhoeae. J Antimicrob Chemother. 2017;72(7):1937–47. 24. Bailey AL, Potter RF, Wallace MA, Johnson C, Dantas G, Burn- ham CA. Genotypic and phenotypic characterization of antimi- crobial resistance in Neisseria gonorrhoeae: a cross-sectional study of isolates recovered from routine urine cultures in a high- incidence setting. MSphere. 2019;4(4):e00373-19. 25. Didelot X, Dordel J, Whittles LK, Collins C, Bilek N, Bishop CJ, et al. Genomic analysis and comparison of two gonorrhea out- breaks. MBio. 2016;7(3):e00525-16. 8. Jolley KA, Maiden MC. BIGSdb: Scalable analysis of bacterial genome variation at the population level. BMC Bioinformatics. 2010;11(1):595. 9. Unemo M, Golparian D, Nicholas R, Ohnishi M, Gallay A, Sed- naoui P. High-level cefixime- and ceftriaxone-resistant Neisseria gonorrhoeae in France: novel penA mosaic allele in a successful international clone causes treatment failure. Antimicrob Agents Chemother. 2012;56(3):1273–80. 10. Harris SR, Cole MJ, Spiteri G, Sanchez-Buso L, Golparian D, Jacobsson S, et al. Public health surveillance of multidrug-resist- ant clones of Neisseria gonorrhoeae in Europe: a genomic survey. Lancet Infect Dis. 2018;18(7):758–68. 11. Williamson DA, Chow EPF, Gorrie CL, Seemann T, Ingle DJ, Hig- gins N, et al. Bridging of Neisseria gonorrhoeae lineages across sexual networks in the HIV pre-exposure prophylaxis era. Nat Commun. 2019;10(1):3988. 12. M100–S26 performance standards for antimicrobial susceptibil- ity testing. 26th edition. Annapolis Junction (MD): Clinical and Laboratory Standards Institute; 2016. 13. Argimon S, Masim MAL, Gayeta JM, Lagrada ML, Macaranas PKV, Cohen V, et al. Integrating whole-genome sequencing with- in the National Antimicrobial Resistance Surveillance Program in the Philippines. Nat Commun. 2020;11(1):2719. 14. Page AJ, De Silva N, Hunt M, Quail MA, Parkhill J, Harris SR, et al. Robust high-throughput prokaryote de novo assembly and improvement pipeline for Illumina data. Microb Genom. 2016;2(8):e000083. 15. Sánchez-Busó L, Yeats CA, Taylor B, Goater R, Underwood A, Abu- dahab K, et al. A community-driven resource for genomic surveillance of Neisseria gonorrhoeae at Pathogenwatch. bioRxiv. 2020. 16. Piekarowicz A, Klyz A, Majchrzak M, Adamczyk-Poplaws- ka M, Maugel TK, Stein DC. Characterization of the dsDNA prophage sequences in the genome of Neisseria gonorrhoeae and visualization of productive bacteriophage. BMC Microbiol. 2007;7(1):66.