Built-with-Claude · Life Sciences · Researcher track

Selective Th2 Suppressor Atlas

Knockdowns that take the Th2 program down without flipping CD4+ T cells toward Th1 — the allergy-relevant target the paper’s single Th2−Th1 axis structurally can’t resolve.

+0.54
Arm correlation
confound = global magnitude
127
Selective (2D gate)
Th2↓, Th1 flat
10
True Th1-skewers
rare — confirms selectivity
21
High-confidence genes
KD-confirmed + footprint
Calibration. GATA3 (Th2 master regulator) is the #1 selective hit and ranks #1 in the scored shortlist. The known-Th2 list is directionally mixed: activators (GATA3, RARA, STAT6) recover as suppressors, while NuRD/chromatin members (MTA2, CHD4, SETDB1) appear high on the Th1 arm — though literature review shows CHD4/SETDB1 are bimodal, so this is a soft calibration, not a clean dichotomy.
Honest status (read this first). This is a rigorous re-analysis + reasoning layer, not a novel method. Independent review (below) flagged real limitations:

The 2D arm plane

headline result
Th2 vs Th1 arm plane
Every point is one perturbation×condition. Green = selective Th2 suppressor (Th2 arm down, Th1 flat); red = Th1-skewer. The cloud tilts along the diagonal (positive arm correlation) — the dominant confound is global dampening, which the Th1 arm controls for.

Context-dependence

the novel cut

Selectivity is strongly stimulation-dependent — something the paper’s per-condition 1D coefficients can’t express as selectivity. Most selective suppressors act in only one context (GATA3 itself is selective at Rest/Stim8hr but not Stim48hr).

Rest-only
60
Stim-induced
49
Mixed
7
Constitutive
1

Scored shortlist

real vs false-positive
GeneScoreTh2 armTh2 residContextClassDonor✓
GATA3 known0.93-3.70-3.52Stim8hrTF✓
AHSA1 novel0.90-2.11-2.53Restother✓
ARNT novel0.83-2.51-2.59Stim48hrTF✓
RAB21 novel0.83-2.76-2.94Restother·
ELAVL1 novel0.81-2.68-2.35Stim48hrother✓
GPATCH8 novel0.79-2.18-1.94Restother✓
HNRNPA2B1 novel0.77-2.31-2.42Restother✓
STK35 noveldruggable0.77-2.11-1.82Stim8hrkinase✓
UFC1 novel0.77-2.41-2.15Restother✓
ARL1 novel0.76-2.29-1.84Stim48hrother✓
MIA2 novel0.74-2.01-1.71Restother✓
KDM2A novel0.72-2.22-1.72Stim48hrTF✓
SSR1 novel0.72-2.58-2.13Restother·
RARA known0.71-2.04-2.54Stim8hrTF·
AMBRA1 novel0.70-2.57-2.23Stim48hrother✓
STAT5A novel0.70-2.29-2.00Stim8hrTF·
COMMD1 novel0.70-2.22-2.00Stim48hrother·
IL2RA novel0.69-2.11-1.89Stim8hrother·
IFT172 novel0.66-2.06-1.70Stim48hrother✓
SRD5A3 novel0.65-2.02-2.41Restother✓
MPG novel0.64-2.21-2.03Stim48hrother✓

Th2 resid = Th2 suppression after regressing out the Th1 (global-magnitude) arm. Donor✓ = cross-donor reproducibility confirmed.

Mechanism dossiers

agentic layer · investigate → verify

Adversarial review of the top 10 candidates kept AHSA1, ARNT, ELAVL1, GPATCH8 as the non-artifact leads worth following up; 4 were flagged as likely global-effect false positives. The mechanism layer is doing the triage the numeric score can’t — surfacing which high-scoring hits are real Th2 biology vs. sick-cell artifacts.

