BackgroundWhole genome sequencing (WGS) of Neisseria gonorrhoeae isolates combined with epidemiological and phenotypic data provides better understanding of population dynamics.AimThe objective of this study was to investigate the molecular epidemiology of N. gonorrhoeae isolates from three centres in Spain and determine associations of antimicrobial resistance.MethodsGenetic characterisation was performed in 170 N. gonorrhoeae isolates. WGS was carried out with the HiSeq platform (Illumina). Genome assemblies were submitted to the PubMLST Neisseria database website to determine NG‐MAST, MLST and NG‐STAR. Antimicrobial resistance genes and point mutations were identified with PubMLST. Phylogenomic comparison was based on whole‐genome single nucleotide polymorphism analysis.ResultsTwenty‐six MLST, 49 NG‐MAST and 41 NG‐STAR sequence types were detected, the most prevalent being MLST‐ST9363 (27.1%), NG‐MAST ST569 (12.4%) and NG‐STAR ST193 (14.7%). Phylogenetic analysis identified 13 clusters comprising 69% of the isolates, with two of note: one involved cefixime‐resistant isolates from Barcelona presenting a mosaic penA X and belonging to MLST‐ST7363 and the other involved azithromycin‐resistant isolates from Mallorca that possessed the C2611T mutation in the four 23S rRNA alleles belonging to MLST‐ST1901.ConclusionThe population of N. gonorrhoeae is quite heterogeneous in Spain. Our results agree with previous data published in Europe, albeit with some differences in distribution between regions. This study describes the circulation of two gonococcal populations with a specific resistance profile and sequence type in a specific geographic area. WGS is an effective tool for epidemiological surveillance of gonococcal infection and detection of resistance genes.