Oncofusion

start / graph / ana-depmap-kmt2a-selective-scan-full

Analysismoderate confidence 3 of 5: One strong study, or several smaller consistent ones.king

The whole protein-coding genome: menin is still first, and a metabolic pathway appears (run)

also called full-genome selectivity scan · 17,728-gene scan · haem pathway signal

The earlier scan covered a quarter of our genes and put menin first; this one covers all of them and menin is still first, so nothing hidden in the rest outranks it. The new result is that four genes from a single chemical pathway — the one cells use to build haem — all rank very high, a pattern that was not visible before.

The same thing, for a scientist

🟠 **ARM-COMPOSITION CAVEAT ADDED 2026-09-24 (cycle 140).** This entry's computation rests on `analyses/depmap-kmt2a-selective-scan-full`, and that scan **fixes both comparison arms from the first gene that returns screen data** (`run.py:170`, `if rear is None:`; line 174 prints *"Groups fixed from first gene with data"*). Because DepMap screens each gene in a gene-specific set of lines, **5 of 16 eligible KMT2A-rearranged lines are excluded from every gene comparison, including genes where they were screened** - 16 resolve under the scan's own query, 11 were used. **This does not show any conclusion here is wrong.** Dropping 5 of 16 costs power, and it biases a ranking only if the dropped lines differ systematically from the kept ones, **which nobody has tested.** The named one-line fix is to derive the arms from each gene's own screen list. Until the scan is re-run, selectivity ranks taken from it carry an unquantified arm-composition caveat. Found by `analyses/depmap-scan-line-identities` on 2026-09-21 and unreported until now; see `ana-the-depmap-selectivity-scan-fixes-its-arms-from-the-first-gene-and-drops-five-of-sixteen-lines-from-every-comparison`. --- Extends the 4,589-gene scan to the full HGNC protein-coding set. 17,728 genes scored; **MEN1 ranks 1**, FLT3 2, PSIP1 7, MEF2C 8, KMT2A 11. The seven cycle-22 strikes occupy the 66th to 98th percentiles. New at genome scale: **four haem-biosynthesis and iron-handling genes in the top 47** — ALAD 14, IREB2 15, UROS 41, ALAS1 47 — all selective, none pan-essential. That pathway-level pattern contradicts the single-gene explanation this project adopted for UROS one day earlier.

Go deeper
The question this settles

The 4,589-gene scan left an obvious open question, recorded in its own entry: "Scanning the remaining ~14,000 protein-coding genes would cost roughly four hours of API time. Is anything above MEN1 hiding there?"

No. Across 17,728 genes — the full HGNC protein-coding set with DepMap data — MEN1 is the single most selectively required gene in KMT2A-rearranged leukemia lines, at the 0.006th percentile 1.

That is the strongest form this project's central premise can take from this dataset. It is also the form most exposed to being wrong: MEN1 competed against every measurable gene, unflagged, and a poor rank would have been unmissable.

The ranking holds, and so does the strike list
RankGenedifferencepan-essentiality
1MEN1-0.866-0.363the premise
2FLT3-0.843-0.147partly FLT3-ITD lines
3SLC5A3-0.740-0.163novel; myo-inositol transporter
4DNAJC9-0.645-1.393pan-essential, discard
5CCND3-0.623-0.348independently confirmed this cycle
7PSIP1 (LEDGF)-0.565+0.028the tether
8MEF2C-0.560-0.039
9MED13L-0.551-0.156novel; Mediator
11KMT2A-0.541-0.206
12IRF8-0.539+0.037novel; myeloid TF, highly selective
13CDK6-0.527-0.632validates the CCND3 finding
14ALAD-0.500-0.088haem biosynthesis
15IREB2-0.494-0.215iron regulation
16MED23-0.483-0.210novel; Mediator
17MYB-0.475-0.384still unresolved
23JMJD1C-0.454-0.012a published MLL-fusion dependency
41UROS-0.387-0.198haem biosynthesis
47ALAS1-0.378-0.119haem biosynthesis, rate-limiting

