Clinical and laboratory-induced colistin-resistance mechanisms in Acinetobacter baumannii.

The increasing incidence and emergence of multi-drug resistant (MDR) Acinetobacter baumannii has become a major global health concern. Colistin is a historic antimicrobial that has become commonly used as a treatment for MDR A. baumannii infections. The increase in colistin usage has been mirrored by an increase in colistin resistance. We aimed to identify the mechanisms associated with colistin resistance in A. baumannii using multiple high-throughput-sequencing technologies, including transposon-directed insertion site sequencing (TraDIS), RNA sequencing (RNAseq) and whole-genome sequencing (WGS) to investigate the genotypic changes of colistin resistance in A. baumannii. Using TraDIS, we found that genes involved in drug efflux (adeIJK), and phospholipid (mlaC, mlaF and mlaD) and lipooligosaccharide synthesis (lpxC and lpsO) were required for survival in sub-inhibitory concentrations of colistin. Transcriptomic (RNAseq) analysis revealed that expression of genes encoding efflux proteins (adeI, adeC, emrB, mexB and macAB) was enhanced in in vitro generated colistin-resistant strains. WGS of these organisms identified disruptions in genes involved in lipid A (lpxC) and phospholipid synthesis (mlaA), and in the baeS/R two-component system (TCS). We additionally found that mutations in the pmrB TCS genes were the primary colistin-resistance-associated mechanisms in three Vietnamese clinical colistin-resistant A. baumannii strains. Our results outline the entire range of mechanisms employed in A. baumannii for resistance against colistin, including drug extrusion and the loss of lipid A moieties by gene disruption or modification.

Morphological, genomic and transcriptomic responses of Klebsiella pneumoniae to the last-line antibiotic colistin.

Colistin remains one of the few antibiotics effective against multi-drug resistant (MDR) hospital pathogens, such as Klebsiella pneumoniae. Yet resistance to this last-line drug is rapidly increasing. Characterized mechanisms of colR in K. pneumoniae are largely due to chromosomal mutations in two-component regulators, although a plasmid-mediated colR mechanism has recently been uncovered. However, the effects of intrinsic colistin resistance are yet to be characterized on a whole-genome level. Here, we used a genomics-based approach to understand the mechanisms of adaptive colR acquisition in K. pneumoniae. In controlled directed-evolution experiments we observed two distinct paths to colistin resistance acquisition. Whole genome sequencing identified mutations in two colistin resistance genes: in the known colR regulator phoQ which became fixed in the population and resulted in a single amino acid change, and unstable minority variants in the recently described two-component sensor crrB. Through RNAseq and microscopy, we reveal the broad range of effects that colistin exposure has on the cell. This study is the first to use genomics to identify a population of minority variants with mutations in a colR gene in K. pneumoniae.

Genomic Analysis and Comparison of Two Gonorrhea Outbreaks.

UNLABELLED: Gonorrhea is a sexually transmitted disease causing growing concern, with a substantial increase in reported incidence over the past few years in the United Kingdom and rising levels of resistance to a wide range of antibiotics. Understanding its epidemiology is therefore of major biomedical importance, not only on a population scale but also at the level of direct transmission. However, the molecular typing techniques traditionally used for gonorrhea infections do not provide sufficient resolution to investigate such fine-scale patterns. Here we sequenced the genomes of 237 isolates from two local collections of isolates from Sheffield and London, each of which was resolved into a single type using traditional methods. The two data sets were selected to have different epidemiological properties: the Sheffield data were collected over 6 years from a predominantly heterosexual population, whereas the London data were gathered within half a year and strongly associated with men who have sex with men. Based on contact tracing information between individuals in Sheffield, we found that transmission is associated with a median time to most recent common ancestor of 3.4 months, with an upper bound of 8 months, which we used as a criterion to identify likely transmission links in both data sets. In London, we found that transmission happened predominantly between individuals of similar age, sexual orientation, and location and also with the same HIV serostatus, which may reflect serosorting and associated risk behaviors. Comparison of the two data sets suggests that the London epidemic involved about ten times more cases than the Sheffield outbreak. IMPORTANCE: The recent increases in gonorrhea incidence and antibiotic resistance are cause for public health concern. Successful intervention requires a better understanding of transmission patterns, which is not uncovered by traditional molecular epidemiology techniques. Here we studied two outbreaks that took place in Sheffield and London, United Kingdom. We show that whole-genome sequencing provides the resolution to investigate direct gonorrhea transmission between infected individuals. Combining genome sequencing with rich epidemiological information about infected individuals reveals the importance of several transmission routes and risk factors, which can be used to design better control measures.

