A Genomic Perspective on the Near-term Impact of Doxycycline Post-exposure Prophylaxis on Neisseria gonorrhoeae Antimicrobial Resistance

Pre-existing tetracycline resistance in Neisseria gonorrhoeae limits the effectiveness of post-exposure prophylaxis (PEP) with doxycycline against gonorrhea, and selection for tetracycline resistance may influence prevalence of multi-drug resistant strains. Using genomic and antimicrobial susceptibility data from N. gonorrhoeae, we assessed the near-term impact of doxycycline PEP on N. gonorrhoeae resistance.

Whole-Genome Sequencing to Predict Antimicrobial Susceptibility Profiles in Neisseria gonorrhoeae

BACKGROUND

Neisseria gonorrhoeae is a major public health problem due to increasing incidence and antimicrobial resistance. Genetic markers of reduced susceptibility have been identified; the extent to which those are representative of global antimicrobial resistance is unknown. We evaluated the performance of whole-genome sequencing (WGS) used to predict susceptibility to ciprofloxacin and other antimicrobials using a global collection of N. gonorrhoeae isolates.

METHODS

Susceptibility testing of common antimicrobials and the recently developed zolifodacin was performed using agar dilution to determine minimum inhibitory concentrations (MICs). We identified resistance alleles at loci known to contribute to antimicrobial resistance in N. gonorrhoeae from WGS data. We tested the ability of each locus to predict antimicrobial susceptibility.

RESULTS

A total of 481 N. gonorrhoeae isolates, collected between 2004 and 2019 and making up 457 unique genomes, were sourced from 5 countries. All isolates with demonstrated susceptibility to ciprofloxacin (MIC ≤0.06 μg/mL) had a wild-type gyrA codon 91. Multilocus approaches were needed to predict susceptibility to other antimicrobials. All isolates were susceptible to zoliflodacin, defined by an MIC ≤0.25 μg/mL.

CONCLUSIONS

Single marker prediction can be used to inform ciprofloxacin treatment of N. gonorrhoeae infection. A combination of molecular markers may be needed to determine susceptibility for other antimicrobials.

Neisseria gonorrhoeae diagnostic escape from a gyrA-based test for ciprofloxacin susceptibility and the effect on zoliflodacin resistance: a bacterial genetics and experimental evolution study

BACKGROUND

The aetiological bacterial agent of gonorrhoea, Neisseria gonorrhoeae, has become resistant to each of the first-line antibiotics used to treat it, including ciprofloxacin. One diagnostic approach to identify ciprofloxacin-susceptible isolates is to determine codon 91 in the gene encoding the A subunit of DNA gyrase, gyrA, where coding for the wild-type serine (gyrA91S) is associated with ciprofloxacin susceptibility and phenylalanine (gyrA91F) with resistance. The aim of this study was to investigate the possibility of diagnostic escape from gyrA susceptibility testing.

METHODS

We used bacterial genetics to introduce pairwise substitutions in GyrA positions 91 (S or F) and 95 (D, G, or N), which is a second site in GyrA associated with ciprofloxacin resistance, into five clinical isolates of N gonorrhoeae. All five isolates encoded GyrA S91F, an additional substitution in GyrA at position 95, substitutions in ParC that are known to cause an increased minimum inhibitory concentration (MIC) to ciprofloxacin, and GyrB 429D, which is associated with susceptibility to zoliflodacin (a spiropyrimidinetrione-class antibiotic in phase 3 trials for treatment of gonorrhoea). We evolved these isolates to assess for the existence of pathways to ciprofloxacin resistance (MIC ≥1 μg/mL) and measured MICs for ciprofloxacin and zoliflodacin. In parallel, we searched metagenomic data for 11 355 N gonorrhoeae clinical isolates with reported ciprofloxacin MICs that were publicly available from the European Nucleotide Archive for strains that would be identified as susceptible by gyrA codon 91-based assays.

FINDINGS

Three clinical isolates of N gonorrhoeae with substitutions in GyrA position 95 associated with resistance (G or N) maintained intermediate ciprofloxacin MICs (0·125-0·5 μg/mL), which has been associated with treatment failure, despite reversion of GyrA position 91 from phenylalanine to serine. From an in-silico analysis of the 11 355 genomes from N gonorrhoeae clinical isolates, we identified 30 isolates with gyrA codon 91 encoding a serine and a ciprofloxacin resistance-associated mutation at codon 95. The reported MICs for these isolates varied from 0·023 μg/mL to 0·25 μg/mL, including four with intermediate ciprofloxacin MICs (associated with substantially increased risk of treatment failure). Finally, through experimental evolution, one clinical isolate of N gonorrhoeae bearing GyrA 91S acquired ciprofloxacin resistance through mutations in the gene encoding for the B subunit of DNA gyrase (gyrB) that also conferred reduced susceptibility to zoliflodacin (ie, MIC ≥2 μg/mL).

INTERPRETATION

Diagnostic escape from gyrA codon 91 diagnostics could occur through either reversion of the gyrA allele or expansion of circulating lineages. N gonorrhoeae genomic surveillance efforts might benefit from including gyrB, given its potential for contributing to ciprofloxacin and zoliflodacin resistance, and diagnostic strategies that reduce the likelihood of escape, such as the incorporation of multiple target sites, should be investigated. Diagnostics that guide antibiotic therapy can have unintended consequences, including novel resistance determinants and antibiotic cross-resistance.

