Identification of Klebsiella capsule synthesis loci from whole genome data.

Klebsiella pneumoniae is a growing cause of healthcare-associated infections for which multi-drug resistance is a concern. Its polysaccharide capsule is a major virulence determinant and epidemiological marker. However, little is known about capsule epidemiology since serological typing is not widely accessible and many isolates are serologically non-typeable. Molecular typing techniques provide useful insights, but existing methods fail to take full advantage of the information in whole genome sequences. We investigated the diversity of the capsule synthesis loci (K-loci) among 2503 K. pneumoniae genomes. We incorporated analyses of full-length K-locus nucleotide sequences and also clustered protein-encoding sequences to identify, annotate and compare K-locus structures. We propose a standardized nomenclature for K-loci and present a curated reference database. A total of 134 distinct K-loci were identified, including 31 novel types. Comparative analyses indicated 508 unique protein-encoding gene clusters that appear to reassort via homologous recombination. Extensive intra- and inter-locus nucleotide diversity was detected among the wzi and wzc genes, indicating that current molecular typing schemes based on these genes are inadequate. As a solution, we introduce Kaptive, a novel software tool that automates the process of identifying K-loci based on full locus information extracted from whole genome sequences (https://github.com/katholt/Kaptive). This work highlights the extensive diversity of Klebsiella K-loci and the proteins that they encode. The nomenclature, reference database and novel typing method presented here will become essential resources for genomic surveillance and epidemiological investigations of this pathogen

Unicycler: Resolving bacterial genome assemblies from short and long sequencing reads.

The Illumina DNA sequencing platform generates accurate but short reads, which can be used to produce accurate but fragmented genome assemblies. Pacific Biosciences and Oxford Nanopore Technologies DNA sequencing platforms generate long reads that can produce complete genome assemblies, but the sequencing is more expensive and error-prone. There is significant interest in combining data from these complementary sequencing technologies to generate more accurate "hybrid" assemblies. However, few tools exist that truly leverage the benefits of both types of data, namely the accuracy of short reads and the structural resolving power of long reads. Here we present Unicycler, a new tool for assembling bacterial genomes from a combination of short and long reads, which produces assemblies that are accurate, complete and cost-effective. Unicycler builds an initial assembly graph from short reads using the de novo assembler SPAdes and then simplifies the graph using information from short and long reads. Unicycler uses a novel semi-global aligner to align long reads to the assembly graph. Tests on both synthetic and real reads show Unicycler can assemble larger contigs with fewer misassemblies than other hybrid assemblers, even when long-read depth and accuracy are low. Unicycler is open source (GPLv3) and available at github.com/rrwick/Unicycler

Completing bacterial genome assemblies with multiplex MinION sequencing

AbstractIllumina sequencing platforms have enabled widespread bacterial whole genome sequencing. While Illumina data is appropriate for many analyses, its short read length limits its ability to resolve genomic structure. This has major implications for tracking the spread of mobile genetic elements, including those which carry antimicrobial resistance determinants. Fully resolving a bacterial genome requires long-read sequencing such as those generated by Oxford Nanopore Technologies (ONT) platforms. Here we describe our use of the ONT MinION to sequence 12 isolates of Klebsiella pneumoniae on a single flow cell. We assembled each genome using a combination of ONT reads and previously available Illumina reads, and little to no manual intervention was needed to achieve fully resolved assemblies using the Unicycler hybrid assembler. Assembling only ONT reads with Canu was less effective, resulting in fewer resolved genomes and higher error rates even following error correction with Nanopolish. We demonstrate that multiplexed ONT sequencing is a valuable tool for high-throughput bacterial genome finishing. Specifically, we advocate the use of Illumina sequencing as a first analysis step, followed by ONT reads as needed to resolve genomic structure.Data summarySequence read files for all 12 isolates have been deposited in SRA, accessible through these NCBI BioSample accession numbers: SAMEA3357010, SAMEA3357043, SAMN07211279, SAMN07211280, SAMEA3357223, SAMEA3357193, SAMEA3357346, SAMEA3357374, SAMEA3357320, SAMN07211281, SAMN07211282, SAMEA3357405.A full list of SRA run accession numbers (both Illumina reads and ONT reads) for these samples are available in Table S1.Assemblies and sequencing reads corresponding to each stage of processing and analysis are provided in the following figshare project: https://figshare.com/projects/Completing_bacterial_genome_assemblies_with_multiplex_MinION_sequencing/23068Source code is provided in the following public GitHub repositories: https://github.com/rrwick/Bacterial-genome-assemblies-with-multiplex-MinION-sequencinghttps://github.com/rrwick/Porechophttps://github.com/rrwick/Fast5-to-FastqImpact StatementLike many research and public health laboratories, we frequently perform large-scale bacterial comparative genomics studies using Illumina sequencing, which assays gene content and provides the high-confidence variant calls needed for phylogenomics and transmission studies. However, problems often arise with resolving genome assemblies, particularly around regions that matter most to our research, such as mobile genetic elements encoding antibiotic resistance or virulence genes. These complexities can often be resolved by long sequence reads generated with PacBio or Oxford Nanopore Technologies (ONT) platforms. While effective, this has proven difficult to scale, due to the relatively high costs of generating long reads and the manual intervention required for assembly. Here we demonstrate the use of barcoded ONT libraries sequenced in multiplex on a single ONT MinION flow cell, coupled with hybrid assembly using Unicycler, to resolve 12 large bacterial genomes. Minor manual intervention was required to fully resolve small plasmids in five isolates, which we found to be underrepresented in ONT data. Cost per sample for the ONT sequencing was equivalent to Illumina sequencing, and there is potential for significant savings by multiplexing more samples on the ONT run. This approach paves the way for high-throughput and cost-effective generation of completely resolved bacterial genomes to become widely accessible.</jats:sec

Genetic diversity, mobilisation and spread of the yersiniabactin-encoding mobile element ICEKp in Klebsiella pneumoniae populations.

Mobile genetic elements (MGEs) that frequently transfer within and between bacterial species play a critical role in bacterial evolution, and often carry key accessory genes that associate with a bacteria's ability to cause disease. MGEs carrying antimicrobial resistance (AMR) and/or virulence determinants are common in the opportunistic pathogen Klebsiella pneumoniae, which is a leading cause of highly drug-resistant infections in hospitals. Well-characterised virulence determinants in K. pneumoniae include the polyketide synthesis loci ybt and clb (also known as pks), encoding the iron-scavenging siderophore yersiniabactin and genotoxin colibactin, respectively. These loci are located within an MGE called ICEKp, which is the most common virulence-associated MGE of K. pneumoniae, providing a mechanism for these virulence factors to spread within the population. Here we apply population genomics to investigate the prevalence, evolution and mobility of ybt and clb in K. pneumoniae populations through comparative analysis of 2498 whole-genome sequences. The ybt locus was detected in 40 % of K. pneumoniae genomes, particularly amongst those associated with invasive infections. We identified 17 distinct ybt lineages and 3 clb lineages, each associated with one of 14 different structural variants of ICEKp. Comparison with the wider population of the family Enterobacteriaceae revealed occasional ICEKp acquisition by other members. The clb locus was present in 14 % of all K. pneumoniae and 38.4 % of ybt+ genomes. Hundreds of independent ICEKp integration events were detected affecting hundreds of phylogenetically distinct K. pneumoniae lineages, including at least 19 in the globally-disseminated carbapenem-resistant clone CG258. A novel plasmid-encoded form of ybt was also identified, representing a new mechanism for ybt dispersal in K. pneumoniae populations. These data indicate that MGEs carrying ybt and clb circulate freely in the K. pneumoniae population, including among multidrug-resistant strains, and should be considered a target for genomic surveillance along with AMR determinants

Gastrointestinal Carriage Is a Major Reservoir of Klebsiella pneumoniae Infection in Intensive Care Patients.

BACKGROUND: Klebsiella pneumoniae is an opportunistic pathogen and leading cause of hospital-associated infections. Intensive care unit (ICU) patients are particularly at risk. Klebsiella pneumoniae is part of the healthy human microbiome, providing a potential reservoir for infection. However, the frequency of gut colonization and its contribution to infections are not well characterized. METHODS: We conducted a 1-year prospective cohort study in which 498 ICU patients were screened for rectal and throat carriage of K. pneumoniae shortly after admission. Klebsiella pneumoniae isolated from screening swabs and clinical diagnostic samples were characterized using whole genome sequencing and combined with epidemiological data to identify likely transmission events. RESULTS: Klebsiella pneumoniae carriage frequencies were estimated at 6% (95% confidence interval [CI], 3%-8%) among ICU patients admitted direct from the community, and 19% (95% CI, 14%-51%) among those with recent healthcare contact. Gut colonization on admission was significantly associated with subsequent infection (infection risk 16% vs 3%, odds ratio [OR] = 6.9, P < .001), and genome data indicated matching carriage and infection isolates in 80% of isolate pairs. Five likely transmission chains were identified, responsible for 12% of K. pneumoniae infections in ICU. In sum, 49% of K. pneumoniae infections were caused by the patients' own unique strain, and 48% of screened patients with infections were positive for prior colonization. CONCLUSIONS: These data confirm K. pneumoniae colonization is a significant risk factor for infection in ICU, and indicate ~50% of K. pneumoniae infections result from patients' own microbiota. Screening for colonization on admission could limit risk of infection in the colonized patient and others

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-11). 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 induces host metabolic stress that promotes tolerance to pulmonary infection

