Campaign run · Alzheimer's disease
alzheimer-rorb-ec-vulnerability-r1
What cell-intrinsic transcriptional programs — specifically excitability regulation, calcium handling, proteostasis capacity, and mitochondrial gene expression — distinguish selectively vulnerable RORB+ entorhinal cortex (EC) layer II excitatory neurons from resilient neurons in the same region, and can these molecular features predict inter-individual variation in tau accumulation rate across Braak stages using existing OPEN snRNA-seq atlases?
Status
- How far it got
- Completed the pipeline
Tested at least one hypothesis and reached the multi-agent review stage. This shows the run went through the whole pipeline, not that its work passed a quality bar.
- Hypotheses
8 nominated; 6 evaluated; 2 not tested because the run did not execute them.
- Recorded status
Done. The run recorded a final status of done.
- Census class
Success. Census class is a label from an earlier, process-level census of the campaign. It does not map one-to-one onto the stages above. How the classes compare.
Research question
The sentence under the run id is quoted from this brief.
Full research brief given to the run
Direction d1 from the Alzheimer's-mechanism IDEATE run.
What cell-intrinsic transcriptional programs — specifically excitability regulation, calcium handling, proteostasis capacity, and mitochondrial gene expression — distinguish selectively vulnerable RORB+ entorhinal cortex (EC) layer II excitatory neurons from resilient neurons in the same region, and can these molecular features predict inter-individual variation in tau accumulation rate across Braak stages using existing OPEN snRNA-seq atlases?
EC layer II neurons are the first to accumulate phospho-tau and die, defining Braak stage I–II. Leng et al. (2021, Nat Neurosci) showed RORB expression marks these vulnerable neurons with ~50% abundance loss at early stages. WHY RORB+ neurons are vulnerable — high intrinsic excitability, low mitochondrial capacity, limited proteostasis reserve, or a combination — remains unresolved. This OmiCairn run re-analyses open EC snRNA-seq atlases to test candidate cell-intrinsic vulnerability mechanisms by contrasting RORB+ vs. RORB- excitatory neurons within region and across Braak stage.
Hypotheses (pre-registered, lifted from d1):
H1 (primary — calcium-handling / excitability signature): RORB+ EC layer II excitatory neurons express significantly higher levels of voltage-gated calcium-channel subunits (CACNA1C, CACNA1D) and lower levels of calcium-buffering genes (CALB1, PVALB) relative to resilient RORB- EC excitatory neurons, consistent with a calcium-/excitability-driven basis for their selective tau vulnerability. Decision rule: per-donor pseudobulk DE within EC excitatory neurons, RORB+ vs RORB-; report log2FC + Benjamini–Hochberg FDR + 95% bootstrap CI for each gene. H1 passes if CACNA1C and/or CACNA1D are up (log2FC > 0, FDR < 0.05) AND CALB1 and/or PVALB are down (log2FC < 0, FDR < 0.05) in the primary atlas, with the directional sign reproduced in the Braak-staged secondary atlas.
H2 (secondary — integrated stress response / SASP precursor trajectory): Across Braak stages I–IV in EC snRNA-seq data, RORB+ neurons show a progressive shift toward an integrated stress response (ATF4-high; phospho-eIF2α surrogate genes — e.g. ATF4, DDIT3/CHOP, PPP1R15A, ATF3, TRIB3 — upregulated) and a senescence-associated secretory (SASP) precursor state, detectable before frank neuronal loss. Decision rule: ISR/SASP module score per RORB+ cell regressed on Braak stage (ordinal); report slope + p-value + 95% CI. H2 passes if the ISR module score increases monotonically with Braak stage in RORB+ neurons (positive slope, p < 0.05) and the increase precedes the inflection point of RORB+ abundance loss. Report the RORB+ abundance-vs-Braak curve alongside.
A negative or partial result is fully informative. If the calcium / excitability contrast does not separate RORB+ from RORB- neurons, or if the ISR shift is absent, that finding must be reported as the headline rather than re-tuning clustering or gene panels until H1/H2 pass. Report which decomposition (excitability vs proteostasis vs mitochondrial) carries the signal.
The research question comes from the topic's current research brief. All attempts at a topic share it; earlier attempts may have started from an earlier wording.
Hypotheses and verdicts
Each hypothesis was written with its decision rule before any analysis (or, for on-the-fly hypotheses, recorded during execution). AI adjudication agents applied that rule to the run's latest outputs and quoted the evidence verbatim. Where the data did not meet the rule, the verdict follows the data, even if the run's own label or report claimed support.
| ID | Hypothesis and decision rule | Verdict |
|---|---|---|
| H1Pre-registered | RORB+ vs RORB- EC excitatory neurons: CACNA1C or CACNA1D up AND CALB1 or PVALB down (BH-FDR<0.05), sign reproduced in second atlas. | Not supportedIn the explorer |
Evidence (verbatim from the run outputs) H1 verdict: NOT SUPPORTED. No gene reaches BH-FDR < 0.05. Adjudication note Latest (008, GSE147528 primary, 2,009 RORB+ cells/10 donors): CACNA1C/D opposite sign, CIs include 0; PVALB log2FC -1.18 CI excludes 0 but FDR 0.865 (t-test NaN issue). Report calls it a powered negative. Atlas roles swapped vs plan (GSE138852 had only 16 RORB+ cells). | ||
| H2Pre-registered | ISR module score in RORB+ EC neurons increases with Braak stage (slope>0, p<0.05, CI excl. 0), preceding RORB+ loss. | Not supportedIn the explorer |
Evidence (verbatim from the run outputs) failing the pre-registered H2 criterion of a positive slope at p < 0.05 Adjudication note Slope -0.006/Braak (CI [-0.018, +0.002], p=0.213), opposite direction; n=10 donors, Braak 0/2/6 only (plan: borderline power). | ||
| H3Pre-registered | Exploratory: one vulnerability program (excitability, ISR/proteostasis, mitochondrial) separates RORB+ from RORB- more strongly than others (Cohen's d, >0.3 SD). | InconclusiveIn the explorer |
Evidence (verbatim from the run outputs) An exploratory excitability gene module score is positive in direction (d = +0.269) but spans zero at the donor level of inference (approximate CI [−0.362, +0.900]) Adjudication note Excitability had the largest |d| but below the 0.3 SD benchmark; cell-level bootstrap CIs anti-conservative, donor-level CI spans 0. Knowledge-graph node ('excitability module higher') = supported (discrepancy). | ||
| H4Pre-registered | Exploratory: RORB+ EC excitatory subcluster recoverable at Leiden 0.6-1.2 in both GSE138852 and GSE147528 with >80% consistent assignment. | Not supportedIn the explorer |
Evidence (verbatim from the run outputs) This is a Tier-A negative finding: the primary atlas does not contain enough RORB+ neurons for valid pseudobulk DE Adjudication note GSE138852: 16 RORB+ nuclei, 1/6 pools >=5 cells (pre-registered Tier-A negative). Leiden resolution sweep never run (OMP crash); threshold-based annotation used. Knowledge-graph status = refuted. | ||
| OTF1Proposed on the fly | ISR module is LOWER in RORB+ than RORB- EC neurons (at all Braak stages), counter to the ISR-driven vulnerability model. | InconclusiveIn the explorer |
Evidence (verbatim from the run outputs) The ISR module score was lower in RORB+ than RORB- neurons (d = −0.056, cell-level CI [−0.100, −0.008]; approximate donor-level CI [−0.678, +0.566]; Fig. 5a) Adjudication note Knowledge-graph status = supported (cell-level CI); donor-level CI spans 0. 'At all Braak stages' not tested stage-stratified. | ||
| OTF2Proposed on the fly | Re-analysis restricted to GSE147528 (2,009 RORB+ cells) will have power for H1; test whether CACNA1C/D reach FDR<0.05 there. | Not supportedIn the explorer |
Evidence (verbatim from the run outputs) CACNA1C and CACNA1D show small NEGATIVE log2FC (opposite the H1 prediction of upregulation in RORB+). Adjudication note Acted on in EXECUTE stage 8 (atlas switch); CACNA1C/D FDR=1.00. | ||
| OTF3Proposed on the fly | Elevated excitability module in RORB+ neurons is driven by non-calcium genes (SCN1A, SCN8A, HCN1). | Not runIn the explorer |
Evidence (verbatim from the run outputs) per-gene identification of driver genes within the module has not been performed Adjudication note Data available; not run. | ||
| OTF4Proposed on the fly | PVALB downregulation in RORB+ neurons (log2FC -1.18) reaches BH-FDR<0.05 under pydeseq2. | Not runIn the explorer |
Evidence (verbatim from the run outputs) For the pseudobulk DE, use pydeseq2 (once OMP issues are resolved) which handles sparse genes more robustly than the simple t-test. Adjudication note Technical: pydeseq2 not run because of native OMP/BLAS library crashes on the analysis platform. | ||