Ten-year longitudinal molecular epidemiology study of Escherichia coli and Klebsiella species bloodstream infections in Oxfordshire, UK

Longitudinal analysis within one hospital in sub-Saharan Africa over 20 years reveals repeated replacements of dominant clones of Klebsiella pneumoniae and stresses the importance to include temporal patterns for vaccine design considerations

BACKGROUND

Infections caused by multidrug-resistant gram-negative bacteria present a severe threat to global public health. The WHO defines drug-resistant Klebsiella pneumoniae as a priority pathogen for which alternative treatments are needed given the limited treatment options and the rapid acquisition of novel resistance mechanisms by this species. Longitudinal descriptions of genomic epidemiology of Klebsiella pneumoniae can inform management strategies but data from sub-Saharan Africa are lacking.

METHODS

We present a longitudinal analysis of all invasive K. pneumoniae isolates from a single hospital in Blantyre, Malawi, southern Africa, from 1998 to 2020, combining clinical data with genome sequence analysis of the isolates.

RESULTS

We show that after a dramatic increase in the number of infections from 2016 K. pneumoniae becomes hyperendemic, driven by an increase in neonatal infections. Genomic data show repeated waves of clonal expansion of different, often ward-restricted, lineages, suggestive of hospital-associated transmission. We describe temporal trends in resistance and surface antigens, of relevance for vaccine development.

CONCLUSIONS

Our data highlight a clear need for new interventions to prevent rather than treat K. pneumoniae infections in our setting. Whilst one option may be a vaccine, the majority of cases could be avoided by an increased focus on and investment in infection prevention and control measures, which would reduce all healthcare-associated infections and not just one.

Risk of death in Klebsiella pneumoniae bloodstream infections is associated with specific phylogenetic lineages

BACKGROUND

Klebsiella pneumoniae species complex (KpSC) bloodstream infections (BSIs) are associated with considerable morbidity and mortality, particularly in elderly and multimorbid patients. Multidrug-resistant (MDR) strains have been associated with poorer outcome. However, the clinical impact of KpSC phylogenetic lineages on BSI outcome is unclear.

METHODS

In an 18-month nationwide Norwegian prospective study of KpSC BSI episodes in adults, we used whole-genome sequencing to describe the molecular epidemiology of KpSC, and multivariable Cox regression analysis including clinical data to determine adjusted hazard ratios (aHR) for death associated with specific genomic lineages.

FINDINGS

We included 1078 BSI episodes and 1082 bacterial isolates from 1055 patients. The overall 30-day case-fatality rate (CFR) was 12.5%. Median patient age was 73.4, 61.7% of patients were male. Median Charlson comorbidity score was 3. Klebsiella pneumoniae sensu stricto (Kp) (79.3%, n = 858/1082) and K. variicola (15.7%, n = 170/1082) were the dominating phylogroups. Global MDR-associated Kp clonal groups (CGs) were prevalent (25.0%, n = 270/1082) but 78.9% (n = 213/270) were not MDR, and 53.7% (n = 145/270) were community acquired. The major findings were increased risk for death within 30 days in monomicrobial BSIs caused by K. variicola (CFR 16.9%, n = 21; aHR 1.86, CI 1.10-3.17, p = 0.02), and global MDR-associated Kp CGs (CFR 17.0%, n = 36; aHR 1.52, CI 0.98-2.38, p = 0.06) compared to Kp CGs not associated with MDR (CFR 10.1%, n = 46).

CONCLUSION

Bacterial traits, beyond antimicrobial resistance, have a major impact on the clinical outcome of KpSC BSIs. The global spread of MDR-associated Kp CGs is driven by other mechanisms than antibiotic selection alone. Further insights into virulence determinants, and their association with phylogenetic lineages are needed to better understand the epidemiology of KpSC infection and clinical outcome.

Evaluation of the accuracy of bacterial genome reconstruction with Oxford Nanopore R10.4.1 long-read-only sequencing

Whole-genome reconstruction of bacterial pathogens has become an important tool for tracking transmission and antimicrobial resistance gene spread, but highly accurate and complete assemblies have largely only historically been achievable using hybrid long- and short-read sequencing. We previously found the Oxford Nanopore Technologies (ONT) R10.4/kit12 flowcell/chemistry produced improved assemblies over the R9.4.1/kit10 combination, however long-read only assemblies contained more errors compared to Illumina-ONT hybrid assemblies. ONT have since released an R10.4.1/kit14 flowcell/chemistry upgrade and recommended the use of Bovine Serum Albumin (BSA) during library preparation, both of which reportedly increase accuracy and yield. They have also released updated basecallers trained using native bacterial DNA containing methylation sites intended to fix systematic basecalling errors, including common adenosine (A) to guanine (G) and cytosine (C) to thymine (T) substitutions. To evaluate these improvements, we successfully sequenced four bacterial reference strains, namely Escherichia coli, Klebsiella pneumoniae, Pseudomonas aeruginosa and Staphylococcus aureus, and nine genetically diverse E. coli bloodstream infection-associated isolates from different phylogroups and sequence types, both with and without BSA. These sequences were de novo assembled and compared against Illumina-corrected reference genomes. In this small evaluation of 13 isolates we found that nanopore long-read-only R10.4.1/kit 14 assemblies with updated basecallers trained using bacterial methylated DNA produce accurate assemblies with ≥40×depth, sufficient to be cost-effective compared with hybrid ONT/Illumina sequencing in our setting.

