Search

Nitrofurantoin is a commonly used chemotherapeutic agent in the treatment of uncomplicated urinary tract infections caused by the problematic multidrug resistant Gram-negative pathogen Klebsiella pneumoniae. The present study aims to elucidate the mechanism of nitrofurantoin action and high-level resistance in K. pneumoniae using whole-genome sequencing (WGS), qPCR analysis, mutation structural modeling and untargeted metabolomic analysis. WGS profiling of evolved highly resistant mutants (nitrofurantoin minimum inhibitory concentrations > 256 mg/L) revealed modified expression of several genes related to membrane transport (porin ompK36 and efflux pump regulator oqxR) and nitroreductase activity (ribC and nfsB, involved in nitrofurantoin reduction). Untargeted metabolomics analysis of total metabolites extracted at 1 and 4 h post-nitrofurantoin treatment revealed that exposure to the drug caused a delayed effect on the metabolome which was most pronounced after 4 h. Pathway enrichment analysis illustrated that several complex interrelated metabolic pathways related to nitrofurantoin bacterial killing (aminoacyl-tRNA biosynthesis, purine metabolism, central carbohydrate metabolism, and pantothenate and CoA biosynthesis) and the development of nitrofurantoin resistance (riboflavin metabolism) were significantly perturbed. This study highlights for the first time the key role of efflux pump regulator oqxR in nitrofurantoin resistance and reveals global metabolome perturbations in response to nitrofurantoin, in K. pneumoniae.

Linear plasmids are extrachromosomal DNA elements that have been found in a small number of bacterial species. To date, the only linear plasmids described in the family Enterobacteriaceae belong to Salmonella, first found in Salmonella enterica Typhi. Here, we describe a collection of 12 isolates of the Klebsiella pneumoniae species complex in which we identified linear plasmids. Screening of assembly graphs assembled from public read sets identified linear plasmid structures in a further 13 K. pneumoniae species complex genomes. We used these 25 linear plasmid sequences to query all bacterial genome assemblies in the National Center for Biotechnology Information database, and discovered an additional 61 linear plasmid sequences in a variety of Enterobacteriaceae species. Gene content analysis divided these plasmids into five distinct phylogroups, with very few genes shared across more than two phylogroups. The majority of linear plasmid-encoded genes are of unknown function; however, each phylogroup carried its own unique toxin-antitoxin system and genes with homology to those encoding the ParAB plasmid stability system. Passage in vitro of the 12 linear plasmid-carrying Klebsiella isolates in our collection (which include representatives of all five phylogroups) indicated that these linear plasmids can be stably maintained, and our data suggest they can transmit between K. pneumoniae strains (including members of globally disseminated multidrug-resistant clones) and also between diverse Enterobacteriaceae species. The linear plasmid sequences, and representative isolates harbouring them, are made available as a resource to facilitate future studies on the evolution and function of these novel plasmids.

The Klebsiella pneumoniae species complex (KpSC) is a set of seven Klebsiella taxa that are found in a variety of niches and are an important cause of opportunistic health care-associated infections in humans. Because of increasing rates of multi-drug resistance within the KpSC, there is a growing interest in better understanding the biology and metabolism of these organisms to inform novel control strategies. We collated 37 sequenced KpSC isolates isolated from a variety of niches, representing all seven taxa. We generated strain-specific genome-scale metabolic models (GEMs) for all 37 isolates and simulated growth phenotypes on 511 distinct carbon, nitrogen, sulfur, and phosphorus substrates. Models were curated and their accuracy was assessed using matched phenotypic growth data for 94 substrates (median accuracy of 96%). We explored species-specific growth capabilities and examined the impact of all possible single gene deletions using growth simulations in 145 core carbon substrates. These analyses revealed multiple strain-specific differences, within and between species, and highlight the importance of selecting a diverse range of strains when exploring KpSC metabolism. This diverse set of highly accurate GEMs could be used to inform novel drug design, enhance genomic analyses, and identify novel virulence and resistance determinants. We envisage that these 37 curated strain-specific GEMs, covering all seven taxa of the KpSC, provide a valuable resource to the Klebsiella research community.