The cycle-22 strikes hold at genome scale, having been reached without any knowledge of them: CDC73 11,645 (65.7th percentile), KDM4C 14,509, RBBP5 14,622, MLLT3 14,902, ASH2L 17,248, KDM6A 17,398 and CTR9 17,413 of 17,728 (98.2nd percentile) — every one with a positive difference, i.e. marginally less required in rearranged lines. DPY30 sits at 2,590 (14.6th), unremarkable for the sixth independent design.

JMJD1C at rank 23 is a useful check. It is a documented MLL-fusion dependency the scan was not told about, and its near-zero pan-essentiality (-0.012) means it is selective rather than generally required. Together with ZFP64 in the smaller scan, that is two known mechanisms recovered from dependency data alone.

The haem cluster, and why it overturns a conclusion this project reached yesterday1 deeper

Four genes in iron and haem handling rank in the top 47 of 17,728 — the top 0.27% — and none is pan-essential:

Generankpathway positiondifferencepan-essentiality
ALAD14haem biosynthesis, enzyme 2-0.500-0.088
IREB215iron-responsive element-binding protein 2-0.494-0.215
UROS41haem biosynthesis, enzyme 4-0.387-0.198
ALAS147haem biosynthesis, enzyme 1, rate-limiting-0.378-0.119

All four accessions were resolved by exact gene symbol and the primary gene name read back 2.

This contradicts the explanation this project adopted one day earlier. UROS appeared at rank 15 in the smaller scan; a literature search found nothing connecting it to leukemia, and the beat recorded a deliberate clean negative with a mechanistic explanation: UROS loss causes accumulation of porphyrin isomer I, the toxic-metabolite mechanism of congenital erythropoietic porphyria, which would harm any cell with upstream haem flux regardless of the fusion.

That was a careful, well-reasoned conclusion on the evidence available, and the pathway view breaks it. The gain-of-toxic-metabolite model makes a directional prediction: reduce flux into the pathway and the toxicity should fall. ALAS1 is the first and rate-limiting enzyme — knocking it out reduces flux, so under that model ALAS1 loss should be neutral or protective. Instead ALAS1 loss is selectively deleterious, and so is ALAD's at enzyme 2. Blocking the tap does not rescue; it hurts more.

The simpler reading is that KMT2A-rearranged lines have an elevated requirement for haem synthesis itself. This project's own graph already holds an unexplored thread in that direction: SETD1A controls pause release at haem genes. That was noted and set aside as pointing the wrong way. With four pathway members ranking, it deserves a proper look.

What would distinguish the two explanations

The two models diverge cleanly, and the test is cheap:

  • Elevated requirement: haem-pathway genes should be more expressed in rearranged lines, and the dependency should scale with flux. A succinylacetone or ALAS1-inhibitor dose-response in rearranged versus non-rearranged lines would show it.
  • Toxic intermediate: only enzymes downstream of the branch point should matter, and upstream loss should rescue. This is already contradicted by ALAS1 and ALAD ranking high.
  • Confounded by erythroid identity: the strongest alternative. If any rearranged lines have erythroid or megakaryocytic character they will carry high haem demand for reasons unrelated to the fusion. The 11 lines are 7 AML and 4 B-ALL by annotation, and checking their lineage subtype is the first thing to do — it is free, and if it explains the signal the rest is moot.

Until that check is done this is a pattern, not a finding.

What has not changed

Eleven lines, none infant-derived. The panel is unchanged and remains the binding limitation on every conclusion here.

Same statistic as the specificity analysis. This is not methodologically independent of it; what is new is the absence of a prior and the complete ranking context.

Pan-essential contaminants persist. DNAJC9 at 4, TINF2 at 6, RPL13A at 19 and MCL1 at 18 are more depleted in rearranged lines because the panel is more sensitive, not because of the fusion. CDK6 at 13 sits on the boundary at -0.632 and should be read with that in mind even though its biology is supported elsewhere.