Whole-Genome Sequencing for Routine Pathogen Surveillance in Public Health: a Population Snapshot of Invasive Staphylococcus aureus in Europe.

UNLABELLED: The implementation of routine whole-genome sequencing (WGS) promises to transform our ability to monitor the emergence and spread of bacterial pathogens. Here we combined WGS data from 308 invasive Staphylococcus aureus isolates corresponding to a pan-European population snapshot, with epidemiological and resistance data. Geospatial visualization of the data is made possible by a generic software tool designed for public health purposes that is available at the project URL (http://www.microreact.org/project/EkUvg9uY?tt=rc). Our analysis demonstrates that high-risk clones can be identified on the basis of population level properties such as clonal relatedness, abundance, and spatial structuring and by inferring virulence and resistance properties on the basis of gene content. We also show that in silico predictions of antibiotic resistance profiles are at least as reliable as phenotypic testing. We argue that this work provides a comprehensive road map illustrating the three vital components for future molecular epidemiological surveillance: (i) large-scale structured surveys, (ii) WGS, and (iii) community-oriented database infrastructure and analysis tools. IMPORTANCE: The spread of antibiotic-resistant bacteria is a public health emergency of global concern, threatening medical intervention at every level of health care delivery. Several recent studies have demonstrated the promise of routine whole-genome sequencing (WGS) of bacterial pathogens for epidemiological surveillance, outbreak detection, and infection control. However, as this technology becomes more widely adopted, the key challenges of generating representative national and international data sets and the development of bioinformatic tools to manage and interpret the data become increasingly pertinent. This study provides a road map for the integration of WGS data into routine pathogen surveillance. We emphasize the importance of large-scale routine surveys to provide the population context for more targeted or localized investigation and the development of open-access bioinformatic tools to provide the means to combine and compare independently generated data with publicly available data sets.

Identifying the effect of patient sharing on between-hospital genetic differentiation of methicillin-resistant Staphylococcus aureus.

BACKGROUND: Methicillin-resistant Staphylococcus aureus (MRSA) is one of the most common healthcare-associated pathogens. To examine the role of inter-hospital patient sharing on MRSA transmission, a previous study collected 2,214 samples from 30 hospitals in Orange County, California and showed by spa typing that genetic differentiation decreased significantly with increased patient sharing. In the current study, we focused on the 986 samples with spa type t008 from the same population. METHODS: We used genome sequencing to determine the effect of patient sharing on genetic differentiation between hospitals. Genetic differentiation was measured by between-hospital genetic diversity, F ST , and the proportion of nearly identical isolates between hospitals. RESULTS: Surprisingly, we found very similar genetic diversity within and between hospitals, and no significant association between patient sharing and genetic differentiation measured by F ST . However, in contrast to F ST , there was a significant association between patient sharing and the proportion of nearly identical isolates between hospitals. We propose that the proportion of nearly identical isolates is more powerful at determining transmission dynamics than traditional estimators of genetic differentiation (F ST ) when gene flow between populations is high, since it is more responsive to recent transmission events. Our hypothesis was supported by the results from coalescent simulations. CONCLUSIONS: Our results suggested that there was a high level of gene flow between hospitals facilitated by patient sharing, and that the proportion of nearly identical isolates is more sensitive to population structure than F ST when gene flow is high.

Evolutionary Trade-Offs Underlie the Multi-faceted Virulence of Staphylococcus aureus.

Bacterial virulence is a multifaceted trait where the interactions between pathogen and host factors affect the severity and outcome of the infection. Toxin secretion is central to the biology of many bacterial pathogens and is widely accepted as playing a crucial role in disease pathology. To understand the relationship between toxicity and bacterial virulence in greater depth, we studied two sequenced collections of the major human pathogen Staphylococcus aureus and found an unexpected inverse correlation between bacterial toxicity and disease severity. By applying a functional genomics approach, we identified several novel toxicity-affecting loci responsible for the wide range in toxic phenotypes observed within these collections. To understand the apparent higher propensity of low toxicity isolates to cause bacteraemia, we performed several functional assays, and our findings suggest that within-host fitness differences between high- and low-toxicity isolates in human serum is a contributing factor. As invasive infections, such as bacteraemia, limit the opportunities for onward transmission, highly toxic strains could gain an additional between-host fitness advantage, potentially contributing to the maintenance of toxicity at the population level. Our results clearly demonstrate how evolutionary trade-offs between toxicity, relative fitness, and transmissibility are critical for understanding the multifaceted nature of bacterial virulence.

Whole genome sequence typing and microarray profiling of nasal and blood stream methicillin-resistant Staphylococcus aureus isolates: Clues to phylogeny and invasiveness.

Hospital-associated methicillin-resistant Staphylococcus aureus (MRSA) infections are frequently caused by predominant clusters of closely related isolates that cannot be discriminated by conventional diagnostic typing methods. Whole genome sequencing (WGS) and DNA microarray (MA) now allow for better discrimination within a prevalent clonal complex (CC). This single center exploratory study aims to distinguish invasive (blood stream infection) and non-invasive (nasal colonization) MRSA isolates of the same CC5 into phylogenetic- and virulence-associated genotypic subgroups by WGS and MA. A cohort of twelve blood stream and fifteen nasal MRSA isolates of CC5 (spa-types t003 and t504) was selected. Isolates were propagated at the same period of time from unrelated patients treated at the University of Saarland Medical Center, Germany. Rooted phylotyping based on WGS with core-genome single nucleotide polymorphism (SNP) analysis revealed two local clusters of closely related CC5 subgroups (t504 and Clade1 t003) which were separated from other local t003 isolates and from unrelated CC5 MRSA reference isolates of German origin. Phylogenetic subtyping was not associated with invasiveness when comparing blood stream and nasal isolates. Clustering based on MA profiles was not concordant with WGS phylotyping, but MA profiles may identify subgroups of isolates with nasal and blood stream origin. Among the new putative virulence associated genes identified by WGS, the strongest association with blood stream infections was shown for ebhB mutants. Analysis of the core-genome together with the accessory genome enables subtyping of closely related MRSA isolates according to phylogeny and presumably also to the potential virulence capacity of isolates.

Genetic Stabilization of the Drug-Resistant PMEN1 Pneumococcus Lineage by Its Distinctive DpnIII Restriction-Modification System.

UNLABELLED: The human pathogen Streptococcus pneumoniae (pneumococcus) exhibits a high degree of genomic diversity and plasticity. Isolates with high genomic similarity are grouped into lineages that undergo homologous recombination at variable rates. PMEN1 is a pandemic, multidrug-resistant lineage. Heterologous gene exchange between PMEN1 and non-PMEN1 isolates is directional, with extensive gene transfer from PMEN1 strains and only modest transfer into PMEN1 strains. Restriction-modification (R-M) systems can restrict horizontal gene transfer, yet most pneumococcal strains code for either the DpnI or DpnII R-M system and neither limits homologous recombination. Our comparative genomic analysis revealed that PMEN1 isolates code for DpnIII, a third R-M system syntenic to the other Dpn systems. Characterization of DpnIII demonstrated that the endonuclease cleaves unmethylated double-stranded DNA at the tetramer sequence 5' GATC 3', and the cognate methylase is a C5 cytosine-specific DNA methylase. We show that DpnIII decreases the frequency of recombination under in vitro conditions, such that the number of transformants is lower for strains transformed with unmethylated DNA than in those transformed with cognately methylated DNA. Furthermore, we have identified two PMEN1 isolates where the DpnIII endonuclease is disrupted, and phylogenetic work by Croucher and colleagues suggests that these strains have accumulated genomic differences at a higher rate than other PMEN1 strains. We propose that the R-M locus is a major determinant of genetic acquisition; the resident R-M system governs the extent of genome plasticity. IMPORTANCE: Pneumococcus is one of the most important community-acquired bacterial pathogens. Pneumococcal strains can develop resistance to antibiotics and to serotype vaccines by acquiring genes from other strains or species. Thus, genomic plasticity is associated with strain adaptability and pneumococcal success. PMEN1 is a widespread and multidrug-resistant highly pathogenic pneumococcal lineage, which has evolved over the past century and displays a relatively stable genome. In this study, we characterize DpnIII, a restriction-modification (R-M) system that limits recombination. DpnIII is encountered in the PMEN1 lineage, where it replaces other R-M systems that do not decrease plasticity. Our hypothesis is that this genomic region, where different pneumococcal lineages code for variable R-M systems, plays a role in the fine-tuning of the extent of genomic plasticity. It is possible that well-adapted lineages such as PMEN1 have a mechanism to increase genomic stability, rather than foster genomic plasticity.

Novel determinants of antibiotic resistance: identification of mutated loci in highly methicillin-resistant subpopulations of methicillin-resistant Staphylococcus aureus.

We identified mutated genes in highly resistant subpopulations of methicillin-resistant Staphylococcus aureus (MRSA) that are most likely responsible for the historic failure of the ß-lactam family of antibiotics as therapeutic agents against these important pathogens. Such subpopulations are produced during growth of most clinical MRSA strains, including the four historically early MRSA isolates studied here. Chromosomal DNA was prepared from the highly resistant cells along with DNA from the majority of cells (poorly resistant cells) followed by full genome sequencing. In the highly resistant cells, mutations were identified in 3 intergenic sequences and 27 genes representing a wide range of functional categories. A common feature of these mutations appears to be their capacity to induce high-level ß-lactam resistance and increased amounts of the resistance protein PBP2A in the bacteria. The observations fit a recently described model in which the ultimate controlling factor of the phenotypic expression of ß-lactam resistance in MRSA is a RelA-mediated stringent response. IMPORTANCE It has been well established that the level of antibiotic resistance (i.e., minimum concentration of a ß-lactam antibiotic needed to inhibit growth) of a methicillin-resistant Staphylococcus aureus (MRSA) strain depends on the transcription and translation of the resistance protein PBP2A. Here we describe mutated loci in an additional novel set of genetic determinants that appear to be essential for the unusually high resistance levels typical of subpopulations of staphylococci that are produced with unique low frequency in most MRSA clinical isolates. We propose that mutations in these determinants can trigger induction of the stringent stress response which was recently shown to cause increased transcription/translation of the resistance protein PBP2A in parallel with the increased level of resistance.

Molecular tracing of the emergence, diversification, and transmission of S. aureus sequence type 8 in a New York community.

During the last 2 decades, community-associated methicillin-resistant Staphylococcus aureus (CA-MRSA) strains have dramatically increased the global burden of S. aureus infections. The pandemic sequence type (ST)8/pulsed-field gel type USA300 is the dominant CA-MRSA clone in the United States, but its evolutionary history and basis for biological success are incompletely understood. Here, we use whole-genome sequencing of 387 ST8 isolates drawn from an epidemiological network of CA-MRSA infections and colonizations in northern Manhattan to explore short-term evolution and transmission patterns. Phylogenetic analysis predicted that USA300 diverged from a most common recent ancestor around 1993. We found evidence for multiple introductions of USA300 and reconstructed the phylogeographic spread of isolates across neighborhoods. Using pair-wise single-nucleotide polymorphism distances as a measure of genetic relatedness between isolates, we observed that most USA300 isolates had become endemic in households, indicating their critical role as reservoirs for transmission and diversification. Using the maximum single-nucleotide polymorphism variability of isolates from within households as a threshold, we identified several possible transmission networks beyond households. Our study also revealed the evolution of a fluoroquinolone-resistant subpopulation in the mid-1990s and its subsequent expansion at a time of high-frequency outpatient antibiotic use. This high-resolution phylogenetic analysis of ST8 has documented the genomic changes associated with USA300 evolution and how some of its recent evolution has been shaped by antibiotic use. By integrating whole-genome sequencing with detailed epidemiological analyses, our study provides an important framework for delineating the full diversity and spread of USA300 and other emerging pathogens in large urban community populations.

Genomic epidemiology of Neisseria gonorrhoeae with reduced susceptibility to cefixime in the USA: a retrospective observational study.

BACKGROUND: The emergence of Neisseria gonorrhoeae with decreased susceptibility to extended spectrum cephalosporins raises the prospect of untreatable gonorrhoea. In the absence of new treatments, efforts to slow the increasing incidence of resistant gonococcus require insight into the factors that contribute to its emergence and spread. We assessed the relatedness between isolates in the USA and reconstructed likely spread of lineages through different sexual networks. METHODS: We sequenced the genomes of 236 isolates of N gonorrhoeae collected by the Centers for Disease Control and Prevention's Gonococcal Isolate Surveillance Project (GISP) from sentinel public sexually transmitted disease clinics in the USA, including 118 (97%) of the isolates from 2009-10 in GISP with reduced susceptibility to cefixime (cef(RS)) and 118 cefixime-susceptible isolates from GISP matched as closely as possible by location, collection date, and sexual orientation. We assessed the association between antimicrobial resistance genotype and phenotype and correlated phylogenetic clustering with location and sexual orientation. FINDINGS: Mosaic penA XXXIV had a high positive predictive value for cef(RS). We found that two of the 118 cef(RS) isolates lacked a mosaic penA allele, and rechecking showed that these two were susceptible to cefixime. Of the 116 remaining cef(RS) isolates, 114 (98%) fell into two distinct lineages that have independently acquired mosaic penA allele XXXIV. A major lineage of cef(RS) strains spread eastward, predominantly through a sexual network of men who have sex with men. Eight of nine inferred transitions between sexual networks were introductions from men who have sex with men into the heterosexual population. INTERPRETATION: Genomic methods might aid efforts to slow the spread of antibiotic-resistant N gonorrhoeae through augmentation of gonococcal outbreak surveillance and identification of populations that could benefit from increased screening for asymptomatic infections.

Dominant role of nucleotide substitution in the diversification of serotype 3 pneumococci over decades and during a single infection.

Streptococcus pneumoniae of serotype 3 possess a mucoid capsule and cause disease associated with high mortality rates relative to other pneumococci. Phylogenetic analysis of a complete reference genome and 81 draft sequences from clonal complex 180, the predominant serotype 3 clone in much of the world, found most sampled isolates belonged to a clade affected by few diversifying recombinations. However, other isolates indicate significant genetic variation has accumulated over the clonal complex's entire history. Two closely related genomes, one from the blood and another from the cerebrospinal fluid, were obtained from a patient with meningitis. The pair differed in their behaviour in a mouse model of disease and in their susceptibility to antimicrobials, with at least some of these changes attributable to a mutation that up-regulated the patAB efflux pump. This indicates clinically important phenotypic variation can accumulate rapidly through small alterations to the genotype.

A timescale for evolution, population expansion, and spatial spread of an emerging clone of methicillin-resistant Staphylococcus aureus.

Due to the lack of fossil evidence, the timescales of bacterial evolution are largely unknown. The speed with which genetic change accumulates in populations of pathogenic bacteria, however, is a key parameter that is crucial for understanding the emergence of traits such as increased virulence or antibiotic resistance, together with the forces driving pathogen spread. Methicillin-resistant Staphylococcus aureus (MRSA) is a common cause of hospital-acquired infections. We have investigated an MRSA strain (ST225) that is highly prevalent in hospitals in Central Europe. By using mutation discovery at 269 genetic loci (118,804 basepairs) within an international isolate collection, we ascertained extremely low diversity among European ST225 isolates, indicating that a recent population bottleneck had preceded the expansion of this clone. In contrast, US isolates were more divergent, suggesting they represent the ancestral population. While diversity was low, however, our results demonstrate that the short-term evolutionary rate in this natural population of MRSA resulted in the accumulation of measurable DNA sequence variation within two decades, which we could exploit to reconstruct its recent demographic history and the spatiotemporal dynamics of spread. By applying Bayesian coalescent methods on DNA sequences serially sampled through time, we estimated that ST225 had diverged since approximately 1990 (1987 to 1994), and that expansion of the European clade began in 1995 (1991 to 1999), several years before the new clone was recognized. Demographic analysis based on DNA sequence variation indicated a sharp increase of bacterial population size from 2001 to 2004, which is concordant with the reported prevalence of this strain in several European countries. A detailed ancestry-based reconstruction of the spatiotemporal dispersal dynamics suggested a pattern of frequent transmission of the ST225 clone among hospitals within Central Europe. In addition, comparative genomics indicated complex bacteriophage dynamics.

A multilocus phylogeny of Malagasy scincid lizards elucidates the relationships of the fossorial genera Androngo and Cryptoscincus.

A phylogeny for 29 species of scincine lizards from Madagascar, based on 3693 bp of six mitochondrial and five nuclear genes, revealed multiple parallel evolution of adaptations for a burrowing life, and unexpected relationships of the monotypic genera Androngo and Cryptoscincus. Androngo trivittatus was sister to Pygomeles braconnieri, and Cryptoscincus minimus was deeply nested within the genus Paracontias, all of these being fossorial taxa of elongated bodies and partly or fully reduced limbs. To account for these results, we place Cryptoscincus as a junior synonym of Paracontias, and discuss possible taxonomic consequences that may affect the status of Androngo, once additional data become available.