IS-capades of Klebsiella pneumoniae: Insertion sequences drive metabolic loss in obscure sub-lineages

IS-capades of Klebsiella pneumoniae: Insertion sequences drive metabolic loss in obscure sub-lineages | bioRxiv /* */ /* */ <!-- <!-- /*! * yepnope1.5.4 * (c) WTFPL, GPLv2 */ (function(a,b,c){function d(a){return"[object Function]"==o.call(a)}function e(a){return"string"==typeof a}function f(){}function g(a){return!a||"loaded"==a||"complete"==a||"uninitialized"==a}function h(){var a=p.shift();q=1,a?a.t?m(function(){("c"==a.t?B.injectCss:B.injectJs)(a.s,0,a.a,a.x,a.e,1)},0):(a(),h()):q=0}function i(a,c,d,e,f,i,j){function k(b){if(!o&&g(l.readyState)&&(u.r=o=1,!q&&h(),l.onload=l.onreadystatechange=null,b)){"img"!=a&&m(function(){t.removeChild(l)},50);for(var d in y[c])y[c].hasOwnProperty(d)&&y[c][d].onload()}}var j=j||B.errorTimeout,l=b.createElement(a),o=0,r=0,u={t:d,s:c,e:f,a:i,x:j};1===y[c]&&(r=1,y[c]=[]),"object"==a?l.data=c:(l.src=c,l.type=a),l.width=l.height="0",l.onerror=l.onload=l.onreadystatechange=function(){k.call(this,r)},p.splice(e,0,u),"img"!=a&&(r||2===y[c]?(t.insertBefore(l,s?null:n),m(k,j)):y[c].push(l))}function j(a,b,c,d,f){return q=0,b=b||"j",e(a)?i("c"==b?v:u,a,b,this.i++,c,d,f):(p.splice(this.i++,0,a),1==p.length&&h()),this}function k(){var a=B;return a.loader={load:j,i:0},a}var l=b.documentElement,m=a.setTimeout,n=b.getElementsByTagName("script")[0],o={}.toString,p=[],q=0,r="MozAppearance"in l.style,s=r&&!!b.createRange().compareNode,t=s?l:n.parentNode,l=a.opera&&"[object Opera]"==o.call(a.opera),l=!!b.attachEvent&&!l,u=r?"object":l?"script":"img",v=l?"script":u,w=Array.isArray||function(a){return"[object Array]"==o.call(a)},x=[],y={},z={timeout:function(a,b){return b.length&&(a.timeout=b[0]),a}},A,B;B=function(a){function b(a){var a=a.split("!"),b=x.length,c=a.pop(),d=a.length,c={url:c,origUrl:c,prefixes:a},e,f,g;for(f=0;f<d;f++)g=a[f].split("="),(e=z[g.shift()])&&(c=e(c,g));for(f=0;f<b;f++)c=x[f](c);return c}function g(a,e,f,g,h){var i=b(a),j=i.autoCallback;i.url.split(".").pop().split("?").shift(),i.bypass||(e&&(e=d(e)?e:e[a]||e[g]||e[a.split("/").pop().split("?")[0]]),i.instead?i.instead(a,e,f,g,h):(y[i.url]?i.noexec=!0:y[i.url]=1,f.load(i.url,i.forceCSS||!i.forceJS&&"css"==i.url.split(".").pop().split("?").shift()?"c":c,i.noexec,i.attrs,i.timeout),(d(e)||d(j))&&f.load(function(){k(),e&&e(i.origUrl,h,g),j&&j(i.origUrl,h,g),y[i.url]=2})))}function h(a,b){function c(a,c){if(a){if(e(a))c||(j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}),g(a,j,b,0,h);else if(Object(a)===a)for(n in m=function(){var b=0,c;for(c in a)a.hasOwnProperty(c)&&b++;return b}(),a)a.hasOwnProperty(n)&&(!c&&!--m&&(d(j)?j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}:j[n]=function(a){return function(){var b=[].slice.call(arguments);a&&a.apply(this,b),l()}}(k[n])),g(a[n],j,b,n,h))}else!c&&l()}var h=!!a.test,i=a.load||a.both,j=a.callback||f,k=j,l=a.complete||f,m,n;c(h?a.yep:a.nope,!!i),i&&c(i)}var i,j,l=this.yepnope.loader;if(e(a))g(a,0,l,0);else if(w(a))for(i=0;i (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0];var j=d.createElement(s);var dl=l!='dataLayer'?'&l='+l:'';j.src='//www.googletagmanager.com/gtm.js?id='+i+dl;j.type='text/javascript';j.async=true;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-M677548'); Skip to main content Home About Submit ALERTS / RSS Search for this keyword Advanced Search New Results IS-capades of Klebsiella pneumoniae : Insertion sequences drive metabolic loss in obscure sub-lineages Ben Vezina , Claire White , Helena B. Cooper , View ORCID Profile Kathryn E. Holt , View ORCID Profile Jane Hawkey , View ORCID Profile Kelly L. Wyres , Margaret M. C. Lam doi: https://doi.org/10.1101/2025.07.24.666535 Ben Vezina 1 Department of Infectious Diseases, School of Translational Medicine, Monash University , Melbourne, Victoria, Australia Find this author on Google Scholar Find this author on PubMed Search for this author on this site For correspondence: benjamin.vezina{at}monash.edu margaret.lam{at}monash.edu Claire White 1 Department of Infectious Diseases, School of Translational Medicine, Monash University , Melbourne, Victoria, Australia Find this author on Google Scholar Find this author on PubMed Search for this author on this site Helena B. Cooper 1 Department of Infectious Diseases, School of Translational Medicine, Monash University , Melbourne, Victoria, Australia 2 Centre to Impact AMR, Monash University , Clayton, Victoria, Australia Find this author on Google Scholar Find this author on PubMed Search for this author on this site Kathryn E. Holt 1 Department of Infectious Diseases, School of Translational Medicine, Monash University , Melbourne, Victoria, Australia 3 Department of Infection Biology, London School of Hygiene and Tropical Medicine , London, UK Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Kathryn E. Holt Jane Hawkey 1 Department of Infectious Diseases, School of Translational Medicine, Monash University , Melbourne, Victoria, Australia 2 Centre to Impact AMR, Monash University , Clayton, Victoria, Australia Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Jane Hawkey Kelly L. Wyres 1 Department of Infectious Diseases, School of Translational Medicine, Monash University , Melbourne, Victoria, Australia 2 Centre to Impact AMR, Monash University , Clayton, Victoria, Australia 3 Department of Infection Biology, London School of Hygiene and Tropical Medicine , London, UK Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Kelly L. Wyres Margaret M. C. Lam 1 Department of Infectious Diseases, School of Translational Medicine, Monash University , Melbourne, Victoria, Australia 2 Centre to Impact AMR, Monash University , Clayton, Victoria, Australia Find this author on Google Scholar Find this author on PubMed Search for this author on this site For correspondence: benjamin.vezina{at}monash.edu margaret.lam{at}monash.edu Abstract Full Text Info/History Metrics Supplementary material Preview PDF Abstract Introduction Klebsiella pneumoniae is an opportunistic pathogen which causes a wide spectrum of infections within healthcare settings and the community. Four K. pneumoniae sub-lineages defined with cgMLST/LINcodes are known to cause distinct infections of the nasal and/or upper respiratory passages: SL91 and SL10031 (also referred to as subspecies ozaenae ), SL10032 (subspecies rhinoscleromatis ) and SL82. These sub-lineages have also demonstrated reduced carbon source utilisation, which in other species has been linked with high loads of insertion sequences (IS). Methods We performed comparative genomics, analysed IS and constructed genome-scale metabolic models for available public sequences from these four sub-lineages and compared them to other sub-lineages from the wider K. pneumoniae population. Results The four focal sub-lineages displayed significantly higher IS loads (median range 88 to 120 per genome) compared to other K. pneumoniae sub-lineages (median range 12 to 73). Notably, each K. pneumoniae sub-lineage had unique IS profiles, consistent with distinct evolutionary trajectories of IS acquisition and expansion. Across sub-lineages, higher IS loads were inversely associated with the number of metabolic model genes per genome (R 2 = 0.16, p <0.001), as well as predicted aerobic substrate utilisation for phosphorus sources (R 2 = 0.39, p<0.001) as per a second-degree polynomial regression model (n = 1,664 genomes). Additionally, the four IS-dense sub-lineages displayed a combination of convergent sub-lineage-specific substrate utilisation losses including parallel loss of 3-Phospho-D-glycerate, D-Glycerate-2-phosphate, Phosphoenolpyruvate utilisation as carbon/phosphorus sources. Finally, inspection of IS insertion sites demonstrated frequent and non-destructive insertion next to transcriptional, carbohydrate and amino acid metabolism genes. Conclusions IS loads were significantly and inversely associated with metabolic substrate usage within K. pneumoniae , whereby sub-lineages that had higher numbers of IS also had reduced metabolic capacity. We hypothesise that an insertional tolerance model explains these findings, whereby IS can only insert into “metabolically-tolerable” sites for the individual cell and any impacts on metabolism are not detrimental for survival. Introduction Klebsiella pneumoniae is an opportunistic pathogen which causes a diverse range of infections within healthcare settings such as pneumonia, bloodstream, surgical site and urinary tract infections. It is also a causative agent of infections in the community, associated with a unique subset of K. pneumoniae sub-lineages distinct from those that cause healthcare-associated infections ( 1 ). As a species, K. pneumoniae exhibits a remarkable amount of genome diversity, with hundreds of unique ‘deep-branching’ sub-lineages and significant variation in its accessory genome content ( 1 , 2 ). This diversity is likely an important driver of the different lifestyles and variable virulence profiles of K. pneumoniae . Two sub-lineages of K. pneumoniae , clonal groups CG67 and CG90/91 based on 7-gene MLST and formally designated as subspecies rhinoscleromatis and ozaenae , respectively, cause rare but unique chronic infections typically associated with the nasal and/or upper respiratory passages. Whole genome sequencing has since clarified that both are divergent sub-lineages of K. pneumoniae rather than separate subspecies as previously characterised Error! Reference source not found .( 3 ). The K. pneumoniae LINcode taxonomic scheme (based on a 629-loci core genome MLST) assigns K. pneumoniae into >705 discrete sub-lineages ( 4 ). Our sub-lineages of interest are SL10032 (previously rhinoscleromatis or CG67) and SL10031/SL91 (previously ozeanae or CG90/CG91). SL10032 strains cause a chronic granulomatous disease described as rhinoscleroma, whereas SL91 and SL10031 typically cause atrophic rhinitis or ozena. There have been at least 63 cases of rhinoscleromatis ( 5 - 9 ) and 45 ozaenae ( 8 - 22 ) infections documented in modern medical literature dating back to 1944, although there appears to be archaeological evidence of rhinoscleroma as far back as 300-600 AD from Maya (modern day Guatemala) ( 23 ). While these sub-lineages are usually associated with the nasopharyngeal sites, there is strong historical proof that they are able to colonise other body sites including the urinary tract, soft tissue, blood, cerebral spinal fluid and gastrointestinal tract ( 11 , 12 , 18 , 21 , 24 , 25 ). Additionally, these sub-lineages have also been found in non-human hosts, isolated from cockroaches, cattle and chicken meat ( 13 , 21 , 22 ). Aside from their characteristic clinical presentations, these sub-lineages display a reduced metabolic capacity in biochemical tests compared to other K. pneumoniae sub-lineages ( 3 , 26 ), explaining why they were originally considered distinct sub-species. These sub-lineages, along with SL82 (ST82 or CG82 based on 7-gene MLST), have demonstrated reduced phenotypic carbon source utilisation ( 3 ). A link between this metabolic reduction and niche/pathogenic lifestyle adaptation has been posited ( 3 ). This phenomenon has also been observed in other bacterial pathogens such as Shigella species and Bordetella pertussis ( 27 ). Shigella in particular is associated with the accumulation of insertion sequences (IS) ( 28 ). Reduced carbon utilisation is also seen in SL82 ( 3 ), which along with ozaenae were unique amongst K. pneumoniae in lacking the mrkD type III fimbrial adhesin. SL82 expressed the virulence-associated K1 capsule serotype. Similar to the other sub- lineages of interest, SL82 was noted to be strongly associated with respiratory infections (8/11 BioSample IDs with sample metadata). In this study, we performed a systematic analysis on a dereplicated collection of 1,664 completed and 210 draft K. pneumoniae genomes, to investigate the prevalence and impacts of IS on metabolic capacity in these four sub-lineages of interest. We hypothesised that, similar to Shigella , IS caused genome degradation and the loss of metabolic capabilities in these sub-lineages. Methods Genome acquisition, assembly and annotation Complete K. pneumoniae genome assemblies (n=2,302) were initially downloaded from NCBI RefSeq (accessed on 04/07/2024, Table S1 ). To further expand genome numbers and improve robustness of population-level inferences, 210 dereplicated NCBI BioSample IDs of target sub-lineages were obtained from BIGSdb ( 29 ) by selectively downloading genomes which matched the 4 sub-lineages of interest and their associated Sequence Types when no LINcode was determined, including SL10032 (ST67, ST3818, ST3819), SL91 (ST91, ST381, ST3766, ST3768), SL10031 (ST91) and SL82 (ST82, ST3764). Where available, short-read sequence reads were downloaded from SRA (n=201), otherwise assembled sequences were downloaded from Genbank (n=20) (accessed on 29/10/2024, Table S1 ). Short-read sequence data were assembled with Unicycler v0.5.1 ( 30 ). Low quality genomes were removed if they had >300 graphical fragment assembly dead ends, or in absence of assembly graphs, an N50 of <65,000. This threshold was more lenient than previously defined quality control metrics ( 31 ), to account for IS fragmenting short read draft assemblies, resulting in larger contig numbers and dead end counts. Remaining genomes (n=2,108) were then dereplicated using Assembly-Dereplicator v0.1.0 ( 32 ) using the following specifications: ‘--threshold 0.0003’. This resulted in a collective dataset of 1,874 dereplicated genomes. All genomes were annotated using Bakta v1.8.1 ( 33 ). Lineage assignment and genotyping Dereplicated genomes were uploaded to Pathogenwatch ( https://pathogen.watch/ ) for genotyping; namely to confirm species and assign sub-lineage, clonal groups and LIN codes ( 4 , 34 ). A neighbour-joining tree was generated with PopPUNK v2.4.0 ( 35 ) as follows. The ‘create-db’ function was used with the following options: ‘--sketch-size 1000000 --min-k 15 -- max-k 29 --qc-filter prune’. The ‘fit-model’ function was subsequently used with the following options: ‘bgmm --ranks 1,2,3,5 --graph-weights --K 3’. Additionally, the ‘poppunk_visualise’ function was used, with the ‘--distances’ and ‘--previous-clustering’ parameters, to output a neighbour joining tree. Kaptive v3.1.0 ( 36 , 37 ) was used to identify capsule synthesis loci (KL), and genotyping was performed with Kleborate v3.1.3 ( 38 ). Metabolic model construction Metabolic models were constructed using Bactabolize v1.0.3 ( 31 ) with the Kp SC pan v2.0.1 model ( 39 ) and the --draft_model command. Growth across 1,278 conditions was simulated using the --fba command on M9 minimal media, where positive growth was defined at a biomass threshold of ≥0.0001 as previously described ( 40 ). Insertion Sequence analysis Plasmids were filtered out from the complete genome sequences using seqtk v1.3 ( 41 ), with the following command: “seq -L 4500000”. Genomes were rotated to dnaA at base pair position 1 using rotate v1.0 ( 42 ) with the following command: “-s gtgtcactttcgctttggcagcagtgtcttgcccgattgcaggatgagtt -m 5”, to facilitate comparison of chromosomal synteny. Rotated, plasmid-free genomes (i.e. the chromosome) were used for the remaining analyses. ISEScan v1.7.2.3 ( 43 ) was used to identify IS elements, which also identifies novel IS elements with divergence from known IS in the ISEScan database. For comparisons within SL82, SKA v1.0.0 ( 44 ) was used to determine single nucleotide variants (SNV) using ska fasta followed by ska distance -s 25. Within sub-lineages, specific phylogenetic minimum evolution trees were constructed using SNV distances inferred from SKA, optimized via Nearest Neighbour Interchange (NNI) with the fastme.bal command from R package ape v5.8 ( 45 ), then midpoint rooted via midpoint.root from phytools v2.3-0 ( 46 ). All-vs-all whole genome alignments were performed using MUMmer4 v4.0.0 ( 47 ) using “nucmer -p”, then coordinates extracted using “show-coords -d -l -T”. Genome locus comparisons were performed by extracting gene coordinates via slice_multi_gbk.py ( https://github.com/bananabenana/slice_multi_gbk ) and Clinker v0.0.31 ( 48 ) was used for visualisation. Statistical analysis and code R v4.4 ( 49 ) and RStudio v2024.04.2+764 ( 50 ) were used for statistical analysis and visualisation. R packages tidyverse v2.0,0 ( 51 ), ggtree v3.12.0 ( 52 ), colorspace v2.1-0 ( 53 ), ggpmisc v0.6.0 ( 54 ), ggpubr v0.6.0 ( 55 ), ggh4x v0.2.8 ( 56 ), rstatix v0.7.2 ( 57 ), gggenomes v1.0.0 ( 58 ) and patchwork v1.2.0 ( 59 ) were used. All R code can be found at Figshare (doi: 10.6084/m9.figshare.28341917). Linear regression was performed between IS loads and the ratio of nucleotides assigned to IS vs total chromosomal sequence. IS loads between sub-lineages were compared using a Kruskal-Wallis test with Holm adjustment for multiple comparisons, followed by a Dunn’s post hoc test with Holm adjustment for multiple comparisons. Results Dataset description To explore the impact of IS in K. pneumoniae , we utilised two dereplicated datasets: i) a collection of complete genomes (n=1,664) used to quantify the exact numbers of IS within chromosomal sequences and phylogenetic clusters; and ii) an expanded collective dataset (n=1,874) which included targeted inclusion of draft assemblies. Throughout the manuscript, we will explicitly state which dataset was used. We found that 26 draft assemblies failed Pathogenwatch’s quality control measures ( Table S1 ), as they had >500 contigs, despite having high levels of genome completeness (≤300 graphical fragment assembly dead ends). We suspected high contig numbers may be caused by large numbers of IS and these assemblies were kept for analysis. Significant diversity was observed in the expanded dataset (n=1,874), which was made up of 244 sub-lineages (based on cgMLST). To ensure appropriate population-level inferences, only sub-lineages with n≥3 genomes were included for further analysis, resulting in a dataset of n=1,441 genomes representing 65 sub-lineages. The majority of genomes were from sub-lineages associated with known multi-drug resistant (n=1,013 across 14 MDR SLs; 61.69%) or hypervirulent clones (n=129 from 4 hypervirulent SLs; 7.86%) as defined by Wyres et al ( 1 ) ( Table S1 ). Some sub-lineages display unusually high IS loads Plasmid-free complete genomes (n=1,664) were screened for IS, and the numbers and types (i.e. IS families) were compared across the 65 sub-lineages with n≥3 genomes ( Fig. 1 ). As expected, linear regression indicated a strong, positive correlation (R 2 = 0.89, p<0.001) between IS loads and the ratio of nucleotides assigned to IS vs total chromosomal sequence ( Fig. S1 ). Across each sub-lineage, the number of IS per genome (i.e. IS-load) varied from 10 to 131 (median 31). A total of 132 unique IS from 23 distinct IS families were detected across this dataset. Sub-lineages contained 3 to 20 distinct IS families, of which only three IS families (IS 3 , IS 21 and IS NCY ) were found in every sub-lineage. Lineage-specific IS patterns were identified ( Fig. S2, Table S2 ). Aside from five novel IS (Novel 20, 24, 189, 272 and 348, which were detected in sub-lineages SL43, SL91, SL107, SL147 and SL17, respectively. No other IS families were restricted to a single sub-lineage. An additional 19 IS were rarely detected across the population (i.e. found in <10 sub-lineages). Download figure Open in new tab Fig. 1: Sub-lineages display significantly different IS loads Distribution of Insertion Sequences (IS) within 65 K. pneumoniae sub-lineages containing n≥3 complete genomes (n=1,441). The top panel shows the IS per genome while bottom panel shows the breakdown of IS families. Sub-lineages are coloured by their relevant global clone information ( 60 ), where each point represents a single genome. ‘<50 occurrences’ refers to IS families which were not shown in this figure due to their low abundance. Significance not shown for brevity - raw data and results of Dunn’s post hoc test in Table S2 . Sub-lineages SL82, SL91-ozaenae, SL10031-ozaenae and SL10032-rhinoscleromatis carried the highest IS loads, with median 88 to 120 IS per genome compared to other SLs (median 12 to 73, Fig. 1 ). Kruskal-Wallis with Holm error correction indicated sub-lineages had significantly different IS loads (p < 0.0001). This was followed by a Dunn’s post hoc test with Holm error correction, showing SL82 had significantly more IS per genome than 37 other sub-lineages (p-values <0.05), while SL91-ozaenae, SL10031-ozaenae and SL10032-rhinoscleromatis had significantly higher IS loads than 25, 5 and 6 other sub-lineages, respectively ( Table S2 ). There were no significant differences between SL0032-rhinoscleromatis (median 120 ± 2 IQR), SL91-ozaenae (median 113 ± 21.5 IQR), SL82 (median 112 ± 7 IQR) or SL10031-ozanae (median 88 ± 1). We focused on these four IS-dense sub-lineages for the remaining analyses. The complete genome set was supplemented with additional short-read assemblies that could be assigned LINcodes, resulting in a total of n=1,874 genomes, which included n=4 SL10032-rhinoscleromatis, n=15 SL91-ozaenae, n=5 SL10031-ozaenae, and n=36 SL82. This statistical analysis also revealed that several other sub-lineages displayed significantly higher IS loads compared to others ( Table S2 ). This included genomes that corresponded to MDR (SL258, SL15, SL395), hypervirulent (SL65) and common clones (SL34, SL3010), although there appeared to be within-clone variation in many sub-lineages. For example, the number of IS per genome within sub-lineage SL258 ranged from 22 to 80. Within sub-lineage variability in IS loads was typically larger for sub-lineages represented by more genomes. With the exception of our four sub-lineages of interest, the IS loads were generally similar across other clones regardless of MDR or hypervirulence status. Aside from isolates of unknown sampling site (n=42, 70% of genomes from the four sub-lineages of interest), the remaining isolates were recovered from human nasopharyngeal sites including sputum (n=8), nasopharynx (n=6), maxillary sinuses (n=2), pleural cavity (n=1) and throat (n=1) ( Table S1 ). One additional isolate was isolated from blood (n=1). While source metadata was missing for many of the SL10032-rhinoscleromatis and SL91-ozanae isolates, sampling dates indicated collection between 1920 to 1952. Many samples collected during this period were generally of clinical origin. While there is likely a sequencing bias for clinical isolates, these sub-lineages have previously been isolated from a variety of non-human sources including chicken meat, cattle and cockroaches ( 13 , 21 , 22 ), although no such isolates are currently represented in public genome repositories. IS profiles varied within IS-dense sub-lineages The four focal IS-dense sub-lineages contained 16 (SL91-ozaenae), 15 (SL82), 13 (SL10031-ozaenae) and 10 (SL10032-rhinoscleromatis) IS families ( Fig. 2 ). Of these, seven IS families (IS 1 , IS L3 , IS 21 , IS 66 , IS 256 , IS NCY and IS 200 /IS 605 ) were detected across all four lineages ( Fig. S2 ). The IS profile of SL10032-rhinoscleromatis was the most distinct from the others while the closely-related ozaenae sub-lineages and SL82 shared highly similar IS profiles, consistent with their divergence from a more recent common ancestor. Download figure Open in new tab Fig. 2: IS family prevalence differ between K. pneumoniae sub-lineages Distribution of Insertion Sequence (IS) families across the four most IS-dense sub-lineages, aligned against a neighbour-joining tree generated from core k-mers, whereby each tip represents a unique genome and is coloured by the sample type. Only complete genomes with plasmids filtered out were used in this analysis. Most notable is the considerable expansation of IS 1 in SL82 (mean 71.53±6.31 SD copies per genome), SL91-ozaenae (70.14±7.9) and SL10031-ozaenae (51.3±3.5) genomes, while it was conversely rarer in SL10032-rhinoscleromatis; (1±0). IS 1 expansion (mean ≥10 copies) was not specific to these sub-lineage and was found in six others including hypervirulent SL25 and other clones SL10022 and SL383. There were, however, notable differences between the four focal sub-lineages and which IS families had undergone significant expansions. For example, expansion of IS 3 was observed in SL10032-rhinoscleromatis (mean of 53±1 SD copies per genome) and SL82 (12.8±3.2), while SL91-ozanae and SL10031-ozanae carried 3.7±1 or no IS 3 , respectively. IS L3 was highly expanded in SL10032-rhinoscleromatis (24±1), SL91-ozaenae (10.57±10.81), SL10031-ozaenae (9±2), but was largely absent in SL82 (1±0). In SL10032-rhinoscleromatis, IS 3 , IS L3 and IS 21 comprised the largest proportions of IS (mean 53±1, 24±1 and 22.7±0.6 SD copies per genome, respectively), while carrying fewer IS 1 (mean 1±0 copies). IS load impacts strain metabolism and redundancy We next analysed the impact of IS on metabolism across the entire collective dataset (complete and short read-assembled genomes), although IS loads of short-read assemblies were not used in correlative measures. Metabolic models were built for all sub-lineages with n≥3 isolates, resulting in 1,874 genomes across 70 sub-lineages. Of these, 9 sub-lineages required minimal gap filling to simulate growth on M9 media + D-glucose (mean 2.89±2.33 SD gap-filled reactions), with SL34 (n=4) being the most prevalent lineage, caused by poorer draft genome assemblies. This is a standard process which improves simulated growth prediction accuracy ( 31 ). Our model predictions were able to simulate 61/99 previously generated biochemical tests ( 3 ) and models are unable simulate the remaining substrates. As such, we compared metabolic modelling growth predictions to these 61 phenotypic substrates to generate accuracy statistics. As many of the biochemically-tested isolates do not have corresponding genome data, we compared model accuracy on a per-sub-lineage rather than individual basis to account for this ( 40 ). F1 scores varied, ranging from 0.72 to 0.83 ( Table 1 ). The inaccuracies were largely driven by false positives, where the models predicted growth due to presence of intact genes, while phenotypic tests predicted no growth, as seen by the high recall values ranging from 0.88 to 0.94. View this table: View inline View popup Download powerpoint Table 1: Summary of accuracy metrics of simulated metabolic model growth predictions. A total of 61 compatible phenotypic biochemical tests from ( 3 ) were used. Full data and confusion matrix can be found in Table S3 . We next examined the impact of IS load on strain metabolism. Analysis of only the completed genomes demonstrated a clear inverse relationship between IS load per genome and metabolic capacity, whereby high-IS sub-lineages displayed reduced metabolism ( Fig. 3 ). Fitting a squared polynomial to the data, the number of model genes was inversely associated with IS load (R 2 =0.16, p <0.001) and slightly less with model reactions (R 2 =0.16, p <0.001). For substrate usage, only aerobic use of phosphorus sources displayed a moderate inverse relationship (R 2 =0.39, p <0.001). The association between IS loads and reduced anaerobic usage of substrates was weaker compared to aerobic substrate usage, likely due to the reduced metabolic capacity of K. pneumoniae as a whole under anaerobic conditions ( 40 ). Download figure Open in new tab Fig. 3: Number of IS per genome is inversely related to loss of metabolic substrate usage. Number of insertion sequences (IS) per genome plotted against key metabolic model metrics, with a squared polynomial model fitted to data. Each point represents a completed genome. The metabolism of some over-represented sub-lineages such as MDR global clone SL258 (n=427) were unaffected by IS load (R 2 <0.01, p <0.263 to <0.96) despite a wide IS range (19–89) ( Fig. S3 ). To control for this over-representation effect, we calculated mean values with each sub-lineage and re-examined the relationships ( Fig. S4 ). Increasing IS load was now weakly associated with reduced aerobic use of carbon (R 2 =0.23, p <0.001), nitrogen (R 2 =0.13, p <0.001) and moderately associated with phosphorus (R 2 =0.42, p <0.001) sources using a squared polynomial model. The reduced metabolic capacity of sub-lineages with higher IS-loads was consistent with the similar inverse relationships of total number of metabolic model genes (R 2 =0.31, p<0.001) and model reactions (R 2 =0.25, p<0.001). There was no correlation between IS load and the number of pseudogenes (R 2 < 0.01, P = 0.775), which was also previously reported in Shigella species ( 28 ). A total of 232 reference model genes were missing across the complete genomes from SL82, SL91-ozaenae, SL10031-ozaenae and SL10032-rhinoscleromatis. These were a combination of core and variable genes across the dataset generally. Flanking sequences of up to 50bp around IS sites were matched against missing model genes to determine if IS had interrupted them directly. Only five genes across 16 total occurrences showed evidence of direct IS interruption. This meant that 1.5% of missing model genes can be directly attributed to a known IS insertion. For example, the loss of kpnE , an essential component of the KpnEF spermidine antiporter was absent in a single SL82 genome (accession GCF_900452625, coordinates 2591281:2597710 on contig 1), caused by insertion of IS 1 ( Fig. 4 ). Download figure Open in new tab Fig. 4: IS 1 interruption of the kpnE gene required for the KpnEF spermidine antiporter in an SL82 genome (accession GCF_900452625.1). Metabolic evolutionary parallelism Aside from their loss of total substrate usage, these sub-lineages displayed distinct but convergent metabolic trajectories. Previously, 616 core metabolic traits were found to be conserved across the 48 most prevalent K. pneumoniae sub-lineages from a dataset comprising large-scale studies (4,621 genomes; n≥15 genomes in each) ( 40 ). By raw counts, each of the genomes from the four focal sub-lineages in our analysis appear to have completely lost the ability to utilise 15 to 20 core K. pneumoniae substrate growth conditions, while intermediately losing 12 to 49 ( Fig. 5 , Table S3 ). Download figure Open in new tab Fig. 5: Sub-lineages display unique and convergent metabolic losses Venn diagrams showing number of K. pneumoniae core substrate usages intermediately (usage within >0% & <95%) or completely lost (usage within 0%) within sub-lineages. There was evidence of convergent substrate utilisation losses whereby all four sub-lineages have lost the ability to use 3-Phospho-D-glycerate, D-Glycerate-2-phosphate, Phosphoenolpyruvate as carbon or phosphorus sources. This was due to the loss of a periplasm antiporter pgtP and associated operon pgtABC , responsible for the import of these substrates from the periplasm to the cell. Due to the complete contextual loss of this operon and numerous intra-genomic rearrangements, it is not clear if this loss was IS-mediated. As expected based on their genetic relatedness, SL91-ozaenae and SL10031-ozaenae displayed 10 conserved substrate losses, likely as a result of common ancestry. Despite similar IS profiles between SL82 and the two ozaenae sub-lineages, they displayed little overlap of metabolic losses. Sub-lineage-specific complete metabolic losses included SL82’s loss of 3-hydroxyphenylacetic acid usage as a carbon source, SL10032-rhinoscleromatis’s loss of 5-Aminopentanoate, L-Lysine and 4-Aminobutanoate usage as carbon and nitrogen sources and SL91-ozaenae/SL10031-ozaenae’s loss of (R)-3-Hydroxybutanoate, Cytosine and Myo-Inositol hexakisphosphate usage as carbon and nitrogen sources ( Table S3 ). Siderophore virulence genes associated with iron uptake were common within these sub-lineages, with 29/60 genomes displaying Kleborate virulence scores ≥3 (i.e. carry yersiniabactin and aerobactin) ( 38 ) ( Table S1 ). Yersiniabactin was found in 3/4 SL10032-rhinoscleromatis genomes (lineage: ybt 11), 13/15 SL91-ozaenae (lineage: ybt 10), 2/5 SL10031-ozanae (lineage: ybt 16), although was truncated in 12/18 cases across this dataset. Aerobactin was found in all SL10032-rhinoscleromatis (lineage: iuc 4), 6/15 SL91-ozaenae (lineage: iuc 2A), 2/5 SL10031-ozaenae (lineage: iuc 2A) and 17/36 SL82 (lineage: iuc 2A), with only one truncated case in SL10032-rhinoscleromatis. Salmochelin (encoded by iro ) was found only in 17/36 SL82. Previously reported mrkD (type 3 fimbriae adhesin) ( 3 ) was completely absent in all focal sub-lineages with the exception of SL10032-rhinoscleromatis (present in 5/5 genomes). The rmpADC locus which controls capsule expression and hypermucoviscosity was also commonly detected in genomes from these four sub-lineages. It was found in 32/36 SL82 genomes (lineage: rmp 2A, truncated in 26 instances), all SL10032-rhinoscleromatis genomes (lineage: rmp 4, truncated in 1 instance), and 1/15 SL91-ozanae (lineage: rmp 2A, truncated). These virulence gene truncations were associated with the ends of assembly contigs, likely caused by IS insertions. All SL82, SL10031-ozaenae and SL10032-rhinoscleromatis isolates were typeable and displayed homogenous capsular loci (KL1, KL5 and KL3, respectively), consistent with previous reports ( 3 ) based on capsule typing. In contrast, SL91-ozaenae isolates displayed varying K loci with closest matches to KL4 (n=13), KL6 (n=1), KL1 (n=1), with all but the KL6 genome being untypeable – again due to capsule gene truncations and fragmented assemblies, likely caused by IS insertions within the capsule loci, which has been previously reported ( 63 ). As public assemblies lacked reads, it was not possible to leverage assembly graphs to determine which IS caused contig breaks – though in complete genomes such as GCF_014218685.1, IS 1 was responsible for the capsule null predictions. IS variation within sub-lineages confers strain-level metabolic diversity To study the localised impact of IS further, we analysed all available SL82 genomes (n=17) as this sub-lineage had the largest number of complete genomes. After re-orientating the genomes to the same starting position, we mapped IS across them and inspected IS insertion sites and chromosomal synteny ( Fig. 6 ). While many of these genomes within this sub-lineage shared similar IS profiles and IS density, no two genomes were identical with respect to IS insertions, highlighting the dynamic nature of IS and consequent losses or gains of DNA segments. Additionally, chromosomal inversions and rearrangements were common, which were often flanked by IS, reasonably explaining the rearrangements occurring at homologous, repeat sites. Download figure Open in new tab Fig. 6: Insertion Sequences vary locationally within SL82. Comparison of insertion sequences (IS) and homologous blocks between SL82 complete chromosomes. Tree shown is a SNV distance tree showing 17 genomes. Grey links between genomes represent homologous sequence blocks as aligned by nucmer with standard parameters. Coloured lines show IS family. All genomes start at dnaA . The genomes within SL82 were also predicted to have small but varying predicted aerobic metabolism across 24 substrates, with 13/17 genomes showing loss of usage for at least one substrate. The most variable substrates were 4-Hydroxybenzoate (usage observed in 7/17 genomes), then L-alanine-D-glutamate-meso-2,6-diaminoheptanedioate (13/17 genomes), and L-Galactonate (15/17 genomes). The remaining 21 substrates were absent only in a single genome (16/17). Discussion Reduced metabolic versatility has previously been documented in select sub-lineages of K. pneumoniae : SL10032-rhinoscleromatis, SL10031-ozaenae, SL91-ozaenae and SL82 ( 3 ). This is thought to be an important driver of their unique ecological and pathogenic profiles, particularly as these strains are often isolated from the upper respiratory tract but are by no means restricted to these sites or human host ( 13 , 21 , 22 ). In this study, application of genome-scale metabolic modelling supported significant loss of substrate usage across SL82, SL10031-ozaenae, SL91-ozaenae and SL10032-rhinoscleromatis compared to other K. pneumoniae sub-lineages. In particular, these sub-lineages appear to have lost the ability to utilise gut microbe/human metabolites for energy including 3-Phospho-D-glycerate, D-glycerate-2-phosphate, phosphoenolpyruvate, 3-hydroxyphenylacetic acid and 5-aminopentanoate ( 64 - 68 ), along with myo-Inositol hexakisphosphate found in plant tissues and human diet ( 70 ) and human metabolites (R)-3-hydroxybutanoate (liver) ( 71 ) and 4-Aminobutanoate (neurotransmitter) ( 72 ). These hint that these sub-lineages have deviated from other K. pneumoniae sub-lineages that readily colonise the gut prior to infection ( 69 ). We demonstrate for the first time that these four sub-lineages also have significantly higher IS loads ( Fig. 1 ), which has been linked to streamlining of metabolic profiles in Bordetella pertussis and various Shigella species ( 27 , 28 ). IS has been previously quantified within K. pneumoniae SL258 (ST258) ( 73 ), which found a stable repertoire of IS elements and small within-lineage variation (specifically, an increase in IS load in one subclade of SL258, referred to as ST258B). This is consistent with the IS profiles and within-sub-lineage variation that were observed in our study ( Fig. 1 , Fig. 2 , Table S2 ), such as intermediate rather than complete loss of core K. pneumoniae substrates within a sub-lineage ( Fig. 5 ). Notably, the lack of impact IS had on the metabolism of SL258 differed from that of the four focal sub-lineages. This indicated that K. pneumoniae sub-lineages have different relationships and IS-mediated evolutionary trajectories. This analysis was limited by the sampling biases of public data, and the relative rarity of these particular sub-lineages. Direct IS interruptions of metabolic genes were detected in only a few instances (only 5/323 examples were found). We hypothesize that IS may have been responsible for the initial disruption of other key metabolic genes, whose loss was selectively tolerated at the time of transposition, leading to subsequent genetic degradation of the associated metabolic loci over time. The inverse relationship between IS loads and substrate usage, primarily aerobic carbon and phosphorus usage ( Fig. 3 , Fig. S4 ) is consistent with this hypothesis. Insertion of an IS may have caused the initial interruption of the gene but have been subsequently purged from the genome over time -either via recombination, rearrangement or genome degradation. Additionally, higher IS loads potentiate a higher rate of genomic rearrangement ( 74 ), deleting genes without direct insertion. This highlighted that DNA sequence analysis of isolated pathogens cannot capture the in vivo conditions or genetic history which shaped the organism even when historical isolates are used. Hawkey et al . noted significant genome degradation from IS insertions in Shigella species which caused metabolic usage losses ( 28 ). We did not observe any correlation between IS load and genome length in our dataset, contrasting to Shigella where the two appeared to be related. The metabolic gene losses in both Shigella and the four K. pneumoniae sub-lineages of interest provide insight into the cellular processes required for these pathogens to colonise and progress to invasive human disease. Our data show that while these sub-lineages displayed loss of some K. pneumoniae -specific core traits, they retained use of 552 to 587 core metabolic traits. Given conservation of these core metabolic traits across all K. pneumoniae sub-lineages, these are likely key to the species’ ability to colonise and/or cause infections. As it currently stands, the metabolic models have broad agreement with previously published biochemical tests ( 3 , 26 ) ( Table 1 ) whereby these IS-dense sub-lineages exhibited considerable loss of metabolic potential ( Fig. 5 ). The decaying, polynomial relationship between increased number of IS per genome and loss of substrate growth conditions, as well as metabolic genes and reactions ( Fig. 3 ), demonstrates this effect. Lower accuracy scores were observed, ( Table 1 ) and can be explained by various reasons. False positives are not always indicative of model inaccuracies, but usually indicate the presence of metabolic genes which failed to be expressed in the experimental condition due to gene regulation ( 40 , 75 ), or phenotypically-tested isolates do not have intact copies of key genes. Discrepancies likely stem from strain choice used in biochemical tests, which may not be reflective of the larger population. In some cases, false-positives could be explained by IS interrupting the promoter/regulatory regions required for this expression to occur. For example, the genes required for L-histidine, L-tyrosine, ethanolamine, quinate, 2-Ketoglutarate and putrescine utilisation are all present in the genomes but phenotypically showed no growth. Alternatively, if these false positives did not arise from gene regulation issues, there may be incorrect assumptions being made in the metabolic models, such as over-assignment of metabolic genes. A metabolic model gene is considered ‘present’ if there is ≥25% bi-directional coverage (standard value) ( 31 , 39 , 76 ). This would indicate our model-based analysis is likely underrepresenting the impact of IS reducing the metabolic capacity of isolates. It is also possible that the four focal IS-dense sub-lineages may contain additional novel or specialised metabolism that are not accounted for in the current model. Insertion of IS not only leads to genomic rearrangements as had been observed for SL82 ( Fig. 6 ) but have also been shown to impact expression of neighbouring genes via their promoters. One prime example is IS 1 ( 61 , 62 ). Considering the vast numbers of IS found across these IS-dense sub-lineages and their close proximity to many regulatory, amino acid and carbohydrate metabolism loci ( Fig. S2 ), it may be that IS are facilitating metabolic specialisation via provision of a selective advantage for organisms by transcriptionally upregulating remaining metabolic traits in addition to purifying unneeded metabolic traits. For example, a copy of IS 1 was found directly upstream of the cra catabolite repressor/activator gene in an SL82 genome (accession GCF_900451215) that regulates the acetolactate ilvIN locus responsible for amino acid biosynthesis from pyruvate. Increasing expression of cra may provide a selective advantage for SL82 and balance the effects arising from usage loss of 17 core substrates. The proximity of IS to metabolic genes could be interpreted as: i) IS ferrying these associated metabolic genes with them during insertion; or ii) IS can only insert into metabolically-tolerable sites for the individual cell, where only metabolism that is not essential for survival can be selectively purified. This insertional tolerance model, combined with potential upregulation of these operons is a compelling hypothesis, as this would improve the fitness of an individual strain via metabolic specialisation. This is consistent with the limited examples of directly interrupted, vestigial metabolic loci observed in our study. Future work could entail using transcriptional data to study the impact of insertions within intergenic regions ( 62 ). Competing interests The author(s) declare no competing interests. Data availability All data used in this study is available as supplemental material (Figs. S1-4, Tables S1-4). Additionally, all analysis code is available at Figshare (doi: 10.6084/m9.figshare.28341917). Author contributions Conceptualization: KEH, KLW, JH, MMCL Data Curation: BV, CW Formal Analysis: BV, CW, HBC, MMCL Funding Acquisition: MMCL Methodology: BV Project Administration: KEH, KLW, JH, MMCL Resources: KEH, KLW, JH, MMCL Supervision: KEH, KLW, JH, MMCL Writing – Original Draft Preparation: BV, KLW, MMCL All co-authors reviewed and approved the submitted manuscript. Funding MMCL is supported by an Australian National Health and Medical Research Council Investigator Grant [APP2009163]. JEH is supported by an Australian National Health and Medical Research Council Investigator Grant [APP2034741]. Supplemental material Table S1: Metadata and strain information of all genomes used in this study. Includes BioSample and accession numbers. Table S2: Results of Insertion Sequence analysis on genomes and statistical analyses Table S3: Results of substrate usage predictions using metabolic models and comparisons to phenotypic tests Fig. S1: Scatterplot of number of IS per genome compared to the IS bp:chromosomal bp ratio per genome. Linear regression model fitted to data. Fig. S2: IS families present in each sub-lineage. Fig. S3: Scatterplot comparing IS load and metabolic model features of SL258. Squared polynomial model fitted. Fig. S4: Scatterplot comparing IS load and metabolic model features of all sub-lineages. Squared polynomial model fitted. Acknowledgements This research/work was supported by Monash eResearch capabilities, including M3 and Research Data Storage. We thank the Institut Pasteur teams for the curation and maintenance of BIGSdb-Pasteur databases at http://bigsdb.pasteur.fr/ Funder Information Declared National Health and Medical Research Council, https://ror.org/011kf5r70 , APP2009163 , APP2034741 References 1. ↵ Wyres KL , Lam MMC , Holt KE . Population genomics of Klebsiella pneumoniae . Nature Reviews Microbiology . 2020 ; 18 ( 6 ): 344 – 59 . OpenUrl CrossRef PubMed 2. ↵ Holt KE , Wertheim H , Zadoks RN , Baker S , Whitehouse CA , Dance D , et al. Genomic analysis of diversity, population structure, virulence, and antimicrobial resistance Klebsiella pneumoniae, an urgent threat to public health . Proceedings of the National Academy of Sciences . 2015 ; 112 ( 27 ): E3574 . OpenUrl Abstract / FREE Full Text 3. ↵ Brisse S , Fevre C , Passet V , Issenhuth-Jeanjean S , Tournebize R , Diancourt L , et al. Virulent Clones of Klebsiella pneumoniae: Identification and Evolutionary Scenario Based on Genomic and Phenotypic Characterization . PLOS ONE . 2009 ; 4 ( 3 ): e4982 . OpenUrl CrossRef PubMed 4. ↵ Hennart M , Guglielmini J , Bridel S , Maiden MCJ , Jolley KA , Criscuolo A , et al. A Dual Barcoding Approach to Bacterial Strain Nomenclature: Genomic Taxonomy of Klebsiella pneumoniae Strains . Molecular Biology and Evolution . 2022 ; 39 ( 7 ): msac135 . OpenUrl CrossRef PubMed 5. ↵ Merbouh M , El Aidouni G , Housni B. Invasive Klebsiella rhinoscleromatis infection leading to ARDS, septic shock, and death: A rare case report . Radiol Case Rep . 2023 ; 18 ( 8 ): 2574 – 6 . OpenUrl PubMed 6. Andraca R , Edson RS , Kern EB . Rhinoscleroma: A Growing Concern in the United States? Mayo Clinic Experience . Mayo Clinic Proceedings . 1993 ; 68 ( 12 ): 1151 – 7 . OpenUrl CrossRef PubMed Web of Science 7. Kumade E , Furusyo N , Takeshima N , Kishihara Y , Mitsumoto-Kaseida F , Etoh Y , et al. A case of lobar pneumonia and sepsis with death caused by invasive Klebsiella rhinoscleromatis infection . Journal of Infection and Chemotherapy . 2016 ; 22 ( 10 ): 707 – 11 . OpenUrl PubMed 8. ↵ Malowany MS , Chester B , Allerhand J. Isolation and Microbiologic Differentiation of Klebsiella rhinoscleromatis and Klebsiella ozaenae in Cases of Chronic Rhinitis . American Journal of Clinical Pathology . 1972 ; 58 ( 5 ): 550 – 3 . OpenUrl CrossRef PubMed 9. ↵ Berger SA , Pollock AA , Richmond AS . Isolation of Klebsiella ozaenae and Klebsiella rhinoscleromatis in a General Hospital . American Journal of Clinical Pathology . 1977 ; 67 ( 5 ): 499 – 502 . OpenUrl CrossRef PubMed 10. De Champs C , Vellin JF , Diancourt L , Brisse S , Kemeny JL , Gilain L , et al. Laryngeal Scleroma Associated with Klebsiella pneumoniae subsp . ozaenae. Journal of Clinical Microbiology . 2005 ; 43 ( 11 ): 5811 – 3 . OpenUrl PubMed 11. ↵ Tang LM , Chen ST . Klebsiella ozaenae meningitis: report of two cases and review of the literature . Infection . 1994 ; 22 ( 1 ): 58 – 61 . OpenUrl CrossRef PubMed 12. ↵ Goldstein EJ , Lewis RP , Martin WJ , Edelstein PH . Infections caused by Klebsiella ozaenae: a changing disease spectrum . J Clin Microbiol . 1978 ; 8 ( 4 ): 413 – 8 . OpenUrl Abstract / FREE Full Text 13. ↵ Stefani S , Giovanelli I , Anacarso I , Condò C , Messi P , de Niederhäusern S , et al. Prevalence and characterization of extended-spectrum β-lactamase-producing Enterobacteriaceae in food-producing animals in Northern Italy . New Microbiol . 2014 ; 37 ( 4 ): 551 – 5 . OpenUrl PubMed 14. AlJindan R. Epidemiological characteristics of Klebsiella ozaenae infection and its antibiotic susceptibility: Experience of a tertiary care hospital in the Eastern Province of Saudi Arabia . J Family Community Med . 2024 ; 31 ( 2 ): 148 – 52 . OpenUrl PubMed 15. Helmy AK , Sidkey NM , El-Badawy RE , Hegazi AG . Emergence of microbial infections in some hospitals of Cairo, Egypt: studying their corresponding antimicrobial resistance profiles . BMC Infectious Diseases . 2023 ; 23 ( 1 ): 424 . OpenUrl PubMed 16. Livermore DM , Yuan M. Antibiotic resistance and production of extended-spectrum beta-lactamases amongst Klebsiella spp . from intensive care units in Europe. J Antimicrob Chemother . 1996 ; 38 ( 3 ): 409 – 24 . OpenUrl PubMed 17. Korvick JA , Bryan CS , Farber B , Beam TR , Schenfeld L , Muder RR , et al. Prospective observational study of Klebsiella bacteremia in 230 patients: outcome for antibiotic combinations versus monotherapy . Antimicrobial Agents and Chemotherapy . 1992 ; 36 ( 12 ): 2639 – 44 . OpenUrl Abstract / FREE Full Text 18. ↵ Tula A , Mikru A , Alemayehu T , Dobo B. Bacterial Profile and Antibiotic Susceptibility Pattern of Urinary Tract Infection among Pregnant Women Attending Antenatal Care at a Tertiary Care Hospital in Southern Ethiopia . Can J Infect Dis Med Microbiol . 2020 ; 2020 : 5321276 . OpenUrl PubMed 19. Tadesse B , Shimelis T , Worku M. Bacterial profile and antibacterial susceptibility of otitis media among pediatric patients in Hawassa, Southern Ethiopia: cross-sectional study . BMC Pediatrics . 2019 ; 19 ( 1 ): 398 . OpenUrl PubMed 20. Regassa BT , Tosisa W , Eshetu D , Beyene D , Abdeta A , Negeri AA , et al. Antimicrobial resistance profiles of bacterial isolates from clinical specimens referred to Ethiopian Public Health Institute: analysis of 5-year data . BMC Infectious Diseases . 2023 ; 23 ( 1 ): 798 . OpenUrl PubMed 21. ↵ Fielding BC , Mnabisa A , Gouws PA , Morris T. Basic researchAntimicrobial-resistant Klebsiella species isolated from free-range chicken samples in an informal settlement . Archives of Medical Science . 2012 ; 8 ( 1 ): 39 – 42 . OpenUrl PubMed 22. ↵ Oothuman P , Jeffery J , Aziz AHA , Bakar EA , Jegathesan M. Bacterial pathogens isolated from cockroaches trapped from paediatric wards in peninsular Malaysia . Transactions of The Royal Society of Tropical Medicine and Hygiene . 1989 ; 83 ( 1 ): 133 – 5 . OpenUrl CrossRef PubMed 23. ↵ Goldman L. Pre-Columbian Rhinoscleroma . Archives of Dermatology . 1979 ; 115 ( 1 ): 106 – 7 . OpenUrl PubMed 24. Jeffery J , Sulaiman S , Oothuman P , Vellayan S , Zainol-Ariffin P , Paramaswaran S , et al. Domiciliary cockroaches found in restaurants in five zones of Kuala Lumpur Federal Territory, peninsular Malaysia . Trop Biomed . 2012 ; 29 ( 1 ): 180 – 6 . OpenUrl PubMed 25. ↵ Gaafar HA , Gaafar AH , Nour YA . Rhinoscleroma: An updated experience through the last 10 years . Acta Oto-Laryngologica . 2011 ; 131 ( 4 ): 440 – 6 . OpenUrl CrossRef PubMed 26. ↵ Imhoff J. Bergey’s Manual® of Systematic Bacteriology . 2005 . p. 587 – 850 . 27. ↵ Belcher T , Dubois V , Rivera-Millot A , Locht C , Jacob-Dubuisson F. Pathogenicity and virulence of Bordetella pertussis and its adaptation to its strictly human host . Virulence . 2021 ; 12 ( 1 ): 2608 – 32 . OpenUrl CrossRef PubMed 28. ↵ Hawkey J , Monk JM , Billman-Jacobe H , Palsson B , Holt KE . Impact of insertion sequences on convergent evolution of Shigella species . PLOS Genetics . 2020 ; 16 ( 7 ): e1008931 . OpenUrl 29. ↵ Jolley KA , Maiden MCJ . BIGSdb: Scalable analysis of bacterial genome variation at the population level . BMC Bioinformatics . 2010 ; 11 ( 1 ): 595 . OpenUrl CrossRef PubMed 30. ↵ Wick RR , Judd LM , Gorrie CL , Holt KE . Unicycler: Resolving bacterial genome assemblies from short and long sequencing reads . PLOS Computational Biology . 2017 ; 13 ( 6 ): e1005595 . OpenUrl CrossRef 31. ↵ Vezina B , Watts SC , Hawkey J , Cooper HB , Judd LM , Jenney AWJ , et al. Bactabolize is a tool for high-throughput generation of bacterial strain-specific metabolic models . eLife . 2023 ; 12 : RP87406 . OpenUrl CrossRef PubMed 32. ↵ Wick RR , Holt KE . Assembly Dereplicator . 0.1.0 ed: https://github.com/rrwick/Assembly-Dereplicator ; 2019 . 33. ↵ Schwengers O , Jelonek L , Dieckmann MA , Beyvers S , Blom J , Goesmann A. Bakta: rapid and standardized annotation of bacterial genomes via alignment-free sequence identification . Microbial Genomics . 2021 ; 7 ( 11 ). 34. ↵ Argimón S , David S , Underwood A , Abrudan M , Wheeler NE , Kekre M , et al. Rapid Genomic Characterization and Global Surveillance of Klebsiella Using Pathogenwatch . Clin Infect Dis . 2021 ; 73 ( Suppl_4 ): S325 – s35 . OpenUrl CrossRef PubMed 35. ↵ Lees JA , Harris SR , Tonkin-Hill G , Gladstone RA , Lo SW , Weiser JN , et al. Fast and flexible bacterial genomic epidemiology with PopPUNK . Genome Research . 2019 ; 29 ( 2 ): 304 – 16 . OpenUrl Abstract / FREE Full Text 36. ↵ Lam MMC , Wick RR , Judd LM , Holt KE , Wyres KL . Kaptive 2.0: updated capsule and lipopolysaccharide locus typing for the Klebsiella pneumoniae species complex . Microbial Genomics . 2022 ; 8 ( 3 ). 37. ↵ Wyres KL , Wick RR , Gorrie C , Jenney A , Follador R , Thomson NR , et al. Identification of Klebsiella capsule synthesis loci from whole genome data . Microbial Genomics . 2016 ; 2 ( 12 ). 38. ↵ Lam MMC , Wick RR , Watts SC , Cerdeira LT , Wyres KL , Holt KE . A genomic surveillance framework and genotyping tool for Klebsiella pneumoniae and its related species complex . Nature Communications . 2021 ; 12 ( 1 ): 4188 . OpenUrl PubMed 39. ↵ Cooper HB , Vezina B , Hawkey J , Passet V , López-Fernández S , Monk JM , et al. A validated pangenome-scale metabolic model for the Klebsiella pneumoniae species complex . Microbial Genomics . 2024 ; 10 ( 2 ). 40. ↵ Vezina B , Cooper HB , Rethoret-Pasty M , Brisse S , Monk JM , Holt KE , et al. A metabolic atlas of the Klebsiella pneumoniae species complex reveals lineage-specific metabolism that supports persistent co-existence of diverse lineages . bioRxiv . 2024 :2024.07.24.605038. 41. ↵ Li H. seqtk . 2018 . 42. ↵ Durbin R , De Sanctis B , Blumer M. Rotate: A command-line program to rotate circular DNA sequences to start at a given position or string . Wellcome Open Res . 2023 ; 8 : 401 . OpenUrl PubMed 43. ↵ Xie Z , Tang H. ISEScan: automated identification of insertion sequence elements in prokaryotic genomes . Bioinformatics . 2017 ; 33 ( 21 ): 3340 – 7 . OpenUrl CrossRef PubMed 44. ↵ Harris SR . SKA: Split Kmer Analysis Toolkit for Bacterial Genomic Epidemiology . bioRxiv . 2018 : 453142 . 45. ↵ Paradis E , Claude J , Strimmer K. APE: Analyses of Phylogenetics and Evolution in R language . Bioinformatics . 2004 ; 20 ( 2 ): 289 – 90 . OpenUrl CrossRef PubMed Web of Science 46. ↵ Revell LJ . phytools 2.0: an updated R ecosystem for phylogenetic comparative methods (and other things) . PeerJ . 2024 ; 12 : e16505 . OpenUrl CrossRef PubMed 47. ↵ Marçais G , Delcher AL , Phillippy AM , Coston R , Salzberg SL , Zimin A. MUMmer4: A fast and versatile genome alignment system . PLOS Computational Biology . 2018 ; 14 ( 1 ): e1005944 . OpenUrl CrossRef 48. ↵ Gilchrist CLM , Chooi Y-H. clinker & clustermap.js: automatic generation of gene cluster comparison figures . Bioinformatics . 2021 ; 37 ( 16 ): 2473 – 5 . OpenUrl CrossRef PubMed 49. ↵ R-Core-Team . R: A Language and Environment for Statistical Computing . R Foundation for Statistical Computing ; 2020 . 50. ↵ RStudio-Team . RStudio: Integrated Development for R. Boston, MA RStudio ; 2020 . 51. ↵ Wickham H , Averick M , Bryan J , Chang W , McGowan LDA , François R , et al. Welcome to the Tidyverse . Journal of Open Source Software . 2019 ; 4 ( 43 ): 1686 . OpenUrl CrossRef 52. ↵ Yu G , Smith DK , Zhu H , Guan Y , Lam TT-Y. ggtree: an r package for visualization and annotation of phylogenetic trees with their covariates and other associated data . Methods in Ecology and Evolution . 2017 ; 8 ( 1 ): 28 – 36 . OpenUrl CrossRef 53. ↵ Zeileis A , Fisher JC , Hornik K , Ihaka R , McWhite CD , Murrell P , et al. colorspace: A Toolbox for Manipulating and Assessing Colors and Palettes . Journal of Statistical Software . 2020 ; 96 ( 1 ): 1 – 49 . OpenUrl 54. ↵ Aphalo PJ , Slowikowski K , Mouksassi S. ggpmisc: Miscellaneous Extensions to ‘ggplot2’. 0.4.1 ed 2021 . 55. ↵ Kassambara A. ggpubr: ‘ggplot2’ Based Publication Ready Plots. 0.5.0 ed 2022 . 56. ↵ Brand Tvd . ggh4x: Hacks for ‘ggplot2’ . 2024 . 57. ↵ Kassambara A. rstatix: Pipe-Friendly Framework for Basic Statistical Tests . 2023 . 58. ↵ Hackl T , Ankenbrand MJ , Adrichem Bv . gggenomes: A Grammar of Graphics for Comparative Genomics. 1.0.0 ed 2024 . 59. ↵ Pedersen TL . patchwork: The Composer of Plots. 1.1.1 ed 2020 . 60. ↵ Wyres KL , Wick RR , Judd LM , Froumine R , Tokolyi A , Gorrie CL , et al. Distinct evolutionary dynamics of horizontal gene transfer in drug resistant and virulent clones of Klebsiella pneumoniae . PLOS Genetics . 2019 ; 15 ( 4 ): e1008114 . OpenUrl PubMed 61. ↵ Amman F , D’Halluin A , Antoine R , Huot L , Bibova I , Keidel K , et al. Primary transcriptome analysis reveals importance of IS elements for the shaping of the transcriptional landscape of Bordetella pertussis . RNA Biology . 2018 ; 15 ( 7 ): 967 – 75 . OpenUrl PubMed 62. ↵ Olliver A , Vallé M , Chaslus-Dancla E , Cloeckaert A. Overexpression of the multidrug efflux operon acrEF by insertional activation with IS1 or IS10 elements in Salmonella enterica serovar typhimurium DT204 acrB mutants selected with fluoroquinolones . Antimicrob Agents Chemother . 2005 ; 49 ( 1 ): 289 – 301 . OpenUrl Abstract / FREE Full Text 63. ↵ Wei D-W , Song Y , Li Y , Zhang G , Chen Q , Wu L , et al. Insertion sequences accelerate genomic convergence of multidrug resistance and hypervirulence in Klebsiella pneumoniae via capsular phase variation . Genome Medicine . 2025 ; 17 ( 1 ): 45 . OpenUrl PubMed 64. ↵ Ho PC , Bihuniak JD , Macintyre AN , Staron M , Liu X , Amezquita R , et al. Phosphoenolpyruvate Is a Metabolic Checkpoint of Anti-tumor T Cell Responses . Cell . 2015 ; 162 ( 6 ): 1217 – 28 . OpenUrl CrossRef PubMed 65. Zhou N , Fan Y , Wang X , Wang J , Wu H. Acute enteric-coated sodium bicarbonate has negligible effect on anaerobic performance but affects metabolomics and attenuates the gastrointestinal response . Frontiers in Physiology . 2022 ;Volume 13 - 2022. 66. Zhu B , Chen X , Zhang T , Zhang Q , Fu K , Hua J , et al. Interactions between intestinal microbiota and metabolites in zebrafish larvae exposed to polystyrene nanoplastics: Implications for intestinal health and glycolipid metabolism . J Hazard Mater . 2024 ; 472 : 134478 . OpenUrl PubMed 67. Jin Z , Yang Y , Cao Y , Wen Q , Xi Y , Cheng J , et al. The gut metabolite 3hydroxyphenylacetic acid rejuvenates spermatogenic dysfunction in aged mice through GPX4-mediated ferroptosis . Microbiome . 2023 ; 11 ( 1 ): 212 . OpenUrl PubMed 68. ↵ Lee Y , Khan A , Hong S , Jee SH , Park YH . A metabolomic study on high-risk stroke patients determines low levels of serum lysine metabolites: a retrospective cohort study . Molecular BioSystems . 2017 ; 13 ( 6 ): 1109 – 20 . OpenUrl PubMed 69. ↵ Gorrie CL , Mirčeta M , Wick RR , Edwards DJ , Thomson NR , Strugnell RA , et al. Gastrointestinal Carriage Is a Major Reservoir of Klebsiella pneumoniae Infection in Intensive Care Patients . Clinical Infectious Diseases . 2017 ; 65 ( 2 ): 208 – 15 . OpenUrl CrossRef PubMed 70. ↵ Raboy V. myo-Inositol-1,2,3,4,5,6-hexakisphosphate . Phytochemistry . 2003 ; 64 ( 6 ): 1033 – 43 . OpenUrl CrossRef PubMed Web of Science 71. ↵ Clarke K , Tchabanenko K , Pawlosky R , Carter E , Todd King M , Musa-Veloso K , et al. Kinetics, safety and tolerability of (R)-3-hydroxybutyl (R)-3-hydroxybutyrate in healthy adult subjects . Regul Toxicol Pharmacol . 2012 ; 63 ( 3 ): 401 – 8 . OpenUrl CrossRef PubMed 72. ↵ Hoffman Susan S , Liang D , Hood Robert B , Tan Y , Terrell Metrecia L , Marder ME , et al. Assessing Metabolic Differences Associated with Exposure to Polybrominated Biphenyl and Polychlorinated Biphenyls in the Michigan PBB Registry . Environmental Health Perspectives . 131 ( 10 ): 107005 . 73. ↵ Adams MD , Bishop B , Wright MS . Quantitative assessment of insertion sequence impact on bacterial genome architecture . Microbial Genomics . 2016 ; 2 ( 7 ). 74. ↵ Siguier P , Gourbeyre E , Varani A , Ton-Hoang B , Chandler M. Everyman’s Guide to Bacterial Insertion Sequences . Microbiol Spectr . 2015 ; 3 ( 2 ): Mdna3 - 0030 -2014. OpenUrl 75. ↵ Ibarra RU , Edwards JS , Palsson BO . Escherichia coli K-12 undergoes adaptive evolution to achieve in silico predicted optimal growth . Nature . 2002 ; 420 ( 6912 ): 186 – 9 . OpenUrl CrossRef PubMed Web of Science 76. ↵ Norsigian CJ , Fang X , Seif Y , Monk JM , Palsson BO . A workflow for generating multi-strain genome-scale metabolic models of prokaryotes . Nature Protocols . 2020 ; 15 ( 1 ): 1 – 14 . OpenUrl PubMed View the discussion thread. Back to top Previous Next Posted July 27, 2025. Download PDF Supplementary Material Email Thank you for your interest in spreading the word about bioRxiv. NOTE: Your email address is requested solely to identify you as the sender of this article. Your Email * Your Name * Send To * Enter multiple addresses on separate lines or separate them with commas. You are going to email the following IS-capades of Klebsiella pneumoniae: Insertion sequences drive metabolic loss in obscure sub-lineages Message Subject (Your Name) has forwarded a page to you from bioRxiv Message Body (Your Name) thought you would like to see this page from the bioRxiv website. Your Personal Message CAPTCHA This question is for testing whether or not you are a human visitor and to prevent automated spam submissions. Share IS-capades of Klebsiella pneumoniae : Insertion sequences drive metabolic loss in obscure sub-lineages Ben Vezina , Claire White , Helena B. Cooper , Kathryn E. Holt , Jane Hawkey , Kelly L. Wyres , Margaret M. C. Lam bioRxiv 2025.07.24.666535; doi: https://doi.org/10.1101/2025.07.24.666535 Share This Article: Copy Citation Tools IS-capades of Klebsiella pneumoniae : Insertion sequences drive metabolic loss in obscure sub-lineages Ben Vezina , Claire White , Helena B. Cooper , Kathryn E. Holt , Jane Hawkey , Kelly L. Wyres , Margaret M. C. Lam bioRxiv 2025.07.24.666535; doi: https://doi.org/10.1101/2025.07.24.666535 Citation Manager Formats BibTeX Bookends EasyBib EndNote (tagged) EndNote 8 (xml) Medlars Mendeley Papers RefWorks Tagged Ref Manager RIS Zotero Tweet Widget Facebook Like Google Plus One Subject Area Genomics Subject Areas All Articles Animal Behavior and Cognition (7617) Biochemistry (17633) Bioengineering (13856) Bioinformatics (41841) Biophysics (21399) Cancer Biology (18529) Cell Biology (25422) Clinical Trials (138) Developmental Biology (13352) Ecology (19860) Epidemiology (2067) Evolutionary Biology (24281) Genetics (15582) Genomics (22461) Immunology (17700) Microbiology (40295) Molecular Biology (17140) Neuroscience (88413) Paleontology (666) Pathology (2823) Pharmacology and Toxicology (4813) Physiology (7632) Plant Biology (15107) Scientific Communication and Education (2042) Synthetic Biology (4284) Systems Biology (9808) Zoology (2267)