AHSA1uncertain
Loss of AHSA1 (Aha1), the sole mammalian activator of HSP90 ATPase, impairs HSP90 chaperoning of the metastable Th2-defining clients STAT6 and GATA3, collapsing the GATA3/IL4/IL5/IL13 program while leaving the less chaperone-addicted T-bet/IFNG axis intact.
cites: partially_supported
ARNTuncertain
ARNT (HIF-1β) is the obligate dimerization partner for HIF-1α/HIF-2α, and because stimulation-induced HIF signaling promotes the GATA3/IL-4 Th2 program while Th1/T-bet is relatively HIF-independent, ARNT knockdown collapses HIF output and selectively lowers Th2 without inducing Th1.
cites: partially_supported
ELAVL1uncertain
Knockdown of ELAVL1 (HuR) removes AU-rich-element-mediated stabilization of GATA3, IL4, and IL13 mRNAs, collapsing the Th2 program post-transcriptionally without needing to induce the Th1/T-bet axis.
cites: partially_supported
GPATCH8uncertain
As a branchpoint splicing quality-control factor, GPATCH8 loss could selectively impair productive splicing/expression of a Th2-program transcript (e.g., GATA3 or an upstream STAT6/IL4-axis component), lowering Th2 output without engaging the T-bet/Th1 axis.
cites: supported
MPGunlikely
Loss of MPG (AAG) may selectively blunt the Th2 program by disrupting base-excision-repair/5hmC-reader activity at the demethylation-dependent IL4/IL5/IL13 Th2 cytokine locus, rather than by driving Th1.
cites: supported
STK35unlikely
Losing STK35 could lower the Th2 program indirectly by blunting the AKT/glycolysis-driven proliferative and metabolic burst that GATA3/IL4 induction disproportionately depends on, rather than by acting as a dedicated Th2 transcriptional regulator.
druggablecites: partially_supported
KDM2Alikely artifact
Because KDM2A is tightly coupled to PI3K/mTOR/SGK1-driven rRNA transcription, ribosome biogenesis and proliferation, its knockdown most likely dampens the highly proliferation/autocrine-amplified Th2 program more than the low-baseline Th1 program in a Th2-polarizing assay, rather than acting as a bona fide selective GATA3/IL4 regulator.
cites: supported
MIA2likely artifact
If MIA2 knockdown selectively lowered the Th2 program, it would most plausibly be by degrading the ER-export/secretory capacity that autocrine IL-4/STAT6/GATA3 feedback (and the high secretory load of Th2 cells) depends on, rather than by any Th2-specific transcriptional role.
cites: partially_supported
RAB21likely artifact
Loss of RAB21 could preferentially blunt the highly glycolysis/mTOR-dependent Th2 (GATA3/IL4) program by mis-sorting surface GLUT1 and impairing receptor recycling, without inducing T-bet/Th1.
cites: supported
SSR1likely artifact
Loss of SSR1 impairs ER translocation/secretion of the heavily secreted Th2 cytokines and could collapse the autocrine IL-4→STAT6→GATA3 positive-feedback loop, indirectly lowering the Th2 mRNA program more than the less autocrine-dependent Th1 program.
cites: partially_supported

Known-regulator recall

calibration
recall overlay
Known Th2 regulators (red) on the plane. Activators land in/near the selective box; repressor-complex members sit high (Th1 induced), correctly excluded.

Independent review

6-agent literature + methods audit

Each dimension was audited by a separate agent with web-sourced citations. Verdicts and the top concern/fix per dimension:

