Oncofusion

start / graph / ana-depmap-kmt2a-rearranged-specificity

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

The assumption the whole project rests on, tested at last: menin dependency is fusion-specific (run)

also called KMT2A-r specificity · fusion-specific menin dependency · Cellosaurus rearrangement split

Everything this project argues assumes leukemia cells need menin because they carry the fusion, and nobody here had ever checked it against our own data. Splitting the cell lines by whether they carry the fusion, the ones that do need menin far more — beyond what any of 918 randomly chosen genes would predict.

The same thing, for a scientist

Every co-dependency analysis in this project correlated candidates against MEN1 across leukemia lines without testing whether MEN1 dependency tracks the rearrangement. Cellosaurus gene-fusion annotations cross-referenced to DepMap ids give the split: 11 KMT2A-rearranged versus 43 other leukemia lines. Mean MEN1 gene effect is -1.263 versus -0.397, a difference more extreme than all 918 random genes tested, permutation p = 0.0002, Holm 0.012. It holds within AML alone (7 versus 25, difference -0.890). KMT2A, PSIP1, MEF2C, BCL2, MEIS1 and KAT7 also separate; DPY30 does not, for the fourth time.

Go deeper
Why this had never been tested

The project's therapeutic argument is that menin is druggable because the fusion depends on it. Four analyses here have correlated candidate genes against MEN1 gene effect across leukemia lines. None of them asked whether MEN1 dependency is higher in the lines that actually carry a KMT2A rearrangement — and a correlation across a mixed panel can arise with no fusion-specific effect at all.

The obstacle was annotation, not interest. The Open Targets DepMap API returns a disease label per screen (acute myeloid leukemia, b-cell acute lymphoblastic leukemia) and nothing about rearrangement status, and the DepMap bulk download is captcha-blocked, which this project documented in an earlier dataset audit.

Cellosaurus has it. Cell lines carry sequence-variation records with variation-type: "Gene fusion" and HGNC-linked gene labels, and many are cross-referenced to DepMap ACH ids. Querying var:KMT2A AND dr:DepMap returns 16 DepMap lines with an annotated KMT2A fusion — MV4-11, MOLM-13, MOLM-14, THP-1, RS4;11, SEM, SEMK2, KOPN-8, NOMO-1, Mono-Mac-6, ML-2, SHI-1, OCI-AML-2, HB11;19, ALL-PO and Karpas-45 — of which five carry KMT2A::AFF1, the fusion this project is actually about. Eleven overlap the leukemia panel with MEN1 data 1.

The result3 deeper

Gene effect is negative when a gene is required, so a more negative difference means more required in the rearranged lines. Each difference is tested both by label permutation and against a null of the same statistic computed on 918 random genes, because plenty of genes differ between any two groups of cell lines.

Genemean KMT2A-r (n=11)mean other (n=43)differencepercentile of random genesHolmAML-only difference
MEN1-1.263-0.397-0.8660.00.012-0.890
FLT3-0.962-0.119-0.8430.00.012-0.754
PSIP1 (LEDGF)-0.553+0.012-0.5650.00.012-0.665
MEF2C-0.843-0.283-0.5600.00.012-0.510
KMT2A-0.811-0.270-0.5410.00.012-0.581
BCL2-0.666-0.270-0.3960.00.012-0.387
MEIS1-0.429-0.123-0.3060.220.022-0.198
KAT7-0.650-0.354-0.2960.220.022-0.290
DOT1L-0.604-0.373-0.2310.650.039-0.257
MLLT1-0.318-0.168-0.1492.180.072-0.077
DPY30-0.510-0.447-0.06314.40.259-0.131
ACTB (control)-0.949-0.768-0.1821.42—-0.366
GAPDH (control)-1.730-1.785+0.05585.8—+0.025

MEN1, FLT3, PSIP1, MEF2C, KMT2A and BCL2 are more extreme than every one of the 918 random genes. That is the headline, and it is the first direct evidence in this project, from this project's own data, that the menin dependency is a property of the rearrangement rather than of leukemia in general.

