P-202 · IDWeek 2026 · supplementary handout
This handout carries the material that did not fit on the poster: the itemised corrections to the published abstract, the node support values and gene tables in full, the checks that support the re-analysis, and the isolation source of every comparison genome as originally recorded.
A subset of the bloodstream isolates had been analysed with an incorrect input file, and two datasets had been mixed. All 218 genomes were re-analysed in a single uniform pipeline. The direction of every principal finding is unchanged.
| Item | Abstract | Re-analysis |
|---|---|---|
| Public genomes | 192 | 194 |
| Human, non-bloodstream, n | 31 † | 29 |
| Bloodstream, n | 24 / 25 | 25 |
| Core genes | 538 | 1,215 |
| Total genes | 30,454 | 23,512 |
| nagK | 100% | 92% |
| nagI | 96% | 100% |
| nanJ | 96% | 100% |
| Statistics | p, uncorrected | q, BH–FDR |
| nagL | 96% vs 45% | withdrawn — 56% (14/25) vs 45% (13/29) q = 0.958 |
| pfoA | 92% vs 52% | withdrawn — 80% (20/25) vs 52% (15/29) q = 0.102 |
| tetA(P) | 8% vs 61% | withdrawn — 56% (14/25) vs 66% (19/29) OR 0.68, p = 0.58 |
Denominators are those of the frozen re-analysis, as on the poster: bloodstream 25, human non-bloodstream 29. nagL and pfoA carry a Benjamini–Hochberg q across the 18 virulence genes tested; tetA(P) is a resistance gene outside that family, so its Fisher p is uncorrected. † Used in the abstract's calculations but not printed there.
The abstract's tetA(P) figure of 8% was a transcribed odds ratio of 0.086. Its direction is not reversed — carriage is lower in the bloodstream group in both — but the effect collapses to an odds ratio of 0.68 and is not significant.
IDWeek permits errors in an accepted abstract to be indicated during the presentation (Abstract Changes / Edits). The abstract is frozen as submitted; this table is the correction.
The core-genome definition is not a knife edge. Relaxing the threshold from 100% to 80% of genomes moves the count from 822 to 2,057, and the 99% value used throughout the poster sits in the middle of that range.
| Present in ≥ | 100% | 99% | 98% | 95% | 90% | 80% |
|---|---|---|---|---|---|---|
| Core genes recovered | 822 | 1,215 | 1,376 | 1,563 | 1,685 | 2,057 |
Computed from the Roary gene_presence_absence.Rtab for the 210-genome set (23,512 genes total). The 99% column is the value reported on the poster.
The poster prints Figure 1 as a cladogram without support values, so they are given here. The tree has 210 tips and 205 internal nodes carrying support.
| Support | Median | Mean | Minimum | Nodes at or above threshold |
|---|---|---|---|---|
| SH-aLRT | 100.0 | 84.5 | 0 | 168 / 205 at ≥ 80 |
| UFBoot | 100.0 | 90.7 | 6 | 161 / 205 at ≥ 95 |
154 of 205 internal nodes (75%) meet both conventional thresholds (SH-aLRT ≥ 80 and UFBoot ≥ 95). Branch support is therefore strong for three quarters of the tree; the poorly supported nodes lie within the densely sampled, shallow parts of phylogroup III, which is consistent with the reticulation that Figure 2 shows directly.
The poster states that the tree was built from 56,596 SNP sites in a 1,066,463 bp core alignment. IQ-TREE reports a total of 1,020,921 sites, because constant sites were supplied as counts of unambiguous A/C/G/T only. The difference is accounted for here so the two figures can be reconciled.
| Component | Sites | |
|---|---|---|
| Roary core gene alignment | 1,066,463 | 210 genomes |
| SNP sites passed to IQ-TREE | 56,596 | snp-sites output |
| Constant sites supplied | 964,325 | as -fconst, so branch lengths are not inflated |
| Total modelled | 1,020,921 | as reported by IQ-TREE |
| Gap or ambiguous columns | 45,542 | neither SNP nor unambiguous constant |
All 18 genes tested, not only the seven that reached significance. Genes detected in neither group are not listed and were not tested.
| Gene | Bloodstream (25) | Human non-blood (29) | OR | p | q |
|---|---|---|---|---|---|
| nagI | 100% 25/25 | 59% 17/29 | 36.43 | 0.000166 | 0.00245 |
| nagH | 100% 25/25 | 62% 18/29 | 31.70 | 0.000408 | 0.00245 |
| nanJ | 100% 25/25 | 62% 18/29 | 31.70 | 0.000408 | 0.00245 |
| nagJ | 100% 25/25 | 66% 19/29 | 27.46 | 0.000974 | 0.00438 |
| nanI | 100% 25/25 | 69% 20/29 | 23.63 | 0.00227 | 0.00816 |
| cpe | 4% 1/25 | 38% 11/29 | 0.10 | 0.00306 | 0.00919 |
| nagK | 92% 23/25 | 59% 17/29 | 6.71 | 0.00604 | 0.0155 |
| pfoA | 80% 20/25 | 52% 15/29 | 3.49 | 0.0455 | 0.102 |
| iap | 0% 0/25 | 7% 2/29 | 0.22 | 0.493 | 0.888 |
| ibp | 0% 0/25 | 7% 2/29 | 0.22 | 0.493 | 0.888 |
| nagL | 56% 14/25 | 45% 13/29 | 1.54 | 0.586 | 0.958 |
| cloSI | 100% 25/25 | 100% 29/29 | 0.86 | 1 | 1 |
| colA | 96% 24/25 | 97% 28/29 | 0.86 | 1 | 1 |
| cpb | 0% 0/25 | 3% 1/29 | 0.37 | 1 | 1 |
| cpb2 | 0% 0/25 | 3% 1/29 | 0.37 | 1 | 1 |
| ext | 0% 0/25 | 3% 1/29 | 0.37 | 1 | 1 |
| nanH | 100% 25/25 | 100% 29/29 | 0.86 | 1 | 1 |
| plc | 100% 25/25 | 100% 29/29 | 0.86 | 1 | 1 |
Fisher exact test, two-sided; Benjamini–Hochberg across all 18 genes. Odds ratios use a 0.5 continuity correction, so a zero cell still gives a finite value. Shaded rows are the seven with q < 0.05 reported on the poster.
netE, netF, tpeL and netB were found in no bloodstream genome and no human non-bloodstream genome, but readily in animal-source genomes. Because these genes are not uniformly represented across annotation sources, each row states which call it comes from, and both available calls are given where they differ.
| Gene | Call from | Blood | Human | Animal | Food | Env. |
|---|---|---|---|---|---|---|
| netE | Bakta | 0% 0 | 0% 0 | 19% 27 | 0% 0 | 0% 0 |
| netF | Bakta | 0% 0 | 0% 0 | 20% 28 | 0% 0 | 0% 0 |
| tpeL | Bakta | 0% 0 | 0% 0 | 12% 17 | 0% 0 | 0% 0 |
| netB | ABRicate/VFDB | 0% 0 | 0% 0 | 36% 51 | 0% 0 | 0% 0 |
| netB | Bakta | 0% 0 | 0% 0 | 16% 23 | 0% 0 | 0% 0 |
| cpe | Bakta | 0% 0 | 38% 11 | 26% 36 | 94% 16 | 0% 0 |
| cpe | ABRicate/VFDB | 4% 1 | 38% 11 | 26% 37 | 94% 16 | 0% 0 |
| cpe | toxinotype F (reported) | 4% 1 | 28% 8 | 24% 34 | 94% 16 | 0% 0 |
netE, netF and tpeL are not represented in VFDB and were assessed from Bakta annotation; netB is reported on the poster from ABRicate/VFDB. The two sources disagree for netB: 51 animal genomes by VFDB against 23 by Bakta. For the four enteric-disease toxin genes the absence from the bloodstream and human groups is nonetheless unanimous across both sources — every cell is zero either way — so the poster's conclusion does not turn on this choice. cpe is listed for context and is the one gene here where the sources differ inside the bloodstream group (VFDB 1 of 25, Bakta 0 of 25; the poster reports the VFDB call throughout). Toxinotype F is the reported typing of the source studies, an independent assignment rather than a gene call, and is given for comparison. Denominators are the 218 genomes used for gene detection; the 8 CheckM2 failures (7 animal, 1 food) carry none of these genes, so including them changes no numerator.
All 15 resistance genes detected across the 218 genomes. The poster has room for 10; the 5 omitted from it are lnu(D), mef(A), aadE, erm(A), erm(B), each found in one or two animal-source genomes and in no bloodstream isolate.
| Gene | Blood (25) | Human (29) | Animal (141) | Food (17) | Env. (6) |
|---|---|---|---|---|---|
| tetA(P) | 56% 14 | 66% 19 | 65% 92 | 41% 7 | 83% 5 |
| tetB(P) | 48% 12 | 28% 8 | 37% 52 | 0% 0 | 83% 5 |
| ant(6)-Ib | 8% 2 | 0% 0 | 5% 7 | 0% 0 | 0% 0 |
| erm(Q) | 12% 3 | 3% 1 | 1% 2 | 0% 0 | 50% 3 |
| tet(44) | 8% 2 | 0% 0 | 5% 7 | 0% 0 | 0% 0 |
| lnu(P) | 8% 2 | 3% 1 | 4% 5 | 0% 0 | 0% 0 |
| aph(2'')-Ii | 8% 2 | 0% 0 | 4% 5 | 0% 0 | 0% 0 |
| lnu(D) — not on poster | 0% 0 | 0% 0 | 1% 2 | 0% 0 | 0% 0 |
| mef(A) — not on poster | 0% 0 | 0% 0 | 1% 2 | 0% 0 | 0% 0 |
| aadE — not on poster | 0% 0 | 0% 0 | 1% 1 | 0% 0 | 0% 0 |
| erm(A) — not on poster | 0% 0 | 0% 0 | 1% 1 | 0% 0 | 0% 0 |
| erm(B) — not on poster | 0% 0 | 0% 0 | 1% 1 | 0% 0 | 0% 0 |
| fexA | 0% 0 | 0% 0 | 1% 1 | 0% 0 | 0% 0 |
| optrA | 0% 0 | 0% 0 | 1% 1 | 0% 0 | 0% 0 |
| tet(M) | 0% 0 | 0% 0 | 1% 1 | 0% 0 | 0% 0 |
ABRicate against NCBI AMRFinderPlus, one uniform run over all 218 genomes, collapsed to presence/absence per genome (a gene may hit a genome more than once). optrA and fexA co-occur in a single animal genome (2C45, chicken caecum, China 2018, which also carries erm(A)); tet(M) is in one other (CP-41, peacock, China 1984). Both lie in phylogroup III, the phylogroup that contains every bloodstream isolate of this study. Tetracycline efflux genes are widespread across every source.
| Toxinotype | Blood | Human | Animal | Food | Env. |
|---|---|---|---|---|---|
| A | 96% 24 | 59% 17 | 43% 60 | 6% 1 | 100% 6 |
| B | 0% 0 | 0% 0 | 2% 3 | 0% 0 | 0% 0 |
| C | 0% 0 | 3% 1 | 10% 14 | 0% 0 | 0% 0 |
| D | 0% 0 | 3% 1 | 4% 6 | 0% 0 | 0% 0 |
| E | 0% 0 | 7% 2 | 1% 2 | 0% 0 | 0% 0 |
| F | 4% 1 | 28% 8 | 24% 34 | 94% 16 | 0% 0 |
| G | 0% 0 | 0% 0 | 16% 22 | 0% 0 | 0% 0 |
| Total | 25 | 29 | 141 | 17 | 6 |
Toxinotypes follow the scheme of Rood et al.6; those of the comparison genomes are as reported by Abdel-Glil et al.1 Every genome carries an assignment, so the columns sum to the group sizes. These are reported typings, not gene calls; see 5.2 for the cpe detections, which give a slightly different count.
Every source string exactly as recorded by the reference studies, for all 193 comparison genomes. The poster collapses these into five groups; the original strings are given here because several carry anatomical or clinical detail that the grouping discards.
| Group | Source as recorded | n |
|---|---|---|
| Human | Human (stool) | 22 |
| Human | 4 | |
| Human (gut) | 2 | |
| Human (vagina) | 1 | |
| Animal | Bird | 44 |
| Dog | 16 | |
| Horse | 16 | |
| Sheep | 14 | |
| Pig | 13 | |
| Cattle | 5 | |
| Rabbit | 5 | |
| Bird (chicken intestine) | 3 | |
| Bird (caecum) | 2 | |
| Camel | 2 | |
| Deer | 2 | |
| Llama (ileum) | 2 | |
| Llama (jejunum) | 2 | |
| Mink | 2 | |
| Bird (chicken) | 1 | |
| Bird (pigeon) | 1 | |
| Bird (turkey intestine) | 1 | |
| Cattle (calf intestinal content) | 1 | |
| Mouse (gut) | 1 | |
| Mule | 1 | |
| Mummified animal (Intestine of a frozen mummy - Bison priscus) | 1 | |
| Mummified animal (remains of an ancient "Tumat" puppy - Canis lupus familiaris) | 1 | |
| Peacock | 1 | |
| Pig (feces) | 1 | |
| Pig (gut) | 1 | |
| Pig (intestine) | 1 | |
| Sheep (lamb intestine) | 1 | |
| Food | Restaurant | 7 |
| Food (salted beef) | 2 | |
| Food (beef) | 1 | |
| Food (braised heart) | 1 | |
| Food (meat rissole) | 1 | |
| Food (minced beef) | 1 | |
| Food (rissoles) | 1 | |
| Food (roast meat) | 1 | |
| Food (steamed Lamb) | 1 | |
| Food (stewed steak) | 1 | |
| Environment | Environment (river water) | 3 |
| Environment (Aquarium water) | 1 | |
| Environment (landfill leachate sludge) | 1 | |
| Environment (soil) | 1 |
Every genome in the comparison set was matched to the reference study by accession or by strain name (193 of 193 matched). Four human genomes are recorded only as “Human” with no anatomical site; three of these carry a recorded enteric disease (food poisoning, sporadic diarrhoea, enteritis necroticans / Darmbrand). The one gall-bladder isolate of that collection, NCTC 2544 (ERR1456745), exists only as raw reads with no public assembly and could not be included. The bloodstream group is 25 genomes: 24 from this study and one public genome (2023/00056), which are not part of this table.
| Step | Genomes |
|---|---|
| Public genomes compiled | 204 |
| − unassignable isolation source (cannot be excluded as bloodstream) | −8 |
| − redundant (MD5 of sequence content, not name) | −2 |
| Public genomes analysed | 194 |
| + bloodstream isolates, this study | +24 |
| Gene detection and statistics | 218 |
| − CheckM2 fail (completeness < 95% or contamination > 5%; all animal or food) | −8 |
| Pan-genome and phylogeny | 210 |
| Step | Tool |
|---|---|
| Assembly | shovill (SPAdes), --depth 150 --minlen 200 |
| Quality | CheckM2 |
| Annotation | Prokka — one version, identical settings, all 218 genomes |
| Pan-genome | Roary -e --mafft; default -i 95, paralogues split |
| Phylogeny | snp-sites → IQ-TREE 2, GTR+F+I+R7, UFBoot 1000 + SH-aLRT 1000, constant sites supplied |
| Network | SplitsTree4 4.19.2, NeighborNet — 863 splits, fit 99.95% |
| Gene detection | ABRicate against VFDB and NCBI AMRFinderPlus, one uniform run over all 218 |
| netE, netF, tpeL | Not represented in VFDB → Bakta v1.12.0, database v6.0 (full), likewise all 218 |
| Statistics | Fisher exact test, Benjamini–Hochberg correction across the 18 genes detected; q < 0.05 |
Seven of the 218 genomes were annotated with Bakta v1.11.4 against the same reference database; annotation output was verified identical to v1.12.0 across 14,427 CDS (coordinates, protein sequences, gene symbols, products and cross-references all matched exactly).
The single public bloodstream genome is strain 2023/00056 (GCA_963424335.1, BioProject PRJEB65712, BioSample SAMEA114318975), a binary-enterotoxin-producing isolate from a fatal case of pneumonia and septic shock reported by Ben Saïd et al.3 It falls within phylogroup II, not phylogroup III.
The tables in this handout are generated from the frozen analysis by scripts/make_supplement_tables.py and are available as tab-separated files alongside this page, under data/.
Disclosures. S.H. has received honoraria from Denka, Eiken Chemical, KYORIN Pharmaceutical, Meiji Seika Pharma, MSD, Pfizer, Shionogi and Ushio. All other authors report nothing to disclose.
Contact. Aiko Okazaki, MD, PhD · Project Assistant Professor, UTOPIA Center, The University of Tokyo · okazaki-aiko321@g.ecc.u-tokyo.ac.jp
Affiliations. 1 UTOPIA, The University of Tokyo · 2 The University of Tokyo Hospital · 3 Juntendo University · 4 The University of Tokyo · 5 Toho University School of Medicine · 6 Institute of Science Tokyo