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.

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.

Impact of insertion sequences on convergent evolution of Shigella species.

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.

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.

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.

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.

Salmonella enterica Serovar Hvittingfoss in Bar-Tailed Godwits (Limosa lapponica) from Roebuck Bay, Northwestern Australia.

Salmonella enterica serovar Hvittingfoss is an important foodborne serotype of Salmonella, being detected in many countries where surveillance is conducted. Outbreaks can occur, and there was a recent multistate foodborne outbreak in Australia. S Hvittingfoss can be found in animal populations, though a definitive animal host has not been established. Six species of birds were sampled at Roebuck Bay, a designated Ramsar site in northwestern Australia, resulting in 326 cloacal swabs for bacterial culture. Among a single flock of 63 bar-tailed godwits (Limosa lapponica menzbieri) caught at Wader Spit, Roebuck Bay, in 2018, 17 (27%) were culture positive for Salmonella All other birds were negative for Salmonella The isolates were identified as Salmonella enterica serovar Hvittingfoss. Phylogenetic analysis revealed a close relationship between isolates collected from godwits and the S Hvittingfoss strain responsible for a 2016 multistate foodborne outbreak originating from tainted cantaloupes (rock melons) in Australia. While it is not possible to determine how this strain of S Hvittingfoss was introduced into the bar-tailed godwits, these findings show that wild Australian birds are capable of carrying Salmonella strains of public health importance.IMPORTANCESalmonella is a zoonotic pathogen that causes gastroenteritis and other disease presentations in both humans and animals. Serovars of S. enterica commonly cause foodborne disease in Australia and globally. In 2016-2017, S Hvittingfoss was responsible for an outbreak that resulted in 110 clinically confirmed human cases throughout Australia. The origin of the contamination that led to the outbreak was never definitively established. Here, we identify a migratory shorebird, the bar-tailed godwit, as an animal reservoir of S Hvittingfoss. These birds were sampled in northwestern Australia during their nonbreeding period. The presence of a genetically similar S Hvittingfoss strain circulating in a wild bird population, 2 years after the 2016-2017 outbreak and ∼1,500 km from the suspected source of the outbreak, demonstrates a potentially unidentified environmental reservoir of S Hvittingfoss. While the birds cannot be implicated in the outbreak that occurred 2 years prior, this study does demonstrate the potential role for wild birds in the transmission of this important foodborne pathogen.

Emergence and rapid global dissemination of CTX-M-15-associated Klebsiella pneumoniae strain ST307.

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.

Inferring demographic parameters in bacterial genomic data using Bayesian and hybrid phylogenetic methods.

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 parameters in a fraction of the time required for fully Bayesian analyses. As such, they are a valuable alternative in outbreak studies involving a large number of isolates.

ISMapper: identifying transposase insertion sites in bacterial genomes from short read sequence data.

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 identifying IS insertion sites directly from short read data, with a high degree of accuracy demonstrated across a wide range of bacteria.

Multiple Genetic Mutations Associated with Polymyxin Resistance in Acinetobacter baumannii.

We studied polymyxin B resistance in 10 pairs of clinical Acinetobacter baumannii isolates, two of which had developed polymyxin B resistance in vivo. All polymyxin B-resistant isolates had lower growth rates than and substitution mutations in the lpx or pmrB gene compared to their parent isolates. There were significant differences in terms of antibiotic susceptibility and genetic determinants of resistance in A. baumannii isolates that had developed polymyxin B resistance in vivo compared to isolates that had developed polymyxin B resistance in vitro.

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.

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.

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.

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.

Evidence of microevolution of Salmonella Typhimurium during a series of egg-associated outbreaks linked to a single chicken farm.

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, in place SL1344's P2 phage SopEϕ. DT135 lacks P2 prophage. Two additional plasmids were identified in the S. Typhimurium 135@ isolates, pSTM2 and pSTM7. Both plasmids were IncI1, but phylogenetic analysis of the plasmids and their bacterial hosts shows these plasmids are genetically distinct and result from independent plasmid acquisition events. CONCLUSIONS: This study provides a high-resolution insight into short-term microevolution of the important human pathogen S. Typhimurium. It indicates that purifying selection occurs rapidly in this population (≤ 6 years) and then declines, and provides an estimate for the short-term substitution rate. The latter is likely to be more relevant for foodborne outbreak investigation than previous estimates based on longer time scales.

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.

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.

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.

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.