Skip to content

Metabolite → gene links: transporter links via KEGG BRITE ko02000 - #89

Merged
sr320 merged 4 commits into
mainfrom
metabolite-links/transporters
Oct 17, 2026
Merged

sr320 merged 4 commits into
mainfrom
metabolite-links/transporters

Conversation

@sr320

@sr320 sr320 commented Oct 17, 2026

Copy link
Copy Markdown
Owner

Three commits from the claude/recursing-varahamihira-f8dfa9 session that landed after #73 (metabolite_linked_convergence tier) had already merged, so they never reached main. They're moved to a fresh branch, with current main merged in (no conflicts).

Commits

  1. 033d3c1 Add transporter metabolite → gene links via KEGG BRITE ko02000
    • KEGG BRITE ko02000 classifies transporter KOs into families but records no substrates. The KO → compound step is therefore curated in data/mappings/transporter_substrates.tsv: 43 rows, each with substrate_basis, a human ortholog and a note. For example: GLUT/SGLT → glucose, MCT/SMCT → lactate/pyruvate, NaDC/DIC → succinate, BGT1 → betaine, system A/ASCT/PAT → alanine.
    • The builder fetches br:ko02000, rejects curated KOs that aren't in it, and records each KO's BRITE family path.
    • Transporter links are always inferred. The hub cap is applied per link type, so transporter links never strip a compound's enzyme-link credit.
    • The link table gains link_basis and brite_family. Evidence cards show link type, basis and family, with a note on why transporter links rest on curation.
    • First local build: 191 enzyme + 58 transporter links. (S)-lactate, which has no L-LDH enzyme link in crg, now links through MCT/SMCT.
  2. a6e1646 Demo: glucose ↔ simulated glut1 transporter link. Adds a fabricated glut1 gene (LOC105342003) to the demo crosswalk and GIGAS_HEAT01, plus a transporter row with the real GLUT1 KO K07299 and its BRITE path. The demo now shows both link types reaching metabolite_linked_convergence.
  3. 8c9abe2 Demo: move GIGAS_GROW06 glucose off the adj. p = 0.05 boundary. Glucose sat exactly at padj 0.05, so the demo case depended on the gate using <=. Its p-value goes from 0.021 to 0.006 and its padj from 0.05 to 0.02 (simulated data). A test checks that every demo link-convergence feature has adj. p ≤ 0.04.

Checked locally (after merging current main)

  • ruff clean; pytest 730 passed
  • CI demo pipeline: register → harmonize → meta-analyze → build-evidence-cards → top-candidates completes. The glucose evidence card shows the transporter link.

All demo changes are to simulated: true data and the synthetic demo crosswalk; no real study is affected.

🤖 Generated with Claude Code

sr320 and others added 4 commits October 17, 2026 06:18
KEGG BRITE ko02000 classifies transporter KOs into families but records no
substrate compounds, so the KO -> compound step is curated:

- data/mappings/transporter_substrates.tsv: 43 curated rows (GLUT/SGLT and
  SLC60A2 -> glucose; MCT1-4 and SMCT -> lactate/pyruvate; NaDC and DIC ->
  succinate; BGT1 -> betaine; system A, ASCT and PAT -> alanine), each with
  substrate_basis (ko_definition | substrate_class), human ortholog and note.
- Builder fetches get/br:ko02000, rejects curated KOs not in it, records each
  KO's BRITE family path, and lists curated KOs with no crg gene.
- link_type transporter: always 'inferred'; hub cap applied per link type so
  transporter links never strip a compound's enzyme-link credit.
- Link table gains link_basis and brite_family columns (demo table and its
  checksum updated; demo results unchanged).
- Evidence cards show link type, basis and BRITE family, with a paragraph on
  why transporter links rest on curation.
- First local build: 191 enzyme + 58 transporter links; (S)-lactate, which has
  no L-LDH enzyme link in crg, now links through MCT/SMCT transporters.

Co-Authored-By: Claude Opus 5.5 <noreply@anthropic.com>
Add a fabricated glut1 gene (LOC105342003) to the demo crosswalk and the
GIGAS_HEAT01 chronic_heat_vs_control table, and a transporter row to the
simulated link table (real GLUT1 KO K07299 and its real KEGG BRITE ko02000
family path; fabricated gene ID). GIGAS_GROW06 glucose and glut1 now reach
metabolite_linked_convergence through a transporter link, so the demo shows
both link types. Demo provenance checksum, tests and docs updated.

Co-Authored-By: Claude Opus 5.5 <noreply@anthropic.com>
Glucose sat exactly at padj 0.05, so the glucose <-> glut1 transporter demo
case depended on the gate using <= rather than <. Set p 0.021 -> 0.006 and
padj 0.05 -> 0.02, and add a test that every demo link-convergence feature
has adj. p <= 0.04.

Co-Authored-By: Claude Opus 5.5 <noreply@anthropic.com>
@sr320
sr320 merged commit 2361bdd into main Oct 17, 2026
6 checks passed
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

1 participant