Whole-genome characterisation of Escherichia coli isolates from patients with bacteraemia presenting with sepsis or septic shock in Spain: a multicentre cross-sectional study

BACKGROUND

Escherichia coli is the most frequent cause of bloodstream infections (BSIs). About one-third of patients with BSIs due to E coli develop sepsis or shock. The objective of this study is to characterise the microbiological features of E coli blood isolates causing sepsis or septic shock to provide exploratory information for future diagnostic, preventive, or therapeutic interventions.

METHODS

E coli blood isolates from a multicentre cross-sectional study of patients older than 14 years presenting with sepsis or septic shock (according to the Third International Consensus Definitions for Sepsis and Septic Shock criteria) from hospitals in Spain between Oct 4, 2016, and Oct 15, 2017, were studied by whole-genome sequencing. Phylogroups, sequence types (STs), serotype, FimH types, antimicrobial resistance (AMR) genes, pathogenicity islands, and virulence factors were identified. Susceptibility testing was performed by broth microdilution. The main outcome of this study was the characterisation of the E coli blood isolates in terms of population structure by phylogroups, groups (group 1: phylogroups B2, F, and G; group 2: A, B1, and C; group 3: D), and STs and distribution by geographical location and bloodstream infection source. Other outcomes were virulence score and prevalence of virulence-associated genes, pathogenicity islands, AMR, and AMR-associated genes. Frequencies were compared using χ² or Fisher's exact tests, and continuous variables using the Mann-Whitney test, with Bonferroni correction for multiple comparisons.

FINDINGS