K. pneumoniae sequence type 258 (Kp ST258) is a major cause of healthcare-associated pneumonia. However, it remains unclear how it causes protracted courses of infection in spite of its expression of immunostimulatory lipopolysaccharide, which should activate a brisk inflammatory response and bacterial clearance. We predicted that the metabolic stress induced by the bacteria in the host cells shapes an immune response that tolerates infection. We combined in situ metabolic imaging and transcriptional analyses to demonstrate that Kp ST258 activates host glutaminolysis and fatty acid oxidation. This response creates an oxidant-rich microenvironment conducive to the accumulation of anti-inflammatory myeloid cells. In this setting, metabolically active Kp ST258 elicits a disease-tolerant immune response. The bacteria, in turn, adapt to airway oxidants by upregulating the type VI secretion system, which is highly conserved across ST258 strains worldwide. Thus, much of the global success of Kp ST258 in hospital settings can be explained by the metabolic activity provoked in the host that promotes disease tolerance. Keywords: immunometabolism, Klebsiella pneumoniae, immunosuppressive, anti-inflammatory, itaconate, Type Six Secretion Syste

Simulated hybrid assemblies: Read length and accuracy.

<p>NGA50 values segregated by read length and read accuracy. These plots summarise results across all reference genomes and replicate tests, but only include the tests of 8x long-read depth. For read lengths, the <i>p</i>-value is from a two-tailed <i>t</i>-test. For read accuracies, the <i>p</i>-value is from a one-way ANOVA test.</p

Simulated hybrid assemblies: Read length and accuracy.

<p>NGA50 values segregated by read length and read accuracy. These plots summarise results across all reference genomes and replicate tests, but only include the tests of 8x long-read depth. For read lengths, the <i>p</i>-value is from a two-tailed <i>t</i>-test. For read accuracies, the <i>p</i>-value is from a one-way ANOVA test.</p