A rank is a lead, not a target. The scan-hits beat found that a rank buys about one afternoon of literature search, and that roughly half the time that search pays.

Facts and figures
scriptanalyses/depmap-kmt2a-selective-scan-full/run.py
statusrun
run date2026-09-18
answers open question fromana-depmap-kmt2a-selective-scan
seed20260918
n scanned17728
n kmt2a rearranged11
n other leukemia43
men1 rank1
men1 percentile0.006
men1 pan essentiality-0.3631
overturnsthe single-gene 'gain-of-toxic-metabolite' explanation adopted for UROS in ana-uros-kmt2a-literature-search-negative, which predicts that reducing upstream flux is protective; ALAS1 and ALAD loss are instead selectively deleterious
figureanalyses/depmap-kmt2a-selective-scan-full/results.md
arm composition caveat cycle 140source scan fixes both arms from the first gene with data (run.py:170); 5 of 16 eligible KMT2A-rearranged lines are excluded from every comparison. Effect on this entry's conclusion untested.
haem iron cluster
ALAD14.0ALAD: 14IREB215.0IREB2: 15UROS41.0UROS: 41ALAS147.0ALAS1: 47
struck gene ranks
CDC7311,645CDC73: 11645KDM4C14,509KDM4C: 14509RBBP514,622RBBP5: 14622MLLT314,902MLLT3: 14902ASH2L17,248ASH2L: 17248KDM6A17,398KDM6A: 17398CTR917,413CTR9: 17413
notable low ranks
DOT1L184.0DOT1L: 184DPY302,590DPY30: 2590HOXA93,649HOXA9: 3649ELL23,697ELL2: 3697
haem iron cluster — the numbers behind the chart
{
  "ALAD": 14,
  "IREB2": 15,
  "UROS": 41,
  "ALAS1": 47
}
novel selective top20 no entry
[
  "SLC5A3",
  "MED13L",
  "IRF8",
  "MED23",
  "MBNL1"
]
recovered known mechanisms
{
  "JMJD1C": 23,
  "CDK6": 13
}
struck gene ranks — the numbers behind the chart
{
  "CDC73": 11645,
  "KDM4C": 14509,
  "RBBP5": 14622,
  "MLLT3": 14902,
  "ASH2L": 17248,
  "KDM6A": 17398,
  "CTR9": 17413
}
notable low ranks — the numbers behind the chart
{
  "DPY30": 2590,
  "HOXA9": 3649,
  "ELL2": 3697,
  "DOT1L": 184
}
pan essential in top20
[
  "DNAJC9",
  "TINF2",
  "MCL1",
  "RPL13A"
]
How this could bite us back
  • MEN1 first out of the entire protein-coding genome is the strongest support this project's premise has, and the same panel that produced it puts four pan-essential housekeeping genes in the top 20. The ranking is exact; what it means biologically is not.
  • The haem cluster contradicts a conclusion one of this project's own beats reached carefully a day earlier. That beat was not sloppy — it had one gene and no pathway context. The lesson is that a clean negative on a single gene is not a clean negative on its pathway, and the correction belongs to the method, not to that agent.
  • Nothing here is a target. SLC5A3, MED13L, IRF8 and the haem enzymes are leads with no druggability assessment, and the cheapest explanation for the haem cluster is a lineage confound that nobody has excluded yet.
  • Eleven cell lines, none from an infant, still stand behind every number in this entry.
