Molecular characterization and descriptive analysis of carbapenemase-producing Gram-negative rod infections in Bogota, Colombia.

In this study, the genetic differences and clinical impact of the carbapenemase-encoding genes among the community and healthcare-acquired infections were assessed. This retrospective, multicenter cohort study was conducted in Colombia and included patients infected with carbapenem-resistant Gram-negative rods between 2017 and 2021. Carbapenem resistance was identified by Vitek, and carbapenemase-encoding genes were identified by whole-genome sequencing (WGS) to classify the alleles and sequence types (STs). Descriptive statistics were used to determine the association of any pathogen or gene with clinical outcomes. A total of 248 patients were included, of which only 0.8% (2/248) had community-acquired infections. Regarding the identified bacteria, the most prevalent pathogens were Pseudomonas aeruginosa and Klebsiella pneumoniae. In the WGS analysis, 228 isolates passed all the quality criteria and were analyzed. The principal carbapenemase-encoding gene was blaKPC, specifically blaKPC-2 [38.6% (88/228)] and blaKPC-3 [36.4% (83/228)]. These were frequently detected in co-concurrence with blaVIM-2 and blaNDM-1 in healthcare-acquired infections. Notably, the only identified allele among community-acquired infections was blaKPC-3 [50.0% (1/2)]. In reference to the STs, 78 were identified, of which Pseudomonas aeruginosa ST111 was mainly related to blaKPC-3. Klebsiella pneumoniae ST512, ST258, ST14, and ST1082 were exclusively associated with blaKPC-3. Finally, no particular carbapenemase-encoding gene was associated with worse clinical outcomes. The most identified genes in carbapenemase-producing Gram-negative rods were blaKPC-2 and blaKPC-3, both related to gene co-occurrence and diverse STs in the healthcare environment. Patients had several systemic complications and poor clinical outcomes that were not associated with a particular gene.IMPORTANCEAntimicrobial resistance is a pandemic and a worldwide public health problem, especially carbapenem resistance in low- and middle-income countries. Limited data regarding the molecular characteristics and clinical outcomes of patients infected with these bacteria are available. Thus, our study described the carbapenemase-encoding genes among community- and healthcare-acquired infections. Notably, the co-occurrence of carbapenemase-encoding genes was frequently identified. We also found 78 distinct sequence types, of which two were novel Pseudomonas aeruginosa, which could represent challenges in treating these infections. Our study shows that in low and middle-income countries, such as Colombia, the burden of carbapenem resistance in Gram-negative rods is a concern for public health, and regardless of the allele, these infections are associated with poor clinical outcomes. Thus, studies assessing local epidemiology, prevention strategies (including trials), and underpinning genetic mechanisms are urgently needed, especially in low and middle-income countries.

A validated pangenome-scale metabolic model for the Klebsiella pneumoniae species complex.

The Klebsiella pneumoniae species complex (KpSC) is a major source of nosocomial infections globally with high rates of resistance to antimicrobials. Consequently, there is growing interest in understanding virulence factors and their association with cellular metabolic processes for developing novel anti-KpSC therapeutics. Phenotypic assays have revealed metabolic diversity within the KpSC, but metabolism research has been neglected due to experiments being difficult and cost-intensive. Genome-scale metabolic models (GSMMs) represent a rapid and scalable in silico approach for exploring metabolic diversity, which compile genomic and biochemical data to reconstruct the metabolic network of an organism. Here we use a diverse collection of 507 KpSC isolates, including representatives of globally distributed clinically relevant lineages, to construct the most comprehensive KpSC pan-metabolic model to date, KpSC pan v2. Candidate metabolic reactions were identified using gene orthology to known metabolic genes, prior to manual curation via extensive literature and database searches. The final model comprised a total of 3550 reactions, 2403 genes and can simulate growth on 360 unique substrates. We used KpSC pan v2 as a reference to derive strain-specific GSMMs for all 507 KpSC isolates, and compared these to GSMMs generated using a prior KpSC pan-reference (KpSC pan v1) and two single-strain references. We show that KpSC pan v2 includes a greater proportion of accessory reactions (8.8 %) than KpSC pan v1 (2.5 %). GSMMs derived from KpSC pan v2 also generate more accurate growth predictions, with high median accuracies of 95.4 % (aerobic, n=37 isolates) and 78.8 % (anaerobic, n=36 isolates) for 124 matched carbon substrates. KpSC pan v2 is freely available at https://github.com/kelwyres/KpSC-pan-metabolic-model, representing a valuable resource for the scientific community, both as a source of curated metabolic information and as a reference to derive accurate strain-specific GSMMs. The latter can be used to investigate the relationship between KpSC metabolism and traits of interest, such as reservoirs, epidemiology, drug resistance or virulence, and ultimately to inform novel KpSC control strategies.