Species loss within a microbial community can increase resource availability and spur adaptive evolution. Environmental shifts that cause species loss or fluctuations in community composition are expected to become more common, so it is important to understand the evolutionary forces that shape the stability and function of the emergent community. Here we study experimental cultures of a simple, ecologically stable community of Saccharomyces cerevisiae and Lactobacillus plantarum, in order to understand how the presence or absence of a species impacts coexistence over evolutionary timescales. We found that evolution in coculture led to drastically altered evolutionary outcomes for L. plantarum, but not S. cerevisiae. Both monoculture- and co-culture-evolved L. plantarum evolved dozens of mutations over 925 generations of evolution, but only L. plantarum that had evolved in isolation from S. cerevisiae lost the capacity to coexist with S. cerevisiae. We find that the evolutionary loss of ecological stability corresponds with fitness differences between monoculture-evolved L. plantarum and S. cerevisiae and genetic changes that repeatedly evolve across the replicate populations of L. plantarum. This work shows how coevolution within a community can prevent destabilising evolution in individual species, thereby preserving ecological diversity and stability, despite rapid adaptation.

BACKGROUND: Infections caused by Klebsiella oxytoca are the second most common cause of Klebsiella infections in humans. Most studies have focused on K. oxytoca outbreaks and few have examined the broader clinical context of K. oxytoca. METHODS: Here, we collected all clinical isolates identified as K. oxytoca in a hospital microbiological diagnostic lab across a 15-month period (n = 239). Whole genome sequencing was performed on a subset of 92 isolates (all invasive, third-generation cephalosporin resistant (3GCR) and non-urinary isolates collected > 48 h after admission), including long-read sequencing on a further six isolates with extended-spectrum beta-lactamase or carbapenemase genes. RESULTS: The majority of isolates were sensitive to antimicrobials, however 22 isolates were 3GCR, of which five were also carbapenem resistant. Genomic analyses showed those identified as K. oxytoca by the clinical laboratory actually encompassed four distinct species (K. oxytoca, Klebsiella michiganensis, Klebsiella grimontii and Klebsiella pasteurii), referred to as the K. oxytoca species complex (KoSC). There was significant diversity within the population, with only 10/67 multi-locus sequence types (STs) represented by more than one isolate. Strain transmission was rare, with only one likely event identified. Six isolates had extended spectrum beta-lactamase (blaSHV-12 and/or blaCTX-M-9) or carbapenemase (blaIMP-4) genes. One pair of K. michiganensis and K. pasteurii genomes carried identical blaIMP-4 IncL/M plasmids, indicative of plasmid transmission. CONCLUSION: Whilst antimicrobial resistance was rare, the resistance plasmids were similar to those found in other Enterobacterales, demonstrating that KoSC has access to the same plasmid reservoir and thus there is potential for multi-drug resistance. Further genomic studies are required to improve our understanding of the KoSC population and facilitate investigation into the attributes of successful nosocomial isolates.

BACKGROUND: Resistance to third-generation cephalosporins, often mediated by extended-spectrum beta-lactamases (ESBLs), is a considerable issue in hospital-associated infections as few drugs remain for treatment. ESBL genes are often located on large plasmids that transfer horizontally between strains and species of Enterobacteriaceae and frequently confer resistance to additional drug classes. Whilst plasmid transmission is recognised to occur in the hospital setting, the frequency and impact of plasmid transmission on infection burden, compared to ESBL + strain transmission, is not well understood. METHODS: We sequenced the genomes of clinical and carriage isolates of Klebsiella pneumoniae species complex from a year-long hospital surveillance study to investigate ESBL burden and plasmid transmission in an Australian hospital. Long-term persistence of a key transmitted ESBL + plasmid was investigated via sequencing of ceftriaxone-resistant isolates during 4 years of follow-up, beginning 3 years after the initial study. RESULTS: We found 25 distinct ESBL plasmids. We identified one plasmid, which we called Plasmid A, that carried blaCTX-M-15 in an IncF backbone similar to pKPN-307. Plasmid A was transmitted at least four times into different Klebsiella species/lineages and was responsible for half of all ESBL episodes during the initial 1-year study period. Three of the Plasmid A-positive strains persisted locally 3-6 years later, and Plasmid A was detected in two additional strain backgrounds. Overall Plasmid A accounted for 21% of ESBL + infections in the follow-up period. CONCLUSIONS: Here, we systematically surveyed ESBL strain and plasmid transmission over 1 year in a single hospital network. Whilst ESBL plasmid transmission events were rare in this setting, they had a significant and sustained impact on the burden of ceftriaxone-resistant and multidrug-resistant infections. If onward transmission of Plasmid A-carrying strains could have been prevented, this m