It is not the lineage contrast. The rearranged set is part B-ALL and part AML, so a raw split partly recapitulates the myeloid-versus-lymphoid difference that ana-depmap-within-subtype removed. Restricting to AML only (7 rearranged versus 25 other) the MEN1 difference is -0.890, slightly larger, at p = 0.0011.

The control that did not behave, and what it costs

ACTB sits at the 1.42nd percentile of the random-gene null, meaning it too is nominally more required in the rearranged lines. Its permutation p is 0.16, so it is not significant by that test, and GAPDH goes the other way at the 86th percentile. But ACTB should be at neither extreme.

The most likely explanation is a mild global offset: the rearranged lines may be somewhat more sensitive across the board, through growth rate, screen quality or depth. If so, every difference in the table is inflated by roughly that amount.

This does not overturn the result, because the top six are at the 0.0 percentile — beyond all 918 random genes, and far beyond where ACTB sits. But it does mean the magnitudes should be read as upper bounds, and it is the reason to prefer the percentile column over the raw difference.

FLT3 nearly matches MEN1, and that is a warning not a finding

FLT3's difference (-0.843) is almost as large as MEN1's. This is very unlikely to be about the fusion. Several of the 11 rearranged lines — MV4-11, MOLM-13 and MOLM-14 among them — are FLT3-ITD lines, and FLT3-ITD is a co-occurring lesion, not a consequence of the rearrangement.

So the rearranged group differs from the rest in more than one way, and this design cannot separate the fusion from what travels with it. For MEN1 there is independent mechanistic reason to attribute the dependency to the fusion; for FLT3 there is an obvious alternative. A reader who takes the whole column as 'fusion-specific dependencies' would be over-reading it, and the same caution applies in weaker form to BCL2 and MEF2C.

DPY30 fails a fourth independent design

DPY30 sits at the 14th percentile with a difference of -0.063 and does not survive correction. It has now failed the variance-matched null, the deep-pool null, the within-subtype test and this one, while passing only the permutation test that asks whether its correlation differs from zero.

Four designs agreeing is worth more than any one of them. DPY30 should be treated as a gene whose correlation with MEN1 is real and unremarkable, and entries resting on it being the distinctive WRAD subunit need revision.

Robustness: dropping the FLT3-ITD lines, which partly corrects the caveat above

The section above warned that MV4-11, MOLM-13 and MOLM-14 are FLT3-ITD lines and that FLT3's near-equal difference is therefore probably a passenger of the group rather than of the fusion. That caveat was worth testing rather than merely stating, so it was. Excluding those three leaves 8 rearranged lines — SEMK2, HB11;19, RS4;11, SEM, THP-1, NOMO-1, OCI-AML-2, SHI-1 — against the same 43.

Genedifference, all 11difference, 8 non-ITDpermutation p
MEN1-0.866-0.7500.0001
FLT3-0.843-0.6210.0001
MEF2C-0.560-0.5740.0147
PSIP1 (LEDGF)-0.565-0.5260.0002
KMT2A-0.541-0.4620.0023
BCL2-0.396-0.3710.0044
DOT1L-0.231-0.2100.177
ACTB (control)-0.182-0.1940.191

MEN1 survives cleanly. It keeps 87% of its difference on eight lines at p = 0.0001, so the FLT3-ITD lines were not carrying the result.

The caveat about FLT3 was half right and should be corrected rather than left standing. FLT3's difference does fall, from -0.843 to -0.621, so the ITD lines contribute. But it stays highly significant without them, which means FLT3 is genuinely more required in KMT2A-rearranged lines independent of ITD status. Writing it off as an ITD artefact, as the section above did, was too quick.