Genomic investigation of multispecies and multivariant blaNDM outbreak reveals key role of horizontal plasmid transmission.

OBJECTIVES: New Delhi metallo-ß-lactamases (NDMs) are major contributors to the spread of carbapenem resistance globally. In Australia, NDMs were previously associated with international travel, but from 2019 we noted increasing incidence of NDM-positive clinical isolates. We investigated the clinical and genomic epidemiology of NDM carriage at a tertiary-care Australian hospital from 2016 to 2021. METHODS: We identified 49 patients with 84 NDM-carrying isolates in an institutional database, and we collected clinical data from electronic medical record. Short- and long-read whole genome sequencing was performed on all isolates. Completed genome assemblies were used to assess the genetic setting of blaNDM genes and to compare NDM plasmids. RESULTS: Of 49 patients, 38 (78%) were identified in 2019-2021 and only 11 (29%) of 38 reported prior travel, compared with 9 (82%) of 11 in 2016-2018 (P = .037). In patients with NDM infection, the crude 7-day mortality rate was 0% and the 30-day mortality rate was 14% (2 of 14 patients). NDMs were noted in 41 bacterial strains (ie, species and sequence type combinations). Across 13 plasmid groups, 4 NDM variants were detected: blaNDM-1, blaNDM-4, blaNDM-5, and blaNDM-7. We noted a change from a diverse NDM plasmid repertoire in 2016-2018 to the emergence of conserved blaNDM-1 IncN and blaNDM-7 IncX3 epidemic plasmids, with interstrain spread in 2019-2021. These plasmids were noted in 19 (50%) of 38 patients and 35 (51%) of 68 genomes in 2019-2021. CONCLUSIONS: Increased NDM case numbers were due to local circulation of 2 epidemic plasmids with extensive interstrain transfer. Our findings underscore the challenges of outbreak detection when horizontal transmission of plasmids is the primary mode of spread.

High-level nitrofurantoin resistance in a clinical isolate of Klebsiella pneumoniae: a comparative genomics and metabolomics analysis.

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.IMPORTANCEA quest for novel antibiotics and revitalizing older ones (such as nitrofurantoin) for treatment of difficult-to-treat Gram-negative bacterial infections has become increasingly popular. The precise antibacterial activity of nitrofurantoin is still not fully understood. Furthermore, although the prevalence of nitrofurantoin resistance remains low currently, the drug's fast-growing consumption worldwide highlights the need to comprehend the emerging resistance mechanisms. Here, we used multidisciplinary techniques to discern the exact mechanism of nitrofurantoin action and high-level resistance in Klebsiella pneumoniae, a common cause of urinary tract infections for which nitrofurantoin is the recommended treatment. We found that the expression of multiple genes related to membrane transport (including active efflux and passive diffusion of drug molecules) and nitroreductase activity was modified in nitrofurantoin-resistant strains, including oqxR, the transcriptional regulator of the oqxAB efflux pump. Furthermore, complex interconnected metabolic pathways that potentially govern the nitrofurantoin-killing mechanisms (e.g., aminoacyl-tRNA biosynthesis) and nitrofurantoin resistance (riboflavin metabolism) were significantly inhibited following nitrofurantoin treatment. Our study could help inform the improvement of nitrofuran derivatives, the development of new pharmacophores, or drug combinations to support the resurgence of nitrofurantoin in the management of multidrug resistant K. pneumouniae infection.

