A Whole-genome Sequencing Analysis of Neisseria gonorrhoeae Isolates in China: An Observational Study

Background: Tracking the spread of the Neisseria gonorrhoeae strains with decreased susceptibility or resistance to cephalosporins is a major priority for global surveillance programmes. Whole-genome sequencing (WGS) has been widely used by increasing countries in North America, Europe, and Pacific to determine the decreased susceptible or resistance determinants of Neisseria gonorrhoeae, track the spread of these determinants throughout the gonococcal population at national or regional level. However, no studies to date have examined the genomic epidemiology of gonorrhea in Asia where the antimicrobial resistant strains of Neisseria gonorrhoeae appears to have emerged before disseminating the strains globally. Methods: We obtained clinical isolates and data from the China Gonococcal Resistance Surveillance Programme (China-GRSP) from 2012 to 2013. We sequenced the genomes of 435 clinical isolates of Neisseria gonorrhoeae, including 112 (25.6%) isolates with decreased susceptibility to ceftriaxone (Cfx-DS). We assessed the association between antimicrobial resistance genotype and phenotype. We also compared our data with the whole genome data of the isolates from the USA and the UK in the GenBank. Findings: The most prevalent MLST STs in our gonococcal population were MLST ST7827 (n = 74), followed by ST7365 (n = 58), ST1600 (n = 38), ST7367 (n = 35), and ST7363 (n = 29). MLST ST1901 which was reported as the predominant ST in the US was not found in our population. A total of 2512 strains, including additional 2077 published NG strains, were further included for phylogenetic analysis. It generated two distinct lineages - lineage 1 and lineage 2. Analysis of MLST ST1901 in the database indicate that most of MLST ST1901 isolates in the lineage2.6 were Cfx-DS isolates while all isolates in the lineage 2.1 were sensitive to ceftriaxone (77/110 vs. 0/13; p