Foundational immunologywell supported
Concern: 'Master regulators = GATA3 + STAT6' slightly conflates two roles: GATA3 is the lineage-defining transcription factor, whereas STAT6 is the upstream signal transducer of the IL-4 cue. GATA3 can drive Th2 commitment STAT6-independently via autoactivation (Stat6-deficient cells reconstituted by GATA3), so they are not interchangeable co-master TFs.
Fix: Frame STAT6 as the upstream IL-4 signal transducer feeding GATA3, not as a co-equal 'master regulator' transcription factor with GATA3.
Known-regulator direction reframemostly supported
Concern: The clean 'activator vs repressor' dichotomy oversimplifies the NuRD biology. CHD4 is explicitly BIMODAL (Hosokawa et al., PNAS 2013): the Gata3/Chd4/p300 complex ACTIVATES Th2 cytokines (IL-4/IL-5/IL-13) while the Gata3/Chd4-NuRD complex REPRESSES Tbx21. Chd4 knockdown decreases IL-4+ cells (18.3%->7.9%) AND increases IFN-g+ cells (3.7%->13.8%). So CHD4 KD suppresses Th2 too - it is not a pure Th1-repressor, and could legitimately score on the Th2-down arm.
Fix: TRAF3 is MISCLASSIFIED. The reframe lists TRAF3 as a negative regulator whose knockdown should induce Th1, but the direct T-cell evidence (Lu/Bishop lab, J Biol Chem 2024) shows TRAF3 ENHANCES STAT6 activation and PROMOTES Th2 skewing, and TRAF3 deficiency IMPAIRS Th2 (and Th1) differentiation. TRAF3 is therefore a POSITIVE Th2 regulator - its knockdown should LOWER Th2, behaving like the activators (GATA3/STAT6), not the repressors. The 'negative regulator' label likely derives from TRAF3's role in B-cell/alternative-NF-kB signaling, which does not transfer to the T-cell Th2 program.
Novelty & prior artproblematic
Concern: Novelty is low. The specific combination we pitched as our angle - split the Ota 2021 Th2-vs-Th1 signature into arms, score every genome-scale CRISPRi perturbation on both arms, decompose in 2D, and land on GATA3 - is a near-duplicate of the published Fig 4 / pert2state analysis in the very paper we are analyzing, using the identical discovery signature. Anyone reading their supplement + code would see this.
Fix: Drop or heavily soften the novelty claim. We should not present 'split the Th2-vs-Th1 signature into arms and score both' or the 2D decomposition as our idea - Zhu & Dann et al. already do exactly this with the same Ota 2021 signature (score_th2/score_th1, Supplementary Fig 17). Position our work as a re-analysis/reframing, at most contributing an explicit selective-suppression flatness criterion, not a new method.
Lead: ARNTmixed
Concern: Context mismatch: the HIF2alpha-Th2 link is established in HYPOXIC / tissue / intestinal-helminth settings and is hypoxia- or genetic-activation-dependent. Whether the ARNT/HIF2alpha axis meaningfully drives GATA3/IL4 in a 48hr in-vitro Stim of normoxic peripheral blood CD4+ T cells (typical Perturb-seq culture at ~21% O2, no HIF stabilization) is unproven and is the main reason to down-weight the lead.
Fix: Frame the mechanistic hypothesis as ARNT-via-HIF2alpha (EPAS1), not ARNT-via-AhR or generic HIF1alpha - HIF2alpha is the Th2-relevant partner and gives the cleanest, direction-correct rationale. Any writeup should name EPAS1/HIF2alpha as the presumed effector.
Lead: ELAVL1mostly supported
Concern: Global-artifact risk is real and non-trivial: HuR is a pan-cellular mRNA-stability hub with hundreds of ARE targets. Full/strong HuR loss impairs core T-cell activation, proliferation, and survival, and also stabilizes Foxp3 (Treg), IL-2, and Th17 program mRNAs. Strong CRISPRi knockdown could push cells toward the 'sick-cell' global-dampening regime our own analysis flagged as the dominant confound (Th2 and Th1 arms +0.5 correlated).
Fix: Do not describe HuR/ELAVL1 as a Th2-specific factor; it is a global ARE-binding mRNA stabilizer that happens to have a disproportionately characterized effect on the GATA3/IL4/IL13 regulon. The lead is plausible but the mechanism is 'global regulator with a strong Th2 sub-module,' not 'selective Th2 gene.'
Methods red-teammixed
Concern: SCORING IS SELF-CONTAINED, NOT COMPETITIVE: the arm score is an unweighted mean of per-gene DE z-scores over signature genes. This measures absolute movement of the signature, so any global transcriptional shift (sick cell) moves both arms together. The +0.5 correlation is therefore largely manufactured by the scoring choice, not a fact of nature to be regressed away afterward. Competitive/rank-based scores (AUCell, singscore, fgsea/CAMERA-style) normalize against the rest of the genome and would remove most of this confound at the scoring step.
Fix: Move the magnitude correction from a post-hoc regressor to the scoring step: replace the mean-of-z self-contained score with a competitive/rank-based enrichment (AUCell, singscore, or an fgsea/CAMERA/ROAST-style rotation test that normalizes against matched random gene sets). This removes most of the +0.5 confound a priori instead of trusting a noisy Th1-proxy regression.