Bactabolize is a tool for high-throughput generation of bacterial strain-specific metabolic models.

Metabolic capacity can vary substantially within a bacterial species, leading to ecological niche separation, as well as differences in virulence and antimicrobial susceptibility. Genome-scale metabolic models are useful tools for studying the metabolic potential of individuals, and with the rapid expansion of genomic sequencing there is a wealth of data that can be leveraged for comparative analysis. However, there exist few tools to construct strain-specific metabolic models at scale. Here, we describe Bactabolize, a reference-based tool which rapidly produces strain-specific metabolic models and growth phenotype predictions. We describe a pan reference model for the priority antimicrobial-resistant pathogen, Klebsiella pneumoniae, and a quality control framework for using draft genome assemblies as input for Bactabolize. The Bactabolize-derived model for K. pneumoniae reference strain KPPR1 performed comparatively or better than currently available automated approaches CarveMe and gapseq across 507 substrate and 2317 knockout mutant growth predictions. Novel draft genomes passing our systematically defined quality control criteria resulted in models with a high degree of completeness (≥99% genes and reactions captured compared to models derived from matched complete genomes) and high accuracy (mean 0.97, n=10). We anticipate the tools and framework described herein will facilitate large-scale metabolic modelling analyses that broaden our understanding of diversity within bacterial species and inform novel control strategies for priority pathogens.

Genomic dissection of endemic carbapenem resistance reveals metallo-beta-lactamase dissemination through clonal, plasmid and integron transfer.

Infections caused by metallo-beta-lactamase-producing organisms (MBLs) are a global health threat. Our understanding of transmission dynamics and how MBLs establish endemicity remains limited. We analysed two decades of blaIMP-4 evolution in a hospital using sequence data from 270 clinical and environmental isolates (including 169 completed genomes) and identified the blaIMP-4 gene across 7 Gram-negative genera, 68 bacterial strains and 7 distinct plasmid types. We showed how an initial multi-species outbreak of conserved IncC plasmids (95 genomes across 37 strains) allowed endemicity to be established through the ability of blaIMP-4 to disseminate in successful strain-genetic setting pairs we termed propagators, in particular Serratia marcescens and Enterobacter hormaechei. From this reservoir, blaIMP-4 persisted through diversification of genetic settings that resulted from transfer of blaIMP-4 plasmids between bacterial hosts and of the integron carrying blaIMP-4 between plasmids. Our findings provide a framework for understanding endemicity and spread of MBLs and may have broader applicability to other carbapenemase-producing organisms.

DksA is a conserved master regulator of stress response in Acinetobacter baumannii.

Coordination of bacterial stress response mechanisms is critical for long-term survival in harsh environments for successful host infection. The general and specific stress responses of well-studied Gram-negative pathogens like Escherichia coli are controlled by alternative sigma factors, archetypically RpoS. The deadly hospital pathogen Acinetobacter baumannii is notoriously resistant to environmental stresses, yet it lacks RpoS, and the molecular mechanisms driving this incredible stress tolerance remain poorly defined. Here, using functional genomics, we identified the transcriptional regulator DksA as a master regulator for broad stress protection and virulence in A. baumannii. Transcriptomics, phenomics and in vivo animal studies revealed that DksA controls ribosomal protein expression, metabolism, mutation rates, desiccation, antibiotic resistance, and host colonization in a niche-specific manner. Phylogenetically, DksA was highly conserved and well-distributed across Gammaproteobacteria, with 96.6% containing DksA, spanning 88 families. This study lays the groundwork for understanding DksA as a major regulator of general stress response and virulence in this important pathogen.

Within-patient and global evolutionary dynamics of Klebsiella pneumoniae ST17.

