Built with Claude · Life Sciences

From a 22-million-cell CD4⁺ T-cell screen to human-genetics-validated drug targets

About 90% of drugs fail on the wrong target. We built a pipeline that makes a cell screen prove itself against real patient genomes — and it works: it re-finds known biology (including IL23R→IBD, the right gene in the right disease) and returns honest "no" answers for the novel lupus hypotheses. Positive controls fire; the pipeline can say "no".

We take every gene from a huge CD4⁺ T-cell screen and ask one simple question for each: does human genetics agree that changing this gene would help a patient? Each gene goes through five checks — cell biology → can a drug hit it → disease genetics → large biobanks → real patient genomes from All of Us. Only genes that pass every check are trusted. The result: the method re-finds known biology (IL23R → IBD), gives honest "no" answers for the lupus ideas, and flags one gene that it then corrects itself on.

IL23R → IBD validated · OR 0.73, P 4.4×10⁻⁵ ATP2B1 surfaced · then down-weighted (cardiometabolic) KDM8 exciting lead · deflated by our own nearest-gene test reusable on any perturbation screen
Dataset: Genome-scale CD4⁺ T-cell Perturb-seq — 22 M cells, 4 donors, CRISPRi knockdown across three activation states (Rest / Stim 8 h / Stim 48 h). Marson Lab, Gladstone Institutes, via the CZI Virtual Cells Platform (Dann et al. 2025).  ·  Built end-to-end in Claude Science — every figure, table and stat is a reproducible artifact.

The result in one look

The funnel, with the actual counts: ~4,000 genes narrow to a handful, and only what clears the final patient-genome gate is trusted. The verdicts below are what came out the other end. (How each gate works is detailed further down.)

