PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of microgenLink to Publisher's site
 
Microb Genom. 2016 December; 2(12): e000102.
Published online 2016 December 12. doi:  10.1099/mgen.0.000102
PMCID: PMC5359410

Identification of Klebsiella capsule synthesis loci from whole genome data

Abstract

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.

Keywords: Klebsiella capsule K-locus genomic surveillance

Data Summary

1. Genome data generated and/or analysed in this work are available in the European Nucleotide Archive and/or PATRIC genome database, individual accession numbers are listed in Table S1.

2. Novel K-locus nucleotide sequences have been deposited in GenBank: accession numbers LT603702–LT603735 (https://www.ncbi.nlm.nih.gov/).

3. The Kaptive source code, along with the curated K-locus nucleotide and annotation databases have been deposited in GitHub – doi: 10.5281/zenodo.55773 (https://github.com/katholt/Kaptive).

Impact Statement

Klebsiella pneumoniae is a major cause of healthcare-associated infections and an urgent public-health threat for which robust epidemiological surveillance is paramount. These bacteria produce polysaccharide capsules that are important virulence factors, as well as epidemiological markers. Seventy-seven distinct capsule types (K-types) were defined by phenotypic studies in the 1950s–1970s, but the true extent of capsule diversity remains unknown. The increasing availability of whole genome sequences provides an unprecedented opportunity to explore capsule diversity, and here we report our study of the capsule synthesis loci (K-loci) in >2500 Klebsiella genomes. We identify a total of 134 distinct K-loci and show that they are extremely diverse, suggesting they encode at least 134 distinct K-types and are subject to unknown, diversifying selective pressures. Furthermore, we present a curated reference database and a new tool for the identification of K-loci from genome sequences, which will greatly assist epidemiological surveillance for K. pneumoniae, and other bacterial pathogens for which capsule epidemiology has been shown to be important.

Introduction

Klebsiella pneumoniae and its close relatives, Klebsiella variicola and Klebsiella quasipneumoniae, are opportunistic pathogens recognized as a significant threat to global health. Antimicrobial resistance, particularly multi-drug resistance and resistance to the carbapenems, is a major concern. Notably, there are a number of globally distributed, multi-drug resistant clones that cause outbreaks of healthcare-associated infections (Munoz-Price et al., 2013; The et al., 2015).

In order to control the emerging threat of K. pneumoniae sensu stricto, K. variicola and K. quasipneumoniae (hereafter collectively called Kp), there is an urgent requirement for genome-based surveillance. Recent advances in understanding population structure (Bialek-Davenet et al., 2014; Holt et al., 2015) highlight immense genomic diversity and provide a framework for tracking this pathogen. Useful strategies involve analyses of lineages or multi-locus sequence types in combination with resistance and virulence gene characterization (Bialek-Davenet et al., 2014), or phylogenetic analysis for outbreak investigation (Snitkin et al., 2012; The et al., 2015). However, reliable methods for tracking Kp capsular variation are currently lacking.

The polysaccharide capsule is the outermost layer of the Kp cell, protecting the bacterium from desiccation, phage and protist predation (March et al., 2013; Whitfield, 2006). The capsule is also a key virulence determinant. In contrast to capsulated strains, isogenic non-capsulated strains are unable to cause disease in murine infection models (Cortés et al., 2002; Lawlor et al., 2005). In addition, the capsule has been shown to suppress the host inflammatory response (Yoshida et al., 2000), provide resistance to antimicrobial immune-peptides (Campos et al., 2004), complement-mediated killing (Merino et al., 1992) and phagocytosis (Evrard et al., 2010; Lee et al., 2014; March et al., 2013). There are 77 immunologically distinct Klebsiella capsule types (K-types) defined by serology (Edmunds, 1954; Edwards & Fife, 1952; Ørskov & Fife-Asbury, 1977). However, serological typing requires specialist techniques and reagents not available to most microbiology laboratories, so is very rarely applied. Furthermore, 10–70 % of Kp isolates are serologically non-typeable, because either they express a novel capsule (common for clinical isolates) or they are non-capsulated (Cryz et al., 1986; Jenney et al., 2006; Tsay et al., 2002).

Kp employ a Wzy-dependent capsule synthesis process (Rahn et al., 1999; Whitfield, 2006) for which the associated genes are in the capsule synthesis locus (K-locus), which is 10–30 kbp in length (Arakawa et al., 1995; Chuang et al., 2006; Fevre et al., 2011; Pan et al., 2008, 2015; Shu et al., 2009). The K-locus includes a set of common genes in the terminal regions that encode the core capsule biosynthesis machinery (e.g. galF, wzi, wza, wzb, wzc, gnd and ugd). The central region is highly variable, encoding the capsule-specific sugar synthesis, processing and export proteins, plus the core assembly components Wzx (flippase) and Wzy (capsule repeat unit polymerase) (Pan et al., 2015).

K-locus nucleotide sequences and annotations are now available for a large number of Kp isolates, including the 77 K-type reference strains (Chuang et al., 2006; Deleo et al., 2014; Follador et al., 2016; Pan et al., 2008, 2015; Shu et al., 2009; The et al., 2015; Wyres et al., 2015). Serological K-types are generally defined by distinct sets of protein-encoding genes in the variable central region of the K-locus; however, two types (K22 and K37) are distinguished by a point mutation resulting in a premature stop codon that affects acetyltransferase function (Pan et al., 2015).

A number of molecular K-typing schemes have been developed, including RFLP (‘C-typing’) (Brisse et al., 2004), wzi and wzc typing (Brisse et al., 2013; Pan et al., 2013), and capsule-specific wzy PCR-based typing (Pan et al., 2008; Yu et al., 2007). These methods are less technically challenging and more discriminatory than serological techniques (Brisse et al., 2004, 2013; Pan et al., 2013). None have been widely adopted, although the single gene wzi and wzc typing schemes are gaining traction in the high-throughput sequencing era (Bialek-Davenet et al., 2014; Bowers et al., 2016; Zhou et al., 2016). Within the wzi scheme, unique alleles are associated with specific K-types (Brisse et al., 2013). Within the wzc scheme, K-types are assigned based on the level of wzc nucleotide similarity to a reference sequence, with a threshold of 94 % (Pan et al., 2013). Regardless of the method, a substantial proportion of isolates remain non-typeable; consequently, the true extent of Kp capsule diversity remains unknown.

Here, we report the K-loci from a collection of 2503 Kp. We have identified 31 novel K-loci, and have provided evidence that limited diversity remains to be discovered. We have defined a standardized nomenclature, provided a curated K-locus reference database and introduced Kaptive, a tool for rapid identification of reference K-loci from genome data, which will greatly facilitate surveillance efforts and evolutionary investigations of this important pathogen.

Methods

Sequences.

We obtained a total of 2600 Kp genomes (2021 publicly available and 579 novel genomes from a diverse set collected in Australia). Sequence reads were generated locally or obtained from the European Nucleotide Archive (accession numbers are listed in Table S1, available with the online Supplementary Material); 916 genomes that were publicly available as assembled contigs only were downloaded from PATRIC (Wattam et al., 2014) and the NCTC3000 Project (Wellcome Trust Sanger Institute – http://www.sanger.ac.uk/resources/downloads/bacteria/nctc/). For isolates sequenced in this study (n=579), DNA was extracted and libraries prepared using the Nextera XT 96 barcode DNA kit and 125 bp paired-end reads were generated on the Illumina HiSeq 2500 platform.

All paired-end read sets were filtered for a mean Phred quality score ≥30, then assembled de novo using SPAdes v3 (Bankevich et al., 2012). Genomes were excluded from the study if they were duplicate samples, or if there was evidence of contamination or mixed culture measured by: (i) <50 % reads mapping to the NTUH-K2044 reference chromosome (accession number: AP006725.1); (ii) the ratio of heterozygous/homozygous single nucleotide polymorphism (SNP)calls compared to the reference chromosome exceeding 20 %; (iii) the total assembly length being >6.5 Mb, or >6.0 Mb with evidence of >1 % non-Klebsiella read contamination as determined by MetaPhlAn (Segata et al., 2012); or (v) the assembly being low quality, i.e. total length <5 Mb.

Existing high-quality K-locus reference sequences.

K-locus nucleotide sequences for each of the 77 K-type references and 2 serologically non-typeable strains published elsewhere (Arakawa et al., 1995; Chuang et al., 2006; Fevre et al., 2011; Pan et al., 2008, 2015; Shu et al., 2009) were obtained from GenBank or directly from the authors (accession numbers are shown in Table S2). A total of 12 additional K-locus sequences had been published prior to the K-type references (Chen et al., 2014; Deleo et al., 2014; The et al., 2015; Wyres et al., 2015); we have compared these loci to those of the 77 K-type references (Arakawa et al., 1995; Chuang et al., 2006; Fevre et al., 2011; Pan et al., 2008, 2015; Shu et al., 2009) and identified 7 that were novel. These 7 novel loci plus 17 distinct loci described in our recent survey (Follador et al., 2016) were added to the non-redundant list of K-locus reference sequences, resulting in a total of 103 loci (see Table S2).

Identification of novel K-loci.

In order to identify novel K-loci, we first classified each genome by similarity to previously known loci. blastn (Camacho et al., 2009) was used to search each genome assembly for sequences with similarity to those of annotated K-locus coding sequences (CDSs) usually located between galF and ugd inclusive (minimum coverage 80 %, minimum identity 50 %). Transposase CDSs present in the published K-locus reference sequences were excluded from this analysis since they are not K-locus specific. Up to three missing CDSs were tolerated for K-locus assignment, to allow for assembly problems and insertion sequence (IS) insertions (see the Supplementary Methods and Fig. S1). This approach successfully distinguished the 77 K-type reference loci (with the exception of K22 and K37).

Genomes that could not be assigned a K-locus were investigated further: blastn was used to identify the galF and ugd genes within the assembly, and single contig loci were extracted. The assembly graph viewer Bandage (Wick et al., 2015) was used to identify K-loci that did not assemble on a single contig or where galF and/or ugd where missing. Loci were clustered, with identity and coverage thresholds of 90 %, using cd-hit-est (Fu et al., 2012; Li & Godzik, 2006). A representative sequence from each cluster was annotated with prokka (Seemann, 2014), using all proteins in the 77 reference K-type loci as the primary annotation database. Novel K-locus sequences were deposited in GenBank (accession numbers LT603702–LT603735; also included in the Kaptive database at https://github.com/katholt/Kaptive/).

Recombination in K-loci was investigated by aligning nucleotide sequences for the eight common genes (extracted from the reference annotations) using muscle (Edgar, 2004). This generated a 9944 bp concatenated gene alignment that was used as input for maximum likelihood (ML) phylogenetic inference with RAxML (Stamatakis, 2006) (best scoring tree from five runs each of 1000 bootstrap replicates with gamma model of rate heterogeneity), and recombination analysis using ClonalFrameML (Didelot & Wilson, 2015) (run for 1000 simulations and using the ML phylogeny as the starting tree).

Amino acid clustering.

Predicted amino acid sequences of all annotated K-locus coding regions were translated from the DNA sequences using BioPython and clustered with cd-hit (Fu et al., 2012; Li & Godzik, 2006) (90, 80, 70, 60, 50, 40 % identity). We explored the co-occurrence of predicted protein clusters present in three or more K-loci each (excluding the common proteins and the initiating glycosyltransferases, WbaP and WcaJ, n=115 clusters for analysis): pairwise Jaccard similarity scores were calculated as J (A, B)=A∩B/A[union or logical sum]B and were used to draw a weighted edge graph with the igraph R package v 1.0.1 (Csardi & Nepusz, 2006). A weight threshold was determined empirically as 0.61 and all edges for which J<0.61 were removed.

wzc and wzi nucleotide sequence determination.

We used srst2 (Inouye et al., 2014) to determine wzi alleles defined in the Kp BIGSdb (Institut Pasteur Klebsiella pneumoniae BIGSdb). blastn was used to determine alleles in genomes available only as assemblies. Novel alleles were submitted to the Kp BIGSdb for official designation. wzc sequences were extracted from genome assemblies by blastn search against a database of published alleles (Pan et al., 2013). Sequences were aligned with muscle (Edgar, 2004) and pairwise nucleotide divergences calculated.

Kaptive, a tool for identification of K-loci in genome data.

We developed an extended procedure for identification and assessment of full-length K-loci among bacterial genomes based on blast analysis of assemblies. The procedure has been automated and is implemented in a freely available open source software tool, Kaptive (https://github.com/katholt/Kaptive). For full details see Supplementary Methods and Results.

Comparison of Kaptive, wzi and wzc typing results.

Kaptive, wzi and wzc typing were applied to 86 genomes that had matched serological typing information available (Table S3). wzi alleles were determined as described above, and used to predict serotypes by comparison to the Kp BIGSdb. wzc sequences were extracted as above and genomes were assigned to serotypes if the sequence was <6 % divergent from the corresponding reference (Pan et al., 2013).

Results

Identification of novel K-loci

We investigated 2503 genomes that passed our quality-control standards (see Methods): 2298 K. pneumoniae sensu stricto, 144 K. variicola, 57 K. quasipneumoniae and 4 unclassified Klebsiella spp. (Tables 1 and S1). Also included were 10 publicly available genomes representing the more distantly related Klebsiella oxytoca (Table S4), as we hypothesized that they may share K-loci with Kp. Isolates had been collected between 1932 and 2014 (Fig. S2a), and from eight geographical regions spanning six continents (Fig. S2b).

Table 1.
Kp genomes investigated in this study

A total of 1371 genomes could be putatively assigned to 63 of the 77 K-type reference loci by the blast screening approach. A further 918 genomes were assigned to 1 of 25 previously published K-loci that are distinct from the K-type reference loci. Among the remaining 213 genomes, 106 were assigned to 29 novel K-loci, bringing the total to 132 (Table S1). Nine genomes harboured deletion or IS variants of known or novel loci (see below). For 93 genomes (3.7 %), no K-locus could be determined; however, we found three or more K-locus-associated genes in all such genomes and hypothesize the lack of assignment was attributable to low read depth and/or fragmented genome assembly rather than complete K-locus deletion.

The K66 and K74 loci were identified in one and four K. oxytoca genomes, respectively. Two novel K-loci were identified from three K. oxytoca genomes, increasing the total known K-loci to 134 (Table S2).

We estimated the extent to which we had captured the repertoire of K-locus diversity in the Kp population (Fig. 1). Rarefaction curves were estimated from: (i) the full genome set for which K-loci were assigned (n=2410); (ii) a ‘non-redundant’ genome set from which highly biased sub-samples such as outbreaks were removed (n=1081); and (iii) genomes from the non-redundant set representing each of the distinct species, K. variicola, K. quasipneumoniae and K. pneumoniae sensu stricto. In comparison to that for the full genome set, the non-redundant curve better represents the true diversity of the Kp population. Note that neither reached the total number of known K-loci, since 13 of the serologically defined K-loci (Pan et al., 2015) were not represented in our 2503 genomes. The rarefaction curves for each of the three Klebsiella species within the non-redundant dataset were similar to one another, indicating comparable levels of capsule diversity within each species (Fig. 1).

Fig. 1.
Rate of discovery of distinct K-loci with increasing genome sample size. Curves indicate the accumulation of distinct K-loci (mean±se) in different genome sets: grey, all genomes (n=2410, excluding K. oxytoca); blue, non-redundant genomes (n=1081, ...

K-locus nomenclature and reference database

We used a standardized K-locus nomenclature based on that proposed for Acinetobacter baumannii (Kenyon & Hall, 2013). Each distinct K-locus was designated as KL (K-locus) and a unique number. The K-type reference K-loci were assigned the same number as the corresponding K-type, e.g. K1 is encoded by the KL1 locus. K-loci for which K-types have not yet been phenotypically defined were assigned identifiers starting from 101 (note KL101 and KL102 correspond to loci previously named KN1 and KN2).

K-loci with IS insertions were distinguished from orthologous IS-free variants by using −1,–2. This nomenclature was consistently applied to the 10 K-type reference K-loci published elsewhere that include one or more ISs (Fevre et al., 2011; Pan et al., 2008, 2015; Shu et al., 2009). Deletion variants derived from a known K-locus were given the suffix -D1, -D2, etc. A complete Klebsiella K-locus reference database is available at https://github.com/katholt/Kaptive (see Table S2 and Supplementary Results for details).

K-locus epidemiology

The genome collection studied here does not represent a systematic sample and, thus, is not appropriate for detailed exploration of epidemiological questions. However, it is of interest to note the following three points. (i) Among the non-redundant genome set (n=1081), 90 % of K-loci identified were represented by just 67 distinct types, suggesting that many K-loci are rare in the Kp population, while others are more common. The five most common K-loci accounted for 20 % of those identified (KL2 5 %, KL1 4 %, KL21 4 %, KL17 4 %, KL30 3 %; see Table S1). (ii) There is evidence that K-types are shared across environmental and source populations. Twenty-eight K-loci were identified among isolates from retail meats (chicken, pork and turkey) for which genomes were originally published by Davis et al. (2015), and we found all (100 %) of these K-loci among human infection isolates in our data set. A total of 37 distinct K-loci were identified for 51 isolates from bovine hosts in the Kp global diversity study (Holt et al., 2015), and we found 35 (94.6 %) of these among human isolates and 12 (32.4 %) among the retail meat isolates. (iii) Our data set included 644 genomes from the globally distributed, multi-drug resistant clonal group 258, which was previously shown to contain extensive K-locus diversity (Bowers et al., 2015; Wyres et al., 2015). We identified a total of 32 distinct K-loci, plus 1 IS and 1 deletion variant amongst clonal group 258 genomes (Table S1). The most common K-loci were KL106 (n=125, 19.4 %) and KL107 (n=437, 67.9 %), which match the loci previously cited in the literature in association with wzi alleles 29 and 154, respectively (Bowers et al., 2015; Deleo et al., 2014; Wyres et al., 2015).

K-locus structures

The novel K-loci identified in this study conformed to the common structure described elsewhere (Fig. S3) (Pan et al., 2015; Rahn et al., 1999; Shu et al., 2009). We also identified six deletion variants: KL5-D1, KL20-D1, KL30-D1, KL62-D1, KL106-D1 and KL107-D1. Each of these was missing several common genes, but the remaining regions were homologous to existing K-loci represented in our genome collection. We suggest that the latter K-loci represent the ancestral forms that have subsequently lost one or more regions through deletion events; thus, generating the variants described here. Isolate NCTC10004 (recorded as serotype K11 in the UK National Collection of Type Cultures) and four other genomes carried K-loci that were nearly identical to the previously published K11 reference sequence (Pan et al., 2015). However, the latter lacked the essential wzx gene plus two other neighbouring genes, and was not identified among any other genomes. We assume the NCTC10004 locus represents the full length KL11 locus and designate the original K11 reference as KL11-D1 (it is unclear whether the original sequenced reference isolate had retained the ability to produce a capsule, since the serological typing was performed decades earlier) (Pan et al., 2015).

In four of the deletion variants, the deleted region was replaced by an IS, which may have mediated the deletion. Of the other IS-related variants, KL157-1 contained an IS903 family IS without an obvious deletion. In addition, we identified two novel IS variants of K-type reference loci (KL15-1 and KL22-1) and five IS-free variants of K-type reference loci (KL3, KL6, KL38, KL57 and KL81), plus one other previously published K-locus (KL103). The KL22-1 locus included a translocation of part of the nearby lipopolysaccharide (LPS) locus to the centre of the K-locus, plus an inversion of the 3′ K-locus region (Fig. 2). The translocated and inverted regions were bound at each end by a copy of ISKpn26.

Fig. 2.
Example K-locus structures and comparisons. CDSs are represented as arrows coloured by predicted function of the protein product and labelled with gene names where known. Grey bars indicate regions of similarity identified by blast comparison, darker ...

We used ClonalFrameML (Didelot & Wilson, 2015) to identify putative recombination events within the common K-locus genes. Analysis of the nucleotide sequences from each of the 134 reference K-loci identified a high number of such events (n=382), which were not distributed equally across the nucleotide alignment (Fig. 3a); rather the genes closest to the central variable region were affected by a greater number of recombination events compared to those at the ends of the locus.

Fig. 3.
Composition and diversity of Klebsiella K-loci. (a) Putative recombination events among the common K-locus genes. Values plotted are the relative number of events per 100 bp window, inferred using ClonalFrameML. (b) Representation of a generalized K-locus ...

Variation in K-locus gene content

A total of 2675 predicted proteins from 134 complete K-loci were clustered using cd-hit (Fu et al., 2012; Li & Godzik, 2006). As the identity threshold was reduced, the number of clusters continued to fall (from 1496 at 90 % identity to 508 at 40 % identity) and showed no signs of stabilizing (Fig. S4). At 40 % identity, which we believe is the lower bound for sensible comparison, the core capsule assembly proteins GalF, Wzi, Wza, Wzb, Gnd and Ugd each formed a single cluster and were present in nearly all loci (Fig. 3b). The Wzc sequences clustered into two groups and each locus encoded one Wzc protein (except KL50). In contrast, Wzx (flippase) clustered into 42 groups and Wzy (capsule repeat unit polymerase) clustered into 83 groups, highlighting the extreme diversity of these proteins compared to other core capsule assembly machinery proteins (Fig. 3c, d).

There were 374 clusters among the remaining proteins, most of which were associated with sugar synthesis and processing (Fig. 3). The initiating sugar transferase proteins, WbaP and WcaJ, were grouped into two clusters. These proteins are considered essential for capsule synthesis. Concordantly, each locus encoded a single protein from one of these two clusters. RmlB, RmlA, RmlD and RmlC, which are associated with rhamnose and typically encoded together in a single operon, were each represented by a single cluster. Similarly, the mannose synthesis and processing proteins, ManC and ManB, were grouped into a single cluster each. The associated operons rmlBADC and manCB were present in 55 and 73 K-loci, respectively (14 contained both operons; Fig. 3). In contrast, 360 of the remaining 366 protein clusters were present in fewer than ten K-loci each (Fig. 3e).

Co-occurrence analysis identified 18 correlated groups of K-locus proteins, each comprising two to five protein clusters (pairwise Jaccard similarity ≥0.61 for all pairs in the group; Fig. 4, Table S5). One group included the four Rml protein clusters; interestingly, this group also included a WcaA glycosyltransferase, which was present in 67.3 % of rmlBADC-containing K-loci and no rmlBADC-negative loci (χ2=70.09, P value <2.2×10−16, two-sided proportion test). Similarly, another group included the ManCB proteins and the putative mannosyl transferase, WbaZ, which was present in 65.8 % of manCB-containing K-loci and one manCB-negative locus (χ2=56.159, P value=6.683×10−14, two-sided proportion test). In addition, several groups included proteins for which the genes were located sequentially in their K-loci (e.g. wckG, wckH and wzx in KL12, KL29 and KL42) consistent with linked gene transfer.

Fig. 4.
Co-occurrence of wzx, wzy and sugar processing genes across reference K-loci. Nodes represent genes, labelled by name and coloured by protein product as in Figs 2 and and3:3: Wzx (red), Wzy (orange), mannose synthesis/processing proteins (dark ...

Diversity of wzc and wzi gene sequences

We sought to explore the utility of the existing wzi and wzc molecular capsule typing methods by characterizing wzi and wzc nucleotide sequence diversity and their association with K-loci. We confidently assigned wzi alleles to 2461 Kp genomes, including 390 distinct alleles, 218 of which were novel. Median pairwise nucleotide divergence was 7 %. Among the non-redundant genome set, there were 54 wzi alleles represented by at least five genomes, and of these 15 (28 %) were associated with more than one K-locus (Table S1). Among the 67 K-loci for which we had ≥5 representative sequences, 64 (95.5 %) were associated with two or more wzi alleles, and there was a general trend towards increasing wzi allelic diversity with increasing K-locus representation (Fig. 5).

Fig. 5.
Allelic diversity of wzi. Within K-locus wzi allelic diversity increases with total K-locus representation. The blue line represents the least-squares regression and grey shading indicates the 95 % confidence interval.

We extracted wzc sequences from 1041 of 1081 genomes in the non-redundant set (Fig. 6). In general, genomes sharing the same K-locus (6262 pairwise observations) showed lower wzc nucleotide divergence than those with different K-loci (491 775 pairwise observations), but the distributions overlapped substantially (Fig. 6). Notably, there were five distinct combinations of K-loci for which one or more pairs harboured wzc sequences that were <6 % divergent [the cut-off for K-type assignment as described in Pan et al. (2013)]: KL1 and KL112, KL9 and KL45, KL15 and KL52, KL30 and KL104, KL40 and KL135. Conversely, two K-loci (KL45, KL112) had more than 25 % wzc nucleotide divergence between representatives of the same K-locus.

Fig. 6.
wzc nucleotide diversity. Barplots showing distribution of pairwise wzc nucleotide divergence for pairs of genomes with the same (light blue) or different (dark blue) K-loci. The inset shows a magnified view of the lower end of the distribution.

Kaptivecapsule locus (K-locus) typing and variant evaluation from genome data

To facilitate easy identification of K-loci from genome assemblies, we developed the software tool Kaptive (see Fig. 7 and Supplementary Methods). We used Kaptive with our primary K-locus reference database (excluding IS and deletion variants) to rapidly type the K-loci in our collection of 2503 Kp genomes, and obtained confident K-locus calls for 2412 genomes (96.4 %, see Supplementary Results for further details).

Fig. 7.
Summary of the Kaptive analysis procedure. Kaptive takes as input a set of annotated reference K-loci in GenBank format and one or more genome assemblies, each as a single fasta file of contigs. Kaptive performs a series of blast searches to identify ...

We compared the K-locus calls from Kaptive, wzc and wzi typing to serological typing results for 86 isolates for which both genome and serology data were available (Holt et al., 2015; Jenney et al., 2006; NCTC3000 Project; Table S3). Five of six isolates that were non-typeable by serological techniques were identified by Kaptive as carrying KL16, KL54, KL81, KL111 and KL149. The KL16, KL54 and KL81 calls were in agreement with wzc and wzi typing results; the other two K-loci were not present in the wzi or wzc schemes and so were not typeable by those methods. Among the 80 serologically typeable isolates, the three molecular methods were generally in agreement with one another, although concordance with recorded phenotypes was quite low (65–74 %, Table S3). Call rates were highest for Kaptive (95 %), followed by wzc (89 %) and wzi typing (75 %).

Discussion

The number of distinct Klebsiella K-loci (now 134) is striking and exceeds that described for K-loci in other bacterial species, such as A. baumannii and Streptococcus pneumoniae. Furthermore, the diversity is an order of magnitude greater than that recently described for Klebsiella LPS, the other major Klebsiella surface antigen (Follador et al., 2016). This suggests that the K-locus is subject to strong diversifying selection. Given that these bacteria are not obligate pathogens and are ubiquitous in non-host-associated environments (Bagley, 1985; Podschun & Ullmann, 1998), the factors driving selection may not be host immune pressures but may include phage and/or protist predation.

Two novel K-loci, plus the KL66 and KL74 loci, were identified from K. oxytoca, a close relative of Kp. Little is known about K. oxytoca capsules, but one previous report also identified several Kp-associated capsules among K. oxytoca isolates (Ishihara et al., 2012). These findings indicate that K. oxytoca is able to exchange genetic material with Kp and, thus, represents a potential reservoir of virulence and other genes.

Our analysis confirms there are strong constraints on the structure of K-loci (Fig. 3). However, our data also reveal the extensive diversity of proteins encoded in the variable central region. The associated genes ranged in frequency from 0.7 to 54.7 % of the K-loci. Among those represented in at least three loci, approximately half co-occurred in groups ranging from two to five genes.

The molecular evolution driving K-locus diversification is not well understood, but likely includes a combination of point mutation, IS-mediated rearrangements and homologous recombination within the locus. In A. baumannii (Holt et al., 2016; Schultz et al., 2016) and S. pneumoniae (Wyres et al., 2013), it has been shown that recombination within the K-locus can drive capsule exchange between distinct clones. We recently speculated that this may also be true for Kp (Wyres et al., 2015) and the recombination analysis presented here supports this theory. The genes closest to the central variable region of the K-locus (i.e. wzb, wzc and gnd) showed evidence of the greatest number of recombination events, consistent with the hypothesis that they act as regions of homology for recombination events that shuffle the central region of the locus.

Prediction of capsule phenotypes from genome data is complex, as capsule expression is highly regulated and involves genes outside the K-locus region (Hsu et al., 2011). However, it is likely that K-loci encoding distinct sets of proteins are associated with distinct phenotypes, as is the case for the vast majority of K-type reference strains (Pan et al., 2015). Therefore, our data suggest that there are at least 134 distinct Klebsiella K-types. This is a lower bound estimate, particularly because our analysis did not capture differences that may arise from point mutations and small-scale insertions or deletions [e.g. in the case of K22 and K37, described by Pan et al. (2015)]. Furthermore, while we did not attempt to thoroughly characterize IS variants, several were apparent. The potential functional impacts of IS insertions likely depend on their location in the locus, but may include up-regulation, loss of capsule production or more subtle changes in sugar structures (Bentley et al., 2006; Hsu et al., 2016; Salter et al., 2012; Uria et al., 2008).

Serological typing of Klebsiella isolates is notoriously difficult and rarely performed. We were able to compare genotypes (whole-locus typing using Kaptive, as well as wzi and wzc typing schemes) with phenotypes for just 86 isolates for which both sequences and serotypes were available. Of the 19 discordant genotype versus phenotype results, 2 were due to deletion variants and were resolved by running Kaptive with the K-locus variants database. Interestingly, one of these isolates was non-typeable by serology, wzi or wzc typing, but recognized as a specific K-locus deletion variant by Kaptive. This highlights a benefit of our whole-locus typing approach; it provides epidemiologically relevant information even when the K-locus is interrupted. Another isolate was serologically typed as K54 but genotyped by Kaptive as KL113, which has sequence homology with KL54 (>84 % nucleotide identity over 76 % of the locus) and may encode a serologically similar or cross-reacting capsule. The other cases of discordance had no obvious explanation and may result from serological typing errors or from mutations arising during subculture (as identified for the K11 reference isolate above), neither of which we were able to check. Some discordance may also be due to unpredictable serological cross-reactions.

Given the problems with serotyping and the comparative robustness and widespread access to genome sequencing, we anticipate that genotyping will remain the preferred method for tracking capsular diversity in Klebsiella. Due to the extensive diversity and potential for ongoing evolution, we strongly advocate for classification based on complete, or near complete K-locus sequences, rather than single genes such as wzi or wzc, which can be misleading due to substitutions and horizontal gene transfer such as that described in this work. Kaptive analyses the full-length K-locus nucleotide sequence and assesses the presence of all K-locus-associated genes by protein blast search; thus, the approach is resilient to spurious results that may arise due to sequence divergence. Furthermore, the information provided allows users to determine confidence in the results and to identify putative novel K-loci or variants of known loci if desired. Along with the curated reference databases, this new tool will greatly facilitate evolutionary investigations and genomic surveillance efforts for Kp and other bacterial pathogens.

Acknowledgements

This work was funded by the National Health and Medical Research Council (NHMRC) of Australia (fellowship no. 1061409 and project no. 1043822 to K. E. H.) and Wellcome Trust (grant no. 098051 to the Wellcome Trust Sanger Institute).

Supplementary Data

626Supplementary File 1:

Supplementary Data

900Supplementary File 2:

Supplementary Data

61Supplementary File 3:

Supplementary Data

101Supplementary File 4:

Supplementary Data

147Supplementary File 5:

Notes

Abbreviations:

CDS
coding sequence
K-locus
capsule synthesis locus
K-type
capsule type
Kp
K. pneumoniae sensu stricto, K. variicola and K. quasipneumoniae
IS
insertion sequence
LPS
lipopolysaccharide
ML
maximum likelihood

References

  • Arakawa Y., Wacharotayankun R., Nagatsuka T., Ito H., Kato N., Ohta M.(1995). Genomic organization of the Klebsiella pneumoniae cps region responsible for serotype K2 capsular polysaccharide synthesis in the virulent strain chedid. J Bacteriol 1771788–1796.10.1128/jb.177.7.1788-1796.1995 [PMC free article] [PubMed] [Cross Ref]
  • Bagley S. T.(1985). Habitat association of Klebsiella species. Infect Contr 652–58. [PubMed]
  • Bankevich A., Nurk S., Antipov D., Gurevich A. A., Dvorkin M., Kulikov A. S., Lesin V. M., Nikolenko S., I, Pham S., et al. (2012). SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol 19455–477.10.1089/cmb.2012.0021 [PMC free article] [PubMed] [Cross Ref]
  • Bentley S. D., Aanensen D. M., Mavroidi A., Saunders D., Rabbinowitsch E., Collins M., Donohoe K., Harris D., Murphy L., et al. (2006). Genetic analysis of the capsular biosynthetic locus from all 90 pneumococcal serotypes. PLoS Genet 2e31. [PubMed]
  • Bialek-Davenet S., Criscuolo A., Ailloud F., Passet V., Jones L., Delannoy-Vieillard A. S., Garin B., Le Hello S., Arlet G., et al. (2014). Genomic definition of hypervirulent and multidrug-resistant Klebsiella pneumoniae clonal groups. Emerg Infect Dis 201812–1820.10.3201/eid2011.140206 [PMC free article] [PubMed] [Cross Ref]
  • Bowers J. R., Kitchel B., Driebe E. M., MacCannell D. R., Roe C., Lemmer D., de Man T., Rasheed J. K., Engelthaler D. M., et al. (2015). Genomic analysis of the emergence and rapid global dissemination of the clonal group 258 Klebsiella pneumoniae pandemic. PLoS One 10e0133727.10.1371/journal.pone.0133727 [PMC free article] [PubMed] [Cross Ref]
  • Bowers J. R., Lemmer D., Sahl J. W., Pearson T., Driebe E. M., Engelthaler D. M., Keim P.(2016). KlebSeq: a diagnostic tool for healthcare surveillance and antimicrobial resistance monitoring of Klebsiella pneumoniae. J Clin Microbiol 542582–2596. [PMC free article] [PubMed]
  • Brisse S., Issenhuth-Jeanjean S., Grimont P. A.(2004). Molecular serotyping of Klebsiella species isolates by restriction of the amplified capsular antigen gene cluster. J Clin Microbiol 423388–3398.10.1128/JCM.42.8.3388-3398.2004 [PMC free article] [PubMed] [Cross Ref]
  • Brisse S., Passet V., Haugaard A. B., Babosan A., Kassis-Chikhani N., Struve C., Decré D.(2013). wzi gene sequencing, a rapid method for determination of capsular type for Klebsiella strains. J Clin Microbiol 514073–4078.10.1128/JCM.01924-13 [PMC free article] [PubMed] [Cross Ref]
  • Camacho C., Coulouris G., Avagyan V., Ma N., Papadopoulos J., Bealer K., Madden T. L.(2009). blast+: architecture and applications. BMC Bioinformatics 10421.10.1186/1471-2105-10-421 [PMC free article] [PubMed] [Cross Ref]
  • Campos M. A., Vargas M. A., Regueiro V., Llompart C. M., Albertí S., Bengoechea J. A., José A.(2004). Capsule polysaccharide mediates bacterial resistance to antimicrobial peptides. Infect Immun 727107–7114.10.1128/IAI.72.12.7107-7114.2004 [PMC free article] [PubMed] [Cross Ref]
  • Chen L., Mathema B., Pitout J. D., DeLeo F. R., Kreiswirth B. N.(2014). Epidemic Klebsiella pneumoniae ST258 is a hybrid strain. MBio 5e01355-14.10.1128/mBio.01355-14 [PMC free article] [PubMed] [Cross Ref]
  • Chuang Y. P., Fang C. T., Lai S. Y., Chang S. C., Wang J. T.(2006). Genetic determinants of capsular serotype K1 of Klebsiella pneumoniae causing primary pyogenic liver abscess. J Infect Dis 193645–654.10.1086/499968 [PubMed] [Cross Ref]
  • Cortés G., Borrell N., Astorza B., Gómez C., Sauleda J., Albertí S.(2002). Molecular analysis of the contribution of the capsular polysaccharide and the lipopolysaccharide O side chain to the virulence of Klebsiella pneumoniae in a murine model of pneumonia. Infect Immun 702583–2590.10.1128/IAI.70.5.2583-2590.2002 [PMC free article] [PubMed] [Cross Ref]
  • Cryz S. J., Mortimer P. M., Mansfield V., Germanier R.(1986). Seroepidemiology of Klebsiella bacteremic isolates and implications for vaccine development. J Clin Microbiol 23687–690. [PMC free article] [PubMed]
  • Csardi G., Nepusz T.(2006). The igraph software package for complex network research. InterJournal, Complex Systems 2006,1695.
  • Davis G. S., Waits K., Nordstrom L., Weaver B., Aziz M., Gauld L., Grande H., Bigler R., Horwinski J., et al. (2015). Intermingled Klebsiella pneumoniae populations between retail meats and human urinary tract infections. Clin Infect Dis 61892–899.10.1093/cid/civ428 [PMC free article] [PubMed] [Cross Ref]
  • Deleo F. R., Chen L., Porcella S. F., Martens C. A., Kobayashi S. D., Porter A. R., Chavda K. D., Jacobs M. R., Mathema B., et al. (2014). Molecular dissection of the evolution of carbapenem-resistant multilocus sequence type 258 Klebsiella pneumoniae. Proc Natl Acad Sci U S A 1114988–4993.10.1073/pnas.1321364111 [PubMed] [Cross Ref]
  • Didelot X., Wilson D. J.(2015). ClonalFrameML: efficient inference of recombination in whole bacterial genomes. PLoS Comp Biol 11e100404110.1371/journal.pcbi.1004041 [PMC free article] [PubMed] [Cross Ref]
  • Edgar R. C.(2004). MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res 321792–1797.10.1093/nar/gkh340 [PMC free article] [PubMed] [Cross Ref]
  • Edmunds P. N.(1954). Further Klebsiella capsule types. J Infect Dis 9465–71.10.1093/infdis/94.1.65 [PubMed] [Cross Ref]
  • Edwards P. R., Fife M. A.(1952). Capsule types of Klebsiella. J Infect Dis 9192–104. [PubMed]
  • Evrard B., Balestrino D., Dosgilbert A., Bouya-Gachancard J. L., Charbonnel N., Forestier C., Tridon A.(2010). Roles of capsule and lipopolysaccharide O antigen in interactions of human monocyte-derived dendritic cells and Klebsiella pneumoniae. Infect Immun 78210–219.10.1128/IAI.00864-09 [PMC free article] [PubMed] [Cross Ref]
  • Fevre C., Passet V., Deletoile A., Barbe V., Frangeul L., Almeida A. S., Sansonetti P., Tournebize R., Brisse S.(2011). PCR-based identification of Klebsiella pneumoniae subsp. rhinoscleromatis, the agent of rhinoscleroma. PLoS Negl Trop Dis 5e1052.10.1371/journal.pntd.0001052 [PMC free article] [PubMed] [Cross Ref]
  • Follador R., Heinz E., Wyres K. L., Ellington M. J., Kowarik M., Holt K. E., Thomson N. R.(2016). The diversity of Klebsiella pneumoniae surface polysaccharides. MGen 2,doi: 10.1099/mgen.0.000073
  • Fu L., Niu B., Zhu Z., Wu S., Li W.(2012). CD-HIT: accelerated for clustering the next-generation sequencing data. Bioinformatics 283150–3152.10.1093/bioinformatics/bts565 [PMC free article] [PubMed] [Cross Ref]
  • Holt K. E., Wertheim H., Zadoks R. N., Baker S., Whitehouse C. A., Dance D., Jenney A., Connor T. R., Hsu L. Y., et al. (2015). Genomic analysis of diversity, population structure, virulence, and antimicrobial resistance in Klebsiella pneumoniae, an urgent threat to public health. Proc Natl Acad Sci U S A 112E3574–E3581.10.1073/pnas.1501049112 [PubMed] [Cross Ref]
  • Holt K., Kenyon J. J., Hamidian M., Schultz M. B., Pickard D. J., Dougan G., Hall R.(2016). Five decades of genome evolution in the globally distributed, extensively antibiotic-resistant Acinetobacter baumannii global clone 1. MGen 2,10.1099/mgen.0.000052.10.1099/mgen.0.000052 [Cross Ref]
  • Hsu C.-R., Lin T.-L., Chen Y.-C., Chou H.-C., Wang J.-T.(2011). The role of Klebsiella pneumoniae rmpA in capsular polysaccharide synthesis and virulence revisited. Microbiology 1573446–3457.10.1099/mic.0.050336-0 [PubMed] [Cross Ref]
  • Hsu C.-R., Liao C.-H., Lin T.-L., Yang H.-R., Yang F.-L., Hsieh P.-F., Wu S.-H., Wang J.-T.(2016). Identification of a capsular variant and characterization of capsular acetylation in Klebsiella pneumoniae PLA-associated type K57. Sci Rep 631946.10.1038/srep31946 [PMC free article] [PubMed] [Cross Ref]
  • Inouye M., Dashnow H., Raven L. A., Schultz M. B., Pope B. J., Tomita T., Zobel J., Holt K. E.(2014). SRST2: rapid genomic surveillance for public health and hospital microbiology labs. Genome Med 690.10.1186/s13073-014-0090-6 [PMC free article] [PubMed] [Cross Ref]
  • Ishihara Y., Yagi T., Mochizuki M., Ohta M.(2012). Capsular types, virulence factors and DNA types of Klebsiella oxytoca strains isolated from blood and bile. Kansenshogaku Zasshi 861221–126. [PubMed]
  • Jenney A. W., Clements A., Farn J. L., Wijburg O. L., McGlinchey A., Spelman D. W., Pitt T. L., Kaufmann M. E., Liolios L., et al. (2006). Seroepidemiology of Klebsiella pneumoniae in an Australian tertiary hospital and its implications for vaccine development. J Clin Microbiol 44102–107.10.1128/JCM.44.1.102-107.2006 [PMC free article] [PubMed] [Cross Ref]
  • Kenyon J. J., Hall R. M.(2013). Variation in the complex carbohydrate biosynthesis loci of Acinetobacter baumannii genomes. PLoS One 8e62160.10.1371/journal.pone.0062160 [PMC free article] [PubMed] [Cross Ref]
  • Lawlor M. S., Hsu J., Rick P. D., Miller V. L.(2005). Identification of Klebsiella pneumoniae virulence determinants using an intranasal infection model. Mol Microbiol 581054–1073.10.1111/j.1365-2958.2005.04918.x [PubMed] [Cross Ref]
  • Lee C. H., Chang C. C., Liu J. W., Chen R. F., Yang K. D.(2014). Sialic acid involved in hypermucoviscosity phenotype of Klebsiella pneumoniae and associated with resistance to neutrophil phagocytosis. Virulence 5673–679.10.4161/viru.32076 [PMC free article] [PubMed] [Cross Ref]
  • Lee R., Molton J. S., Wyres K. L., Gorrie C., Wong J., Hoh C. H., Teo J., Kalimuddin S., Lye D. C., et al. (2016). Differential host susceptibility and bacterial virulence factors driving Klebsiella liver abscess in an ethnically diverse population. Sci Rep 1329316. [PMC free article] [PubMed]
  • Li W., Godzik A.(2006). CD-Hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics 221658–1659.10.1093/bioinformatics/btl158 [PubMed] [Cross Ref]
  • March C., Cano V., Moranta D., Llobet E., Pérez-Gutiérrez C., Tomás J. M., Suárez T., Garmendia J., Bengoechea J. A.(2013). Role of bacterial surface structures on the interaction of Klebsiella pneumoniae with phagocytes. PLoS One 8,e56847. [PMC free article] [PubMed]
  • Merino S., Camprubí S., Albertí S., Benedí V. J., Tomás J. M.(1992). Mechanisms of Klebsiella pneumoniae resistance to complement-mediated killing. Infect Immun 602529–2535. [PMC free article] [PubMed]
  • Munoz-Price L. S., Poirel L., Bonomo R. A., Schwaber M. J., Daikos G. L., Cormican M., Cornaglia G., Garau J., Gniadkowski M., et al. (2013). Clinical epidemiology of the global expansion of Klebsiella pneumoniae carbapenemases. Lancet Infect Dis 13785–796.10.1016/S1473-3099(13)70190-7 [PMC free article] [PubMed] [Cross Ref]
  • Ørskov I. D. A., Fife-Asbury M. A.(1977). New Klebsiella capsular antigen, K82, and the deletion of five of those previously assigned. Int J Syst Bacteriol 27386–387.
  • Pan Y.-J., Fang H.-C., Yang H.-C., Lin T.-L., Hsieh P.-F., Tsai F.-C., Keynan Y., Wang J.-T.(2008). Capsular polysaccharide synthesis regions in Klebsiella pneumoniae serotype K57 and a new capsular serotype. J Clin Microbiol 462231–2240.10.1128/JCM.01716-07 [PMC free article] [PubMed] [Cross Ref]
  • Pan Y.-J., Lin T.-L., Chen Y.-H., Hsu C.-R., Hsieh P.-F., Wu M.-C., Wang J.-T.(2013). Capsular types of Klebsiella pneumoniae revisited by wzc sequencing. PLoS One 8e80670.10.1371/journal.pone.0080670 [PMC free article] [PubMed] [Cross Ref]
  • Pan Y.-J., Lin T.-L., Chen C.-T., Chen Y.-Y., Hsieh P.-F., Hsu C.-R., Wu M.-C., Wang J.-T.(2015). Genetic analysis of capsular polysaccharide synthesis gene clusters in 79 capsular types of Klebsiella spp. Nat Sci Rep 515573. [PMC free article] [PubMed]
  • Podschun R., Ullmann U.(1998). Klebsiella spp. as nosocomial pathogens: epidemiology, taxonomy, typing methods, and pathogenicity factors. Clin Microbiol Rev 11589–603. [PMC free article] [PubMed]
  • Rahn A., Drummelsmith J., Whitfield C.(1999). Conserved organization in the cps gene clusters for expression of Escherichia coli group 1 K antigens: relationship to the colanic acid biosynthesis locus and the cps genes from Klebsiella pneumoniae. J Bacteriol 1812307–2713. [PMC free article] [PubMed]
  • Salter S. J., Gould K. A., Lambertsen L. M., Hanage W. P., Antonio M., Turner P., Hermans P. W. M., Bootsma H. J., O’Brien K. L., et al. (2012). Variation at the capsule locus, cps, of mistyped and non-typable Streptococcus pneumoniae isolates. Microbiology 1581560–1569.10.1099/mic.0.056580-0 [PMC free article] [PubMed] [Cross Ref]
  • Schultz M. B., Thanh D. P., Do Hoan N. T., Wick R. R., Ingle D. J., Hawkey J., Edwards D. J., Kenyon J. J., Lan N. P. H., Campbell J. I.(2016). Repeated local emergence of carbapenem resistant Acinetobacter baumannii in a single hospital ward. MGen 2,doi: 10.1099/mgen.0.000050
  • Seemann T.(2014). prokka: rapid prokaryotic genome annotation. Bioinformatics 302068–2069.10.1093/bioinformatics/btu153 [PubMed] [Cross Ref]
  • Segata N., Waldron L., Ballarini A., Narasimhan V., Jousson O., Huttenhower C.(2012). Metagenomic microbial community profiling using unique clade-specific marker genes. Nat Methods 9811–814.10.1038/nmeth.2066 [PMC free article] [PubMed] [Cross Ref]
  • Shu H. Y., Fung C. P., Liu Y. M., Wu K. M., Chen Y. T., Li L. H., Liu T. T., Kirby R., Tsai S. F.(2009). Genetic diversity of capsular polysaccharide biosynthesis in Klebsiella pneumoniae clinical isolates. Microbiology 1554170–4183.10.1099/mic.0.029017-0 [PubMed] [Cross Ref]
  • Snitkin E. S., Zelazny A. M., Thomas P. J., Stock F., NISC Comparative Sequencing Program, Henderson D. K., Palmore T. N., Segre J. A.(2012). Tracking a hospital outbreak of carbapenem-resistant Klebsiella pneumoniae with whole-genome sequencing. Sci Transl Med 4148ra116.10.1126/scitranslmed.3004129 [PMC free article] [PubMed] [Cross Ref]
  • Stamatakis A.(2006). RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics 222688–2690.10.1093/bioinformatics/btl446 [PubMed] [Cross Ref]
  • Stoesser N., Batty E. M., Eyre D. W., Morgan M., Wyllie D. H., Del Ojo Elias C., Johnson J. R., Walker A. S., Peto T. E. A., Crook D. W.(2013). Predicting antimicrobial susceptibilities for Escherichia coli and Klebsiella pneumoniae isolates using whole genomic sequence data. J Antimicrob Chemother 682234–2244.10.1093/jac/dkt180 [PMC free article] [PubMed] [Cross Ref]
  • Stoesser N., Giess A., Batty E. M., Sheppard A. E., Walker A. S., Wilson D. J., Didelot X., Bashir A., Sebra R., et al. (2014). Genome sequencing of an extended series of NDM-producing Klebsiella pneumoniae isolates from neonatal infections in a Nepali hospital characterizes the extent of community- versus hospital-associated transmission in an endemic setting. Antimicrob Agents Chemother 587347–7357.10.1128/AAC.03900-14 [PMC free article] [PubMed] [Cross Ref]
  • Struve C., Roe C. C., Stegger M., Stahlhut S. G., Hansen D. S., Engelthaler D. M., Andersen P. S., Driebe E. M., Keim P., Krogfelt A.(2015). Mapping the evolution of hypervirulent Klebsiella pneumoniae. MBio 6e00630-15. [PMC free article] [PubMed]
  • The H. C., Karkey A., Thanh D. P., Boinett C. J., Cain A. K., Ellington M., Baker K. S., Dongol S., Thompson C., et al. (2015). A high-resolution genomic analysis of multidrug- resistant hospital outbreaks of Klebsiella pneumoniae. EMBO Molec Med 7227–239. [PMC free article] [PubMed]
  • Tsay R.-W., Siu L. K., Fung C.-P., Chang F.-Y.(2002). Characteristics of bacteremia between community-acquired and nosocomial Klebsiella pneumoniae infection. Arch Intern Med 1621021–1027. [PubMed]
  • Uria M. J., Zhang Q., Li Y., Chan A., Exley R. M., Gollan B., Chan H., Feavers I., Yarwood A., et al. (2008). A generic mechanism in Neisseria meningitidis for enhanced resistance against bactericidal antibodies. J Exp Med 2051423–1434.10.1084/jem.20072577 [PMC free article] [PubMed] [Cross Ref]
  • Wand M. E., Baker K. S., Benthall G., McGregor H., McCowen J. W. I., Deheer-Graham A., Sutton J. M.(2015). Characterization of pre-antibiotic era Klebsiella pneumoniae isolates with respect to antibiotic/disinfectant susceptibility and virulence in Galleria mellonella. Antimicrob Agents Chemother 593966–3972.10.1128/AAC.05009-14 [PMC free article] [PubMed] [Cross Ref]
  • Wattam A. R., Abraham D., Dalay O., Disz T. L., Driscoll T., Gabbard J. L., Gillespie J. J., Gough R., Hix D., et al. (2014). PATRIC, the bacterial bioinformatics database and analysis resource. Nucleic Acids Res 42581–591. [PMC free article] [PubMed]
  • Whitfield C.(2006). Biosynthesis and assembly of capsular polysaccharides in Escherichia coli. Annu Rev Biochem 7539–68.10.1146/annurev.biochem.75.103004.142545 [PubMed] [Cross Ref]
  • Wick R. R., Schultz M. B., Zobel J., Holt K. E.(2015). Bandage: interactive visualization of de novo genome assemblies. Bioinformatics 313350–3352.10.1093/bioinformatics/btv383 [PMC free article] [PubMed] [Cross Ref]
  • Wyres K. L., Lambertsen L. M., Croucher N. J., McGee L., von Gottberg A., Liñares J., Jacobs M. R., Kristinsson K. G., Beall B. W., et al. (2013). Pneumococcal capsular switching: a historical perspective. J Infect Dis 207439–449.10.1093/infdis/jis703 [PMC free article] [PubMed] [Cross Ref]
  • Wyres K. L., Gorrie C., Edwards D. J., Wertheim H. F. L., Hsu L. Y., Van Kinh N., Zadoks R., Baker S., Holt K. E.(2015). Extensive capsule locus variation and large-scale genomic recombination within the Klebsiella pneumoniae clonal group 258. Genome Biol Evol 71267–1279.10.1093/gbe/evv062 [PMC free article] [PubMed] [Cross Ref]
  • Yoshida K., Matsumoto T., Tateda K., Uchida K., Tsujimoto S., Yamaguchi K.(2000). Role of bacterial capsule in local and systemic inflammatory responses of mice during pulmonary infection with Klebsiella pneumoniae. J Med Microbiol 491003–1010.10.1099/0022-1317-49-11-1003 [PubMed] [Cross Ref]
  • Yu W. L., Fung C. P., Ko W. C., Cheng K. C., Lee C. C., Chuang Y. C.(2007). Polymerase chain reaction analysis for detecting capsule serotypes K1 and K2 of Klebsiella pneumoniae causing abscesses of the liver and other sites. J Infect Dis 1951235–1236.10.1086/512686 [PubMed] [Cross Ref]
  • Zhou K., Lokate M., Deurenberg R. H., Tepper M., Arends J. P., Raangs E. G. C., Lo-Ten-Foe J., Grundmann H., Rossen J. W. A., Friedrich A. W.(2016). Use of whole-genome sequencing to trace, control and characterize the regional expansion of extended-spectrum β-lactamase producing ST15 Klebsiella pneumoniae. Sci Rep 620840.10.1038/srep20840 [PMC free article] [PubMed] [Cross Ref]

Data Bibliography

1. Gorrie, C., Jenney A. & Holt, K. E. NCBI BioProject PRJNA351909 (2016).
2. Einsiedel, L. & Holt, K. E. NCBI BioProject PRJNA351911 (2016).
3. Holt K. E. & Hall, R. NCBI BioProject PRJNA356346 (2016).
4. The Wellcome Trust Sanger Institute, European Nucleotide Archive PRJEB6891, pre-publication release.
5. Bialek-Davenet, S., Criscuolo, A., Ailloud, F., Passet, V., Jones, L., Delannoy-Vieillard, A. S., Garin, B., Le Hello, S., Arlet, G. and other authors. Genomic definition of hypervirulent and multidrug-resistant Klebsiella pneumoniae clonal groups. Emerg Infect Dis 20, 1812–1820. European Nucleotide Archive PRJEB6688 (2014). [PMC free article] [PubMed]
6. Bowers, J. R., Kitchel, B., Driebe, E. M., MacCannell, D. R., Roe, C., Lemmer, D., de Man, T., Rasheed, J. K., Engelthaler, D. M. and other authors. Genomic analysis of the emergence and rapid global dissemination of the clonal group 258 Klebsiella pneumoniae pandemic. PLoS One 10, e0133727. European Nucleotide Archive PRJNA252957 (2015). [PMC free article] [PubMed]
7. Davis, G. S., Waits, K., Nordstrom, L., Weaver, B., Aziz, M., Gauld, L., Grande, H., Bigler, R., Horwinski, J. and other authors. Intermingled Klebsiella pneumoniae populations between retail meats and human urinary tract infections. Clin Infect Dis 61, 892–899. European Nucleotide Archive PRJNA289272 (2015). [PMC free article] [PubMed]
8. Deleo, F. R., Chen, L., Porcella, S. F., Martens, C. A., Kobayashi, S. D., Porter, A. R., Chavda, K. D., Jacobs, M. R., Mathema, B. and other authors. Molecular dissection of the evolution of carbapenem-resistant multilocus sequence type 258 Klebsiella pneumoniae. Proc Natl Acad Sci U S A 111, 4988–4993. European Nucleotide PRJNA237670 (2014). [PubMed]
9. Ellington M. J. European Nucleotide Archive PRJEB1271 pre-publication release.
10. Holt, K. E., Wertheim, H., Zadoks, R. N., Baker, S., Whitehouse, C. A., Dance, D., Jenney, A., Connor, T. R., Hsu, L. Y. and other authors. (2015). Genomic analysis of diversity, population structure, virulence, and antimicrobial resistance in Klebsiella pneumoniae, an urgent threat to public health. Proc Natl Acad Sci U S A 112, E3574–E3581. European Nucleotide Archive PRJEB2111 (2015). [PubMed]
11. Lee, I. R., Molton, J. S., Wyres, K. L., Gorrie, C., Wong, J., Hoh, C. H., Teo, J., Kalimuddin, S., Lye, D. C. and other authors. Differential host susceptibility and bacterial virulence factors driving Klebsiella liver abscess in an ethnically diverse population. Sci Rep 6, 29316. NCBI BioProject PRJNA351910 (2016). [PMC free article] [PubMed]
12. Stoesser, N., Batty, E. M., Eyre, D. W., Morgan, M., Wyllie, D. H., Del Ojo Elias, C., Johnson, J. R., Walker, A. S., Peto, T. E. A. & Crook, D. W. Predicting antimicrobial susceptibilities for Escherichia coli and Klebsiella pneumoniae isolates using whole genomic sequence data. J Antimicrob Chemother 68, 2234–2244. European Nucleotide Archive PRJEB1963 (2013). [PMC free article] [PubMed]
13. Stoesser, N., Giess, A., Batty, E. M., Sheppard, A. E., Walker, A. S., Wilson, D. J., Didelot, X., Bashir, A., Sebra, R. and other authors. Genome sequencing of an extended series of NDM-producing Klebsiella pneumoniae isolates from neonatal infections in a Nepali hospital characterizes the extent of community- versus hospital-associated transmission in an endemic setting. Antimicrob Agents Chemother 58, 7347–7357. European Nucleotide Archive PRJNA253300 (2014). [PMC free article] [PubMed]
14. Struve, C., Roe, C. C., Stegger, M., Stahlhut, S. G., Hansen, D. S., Engelthaler, D. M., Andersen, P. S., Driebe, E. M., Keim, P. & Krogfelt, K. A. Mapping the evolution of hypervirulent Klebsiella pneumoniae. MBio 6, e00630-15. European Nucleotide Archive PRJEB7967 (2015). [PMC free article] [PubMed]
15. The, H. C., Karkey, A., Thanh, D. P., Boinett, C. J., Cain, A. K., Ellington, M., Baker, K. S., Dongol, S., Thompson, C. and other authors. A high-resolution genomic analysis of multidrug-resistant hospital outbreaks of Klebsiella pneumoniae. EMBO Molec Med 7, 227–239. European Nucleotide Archive PRJEB1800 (2015). [PMC free article] [PubMed]
16. Wand, M. E., Baker, K. S., Benthall, G., McGregor, H., McCowen, J. W. I., Deheer-Graham, A. & Sutton, J. M. Characterization of pre-antibiotic era Klebsiella pneumoniae isolates with respect to antibiotic/disinfectant susceptibility and virulence in Galleria mellonella. Antimicrob Agents Chemother 59, 3966–3972. European Nucleotide PRJEB3255 (2015). [PMC free article] [PubMed]
17. Wyres, K. L., Wick, R. R., Gorrie, C., Jenney, A., Follador, R., Thomson, N. R. & Holt, K. E. GenBank LT603702–LT603735 (2016).
18. Wick, R. R., Wyres, K. L. & Holt, K. E. Github, DOI: 10.5281/zenodo.55773.

Articles from Microbial Genomics are provided here courtesy of Microbiology Society