Klebsiella pneumoniae sequence type (ST) 17 is a global problem clone that causes multidrug-resistant (MDR) hospital infections worldwide. In 2008-2009, an outbreak of MDR ST17 occurred at a neonatal intensive care unit (NICU) in Stavanger, Norway. Fifty-seven children were colonized. We observed intestinal persistence of ST17 in all of the children for up to two years after hospital discharge. Here, we investigated the within-host evolution of ST17 in 45 of those children during long-term colonization and compared the outbreak with 254 global strains. Ninety-two outbreak-related isolates were whole-genome sequenced. They had capsule locus KL25, O locus O5 and carried yersiniabactin. During within-host colonization ST17 remained stable with few single nucleotide polymorphisms, no acquisition of antimicrobial resistance (AMR) or virulence determinants, and persistent carriage of a bla CTX-M-15-encoding IncFII(K) IncFIB(K) plasmid (pKp2177_1). The global collection included ST17 from 1993 to 2020 from 34 countries, that were from human infection (41.3%), colonization (39.3%) and respiratory specimens (7.3%), from animals (9.3%), and from the environment (2.7%). We estimate that ST17 emerged mid-to-late 19th century (1859, 95 % HPD 1763-1939) and diversified through recombinations of the K and O loci to form several sublineages, with various AMR genes, virulence loci and plasmids. There was limited evidence of persistence of AMR genes in any of these lineages. A globally disseminated sublineage with KL25/O5 accounted for 52.7 % of the genomes. It included a monophyletic subclade that emerged in the mid-1980s, which comprised the Stavanger NICU outbreak and 10 genomes from three other countries, which all carried pKp2177_1. The plasmid was also observed in a KL155/OL101 subclade from the 2000s. Three clonal expansions of ST17 were identified; all were healthcare-associated and carried either yersiniabactin and/or pKp2177_1. To conclude, ST17 is globally disseminated and associated with opportunistic hospital-acquired infections. It contributes to the burden of global MDR infections, but many diverse lineages persist without acquired AMR. We hypothesize that non-human sources and human colonization may play a crucial role for severe infections in vulnerable patients, such as preterm neonates.

A Shigella sonnei clone with extensive drug resistance associated with waterborne outbreaks in China.

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.

ESBL plasmids in Klebsiella pneumoniae: diversity, transmission and contribution to infection burden in the hospital setting.

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 may have reduced the number of opportunities for Plasmid A to transmit and create novel ESBL + strains, as well as reducing overall ESBL infection burden.

Epidemiology and genomic analysis of Klebsiella oxytoca from a single hospital network in Australia.

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.

Attributable mortality and excess length of stay associated with third-generation cephalosporin-resistant Enterobacterales bloodstream infections: a prospective cohort study in Suva, Fiji.

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.

Transmission of Klebsiella strains and plasmids within and between grey-headed flying fox colonies.

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.

Linear plasmids in Klebsiella and other Enterobacteriaceae.

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.

A curated collection of Klebsiella metabolic models reveals variable substrate usage and gene essentiality.

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 interactions constrain adaptation and preserve ecological stability in an experimental microbial 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.

wMel Wolbachia genome remains stable after 7 years in Australian Aedes aegypti field populations.

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.

Pandora: nucleotide-resolution bacterial pan-genomics with reference graphs.

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.

Trycycler: consensus long-read assemblies for bacterial genomes.

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.

Genomic surveillance of antimicrobial resistant bacterial colonisation and infection in intensive care patients.

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 3GCR-GN, ≤3 SNVs for VRE): four VRE, six 3GCR-GN, with epidemiologically linked clusters accounting for 21 and 6% of episodes, respectively (OR 4.3, p = 0.02). CONCLUSIONS: 3GCR-E. coli and VRE were the most common gut colonisers. E. coli was the most common cause of 3GCR-GN infection, but other 3GCR-GN species showed greater risk for infection in colonised patients. Larger studies are warranted to elucidate the relative risks of different colonisers and guide the use of screening in ICU infection control.