Still open
  1. What are the lineage subtypes of the 11 rearranged lines? If any have erythroid or megakaryocytic character the haem cluster is explained and needs no further work. This is free and should be done first.
  2. Are haem-pathway genes more highly expressed in KMT2A-rearranged lines? Expression data would separate an elevated requirement from a structural artefact of the panel.
  3. This graph already records that SETD1A controls pause release at haem genes, noted once and set aside. With four pathway members ranking, does that thread connect?
  4. SLC5A3 ranks 3 with no entry anywhere and no obvious leukemia literature. Is it a real lead or the next UROS?
  5. MED13L and MED23 both rank in the top 16 and MED26 is already in the graph as a SEC link. Is Mediator under-represented in this project's interactome?
  6. DOT1L ranks 184 and HOXA9 3,649. How many other genes in this scan are paralogue-buffered and therefore invisible, and can that be estimated rather than guessed?
Every source behind this entry (2)
  1. E1moderate confidence 3 of 5: One strong study, or several smaller consistent ones.dataset

    Across 17,728 HGNC protein-coding genes with DepMap gene-effect data, ranked by difference in mean effect between 11 Cellosaurus-annotated KMT2A-rearranged and 43 other leukemia lines, MEN1 ranks 1 (difference -0.866, 0.006th percentile, pan-essentiality -0.363). FLT3 ranks 2, SLC5A3 3, CCND3 5, PSIP1 7, MEF2C 8, MED13L 9, KMT2A 11, IRF8 12, CDK6 13, ALAD 14, IREB2 15, MED23 16, MYB 17, JMJD1C 23, UROS 41 and ALAS1 47. The seven genes struck in the cycle-22 target rewrite rank from 11,645 (CDC73) to 17,413 (CTR9), all with positive differences; DPY30 ranks 2,590 and HOXA9 3,649.

    What this evidence cannot show

    Our own computation, not peer reviewed. Eleven rearranged lines, none infant-derived, is the binding limit. Same difference statistic and same line split as ana-depmap-kmt2a-rearranged-specificity, so not methodologically independent of it; what is new is the absence of a prior and the complete ranking. Several top-20 entries are pan-essential genes depleted more in a more sensitive panel (DNAJC9, TINF2, RPL13A, MCL1) and CDK6 sits on that boundary. Cellosaurus fusion annotation may be incomplete, biasing toward the null. DepMap release version is not stated by the API.

    Evidence that argues against this (2)
    • MEN1 ranking first is the result the project most wants and should be held to the strictest reading. The ranking is over a panel of 11 lines whose composition also produces four pan-essential genes in the top 20, so the same panel effect that inflates those inflates MEN1's difference to an unknown degree. The percentile is exact; the biology it implies is not. analysis run by this project 2026
    • HOXA9 ranks 3,649 of 17,728 and DOT1L 184. The canonical downstream effector and a central mechanistic player both rank far below genes with no established connection to this biology, which is a standing reason to treat the ranking as a lead generator rather than a mechanistic statement. Paralogue buffering explains HOXA9 but nothing in the scan flags which other genes are similarly buffered. analysis run by this project 2026
  2. E2established confidence 5 of 5: Multiple independent lines of evidence; treated as settled in the field.dataset

    ALAS1 (P13196) is 5-aminolevulinate synthase, the first and rate-limiting enzyme of haem biosynthesis; ALAD (P13716) is delta-aminolevulinic acid dehydratase, the second enzyme; UROS (P10746) is uroporphyrinogen-III synthase, the fourth; and IREB2 (P48200) is iron-responsive element-binding protein 2. All four accessions were resolved by exact gene-symbol search restricted to reviewed human entries with the primary gene name confirmed in the response.

    What this evidence cannot show

    Pathway position is standard biochemistry and is stated from the UniProt protein names; no pathway database was queried to confirm ordering. Membership of a pathway does not imply the dependency runs through that pathway's function.

    Evidence that argues against this (1)
    • Four genes sharing a pathway can rank highly for a reason unrelated to the pathway. The strongest alternative here is lineage: erythroid or megakaryocytic character in any of the 11 rearranged lines would raise haem demand independently of the fusion. That check has not been done and is free, so this cluster should be treated as a pattern rather than a finding until it is. analysis run by this project 2026
Connected to