Pseudotimecascade takes an scRNA-seq matrix plus a cell ordering (pseudotime or RNA velocity), fits smoothed per-gene curves with VGAM GAMs, applies BH correction, scales fits, and classifies genes by monotonic/non-monotonic temporal patterns. Genes are grouped by pattern and ordered by their first switch point (zero-crossing), then visualized as a cascade heatmap; a sliding-window GO analysis adds stage-specific pathway annotation . The published CSBJ version adds left-censoring via Tobit, alternative switch-point definitions (threshold, max first derivative), and ssGSEA pathway scoring, with two datasets: an HCA bone marrow sample and GSE253355 (12 donors) .
The schematic above reproduces the reported temporal ordering β early hematopoietic genes (CEBPA, SPI1, CD34) declining; GATA1/KLF1 rising at intermediate stages; HBA1/HBA2 late; GATA2 repressed β as described by the authors; exact fitted values were not published in the text, so this is qualitative, not quantitative reproduction .
Strengths: The core insight β treating genes as a coordinated cascade ordered by switch points rather than independently β is practical and well-motivated; the paper shows hierarchical/k-means orderings fail to resolve temporal structure, supporting the added value. Temporal GO (early: cell cycle/remodeling; late: hemoglobin biosynthesis, heme metabolism) aligns with known erythropoiesis biology .
Limitations the paper does not fully resolve: (1) Results inherit all errors of the input ordering β no sensitivity analysis across multiple pseudotime methods (e.g., Monocle vs DPT vs TSCAN) is reported despite this being the dominant uncertainty. (2) GO enrichment is annotation-dependent, and the paper acknowledges no external experimental validation of the inferred cascade hierarchy. (3) Scope is narrow: two bone marrow datasets; generalization to branching trajectories (where switch-point ordering across branches is ambiguous), other tissues, or species remains untested . (4) Statistical novelty is thin β VGAM smoothing plus zero-crossing detection and BH correction are standard components; the contribution is ordering + visualization logic, not new inference. Confidence in the erythroid findings is moderate (consistent with established GATA1/KLF1 biology); confidence in broad generality is low.
Know what changed, what holds up, and what remains uncertain. Every Friday. No ads.