Antimicrobial resistance and heterogeneity of Neisseria gonorrhoeae isolated from patients attending sexually transmitted infection clinics in Lusaka, Zambia.

Antimicrobial resistance (AMR) of Neisseria gonorrhoeae is a threat to public health as strains have developed resistance to antimicrobials available for the treatment of gonorrhea. Whole genome sequencing (WGS) can detect and predict antimicrobial resistance to enhance the control and prevention of gonorrhea. Data on the molecular epidemiology of N. gonorrhoeae is sparse in Zambia. This study aimed to determine the genetic diversity of N. gonorrhoeae isolated from patients attending sexually transmitted infection (STI) clinics in Lusaka, Zambia.A cross-sectional study that sequenced 38 N. gonorrhoeae isolated from 122 patients with gonorrhea from 2019 to 2020 was conducted. The AMR profiles were determined by the E-test, and the DNA was extracted using the NucliSens easyMaG magnetic device. Whole genome sequencing was performed on the Illumina NextSeq550 platform. The Bacterial analysis pipeline (BAP) that is readily available at: https://cge.cbs.dtu.dk/services/CGEpipeline-1.1 was used for the identification of the species, assembling the genome, multi-locus sequence typing (MLST), detection of plasmids and AMR genes. Phylogeny by single nucleotide polymorphisms (SNPs) was determined with the CCphylo dataset.The most frequent STs with 18.4% of isolates each were ST⁷³⁶³, ST¹⁹²¹ and ST¹⁵⁸², followed by ST¹⁵⁸³ (13%), novel ST¹⁷⁰²⁶ (7.9%), ST¹⁵⁸⁸ (7.9%), ST¹⁵⁹⁶ (5.3%), ST¹¹¹⁸¹ (5.3%), ST¹¹⁷⁵⁰ (2.6/%) and ST¹¹²⁴¹ (2.6%) among the 38 genotyped isolates. The blaTeM-1B and tetM (55%) was the most prevalent combination of AMR genes, followed by blaTeM-1B (18.4%), tetM (15.8%), and the combination of blaTeM-1B, ermT, and tetL was 2.6% of the isolates. The AMR phenotypes were predicted in ciprofloxacin, penicillin, tetracycline, azithromycin, and cefixime. The combination of mutations 23.7% was gryA (S91F), parC (E91G), ponA (L421) and rpsJ (V57M), followed by 18.4% in gyrA (S91F), ponA (L421P), rpsJ (V57M), and 18.4% in gyrA (D95G, S91F), ponA (L421P), and rpsJ (V57M). The combinations in gyrA (D95G, S91F) and rpsJ (V57M), and gyrA (D95G, S91F), parC (E91F), ponA (L421P) and rpsJ (V57M) were 13.2% each of the isolates. Plasmid TEM-1 (84.2%), tetM (15.8%), and gonococcal genetic island (GGI) was detected in all isolates.This study revealed remarkable heterogeneity of N. gonorrhoeae with bla_TEM-1, tetM, ponA, gyrA, and parC genes associated with high resistance to penicillin, tetracycline, and ciprofloxacin demanding revision of the standard treatment guidelines and improved antimicrobial stewardship in Zambia.