We analysed 224 isolates: 140 isolates (63%) were included in phylogenetic group 1, 52 (23%) in group 2, and 32 (14%) in group 3. 85 STs were identified, with four comprising 44% (n=98) of the isolates: ST131 (38 [17%]), ST73 (25 [11%]), ST69 (23 [10%]), and ST95 (12 [5%]). No significant differences in phylogroup or ST distribution were found according to geographical areas or source of bloodstream infection, except for ST95, which was more frequent in urinary tract infections than in other sources (11 [9%] of 116 vs 1 [1%] of 108, p=0·0045). Median virulence score was higher in group 1 (median 25·0 [IQR 20·5-29·0) than in group 2 (median 14·5 [9·0-20·0]; p<0·0001) and group 3 (median 21 [16·5-23·0]; p<0·0001); prevalence of several pathogenicity islands was higher in group 1. No significant differences were found between phylogenetic groups in proportions of resistance to antibiotics. ST73 had higher median virulence score (32 [IQR 29-35]) than the other predominant clones (median range 21-28). Some virulence genes and pathogenicity islands were significantly associated with each ST. ST131 isolates had higher prevalence of AMR and a higher proportion of AMR genes, notably blaCTX-M-15 and blaOXA-1.

INTERPRETATION

In this exploratory study, the population structure of E coli causing sepsis or shock was similar to previous studies that included all bacteraemic isolates. Virulence genes, pathogenicity islands, and AMR genes were not randomly distributed among phylogroups or STs. These results provide a comprehensive characterisation of invasive E coli isolates causing severe response syndrome. Future studies are required to determine the contribution of these microbiological factors to severe clinical presentation and worse outcomes in patients with E coli bloodstream infection.

FUNDING

Instituto de Salud Carlos III.

The plasmidome associated with Gram-negative bloodstream infections: A large-scale observational study using complete plasmid assemblies

Plasmids carry genes conferring antimicrobial resistance and other clinically important traits, and contribute to the rapid dissemination of such genes. Previous studies using complete plasmid assemblies, which are essential for reliable inference, have been small and/or limited to plasmids carrying antimicrobial resistance genes (ARGs). In this study, we sequenced 1,880 complete plasmids from 738 isolates from bloodstream infections in Oxfordshire, UK. The bacteria had been originally isolated in 2009 (194 isolates) and 2018 (368 isolates), plus a stratified selection from intervening years (176 isolates). We demonstrate that plasmids are largely, but not entirely, constrained to a single host species, although there is substantial overlap between species of plasmid gene-repertoire. Most ARGs are carried by a relatively small number of plasmid groups with biological features that are predictable. Plasmids carrying ARGs (including those encoding carbapenemases) share a putative 'backbone' of core genes with those carrying no such genes. These findings suggest that future surveillance should, in addition to tracking plasmids currently associated with clinically important genes, focus on identifying and monitoring the dissemination of high-risk plasmid groups with the potential to rapidly acquire and disseminate these genes.

Modulation of multidrug-resistant clone success in Escherichia coli populations: a longitudinal, multi-country, genomic and antibiotic usage cohort study

BACKGROUND

The effect of antibiotic usage on the success of multidrug-resistant (MDR) clones in a population remains unclear. With this genomics-based molecular epidemiology study, we aimed to investigate the contribution of antibiotic use to Escherichia coli clone success, relative to intra-strain competition for colonisation and infection.

METHODS

We sequenced all the available E coli bloodstream infection isolates provided by the British Society for Antimicrobial Chemotherapy (BSAC) from 2012 to 2017 (n=718) and combined these with published data from the UK (2001-11; n=1090) and Norway (2002-17; n=3254). Defined daily dose (DDD) data from the European Centre for Disease Prevention and Control (retrieved on Sept 21, 2021) for major antibiotic classes (β-lactam, tetracycline, macrolide, sulfonamide, quinolone, and non-penicillin β-lactam) were used together with sequence typing, resistance profiling, regression analysis, and non-neutral Wright-Fisher simulation-based modelling to enable systematic comparison of resistance levels, clone success, and antibiotic usage between the UK and Norway.

FINDINGS

Sequence type (ST)73, ST131, ST95, and ST69 accounted for 892 (49·3%) of 1808 isolates in the BSAC collection. In the UK, the proportion of ST69 increased between 2001-10 and 2011-17 (p=0·0004), whereas the proportions of ST73 and ST95 did not vary between periods. ST131 expanded quickly after its emergence in 2003 and its prevalence remained consistent throughout the study period (apart from a brief decrease in 2009-10). The extended-spectrum β-lactamase (ESBL)-carrying, globally disseminated MDR clone ST131-C2 showed overall greater success in the UK (154 [56·8%] of 271 isolates in 2003-17) compared with Norway (51 [18·3%] of 278 isolates in 2002-17; p<0·0001). DDD data indicated higher total use of antimicrobials in the UK, driven mainly by the class of non-penicillin β-lactams, which were used between 2·7-times and 5·1-times more in the UK per annum (ratio mean 3·7 [SD 0·8]). This difference was associated with the higher success of the MDR clone ST131-C2 (pseudo-R2 69·1%). A non-neutral Wright-Fisher model replicated the observed expansion of non-MDR and MDR sequence types under higher DDD regimes.

INTERPRETATION

Our study indicates that resistance profiles of contemporaneously successful clones can vary substantially, warranting caution in the interpretation of correlations between aggregate measures of resistance and antibiotic usage. Our study further suggests that in countries with low-to-moderate use of antibiotics, such as the UK and Norway, the extent of non-penicillin β-lactam use modulates rather than determines the success of widely disseminated MDR ESBL-carrying E coli clones. Detailed understanding of underlying causal drivers of success is important for improved control of resistant pathogens.

FUNDING

Trond Mohn Foundation, Marie Skłodowska-Curie Actions, European Research Council, Royal Society, and Wellcome Trust.

Clinical and Genomic Characterization of a Cohort of Patients With Klebsiella pneumoniae Bloodstream Infection

BACKGROUND

The clinical and microbial factors associated with Klebsiella pneumoniae bloodstream infections (BSIs) are not well characterized. Prior studies have focused on highly resistant or hypervirulent isolates, limiting our understanding of K. pneumoniae strains that commonly cause BSI. We performed a record review and whole-genome sequencing to investigate the clinical characteristics, bacterial diversity, determinants of antimicrobial resistance, and risk factors for in-hospital death in a cohort of patients with K. pneumoniae BSI.

METHODS

We identified 562 patients at Massachusetts General Hospital with K. pneumoniae BSIs between 2016 and 2022. We collected data on comorbid conditions, infection source, clinical outcomes, and antibiotic resistance and performed whole-genome sequencing on 108 sequential BSI isolates from 2021 to 2022.

RESULTS

Intra-abdominal infection was the most common source of infection accounting for 34% of all BSIs. A respiratory tract source accounted for 6% of BSIs but was associated with a higher in-hospital mortality rate (adjusted odds ratio, 5.4 [95% confidence interval, 2.2-12.8]; P < .001 for comparison with other sources). Resistance to the first antibiotic prescribed was also associated with a higher risk of death (adjusted odds ratio, 5.2 [95% confidence interval, 2.2-12.4]; P < .001). BSI isolates were genetically diverse, and no clusters of epidemiologically and genetically linked cases were observed. Virulence factors associated with invasiveness were observed at a low prevalence, although an unexpected association between O-antigen type and the source of infection was found.

CONCLUSIONS

These observations demonstrate the versatility of K. pneumoniae as an opportunistic pathogen and highlight the need for new approaches for surveillance and the rapid identification of patients with invasive antimicrobial-resistant K. pneumoniae infection.

Deep learning and single-cell phenotyping for rapid antimicrobial susceptibility detection in Escherichia coli

The rise of antimicrobial resistance (AMR) is one of the greatest public health challenges, already causing up to 1.2 million deaths annually and rising. Current culture-based turnaround times for bacterial identification in clinical samples and antimicrobial susceptibility testing (AST) are typically 18-24 h. We present a novel proof-of-concept methodological advance in susceptibility testing based on the deep-learning of single-cell specific morphological phenotypes directly associated with antimicrobial susceptibility in Escherichia coli. Our models can reliably (80% single-cell accuracy) classify untreated and treated susceptible cells for a lab-reference fully susceptible E. coli strain, across four antibiotics (ciprofloxacin, gentamicin, rifampicin and co-amoxiclav). For ciprofloxacin, we demonstrate our models reveal significant (p < 0.001) differences between bacterial cell populations affected and unaffected by antibiotic treatment, and show that given treatment with a fixed concentration of 10 mg/L over 30 min these phenotypic effects correlate with clinical susceptibility defined by established clinical breakpoints. Deploying our approach on cell populations from six E. coli strains obtained from human bloodstream infections with varying degrees of ciprofloxacin resistance and treated with a range of ciprofloxacin concentrations, we show single-cell phenotyping has the potential to provide equivalent information to growth-based AST assays, but in as little as 30 min.

Parallel loss of type VI secretion systems in two multi-drug-resistant Escherichia coli lineages

The repeated emergence of multi-drug-resistant (MDR) Escherichia coli clones is a threat to public health globally. In recent work, drug-resistant E. coli were shown to be capable of displacing commensal E. coli in the human gut. Given the rapid colonization observed in travel studies, it is possible that the presence of a type VI secretion system (T6SS) may be responsible for the rapid competitive advantage of drug-resistant E. coli clones. We employed large-scale genomic approaches to investigate this hypothesis. First, we searched for T6SS genes across a curated dataset of over 20 000 genomes representing the full phylogenetic diversity of E. coli. This revealed large, non-phylogenetic variation in the presence of T6SS genes. No association was found between T6SS gene carriage and MDR lineages. However, multiple clades containing MDR clones have lost essential structural T6SS genes. We characterized the T6SS loci of ST410 and ST131 and identified specific recombination and insertion events responsible for the parallel loss of essential T6SS genes in two MDR clones.

The clinical, genomic, and microbiological profile of invasive multi-drug resistant Escherichia coli in a major teaching hospital in the United Kingdom

Escherichia coli is a ubiquitous component of the human gut microbiome, but is also a common pathogen, causing around 40, 000 bloodstream infections (BSI) in the United Kingdom (UK) annually. The number of E. coli BSI has increased over the last decade in the UK, and emerging antimicrobial resistance (AMR) profiles threaten treatment options. Here, we combined clinical, epidemiological, and whole genome sequencing data with high content imaging to characterise over 300 E. coli isolates associated with BSI in a large teaching hospital in the East of England. Overall, only a limited number of sequence types (ST) were responsible for the majority of organisms causing invasive disease. The most abundant (20 % of all isolates) was ST131, of which around 90 % comprised the pandemic O25b:H4 group. ST131-O25b:H4 isolates were frequently multi-drug resistant (MDR), with a high prevalence of extended spectrum β-lactamases (ESBL) and fluoroquinolone resistance. There was no association between AMR phenotypes and the source of E. coli bacteraemia or whether the infection was healthcare-associated. Several clusters of ST131 were genetically similar, potentially suggesting a shared transmission network. However, there was no clear epidemiological associations between these cases, and they included organisms from both healthcare-associated and non-healthcare-associated origins. The majority of ST131 isolates exhibited strong binding with an anti-O25b antibody, raising the possibility of developing rapid diagnostics targeting this pathogen. In summary, our data suggest that a restricted set of MDR E. coli populations can be maintained and spread across both community and healthcare settings in this location, contributing disproportionately to invasive disease and AMR.

Temporal dynamics of genetically heterogeneous extended-spectrum cephalosporin-resistant Escherichia coli bloodstream infections

Extended-spectrum cephalosporin-resistant Escherichia coli (ESC-R-Ec) is an urgent public health threat with sequence type clonal complex 131 (STc131), phylogroup B2 strains being particularly concerning as the dominant cause of ESC-R-Ec infections. To address the paucity of recent ESC-R-Ec molecular epidemiology data in the United States, we used whole-genome sequencing (WGS) to fully characterize a large cohort of invasive ESC-R-Ec at a tertiary care cancer center in Houston, Texas, collected from 2016 to 2020. During the study time frame, there were 1,154 index E. coli bloodstream infections (BSIs) of which 389 (33.7%) were ESC-R-Ec. Using time series analyses, we identified a temporal dynamic of ESC-R-Ec distinct from ESC-susceptible E. coli (ESC-S-Ec), with cases peaking in the last 6 months of the calendar year. WGS of 297 ESC-R-Ec strains revealed that while STc131 strains accounted for ~45% of total BSIs, the proportion of STc131 strains remained stable across the study time frame with infection peaks driven by genetically heterogeneous ESC-R-Ec clonal complexes. bla CTX-M variants accounted for most β-lactamases conferring the ESC-R phenotype (89%; 220/248 index ESC-R-Ec), and amplification of bla CTX-M genes was widely detected in ESC-R-Ec strains, particularly in carbapenem non-susceptible, recurrent BSI strains. Bla CTX-M-55 was significantly enriched within phylogroup A strains, and we identified bla CTX-M-55 plasmid-to-chromosome transmission occurring across non-B2 strains. Our data provide important information regarding the current molecular epidemiology of invasive ESC-R-Ec infections at a large tertiary care cancer center and provide novel insights into the genetic basis of observed temporal variability for these clinically important pathogens. IMPORTANCE Given that E. coli is the leading cause of worldwide ESC-R Enterobacterales infections, we sought to assess the current molecular epidemiology of ESC-R-Ec using a WGS analysis of many BSIs over a 5-year period. We identified fluctuating temporal dynamics of ESC-R-Ec infections, which have also recently been identified in other geographical regions such as Israel. Our WGS data allowed us to visualize the stable nature of STc131 over the study period and demonstrate a limited but genetically diverse group of ESC-R-Ec clonal complexes are detected during infection peaks. Additionally, we provide a widespread assessment of β-lactamase gene copy number in ESC-R-Ec infections and delineate mechanisms by which such amplifications are achieved in a diverse array of ESC-R-Ec strains. These data suggest that serious ESC-R-Ec infections are driven by a diverse array of strains in our cohort and impacted by environmental factors suggesting that community-based monitoring could inform novel preventative measures.

Virulence Factor Genes in Invasive Escherichia coli Are Associated with Clinical Outcomes and Disease Severity in Patients with Sepsis: A Prospective Observational Cohort Study

BACKGROUND

Escherichia coli harbours virulence factors that facilitate the development of bloodstream infections. Studies determining virulence factors in clinical isolates often have limited access to clinical data and lack associations with patient outcome. The goal of this study was to correlate sepsis outcome and virulence factors of clinical E. coli isolates in a large cohort.

METHODS

Patients presenting at the emergency department whose blood cultures were positive for E. coli were prospectively included. Clinical and laboratory parameters were collected at admission. SOFA-score was calculated to determine disease severity. Patient outcomes were in-hospital mortality and ICU admission. Whole genome sequencing was performed for E. coli isolates and virulence genes were detected using the VirulenceFinder database.

RESULTS

In total, 103 E. coli blood isolates were sequenced. Isolates had six to 41 virulence genes present. One virulence gene, kpsMII_K23, a K1 capsule group 2 of E. coli type K23, was significantly more present in isolates of patients who died. kpsMII_K23 and cvaC (Microcin C) were significantly more frequent in isolates of patients who were admitted to the ICU. Fourteen virulence genes (mchB, mchC, papA_fsiA_F16, sat, senB, iucC, iutA, iha, sfaD, cnf1, focG, vat, cldB, and mcmA) significantly differed between patients with and without sepsis.

CONCLUSIONS

Microcins, toxins, and fimbriae were associated with disease severity. Adhesins and iron uptake proteins seemed to be protective. Two genes were associated with worse clinical outcome. These findings contribute to a better understanding of host-pathogen interactions and could help identifying patients most at risk for a worse outcome.

Evolutionary and functional history of the Escherichia coli K1 capsule

Escherichia coli is a leading cause of invasive bacterial infections in humans. Capsule polysaccharide has an important role in bacterial pathogenesis, and the K1 capsule has been firmly established as one of the most potent capsule types in E. coli through its association with severe infections. However, little is known about its distribution, evolution and functions across the E. coli phylogeny, which is fundamental to elucidating its role in the expansion of successful lineages. Using systematic surveys of invasive E. coli isolates, we show that the K1-cps locus is present in a quarter of bloodstream infection isolates and has emerged in at least four different extraintestinal pathogenic E. coli (ExPEC) phylogroups independently in the last 500 years. Phenotypic assessment demonstrates that K1 capsule synthesis enhances E. coli survival in human serum independent of genetic background, and that therapeutic targeting of the K1 capsule re-sensitizes E. coli from distinct genetic backgrounds to human serum. Our study highlights that assessing the evolutionary and functional properties of bacterial virulence factors at population levels is important to better monitor and predict the emergence of virulent clones, and to also inform therapies and preventive medicine to effectively control bacterial infections whilst significantly lowering antibiotic usage.

Integrative omics identifies conserved and pathogen-specific responses of sepsis-causing bacteria

Even in the setting of optimal resuscitation in high-income countries severe sepsis and septic shock have a mortality of 20-40%, with antibiotic resistance dramatically increasing this mortality risk. To develop a reference dataset enabling the identification of common bacterial targets for therapeutic intervention, we applied a standardized genomic, transcriptomic, proteomic and metabolomic technological framework to multiple clinical isolates of four sepsis-causing pathogens: Escherichia coli, Klebsiella pneumoniae species complex, Staphylococcus aureus and Streptococcus pyogenes. Exposure to human serum generated a sepsis molecular signature containing global increases in fatty acid and lipid biosynthesis and metabolism, consistent with cell envelope remodelling and nutrient adaptation for osmoprotection. In addition, acquisition of cholesterol was identified across the bacterial species. This detailed reference dataset has been established as an open resource to support discovery and translational research.

Genomic analysis of the international high-risk clonal lineage Klebsiella pneumoniae sequence type 395

BACKGROUND

Klebsiella pneumoniae, which is frequently associated with hospital- and community-acquired infections, contains multidrug-resistant (MDR), hypervirulent (hv), non-MDR/non-hv as well as convergent representatives. It is known that mostly international high-risk clonal lineages including sequence types (ST) 11, 147, 258, and 307 drive their global spread. ST395, which was first reported in the context of a carbapenemase-associated outbreak in France in 2010, is a less well-characterized, yet emerging clonal lineage.

METHODS

We computationally analyzed a large collection of K. pneumoniae ST395 genomes (n = 297) both sequenced in this study and reported previously. By applying multiple bioinformatics tools, we investigated the core-genome phylogeny and evolution of ST395 as well as distribution of accessory genome elements associated with antibiotic resistance and virulence features.

RESULTS

Clustering of the core-SNP alignment revealed four major clades with eight smaller subclades. The subclades likely evolved through large chromosomal recombination, which involved different K. pneumoniae donors and affected, inter alia, capsule and lipopolysaccharide antigen biosynthesis regions. Most genomes contained acquired resistance genes to extended-spectrum cephalosporins, carbapenems, and other antibiotic classes carried by multiple plasmid types, and many were positive for hypervirulence markers, including the siderophore aerobactin. The detection of "hybrid" resistance and virulence plasmids suggests the occurrence of the convergent ST395 pathotype.

CONCLUSIONS

To the best of our knowledge, this is the first study that investigated a large international collection of K. pneumoniae ST395 genomes and elucidated phylogenetics and detailed genomic characteristics of this emerging high-risk clonal lineage.

Evolution of Virulence, Fitness, and Carbapenem Resistance Transmission in ST23 Hypervirulent Klebsiella pneumoniae with the Capsular Polysaccharide Synthesis Gene wcaJ Inserted via Insertion Sequence Elements

Carbapenem-resistant hypervirulent Klebsiella pneumoniae (CR-hvKP) is recognized as a threat worldwide, but the mechanisms underlying its emergence remain unclear. As most CR-hvKP isolates are not hypermucoviscous, we speculated that the evolution of the capsule might result in the convergence of carbapenem resistance and hypervirulence. Here, 2,096 K. pneumoniae isolates were retrospectively collected to screen the ST23-K1 clone, and hypervirulence was roughly defined as being highly resistant to serum killing. The effect of wcaJ on the capsule, virulence, fitness, and resistance acquisition was further analyzed. The capsule gene wcaJ, inserted by ISKpn26/ISKpn74, was identified via whole-genome sequencing in four hvKP, but not hypermucoviscous, isolates. Uronic acid quantitation results revealed that these isolates produced significantly less capsular polysaccharides than NTUH-K2044. A significant increase in capsular production was observed in wcaJ-complemented isolates and confirmed by transmission electron microscopy. Further, all wcaJ-complemented isolates acquired greater resistance to macrophage phagocytosis, and one representative isolate resulted in a significantly higher mortality rate than the parental isolate in mice, indicating that wcaJ inactivation might compromise virulence. However, isolates with wcaJ interruption demonstrated a lower fitness cost and a high conjugation frequency of the blaKPC-2 plasmid, raising concerns about the emergence of carbapenem resistance in hvKP. IMPORTANCE Klebsiella pneumoniae is one of the most common nosocomial pathogens worldwide, and we speculated that the evolution of the capsule might result in the convergence of carbapenem resistance and hypervirulence of K. pneumoniae. The wcaJ gene was first reported to be interrupted by insertion sequence elements in ST23-K1 hypervirulent Klebsiella pneumoniae, resulting in little capsule synthesis, which plays an important role in virulence. We examined the effect of wcaJ on the capsule, virulence, and fitness. Isolates with wcaJ interruption might compromise virulence and demonstrated a lower fitness cost and a high conjugation frequency of the blaKPC-2 plasmid, highlighting its role as a potential factor facilitating hypervirulence and carbapenem resistance.

Benchmarking taxonomic classifiers with Illumina and Nanopore sequence data for clinical metagenomic diagnostic applications

Culture-independent metagenomic detection of microbial species has the potential to provide rapid and precise real-time diagnostic results. However, it is potentially limited by sequencing and taxonomic classification errors. We use simulated and real-world data to benchmark rates of species misclassification using 100 reference genomes for each of the ten common bloodstream pathogens and six frequent blood-culture contaminants (n=1568, only 68 genomes were available for Micrococcus luteus). Simulating both with and without sequencing error for both the Illumina and Oxford Nanopore platforms, we evaluated commonly used classification tools including Kraken2, Bracken and Centrifuge, utilizing mini (8 GB) and standard (30–50 GB) databases. Bracken with the standard database performed best, the median percentage of reads across both sequencing platforms identified correctly to the species level was 97.8% (IQR 92.7:99.0) [range 5:100]. For Kraken2 with a mini database, a commonly used combination, median species-level identification was 86.4% (IQR 50.5:93.7) [range 4.3:100]. Classification performance varied by species, with Escherichia coli being more challenging to classify correctly (probability of reads being assigned to the correct species: 56.1–96.0%, varying by tool used). Human read misclassification was negligible. By filtering out shorter Nanopore reads we found performance similar or superior to Illumina sequencing, despite higher sequencing error rates. Misclassification was more common when the misclassified species had a higher average nucleotide identity to the true species. Our findings highlight taxonomic misclassification of sequencing data occurs and varies by sequencing and analysis workflow. To account for ‘bioinformatic contamination’ we present a contamination catalogue that can be used in metagenomic pipelines to ensure accurate results that can support clinical decision making.

Prevalence of Escherichia coli ST1193 Causing Intracranial Infection in Changsha, China

ST1193 is an emerging new virulent and resistant clone among Escherichia coli with a tendency to spread rapidly across the globe. However, the prevalence of intracranial infection-causing E. coli ST1193 is rarely reported. This study aimed at determining the prevalence of E. coli ST1193 isolates, causing intracranial infections in Changsha, central China. A total of 28 E. coli isolates were collected from the cerebrospinal fluid of patients with intracranial infection over a four-year period. All isolates were differentiated using multilocus sequence typing (MLST), and phylogenetic grouping, and tested for antibiotic resistance. MLST analysis showed 11 sequence types (ST) among the 28 E. coli isolates. The most prevalent ST was B2-ST1193 (28.6%, 8/28), followed by B2-ST131 (21.4%, 6/28) and F-ST648 (10.7%, 3/28). Of the eight ST1193 isolates, three carried CTX-M-55, and one carried CTX-M-27. All eight ST1193 isolates were resistant to Ciprofloxacin, showing gyrA1AB/parC4A mutations. Two ST1193 isolates carried the aac(6′)-Ib-cr gene. All ST1193 isolates were recovered from infants with meningitis, with a fatal outcome for one three-month-old infant. ST1193 has emerged as the predominant type of E. coli strain causing intracranial infections in Changsha, China. This study highlights the importance of implementing appropriate surveillance measures to prevent the spread of this emerging public health threat.

Nosocomial dissemination of hypervirulent Klebsiella pneumoniae with high-risk clones among children in Shanghai

Objectives

Although hypervirulent Klebsiella pneumoniae (hvKp) is an increasing public health problem, there remains limited epidemiological information regarding hvKp infections in children. Here, we conducted a clinical, molecular and phenotypic surveillance of hvKp strains in a pediatric population.

Methods

Non-repetitive K. pneumoniae (Kp) strains consecutively collected during 2019-2020 were screened for hypervirulence genes (prmpA, prmpA2, iucA, iroB, and peg344) using PCR. Positive strains were further characterized by four phenotypic assays (string test, serum killing assay, siderophore production, Galleria mellonella lethality assay), followed by murine sepsis model to determine virulence in vitro and in vivo. Also, capsular types, sequence types, plasmid replicon types, antimicrobial resistance determinants and susceptibility were analyzed.

Results

A total of 352 isolates were collected, wherein 83 (23.6%) were hypervirulence genes-positive Kp (hgKp). A significant increase in KPC-2-producing KL47-ST11 among hgKp strains was observed, from 5.3% (1/19) in 2019 to 67.6% (25/37) in 2020 (P<.0001), suggesting the potential dissemination of the hybrid virulence and carbapenem-resistance encoding plasmid among children. Further, hgKp isolates were classified into hvKp (n = 27) and hgKp-low virulence (hgKp-Lv) (n = 56) based on virulence phenotypic assays. In hvKp, diverse genetic clones were observed and K1-ST23 or K2-ST25 strains with sensitivity to multiple antibiotics were prevalent (25.9%, 7/27). Compared with hgKp-Lv, hvKp infection had a higher propensity to involve severe pneumonia (22.2% vs. 12.5%) in elder children and significant higher mortality in mice (P = 0.0086). Additionally, either hvKp or hgKp-Lv infections were mostly healthcare-associated and hospital-acquired (74.1% vs. 91.9%).

Conclusions

These data suggest that K1-ST23 and K2-ST25 are high-risk clones of hvKp, and the genetic convergence of virulence and carbapenem-resistance is increasing among children. Control measures are needed to prevent the dissemination in clinical settings.

Klebsiella pneumoniae hijacks the Toll-IL-1R protein SARM1 in a type I IFN-dependent manner to antagonize host immunity

Many bacterial pathogens antagonize host defense responses by translocating effector proteins into cells. It remains an open question how those pathogens not encoding effectors counteract anti-bacterial immunity. Here, we show that Klebsiella pneumoniae exploits the evolutionary conserved innate protein SARM1 to regulate negatively MyD88- and TRIF-governed inflammation, and the activation of the MAP kinases ERK and JNK. SARM1 is required for Klebsiella induction of interleukin-10 (IL-10) by fine-tuning the p38-type I interferon (IFN) axis. SARM1 inhibits the activation of Klebsiella-induced absent in melanoma 2 inflammasome to limit IL-1β production, suppressing further inflammation. Klebsiella exploits type I IFNs to induce SARM1 in a capsule and lipopolysaccharide O-polysaccharide-dependent manner via the TLR4-TRAM-TRIF-IRF3-IFNAR1 pathway. Absence of SARM1 reduces the intracellular survival of K. pneumoniae in macrophages, whereas sarm1-deficient mice control the infection. Altogether, our results illustrate an anti-immunology strategy deployed by a human pathogen. SARM1 inhibition will show a beneficial effect to treat Klebsiella infections.

•

SARM1 is an evolutionary conserved innate immune protein with a TIR domain

•

Klebsiella pneumoniae induces SARM1 to limit inflammation, and to induce IL-10

•

SARM1 inhibits the action of AIM2 inflammasome to limit IL-1β production

•

In vivo, absence of SARM1 facilitates the clearance of Klebsiella pneumoniae

Genomic and Temporal Trends in Canine ExPEC Reflect Those of Human ExPEC

Companion animals and humans are known to share extraintestinal pathogenic Escherichia coli (ExPEC), but the extent of E. coli sequence types (STs) that cause extraintestinal diseases in dogs is not well understood. Here, we generated whole-genome sequences of 377 ExPEC collected by the University of Melbourne Veterinary Hospital from dogs over an 11-year period from 2007 to 2017. Isolates were predominantly from urogenital tract infections (219, 58.1%), but isolates from gastrointestinal specimens (51, 13.5%), general infections (72, 19.1%), and soft tissue infections (34, 9%) were also represented. A diverse collection of 53 STs were identified, with 18 of these including at least five sequences. The five most prevalent STs were ST372 (69, 18.3%), ST73 (31, 8.2%), ST127 (22, 5.8%), ST80 (19, 5.0%), and ST58 (14, 3.7%). Apart from ST372, all of these are prominent human ExPEC STs. Other common ExPEC STs identified included ST12, ST131, ST95, ST141, ST963, ST1193, ST88, and ST38. Virulence gene profiles, antimicrobial resistance carriage, and trends in plasmid carriage for specific STs were generally reflective of those seen in humans. Many of the prominent STs were observed repetitively over an 11-year time span, indicating their persistence in the dogs in the community, which is most likely driven by household sharing of E. coli between humans and their pets. The case of ST372 as a dominant canine lineage observed sporadically in humans is flagged for further investigation. IMPORTANCE Pathogenic E. coli that causes extraintestinal infections (ExPEC) in humans and canines represents a significant burden in hospital and veterinary settings. Despite the obvious interrelationship between dogs and humans favoring both zoonotic and anthropozoonotic infections, whole-genome sequencing projects examining large numbers of canine-origin ExPEC are lacking. In support of anthropozoonosis, we found that most STs from canine infections are dominant human ExPEC STs (e.g., ST73, ST127, ST131) with similar genomic traits, such as plasmid carriage and virulence gene burden. In contrast, we identified ST372 as the dominant canine ST and a sporadic cause of infection in humans, supporting zoonotic transfer. Furthermore, we highlight that, as is the case in humans, STs in canine disease are consistent over time, implicating the gastrointestinal tract as the major community reservoir, which is likely augmented by exposure to human E. coli via shared diet and proximity.

Genomic dissection of Klebsiella pneumoniae infections in hospital patients reveals insights into an opportunistic pathogen

Klebsiella pneumoniae is a major cause of opportunistic healthcare-associated infections, which are increasingly complicated by the presence of extended-spectrum beta-lactamases (ESBLs) and carbapenem resistance. We conducted a year-long prospective surveillance study of K. pneumoniae clinical isolates in hospital patients. Whole-genome sequence (WGS) data reveals a diverse pathogen population, including other species within the K. pneumoniae species complex (18%). Several infections were caused by K. variicola/K. pneumoniae hybrids, one of which shows evidence of nosocomial transmission. A wide range of antimicrobial resistance (AMR) phenotypes are observed, and diverse genetic mechanisms identified (mainly plasmid-borne genes). ESBLs are correlated with presence of other acquired AMR genes (median n = 10). Bacterial genomic features associated with nosocomial onset are ESBLs (OR 2.34, p = 0.015) and rhamnose-positive capsules (OR 3.12, p < 0.001). Virulence plasmid-encoded features (aerobactin, hypermucoidy) are observed at low-prevalence (<3%), mostly in community-onset cases. WGS-confirmed nosocomial transmission is implicated in just 10% of cases, but strongly associated with ESBLs (OR 21, p < 1 × 10⁻¹¹). We estimate 28% risk of onward nosocomial transmission for ESBL-positive strains vs 1.7% for ESBL-negative strains. These data indicate that K. pneumoniae infections in hospitalised patients are due largely to opportunistic infections with diverse strains, with an additional burden from nosocomially-transmitted AMR strains and community-acquired hypervirulent strains.

Klebsiella pneumoniae is an opportunistic pathogen of increasing public health concern due to the prevalence of antimicrobial resistance. Here, the authors provide insight into the resistance profiles, bacterial genome features and virulence genes, in a year-long prospective study of K. pneumoniae clinical isolates.

SARS-CoV-2 genomic surveillance as an evidence-based infection control approach in an offshore petroleum employee population

BACKGROUND

Industrial hygienists (IH) in the oil and gas business instituted an extraordinary number of safety protocols to limit spread of SARS-CoV-2 onto offshore platforms in the Gulf of Mexico. We used genomic surveillance to provide actionable information concerning the efficacy of their efforts.

METHODS

Over 6 months, employees at a single company were serology and PCR tested during a 1-5 day predeployment quarantine and when postdeployment symptoms were reported. From each positive test (n = 49), SARS-CoV-2 genomes were sequenced. Phylogenetic analysis was used to investigate the epidemiology of transmissions.

RESULTS

Genomic surveillance confirmed 2 viral strains were infecting 18 offshore workers. Genomic data combined with epidemiological data suggested that a change in quarantine protocols contributed to these outbreaks. A pre-deployment outbreak involved a WHO variant of interest (Theta) that had infected 4 international workers. Two additional predeployment clusters of infections were identified.

CONCLUSIONS

Our findings support that IH quarantine/testing protocols limited viral transmissions, halted offshore outbreaks, and stopped the spread of a variant of interest. The study demonstrates how genomic data can be used to understand viral transmission dynamics in employee populations and evaluate safety protocols in the offshore oil and gas industry.

Molecular epidemiology and antimicrobial resistance phenotype of paediatric bloodstream infections caused by Gram-negative bacteria

Background

Gram-negative organisms are common causes of bloodstream infection (BSI) during the neonatal period and early childhood. Whilst several large studies have characterised these isolates in adults, equivalent data (particularly incorporating whole genome sequencing) is lacking in the paediatric population.

Methods

We perform an epidemiological and sequencing based analysis of Gram-negative bloodstream infections (327 isolates (296 successfully sequenced) from 287 patients) in children <18 years old between 2008 and 2018 in Oxfordshire, UK.

Results

Here we show that the burden of infection lies predominantly in neonates and that most infections are caused by Escherichia coli, Klebsiella spp. and Enterobacter hormaechei. There is no evidence in our setting that the proportion of antimicrobial resistant isolates is increasing in the paediatric population although we identify some evidence of sub-breakpoint increases in gentamicin resistance. The population structure of E. coli BSI isolates in neonates and children mirrors that in adults with a predominance of STs 131/95/73/69 and the same proportions of O-antigen serotypes. In most cases in our setting there is no evidence of transmission/point-source acquisition and we demonstrate the utility of whole genome sequencing to refute a previously suspected outbreak.

Conclusions

Our findings support continued use of current empirical treatment guidelines and suggest that O-antigen targeted vaccines may have a role in reducing the incidence of neonatal sepsis.