ACTB is reassuring here. On the permutation test it is not significant either with the ITD lines (p = 0.16) or without (p = 0.19). The global-offset worry rests on its position in the random-gene null, not on the within-panel test, and it remains a reason to treat magnitudes as upper bounds rather than a reason to doubt the ranking.

DOT1L does not survive this. At p = 0.177 it is the one candidate whose group difference is not distinguishable from relabelling, consistent with it sitting lowest of the Holm survivors in the main table. DOT1L's case rests on correlation and on mechanism, not on differential dependency between these groups.

The global-offset explanation was mine, it was testable, and it is wrong

The section above proposed that the ACTB control's position at the 1.4th percentile reflected a mild global offset — rearranged lines being more dependency-sensitive across the board, inflating every difference in the table. That was a mechanism, not an observation, so it was run rather than left as a caveat (analyses/depmap-line-normalised-specificity/).

Each line's own sensitivity was estimated as its mean gene effect across the random panel only — excluding the candidates, so they cannot help set the baseline they are then measured against — and subtracted from every value before the comparison was rerun, with the null rebuilt on the normalised statistic.

There is no offset to remove. Mean per-line sensitivity is -0.1551 in the rearranged lines and -0.1515 in the others: a difference of -0.0036, which is negligible beside differences of 0.3 to 0.9.

Generaw differencenormalised difference% retainedpercentile
MEN1-0.866-0.86299.60.0
FLT3-0.843-0.83999.60.0
PSIP1-0.565-0.56199.40.0
MEF2C-0.560-0.55799.40.0
KMT2A-0.541-0.53899.30.0
BCL2-0.396-0.39299.10.0
ACTB (control)-0.182-0.17898.01.28
DPY30-0.063-0.05994.214.3

Two conclusions, and the second is against me.

First, the main result is not an artefact of differential screen sensitivity. MEN1 keeps 99.6% of its difference. That strengthens the finding and the caveat it was carrying can be discharged: magnitudes no longer need to be read as upper bounds on this account.

Second, my explanation for the ACTB anomaly was wrong. ACTB keeps 98% of its difference and stays at the 1.3rd percentile. Whatever puts it there, it is not a global offset. The better reading is the one the percentile obscured: ACTB's permutation p is 0.169, so it is not significant on its own test, and a single gene's rank in a 1,097-gene null is a noisy statistic — roughly 14 random genes sit below it. Ranking one control at the 1.3rd percentile out of thirteen genes tested is unremarkable.

So the anomaly was largely an artefact of how I read the percentile, not of the data. The percentile column is still the right one for the headline genes, which sit beyond every random gene; it is the wrong tool for adjudicating a single borderline control.

What this does not establish

Eleven lines. Every conclusion rests on a group of 11, and 7 in the AML-only test. That is small, and the random-gene null is what makes the comparison interpretable at all.

Cellosaurus annotation is curated from publications and may be incomplete. Any unannotated rearranged line sitting in the 'other' group weakens the contrast, so this bias is conservative — the true difference is if anything larger. But it also means the 'other' group is not cleanly fusion-negative.

Co-occurring lesions are not separated from the fusion, as FLT3 shows plainly.

Cell lines are not infant leukemia. None of these 11 is an infant-derived line, and the population this project exists to serve is infants. This is the same gap every DepMap analysis here has.

Dependency is not druggability. A gene more required in rearranged lines may still have no ligandable pocket and no therapeutic window.

