Skip to main content
  • ASM Journals
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems
  • Log in
  • My alerts
  • My Cart

Main menu

  • Home
  • Articles
    • Latest Articles
    • COVID-19 Research and News from ASM Journals
    • mSphere of Influence: Commentaries from Early Career Microbiologists
    • Archive
  • Topics
    • Applied and Environmental Science
    • Clinical Science and Epidemiology
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Therapeutics and Prevention
  • For Authors
    • Getting Started
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About mSphere
    • Editor in Chief
    • Board of Editors
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
  • ASM Journals
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems

User menu

  • Log in
  • My alerts
  • My Cart

Search

  • Advanced search
mSphere
publisher-logosite-logo

Advanced Search

  • Home
  • Articles
    • Latest Articles
    • COVID-19 Research and News from ASM Journals
    • mSphere of Influence: Commentaries from Early Career Microbiologists
    • Archive
  • Topics
    • Applied and Environmental Science
    • Clinical Science and Epidemiology
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Therapeutics and Prevention
  • For Authors
    • Getting Started
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About mSphere
    • Editor in Chief
    • Board of Editors
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
Research Article | Host-Microbe Biology

Transcriptome Analysis of Neisseria gonorrhoeae during Natural Infection Reveals Differential Expression of Antibiotic Resistance Determinants between Men and Women

Kathleen Nudel, Ryan McClure, Matthew Moreau, Emma Briars, A. Jeanine Abrams, Brian Tjaden, Xiao-Hong Su, David Trees, Peter A. Rice, Paola Massari, Caroline A. Genco
Sarah E. F. D'Orazio, Editor
Kathleen Nudel
aDepartment of Immunology, Tufts University School of Medicine, Boston, Massachusetts, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Kathleen Nudel
Ryan McClure
bBiological Sciences Division, Pacific Northwest National Laboratory, Richland, Washington, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Matthew Moreau
aDepartment of Immunology, Tufts University School of Medicine, Boston, Massachusetts, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Emma Briars
cProgram in Bioinformatics, Boston University, Boston, Massachusetts, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
A. Jeanine Abrams
dCenters for Disease Control and Prevention, Atlanta, Georgia, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Brian Tjaden
eDepartment of Computer Science, Wellesley College, Wellesley, Massachusetts, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Xiao-Hong Su
fNanjing STD Clinic, Institute of Dermatology, Chinese Academy of Medical Sciences and Peking Union Medical College, Nanjing, China
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
David Trees
dCenters for Disease Control and Prevention, Atlanta, Georgia, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Peter A. Rice
gDepartment of Medicine, Division of Infectious Diseases and Immunology, University of Massachusetts Medical School, Worcester, Massachusetts, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Paola Massari
aDepartment of Immunology, Tufts University School of Medicine, Boston, Massachusetts, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Caroline A. Genco
aDepartment of Immunology, Tufts University School of Medicine, Boston, Massachusetts, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Sarah E. F. D'Orazio
University of Kentucky
Roles: Editor
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Alison Criss
University of Virginia
Roles: Solicited external reviewer
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Joseph Dillard
University of Wisconsin—Madison
Roles: Solicited external reviewer
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
DOI: 10.1128/mSphereDirect.00312-18
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

ABSTRACT

Neisseria gonorrhoeae is a bacterial pathogen responsible for the sexually transmitted infection gonorrhea. Emergence of antimicrobial resistance (AMR) of N. gonorrhoeae worldwide has resulted in limited therapeutic choices for this infection. Men who seek treatment often have symptomatic urethritis; in contrast, gonococcal cervicitis in women is usually minimally symptomatic, but may progress to pelvic inflammatory disease. Previously, we reported the first analysis of gonococcal transcriptome expression determined in secretions from women with cervical infection. Here, we defined gonococcal global transcriptional responses in urethral specimens from men with symptomatic urethritis and compared these with transcriptional responses in specimens obtained from women with cervical infections and in vitro-grown N. gonorrhoeae isolates. This is the first comprehensive comparison of gonococcal gene expression in infected men and women. RNA sequencing analysis revealed that 9.4% of gonococcal genes showed increased expression exclusively in men and included genes involved in host immune cell interactions, while 4.3% showed increased expression exclusively in women and included phage-associated genes. Infected men and women displayed comparable antibiotic-resistant genotypes and in vitro phenotypes, but a 4-fold higher expression of the Mtr efflux pump-related genes was observed in men. These results suggest that expression of AMR genes is programed genotypically and also driven by sex-specific environments. Collectively, our results indicate that distinct N. gonorrhoeae gene expression signatures are detected during genital infection in men and women. We propose that therapeutic strategies could target sex-specific differences in expression of antibiotic resistance genes.

IMPORTANCE Recent emergence of antimicrobial resistance of Neisseria gonorrhoeae worldwide has resulted in limited therapeutic choices for treatment of infections caused by this organism. We performed global transcriptomic analysis of N. gonorrhoeae in subjects with gonorrhea who attended a Nanjing, China, sexually transmitted infection (STI) clinic, where antimicrobial resistance of N. gonorrhoeae is high and increasing. We found that N. gonorrhoeae transcriptional responses to infection differed in genital specimens taken from men and women, particularly antibiotic resistance gene expression, which was increased in men. These sex-specific findings may provide a new approach to guide therapeutic interventions and preventive measures that are also sex specific while providing additional insight to address antimicrobial resistance of N. gonorrhoeae.

INTRODUCTION

Neisseria gonorrhoeae is responsible for the sexually transmitted infection (STI) gonorrhea. Rates of gonococcal infection in the United States have increased 46% since 2011, and approximately 450,000 cases of gonorrhea were reported to the U.S. Centers for Disease Control and Prevention in 2015. N. gonorrhoeae is the second most common reportable bacterial STI in the United States (1). Infection with N. gonorrhoeae also contributes significantly to global STI morbidity and is responsible for 78 million cases each year (2). For the past 8 decades, gonorrhea has been treated successfully with antibiotics, despite patterns of resistance, but over the last several years, multidrug-resistant and untreatable strains have begun to emerge worldwide. Resistant strains include isolates from both Eastern and Western Europe (3, 4), Asia (5–7), and recently the United States (8–10). This rise in antibiotic-resistant strains may be linked, in part, to an increase in availability of over-the-counter antibiotics, particularly in many parts of Asia. New strategies for combating this disease are necessary, as evidenced by increased efforts to develop gonococcal vaccines (11, 12) and new drug regimens (7, 13, 14).

In humans, the only natural host for N. gonorrhoeae, symptomatic responses to infection are unique to men and women. While most infected men remain asymptomatic (15), those who develop symptoms often show robust inflammation characterized by purulent urethral discharge accompanied by large numbers of polymorphonuclear leukocytes (PMNs). In women, gonococcal infections are asymptomatic 50 to 80% of the time (16) or are accompanied by nonspecific symptoms such as vaginal discharge (17). The presence of cervical mucus or microscopic PMN counts of ≥10 seen microscopically on an oil immersion field does not correlate well with gonococcal infection in the absence of coinfecting agents (18). Manifestations of gonococcal infection may be driven by environmental factors specifically associated with the genital tracts of men and women: these may include biofilm formation (19, 20); the influence of the microbiome, particularly in women (21–25); and both direct and indirect molecular interactions between gonococci and specific host cells from men or women (26–30).

Gonococcal pathogenesis and host immune responses have been studied principally in vitro and in animal models (31), both of which may not faithfully replicate differences in infection between men and women. Historically, the only human model of gonococcal infection is based on genital experimental challenge of male volunteers, which has provided valuable information on gonococcal virulence factors, but is limited by the duration of the experimental infection and is not applicable to women (32, 33). A better understanding of gonococcal pathogenesis during natural mucosal infection, which is less well studied, is critical to implement biological strategies to treat and prevent infection. Previous work from our group using a gene-specific microarray and reverse transcription-quantitative PCR (qRT-PCR) analysis demonstrated that at least 20 gonococcal iron-regulated genes were expressed during human mucosal infection in men and women (34, 35). Recently, we also examined the complete N. gonorrhoeae transcriptome using cervicovaginal lavage specimens of infected women attending an STI clinic in China where antibiotic-resistant N. gonorrhoeae is prevalent (36, 37). We reported that a large portion of the gonococcal genome is expressed (65% of protein-coding genes) and regulated during infection (305 genes regulated compared to in vitro growth) (38). In the present study, we examined gonococcal transcriptomes from men attending the same clinic and compared bacterial gene expression and regulation in infected men with expression and regulation in infected women. The resulting data revealed sex-specific adaptation of genes involved in host-pathogen interactions and phage and antibiotic resistance.

RESULTS

Description of subjects.The male cohort comprised six male subjects (average age, 34 years) attending the Nanjing (China) Sexually Transmitted Disease (STD) Clinic, who presented with urethral discharge and/or dysuria and were diagnosed with N. gonorrhoeae infection (Table 1). Four men were coinfected with Chlamydia trachomatis and/or Mycoplasma genitalium, three had a history of prior gonococcal infection, and three had self-administered antibiotics prior to their clinic visit (Table 1). In contrast, none of the seven enrolled women had a history of gonococcal infection or had taken antibiotics prior to their clinic visit (see Table S1 in the supplemental material). Two gonococcal strains were isolated from a man and a woman who comprised a single dyad (subjects M-4 and F-6, respectively).

TABLE S1 

Characteristics of female subjects. Download TABLE S1, DOCX file, 0.02 MB.
Copyright © 2018 Nudel et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

View this table:
  • View inline
  • View popup
TABLE 1 

Characteristics of male subjects with gonococcal urethritis

Gonococcal gene expression in the human genital tract.RNA was extracted from male urethral or cervicovaginal lavage specimens, and individual cDNA reads were aligned to N. gonorrhoeae strain FA1090 using Rockhopper (39). Total gene expression was reported in reads per kilobase per million (RPKM). On average, 8.1% of the total RNA from male urethral specimens aligned to the N. gonorrhoeae FA1090 genome and 42.5% to the human genome; 4.3% of total RNA from cervicovaginal lavages aligned to FA1090 and 45% to the human genome.

Expression of 1,725 gonococcal genes (approximately 65% of the gonococcal genome) was detected in male urethral specimens. High expression of a core set of genes that encoded oxidative stress products (e.g., peroxiredoxin/glutaredoxin), housekeeping genes (e.g., the gene coding for elongation factor Tu), outer membrane protein genes (e.g., genes coding for Rmp [Omp3] and porin P.IB), and hypothetical protein genes were identified across male specimens (see Table S3 in Data Set S1 in the supplemental material). A more detailed analysis of the top 200 gonococcal genes (approximately 10% of total genome) that were expressed in all six male urethral secretions (in aggregate) was carried out by determining the ratio of gene enrichment and segregation into broad categories. Despite some variability in RPKM values among subjects, significantly enriched categories included genes coding for pilin and oxidative stress and regulatory proteins (Fig. 1), which reflected a common pattern of high gene expression during natural infection.