OBJECTIVES: There are scant primary clinical data on antimicrobial resistance (AMR) burden from low- and middle-income countries (LMICs). We adapted recent World Health Organization methodology to measure the effect of third-generation cephalosporin resistance (3GC-R) on mortality and excess length of hospital stay in Fiji. METHODS: We conducted a prospective cohort study of inpatients with Enterobacterales bloodstream infections (BSIs) at Colonial War Memorial Hospital, Suva. We used cause-specific Cox proportional hazards models to estimate the effect of 3GC-R on the daily risk (hazard) of in-hospital mortality and being discharged alive (competing risks), and we used multistate modelling to estimate the excess length of hospital stay. RESULTS: From July 2020 to February 2021 we identified 162 consecutive Enterobacterales BSIs; 3GC-R was present in 66 (40.7%). Crude mortality for patients with 3GC-susceptible and 3GC-R BSIs was 16.7% (16/96) and 30.3% (20/66), respectively. 3GC-R was not associated with the in-hospital mortality hazard rate (adjusted hazard ratio [aHR] 1.13, 95% confidence interval [CI] 0.51-2.53) or being discharged alive (aHR 0.99, 95% CI 0.65-1.50), whereas Charlson comorbidity index score (aHR 1.62, 95% CI 1.36-1.93) and Pitt bacteraemia score (aHR 3.57, 95% CI 1.31-9.71) were both associated with an increased hazard rate of in-hospital mortality. 3GC-R was associated with an increased length of stay of 2.6 days (95% CI 2.5-2.8). 3GC-R was more common among hospital-associated infections, but genomics did not identify clonal transmission. CONCLUSION: Patients with Enterobacterales BSIs in Fiji had high mortality. There were high rates of 3GC-R, which was associated with increased hospital length of stay but not with in-hospital mortality.

Antimicrobial resistance of Shigella sonnei has become a global concern. Here, we report a phylogenetic group of S. sonnei with extensive drug resistance, including a combination of multidrug resistance, coresistance to ceftriaxone and azithromycin (cefRaziR), reduced susceptibility to fluoroquinolones, and even colistin resistance (colR). This distinct clone caused six waterborne shigellosis outbreaks in China from 2015 to 2020. We collect 155 outbreak isolates and 152 sporadic isolates. The cefRaziR isolates, including outbreak strains, are mainly distributed in a distinct clade located in global Lineage III. The outbreak strains form a recently derived monophyletic group that may have emerged circa 2010. The cefRaziR and colR phenotypes are attributed to the acquisition of different plasmids, particularly the IncB/O/K/Z plasmid coharboring the blaCTX-M-14, mphA, aac(3)-IId, dfrA17, aadA5, and sul1 genes and the IncI2 plasmid with an mcr-1 gene. Genetic analyses identify 92 accessory genes and 60 single-nucleotide polymorphisms associated with the cefRaziR phenotype. Surveillance of this clone is required to determine its dissemination and threat to global public health.

The grey-headed flying fox (Pteropus poliocephalus) is an endemic Australian fruit bat, known to carry zoonotic pathogens. We recently showed they harbour bacterial pathogen Klebsiella pneumoniae and closely related species in the K. pneumoniae species complex (KpSC); however, the dynamics of KpSC transmission and gene flow within flying fox colonies are poorly understood. High-resolution genome comparisons of 39 KpSC isolates from grey-headed flying foxes identified five putative strain transmission clusters (four intra- and one inter-colony). The instance of inter-colony strain transmission of K. africana was found between two flying fox populations within flying distance, indicating either direct or indirect transmission through a common food/water source. All 11 plasmids identified within the KpSC isolates showed 73% coverage (mean) and ≥95% identity to human-associated KpSC plasmids, indicating gene flow between human clinical and grey-headed flying fox isolates. Along with strain transmission, inter-species horizontal plasmid transmission between K. pneumoniae and Klebsiella africana was also identified within a flying fox colony. Finally, genome-scale metabolic models were generated to predict and compare substrate usage to previously published KpSC models, from human and environmental sources. These models indicated no distinction on the basis of metabolic capabilities. Instead, metabolic capabilities were consistent with population structure and ST/lineage.

We present pandora, a novel pan-genome graph structure and algorithms for identifying variants across the full bacterial pan-genome. As much bacterial adaptability hinges on the accessory genome, methods which analyze SNPs in just the core genome have unsatisfactory limitations. Pandora approximates a sequenced genome as a recombinant of references, detects novel variation and pan-genotypes multiple samples. Using a reference graph of 578 Escherichia coli genomes, we compare 20 diverse isolates. Pandora recovers more rare SNPs than single-reference-based tools, is significantly better than picking the closest RefSeq reference, and provides a stable framework for analyzing diverse samples without reference bias.

