The Nineteen
## What we built / investigated We mined a genome-scale CRISPRi Perturb-seq atlas of about 22 million primary human CD4⁺ T cells to discover **novel drug targets**.
Core trick
Reads a genome-scale CRISPRi Perturb-seq atlas backward from desired T-cell programs to upstream regulators, then validates candidates across donors at single-cell level.
Description
## What we built / investigated
We mined a genome-scale CRISPRi Perturb-seq atlas of about 22 million primary human CD4⁺ T cells to discover **novel drug targets**. (Per the project's dataset documentation, the atlas is the CD4⁺ T-cell Perturb-seq resource from the Marson and Pritchard labs, distributed on the CZI Virtual Cells Platform; this attribution is as recorded in the project files and has not been independently re-verified here.) The core idea is that a CRISPRi knockdown is a genetic model of drug-induced loss of function, so we read the perturbation-to-transcriptome map backwards: starting from a therapeutically desirable T-cell program and asking which upstream regulators control it. Those regulators become drug-target hypotheses. On top of this we built an end-to-end computational-to-translational pipeline: foundation QC and reproducibility filtering, then six orthogonal discovery lenses, then an integrated ranked scorecard, then single-cell validation across donors, and finally structure, genetics, and chemistry-based translational dossiers.
## What we found
The pipeline narrowed roughly 34,000 perturbations down to a short list of **19 single-cell-validated targets**, and this is the main result (full detail in the Appendix table). It splits into two therapeutic groups.
**A TCR-proximal signaling group** (LAT, PLCG1, CD247, CD3E, ZAP70, VAV1, SIK3) contains the highest-ranked hits and mostly points toward suppressing inflammation. It recovers known immune-drug biology on its own, which validates the logic, and it surfaces **SIK3**, a novel kinase that became our structurally validated lead (redocking under 2 Å, actives-versus-decoys enrichment AUC 0.905).
**A drug-naive chromatin and transcriptional group** (a SAGA/Mediator axis plus remodelers: SMARCE1, STAT6, SGF29, MED24, MED12, TADA2B, NSD1, CHD4, SMARCB1, TRIP12, ARNT, SEL1L) is the central novelty. All 12 of these validate at single-cell resolution and form a mechanism distinct from acute TCR signaling (89% cluster recovery, 2.4× smaller state-space displacement). Of the full set of 19, **18 are clinically unprecedented** (only CD3E is already drugged).
The single top opportunity is **STAT6** (rank 5): a boost-immunity target backed by strong asthma and allergy genetics (55 immune GWAS associations), rich chemical matter (552 ChEMBL bioactivities, max pChEMBL 9.15), and small-molecule tractability, yet with no approved drug. The remaining standouts are the drug-naive SAGA/Mediator readers **MED24**, **SGF29**, and **TADA2B**. Two targets (CHD4, SMARCB1) are flagged for essentiality and toxicity risk on weaker knockdown. The complete ranked list, with validation and translational evidence for each target, is given in the Appendix.
## Why it matters
By recovering known targets from scratch, which validates the logic, and then nominating a genetically-supported, single-cell-validated, structurally-tractable set of *novel* regulators, the project turns a 22-million-cell atlas into a short list of de-risked, mechanistically-distinct therapeutic entry points for autoimmune and allergic disease and for immuno-oncology. STAT6 and a previously undrugged chromatin/transcriptional axis head that list. More broadly, the work demonstrates a reusable in-silico funnel that goes from genome-scale perturbation data to prioritized, translation-ready drug targets.
MengQiuchen/perturb-seq-cd4-target-discovery★ 00 forksPythonMITCreated Jul 13, 2026 · last push Jul 13, 2026
README excerpt
# CD4⁺ T-Cell Perturb-seq Drug-Target Discovery Discovery and translational prioritization of **novel druggable regulators of CD4⁺ T-cell programs** from a genome-scale CRISPRi Perturb-seq atlas. A CRISPRi knockdown is read as a genetic model of drug-induced loss of function; this pipeline mines the perturbation→transcriptome map **backwards** — from therapeutically desirable T-cell programs to the upstream regulators that control them — and carries the resulting candidates from genome-scale nomination through single-cell validation to translational dossiers. > **Scope note.** This repository contains the **analysis code, curated result tables, figures, and reports** produced on top of a public dataset. It does **not** redistribute the raw/processed expression matrices (hundreds of GB to ~2.3 TB); those are fetched from the public source (see [Data access](#data-access)). --- ## Key results - **Foundation (Phase A).** From the precomputed DESeq2 layer (33,983 perturbation×context contrasts × 10,282 genes) a reproducible, on-target, off-target-clean **high-confidence set of 12,576 perturbations / 5,728 genes** was defined. Known drug-target biology reproduces: **13/16 positive controls** move their canonical readouts in the expected direction. - **Discovery + integration (Phases B–C).** Four orthogonal functional lenses + two external evidence layers (human genetics, druggability) converged on **301 candidate genes** (≥2 directions; 61 at ≥3; LAT/PLCG1/SENP5 at all 4). **