Escherichia coli ST131 is the most frequently isolated fluoroquinolone-resistant (FQR) E. coli clone worldwide and a major cause of urinary tract and bloodstream infections. Although originally identified through its association with the CTX-M-15 extended-spectrum β-lactamase resistance gene, global genomic epidemiology studies have failed to resolve the geographical and temporal origin of the ST131 ancestor. Here, we developed a framework for the reanalysis of publically available genomes from different countries and used this data set to reconstruct the evolutionary steps that led to the emergence of FQR ST131. Using Bayesian estimation, we show that point mutations in chromosomal genes that confer FQR coincide with the first clinical use of fluoroquinolone in 1986 and illustrate the impact of this pivotal event on the rapid population expansion of ST131 worldwide from an apparent origin in North America. Furthermore, we identify virulence factor acquisition events that predate the development of FQR, suggesting that the gain of virulence-associated genes followed by the tandem development of antibiotic resistance primed the successful global dissemination of ST131.

Escherichia coli ST131 is the most frequently isolated fluoroquinolone-resistant (FQR) E. coli clone worldwide and a major cause of urinary tract and bloodstream infections. Although originally identified through its association with the CTX-M-15 extended-spectrum ␤-lactamase resistance gene, global genomic epidemiology studies have failed to resolve the geographical and temporal origin of the ST131 ancestor. Here, we developed a framework for the reanalysis of publically available genomes from different countries and used this data set to reconstruct the evolutionary steps that led to the emergence of FQR ST131. Using Bayesian estimation, we show that point mutations in chromosomal genes that confer FQR coincide with the first clinical use of fluoroquinolone in 1986 and illustrate the impact of this pivotal event on the rapid population expansion of ST131 worldwide from an apparent origin in North America. Furthermore, we identify virulence factor acquisition events that predate the development of FQR, suggesting that the gain of virulence-associated genes followed by the tandem development of antibiotic resistance primed the successful global dissemination of ST131.IMPORTANCE Escherichia coli sequence type 131 (ST131) is a recently emerged and globally disseminated multidrug-resistant clone frequently associated with human urinary tract and bloodstream infections. In this study, we have used two large publically available genomic data sets to define a number of critical steps in the evolution of this important pathogen. We show that resistance to fluoroquinolones, a class of broad-spectrum antibiotic used extensively in human medicine and veterinary practice, developed in ST131 soon after the introduction of these antibiotics in the United States, most likely in North America. We also mapped the acquisition of several fitness and virulence determinants by ST131 and demonstrate these events occurred prior to the development of fluoroquinolone resistance. Thus, ST131 has emerged by stealth, first acquiring genes associated with an increased capacity to cause human infection, and then gaining a resistance armory that has driven its massive population expansion across the globe.

ObjectivesThis prospective cohort study aimed to determine the natural history and incidence of oropharyngeal gonorrhoea and chlamydia among a cohort of men who have sex with men (MSM) over a 12-week period, and to examine risk factors associated with incident oropharyngeal infections.MethodsMSM either aged ≥18 years and had a diagnosis of oropharyngeal gonorrhoea by nucleic acid amplification test (NAAT) in the past 3 months or aged 18–35 years who were HIV-negative taking pre-exposure prophylaxis (PrEP) were eligible for this study. Enrolled men were followed up for 12 weeks. Oropharyngeal swabs were collected at week 0 (baseline) and week 12 (end of study). Between these time points, weekly saliva specimens and the number of tongue kissing, penile–oral and insertive rimming partners were collected by post. Oropharyngeal swabs and saliva specimens were tested by NAAT for Neisseria gonorrhoeae and Chlamydia trachomatis. Poisson regression was performed to examine the risk factors (weekly number of partners) associated with incident oropharyngeal gonorrhoea.ResultsA total of 100 MSM were recruited. The incidence of oropharyngeal gonorrhoea and chlamydia was 62 (95% CI 37 to 105) and 9 (95% CI 2 to 35)/100 person-years, respectively. The median duration of incident oropharyngeal infection with gonorrhoea was 28 days (IQR=21–36, n=7). The incidence rate ratio (IRR) for oropharyngeal gonorrhoea increased with an increased number of kissing partners (IRR=1.08; 95% CI 1.03 to 1.12) an increased number of penile-oral sex partners (IRR=1.07, 95% CI 1.01 to 1.14) but not with an increased number of insertive rimming partners (IRR=1.11, 95% CI 0.96 to 1.29) or other demographic factors. The IRR and duration of incident oropharyngeal chlamydia were not calculated due to the small number of cases (n=2).ConclusionsMSM have a high incidence of oropharyngeal gonorrhoea and the median duration of infection was less than 3 months.

Escherichia coli Sequence Type (ST)101 is an emerging, multi-drug resistant lineage associated with carbapenem resistance. We recently completed a comprehensive genomics study on mobile genetic elements (MGEs) and their role in blaNDM-1 dissemination within the ST101 lineage. DNA methyltransferases (MTases) are also frequently associated with MGEs, with DNA methylation guiding numerous biological processes including genomic defence against foreign DNA and regulation of gene expression. The availability of Pacific Biosciences Single Molecule Real Time Sequencing data for seven ST101 strains enabled us to investigate the role of DNA methylation on a genome-wide scale (methylome). We defined the methylome of two complete (MS6192 and MS6193) and five draft (MS6194, MS6201, MS6203, MS6204, MS6207) ST101 genomes. Our analysis identified 14 putative MTases and eight N6-methyladenine DNA recognition sites, with one site that has not been described previously. Furthermore, we identified a Type I MTase encoded within a Transposon 7-like Transposon and show its acquisition leads to differences in the methylome between two almost identical isolates. Genomic comparisons with 13 previously published ST101 draft genomes identified variations in MTase distribution, consistent with MGE differences between genomes, highlighting the diversity of active MTases within strains of a single E. coli lineage. It is well established that MGEs can contribute to the evolution of E. coli due to their virulence and resistance gene repertoires. This study emphasises the potential for mobile genetic elements to also enable highly similar bacterial strains to rapidly acquire genome-wide functional differences via changes to the methylome.Impact StatementEscherichia coli ST101 is an emerging human pathogen frequently associated with carbapenem resistance. E. coli ST101 strains carry numerous mobile genetic elements that encode virulence determinants, antimicrobial resistance, and DNA methyltransferases (MTases). In this study we provide the first comprehensive analysis of the genome-wide complement of DNA methylation (methylome) in seven E. coli ST101 genomes. We identified a Transposon carrying a Type I restriction modification system that may lead to functional differences between two almost identical genomes and showed how small recombination events at a single genomic region can lead to global methylome changes across the lineage. We also showed that the distribution of MTases throughout the ST101 lineage was consistent with the presence or absence of mobile genetic elements on which they are encoded. This study shows the diversity of MTases within a single bacterial lineage and shows how strain and lineage-specific methylomes may drive host adaptation.Data SummarySequence data including reads, assemblies and motif summaries have previously been submitted to the National Center for Biotechnology Information (https://www.ncbi.nlm.nih.gov) under the BioProject Accessions: PRJNA580334, PRJNA580336, PRJNA580337, PRJNA580338, PRJNA580339, PRJNA580341 and PRJNA580340 for MS6192, MS6193, MS6194, MS6201, MS6203, MS6204 and MS6207 respectively. All supporting data, code, accessions, and protocols have been provided within the article or through supplementary data files.