Facts and figures
scriptanalyses/depmap-kmt2a-rearranged-specificity/run.py
statusrun
run date2026-09-18
answers open question fromana-depmap-within-subtype
methodCellosaurus var:KMT2A AND dr:DepMap, filtered to variation-type 'Gene fusion' with KMT2A among the HGNC-linked genes; difference in mean gene effect between rearranged and other leukemia lines; 10,000-iteration label permutation plus a 918-gene random null of the same statistic; repeated within AML only; Holm over eleven candidates
seed20260918
n kmt2a rearranged11
n other leukemia43
n kmt2a rearranged aml7
n other aml25
n random usable918
figureanalyses/depmap-kmt2a-rearranged-specificity/results.md
rearranged lines
[
  "SEMK2",
  "HB11;19",
  "RS4;11",
  "SEM",
  "THP-1",
  "NOMO-1",
  "OCI-AML-2",
  "MV4-11",
  "MOLM-13",
  "SHI-1",
  "MOLM-14"
]
kmt2a aff1 lines
[
  "SEMK2",
  "RS4;11",
  "SEM",
  "MV4-11"
]
beats all 918 random genes
[
  "MEN1",
  "FLT3",
  "PSIP1",
  "MEF2C",
  "KMT2A",
  "BCL2"
]
survives holm
[
  "MEN1",
  "FLT3",
  "PSIP1",
  "MEF2C",
  "KMT2A",
  "BCL2",
  "MEIS1",
  "KAT7",
  "DOT1L"
]
does not survive holm
[
  "MLLT1",
  "DPY30"
]
control percentiles
{
  "ACTB": 1.42,
  "GAPDH": 85.84
}
robustness no flt3 itd
{
  "excluded": [
    "MV4-11",
    "MOLM-13",
    "MOLM-14"
  ],
  "n_rearranged": 8,
  "MEN1": {
    "difference": -0.7502,
    "perm_p": 0.0001
  },
  "FLT3": {
    "difference": -0.621,
    "perm_p": 0.0001
  },
  "PSIP1": {
    "difference": -0.526,
    "perm_p": 0.0002
  },
  "KMT2A": {
    "difference": -0.462,
    "perm_p": 0.0023
  },
  "DOT1L": {
    "difference": -0.21,
    "perm_p": 0.17748
  },
  "ACTB": {
    "difference": -0.194,
    "perm_p": 0.19108
  }
}
line normalisation control
{
  "script": "analyses/depmap-line-normalised-specificity/run.py",
  "status": "run",
  "run_date": "2026-09-18",
  "mean_offset_kmt2a_r": -0.1551,
  "mean_offset_other": -0.1515,
  "offset_difference": -0.0036,
  "MEN1_pct_retained": 99.6,
  "ACTB_pct_retained": 98,
  "verdict": "no global offset exists between the groups; the main result is not an artefact of screen sensitivity, and the proposed explanation for the ACTB anomaly is refuted"
}
How this could bite us back
  • This is the result the project most wanted and the one it would be most damaging to over-read. It shows MEN1 dependency tracks the rearrangement in 11 cell lines, and survives dropping the FLT3-ITD lines. It still does not show the fusion *causes* that dependency: the rearranged lines differ from the others in more than one way, and the design cannot separate the fusion from what travels with it.
  • A single control behaving badly is enough to require the magnitudes be read as upper bounds. ACTB at the 1.4th percentile is that control.
  • Eleven lines, none of them infant-derived, standing in for a disease of infants. Every DepMap analysis in this project shares that gap and none of them fixes it.
  • Dependency is not druggability. Nothing here says any of these proteins has a ligandable pocket or a therapeutic window, and menin's own window is challenged elsewhere in this graph on endocrine grounds.
  • A caveat that can be tested should be tested, not carried. This entry shipped with a global-offset caveat that turned out to describe nothing, and carrying it would have understated a real result while leaving a wrong explanation on the record.
Still open
  1. ANSWERED within this entry: excluding the three FLT3-ITD lines, MEN1 keeps 87% of its difference at p = 0.0001. The remaining question is narrower — FLT3 itself stays significant without them, so what makes FLT3 required in non-ITD KMT2A-rearranged lines?
  2. ANSWERED within this entry: there is no global sensitivity offset between the groups (-0.0036), so the main result is not an artefact of it and my proposed explanation for the ACTB anomaly was wrong. ACTB is simply not significant on its own permutation test (p = 0.169).
  3. Five of the 11 lines carry KMT2A::AFF1 specifically. Is the menin dependency equally strong across partners, or does it differ between AFF1, MLLT3 and AFDN fusions? n is too small here but the question is now askable.
  4. Cellosaurus annotates fusions for only a subset of DepMap lines. How many of the 43 'other' leukemia lines are unannotated rather than genuinely fusion-negative, and can Cell Model Passports or Cosmic fill the gap?
  5. DPY30 has now failed four independent designs. Which entries still rest on it being the distinctive WRAD subunit, and should they be rewritten or retired?