DATA SET S1

Supplemental Tables S3 to S7. Table S3 shows the top 10 most highly expressed genes in each male subject. Table S4 contains the composite RPKM values for each gene in men in vivo versus in vitro. Table S5 contains the composite RPKM values for significantly different noncoding RNA genes in men in vivo versus in vitro. Table S6 contains the composite RPKM values for significantly different predictive sRNAs in men in vivo versus in vitro and men in vivo versus women in vivo. Table S7 contains the composite RPKM values for each gene in men versus women. Download DATA SET S1, XLSX file, 0.2 MB.
Copyright © 2018 Nudel et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

FIG 1 
  • Open in new tab
  • Download powerpoint
FIG 1 

Categorization of the top 200 N. gonorrhoeae genes expressed during natural infection in men. Expression of gonococcal genes in 6 infected men was averaged, and the top 200 genes were categorized based on characteristics indicated for N. gonorrhoeae strain FA1090 available in NCBI. Categories are shown on the y axis, and the ratio of enrichment (% of genes of a given functional category in RNA-seq data set/% of genes assigned to that functional category in whole gonococcal genome) is shown on the x axis. For clarity, ribosomal protein, tRNA, and rRNA genes were removed. Asterisks indicate that enrichment was significant (P ≤ 0.05) by Fisher’s exact test.

Gonococcal gene regulation in the male genital tract.Our earlier studies had demonstrated that approximately 17% of gonococcal genes were regulated in vivo during cervical infection in women (38). Gonococcal gene expression in the male genital tract was compared to expression of matched infecting strains grown in vitro by analysis of transcriptome sequencing (RNA-seq) data sets for individual strains and for the composite data sets. Genes that displayed a significant difference in expression under in vivo versus in vitro conditions were selected for further analysis based on the following criteria: a ≥2-fold difference in RPKM value, a q (false-discovery rate) value of ≤0.05, and expression greater than 10 RPKM under at least one condition (see Table S4 in Data Set S1). The majority of genes were expressed similarly in vivo and in vitro, but approximately 4.2% of the total number of genes detected (79 genes) revealed increased expression in vivo compared to that observed in vitro; expression of 5.2% (99 genes) was decreased in vivo compared to that observed in vitro (Fig. 2A).

FIG 2 
  • Open in new tab
  • Download powerpoint
FIG 2 

Comparison of N. gonorrhoeae gene expression in infected men in vivo and isolates grown in vitro. (A) Composite gene expression levels in urethral specimens (n = 6) on the y axis plotted against composite expression of corresponding gonococcal isolates (n = 6) grown in vitro (x axis). Data are shown as log10 expression (RPKM) levels. The genes indicated in light blue had q values of ≤0.05 and ≥2-fold changes in vivo (versus in vitro). (B) Functional enrichment of genes with statistically significant changes in expression under in vivo and in vitro conditions. Genes were differentially expressed if the q value was ≤0.05 and the fold change was ≥2. Black bars represent functional enrichment of genes with increased expression in vivo. Gray bars represent functional enrichment of genes with decreased expression in vivo. Categories are shown on the y axis, and the ratio of enrichment (% of genes of a given functional category in RNA-seq data set/% of genes assigned to that functional category in whole gonococcal genome) is shown on the x axis. Ribosomal protein, tRNA, and rRNA genes were removed for clarity. *, significant enrichment among genes increased in vivo (versus in vitro), and #, significant enrichment among genes decreased in vivo, with P ≤ 0.05, by Fisher’s exact test. (C and D) Expression levels of iron genes (C) and oxidative stress genes (D) under in vivo and in vitro conditions. (Six male specimens [in vivo] and the corresponding 6 isolates [in vitro] are shown in the same color in adjacent columns.) *, q ≤ 0.05.

Regulation of gonococcal gene expression was then examined by functional enrichment analysis (Fig. 2B). Genes involved in iron and oxidative stress pathways and genes encoding hypothetical proteins were among those significantly enriched and were more highly expressed in vivo than in vitro. For example, expression of certain iron-scavenging genes (e.g., fbpA, tbpA, and tbpB) and of the fur-controlled regulator arsR was significantly higher in vivo than in vitro, although variance in arsR in vivo expression was large (Fig. 2C). This observation suggests that the male genital tract is iron depleted (similar to the female genital tract [38]). Genes involved in energy metabolism were also significantly enriched but had decreased expression in vivo (Fig. 2B). Among the downregulated energy metabolism genes was the nuoA–N operon (NGO1737–1751) (40), involved in aerobic respiration; its overall decreased expression in vivo may suggest a mechanism for bacterial survival in an anaerobic environment (41).

Among the enriched oxidative stress genes increased in vivo, we observed a significant increase in a glutaredoxin family protein gene (NGO0031), the thioredoxin I (trx1) gene (NGO0652), and msrAB (NGO2059) (Fig. 2D). The family of glutaredoxin proteins function as glutathione-dependent reductases (42); msrAB plays a role in reducing oxidized methionines to reactive peptides (43) and trx1 has been shown to respond to oxidation and is under control of the Fur regulon (44, 45). While expression of the well-studied OxyR regulon (oxyR, prx, kat, and gor) was not significantly different in in vivo samples compared to cultures grown in vitro, we observed decreased expression of the OxyR repressor in vivo, with increased expression of genes within the OxyR regulon (prx, kat, and gor). This notable trend further suggests that the gonococcus is exposed to reactive oxygen species during mucosal infection in men (Fig. 2D) (46, 47).

Several gonococcal noncoding RNAs were also regulated during infection in men. These included 4 rRNAs and 29 tRNAs, which were expressed at higher levels during infection compared to growth in vitro, except for NGO t31, an arginine tRNA that was expressed more highly during in vitro growth (see Table S5 in Data Set S1). Putative small RNAs (sRNAs: defined as expression representing an intergenic region of the complementary strand of a known open reading frame [ORF] and 30 to 250 nucleotides in length) were also regulated during infection; 29 sRNAs were increased in vivo compared to in vitro, and 21 sRNAs were decreased in vivo (see Table S6 in Data Set S1). Among these sRNAs, several have been described previously (48, 49).

Comparison of gonococcal gene signatures in the male and female genital tracts.Building on gonococcal transcriptome data in infected men and women, we examined whether site-specific, environmental conditions influenced expression of gonococcal genes during infection. A comparative transcriptome analysis was carried out using the 6 urethral and 7 cervicovaginal specimens. We performed principal-component analysis (PCA) to assess whether global in vivo gene expression revealed transcriptional signatures specific for infections in men or women (Fig. 3A). PCA showed a clear distinction between gonococcal genes expressed in specimens from men and women, with high variance among the 7 specimens from women compared to low variance in men, as evidenced by tight clustering when the first two principal components were plotted against each other (Fig. 3A). Gene expression analysis of the composite data sets was also performed; parameters set for defining differential gene regulation were a ≥2-fold change in RPKM, a q value of ≤0.05, and expression greater than 10 RPKM under at least one condition (see Table S5 in Data Set S1). Of 1,879 genes expressed, 176 (9.4%) were increased in men compared to women and 80 (4.3%) were decreased in men, indicating that nearly 14% of gonococcal genes were differentially regulated during gonococcal infection in the male and female genital tracts (see Table S7 in Data Set S1).

FIG 3 
  • Open in new tab
  • Download powerpoint
FIG 3 

Comparison of N. gonorrhoeae gene expression during infection in men and women in vivo. (A) Variance in global expression in each specimen from men and women. The first principal component (PC1 [x axis]) had a variance of ~26%, and the second principal component (PC2 [y axis]) had a variance of ~14%. (B) Functional enrichment of genes with significant changes in expression levels in men and women. Genes were differentially expressed if the q value was ≤0.05 and the fold change was ≥2. Blue bars represent genes with increased expression in vivo in men compared to women. Red bars represent genes with increased expression in vivo in women compared to men. *, significant enrichment among genes increased in vivo in men; #, significant enrichment among genes increased in vivo in women.

A functional enrichment analysis of gene categories performed on male and female specimens revealed that phage-associated genes were significantly enriched in female specimens (Fig. 3B). While female specimens had higher expression of a broad range of phage types (e.g., double-stranded and filamentous single-stranded phages), male specimens had higher expression in the subset of double-stranded DNA prophages (see Fig. S1A in the supplemental material) (50). Other highly expressed categories enriched in male specimens included genes that encoded hypothetical proteins and proteins involved in host interactions (Fig. 3B). Iron and oxidative stress genes that were differentially expressed in male and female specimens included the ferredoxin gene (NGO1859), a TonB receptor protein gene (NGO1205), and fetA (NGO2093). (fetA is regulated by MpeR and induced under low-iron conditions [51].) These genes all manifested increased expression in women (Fig. S1B). In contrast, expression of the oxidative stress gene bfrB was increased in male specimens; other oxidative stress genes were similarly expressed in both men and women (Fig. S1C), suggesting that exposure to oxidative stress occurs in both the male and female genital tracts.

FIG S1 

Phage expression, iron regulation, and oxidative response genes in the male and female genital tracts. (A) Expression levels of phage genes are shown on the x axis and represent male urethral samples (blue) and female cervicovaginal lavage samples (red). (B and C) Expression levels of significantly differentially expressed iron genes (B) and oxidative stress genes (C). *, gene more highly expressed in the male genital tract; #, gene more highly expressed in the female genital tract (q value of <0.05). Red and blue X’s indicate the male and female counterparts of the dyad. Download FIG S1, TIF file, 0.8 MB.
Copyright © 2018 Nudel et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Regulation of several gonococcal sRNAs was also observed. A significant change in expression of 119 sRNAs was observed in specimens obtained from men and women: 28 were expressed at higher levels in male specimens and 91 in female specimens. Previously identified sRNAs showing regulation based on sex included smRNA4, smRNA5, and smRNA8 (49), fnrS (a small RNA induced under anaerobic conditions [44]), and the iron-repressed sRNA encoding NrrF (52) (see Table S6 in Data Set S1). In addition, 45 tRNAs showed differential regulation, with all being more highly expressed in female patients than in male patients (see Table S5 in Data Set S1).

