SciFlow lab notebook · Shared read-only

β-Lactamase Phylogeny — SciFlow

Reconstruct the evolutionary history of the serine β-lactamases — the most clinically important antibiotic-resistance enzymes in Gram-negative bacteria — directly from protein sequences.

Overview Methods & provenance Results Raw artifacts

Run date

2026-05-30 → 2026-05-31

Stages

3 3 completed

Pipeline

mafft → trim_alignment → iqtree

Software

3 container images mafft, trimal, iqtree2

Outputs

26 files 1.8 MB

Overview

Reconstruct the evolutionary history of the serine β-lactamases — the most clinically important antibiotic-resistance enzymes in Gram-negative bacteria — directly from protein sequences. The curated demo FASTA holds Ambler class A families (TEM, SHV, CTX-M, KPC) and Ambler class D (OXA-family) enzymes, rooted on a genuine class C AmpC cephalosporinase outgroup (E. coli K-12 AmpC, UniProt P00811). You get back a rooted maximum-likelihood tree with bootstrap support, model-selection diagnostics, and a clean alignment you can actually trust. Built for AMR surveillance teams, structural biologists mapping active-site evolution, and anyone teaching molecular evolution who wants a real-world example instead of toy data. Tips are labelled by family/class (e.g. TEM__, SHV__, KPC__, OXA__/clD__ for class D, AmpC__ for the class C outgroup) so you can verify directly that each family is recovered as a monophyletic clade — the point of a known-answer demo: confirm the pipeline gets the right answer before running it on novel sequences.

The methodology is the substance, not the topology alone. Alignment uses MAFFT (L-INS-i for accuracy at this size), trimAl removes ambiguously aligned columns under an automated heuristic, ModelFinder picks the best amino-acid substitution model by BIC, and IQ-TREE infers the tree with 1000 ultrafast bootstrap replicates and an SH-aLRT branch test, rooted on the AmpC outgroup. The output panel reports model choice, the SH-aLRT/UFBoot support distribution, alignment length retained after trimming, and class A vs class D clade recovery — because a tree without those numbers is just a picture.

Methods & per-stage provenance

3 stages, 3 completed. Each stage below lists the parameters it ran with, the software the run recorded, and the files it produced.

Not recorded by this run
  • no figures rendered (no PNG/SVG outputs and no Plotly HTML figure could be extracted)
02 Align (MAFFT-LINSI) mafft 43 min completed
input fasta: betalactamases.fastaalgorithm: linsioutput format: fastaadjust direction: true

Verbatim command / script not recorded in this run's stage record.

Container image: staphb/mafft:7.526 Image digest: sha256:a688373aab143e1d3a458440ac378a98e0f720766143803c89117cc3ec12bbb7
  • alignment.fastaFASTA · 82 KB
  • alignment_report.htmlreport · 113 KB
  • input_data_summary.htmlreport · 4 KB
  • mafft.loglog · 120 KB
  • summary.jsondata · 396 B
03 Trim alignment (trimAl -automated1) trim_alignment 43 min completed
alignment: alignment.fastamethod: automated1

Verbatim command / script not recorded in this run's stage record.

Container image: staphb/trimal:1.5.0 Image digest: sha256:dab106310271bbfe652aa07c378e537b2ab68a4537d172facecefb272d0aa096
  • alignment.trimmed.fastaFASTA · 34 KB
  • alignment.untrimmed.fastaFASTA · 82 KB
  • columns_kept.txttext · 1 KB
  • trimal.loglog · 20 KB
  • trimming_report.htmlreport · 567 KB
  • trimming_report.jsondata · 1 KB
04 ML tree (IQ-TREE 2) iqtree 43 min completed
alrt: trueseed: 42model: MFPoutgroup: AmpC__P00811_ECOLIalignment: alignment.trimmed.fastabootstrap: ufbootsequence type: AAbootstrap reps: 1000

Verbatim command / script not recorded in this run's stage record.

Container image: staphb/iqtree2:2.3.6 Image digest: sha256:14d83f760b25cf61cf7ebcf0322f1de166eddf08ae04c8fa5f656cd2de9902a3
  • iqtree.loglog · 35 KB
  • iqtree_report.htmlreport · 6 KB
  • phylogeny.nwkNWK · 7 KB
  • run_metadata.jsondata · 558 B
  • summary.jsondata · 800 B
  • tree.bionjBIONJ · 6 KB
  • tree.ckp.gzGZ · 320 KB
  • tree.contreeCONTREE · 7 KB
  • tree.iqtreeIQTREE · 81 KB
  • tree.loglog · 35 KB
  • tree.mldistMLDIST · 150 KB
  • tree.model.gzGZ · 11 KB
  • tree.splits.nexNEX · 123 KB
  • tree.treefileTREEFILE · 7 KB
  • tree.uniqueseq.phyPHY · 29 KB

Results

No static figure could be embedded in this notebook — this run recorded no PNG/SVG output, and no interactive plot could be exported to a static image. Any figures this run produced are listed among the raw artifacts below and remain viewable in the project itself.

Raw artifacts

26 output files produced across all stages (1.8 MB total). This report is self-contained; structures and trajectories are listed below for reference.