Population genetics of Neisseria gonorrhoeae in a high-prevalence community using a hypervariable outer membrane porB and 13 slowly evolving housekeeping genes

Abstract

Baltimore, Md., is an urban community with a high prevalence of Neisseria gonorrhoeae. Due to partially protective immune responses, introduction of new strains from other host populations, and exposure of N. gonorrhoeae to antibiotics, the phenotypic and genotypic characteristics of the circulating strains can fluctuate over time. Understanding the overall genetic diversity and population structure of N. gonorrhoeae is essential for informing public health interventions to eliminate this pathogen. We studied gonococci population genetics in Baltimore by analyzing a hypervariable and strongly selected outer membrane porB gene and 13 slowly evolving and presumably neutral housekeeping genes (abcZ, adk, aroE, fumC, gdh, glnA, gnd, pdhC, pgm, pilA, ppk, pyrD, and serC) in 204 isolates collected in 1991, 1996, and 2001 from male and female patients of two public sexually transmitted diseases clinics. Genetic diversity (), recombination (C), growth (g), population structure, and adaptive selection under codon-substitution and amino acid property models were estimated and compared between these two gene classes. Estimates of the F(ST) fixation index and the chi(2) test of sequence absolute frequencies revealed significant temporal substructuring for both gene types. Baltimore's N. gonorrhoeae populations have increased since 1991 as indicated by consistent positive values of g. Female patients showed similar or lower levels of and C than male patients. Within the MLST housekeeping genes, levels of and C ranged from 0.001-0.013 and 0.000-0.018, respectively. Overall recombination seems to be the dominant force driving evolution in these populations. All loci showed amino acid sites and physicochemical properties under adaptive (or positive-destabilizing) selection, rejecting the generally assumed hypothesis of stabilizing selection for these MLST genes. Within the porB gene, protein I B showed higher and C values than protein I A. Directional positive selection possibly mediated by the immune system operates to a significant extent in the protein I sequences, as indicated by the distribution of the positively selected sites in the surface-exposed loops. Thirteen amino acid physicochemical properties seem to drive protein evolution of the PI porins in N. gonorrhoeae.