FUNDING

US National Institutes of Health National Institute of Allergy and Infectious Diseases, National Institute of General Medical Sciences, and the Smith Family Foundation.

A genomic perspective on the near-term impact of doxycycline post-exposure prophylaxis on Neisseria gonorrhoeae antimicrobial resistance

Post-exposure prophylaxis with doxycycline (doxyPEP) is being introduced to prevent bacterial sexually transmitted infections (STIs). Pre-existing tetracycline resistance in Neisseria gonorrhoeae limits doxyPEP effectiveness against gonorrhea, and selection for tetracycline resistant lineages may influence prevalence of resistance to other antimicrobials via selection for multi-drug resistant strains. Using genomic and antimicrobial susceptibility data from 5,644 clinical isolates of N. gonorrhoeae , we assessed the near-term impact of doxyPEP on N. gonorrhoeae antimicrobial resistance. We found that the impact on antimicrobial resistance is likely to be influenced by the strength of selection for plasmid-encoded and chromosomally-encoded tetracycline resistance, as isolates with high-level, plasmid-encoded resistance had lower minimum inhibitory concentrations to other antimicrobials compared to isolates with low-level tetracycline resistance. The impact of doxyPEP may differ across demographic groups and geographic regions within the United States due to variation in pre-existing tetracycline resistance.

Machine learning models for Neisseria gonorrhoeae antimicrobial susceptibility tests

Neisseria gonorrhoeae is an urgent public health threat due to the emergence of antibiotic resistance. As most isolates in the United States are susceptible to at least one antibiotic, rapid molecular antimicrobial susceptibility tests (ASTs) would offer the opportunity to tailor antibiotic therapy, thereby expanding treatment options. With genome sequence and antibiotic resistance phenotype data for nearly 20,000 clinical N. gonorrhoeae isolates now available, there is an opportunity to use statistical methods to develop sequence-based diagnostics that predict antibiotic susceptibility from genotype. N. gonorrhoeae, therefore, provides a useful example illustrating how to apply machine learning models to aid in the design of sequence-based ASTs. We present an overview of this framework, which begins with establishing the assay technology, the performance criteria, the population in which the diagnostic will be used, and the clinical goals, and extends to the choices that must be made to arrive at a set of features with the desired properties for predicting susceptibility phenotype from genotype. While we focus on the example of N. gonorrhoeae, the framework generalizes to other organisms for which large-scale genotype and antibiotic resistance data can be combined to aid in diagnostics development.

Sculpting the Bacterial O-Glycoproteome: Functional Analyses of Orthologous Oligosaccharyltransferases with Diverse Targeting Specificities

Protein glycosylation systems are widely recognized in bacteria, including members of the genus Neisseria. In most bacterial species, the molecular mechanisms and evolutionary contexts underpinning target protein selection and the glycan repertoire remain poorly understood. Broad-spectrum O-linked protein glycosylation occurs in all human-associated species groups within the genus Neisseria, but knowledge of their individual glycoprotein repertoires is limited. Interestingly, PilE, the pilin subunit of the type IV pilus (Tfp) colonization factor, is glycosylated in Neisseria gonorrhoeae and Neisseria meningitidis but not in the deeply branching species N. elongata subsp. glycolytica. To examine this in more detail, we assessed PilE glycosylation status across the genus and found that PilEs of commensal clade species are not modified by the gonococcal PglO oligosaccharyltransferase. Experiments using PglO oligosaccharyltransferases from across the genus expressed in N. gonorrhoeae showed that although all were capable of broad-spectrum protein glycosylation, those from a deep-branching group of commensals were unable to support resident PilE glycosylation. Further glycoproteomic analyses of these strains using immunoblotting and mass spectrometry revealed other proteins differentially targeted by otherwise remarkably similar oligosaccharyltransferases. Finally, we generated pglO allelic chimeras that begin to localize PglO protein domains associated with unique substrate targeting activities. These findings reveal previously unappreciated differences within the protein glycosylation systems of highly related bacterial species. We propose that the natural diversity manifest in the neisserial protein substrates and oligosaccharyltransferases has significant potential to inform the structure-function relationships operating in these and related bacterial protein glycosylation systems. IMPORTANCE Although general protein glycosylation systems have been well recognized in prokaryotes, the processes governing their distribution, function, and evolution remain poorly understood. Here, we have begun to address these gaps in knowledge by comparative analyses of broad-spectrum O-linked protein glycosylation manifest in species within the genus Neisseria that strictly colonize humans. Using N. gonorrhoeae as a well-defined model organism in conjunction with comparative genomics, intraspecies gene complementation, and glycoprotein phenotyping, we discovered clear differences in both glycosylation susceptibilities and enzymatic targeting activities of otherwise largely conserved proteins. These findings reveal previously unappreciated differences within the protein glycosylation systems of highly related bacterial species. We propose that the natural diversity manifest within Neisseria species has significant potential to elucidate the structure-function relationships operating in these and related systems and to inform novel approaches to applied glycoengineering strategies.

Rapid adaptation of a complex trait during experimental evolution of Mycobacterium tuberculosis