Data/QC fixes already applied: neighboring-gene-KD filter, corrected donor metric (hits-mean, median 0.81), single-guide penalty.

Hardened re-analysis

answering the review

Directly addressing the methods critique: competitive rank-based scoring (removes the confound at the scoring step), a permutation null → empirical FDR, multi-lineage + proliferation control arms, and two-signature concordance.

+0.54
arm corr (mean-of-z)
the manufactured confound
+0.18
arm corr (competitive)
after rank-based scoring
420
selective, 2-arm
FDR 1.25
41
selective, multi-arm
+Th17/Treg/prolif flat · FDR 5.51
7
hardened survivors
+2-signature concordance
The rigorous re-analysis undercuts the original claim — honestly.
  • Competitive rank scoring collapses the arm correlation +0.54→+0.18 — most of the "global-magnitude confound" was a scoring artifact, as the review predicted.
  • Permutation null: the 2-arm selective set (420) is not enriched over random gene sets of the same size (expected ~525, empirical FDR ≈ 1.25). Genuine selective suppressors are no more common than chance here — a consequence of real Th1/Th2 reciprocity.
  • GATA3 fails under proper scoring: with the matched Th1-vs-Th0 arm, GATA3 KD drives the Th1 arm to +2.8 — it reads as a Th1-skewer, exactly as canonical biology predicts. The earlier "top selective hit" was a mean-of-z + noisy-Ota-arm artifact.
  • Only 7 genes pass the strict multi-arm + 2-signature concordance gate, and none of the earlier 21 survive; with FDR ≥ 1 even these are not credible targets.

Bottom line: the naive selective-Th2-suppressor atlas does not survive rigorous statistics. The real deliverables are the calibrated negative result, the reproducible pipeline, and the agentic review/QC infrastructure — not a target list. The right next step is a functional (IL-4/5/13) readout and control-guide-based FDR before any hit is claimed.

Review independently

self-contained + local files

This page embeds everything (figures, full shortlist, full dossiers, review). For click-through access to the raw files on your machine, open review_index.html at the project root — it links every artifact with relative paths.

  • dashboard.html
  • review_index.html
  • outputs/synthesis.md
  • outputs/review_report.md
  • outputs/real_vs_fp_shortlist.csv
  • outputs/selective_tiered.csv
  • outputs/condition_selectivity_matrix.csv
  • outputs/qc_audit_highconf.csv
  • outputs/plane_plot.svg
  • outputs/recall_overlay.svg
  • outputs/mechanism_dossiers/

Pipeline

2026-07-08
Done
Phase 0 · Setup & data access
Streamed the authoritative DE_stats.h5ad from public S3 via parallel HDF5 range reads (contiguous float64 matrix; only the 2.8 GB zscore layer, not the full 16.8 GB).
Done
Phase 1 · Reconstruct Th2 / Th1 arms
Projected every perturbation×condition onto the two Ota-2021 arm gene sets; z-scored; validated against the teammate’s checkpoint.
arm corr = +0.54 (teammate +0.49) · selective = 127 · skewers = 10 · GATA3 = top hit ✓
Done
Phase 2 · Credibility pack
Cutoff tiering · donor stability · KD self-silencing · known-regulator recall · per-condition split.
Known Th2 activators w/ detectable KD recovered 2/4 (GATA3, RARA; STAT6 trending) · Th1 regulators correctly excluded (0 in box)
Done
Phase 3 · Real-vs-false-positive score
KD-confirmed + real-footprint + magnitude-controlled selectivity (Th2 residual after regressing out the Th1 arm) → ranked shortlist.
GATA3 ranks #1 (0.90), RARA #2 — known positives at the top · high-confidence set = 21 genes
Done
Phase 4 · Agentic mechanism layer
Per top novel hit: agent investigates literature + pathways → mechanistic hypothesis bypassing T-bet → proposed confirming experiment, citations verified.
10 mechanism dossiers generated
Dataset: Zhu & Dann et al. 2025, genome-scale CRISPRi Perturb-seq in primary human CD4+ T cells (Marson lab). Signature: Ota et al. 2021. · Streamed live from the public Virtual Cells S3 bucket.