Expression of antimicrobial resistance genes.Gonococcal antimicrobial resistance (AMR) is a growing concern and is common in strains from the eastern regions of the world, including China. We examined phenotypic and genotypic AMR patterns from our cohort and found that all the isolates exhibited resistance to penicillin, tetracycline, and ciprofloxacin, with one isolate showing resistance to azithromycin (Table 2). One isolate also showed reduced susceptibility to ceftriaxone. Strains isolated from women exhibited similar antimicrobial susceptibility patterns (Table 2). These results were similar to AMR patterns of male urethritis strains reported earlier (2011 to 2012) from the Nanjing STD Clinic (53).

View this table:
  • View inline
  • View popup
  • Download powerpoint
TABLE 2 

Phenotypic resistance and AMR genes identified in male and female N. gonorrhoeae isolates (n = 13) by WGSa

Whole-genome DNA sequencing (WGS) was performed on all available strains (6 male and 6 female) to define strain relatedness and genotypic determinants of antibiotic resistance. A single-nucleotide polymorphism (SNP) distance tree was constructed that used SNP distances (differences) that ranged between 81 and 3,913 separating the genomes (see Fig. S2 in the supplemental material). The male and female dyad (M-4 and F-6) yielded genomic sequences that were separated by only 81 SNPs. Strain NCCP11945 isolated from a subject with gonococcal disease in South Korea (average SNP distance of 2,828 SNPs compared to strain FA1090) shared considerable homology with the isolates obtained from our cohort compared to the commonly used laboratory strain, FA1090 (average distance of 4,212 compared to FA1090).

FIG S2 

Hamming distance tree and SNP distance matrix of N. gonorrhoeae isolates. Tree visualizing Hamming distances of total pairwise SNPs identified across strains. The scale bar shows a distance of 300 SNPs. Download FIG S2, TIF file, 0.9 MB.
Copyright © 2018 Nudel et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

Genomic analyses based on publicly available, curated, antibiotic resistance databases (CARD and GC-MLST) identified multiple AMR genes (Table 2) and SNP-level variants that likely contributed to phenotypic resistance. All gonococcal strains exhibited SNP mutations in gyrA and/or parC (which confer resistance to fluoroquinolones), SNP mutations in penA, ponA, porB, and pilQ (β-lactam resistance), and SNP mutations in rpsJ (tetracycline resistance) (54). Only one isolate (M-5) carried plasmid-based resistance genes (tem-1 and tetM). Interestingly, among the male subjects, M-3 was the only isolate to exhibit azithromycin resistance and a corresponding mutation in the 23S rRNA gene (55). Three isolates had an adenosine deletion in the promoter region of mtrR, the repressor of the mtrCDE efflux pump, a system that exports bactericidal host-derived compounds and hydrophobic antibiotics (β-lactams and macrolides) from the bacterial cell (56–58). Isolates from women exhibited similar genomic patterns of resistance: by and large, the same mutations identified in the male isolates were also present in the female isolates (Table 2). Among 6 SNPs associated with β-lactam resistance in the dyad isolates (subjects M-4 and F-6, who were infected by the same strain), 5 SNPs were shared in the dyad isolates. Two SNPs in gyrA were shared (S91F and D95N), and a Ser SNP in parC [parC(S87I)] was detected only in the male isolate. Finally, while several SNP mutations coincided with phenotypic resistance to azithromycin, ciprofloxacin, and penicillin (Table 2), the single isolate with decreased sensitivity to ceftriaxone (M-3) did not possess the mosaic penA XXXIV allele (59) or any additional curated mutations that distinguished it from the phenotypically sensitive strains (Table 2) (59, 60).

Expression levels of genes associated with antibiotic resistance in men were also compared to expression of corresponding genes in vitro (Fig. 4A and B). Expression of the pilQ gene, associated with antibiotic permeability (14), was significantly decreased by 1.78-fold in vivo compared to in vitro (Fig. 4A). The MisRS two-component regulatory system is important for resistance to host-derived antimicrobial peptides and antibiotics (61), and though not statistically significant, misR had 2.08-fold higher expression in vivo compared to in vitro, supporting its important role in vivo. The expression of the Mtr efflux pump was also closely examined; expression of mtrA, which can influence expression of mtrCDE (62), was significantly decreased by 2.75-fold in vivo (Fig. 4B). Gonococcal strains from male urethral specimens (M-2, M-3, M-4, and M-6) contained an adenosine deletion in the mtrR promoter that was accompanied by the lowest expression of mtrR and highest expression of mtrCDE and mtrF, confirming that the adenosine deletion directly affects expression of the Mtr locus (63) in vivo. Expression of mpeR was also investigated, as it has been shown to regulate the Mtr locus; however, expression was below 10 RPKM in all in vitro and in vivo samples.

FIG 4 
  • Open in new tab
  • Download powerpoint
FIG 4 

Expression of gonococcal antibiotic resistance genes. Expression levels of antibiotic resistance genes (A) and the Mtr efflux pump (B) in the male genital tract (in vivo) and the corresponding strain in CDM (in vitro); genes are identified on the x axis. (Six male specimens [in vivo] and the corresponding 6 isolates [in vitro] are shown in the same color.) *, q ≤ 0.05. (C) Expression levels of antibiotic resistance genes from the male genital tract (blue) or the female genital tract (red). *, higher gene expression in men; #, higher gene expression in women (q ≤ 0.05). Red and blue X’s indicate the male (M-4) and female (F-6) counterparts of the dyad.

Transcriptomes from male urethral specimens were also compared with female cervicovaginal lavage specimens to determine if expression of AMR genes was influenced by sex (Fig. 4C). Because of the adenosine deletion in the mtrR promoter in 4/6 males (described above) and all female isolates (Table 2), expression of mtrR was low in both groups and resulted in significantly higher expression of mtrCDE and mtrF in men, as predicted (63), but not in women. Even within the dyad, where both M-4 and F-6 contained the adenosine deletion in the mtrR promoter, M-4 had approximately 2× higher expression levels of mtrCDE than F-6. Interestingly, the expression of mtrA, an activator of mtrCDE (62), was significantly higher during infection in women compared to men, despite lower expression of mtrCDE in women (Fig. 4C). The farR gene, which negatively controls expression of the farAB operon by directly binding to the farAB promoter (64) to prevent excess expression of gonococcal efflux pumps, also showed significantly higher expression in cervicovaginal lavages. Collectively these results suggest that, although strains isolated from infected men and women exhibit comparable AMR genotypes and in vitro phenotypes, expression of AMR genes is different during infection in the male and female genital tracts.

DISCUSSION

N. gonorrhoeae infects the male and female genital tracts, two very distinct environments in humans. As there is little doubt about the intrinsic tissue, cellular, and molecular differences that define these host environments, it is reasonable to assume that the gonococcus will adapt to these environmental differences during infection (65). Building on our previous results on the gonococcal transcriptome profile during infection in women (37), we extended our analysis to specimens from infected men. In the present study, we not only report a distinct gonococcal transcriptome profile in vivo compared to in vitro in infected men, but also, for the first time, describe a comparison of N. gonorrhoeae gene expression during infection in the male and female genital tracts. Our results also demonstrate distinct gonococcal gene expression signatures in men and women, consistent with the intrinsically different makeups and natures of the two sites of infection. Furthermore, our approach based on both whole-genome and RNA sequencing enabled us to evaluate expression of antibiotic resistance determinants during infection in both the male and female genital tracts.

Currently, the only human model of gonococcal infection is experimental urethral challenge of male volunteers. Those studies have examined gonococcal virulence factors and have provided valuable information on opacity proteins and genes involved in iron acquisition, structural modification of lipooligosaccharides, and a variety of other established virulence factors (32, 66, 67) but are limited by the short duration of experimental infection. The small number of infecting strains used experimentally in men (mostly laboratory-adapted strains FA1090 and MS11mkC [32]) have precluded an in-depth understanding of the range of adaptions employed during gonococcal infections generally. Using a panel of N. gonorrhoeae strains isolated from infected subjects, we compared the gonococcal in vivo transcriptome expressed during infection of the male genital tract to that of the corresponding infecting strains grown in vitro. Our results revealed increased expression of genes involved in oxidative stress and iron scavenging in vivo, accompanied by decreased expression of genes involved in metabolism-associated categories (i.e., energy and amino acid metabolism). In particular, we observed that gonococci expressed high levels of trx1, msrAB, and glutaredoxin genes during infection in male subjects, in agreement with the large PMN influx that occurs during infection in symptomatic men. Our results also demonstrated increased expression of iron-regulated genes (tbpAB [66]), indicating that the male genital tract is an iron-depleted environment. We observed lower expression of genes involved in gonococcal aerobic energy metabolism processes and higher expression of oxidative stress response genes in urethral specimens compared to that of the corresponding infecting strains grown in vitro. It is established that the genital tract is anaerobic or microaerobic (68, 69). Thus, we would expect that a comparison of RNA-seq data sets from urethral and cervicovaginal lavage specimens to RNA-seq data sets from in vitro growth conditions would identify differential regulation of gonococcal genes utilized for aerobic energy production. However, the presence of neutrophils, especially in male subjects during infection, likely requires a corresponding N. gonorrhoeae oxidative defense response. Such immune cells are obviously absent under our in vitro conditions, and thus these oxidative stress pathways are expressed at lower levels. While we recognize that in vitro culture conditions differ compared to the conditions in the human genital tract beyond metabolite content, this study provides a gene-specific view of the response of N. gonorrheoae to infection and leads to a number of hypotheses regarding infection mechanisms that can be explored in more detail in future studies.

While we identified both expected and novel differences in gonococcal gene expression during human infection compared to in vitro growth, the most revealing aspect of our study was the differential gonococcal expression profiles observed in men and women. Our results showed that nearly 14% of gonococcal genes that were expressed in human infection were differentially expressed in men and women. Principal-component analysis demonstrated that gonococcal expression in men and women separated into distinct gene expression signatures, with high variance among the 7 female specimens and low variance among the 6 male specimens. We acknowledge that high variance could be due to the small sample size, but it may also be due to other female-specific environmental factors that could alter gonococcal gene expression, such as menstrual cycle or the microbiome, which exhibits diversity among individuals (70). We also observed differences in the presence and type of coinfecting organisms in infected men and women. Several male subjects had a coinfection of Chlamydia trachomatis or Mycobacterium genitalium, and while some female subjects were also infected with these pathogens, several females additionally carried Ureaplasma urealyticum. In addition, 3/7 of the female subjects had bacterial vaginosis (BV), suggesting differences in their natural microflora in combination with these coinfections. We hypothesize that changes in the microflora and the presence of coinfecting organisms have a direct effect on the transcriptomic response of N. gonorrhoeae, as does the variable host response during coinfection compared to infection with N. gonorrhoeae alone. Both of these factors most likely impact the overall differences observed between samples of the same gender as well as differences observed in the gonococcal transcriptomic response during infection in men and women. When genomic DNA sequence analysis was also considered, there was no clustering among strains derived from men or women (except for strains in the single dyad), which strongly suggests that the observed differences in transcriptomic signatures were driven by exposure of gonococci, generally, to the male or the female genital tract and not by strain relatedness. Furthermore, this is also supported in the single dyad by differences in gene expression manifested by the same strain in each partner, which reflected changes in gene expression that were seen overall in the other men and women. In addition to the small sample size, it is important to also consider other limitations of the study when interpreting sex-specific and growth-specific aspects of N. gonorrhoeae infection. Due to the design of the study, men are recruited because they enter the clinic with a symptomatic response: thus, men may have a similar stage of infection and gonococcal growth stage. This is in contrast to the matched females, who are partners of the recruited men and who were at various stages of infection and therefore gonococcal growth.

