P-202  ·  IDWeek 2026  ·  supplementary handout

Genomic Characterization of Clostridium perfringens Bacteremia Isolates at a Tertiary Care Center in Japan

Aiko Okazaki1, Shu Okugawa2, Satoshi Kitaura2, Mahoko Ikeda3, Shintaro Yanagimoto4, Yusuke Nomura2, Saho Koyano2, Yoshimi Higurashi2, Sohei Harada5, Ryoichi Saito6, Takeya Tsutsumi4

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.

1  Corrections to the published abstract

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.

ItemAbstractRe-analysis
Public genomes192194
Human, non-bloodstream, n31 †29
Bloodstream, n24 / 2525
Core genes5381,215
Total genes30,45423,512
nagK100%92%
nagI96%100%
nanJ96%100%
Statisticsp, uncorrectedq, BH–FDR
nagL96% vs 45%withdrawn — 56% (14/25) vs 45% (13/29) q = 0.958
pfoA92% 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.

2  Robustness checks

3  Core-gene thresholds

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 recovered8221,2151,3761,5631,6852,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.

4  Node support and core alignment

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.

SupportMedianMeanMinimumNodes at or above threshold
SH-aLRT100.084.5 0 168 / 205 at ≥ 80
UFBoot100.090.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.

Alignment accounting

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.

ComponentSites
Roary core gene alignment1,066,463 210 genomes
SNP sites passed to IQ-TREE56,596 snp-sites output
Constant sites supplied964,325 as -fconst, so branch lengths are not inflated
Total modelled1,020,921 as reported by IQ-TREE
Gap or ambiguous columns45,542 neither SNP nor unambiguous constant

5  Gene tables

5.1  Virulence genes, bloodstream vs non-bloodstream

All 18 genes tested, not only the seven that reached significance. Genes detected in neither group are not listed and were not tested.

GeneBloodstream (25)Human non-blood (29) ORpq
nagI100% 25/2559% 17/2936.430.0001660.00245
nagH100% 25/2562% 18/2931.700.0004080.00245
nanJ100% 25/2562% 18/2931.700.0004080.00245
nagJ100% 25/2566% 19/2927.460.0009740.00438
nanI100% 25/2569% 20/2923.630.002270.00816
cpe4% 1/2538% 11/290.100.003060.00919
nagK92% 23/2559% 17/296.710.006040.0155
pfoA80% 20/2552% 15/293.490.04550.102
iap0% 0/257% 2/290.220.4930.888
ibp0% 0/257% 2/290.220.4930.888
nagL56% 14/2545% 13/291.540.5860.958
cloSI100% 25/25100% 29/290.8611
colA96% 24/2597% 28/290.8611
cpb0% 0/253% 1/290.3711
cpb20% 0/253% 1/290.3711
ext0% 0/253% 1/290.3711
nanH100% 25/25100% 29/290.8611
plc100% 25/25100% 29/290.8611

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.

5.2  Enteric toxin genes

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.

GeneCall fromBloodHumanAnimalFoodEnv.
netEBakta0% 00% 019% 270% 00% 0
netFBakta0% 00% 020% 280% 00% 0
tpeLBakta0% 00% 012% 170% 00% 0
netBABRicate/VFDB0% 00% 036% 510% 00% 0
netBBakta0% 00% 016% 230% 00% 0
cpeBakta0% 038% 1126% 3694% 160% 0
cpeABRicate/VFDB4% 138% 1126% 3794% 160% 0
cpetoxinotype F (reported)4% 128% 824% 3494% 160% 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.

5.3  Resistance genes by source

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.

GeneBlood (25)Human (29)Animal (141)Food (17)Env. (6)
tetA(P)56% 1466% 1965% 9241% 783% 5
tetB(P)48% 1228% 837% 520% 083% 5
ant(6)-Ib8% 20% 05% 70% 00% 0
erm(Q)12% 33% 11% 20% 050% 3
tet(44)8% 20% 05% 70% 00% 0
lnu(P)8% 23% 14% 50% 00% 0
aph(2'')-Ii8% 20% 04% 50% 00% 0
lnu(D) — not on poster0% 00% 01% 20% 00% 0
mef(A) — not on poster0% 00% 01% 20% 00% 0
aadE — not on poster0% 00% 01% 10% 00% 0
erm(A) — not on poster0% 00% 01% 10% 00% 0
erm(B) — not on poster0% 00% 01% 10% 00% 0
fexA0% 00% 01% 10% 00% 0
optrA0% 00% 01% 10% 00% 0
tet(M)0% 00% 01% 10% 00% 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.

5.4  Toxinotype by source

ToxinotypeBloodHumanAnimalFoodEnv.
A96% 2459% 1743% 606% 1100% 6
B0% 00% 02% 30% 00% 0
C0% 03% 110% 140% 00% 0
D0% 03% 14% 60% 00% 0
E0% 07% 21% 20% 00% 0
F4% 128% 824% 3494% 160% 0
G0% 00% 016% 220% 00% 0
Total2529141176

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.

6  Isolation sources, as recorded

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.

GroupSource as recordedn
HumanHuman (stool)22
Human4
Human (gut)2
Human (vagina)1
AnimalBird44
Dog16
Horse16
Sheep14
Pig13
Cattle5
Rabbit5
Bird (chicken intestine)3
Bird (caecum)2
Camel2
Deer2
Llama (ileum)2
Llama (jejunum)2
Mink2
Bird (chicken)1
Bird (pigeon)1
Bird (turkey intestine)1
Cattle (calf intestinal content)1
Mouse (gut)1
Mule1
Mummified animal (Intestine of a frozen mummy - Bison priscus)1
Mummified animal (remains of an ancient "Tumat" puppy - Canis lupus familiaris)1
Peacock1
Pig (feces)1
Pig (gut)1
Pig (intestine)1
Sheep (lamb intestine)1
FoodRestaurant7
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
EnvironmentEnvironment (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.

7  Analysis set

StepGenomes
Public genomes compiled204
− unassignable isolation source (cannot be excluded as bloodstream)−8
− redundant (MD5 of sequence content, not name)−2
Public genomes analysed194
+ bloodstream isolates, this study+24
Gene detection and statistics218
− CheckM2 fail (completeness < 95% or contamination > 5%; all animal or food)−8
Pan-genome and phylogeny210

8  Software and data

StepTool
Assemblyshovill (SPAdes), --depth 150 --minlen 200
QualityCheckM2
AnnotationProkka — one version, identical settings, all 218 genomes
Pan-genomeRoary -e --mafft; default -i 95, paralogues split
Phylogenysnp-sites → IQ-TREE 2, GTR+F+I+R7, UFBoot 1000 + SH-aLRT 1000, constant sites supplied
NetworkSplitsTree4 4.19.2, NeighborNet — 863 splits, fit 99.95%
Gene detectionABRicate against VFDB and NCBI AMRFinderPlus, one uniform run over all 218
netE, netF, tpeLNot represented in VFDB → Bakta v1.12.0, database v6.0 (full), likewise all 218
StatisticsFisher 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/.

References

  1. Abdel-Glil MY et al. Sci Rep 2021;11:6756.
  2. Abdel-Glil MY et al. Microbiol Spectr 2021;9:e00533-21.
  3. Ben Saïd L et al. Microorganisms 2024;12:1095.
  4. Okazaki A et al. Int J Infect Dis 2025;151:107358.
  5. Kiu R et al. Nat Microbiol 2023;8:1160–75.
  6. Rood JI et al. Anaerobe 2018;53:5–10.

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