1
Cell biology
Switching the gene off changes the T‑cell activation program
2,771 candidate genes
→
2
Druggable
A drug can physically reach the protein
86 ranked targets
→
3
Disease genetics
Linked to autoimmune disease in human GWAS
shortlist of receptors
→
4
Biobanks
The link repeats in independent national biobanks
FinnGen + GWAS Catalog
→
5
Real patient genomes
Tested in All of Us individual-level DNA
the final verdict
IL23R → IBDpassed all 5 gates · a genuine win
IL12RB2 · IL7Rdirectly tested in All of Us · honest null. IL21R null via a common-variant proxy (its lead SNP wasn't in the cohort).
ATP2B1flagged, then down-weighted · self-correction

Why this matters

About 90% of new drugs fail in clinical trials — most often because the target gene was wrong for the disease. Targets that also have human genetic support are about twice as likely to succeed. A screen shows what a gene does inside a cell. It does not tell you whether changing that gene would help a real patient. That is the gap we close.

The gap we close

The screen gives thousands of candidate genes. We link each one to human disease genetics — including real patient genomes from All of Us. A gene is trusted only when the cell biology and the human genetics point the same way. That overlap is where good drug targets are.

Who it's for

Researchers working on autoimmune disease — lupus, IBD, rheumatoid arthritis, MS — and anyone with a gene screen asking "which of these hits is worth turning into a drug?" The answer here rests on five different kinds of evidence, not one experiment.

Why it's credible

The method re-finds biology we already know is true (PTPN22→RA, TYK2→RA, IL23R→IBD). It also says a clear "no" to the lupus receptor ideas when the data don't support them. A method that can say "no" is one you can trust when it says "yes". Every number can be reproduced in Claude Science.

The dataset

Our starting point is one of the largest single-cell experiments ever run in human immune cells — the genome-scale CD4⁺ T-cell Perturb-seq atlas from the Marson Lab at Gladstone Institutes, shared through the CZI Virtual Cells Platform (Dann et al., 2025).

22 M
single cells measured
~2,000
genes switched off, one at a time
4
human blood donors
3
T-cell activation states

What is in the data — in plain words

CD4⁺ T cells are the "commander" cells of the immune system. When they behave badly, they drive autoimmune diseases like lupus, IBD, rheumatoid arthritis and MS. So they are exactly the right cells to look in for new autoimmune drug targets.

Perturb-seq means: switch off one gene at a time (here using CRISPRi, which quiets a gene instead of cutting it), then read the full gene-activity profile of each single cell. Do this for almost every gene in the genome, and you learn what each gene actually does to the T cell.

Why three activation states matter

The same gene can matter a lot or not at all depending on whether the T cell is resting or activated. The atlas captures three states — Rest, Stim 8 h, and Stim 48 h (8 and 48 hours after activation) — so a target's effect is measured in the state where it is actually used.

We rank genes within each state separately and never average across them. A gene that only acts in activated cells is not hidden by its "quiet" resting behaviour.

We do not re-process the 22 M raw cells. We build on the study's own published supplementary tables — the distilled signal from all cells and all genes — and add the human-genetics layers on top.

What we add beyond the original screen

The Marson dataset already links these clusters to autoimmune genetics. Our contribution is the translational layer on top — a per-gene decision framework, not a repeat of the paper's enrichment analysis. The genetics layer, not the bench score, is what re-orders the list.

Original Perturb-seq paper provides
  • Open Targets autoimmune genes
  • Regulatory clusters & downstream programs
  • State-specific perturbation effects
  • Autoimmune-GWAS enrichment (cluster-level)
What THIS pipeline adds (translational layer)
  • Therapeutic DIRECTION call (inhibit vs activate)
  • Multi-biobank replication (FinnGen + GWAS Catalog)
  • Cross-ancestry check (BioBank Japan, East Asian)
  • All of Us individual-level rare-variant test
  • Druggability + clinical stage + failure mode
  • Novelty-gap discovery + reproducible auto-audit
Ablation: dropping genetics reshuffles the list most
− genetics0.25
− druggability0.82
− bench0.54
top-10 overlap with full ranking (Jaccard)
7 / 20
of our top-20 targets are not in the paper's autoimmune-intersecting gene set
ρ = 0.55
Spearman of our composite vs the pure-bench rank — genetics + druggability substantially re-order it
25%
top-10 overlap after removing the genetics layer — genetics is the decisive re-ranker

A translational decision layer, quantified. 7 of our top-20 targets are absent from the paper's autoimmune-intersecting gene set; our composite correlates only ρ=0.55 with the pure-bench rank; and removing the genetics layer collapses top-10 overlap to 0.25 — genetics is the decisive re-ranker.

How each gate works

The scoreboard above shows what passed; here is the method behind each gate. Every layer is a different kind of evidence, drawn from a different data source — a gene earns trust only by clearing checks that ask genuinely different questions. Cell biology and human genetics are strongest when they agree.

LAYER 1

Bench causality

We switch each gene off (CRISPRi) and measure how the T-cell activation program changes. Genes are ranked in each condition and must repeat.

→
LAYER 2

Druggability

How easy is the gene to hit with a drug? Cell-surface receptors and kinases score higher than harder classes like transporters.

→
LAYER 3

Human GWAS

Is the gene linked to disease in genome-wide studies? We use Open Targets and the GWAS Catalog, and flag known traps (MHC, nearest-gene).

→
LAYER 4

Biobank PheWAS

Does the link repeat in large national biobanks? We check FinnGen and the GWAS Catalog across 14 autoimmune diseases.

→
LAYER 5

All of Us

The strongest test: real patient genomes (All of Us). We test both common and rare gene variants — something no public summary table can give you.

The targets, in 3D

Live AlphaFold structures — rotate, zoom, explore. The 3D model is colored by AlphaFold per-residue confidence (blue = high, orange = low). Each target's pipeline verdict is shown by its tab colour and in the card: green = validated, grey = tested null, orange = flagged.

The All of Us reveal

This was never a lupus-only search. Every target is tested against four autoimmune diseases (lupus, IBD, rheumatoid arthritis, MS) in individual-level data — and screened across 14 in public biobanks. Two known controls (PTPN22, TYK2) come out positive — proof the method works on real patient data. IL23R is confirmed in the right disease (IBD). The three novel receptor ideas come back as a clean, well-checked "no" across all four — lupus included, but not lupus alone.

All of Us validation scorecard
Individual-level validation scorecard. Method controls (blue) confirm the pipeline; IL23R→IBD (green) is a genuine bench→biobank win; the three novel receptors are tested-null (grey). Right: the same IL-12/23 receptor family splits — IL23R validates while its sibling IL12RB2 is null — and coherent rare-variant burden trends.

The moat: discovery and self-correction

We flipped the usual search to find genes with strong genetics but almost no published papers, then tested them where only patient-level data can — rare variants. The result is honest in both directions.

Rare-variant burden and hidden-hit verdict
Hidden-hit burden. The flagship KDM8 (0 autoimmune papers) comes back null — a moat that can say "no". ATP2B1 is the one gene to survive Bonferroni (SLE P 1×10⁻⁶, RA replicates).
ATP2B1 phenome specificity
…then down-weighted, honestly. ATP2B1's genome-wide-significant phenome is dominated by cardiometabolic disease (hypertension P 3×10⁻²⁴; 0 autoimmune phenotypes reach significance). We flag it as a lead needing blood-pressure-conditioned analysis, not an autoimmune target claim.
Round-3 hardened rare-variant test
Hardened rare-variant test. Splitting the mask shows ATP2B1’s protection is a damaging-missense effect (1,469 carriers, OR 0.73, P 9×10⁻⁷), not pLoF (31 carriers, null) — a mechanism a merged mask would hide. Firth regression handles the small IBD/MS cohorts, and for 6 of 8 genes the rare-variant direction agrees with the CRISPRi knockdown sign.

Where each target lands

Instead of collapsing everything into one score, we sort candidates into four honest tiers using a therapeutic-window layer (gnomAD LOEUF), clinical stage, and causal-gene status. This separates clinically supported targets from novel-but-unresolved hypotheses — and flags the ones with strong biology but a poor therapeutic window.

Clinically supported
IL7R
LOEUF 1.06 (tolerant) · OT 0.59 · Phase 2 UC + Phase 2a Sjögren
anti-IL7R in randomized autoimmune trials — framework re-finds an emerging validated target
Mechanistically convergent
IL21R
LOEUF 0.53 (constrained) · FinnGen 49.3 · ATR-107 (immunogenic)
strong bench + genetics; prior antibody failed on immunogenicity, not biology — re-engineerable
IL12RB2
LOEUF 1.10 (tolerant) · FinnGen 32.7 · ustekinumab (pathway)
IL-12/23 pathway clinically validated; direct receptor-subunit targeting still unvalidated
Novel — unresolved
KDM8
LOEUF n/a · FinnGen 35.9 → maps to IL4R · AoU burden null
high-priority hypothesis pending causal-gene resolution — NOT a confirmed target
BCAT2
LOEUF 1.07 · FinnGen 11.5 → maps to FUT2 · no CD4 eQTL
strong bench (z=−91) but disease signal is a neighbor; needs coloc
Strong biology · poor therapeutic window
ZAP70
LOEUF 0.57 (constrained) · core TCR kinase
top bench score, but pan-T essentiality → broad immunosuppression risk
JAK2
LOEUF 0.74 · approved (ruxolitinib/baricitinib)
validated drug class but broad cytokine effects; positive control, not novel

Tiered target cards. IL7R is clinically advancing (Phase 2); IL21R/IL12RB2 are mechanistically convergent; KDM8/BCAT2 are novel-unresolved (pending causal-gene resolution); ZAP70/JAK2 have strong biology but a poor window (constrained LOEUF, pan-T essentiality).

Closing the obvious critiques

Five common objections, answered with data across the six panels below: arbitrary weighting (weight-sensitivity), nearest-gene (answered twice — a stress test and formal GWAS×eQTL colocalization), single-disease (cross-disease pleiotropy), screen-replication (an orthogonal CRISPR screen), and Europe-only (cross-ancestry BioBank Japan).

Composite-weight sensitivity
The ranking is not an artifact of hand-tuned weights. Across 64 composite-weight combinations, IL21R, IL12RB2 and IL7R stay in the top-10 100% of the time — the headline receptors are robust to any reasonable weighting.
Hidden-hit nearest-gene stress test
The nearest-gene guard, turned on our own leads. The four most exciting "hidden hits" (strong public genetics, almost no autoimmune literature) each fail when we ask which gene the lead variant actually maps to: the signal belongs to a well-known autoimmune gene next door — IL4R, GSDMA, FUT2, GCKR — not the screen hit. Our flagship KDM8 (0 autoimmune papers) is really the IL4R allergy locus, and it is also null in All of Us rare-variant burden. A pipeline that deflates its own best-looking leads is one you can trust when it says "yes".
GWAS x eQTL colocalization
GWAS × cis-eQTL colocalization. Gene-specific verdicts: IL21R's disease lead sits near its eQTL peak (partial overlap); IL12RB2's disease signal is 168 kb from the eQTL peak — distinct signals, so the association is not a simple cis-eQTL of the named gene.
Cross-disease pleiotropy
Cross-disease pleiotropy. Perturbation-response programs shared across 14 autoimmune diseases — the shared nodes are the most attractive targets: one mechanism, several indications.
Orthogonal CRISPR-screen replication
Orthogonal bench replication. The TCR-proximal core (ZAP70, ITK, PLCG1) replicates strongly in an independent genome-scale CRISPR screen (Schmidt 2022; lfc −2.8/−1.5/−1.0, FDR<0.001) — a method sanity check. The cytokine receptors are read-out-orthogonal (invisible to a single-cytokine fitness screen), which is the point.
Cross-ancestry replication in BioBank Japan
Cross-ancestry replication. Not just European: for IL7R and IL21R the autoimmune/allergy signal reaches genome-wide significance in the multi-ancestry (BBJ+UKB) meta (IL7R p=5×10⁻²³, IL21R p=2×10⁻¹²), with the East-Asian-only stratum concordant in direction but sub-significant on its own (IL21R BBJ p=3×10⁻⁶). Crucially the effect directions agree across populations (panel b, on the y=x line); IL12RB2’s RA signal is genome-wide significant in the multi-ancestry meta. IL23R R381Q is ~7× rarer in East Asians (its classic frequency story).

Run it on your own data

The pipeline is a config-driven package — bench2biobank. Point it at any perturbation screen (CRISPRi/a, ORF, shRNA, drug; any cell type), map your column names, and get a ranked, druggability-weighted, disease-linked target shortlist plus a genome-scale enrichment test.

# 1. describe YOUR data in a config — just map the column names
perturbation_table:
  path: my_screen/effects.csv
  columns:
    perturbation: gene        # your header for the perturbed gene
    effect: logFC            # signed effect on the readout
    stderr: se
    signature: readout
    condition: timepoint
signature_of_interest: {name: activation, conditions: [8h, 48h], direction: lower}

# 2. run it
$ bench2biobank --config my_config.yaml

[bench2biobank] 3994 perturbations → 2771 candidates → 86 ranked targets
[headline] top-9% effect × disease-linked enrichment: OR=1.47, p=2.36e-03
  ZAP70   composite=0.631   IL21R  composite=0.593   BCAT2  composite=0.546 ...
Disease and druggability layers are optional and degrade gracefully — with neither, you still get a clean per-condition bench ranking. A schema check tells a new user exactly which columns their data needs.