Among the most relevant gonococcal transcriptome differences that emerged in male and female specimens were those in genes encoding antimicrobial resistance determinants. The rise of antibiotic resistance in N. gonorrhoeae is now a global concern (71). Patterns of resistance vary among countries because of population mobility and geographically localized increases in antimicrobial resistance—such as those observed in China (72); these patterns are predictive of future resistance patterns in countries that currently have lower levels of resistance, such as the United States (73, 74). Whole-genome DNA sequencing (WGS) of the gonococcal strains isolated from our male and female cohort revealed similarity of genetic determinants of resistance to strains isolated in the United States (59). In all strains, regardless of the origin of the specimens (men or women), the presence of an adenosine deletion in the mtrR promoter (also detected by WGS) correlated with increased expression of the Mtr regulon. Nevertheless, strains isolated from men showed a significantly higher expression of mtrCDE and mtrF than those from women. While a direct explanation for these results remains speculative, it is possible that the higher expression of mtrCDE observed in men may be due to the readily accessible antibiotics in China; approximately half of the male subjects (but none of the female subjects) in this study had self-administered antibiotics prior to seeking care at the STD clinic. Future studies including a larger cohort may help resolve these differences. An additional explanation for gender-specific differences is that biofilm formation in the female genital tract (and apparent lack of it in the male genital tract) is a contributing factor (19, 20). Biofilms can provide protection against harsh environments, and their absence in the male genital tract may lead to increased expression of antibiotic-resistant efflux pumps as a defense mechanism (75, 76).

N. gonorrhoeae contains phage DNA inserted in the genome: five regions comprising double-stranded DNA (dsDNA) lysogenic phage genomes (50) and four additional regions comprising filamentous phages (77). The role of phage genes in N. gonorrhoeae infection has not been studied extensively. Previous work from our group has identified a phage gene, npr, involved in invasion of female epithelial cells in vitro, which correlates with disease progression in a female mouse model of gonococcal infection (78). Our finding that expression of phage genes was increased in female subjects may also relate to the formation of biofilms in the female genital tract (20). We propose that biofilm formation may facilitate DNA exchange via phages among bacteria, thus contributing to differences in infections of men and women. Additional transcriptional patterns that differentiated organisms obtained during infection in men and women included genes encoding tRNAs with higher expression observed in women compared to men. However, apparent expression of tRNAs may be partially affected by reads aligning to additional microbes present in the genital tract: tRNA genes are highly conserved across species, and the women in the study may have a larger population of microbial organisms. Six of seven gonococcus-infected women reported in this study had BV and/or other bacterial genital infections. Some studies have suggested a role for tRNAs in bacterial infection (79, 80), and it is possible that differences observed in men and women also correlate with the enhanced microbial environment in women. Small RNA expression in vivo was also different; NrrF and FnrS (an sRNA induced under anaerobic conditions) were expressed at higher levels in women than in men.

Collectively, our results provide the first global view of gonococcal gene expression during infection in humans and define gene signatures specific to infections in men and women. We report important differences related to antibiotic resistance and gonococcal pathogenesis that can be extrapolated to improve understanding of gonococcal disease outcomes and potential treatments. Our analysis also highlights shortfalls of studying bacterial infections using in vitro models and systems. It is critical that studies designed to identify targeted therapies for gonococcal infections consider sex-specific differences in gene expression profiles that may impact treatment outcomes. Furthermore, addressing how expression of antimicrobial resistance genes are driven by environmental cues in the male and female genital tracts has important implications for the use of targeted antibiotics.

MATERIALS AND METHODS

Ethics statement.All subjects provided written informed consent in accordance with requirements by Institutional Review Boards from Tufts University, Boston, MA (protocol no. 11710), the University of Massachusetts Medical School, Worcester, MA, and Boston University School of Medicine, Boston, MA, and the Institute of Dermatology, Chinese Academy of Medical Sciences and Peking Union Medical College, Nanjing, China.

Identification of subjects and characterization of urethral infections.Urethral swab specimens were collected from six men attending the Nanjing Sexually Transmitted Disease (STD) Clinic, chosen at random, who were participants in a study of gonococcal transmission (81). Male subjects presented with symptoms and signs of urethritis (dysuria and/or urethral discharge). Gram stains of urethral exudates showed PMNs and Gram-negative intracellular diplococci, criteria that are highly specific (99%) for gonococcal infection (82). Urethral swab specimens from men with urethritis were also inoculated onto Thayer-Martin medium (DL Biotech, China) and cultured in candle jars at 36°C. Gonococcal isolates were identified by colonial morphology, Gram stain, and oxidase testing. Urethral swab specimens were also tested for Chlamydia trachomatis and Mycoplasma genitalium by PCR (83, 84) and Ureaplasma urealyticum and Mycoplasma hominis by culture using the Mycoplasma IST2 test (BioMérieux, France).

Collection of cervicovaginal lavage specimens (n = 4) was described previously (38), and three additional specimens were obtained for the present study (total n = 7). PCR testing using primers directed against the gonococcal porA pseudogene (85) had previously identified N. gonorrhoeae infection in one woman (F-4) whose cervical culture was negative. Women reported that they were monogamous (3/7 were married), having been the sole sex partner of a man with gonococcal urethritis. One woman (F-6) and one man (M-4) were members of a single dyad.