Tuberculosis (TB), caused by Mycobacterium tuberculosis (M. tb), is a leading cause of death due to infectious disease. TB is not traditionally associated with biofilms, but M. tb biofilms are linked with drug and immune tolerance and there is increasing recognition of their contribution to the recalcitrance of TB infections. Here, we used M. tb experimental evolution to investigate this complex phenotype and identify candidate loci controlling biofilm formation. We identified novel candidate loci, adding to our understanding of the genetic architecture underlying M. tb biofilm development. Under selective pressure to grow as a biofilm, regulatory mutations rapidly swept to fixation and were associated with changes in multiple traits, including extracellular matrix production, cell size, and growth rate. Genetic and phenotypic paths to enhanced biofilm growth varied according to the genetic background of the parent strain, suggesting that epistatic interactions are important in M. tb adaptation to changing environments.

A community-driven resource for genomic epidemiology and antimicrobial resistance prediction of Neisseria gonorrhoeae at Pathogenwatch

BACKGROUND

Antimicrobial-resistant (AMR) Neisseria gonorrhoeae is an urgent threat to public health, as strains resistant to at least one of the two last-line antibiotics used in empiric therapy of gonorrhoea, ceftriaxone and azithromycin, have spread internationally. Whole genome sequencing (WGS) data can be used to identify new AMR clones and transmission networks and inform the development of point-of-care tests for antimicrobial susceptibility, novel antimicrobials and vaccines. Community-driven tools that provide an easy access to and analysis of genomic and epidemiological data is the way forward for public health surveillance.

METHODS