Every source behind this entry (1)
  1. E1moderate confidence 3 of 5: One strong study, or several smaller consistent ones.dataset

    Among 54 DepMap leukemia lines with MEN1 gene-effect data, the 11 annotated by Cellosaurus as carrying a KMT2A gene fusion have mean MEN1 gene effect -1.263 versus -0.397 in the other 43, a difference of -0.866 that is more extreme than all 918 random protein-coding genes tested (permutation p = 0.0002; Holm-corrected empirical p = 0.012), and that persists within AML alone (7 versus 25 lines, difference -0.890, p = 0.0011). FLT3, PSIP1, MEF2C, KMT2A and BCL2 also exceed all 918 random genes; MEIS1, KAT7 and DOT1L survive Holm correction; MLLT1 and DPY30 do not. Control genes ACTB and GAPDH sit at the 1.4th and 85.8th percentiles.

    What this evidence cannot show

    Our own computation, not peer reviewed. Only 11 rearranged lines, 7 in the AML-only test. Cellosaurus fusion annotation is curated from publications and may be incomplete, so the 'other' group may contain unannotated rearranged lines — a conservative bias. Co-occurring lesions are not separated from the rearrangement, and FLT3's near-equal difference is almost certainly driven by FLT3-ITD in MV4-11 and the MOLM lines rather than by the fusion. The ACTB control sits at the 1.4th percentile of the random-gene null but is not significant on its own permutation test (p = 0.169); per-line normalisation shows no global sensitivity offset between the groups (-0.0036), so magnitudes do NOT need discounting on that account. No line is infant-derived. DepMap release version is not stated by the API. **The eleven lines are nine individuals, and the index-fusion arm is three people.** Verified by the coordinator against Cellosaurus by accession on 2026-09-20: **SEMK2 (CVCL_S906) is a child record of SEM (CVCL_0095)** — one 5-year-old patient, the same relationship MOLM-13/MOLM-14 already carry. Of the nine individuals only three carry KMT2A::AFF1, and they are **RS4;11 (CVCL_0093, age 32Y, Cellosaurus disease term 'Adult B acute lymphoblastic leukemia')**, **MV4-11 (CVCL_0064, 10Y male, 'Childhood acute monocytic leukemia')** and that one 5-year-old counted twice. Any per-line statistic computed over the eleven over-counts.

    Evidence that argues against this (3)
    • The ACTB control should sit in the middle of the null and sits at the 1.4th percentile, so the rearranged lines appear mildly more dependency-sensitive overall. Some fraction of every difference in this table is therefore a global offset rather than biology. GAPDH going the opposite way (85.8th percentile) argues the offset is small, but two controls disagreeing is a reason for caution, not reassurance. analysis run by this project 2026
    • Correlation-based analyses in this project reached a partly different ordering: DOT1L ranks third by correlation with MEN1 but ninth here, and BCL2 ranks low by correlation but sixth here. The two statistics measure different things — covariation across lines versus difference in mean dependency between groups — and a gene can score well on one and poorly on the other. Neither ordering should be treated as the target ranking. analysis run by this project 2026
    • Effective n is the number this analysis should report and does not. Nine patients for the rearranged arm and three for the index fusion, one of them an adult and one a paediatric monocytic AML, is a different claim from eleven lines - and the confidence intervals, the lineage split and the line-normalisation all assume independence that SEM and SEMK2 do not have. SIB Swiss Institute of Bioinformatics 2026
Connected to