Infection of wMel Wolbachia in Aedes aegypti imparts two signature features that enable its application for biocontrol of dengue. First, the susceptibility of mosquitoes to viruses such as dengue and Zika is reduced. Second, a reproductive manipulation is caused that enables wMel introgression into wild-type mosquito populations. The long-term success of this method relies, in part, on evolution of the wMel genome not compromising the critical features that make it an attractive biocontrol tool. This study compared the wMel Wolbachia genome at the time of initial releases and 1-7 years post-release in Cairns, Australia. Our results show the wMel genome remains highly conserved up to 7 years post-release in gene sequence, content, synteny and structure. This work suggests the wMel genome is stable in its new mosquito host and, therefore, provides reassurance on the potential for wMel to deliver long-term public-health impacts.

Shigella sonnei is the most common agent of shigellosis in high-income countries, and causes a significant disease burden in low- and middle-income countries. Antimicrobial resistance is increasingly common in all settings. Whole genome sequencing (WGS) is increasingly utilised for S. sonnei outbreak investigation and surveillance, but comparison of data between studies and labs is challenging. Here, we present a genomic framework and genotyping scheme for S. sonnei to efficiently identify genotype and resistance determinants from WGS data. The scheme is implemented in the software package Mykrobe and tested on thousands of genomes. Applying this approach to analyse >4,000 S. sonnei isolates sequenced in public health labs in three countries identified several common genotypes associated with increased rates of ciprofloxacin resistance and azithromycin resistance, confirming intercontinental spread of highly-resistant S. sonnei clones and demonstrating the genomic framework can facilitate monitoring the spread of resistant clones, including those that have recently emerged, at local and global scales.