Here we present a public health-focussed scheme for genomic epidemiology of N. gonorrhoeae at Pathogenwatch ( https://pathogen.watch/ngonorrhoeae ). An international advisory group of experts in epidemiology, public health, genetics and genomics of N. gonorrhoeae was convened to inform on the utility of current and future analytics in the platform. We implement backwards compatibility with MLST, NG-MAST and NG-STAR typing schemes as well as an exhaustive library of genetic AMR determinants linked to a genotypic prediction of resistance to eight antibiotics. A collection of over 12,000 N. gonorrhoeae genome sequences from public archives has been quality-checked, assembled and made public together with available metadata for contextualization.

RESULTS

AMR prediction from genome data revealed specificity values over 99% for azithromycin, ciprofloxacin and ceftriaxone and sensitivity values around 99% for benzylpenicillin and tetracycline. A case study using the Pathogenwatch collection of N. gonorrhoeae public genomes showed the global expansion of an azithromycin-resistant lineage carrying a mosaic mtr over at least the last 10 years, emphasising the power of Pathogenwatch to explore and evaluate genomic epidemiology questions of public health concern.

CONCLUSIONS

The N. gonorrhoeae scheme in Pathogenwatch provides customised bioinformatic pipelines guided by expert opinion that can be adapted to public health agencies and departments with little expertise in bioinformatics and lower-resourced settings with internet connection but limited computational infrastructure. The advisory group will assess and identify ongoing public health needs in the field of gonorrhoea, particularly regarding gonococcal AMR, in order to further enhance utility with modified or new analytic methods.

Emergence and evolution of antimicrobial resistance genes and mutations in Neisseria gonorrhoeae

BACKGROUND

Antimicrobial resistance in Neisseria gonorrhoeae is a global health concern. Strains from two internationally circulating sequence types, ST-7363 and ST-1901, have acquired resistance to third-generation cephalosporins, mainly due to mosaic penA alleles. These two STs were first detected in Japan; however, the timeline, mechanism, and process of emergence and spread of these mosaic penA alleles to other countries remain unknown.

METHODS

We studied the evolution of penA alleles by obtaining the complete genomes from three Japanese ST-1901 clinical isolates harboring mosaic penA allele 34 (penA-34) dating from 2005 and generating a phylogenetic representation of 1075 strains sampled from 35 countries. We also sequenced the genomes of 103 Japanese ST-7363 N. gonorrhoeae isolates from 1996 to 2005 and reconstructed a phylogeny including 88 previously sequenced genomes.

RESULTS

Based on an estimate of the time-of-emergence of ST-1901 (harboring mosaic penA-34) and ST-7363 (harboring mosaic penA-10), and > 300 additional genome sequences of Japanese strains representing multiple STs isolated in 1996-2015, we suggest that penA-34 in ST-1901 was generated from penA-10 via recombination with another Neisseria species, followed by recombination with a gonococcal strain harboring wildtype penA-1. Following the acquisition of penA-10 in ST-7363, a dominant sub-lineage rapidly acquired fluoroquinolone resistance mutations at GyrA 95 and ParC 87-88, by independent mutations rather than horizontal gene transfer. Data in the literature suggest that the emergence of these resistance determinants may reflect selection from the standard treatment regimens in Japan at that time.

CONCLUSIONS

Our findings highlight how antibiotic use and recombination across and within Neisseria species intersect in driving the emergence and spread of drug-resistant gonorrhea.

Using rapid point-of-care tests to inform antibiotic choice to mitigate drug resistance in gonorrhoea

Background

The first cases of extensively drug resistant gonorrhoea were recorded in the United Kingdom in 2018. There is a public health need for strategies on how to deploy existing and novel antibiotics to minimise the risk of resistance development. As rapid point-of-care tests (POCTs) to predict susceptibility are coming to clinical use, coupling the introduction of an antibiotic with diagnostics that can slow resistance emergence may offer a novel paradigm for maximising antibiotic benefits. Gepotidacin is a novel antibiotic with known resistance and resistance-predisposing mutations. In particular, a mutation that confers resistance to ciprofloxacin acts as the ‘stepping-stone’ mutation to gepotidacin resistance.

Aim

To investigate how POCTs detecting Neisseria gonorrhoeae resistance mutations for ciprofloxacin and gepotidacin can be used to minimise the risk of resistance development to gepotidacin.

Methods

We use individual-based stochastic simulations to formally investigate the aim.

Results

The level of testing needed to reduce the risk of resistance development depends on the mutation rate under treatment and the prevalence of stepping-stone mutations. A POCT is most effective if the mutation rate under antibiotic treatment is no more than two orders of magnitude above the mutation rate without treatment and the prevalence of stepping-stone mutations is 1–13%.

Conclusion

Mutation frequencies and rates should be considered when estimating the POCT usage required to reduce the risk of resistance development in a given population. Molecular POCTs for resistance mutations and stepping-stone mutations to resistance are likely to become important tools in antibiotic stewardship.

Disseminated Gonococcal Infection Complicated by Prosthetic Joint Infection: Case Report and Genomic and Phylogenetic Analysis

Neisseria gonorrhoeae infections have been increasing globally, with prevalence rising across age groups. In this study, we report a case of disseminated gonococcal infection (DGI) involving a prosthetic joint, and we use whole-genome sequencing to characterize resistance genes, putative virulence factors, and the phylogenetic lineage of the infecting isolate. We review the literature on sequence-based prediction of antibiotic resistance and factors that contribute to risk for DGI. We argue for routine sequencing and reporting of invasive gonococcal infections to aid in determining whether an invasive gonococcal infection is sporadic or part of an outbreak and to accelerate understanding of the genetic features of N gonorrhoeae that contribute to pathogenesis.

Increased power from conditional bacterial genome-wide association identifies macrolide resistance mutations in Neisseria gonorrhoeae

The emergence of resistance to azithromycin complicates treatment of Neisseria gonorrhoeae, the etiologic agent of gonorrhea. Substantial azithromycin resistance remains unexplained after accounting for known resistance mutations. Bacterial genome-wide association studies (GWAS) can identify novel resistance genes but must control for genetic confounders while maintaining power. Here, we show that compared to single-locus GWAS, conducting GWAS conditioned on known resistance mutations reduces the number of false positives and identifies a G70D mutation in the RplD 50S ribosomal protein L4 as significantly associated with increased azithromycin resistance (p-value = 1.08 × 10-11). We experimentally confirm our GWAS results and demonstrate that RplD G70D and other macrolide binding site mutations are prevalent (present in 5.42% of 4850 isolates) and widespread (identified in 21/65 countries across two decades). Overall, our findings demonstrate the utility of conditional associations for improving the performance of microbial GWAS and advance our understanding of the genetic basis of macrolide resistance.

The distribution and spread of susceptible and resistant Neisseria gonorrhoeae across demographic groups in a major metropolitan center

Background

Genomic epidemiology studies of gonorrhea in the United States have primarily focused on national surveillance for antibiotic resistance, and patterns of local transmission between demographic groups of resistant and susceptible strains are unknown.

Methods

We analyzed a convenience sample of genome sequences, antibiotic susceptibility, and patient data from 897 gonococcal isolates cultured at the New York City (NYC) Public Health Laboratory from NYC Department of Health and Mental Hygiene (DOHMH) Sexual Health Clinic (SHC) patients, primarily in 2012–2013. We reconstructed the gonococcal phylogeny, defined transmission clusters using a 10 nonrecombinant single nucleotide polymorphism threshold, tested for clustering of demographic groups, and placed NYC isolates in a global phylogenetic context.

Results

The NYC gonococcal phylogeny reflected global diversity with isolates from 22/23 of the prevalent global lineages (96%). Isolates clustered on the phylogeny by patient sexual behavior (P < .001) and race/ethnicity (P < .001). Minimum inhibitory concentrations were higher across antibiotics in isolates from men who have sex with men compared to heterosexuals (P < .001) and white heterosexuals compared to black heterosexuals (P < .01). In our dataset, all large transmission clusters (≥10 samples) of N. gonorrhoeae were susceptible to ciprofloxacin, ceftriaxone, and azithromycin, and comprised isolates from patients across demographic groups.

Conclusions

All large transmission clusters were susceptible to gonorrhea therapies, suggesting that resistance to empiric therapy was not a main driver of spread, even as risk for resistance varied across demographic groups. Further study of local transmission networks is needed to identify drivers of transmission.

Adaptation to the cervical environment is associated with increased antibiotic susceptibility in Neisseria gonorrhoeae

Neisseria gonorrhoeae is an urgent public health threat due to rapidly increasing incidence and antibiotic resistance. In contrast with the trend of increasing resistance, clinical isolates that have reverted to susceptibility regularly appear, prompting questions about which pressures compete with antibiotics to shape gonococcal evolution. Here, we used genome-wide association to identify loss-of-function (LOF) mutations in the efflux pump mtrCDE operon as a mechanism of increased antibiotic susceptibility and demonstrate that these mutations are overrepresented in cervical relative to urethral isolates. This enrichment holds true for LOF mutations in another efflux pump, farAB, and in urogenitally-adapted versus typical N. meningitidis, providing evidence for a model in which expression of these pumps in the female urogenital tract incurs a fitness cost for pathogenic Neisseria. Overall, our findings highlight the impact of integrating microbial population genomics with host metadata and demonstrate how host environmental pressures can lead to increased antibiotic susceptibility.

Targeted surveillance strategies for efficient detection of novel antibiotic resistance variants

Genotype-based diagnostics for antibiotic resistance represent a promising alternative to empiric therapy, reducing inappropriate antibiotic use. However, because such assays infer resistance based on known genetic markers, their utility will wane with the emergence of novel resistance. Maintenance of these diagnostics will therefore require surveillance to ensure early detection of novel resistance variants, but efficient strategies to do so remain undefined. We evaluate the efficiency of targeted sampling approaches informed by patient and pathogen characteristics in detecting antibiotic resistance and diagnostic escape variants in Neisseria gonorrhoeae, a pathogen associated with a high burden of disease and antibiotic resistance and the development of genotype-based diagnostics. We show that patient characteristic-informed sampling is not a reliable strategy for efficient variant detection. In contrast, sampling informed by pathogen characteristics, such as genomic diversity and genomic background, is significantly more efficient than random sampling in identifying genetic variants associated with resistance and diagnostic escape.

Lateral Gene Transfer Shapes Diversity of Gardnerella spp

Gardnerella spp. are pathognomonic for bacterial vaginosis, which increases the risk of preterm birth and the transmission of sexually transmitted infections. Gardnerella spp. are genetically diverse, comprising what have recently been defined as distinct species with differing functional capacities. Disease associations with Gardnerella spp. are not straightforward: patients with BV are usually infected with multiple species, and Gardnerella spp. are also found in the vaginal microbiome of healthy women. Genome comparisons of Gardnerella spp. show evidence of lateral gene transfer (LGT), but patterns of LGT have not been characterized in detail. Here we sought to define the role of LGT in shaping the genetic structure of Gardnerella spp. We analyzed whole genome sequencing data for 106 Gardnerella strains and used these data for pan genome analysis and to characterize LGT in the core and accessory genomes, over recent and remote timescales. In our diverse sample of Gardnerella strains, we found that both the core and accessory genomes are clearly differentiated in accordance with newly defined species designations. We identified putative competence and pilus assembly genes across most species; we also found them to be differentiated between species. Competence machinery has diverged in parallel with the core genome, with selection against deleterious mutations as a predominant influence on their evolution. By contrast, the virulence factor vaginolysin, which encodes a toxin, appears to be readily exchanged among species. We identified five distinct prophage clusters in Gardnerella genomes, two of which appear to be exchanged between Gardnerella species. Differences among species are apparent in their patterns of LGT, including their exchange with diverse gene pools. Despite frequent LGT and co-localization in the same niche, our results show that Gardnerella spp. are clearly genetically differentiated and yet capable of exchanging specific genetic material. This likely reflects complex interactions within bacterial communities associated with the vaginal microbiome. Our results provide insight into how such interactions evolve and are maintained, allowing these multi-species communities to colonize and invade human tissues and adapt to antibiotics and other stressors.

RNA polymerase mutations cause cephalosporin resistance in clinical Neisseria gonorrhoeae isolates

Increasing Neisseria gonorrhoeae resistance to ceftriaxone, the last antibiotic recommended for empiric gonorrhea treatment, poses an urgent public health threat. However, the genetic basis of reduced susceptibility to ceftriaxone is not completely understood: while most ceftriaxone resistance in clinical isolates is caused by target site mutations in penA, some isolates lack these mutations. We show that penA-independent ceftriaxone resistance has evolved multiple times through distinct mutations in rpoB and rpoD. We identify five mutations in these genes that each increase resistance to ceftriaxone, including one mutation that arose independently in two lineages, and show that clinical isolates from multiple lineages are a single nucleotide change from ceftriaxone resistance. These RNA polymerase mutations cause large-scale transcriptional changes without altering susceptibility to other antibiotics, reducing growth rate, or deranging cell morphology. These results underscore the unexpected diversity of pathways to resistance and the importance of continued surveillance for novel resistance mutations.

Epidemiological and genomic determinants of tuberculosis outbreaks in First Nations communities in Canada

BACKGROUND

In Canada, tuberculosis disproportionately affects foreign-born and First Nations populations. Within First Nations' peoples, a high proportion of cases occur in association with outbreaks. Tuberculosis transmission in the context of outbreaks is thought to result from the convergence of several factors including characteristics of the cases, contacts, the environment, and the pathogen.

METHODS

We examined the epidemiological and genomic determinants of two well-characterized tuberculosis outbreaks attributed to two super-spreaders among First Nations in the province of Alberta. These outbreaks were associated with two distinct DNA fingerprints (restriction fragment-length polymorphisms or RFLPs 0.0142 and 0.0728). We compared outbreak isolates with endemic isolates not spatio-temporarily linked to outbreak cases. We extracted epidemiological variables pertaining to tuberculosis cases and contacts from individual public health records and the provincial tuberculosis registry. We conducted group analyses using parametric and non-parametric statistical tests. We carried out whole-genome sequencing and bioinformatic analysis using validated protocols.

RESULTS

We observed differences between outbreak and endemic groups in the mean number of total and child-aged contacts and the number of contacts with new positive and converted tuberculin skin tests in all group comparisons (p < 0.05). Differences were also detected in the proportion of cases with cavitation on a chest radiograph and the mean number of close contacts in selected group comparisons (p < 0.02). A phylogenetic network analysis of whole-genome sequencing data indicated that most outbreak and endemic strains were closely related to the source case for the 0.0142 fingerprint. For the 0.0728 fingerprint, the source case haplotype was circulating among endemic cases prior to the outbreak. Genetic and temporal distances were not correlated for either RFLP 0.0142 (r² = - 0.05) or RFLP 0.0728 (r² = 0.09) when all isolates were analyzed.

CONCLUSIONS

We found no evidence that endemic strains acquired mutations resulting in their emergence in outbreak form. We conclude that the propagation of these outbreaks was likely driven by the combination of characteristics of the source cases, contacts, and the environment. The role of whole-genome sequencing in understanding mycobacterial evolution and in assisting public health authorities in conducting contact investigations and managing outbreaks is important and expected to grow in the future.

Applications of genomics to slow the spread of multidrug-resistant Neisseria gonorrhoeae

Infections with Neisseria gonorrhoeae, a sexually transmitted pathogen that causes urethritis, cervicitis, and more severe complications, are increasing. Gonorrhea is typically treated with antibiotics; however, N. gonorrhoeae has rapidly acquired resistance to many antibiotic classes, and lineages with reduced susceptibility to the currently recommended therapies are emerging worldwide. In this review, we discuss the contributions of whole genome sequencing (WGS) to our understanding of resistant N. gonorrhoeae. Genomics has illuminated the evolutionary origins and population structure of N. gonorrhoeae and the magnitude of horizontal gene transfer within and between Neisseria species. WGS can be used to predict the susceptibility of N. gonorrhoeae based on known resistance determinants, track the spread of these determinants throughout the N. gonorrhoeae population, and identify novel loci contributing to resistance. WGS has also allowed more detailed epidemiological analysis of transmission of N. gonorrhoeae between individuals and populations than previously used typing methods. Ongoing N. gonorrhoeae genomics will complement other laboratory techniques to understand the biology and evolution of the pathogen, improve diagnostics and treatment in the clinic, and inform public health policies to limit the impact of antibiotic resistance.

Signatures of Selection at Drug Resistance Loci in Mycobacterium tuberculosis

Tuberculosis (TB) is the leading cause of death by an infectious disease, and global TB control efforts are increasingly threatened by drug resistance in Mycobacterium tuberculosis. Unlike most bacteria, where lateral gene transfer is an important mechanism of resistance acquisition, resistant M. tuberculosis arises solely by de novo chromosomal mutation. Using whole-genome sequencing data from two natural populations of M. tuberculosis, we characterized the population genetics of known drug resistance loci using measures of diversity, population differentiation, and convergent evolution. We found resistant subpopulations to be less diverse than susceptible subpopulations, consistent with ongoing transmission of resistant M. tuberculosis. A subset of resistance genes ("sloppy targets") were characterized by high diversity and multiple rare variants; we posit that a large genetic target for resistance and relaxation of purifying selection contribute to high diversity at these loci. For "tight targets" of selection, the path to resistance appeared narrower, evidenced by single favored mutations that arose numerous times in the phylogeny and segregated at markedly different frequencies in resistant and susceptible subpopulations. These results suggest that diverse genetic architectures underlie drug resistance in M. tuberculosis and that combined approaches are needed to identify causal mutations. Extrapolating from patterns observed for well-characterized genes, we identified novel candidate variants involved in resistance. The approach outlined here can be extended to identify resistance variants for new drugs, to investigate the genetic architecture of resistance, and when phenotypic data are available, to find candidate genetic loci underlying other positively selected traits in clonal bacteria. IMPORTANCEMycobacterium tuberculosis, the causative agent of tuberculosis (TB), is a significant burden on global health. Antibiotic treatment imposes strong selective pressure on M. tuberculosis populations. Identifying the mutations that cause drug resistance in M. tuberculosis is important for guiding TB treatment and halting the spread of drug resistance. Whole-genome sequencing (WGS) of M. tuberculosis isolates can be used to identify novel mutations mediating drug resistance and to predict resistance patterns faster than traditional methods of drug susceptibility testing. We have used WGS from natural populations of drug-resistant M. tuberculosis to characterize effects of selection for advantageous mutations on patterns of diversity at genes involved in drug resistance. The methods developed here can be used to identify novel advantageous mutations, including new resistance loci, in M. tuberculosis and other clonal pathogens.

Adaptation in a Fibronectin Binding Autolysin of Staphylococcus saprophyticus

Human-pathogenic bacteria are found in a variety of niches, including free-living, zoonotic, and microbiome environments. Identifying bacterial adaptations that enable invasive disease is an important means of gaining insight into the molecular basis of pathogenesis and understanding pathogen emergence. Staphylococcus saprophyticus, a leading cause of urinary tract infections, can be found in the environment, food, animals, and the human microbiome. We identified a selective sweep in the gene encoding the Aas adhesin, a key virulence factor that binds host fibronectin. We hypothesize that the mutation under selection (aas_2206A>C) facilitates colonization of the urinary tract, an environment where bacteria are subject to strong shearing forces. The mutation appears to have enabled emergence and expansion of a human-pathogenic lineage of S. saprophyticus. These results demonstrate the power of evolutionary genomic approaches in discovering the genetic basis of virulence and emphasize the pleiotropy and adaptability of bacteria occupying diverse niches. IMPORTANCEStaphylococcus saprophyticus is an important cause of urinary tract infections (UTI) in women; such UTI are common, can be severe, and are associated with significant impacts to public health. In addition to being a cause of human UTI, S. saprophyticus can be found in the environment, in food, and associated with animals. After discovering that UTI strains of S. saprophyticus are for the most part closely related to each other, we sought to determine whether these strains are specially adapted to cause disease in humans. We found evidence suggesting that a mutation in the gene aas is advantageous in the context of human infection. We hypothesize that the mutation allows S. saprophyticus to survive better in the human urinary tract. These results show how bacteria found in the environment can evolve to cause disease.

Evolutionary Thrift: Mycobacteria Repurpose Plasmid Diversity during Adaptation of Type VII Secretion Systems

Mycobacteria have a distinct secretion system, termed type VII (T7SS), which is encoded by paralogous chromosomal loci (ESX) and associated with pathogenesis, conjugation, and metal homeostasis. Evolution of paralogous gene families is of interest because duplication is an important mechanism by which novel genes evolve, but there are potential conflicts between adaptive forces that stabilize duplications and those that enable evolution of new functions. Our objective was to delineate the adaptive forces underlying diversification of T7SS. Plasmid-borne ESX were described recently, and we found evidence that the initial duplication and divergence of ESX systems occurred on plasmids and was driven by selection for advantageous mutations. Plasmid conjugation has been linked to T7SS and type IV secretion systems (T4SS) in mycobacteria, and we discovered that T7SS and T4SS genes evolved in concert on the plasmids. We hypothesize that differentiation of plasmid ESX helps to prevent conjugation among cells harboring incompatible plasmids. Plasmid ESX appear to have been repurposed following migration to the chromosome, and there is evidence of positive selection driving further differentiation of chromosomal ESX. We hypothesize that ESX loci were initially stabilized on the chromosome by mediating their own transfer. These results emphasize the diverse adaptive paths underlying evolution of novelty, which in this case involved plasmid duplications, selection for advantageous mutations in the mobile and core genomes, migration of the loci between plasmids and chromosomes, and lateral transfer among chromosomes. We discuss further implications for the choice of model organism to study ESX functions in Mycobacterium tuberculosis.

A molecular portrait of maternal sepsis from Byzantine Troy

Pregnancy complications are poorly represented in the archeological record, despite their importance in contemporary and ancient societies. While excavating a Byzantine cemetery in Troy, we discovered calcified abscesses among a woman’s remains. Scanning electron microscopy of the tissue revealed ‘ghost cells’, resulting from dystrophic calcification, which preserved ancient maternal, fetal and bacterial DNA of a severe infection, likely chorioamnionitis. Gardnerella vaginalis and Staphylococcus saprophyticus dominated the abscesses. Phylogenomic analyses of ancient, historical, and contemporary data showed that G. vaginalis Troy fell within contemporary genetic diversity, whereas S. saprophyticus Troy belongs to a lineage that does not appear to be commonly associated with human disease today. We speculate that the ecology of S. saprophyticus infection may have differed in the ancient world as a result of close contacts between humans and domesticated animals. These results highlight the complex and dynamic interactions with our microbial milieu that underlie severe maternal infections.

DOI:http://dx.doi.org/10.7554/eLife.20983.001

Why and how have some bacteria evolved to cause illness in humans? One way to study bacterial evolution is to search for ancient samples of bacteria and use DNA sequencing technology to investigate how modern bacteria have changed from their ancestors. Understanding the evolution process may help researchers to understand how some bacteria become resistant to the antibiotics designed to kill them.

Complications that occur during pregnancy, including bacterial infections, have long been a major cause of death for women. Now, Devault, Mortimer et al. have been able to sequence the DNA of bacteria found in tissue collected from a woman buried 800 years ago in a cemetery in Troy. Some of the woman’s tissues had been well preserved because they had calcified (probably as the result of infection), which preserved their structure in a mineralized layer. Two mineralized “nodules” in the body appear to be the remains of abscesses. Some of the human DNA in the nodules came from a male, suggesting that the woman was pregnant with a boy and that the abscesses formed in placental tissue.

Sequencing the DNA of the bacteria in the abscess allowed Devault, Mortimer et al. to diagnose the woman’s infection, which was caused by two types of bacteria. One species, called Gardnerella vaginalis, is found in modern pregnancy-related infections. The DNA of the ancient samples was similar to that of modern bacteria. The other bacteria species was an ancient form of Staphylococcus saprophyticus, a type of bacteria that causes urinary tract infections. However, the DNA of the ancient S. saprophyticus bacteria is quite different to that of the bacteria found in modern humans. Instead, their DNA sequence appears more similar to forms of the bacteria that infect currently livestock. As humans lived closely with their livestock at the time the woman lived, her infection may be due to a type of bacteria that passed easily between humans and animals.

Overall, the results suggest that the disease-causing properties of bacteria can arise from a wide range of sources. In addition, Devault, Mortimer et al. have demonstrated that certain types of tissue found in archeological remains are a potential gold mine of information about the evolution of bacteria and other microbes found in the human body.

DOI:http://dx.doi.org/10.7554/eLife.20983.002

Diversity of Mycobacterium tuberculosis across Evolutionary Scales

Tuberculosis (TB) is a global public health emergency. Increasingly drug resistant strains of Mycobacterium tuberculosis (M.tb) continue to emerge and spread, highlighting adaptability of this pathogen. Most studies of M.tb evolution have relied on ‘between-host’ samples, in which each person with TB is represented by a single M.tb isolate. However, individuals with TB commonly harbor populations of M.tb numbering in the billions. Here, we use analyses of M.tb genomic data from within and between hosts to gain insight into influences shaping genetic diversity of this pathogen. We find that the amount of M.tb genetic diversity harbored by individuals with TB can vary dramatically, likely as a function of disease severity. Surprisingly, we did not find an appreciable impact of TB treatment on M.tb diversity. In examining genomic data from M.tb samples within and between hosts with TB, we find that genes involved in the regulation, synthesis, and transportation of immunomodulatory cell envelope lipids appear repeatedly in the extremes of various statistical measures of diversity. Many of these genes have been identified as possible targets of selection in other studies employing different methods and data sets. Taken together, these observations suggest that M.tb cell envelope lipids are targets of selection within hosts. Many of these lipids are specific to pathogenic mycobacteria and, in some cases, human-pathogenic mycobacteria. We speculate that rapid adaptation of cell envelope lipids is facilitated by functional redundancy, flexibility in their metabolism, and their roles mediating interactions with the host.

Tuberculosis (TB) is a grave threat to global public health and is the second leading cause of death due to infectious disease. The causative agent, Mycobacterium tuberculosis (M.tb), has emerged in increasingly drug resistant forms that hamper our efforts to control TB. We need a better understanding of M.tb adaptation to guide development of more effective TB treatment and control strategies. The goal of this study was to gain insight into M.tb evolution within individual patients with TB. We found that TB patients harbor a diverse population of M.tb. We further found evidence to suggest that the bacterial population evolves measurably in response to selection pressures imposed by the environment within hosts. Changes were particularly notable in M.tb genes involved in the regulation, synthesis, and transportation of lipids and glycolipids of the bacterial cell envelope. These findings have important implications for drug and vaccine development, and provide insight into TB host pathogen interactions.

Genomic signatures of distributive conjugal transfer among mycobacteria

Distributive conjugal transfer (DCT) is a newly described mechanism of lateral gene transfer (LGT) that results in a mosaic transconjugant structure, similar to the products of meiosis. We have tested popular LGT detection methods on whole-genome sequence data from experimental DCT transconjugants and used the best performing methods to compare genomic signatures of DCT with those of LGT through natural transformation, conjugative plasmids, and mobile genetic elements (MGE). We found that DCT results in transfer of larger chromosomal segments, that these segments are distributed more broadly around the chromosome, and that a greater proportion of the chromosome is affected by DCT than by other mechanisms of LGT. We used the best performing methods to characterize LGT in Mycobacterium canettii, the mycobacterial species most closely related to Mycobacterium tuberculosis. Patterns of LGT among M. canettii were highly distinctive. Gene flow appeared unidirectional, from lineages with minimal evidence of LGT to isolates with a substantial proportion (6–13%) of sites identified as recombinant. Among M. canettii isolates with evidence of LGT, recombinant fragments were larger and more evenly distributed relative to bacteria that undergo LGT through natural transformation, conjugative plasmids, and MGE. Spatial bias in M. canettii was also unusual in that patterns of recombinant fragment sharing mirrored overall phylogenetic structure. Based on the proportion of recombinant sites, the size of recombinant fragments, their spatial distribution and lack of association with MGE, as well as unidirectionality of DNA transfer, we conclude that DCT is the predominant mechanism of LGT among M. canettii.

A novel nonhuman primate model for influenza transmission

Studies of influenza transmission are necessary to predict the pandemic potential of emerging influenza viruses. Currently, both ferrets and guinea pigs are used in such studies, but these species are distantly related to humans. Nonhuman primates (NHP) share a close phylogenetic relationship with humans and may provide an enhanced means to model the virological and immunological events in influenza virus transmission. Here, for the first time, it was demonstrated that a human influenza virus isolate can productively infect and be transmitted between common marmosets (Callithrix jacchus), a New World monkey species. We inoculated four marmosets with the 2009 pandemic virus A/California/07/2009 (H1N1pdm) and housed each together with a naïve cage mate. We collected bronchoalveolar lavage and nasal wash samples from all animals at regular intervals for three weeks post-inoculation to track virus replication and sequence evolution. The unadapted 2009 H1N1pdm virus replicated to high titers in all four index animals by 1 day post-infection. Infected animals seroconverted and presented human-like symptoms including sneezing, nasal discharge, labored breathing, and lung damage. Transmission occurred in one cohabitating pair. Deep sequencing detected relatively few genetic changes in H1N1pdm viruses replicating in any infected animal. Together our data suggest that human H1N1pdm viruses require little adaptation to replicate and cause disease in marmosets, and that these viruses can be transmitted between animals. Marmosets may therefore be a viable model for studying influenza virus transmission.