Search

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.

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.

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.

Infection of w Mel 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 w Mel introgression into wild-type mosquito populations. The long-term success of this method relies, in part, on evolution of the w Mel genome not compromising the critical features that make it an attractive biocontrol tool. This study compared the w Mel Wolbachia genome at the time of initial releases and 1-7 years post-release in Cairns, Australia. Our results show the w Mel genome remains highly conserved up to 7 years post-release in gene sequence, content, synteny and structure. This work suggests the w Mel genome is stable in its new mosquito host and, therefore, provides reassurance on the potential for w Mel to deliver long-term public-health impacts.

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. These results indicate that most MDR S. enterica serotype Kentucky circulating globally result from the clonal expansion of a single lineage that acquired chromosomal AMR genes 30 years ago, and has continued to diversify and accumulate additional resistances to last-line oral antimicrobials. This article contains data hosted by Microreact.

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.

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.

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.

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.

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.