BACKGROUND: Third-generation cephalosporin-resistant Gram-negatives (3GCR-GN) and vancomycin-resistant enterococci (VRE) are common causes of multi-drug resistant healthcare-associated infections, for which gut colonisation is considered a prerequisite. However, there remains a key knowledge gap about colonisation and infection dynamics in high-risk settings such as the intensive care unit (ICU), thus hampering infection prevention efforts. METHODS: We performed a three-month prospective genomic survey of infecting and gut-colonising 3GCR-GN and VRE among patients admitted to an Australian ICU. Bacteria were isolated from rectal swabs (n = 287 and n = 103 patients ≤2 and > 2 days from admission, respectively) and diagnostic clinical specimens between Dec 2013 and March 2014. Isolates were subjected to Illumina whole-genome sequencing (n = 127 3GCR-GN, n = 41 VRE). Multi-locus sequence types (STs) and antimicrobial resistance determinants were identified from de novo assemblies. Twenty-three isolates were selected for sequencing on the Oxford Nanopore MinION device to generate completed reference genomes (one for each ST isolated from ≥2 patients). Single nucleotide variants (SNVs) were identified by read mapping and variant calling against these references. RESULTS: Among 287 patients screened on admission, 17.4 and 8.4% were colonised by 3GCR-GN and VRE, respectively. Escherichia coli was the most common species (n = 36 episodes, 58.1%) and the most common cause of 3GCR-GN infection. Only two VRE infections were identified. The rate of infection among patients colonised with E. coli was low, but higher than those who were not colonised on admission (n = 2/33, 6% vs n = 4/254, 2%, respectively, p = 0.3). While few patients were colonised with 3GCR- Klebsiella pneumoniae or Pseudomonas aeruginosa on admission (n = 4), all such patients developed infections with the colonising strain. Genomic analyses revealed 10 putative nosocomial transmission clusters (≤20 SNVs for

While long-read sequencing allows for the complete assembly of bacterial genomes, long-read assemblies contain a variety of errors. Here, we present Trycycler, a tool which produces a consensus assembly from multiple input assemblies of the same genome. Benchmarking showed that Trycycler assemblies contained fewer errors than assemblies constructed with a single tool. Post-assembly polishing further reduced errors and Trycycler+polishing assemblies were the most accurate genomes in our study. As Trycycler requires manual intervention, its output is not deterministic. However, we demonstrated that multiple users converge on similar assemblies that are consistently more accurate than those produced by automated assembly tools.

Shigella species are specialised lineages of Escherichia coli that have converged to become human-adapted and cause dysentery by invading human gut epithelial cells. Most studies of Shigella evolution have been restricted to comparisons of single representatives of each species; and population genomic studies of individual Shigella species have focused on genomic variation caused by single nucleotide variants and ignored the contribution of insertion sequences (IS) which are highly prevalent in Shigella genomes. Here, we investigate the distribution and evolutionary dynamics of IS within populations of Shigella dysenteriae Sd1, Shigella sonnei and Shigella flexneri. We find that five IS (IS1, IS2, IS4, IS600 and IS911) have undergone expansion in all Shigella species, creating substantial strain-to-strain variation within each population and contributing to convergent patterns of functional gene loss within and between species. We find that IS expansion and genome degradation are most advanced in S. dysenteriae and least advanced in S. sonnei; and using genome-scale models of metabolism we show that Shigella species display convergent loss of core E. coli metabolic capabilities, with S. sonnei and S. flexneri following a similar trajectory of metabolic streamlining to that of S. dysenteriae. This study highlights the importance of IS to the evolution of Shigella and provides a framework for the investigation of IS dynamics and metabolic reduction in other bacterial species.

Knowledge of mobile genetic elements that capture and disseminate antimicrobial resistance genes between diverse environments, particularly across human-animal boundaries, is key to understanding the role anthropogenic activities have in the evolution of antimicrobial resistance. Plasmids that circulate within the Enterobacteriaceae and the Proteobacteria more broadly are well placed to acquire resistance genes sourced from separate niche environments and provide a platform for smaller mobile elements such as IS26 to assemble these genes into large, complex genomic structures. Here, we characterised two atypical Z/I1 hybrid plasmids, pSTM32-108 and pSTM37-118, hosting antimicrobial resistance and virulence associated genes within endemic pathogen Salmonella enterica serovar Typhimurium 1,4,[5],12:i:-, sourced from Australian swine production facilities during 2013. We showed that the plasmids found in S. Typhimurium 1,4,[5],12:i:- are close relatives of two plasmids identified from Escherichia coli of human and bovine origin in Australia circa 1998. The older plasmids, pO26-CRL125 and pO111-CRL115, encoded a putative serine protease autotransporter and were host to a complex resistance region composed of a hybrid Tn21-Tn1721 mercury resistance transposon and composite IS26 transposon Tn6026. This gave a broad antimicrobial resistance profile keyed towards first generation antimicrobials used in Australian agriculture but also included a class 1 integron hosting the trimethoprim resistance gene dfrA5. Genes encoding resistance to ampicillin, trimethoprim, sulphonamides, streptomycin, aminoglycosides, tetracyclines and mercury were a feature of these plasmids. Phylogenetic analyses showed very little genetic drift in the sequences of these plasmids over the past 15 years; however, some alterations within the complex resistance regions present on each plasmid have led to the loss of various resistance genes, presumably as a result of the activity of IS26. These alterat

Resistance to carbapenem and aminoglycoside antibiotics is a critical problem in Acinetobacter baumannii, particularly when genes conferring resistance are acquired by multiply or extensively resistant members of successful globally distributed clonal complexes, such as global clone 1 (GC1) . Here, we investigate the evolution of an expanding clade of lineage 1 of the GC1 complex via repeated acquisition of carbapenem- and aminoglycoside-resistance genes. Lineage 1 arose in the late 1970s and the Tn6168/OCL3 clade arose in the late 1990s from an ancestor that had already acquired resistance to third-generation cephalosporins and fluoroquinolones. Between 2000 and 2002, two distinct subclades have emerged, and they are distinguishable via the presence of an integrated phage genome in subclade 1 and AbaR4 (carrying the oxa23 carbapenem-resistance gene in Tn2006) at a specific chromosomal location in subclade 2. Part or all of the original resistance gene cluster in the chromosomally located AbaR3 has been lost from some isolates, but plasmids carrying alternate resistance genes have been gained. In one group in subclade 2, the chromosomally located AbGRI3, carrying the armA aminoglycoside-resistance gene, has been acquired from a GC2 isolate and incorporated via homologous recombination. ISAba1 entered the common ancestor of this clade as part of the cephalosporin-resistance transposon Tn6168 and has dispersed differently in each subclade. Members of subclade 1 share an ISAba1 in one specific position in the chromosome and in subclade 2 two different ISAba1 locations are shared. Further shared ISAba1 locations distinguish further divisions, potentially providing simple markers for epidemiological studies.

OBJECTIVES: Recent reports indicate the emergence of a new carbapenemase-producing Klebsiella pneumoniae clone, ST307. We sought to better understand the global epidemiology and evolution of this clone and evaluate its association with antimicrobial resistance (AMR) genes. METHODS: We collated information from the literature and public databases and performed a comparative analysis of 95 ST307 genomes (including 37 that were newly sequenced). RESULTS: We show that ST307 emerged in the mid-1990s (nearly 20 years prior to its first report), is already globally distributed and is intimately associated with a conserved plasmid harbouring the blaCTX-M-15 ESBL gene and several other AMR determinants. CONCLUSIONS: Our findings support the need for enhanced surveillance of this widespread ESBL clone in which carbapenem resistance has occasionally emerged.

Salmonella enterica serotype Kentucky can be a common causative agent of salmonellosis, usually associated with consumption of contaminated poultry. Antimicrobial resistance (AMR) to multiple drugs, including ciprofloxacin, is an emerging problem within this serotype. We used whole-genome sequencing (WGS) to investigate the phylogenetic structure and AMR content of 121 S.enterica serotype Kentucky sequence type 198 isolates from five continents. Population structure was inferred using phylogenomic analysis and whole genomes were compared to investigate changes in gene content, with a focus on acquired AMR genes. Our analysis showed that multidrug-resistant (MDR) S.enterica serotype Kentucky isolates belonged to a single lineage, which we estimate emerged circa 1989 following the acquisition of the AMR-associated Salmonella genomic island (SGI) 1 (variant SGI1-K) conferring resistance to ampicillin, streptomycin, gentamicin, sulfamethoxazole and tetracycline. Phylogeographical analysis indicates this clone emerged in Egypt before disseminating into Northern, Southern and Western Africa, then to the Middle East, Asia and the European Union. The MDR clone has since accumulated various substitution mutations in the quinolone-resistance-determining regions (QRDRs) of DNA gyrase (gyrA) and DNA topoisomerase IV (parC), such that most strains carry three QRDR mutations which together confer resistance to ciprofloxacin. The majority of AMR genes in the S. enterica serotype Kentucky genomes were carried either on plasmids or SGI structures. Remarkably, each genome of the MDR clone carried a different SGI1-K derivative structure; this variation could be attributed to IS26-mediated insertions and deletions, which appear to have hampered previous attempts to trace the clone's evolution using sub-WGS resolution approaches. Several different AMR plasmids were also identified, encoding resistance to chloramphenicol, third-generation cephalosporins, carbapenems and/or azithromycin.

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.

Acinetobacter baumannii is a common causative agent of hospital-acquired infections and a leading cause of infection in burns patients. Carbapenem-resistant A. baumannii is considered a major public-health threat and has been identified by the World Health Organization as the top priority organism requiring new antimicrobials. The most common mechanism for carbapenem resistance in A. baumannii is via horizontal acquisition of carbapenemase genes. In this study, we sampled 20 A. baumannii isolates from a patient with extensive burns, and characterized the evolution of carbapenem resistance over a 45 day period via Illumina and Oxford Nanopore sequencing. All isolates were multidrug resistant, carrying two genomic islands that harboured several antibiotic-resistance genes. Most isolates were genetically identical and represented a single founder genotype. We identified three novel non-synonymous substitutions associated with meropenem resistance: F136L and G288S in AdeB (part of the AdeABC efflux pump) associated with an increase in meropenem MIC to ≥8 µg ml-1; and A515V in FtsI (PBP3, a penicillin-binding protein) associated with a further increase in MIC to 32 µg ml-1. Structural modelling of AdeB and FtsI showed that these mutations affected their drug-binding sites and revealed mechanisms for meropenem resistance. Notably, one of the adeB mutations arose prior to meropenem therapy but following ciprofloxacin therapy, suggesting exposure to one drug whose resistance is mediated by the efflux pump can induce collateral resistance to other drugs to which the bacterium has not yet been exposed.

BACKGROUND: Recent developments in sequencing technologies make it possible to obtain genome sequences from a large number of isolates in a very short time. Bayesian phylogenetic approaches can take advantage of these data by simultaneously inferring the phylogenetic tree, evolutionary timescale, and demographic parameters (such as population growth rates), while naturally integrating uncertainty in all parameters. Despite their desirable properties, Bayesian approaches can be computationally intensive, hindering their use for outbreak investigations involving genome data for a large numbers of pathogen isolates. An alternative to using full Bayesian inference is to use a hybrid approach, where the phylogenetic tree and evolutionary timescale are estimated first using maximum likelihood. Under this hybrid approach, demographic parameters are inferred from estimated trees instead of the sequence data, using maximum likelihood, Bayesian inference, or approximate Bayesian computation. This can vastly reduce the computational burden, but has the disadvantage of ignoring the uncertainty in the phylogenetic tree and evolutionary timescale. RESULTS: We compared the performance of a fully Bayesian and a hybrid method by analysing six whole-genome SNP data sets from a range of bacteria and simulations. The estimates from the two methods were very similar, suggesting that the hybrid method is a valid alternative for very large datasets. However, we also found that congruence between these methods is contingent on the presence of strong temporal structure in the data (i.e. clocklike behaviour), which is typically verified using a date-randomisation test in a Bayesian framework. To reduce the computational burden of this Bayesian test we implemented a date-randomisation test using a rapid maximum likelihood method, which has similar performance to its Bayesian counterpart. CONCLUSIONS: Hybrid approaches can produce reliable inferences of evolutionary timescales and phylodynamic

Recent reports indicate the emergence of a new carbapenemase producing Klebsiella pneumoniae clone, ST307. Here we show that ST307 emerged in the mid-1990s (nearly 20 years prior to its first report), is already globally distributed and is intimately associated with a conserved plasmid harbouring the bla CTX-M-15 extended-spectrum beta-lactamase (ESBL) gene plus other antimicrobial resistance determinants. Our findings support the need for enhanced surveillance of this widespread ESBL clone in which carbapenem resistance is now emerging.

We recently reported a dramatic increase in the prevalence of carbapenem-resistant Acinetobacter baumannii infections in the intensive care unit (ICU) of a Vietnamese hospital. This upsurge was associated with a specific oxa23-positive clone that was identified by multilocus VNTR analysis. Here, we used whole-genome sequence analysis to dissect the emergence of carbapenem-resistant A. baumannii causing ventilator-associated pneumonia (VAP) in the ICU during 2009-2012. To provide historical context and distinguish microevolution from strain introduction, we compared these genomes with those of A. baumannii asymptomatic carriage and VAP isolates from this same ICU collected during 2003-2007. We identified diverse lineages co-circulating over many years. Carbapenem resistance was associated with the presence of oxa23, oxa40, oxa58 and ndm1 genes in multiple lineages. The majority of resistant isolates were oxa23-positive global clone GC2; fine-scale phylogenomic analysis revealed five distinct GC2 sublineages within the ICU that had evolved locally via independent chromosomal insertions of oxa23 transposons. The increase in infections caused by carbapenem-resistant A. baumannii was associated with transposon-mediated transmission of a carbapenemase gene, rather than clonal expansion or spread of a carbapenemase-harbouring plasmid. Additionally, we found evidence of homologous recombination creating diversity within the local GC2 population, including several events resulting in replacement of the capsule locus. We identified likely donors of the imported capsule locus sequences amongst the A. baumannii isolated on the same ward, suggesting that diversification was largely facilitated via reassortment and sharing of genetic material within the localized A. baumannii population.

Estimating the rates at which bacterial genomes evolve is critical to understanding major evolutionary and ecological processes such as disease emergence, long-term host-pathogen associations and short-term transmission patterns. The surge in bacterial genomic data sets provides a new opportunity to estimate these rates and reveal the factors that shape bacterial evolutionary dynamics. For many organisms estimates of evolutionary rate display an inverse association with the time-scale over which the data are sampled. However, this relationship remains unexplored in bacteria due to the difficulty in estimating genome-wide evolutionary rates, which are impacted by the extent of temporal structure in the data and the prevalence of recombination. We collected 36 whole genome sequence data sets from 16 species of bacterial pathogens to systematically estimate and compare their evolutionary rates and assess the extent of temporal structure in the absence of recombination. The majority (28/36) of data sets possessed sufficient clock-like structure to robustly estimate evolutionary rates. However, in some species reliable estimates were not possible even with 'ancient DNA' data sampled over many centuries, suggesting that they evolve very slowly or that they display extensive rate variation among lineages. The robustly estimated evolutionary rates spanned several orders of magnitude, from approximately 10-5 to 10-8 nucleotide substitutions per site year-1. This variation was negatively associated with sampling time, with this relationship best described by an exponential decay curve. To avoid potential estimation biases, such time-dependency should be considered when inferring evolutionary time-scales in bacteria.

Multiply antibiotic-resistant Acinetobacter baumannii isolate D36 was recovered in Australia in 2008 and belongs to a distinct lineage of global clone 1 (GC1). Here, we present the complete 4.13 Mbp genome sequence (chromosome plus 4 plasmids), generated via long read sequencing (PacBio).

BACKGROUND: Insertion sequences (IS) are small transposable elements, commonly found in bacterial genomes. Identifying the location of IS in bacterial genomes can be useful for a variety of purposes including epidemiological tracking and predicting antibiotic resistance. However IS are commonly present in multiple copies in a single genome, which complicates genome assembly and the identification of IS insertion sites. Here we present ISMapper, a mapping-based tool for identification of the site and orientation of IS insertions in bacterial genomes, directly from paired-end short read data. RESULTS: ISMapper was validated using three types of short read data: (i) simulated reads from a variety of species, (ii) Illumina reads from 5 isolates for which finished genome sequences were available for comparison, and (iii) Illumina reads from 7 Acinetobacter baumannii isolates for which predicted IS locations were tested using PCR. A total of 20 genomes, including 13 species and 32 distinct IS, were used for validation. ISMapper correctly identified 97 % of known IS insertions in the analysis of simulated reads, and 98 % in real Illumina reads. Subsampling of real Illumina reads to lower depths indicated ISMapper was able to correctly detect insertions for average genome-wide read depths >20x, although read depths >50x were required to obtain confident calls that were highly-supported by evidence from reads. All ISAba1 insertions identified by ISMapper in the A. baumannii genomes were confirmed by PCR. In each A. baumannii genome, ISMapper successfully identified an IS insertion upstream of the ampC beta-lactamase that could explain phenotypic resistance to third-generation cephalosporins. The utility of ISMapper was further demonstrated by profiling genome-wide IS6110 insertions in 138 publicly available Mycobacterium tuberculosis genomes, revealing lineage-specific insertions and multiple insertion hotspots. CONCLUSIONS: ISMapper provides a rapid and robust method for id

Acinetobacter baumannii isolate A1 was recovered in the United Kingdom in 1982 and belongs to global clone 1 (GC1). Here, we present its complete 3.91-Mbp genome sequence, generated via a combination of short-read sequencing (Illumina), long-read sequencing (PacBio), and manual finishing.

Multiply antibiotic-resistant Acinetobacter baumannii isolate D36 was recovered in Australia in 2008 and belongs to a distinct lineage of global clone 1 (GC1). Here, we present the complete 4.13 Mbp genome sequence (chromosome plus 4 plasmids), generated via long read sequencing (PacBio).

BACKGROUND: The bacterium Salmonella enterica serovar Typhimurium (S. Typhimurium) is one of the most frequent causes of foodborne outbreaks of gastroenteritis. Between 2005-2008 a series of S. Typhimurium outbreaks occurred in Tasmania, Australia, that were all traced to eggs originating from a single chicken farm. We sequenced the genomes of 12 isolates linked to these outbreaks, in order to investigate the microevolution of a pathogenic S. Typhimurium clone in a natural, spatiotemporally restricted population. RESULTS: The isolates, which shared a phage type similar to DT135 known locally as 135@ or 135a, formed a clade within the S. Typhimurium population with close similarity to the reference genome SL1334 (160 single nucleotide polymorphisms, or SNPs). Ten of the isolates belonged to a single clone (<23 SNPs between isolate pairs) which likely represents the population of S. Typhimurium circulating at the chicken farm; the other two were from sporadic cases and were genetically distinct from this clone. Divergence dating indicated that all 12 isolates diverged from a common ancestor in the mid 1990 s, and the clone began to diversify in 2003-2004. This clone spilled out into the human population several times between 2005-2008, during which time it continued to accumulate SNPs at a constant rate of 3-5 SNPs per year or 1x10-6 substitutions site-1 year-1, faster than the longer-term (~50 year) rates estimated previously for S. Typhimurium. Our data suggest that roughly half of non-synonymous substitutions are rapidly removed from the S. Typhimurium population, after which purifying selection is no longer important and the remaining substitutions become fixed in the population. The S. Typhimurium 135@ isolates were nearly identical to SL1344 in terms of gene content and virulence plasmids. Their phage contents were close to SL1344, except that they carried a different variant of Gifsy-1, lacked the P2 remnant found in SL1344 and carried a novel P2 phage, P2-Hawk