Antimicrobial testing of N. gonorrhoeae.N. gonorrhoeae strains from urethral swab and cervicovaginal lavage specimens were identified and processed as previously described (38). Antimicrobial sensitivity testing of N. gonorrhoeae isolates was determined for penicillin, tetracycline, spectinomycin, azithromycin, ceftriaxone, cefixime, and ciprofloxacin at the Nanjing STD Clinic. Mean inhibitory concentrations (MICs) were determined by the agar dilution method, and current Clinical and Laboratory Standards Institute cutoffs were used, except for azithromycin, in which we followed the European Committee on Antimicrobial Susceptibility Testing (http://www.eucast.org) (86). A single isolate (from M-2) that did not grow in agar dilution media or GC agar had MICs determined on chocolate agar by Etest: thus, MIC values may be one dilution lower due to enriched agar (87).

RNA isolation.Total RNA was isolated from urethral swab and cervicovaginal lavage specimens using TRIzol (Invitrogen) as described previously (38). Briefly, specimens were washed twice with 70% ethanol and treated with Turbo DNase (Ambion), the RiboMinus kit (Invitrogen) (to deplete eukaryotic rRNA), and the Microbe Express kit (Ambion) (to deplete bacterial rRNA) using diethylpyrocarbonate (DEPC)-treated EDTA-free water. Gonococcal isolates that corresponded to each specimen were grown overnight on chocolate agar plates at 37°C in 5% CO2 prior to inoculation in chemically defined medium (CDM) containing ferric nitrate (100 µM). Liquid cultures were inoculated at an optical density at 600 nm (OD600) of 0.1, harvested after 3 h, and pelleted, and RNA was extracted as described above.

RNA sequencing and analysis.RNA sequencing was performed as previously described (38). cDNA libraries were prepared with TruSeq RNA preparation kit and sequenced on an Illumina HiSeq 2500 instrument using high-output V3 chemistry in single-read 100 formats. Analysis of RNA-seq data was performed using the Rockhopper program, aligning the reads to the NCBI genome sequence of N. gonorrhoeae strain FA1090 and then reporting the RNA-seq data in reads per kilobase per million (RPKM) (37). The cutoff for gene expression was set at an RPKM value of at least 10. A composite set of RPKM values was derived from RNA-seq data from the six male urethral samples and another from the three newly collected cervicovaginal samples added to the four reported earlier (in the Gene Expression Omnibus database under GenBank accession no. GSE71151) (38). These are referred to as “in vivo” RNA-seq data. Data from gonococcal strains isolated from each subject and cultured in vitro were also obtained (referred to as “in vitro” data).

Statistical analysis was performed on pooled data from each category of specimens: in vivo urethral, in vivo cervicovaginal lavage, and in vitro. Genes with ≥2-fold change in expression and a q value of ≤0.05 in each category were deemed as regulated. Gene assignments and categorization were performed using the KEGG and NCBI RefSeq public databases. Functional enrichment ratios were determined by dividing the percentage of genes in a functional category for a particular set of genes (e.g., highly expressed genes or differentially regulated genes) by the percentage of genes in the same category present in the gonococcal genome as a whole. This ratio was used to interpret enrichment of a particular set of genes that possessed a specific function. Fisher’s exact test was used to assign a P value to the enrichment. Principal-component analysis (PCA) (88) was performed in R using singular value decomposition on the scaled and centered data set of RPKM values for all genes in each subject, and the resulting data were visualized by sex on a colored plot.

Whole-genome sequencing and analysis.N. gonorrhoeae isolates were grown on GC agar plates supplemented with 1% IsoVitaleX and 0.6% fetal bovine serum at 37°C and 5% CO2 for 16 to 18 h. Genomic DNA was extracted from a 10-µl inoculum of each bacterial culture using a 5Prime DNA extraction kit (5Prime, San Francisco, CA), following the manufacturer’s recommendations with slight modifications. Whole-genome sequencing (WGS) was conducted using a PacBio RSII platform (Pacific Biosciences, Menlo Park, CA) with P5-C3 chemistry. De novo genome assembly was conducted using the hierarchical genome assembly process workflow (HGAP3; SMRTAnalysis 2.3.0), which included consensus polishing using Quiver (89). Single-nucleotide polymorphisms (SNPs) across sequenced isolates were called de novo-assembled contigs using kSNP3 without a reference strain (90). Concatenated SNPs were used to construct phylogenetic trees using the approximate-maximum-likelihood-based approach; trees were visualized in FigTree (http://tree.bio.ed.ac.uk/software/figtree/). Antibiotic resistance elements were identified using the Comprehensive Antibiotic Resistance Database (91) and the PubMLST Neisseria database (https://pubmlst.org/neisseria) (92).

Data availability.New RNA-seq data from this study have been deposited in the Gene Expression Omnibus (GSE113290). Genomic data from isolates have been deposited into NCBI under project no. PRJNA329501 (see Table S2 in the supplemental material).

TABLE S2 

Deposited genomic data. Download TABLE S2, DOCX file, 0.01 MB.
Copyright © 2018 Nudel et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

ACKNOWLEDGMENTS

This work was supported by NIH/NIAID grants U19AI084048 (P. A. Rice and C. A. Genco) and R01AI116969 (C. A. Genco). PNNL is operated for the DOE by Battelle Memorial Institute under contract DE-AC05-76RLO 1830.

The conclusions, findings, and opinions expressed by the authors do not necessarily reflect the official position of the U.S. Centers for Disease Control and Prevention (CDC).

FOOTNOTES

    • Received June 6, 2018.
    • Accepted June 6, 2018.
  • This paper was submitted via the mSphereDirect™ pathway.

  • Copyright © 2018 Nudel et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

REFERENCES

  1. 1.↵
    1. Centers for Disease Control and Prevention
    . 2016. 2015 sexually transmitted disease surveillance. Centers for Disease Control and Prevention, Atlanta, GA.
  2. 2.↵
    1. Newman L,
    2. Rowley J,
    3. Vander Hoorn S,
    4. Wijesooriya NS,
    5. Unemo M,
    6. Low N,
    7. Stevens G,
    8. Gottlieb S,
    9. Kiarie J,
    10. Temmerman M
    . 2015. Global estimates of the prevalence and incidence of four curable sexually transmitted infections in 2012 based on systematic review and global reporting. PLoS One 10:e0143304. doi:10.1371/journal.pone.0143304.
    OpenUrlCrossRefPubMed
  3. 3.↵
    1. Brunner A,
    2. Nemes-Nikodem E,
    3. Jeney C,
    4. Szabo D,
    5. Marschalko M,
    6. Karpati S,
    7. Ostorhazi E
    . 2016. Emerging azithromycin-resistance among the Neisseria gonorrhoeae strains isolated in Hungary. Ann Clin Microbiol Antimicrob 15:53. doi:10.1186/s12941-016-0166-9.
    OpenUrlCrossRef
  4. 4.↵
    1. Unemo M,
    2. Golparian D,
    3. Nicholas R,
    4. Ohnishi M,
    5. Gallay A,
    6. Sednaoui P
    . 2012. 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 56:1273–1280. doi:10.1128/AAC.05760-11.
    OpenUrlAbstract/FREE Full Text
  5. 5.↵
    1. Cheng CW,
    2. Li LH,
    3. Su CY,
    4. Li SY,
    5. Yen MY
    . 2016. Changes in the six most common sequence types of Neisseria gonorrhoeae, including ST4378, identified by surveillance of antimicrobial resistance in northern Taiwan from 2006 to 2013. J Microbiol Immunol Infect 49:708–716. doi:10.1016/j.jmii.2014.08.016.
    OpenUrlCrossRef
  6. 6.↵
    1. Shimuta K,
    2. Unemo M,
    3. Nakayama S,
    4. Morita-Ishihara T,
    5. Dorin M,
    6. Kawahata T,
    7. Ohnishi M, Antibiotic-Resistant Gonorrhea Study Group
    . 2013. Antimicrobial resistance and molecular typing of Neisseria gonorrhoeae isolates in Kyoto and Osaka, Japan, 2010 to 2012: intensified surveillance after identification of the first strain (H041) with high-level ceftriaxone resistance. Antimicrob Agents Chemother 57:5225–5232. doi:10.1128/AAC.01295-13.
    OpenUrlAbstract/FREE Full Text
  7. 7.↵
    1. Su XH,
    2. Wang BX,
    3. Le WJ,
    4. Liu YR,
    5. Wan C,
    6. Li S,
    7. Alm RA,
    8. Mueller JP,
    9. Rice PA
    . 2016. Multidrug-resistant Neisseria gonorrhoeae isolates from Nanjing, China, are sensitive to killing by a novel DNA gyrase inhibitor, ETX0914 (AZD0914). Antimicrob Agents Chemother 60:621–623. doi:10.1128/AAC.01211-15.
    OpenUrlAbstract/FREE Full Text
  8. 8.↵
    1. Gose S,
    2. Nguyen D,
    3. Lowenberg D,
    4. Samuel M,
    5. Bauer H,
    6. Pandori M
    . 2013. Neisseria gonorrhoeae and extended-spectrum cephalosporins in California: surveillance and molecular detection of mosaic penA. BMC Infect Dis 13:570. doi:10.1186/1471-2334-13-570.
    OpenUrlCrossRefPubMed
  9. 9.↵
    1. Nguyen D,
    2. Gose S,
    3. Castro L,
    4. Chung K,
    5. Bernstein K,
    6. Samuel M,
    7. Bauer H,
    8. Pandori M
    . 2014. Neisseria gonorrhoeae strain with reduced susceptibilities to extended-spectrum cephalosporins. Emerg Infect Dis 20:1211–1213. doi:10.3201/eid2007.131396.
    OpenUrlCrossRefPubMed
  10. 10.↵
    1. Wu A,
    2. Buono S,
    3. Katz KA,
    4. Pandori MW
    . 2011. Clinical Neisseria gonorrhoeae isolates in the United States with resistance to azithromycin possess mutations in all 23S rRNA alleles and the mtrR coding region. Microb Drug Resist 17:425–427. doi:10.1089/mdr.2010.0199.
    OpenUrlCrossRefPubMedWeb of Science
  11. 11.↵
    1. Cash DR,
    2. Noinaj N,
    3. Buchanan SK,
    4. Cornelissen CN
    . 2015. Beyond the crystal structure: insight into the function and vaccine potential of TbpA expressed by Neisseria gonorrhoeae. Infect Immun 83:4438–4449. doi:10.1128/IAI.00762-15.
    OpenUrlAbstract/FREE Full Text
  12. 12.↵
    1. Piekarowicz A,
    2. Kłyż A,
    3. Majchrzak M,
    4. Stein DC
    . 2016. Oral immunization of rabbits with S. enterica typhimurium expressing Neisseria gonorrhoeae filamentous phage Phi6 induces bactericidal antibodies against N. gonorrhoeae. Sci Rep 6:22549. doi:10.1038/srep22549.
    OpenUrlCrossRef
  13. 13.↵
    1. Cybulska P,
    2. Thakur SD,
    3. Foster BC,
    4. Scott IM,
    5. Leduc RI,
    6. Arnason JT,
    7. Dillon JA
    . 2011. Extracts of Canadian First Nations medicinal plants, used as natural products, inhibit Neisseria gonorrhoeae isolates with different antibiotic resistance profiles. Sex Transm Dis 38:667–671. doi:10.1097/OLQ.0b013e31820cb166.
    OpenUrlCrossRefPubMed
  14. 14.↵
    1. Nandi S,
    2. Swanson S,
    3. Tomberg J,
    4. Nicholas RA
    . 2015. Diffusion of antibiotics through the PilQ secretin in Neisseria gonorrhoeae occurs through the immature, sodium dodecyl sulfate-labile form. J Bacteriol 197:1308–1321. doi:10.1128/JB.02628-14.
    OpenUrlAbstract/FREE Full Text
  15. 15.↵
    1. Detels R,
    2. Green AM,
    3. Klausner JD,
    4. Katzenstein D,
    5. Gaydos C,
    6. Handsfield H,
    7. Pequegnat W,
    8. Mayer K,
    9. Hartwell TD,
    10. Quinn TC
    . 2011. The incidence and correlates of symptomatic and asymptomatic Chlamydia trachomatis and Neisseria gonorrhoeae infections in selected populations in five countries. Sex Transm Dis 38:503–509.
    OpenUrlPubMed
  16. 16.↵
    1. Hook EW,
    2. Holmes KK
    . 1985. Gonococcal infections. Ann Intern Med 102:229–243. doi:10.7326/0003-4819-102-2-229.
    OpenUrlCrossRefPubMedWeb of Science
  17. 17.↵
    1. McCormack WM,
    2. Stumacher RJ,
    3. Johnson K,
    4. Donner A
    . 1977. Clinical spectrum of gonococcal infection in women. Lancet i:1182–1185. doi:10.1016/S0140-6736(77)92720-9.
    OpenUrlCrossRef
  18. 18.↵
    1. Brunham RC,
    2. Paavonen J,
    3. Stevens CE,
    4. Kiviat N,
    5. Kuo CC,
    6. Critchlow CW,
    7. Holmes KK
    . 1984. Mucopurulent cervicitis—the ignored counterpart in women of urethritis in men. N Engl J Med 311:1–6. doi:10.1056/NEJM198407053110101.
    OpenUrlCrossRefPubMedWeb of Science
  19. 19.↵
    1. Greiner LL,
    2. Edwards JL,
    3. Shao J,
    4. Rabinak C,
    5. Entz D,
    6. Apicella MA
    . 2005. Biofilm formation by Neisseria gonorrhoeae. Infect Immun 73:1964–1970. doi:10.1128/IAI.73.4.1964-1970.2005.
    OpenUrlAbstract/FREE Full Text
  20. 20.↵
    1. Steichen CT,
    2. Shao JQ,
    3. Ketterer MR,
    4. Apicella MA
    . 2008. Gonococcal cervicitis: a role for biofilm in pathogenesis. J Infect Dis 198:1856–1861. doi:10.1086/593336.
    OpenUrlCrossRefPubMed
  21. 21.↵
    1. Nikolaitchouk N,
    2. Andersch B,
    3. Falsen E,
    4. Strömbeck L,
    5. Mattsby-Baltzer I
    . 2008. The lower genital tract microbiota in relation to cytokine-, SLPI- and endotoxin levels: application of checkerboard DNA-DNA hybridization (CDH). APMIS 116:263–277. doi:10.1111/j.1600-0463.2008.00808.x.
    OpenUrlCrossRefPubMedWeb of Science
  22. 22.↵
    1. Graver MA,
    2. Wade JJ
    . 2011. The role of acidification in the inhibition of Neisseria gonorrhoeae by vaginal lactobacilli during anaerobic growth. Ann Clin Microbiol Antimicrob 10:8. doi:10.1186/1476-0711-10-8.
    OpenUrlCrossRefPubMed
  23. 23.↵
    1. Zheng HY,
    2. Alcorn TM,
    3. Cohen MS
    . 1994. Effects of H2O2-producing lactobacilli on Neisseria gonorrhoeae growth and catalase activity. J Infect Dis 170:1209–1215. doi:10.1093/infdis/170.5.1209.
    OpenUrlCrossRefPubMedWeb of Science
  24. 24.↵
    1. Montagnini Spaine D,
    2. Mamizuka EM,
    3. Pereira Cedenho A,
    4. Srougi M
    . 2000. Microbiologic aerobic studies on normal male urethra. Urology 56:207–210. doi:10.1016/S0090-4295(00)00615-4.
    OpenUrlCrossRefPubMed
  25. 25.↵
    1. Ketterer MR,
    2. Rice PA,
    3. Gulati S,
    4. Kiel S,
    5. Byerly L,
    6. Fortenberry JD,
    7. Soper DE,
    8. Apicella MA
    . 2016. Desialylation of Neisseria gonorrhoeae lipooligosaccharide by cervicovaginal microbiome sialidases: the potential for enhancing infectivity in men. J Infect Dis 214:1621–1628. doi:10.1093/infdis/jiw329.
    OpenUrlCrossRefPubMed
  26. 26.↵
    1. Edwards JL,
    2. Brown EJ,
    3. Ault KA,
    4. Apicella MA
    . 2001. The role of complement receptor 3 (CR3) in Neisseria gonorrhoeae infection of human cervical epithelia. Cell Microbiol 3:611–622. doi:10.1046/j.1462-5822.2001.00140.x.
    OpenUrlCrossRefPubMed
  27. 27.↵
    1. Edwards JL,
    2. Brown EJ,
    3. Uk-Nham S,
    4. Cannon JG,
    5. Blake MS,
    6. Apicella MA
    . 2002. A co-operative interaction between Neisseria gonorrhoeae and complement receptor 3 mediates infection of primary cervical epithelial cells. Cell Microbiol 4:571–584. doi:10.1046/j.1462-5822.2002.t01-1-00215.x.
    OpenUrlCrossRefPubMed
  28. 28.↵
    1. Edwards JL,
    2. Butler EK
    . 2011. The pathobiology of Neisseria gonorrhoeae lower female genital tract infection. Front Microbiol 2:102. doi:10.3389/fmicb.2011.00102.
    OpenUrlCrossRefPubMed
  29. 29.↵
    1. Caron E,
    2. Hall A
    . 1998. Identification of two distinct mechanisms of phagocytosis controlled by different Rho GTPases. Science 282:1717–1721. doi:10.1126/science.282.5394.1717.
    OpenUrlAbstract/FREE Full Text
  30. 30.↵
    1. Quillin SJ,
    2. Seifert HS
    . 2018. Neisseria gonorrhoeae host adaptation and pathogenesis. Nat Rev Microbiol 16:226–240. doi:10.1038/nrmicro.2017.169.
    OpenUrlCrossRef
  31. 31.↵
    1. Jerse AE,
    2. Deal CD
    . 2013. Vaccine research for gonococcal infections: where are we? Sex Transm Infect 89(Suppl 4):iv63–iv68. doi:10.1136/sextrans-2013-051225.
    OpenUrlAbstract/FREE Full Text
  32. 32.↵
    1. Hobbs MM,
    2. Sparling PF,
    3. Cohen MS,
    4. Shafer WM,
    5. Deal CD,
    6. Jerse AE
    . 2011. Experimental gonococcal infection in male volunteers: cumulative experience with Neisseria gonorrhoeae strains FA1090 and MS11mkC. Front Microbiol 2:123. doi:10.3389/fmicb.2011.00123.
    OpenUrlCrossRefPubMed
  33. 33.↵
    1. Handsfield HH,
    2. Lipman TO,
    3. Harnisch JP,
    4. Tronca E,
    5. Holmes KK
    . 1974. Asymptomatic gonorrhea in men. Diagnosis, natural course, prevalence and significance. N Engl J Med 290:117–123. doi:10.1056/NEJM197401172900301.
    OpenUrlCrossRefPubMedWeb of Science
  34. 34.↵
    1. Agarwal S,
    2. King CA,
    3. Klein EK,
    4. Soper DE,
    5. Rice PA,
    6. Wetzler LM,
    7. Genco CA
    . 2005. The gonococcal Fur-regulated tbpA and tbpB genes are expressed during natural mucosal gonococcal infection. Infect Immun 73:4281–4287. doi:10.1128/IAI.73.7.4281-4287.2005.
    OpenUrlAbstract/FREE Full Text
  35. 35.↵
    1. Agarwal S,
    2. Sebastian S,
    3. Szmigielski B,
    4. Rice PA,
    5. Genco CA
    . 2008. Expression of the gonococcal global regulatory protein Fur and genes encompassing the Fur and iron regulon during in vitro and in vivo infection in women. J Bacteriol 190:3129–3139. doi:10.1128/JB.01830-07.
    OpenUrlAbstract/FREE Full Text
  36. 36.↵
    1. Su X,
    2. Jiang F,
    3. Ge Q,
    4. Dai X,
    5. Sun H,
    6. Ye S
    . 2007. Surveillance of antimicrobial susceptibilitites in Neisseria gonorrhoeae in Nanjing, China, 1999–2006. Sex Transm Dis 34:995–999.
    OpenUrlPubMed
  37. 37.↵
    1. Wan C,
    2. Li Y,
    3. Le W,
    4. Liu Y,
    5. Li S,
    6. Wang BX,
    7. Rice PA,
    8. Su X-H
    . 2018. Increasing resistance to azithromycin of Neisseria gonorrhoeae in eastern Chinese cities: mechanisms and genetic diversity of resistant Nanjing isolates. Antimicrob Agents Chemother 62:e02499-17. doi:10.1128/AAC.02499-17.
    OpenUrlAbstract/FREE Full Text
  38. 38.↵
    1. McClure R,
    2. Nudel K,
    3. Massari P,
    4. Tjaden B,
    5. Su X,
    6. Rice PA,
    7. Genco CA
    . 2015. The gonococcal transcriptome during infection of the lower genital tract in women. PLoS One 10:e0133982. doi:10.1371/journal.pone.0133982.
    OpenUrlCrossRefPubMed
  39. 39.↵
    1. McClure R,
    2. Balasubramanian D,
    3. Sun Y,
    4. Bobrovskyy M,
    5. Sumby P,
    6. Genco CA,
    7. Vanderpool CK,
    8. Tjaden B
    . 2013. Computational analysis of bacterial RNA-Seq data. Nucleic Acids Res 41:e140. doi:10.1093/nar/gkt444.
    OpenUrlCrossRefPubMed
  40. 40.↵
    1. Wackwitz B,
    2. Bongaerts J,
    3. Goodman SD,
    4. Unden G
    . 1999. Growth phase-dependent regulation of nuoA-N expression in Escherichia coli K-12 by the Fis protein: upstream binding sites and bioenergetic significance. Mol Gen Genet 262:876–883. doi:10.1007/s004380051153.
    OpenUrlCrossRefPubMedWeb of Science
  41. 41.↵
    1. Seib KL,
    2. Wu HJ,
    3. Kidd SP,
    4. Apicella MA,
    5. Jennings MP,
    6. McEwan AG
    . 2006. Defenses against oxidative stress in Neisseria gonorrhoeae: a system tailored for a challenging environment. Microbiol Mol Biol Rev 70:344–361. doi:10.1128/MMBR.00044-05.
    OpenUrlAbstract/FREE Full Text
  42. 42.↵
    1. Rouhier N,
    2. Jacquot JP
    . 2003. Molecular and catalytic properties of a peroxiredoxin-glutaredoxin hybrid from Neisseria meningitidis. FEBS Lett 554:149–153. doi:10.1016/S0014-5793(03)01156-6.
    OpenUrlCrossRefPubMed
  43. 43.↵
    1. Antoine M,
    2. Gand A,
    3. Boschi-Muller S,
    4. Branlant G
    . 2006. Characterization of the amino acids from Neisseria meningitidis MsrA involved in the chemical catalysis of the methionine sulfoxide reduction step. J Biol Chem 281:39062–39070. doi:10.1074/jbc.M608844200.
    OpenUrlAbstract/FREE Full Text
  44. 44.↵
    1. Isabella VM,
    2. Clark VL
    . 2011. Deep sequencing-based analysis of the anaerobic stimulon in Neisseria gonorrhoeae. BMC Genomics 12:51. doi:10.1186/1471-2164-12-51.
    OpenUrlCrossRefPubMed
  45. 45.↵
    1. Yu C,
    2. Genco CA
    . 2012. Fur-mediated activation of gene transcription in the human pathogen Neisseria gonorrhoeae. J Bacteriol 194:1730–1742. doi:10.1128/JB.06176-11.
    OpenUrlAbstract/FREE Full Text
  46. 46.↵
    1. Stohl EA,
    2. Criss AK,
    3. Seifert HS
    . 2005. The transcriptome response of Neisseria gonorrhoeae to hydrogen peroxide reveals genes with previously uncharacterized roles in oxidative damage protection. Mol Microbiol 58:520–532. doi:10.1111/j.1365-2958.2005.04839.x.
    OpenUrlCrossRefPubMedWeb of Science
  47. 47.↵
    1. Seib KL,
    2. Wu HJ,
    3. Srikhanta YN,
    4. Edwards JL,
    5. Falsetta ML,
    6. Hamilton AJ,
    7. Maguire TL,
    8. Grimmond SM,
    9. Apicella MA,
    10. McEwan AG,
    11. Jennings MP
    . 2007. Characterization of the OxyR regulon of Neisseria gonorrhoeae. Mol Microbiol 63:54–68. doi:10.1111/j.1365-2958.2006.05478.x.
    OpenUrlCrossRefPubMedWeb of Science
  48. 48.↵
    1. Jackson LA,
    2. Day M,
    3. Allen J,
    4. Scott E, II,
    5. Dyer DW
    . 2017. Iron-regulated small RNA expression as Neisseria gonorrhoeae FA 1090 transitions into stationary phase growth. BMC Genomics 18:317. doi:10.1186/s12864-017-3684-8.
    OpenUrlCrossRef
  49. 49.↵
    1. McClure R,
    2. Tjaden B,
    3. Genco C
    . 2014. Identification of sRNAs expressed by the human pathogen Neisseria gonorrhoeae under disparate growth conditions. Front Microbiol 5:456. doi:10.3389/fmicb.2014.00456.
    OpenUrlCrossRef
  50. 50.↵
    1. Piekarowicz A,
    2. Kłyz A,
    3. Majchrzak M,
    4. Adamczyk-Popławska M,
    5. Maugel TK,
    6. Stein DC
    . 2007. Characterization of the dsDNA prophage sequences in the genome of Neisseria gonorrhoeae and visualization of productive bacteriophage. BMC Microbiol 7:66. doi:10.1186/1471-2180-7-66.
    OpenUrlCrossRefPubMed
  51. 51.↵
    1. Hollander A,
    2. Mercante AD,
    3. Shafer WM,
    4. Cornelissen CN
    . 2011. The iron-repressed, AraC-like regulator MpeR activates expression of fetA in Neisseria gonorrhoeae. Infect Immun 79:4764–4776. doi:10.1128/IAI.05806-11.
    OpenUrlAbstract/FREE Full Text
  52. 52.↵
    1. Ducey TF,
    2. Jackson L,
    3. Orvis J,
    4. Dyer DW
    . 2009. Transcript analysis of nrrF, a Fur repressed sRNA of Neisseria gonorrhoeae. Microb Pathog 46:166–170. doi:10.1016/j.micpath.2008.12.003.
    OpenUrlCrossRefPubMed
  53. 53.↵
    1. Li S,
    2. Su XH,
    3. Le WJ,
    4. Jiang FX,
    5. Wang BX,
    6. Rice PA
    . 2014. Antimicrobial susceptibility of Neisseria gonorrhoeae isolates from symptomatic men attending the Nanjing sexually transmitted diseases clinic (2011–2012): genetic characteristics of isolates with reduced sensitivity to ceftriaxone. BMC Infect Dis 14:622. doi:10.1186/s12879-014-0622-0.
    OpenUrlCrossRefPubMed
  54. 54.↵
    1. Ilina EN,
    2. Vereshchagin VA,
    3. Borovskaya AD,
    4. Malakhova MV,
    5. Sidorenko SV,
    6. Al-Khafaji NC,
    7. Kubanova AA,
    8. Govorun VM
    . 2008. Relation between genetic markers of drug resistance and susceptibility profile of clinical Neisseria gonorrhoeae strains. Antimicrob Agents Chemother 52:2175–2182. doi:10.1128/AAC.01420-07.
    OpenUrlAbstract/FREE Full Text
  55. 55.↵
    1. Matejková P,
    2. Flasarová M,
    3. Zákoucká H,
    4. Borek M,
    5. Kremenová S,
    6. Arenberger P,
    7. Woznicová V,
    8. Weinstock GM,
    9. Smajs D
    . 2009. Macrolide treatment failure in a case of secondary syphilis: a novel A2059G mutation in the 23S rRNA gene of Treponema pallidum subsp. pallidum. J Med Microbiol 58:832–836. doi:10.1099/jmm.0.007542-0.
    OpenUrlCrossRefPubMedWeb of Science
  56. 56.↵
    1. Hagman KE,
    2. Shafer WM
    . 1995. Transcriptional control of the mtr efflux system of Neisseria gonorrhoeae. J Bacteriol 177:4162–4165. doi:10.1128/jb.177.14.4162-4165.1995.
    OpenUrlAbstract/FREE Full Text
  57. 57.↵
    1. Jerse AE,
    2. Sharma ND,
    3. Simms AN,
    4. Crow ET,
    5. Snyder LA,
    6. Shafer WM
    . 2003. A gonococcal efflux pump system enhances bacterial survival in a female mouse model of genital tract infection. Infect Immun 71:5576–5582. doi:10.1128/IAI.71.10.5576-5582.2003.
    OpenUrlAbstract/FREE Full Text
  58. 58.↵
    1. Shafer WM,
    2. Qu X,
    3. Waring AJ,
    4. Lehrer RI
    . 1998. Modulation of Neisseria gonorrhoeae susceptibility to vertebrate antibacterial peptides due to a member of the resistance/nodulation/division efflux pump family. Proc Natl Acad Sci U S A 95:1829–1833. doi:10.1073/pnas.95.4.1829.
    OpenUrlAbstract/FREE Full Text
  59. 59.↵
    1. Grad YH,
    2. Harris SR,
    3. Kirkcaldy RD,
    4. Green AG,
    5. Marks DS,
    6. Bentley SD,
    7. Trees D,
    8. Lipsitch M
    . 2016. Genomic epidemiology of gonococcal resistance to extended-spectrum cephalosporins, macrolides, and fluoroquinolones in the United States, 2000–2013. J Infect Dis 214:1579–1587. doi:10.1093/infdis/jiw420.
    OpenUrlCrossRefPubMed
  60. 60.↵
    1. Seike K,
    2. Yasuda M,
    3. Hatazaki K,
    4. Mizutani K,
    5. Yuhara K,
    6. Ito Y,
    7. Fujimoto Y,
    8. Ito S,
    9. Tsuchiya T,
    10. Yokoi S,
    11. Nakano M,
    12. Deguchi T
    . 2016. Novel penA mutations identified in Neisseria gonorrhoeae with decreased susceptibility to ceftriaxone isolated between 2000 and 2014 in Japan. J Antimicrob Chemother 71:2466–2470. doi:10.1093/jac/dkw161.
    OpenUrlCrossRefPubMed
  61. 61.↵
    1. Kandler JL,
    2. Holley CL,
    3. Reimche JL,
    4. Dhulipala V,
    5. Balthazar JT,
    6. Muszyński A,
    7. Carlson RW,
    8. Shafer WM
    . 2016. The MisR response regulator is necessary for intrinsic cationic antimicrobial peptide and aminoglycoside resistance in Neisseria gonorrhoeae. Antimicrob Agents Chemother 60:4690–4700. doi:10.1128/AAC.00823-16.
    OpenUrlAbstract/FREE Full Text
  62. 62.↵
    1. Rouquette C,
    2. Harmon JB,
    3. Shafer WM
    . 1999. Induction of the mtrCDE-encoded efflux pump system of Neisseria gonorrhoeae requires MtrA, an AraC-like protein. Mol Microbiol 33:651–658. doi:10.1046/j.1365-2958.1999.01517.x.
    OpenUrlCrossRefPubMedWeb of Science
  63. 63.↵
    1. Veal WL,
    2. Nicholas RA,
    3. Shafer WM
    . 2002. Overexpression of the MtrC-MtrD-MtrE efflux pump due to an mtrR mutation is required for chromosomally mediated penicillin resistance in Neisseria gonorrhoeae. J Bacteriol 184:5619–5624. doi:10.1128/JB.184.20.5619-5624.2002.
    OpenUrlAbstract/FREE Full Text
  64. 64.↵
    1. Lee EH,
    2. Rouquette-Loughlin C,
    3. Folster JP,
    4. Shafer WM
    . 2003. FarR regulates the farAB-encoded efflux pump of Neisseria gonorrhoeae via an MtrR regulatory mechanism. J Bacteriol 185:7145–7152. doi:10.1128/JB.185.24.7145-7152.2003.
    OpenUrlAbstract/FREE Full Text
  65. 65.↵
    1. Edwards JL,
    2. Apicella MA
    . 2004. The molecular mechanisms used by Neisseria gonorrhoeae to initiate infection differ between men and women. Clin Microbiol Rev 17:965–981. doi:10.1128/CMR.17.4.965-981.2004.
    OpenUrlAbstract/FREE Full Text
  66. 66.↵
    1. Cornelissen CN,
    2. Kelley M,
    3. Hobbs MM,
    4. Anderson JE,
    5. Cannon JG,
    6. Cohen MS,
    7. Sparling PF
    . 1998. The transferrin receptor expressed by gonococcal strain FA1090 is required for the experimental infection of human male volunteers. Mol Microbiol 27:611–616. doi:10.1046/j.1365-2958.1998.00710.x.
    OpenUrlCrossRefPubMedWeb of Science
  67. 67.↵
    1. Hobbs MM,
    2. Anderson JE,
    3. Balthazar JT,
    4. Kandler JL,
    5. Carlson RW,
    6. Ganguly J,
    7. Begum AA,
    8. Duncan JA,
    9. Lin JT,
    10. Sparling PF,
    11. Jerse AE,
    12. Shafer WM
    . 2013. Lipid A’s structure mediates Neisseria gonorrhoeae fitness during experimental infection of mice and men. mBio 4:e00892-13. doi:10.1128/mBio.00892-13.
    OpenUrlAbstract/FREE Full Text
  68. 68.↵
    1. Larsen B,
    2. Monif GR
    . 2001. Understanding the bacterial flora of the female genital tract. Clin Infect Dis 32:e69–e77. doi:10.1086/318710.
    OpenUrlCrossRefPubMedWeb of Science
  69. 69.↵
    1. Ma B,
    2. Forney LJ,
    3. Ravel J
    . 2012. Vaginal microbiome: rethinking health and disease. Annu Rev Microbiol 66:371–389. doi:10.1146/annurev-micro-092611-150157.
    OpenUrlCrossRefPubMedWeb of Science
  70. 70.↵
    1. Ravel J,
    2. Gajer P,
    3. Abdo Z,
    4. Schneider GM,
    5. Koenig SS,
    6. McCulle SL,
    7. Karlebach S,
    8. Gorle R,
    9. Russell J,
    10. Tacket CO,
    11. Brotman RM,
    12. Davis CC,
    13. Ault K,
    14. Peralta L,
    15. Forney LJ
    . 2011. Vaginal microbiome of reproductive-age women. Proc Natl Acad Sci U S A 108:4680–4687. doi:10.1073/pnas.1002611107.
    OpenUrlAbstract/FREE Full Text
  71. 71.↵
    1. Unemo M,
    2. del Rio C,
    3. Shafer WM
    . 10 June 2016. Antimicrobial resistance expressed by Neisseria gonorrhoeae: a major global public health problem in the 21st century. Microbiol Spectr doi:10.1128/microbiolspec.EI10-0009-2015.
    OpenUrlCrossRef
  72. 72.↵
    1. Gu WM,
    2. Chen Y,
    3. Yang Y,
    4. Wu L,
    5. Hu WZ,
    6. Jin YL
    . 2014. Twenty-five-year changing pattern of gonococcal antimicrobial susceptibility in Shanghai: surveillance and its impact on treatment guidelines. BMC Infect Dis 14:731. doi:10.1186/s12879-014-0731-9.
    OpenUrlCrossRefPubMed
  73. 73.↵
    1. Wi T,
    2. Lahra MM,
    3. Ndowa F,
    4. Bala M,
    5. Dillon J-AR,
    6. Ramon-Pardo P,
    7. Eremin SR,
    8. Bolan G,
    9. Unemo M
    . 2017. Antimicrobial resistance in Neisseria gonorrhoeae: global surveillance and a call for international collaborative action. PLoS Med 14:e1002344. doi:10.1371/journal.pmed.1002344.
    OpenUrlCrossRef
  74. 74.↵
    1. MacPherson DW,
    2. Gushulak BD,
    3. Baine WB,
    4. Bala S,
    5. Gubbins PO,
    6. Holtom P,
    7. Segarra-Newnham M
    . 2009. Population mobility, globalization, and antimicrobial drug resistance. Emerg Infect Dis 15:1727–1732. doi:10.3201/eid1511.090419.
    OpenUrlCrossRefPubMed
  75. 75.↵
    1. Falsetta ML,
    2. Steichen CT,
    3. McEwan AG,
    4. Cho C,
    5. Ketterer M,
    6. Shao J,
    7. Hunt J,
    8. Jennings MP,
    9. Apicella MA
    . 2011. The composition and metabolic phenotype of Neisseria gonorrhoeae biofilms. Front Microbiol 2:75. doi:10.3389/fmicb.2011.00075.
    OpenUrlCrossRefPubMed
  76. 76.↵
    1. Goytia M,
    2. Dhulipala VL,
    3. Shafer WM
    . 2013. Spermine impairs biofilm formation by Neisseria gonorrhoeae. FEMS Microbiol Lett 343:64–69. doi:10.1111/1574-6968.12130.
    OpenUrlCrossRefPubMed
  77. 77.↵
    1. Piekarowicz A,
    2. Majchrzak M,
    3. Kłyz A,
    4. Adamczyk-Popławska M
    . 2006. Analysis of the filamentous bacteriophage genomes integrated into Neisseria gonorrhoeae FA1090 chromosome. Pol J Microbiol 55:251–260.
    OpenUrlCrossRefPubMed
  78. 78.↵
    1. Daou N,
    2. Yu C,
    3. McClure R,
    4. Gudino C,
    5. Reed GW,
    6. Genco CA
    . 2013. Neisseria prophage repressor implicated in gonococcal pathogenesis. Infect Immun 81:3652–3661. doi:10.1128/IAI.00298-13.
    OpenUrlAbstract/FREE Full Text
  79. 79.↵
    1. Klein S,
    2. Lorenzo C,
    3. Hoffmann S,
    4. Walther JM,
    5. Storbeck S,
    6. Piekarski T,
    7. Tindall BJ,
    8. Wray V,
    9. Nimtz M,
    10. Moser J
    . 2009. Adaptation of Pseudomonas aeruginosa to various conditions includes tRNA-dependent formation of alanyl-phosphatidylglycerol. Mol Microbiol 71:551–565. doi:10.1111/j.1365-2958.2008.06562.x.
    OpenUrlCrossRefPubMed
  80. 80.↵
    1. Shepherd J,
    2. Ibba M
    . 2013. Direction of aminoacylated transfer RNAs into antibiotic synthesis and peptidoglycan-mediated antibiotic resistance. FEBS Lett 587:2895–2904. doi:10.1016/j.febslet.2013.07.036.
    OpenUrlCrossRefPubMedWeb of Science
  81. 81.↵
    1. Lin JS,
    2. Donegan SP,
    3. Heeren TC,
    4. Greenberg M,
    5. Flaherty EE,
    6. Haivanis R,
    7. Su XH,
    8. Dean D,
    9. Newhall WJ,
    10. Knapp JS,
    11. Sarafian SK,
    12. Rice RJ,
    13. Morse SA,
    14. Rice PA
    . 1998. Transmission of Chlamydia trachomatis and Neisseria gonorrhoeae among men with urethritis and their female sex partners. J Infect Dis 178:1707–1712. doi:10.1086/314485.
    OpenUrlCrossRefPubMedWeb of Science
  82. 82.↵
    1. Sherrard J,
    2. Barlow D
    . 1996. Gonorrhoea in men: clinical and diagnostic aspects. Genitourin Med 72:422–426. doi:10.1136/sti.72.6.422.
    OpenUrlAbstract/FREE Full Text
  83. 83.↵
    1. Jensen JS,
    2. Uldum SA,
    3. Søndergård-Andersen J,
    4. Vuust J,
    5. Lind K
    . 1991. Polymerase chain reaction for detection of Mycoplasma genitalium in clinical samples. J Clin Microbiol 29:46–50.
    OpenUrlAbstract/FREE Full Text
  84. 84.↵
    1. Cheng G,
    2. He Y,
    3. Zhou X,
    4. Li H
    . 2002. Development of Chlamydia trachomatis fluorescence PCR diagnostic kit and its clinical trial. Acad J Sun Yat-Sen 23:397–400.
    OpenUrl
  85. 85.↵
    1. Feavers IM,
    2. Maiden MC
    . 1998. A gonococcal porA pseudogene: implications for understanding the evolution and pathogenicity of Neisseria gonorrhoeae. Mol Microbiol 30:647–656. doi:10.1046/j.1365-2958.1998.01101.x.
    OpenUrlCrossRefPubMedWeb of Science
  86. 86.↵
    1. CLSI
    . 2017. Performance standards for antimicrobial susceptibility testing, 27th ed. Clinical and Laboratory Standards Institute, Wayne, PA.
  87. 87.↵
    1. Liao CH,
    2. Lai CC,
    3. Hsu MS,
    4. Chu FY,
    5. Wu MY,
    6. Huang YT,
    7. Hsueh PR
    . 2010. Antimicrobial susceptibility of Neisseria gonorrhoeae isolates determined by the agar dilution, disk diffusion and Etest methods: comparison of results using GC agar and chocolate agar. Int J Antimicrob Agents 35:457–460. doi:10.1016/j.ijantimicag.2010.01.007.
    OpenUrlCrossRefPubMed
  88. 88.↵
    1. Bibby J,
    2. Kent J,
    3. Mardia K
    . 1979. Multivariate analysis. Academic Press, London, United Kingdom.
  89. 89.↵
    1. Chin CS,
    2. Alexander DH,
    3. Marks P,
    4. Klammer AA,
    5. Drake J,
    6. Heiner C,
    7. Clum A,
    8. Copeland A,
    9. Huddleston J,
    10. Eichler EE,
    11. Turner SW,
    12. Korlach J
    . 2013. Nonhybrid, finished microbial genome assemblies from long-read SMRT sequencing data. Nat Methods 10:563–569. doi:10.1038/nmeth.2474.
    OpenUrlCrossRefPubMedWeb of Science
  90. 90.↵
    1. Gardner SN,
    2. Slezak T,
    3. Hall BG
    . 2015. kSNP3.0: SNP detection and phylogenetic analysis of genomes without genome alignment or reference genome. Bioinformatics 31:2877–2878. doi:10.1093/bioinformatics/btv271.
    OpenUrlCrossRefPubMed
  91. 91.↵
    1. Jia B,
    2. Raphenya AR,
    3. Alcock B,
    4. Waglechner N,
    5. Guo P,
    6. Tsang KK,
    7. Lago BA,
    8. Dave BM,
    9. Pereira S,
    10. Sharma AN,
    11. Doshi S,
    12. Courtot M,
    13. Lo R,
    14. Williams LE,
    15. Frye JG,
    16. Elsayegh T,
    17. Sardar D,
    18. Westman EL,
    19. Pawlowski AC,
    20. Johnson TA,
    21. Brinkman FS,
    22. Wright GD,
    23. McArthur AG
    . 2017. CARD 2017: expansion and model-centric curation of the comprehensive antibiotic resistance database. Nucleic Acids Res 45:D566–D573. doi:10.1093/nar/gkw1004.
    OpenUrlCrossRefPubMed
  92. 92.↵
    1. Harrison OB,
    2. Clemence M,
    3. Dillard JP,
    4. Tang CM,
    5. Trees D,
    6. Grad YH,
    7. Maiden MC
    . 2016. Genomic analyses of Neisseria gonorrhoeae reveal an association of the gonococcal genetic island with antimicrobial resistance. J Infect 73:578–587. doi:10.1016/j.jinf.2016.08.010.
    OpenUrlCrossRef
PreviousNext
Back to top
Download PDF
Citation Tools
Transcriptome Analysis of Neisseria gonorrhoeae during Natural Infection Reveals Differential Expression of Antibiotic Resistance Determinants between Men and Women
Kathleen Nudel, Ryan McClure, Matthew Moreau, Emma Briars, A. Jeanine Abrams, Brian Tjaden, Xiao-Hong Su, David Trees, Peter A. Rice, Paola Massari, Caroline A. Genco
mSphere Jun 2018, 3 (3) e00312-18; DOI: 10.1128/mSphereDirect.00312-18

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Print
Alerts
Sign In to Email Alerts with your Email Address
Email

Thank you for sharing this mSphere article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Transcriptome Analysis of Neisseria gonorrhoeae during Natural Infection Reveals Differential Expression of Antibiotic Resistance Determinants between Men and Women
(Your Name) has forwarded a page to you from mSphere
(Your Name) thought you would be interested in this article in mSphere.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
Transcriptome Analysis of Neisseria gonorrhoeae during Natural Infection Reveals Differential Expression of Antibiotic Resistance Determinants between Men and Women
Kathleen Nudel, Ryan McClure, Matthew Moreau, Emma Briars, A. Jeanine Abrams, Brian Tjaden, Xiao-Hong Su, David Trees, Peter A. Rice, Paola Massari, Caroline A. Genco
mSphere Jun 2018, 3 (3) e00312-18; DOI: 10.1128/mSphereDirect.00312-18
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Top
  • Article
    • ABSTRACT
    • INTRODUCTION
    • RESULTS
    • DISCUSSION
    • MATERIALS AND METHODS
    • ACKNOWLEDGMENTS
    • FOOTNOTES
    • REFERENCES
  • Figures & Data
  • Info & Metrics
  • PDF

KEYWORDS

Neisseria gonorrhoeae
RNA-seq
antibiotic resistance
human mucosal infection

Related Articles

Cited By...

About

  • About mSphere
  • Board of Editors
  • Policies
  • For Reviewers
  • For the Media
  • Embargo Policy
  • For Librarians
  • For Advertisers
  • Alerts
  • RSS
  • FAQ
  • Permissions
  • Journal Announcements

Authors

  • ASM Author Center
  • Submit a Manuscript
  • Author Warranty
  • Types of Articles
  • Getting Started
  • Ethics
  • Contact Us

Follow #mSphereJ

@ASMicrobiology

       

 

Website feedback

ASM Journals

ASM journals are the most prominent publications in the field, delivering up-to-date and authoritative coverage of both basic and clinical microbiology.

About ASM | Contact Us | Press Room

 

ASM is a member of

Scientific Society Publisher Alliance

 

American Society for Microbiology
1752 N St. NW
Washington, DC 20036
Phone: (202) 737-3600

Copyright © 2021 American Society for Microbiology | Privacy Policy | Website feedback

Online ISSN: 2379-5042