pith. sign in

arxiv: 2606.00568 · v2 · pith:WHR6S5PVnew · submitted 2026-05-30 · 💻 cs.LG · q-bio.GN

On the Recoverability of Causal Relations from Bulk Gene Expression Data

Pith reviewed 2026-06-28 19:15 UTC · model grok-4.3

classification 💻 cs.LG q-bio.GN
keywords causal discoverygene expressionbulk dataaggregationrecoverabilitylinearitystructural equations
0
0 comments X

The pith

Causal relations are recoverable from bulk gene expression data only under linear aggregation and affine structural equations.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper defines recoverability of causal gene relations from aggregated bulk data via two consistency properties. It derives necessary and sufficient conditions showing that both properties survive aggregation only when the aggregation step is linear, such as summing or averaging expression values, and the underlying equations linking genes are affine. Checks across four bulk and four single-cell datasets indicate that pairwise regulatory functions commonly depart from linearity. The results therefore advise against treating causal recovery from bulk expression data as routine without extra assumptions that enforce the required linearity.

Core claim

Recoverability under aggregation, defined by preservation of functional-form consistency and conditional-independence consistency, holds if and only if the aggregation is linear and the structural equations are affine.

What carries the argument

Functional-form consistency and conditional-independence consistency under aggregation.

If this is right

  • Nonlinear aggregation functions destroy both consistency properties.
  • Non-affine structural equations break recoverability even under linear aggregation.
  • Only sum or mean aggregation paired with affine equations preserves the original causal relations.
  • Empirical gene-expression data supplies limited support for the linearity conditions.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Causal methods applied to bulk data may need explicit linearity diagnostics or corrections before use.
  • Single-cell assays could avoid the aggregation barrier if their own measurement models satisfy the same conditions.
  • Alternative aggregation models that restore consistency without forcing linearity could be tested on the same datasets.

Load-bearing premise

Estimated pairwise regulatory functions from the eight datasets reliably indicate deviations from linearity.

What would settle it

A new dataset in which every estimated pairwise regulatory function is exactly linear would confirm that the linearity conditions can hold in practice.

Figures

Figures reproduced from arXiv: 2606.00568 by Boyang Sun, Gongxu Luo, Kun Zhang.

Figure 1
Figure 1. Figure 1: Illustration of aggregated bulk gene expression data. (a) The causal relationships among [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Illustration of the aggregated nature of bulk expression after standard normalization. [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: The three fundamental causal structure where we can only ob￾serve white nodes. of X and Z. Under the faithfulness assumption and the Markov condition, conditioning on Y (or any of its descen￾dants) opens the path between X and Z, rendering them conditionally dependent. While Y¯ is the descendant of Y , the observed variables satisfy X¯ ̸⊥⊥ Z¯|Y¯ , which is consis￾tent with the true underlying latent relati… view at source ↗
Figure 4
Figure 4. Figure 4: Distribution of inverse smooth edf and signal-weighted linearity score across microarray (a) & (b) and bulk RNA-seq (c) & (d) datasets. GSE39582 (585 samples) [34] and GSE20142 (1240 samples) [35], and two bulk RNA-seq datasets, GSE57945 (260 samples) [36] and GSE245006 (518 samples) [37]. The use of bulk data is motivated by the principle above: if the underlying regulatory relations are predominantly aff… view at source ↗
Figure 5
Figure 5. Figure 5: Distribution of inverse smooth edf and signal-weighted linearity score across single-cell datasets. where s(·) is a smooth function estimated from the data. The effective degrees of freedom, denoted by edf c j→l , measure the flexibility used by the fitted smooth function. An edf close to 1 indicates that the fitted function is approximately linear, whereas a larger edf indicates stronger nonlinearity. We … view at source ↗
read the original abstract

Bulk gene expression profiling, which aggregates pooled RNA across cells within a biological sample, remains important in the single-cell era because it is typically less noisy, more sensitive, and more cost-effective than single-cell assays. Accordingly, a growing body of computational methods seeks to recover causal relations among genes from bulk expression data. However, aggregation is a lossy, non-invertible coarsening of the underlying cellular system, and it remains unclear whether and under what conditions causal relations are recoverable from aggregated bulk gene expression data. To answer this, we formalize recoverability under aggregation through two notions of consistency: functional-form consistency and conditional-independence consistency. We then derive necessary and sufficient conditions for recoverability, showing that these properties are preserved only under linear aggregations (e.g., sum/mean) coupled with affine structural equations. To assess the practical plausibility of these conditions, analyses of four bulk and four single-cell gene expression datasets further reveal that the estimated pairwise regulatory functions among genes deviate from linearity in both data types, providing limited empirical support for the linearity assumptions required for recoverability. Together, these results caution against recovering causal relations from aggregated bulk expression data without strong additional assumptions.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper formalizes recoverability of causal relations from bulk (aggregated) gene expression data via two consistency notions—functional-form consistency and conditional-independence consistency. It derives necessary and sufficient conditions under which these consistencies are preserved, concluding that recoverability holds only for linear aggregations (e.g., sum or mean) paired with affine structural equations. Analyses of four bulk and four single-cell datasets are used to estimate pairwise regulatory functions, which are reported to deviate from linearity in both data types, leading to the conclusion of limited empirical support for the required assumptions and a caution against causal recovery from bulk data without strong additional assumptions.

Significance. If the derivation is correct, the work supplies a precise theoretical boundary on when causal structure can survive aggregation, which is valuable given the continued use of bulk profiling. The explicit linkage of recoverability to linear aggregation plus affine equations, together with the multi-dataset empirical check, offers a falsifiable framing that could guide future method development. The absence of free parameters in the consistency definitions is a strength.

major comments (2)
  1. [§3 (derivation of nec/suff conditions)] Theoretical derivation (likely §3): the manuscript states that necessary and sufficient conditions for recoverability are derived, yet supplies no lemmas, proof steps, or verification that non-linear aggregations necessarily violate functional-form consistency while linear ones preserve both consistencies. Without these details the central claim cannot be assessed.
  2. [§4–5 (dataset analyses and linearity checks)] Empirical section (likely §4–5): the claim that estimated pairwise regulatory functions deviate from linearity (and thus violate the affine requirement) rests on unspecified regression procedures, absence of reported statistical tests for linearity, and no controls for unmeasured confounding or aggregation-induced bias across the eight datasets. This directly weakens the auxiliary claim that real data provide limited support for recoverability.
minor comments (2)
  1. [§2] Notation for the two consistency notions is introduced without an explicit comparison table; a side-by-side definition would improve readability.
  2. [Figures 3–6] Figure captions for the regulatory-function plots do not state the exact functional form fitted or the number of gene pairs examined per dataset.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their thoughtful review and constructive comments. We address each major comment below.

read point-by-point responses
  1. Referee: [§3 (derivation of nec/suff conditions)] Theoretical derivation (likely §3): the manuscript states that necessary and sufficient conditions for recoverability are derived, yet supplies no lemmas, proof steps, or verification that non-linear aggregations necessarily violate functional-form consistency while linear ones preserve both consistencies. Without these details the central claim cannot be assessed.

    Authors: We agree that the detailed proof steps and lemmas were not included in the submitted manuscript, which limits the assessability of the central theoretical claim. In the revised manuscript, we will provide the full derivation in an appendix or dedicated section, including lemmas establishing the necessary and sufficient conditions for both consistency notions under linear vs. non-linear aggregations paired with affine structural equations. revision: yes

  2. Referee: [§4–5 (dataset analyses and linearity checks)] Empirical section (likely §4–5): the claim that estimated pairwise regulatory functions deviate from linearity (and thus violate the affine requirement) rests on unspecified regression procedures, absence of reported statistical tests for linearity, and no controls for unmeasured confounding or aggregation-induced bias across the eight datasets. This directly weakens the auxiliary claim that real data provide limited support for recoverability.

    Authors: We acknowledge the need for greater transparency in the empirical section. We will revise to explicitly describe the regression procedures employed for estimating pairwise regulatory functions, report statistical tests for deviations from linearity, and address potential issues of unmeasured confounding and aggregation bias, including any limitations. These additions will strengthen the presentation while preserving the observed deviations from linearity in the analyzed datasets. revision: yes

Circularity Check

0 steps flagged

Derivation of recoverability conditions proceeds directly from explicit consistency definitions without reduction to inputs or self-citations.

full rationale

The paper first defines functional-form consistency and conditional-independence consistency, then states necessary and sufficient conditions for recoverability under aggregation. These conditions are presented as logical consequences of the definitions (linear aggregation + affine equations preserve the consistencies). No step reduces a derived quantity to a fitted parameter, renames an input, or relies on a load-bearing self-citation whose content is unverified. The empirical dataset analyses are auxiliary and concern plausibility of the affine assumption; they do not enter the theoretical derivation. The derivation chain is therefore self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on domain assumptions about the form of structural causal models for gene expression and the nature of aggregation as a coarsening operation.

axioms (2)
  • domain assumption Gene expression data arises from a structural causal model
    Required to define causal relations and the consistency properties
  • domain assumption Bulk data results from a non-invertible aggregation of single-cell expressions
    Stated as the core setup for the recoverability question

pith-pipeline@v0.9.1-grok · 5740 in / 1221 out tokens · 39963 ms · 2026-06-28T19:15:13.444895+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

42 extracted references

  1. [1]

    Properties of developmental gene regulatory networks

    Eric H Davidson and Michael S Levine. Properties of developmental gene regulatory networks. Proceedings of the National Academy of Sciences, 105(51):20063–20066, 2008

  2. [2]

    recount3: summaries and queries for large-scale rna-seq expression and splicing.Genome biology, 22(1):323, 2021

    Christopher Wilks, Shijie C Zheng, Feng Yong Chen, Rone Charles, Brad Solomon, Jonathan P Ling, Eddie Luidy Imada, David Zhang, Lance Joseph, Jeffrey T Leek, et al. recount3: summaries and queries for large-scale rna-seq expression and splicing.Genome biology, 22(1):323, 2021

  3. [3]

    Perturbatlas: a comprehensive atlas of public genetic perturbation bulk rna-seq datasets.Nucleic Acids Research, 53(D1):D1112–D1119, 2025

    Yiming Zhang, Ting Zhang, Gaoxia Yang, Zhenzhong Pan, Min Tang, Yue Wen, Ping He, Yuan Wang, and Ran Zhou. Perturbatlas: a comprehensive atlas of public genetic perturbation bulk rna-seq datasets.Nucleic Acids Research, 53(D1):D1112–D1119, 2025

  4. [4]

    Aracne: an algorithm for the reconstruction of gene regulatory networks in a mammalian cellular context.BMC bioinformatics, 7(Suppl 1):S7, 2006

    Adam A Margolin, Ilya Nemenman, Katia Basso, Chris Wiggins, Gustavo Stolovitzky, Ric- cardo Dalla Favera, and Andrea Califano. Aracne: an algorithm for the reconstruction of gene regulatory networks in a mammalian cellular context.BMC bioinformatics, 7(Suppl 1):S7, 2006

  5. [5]

    Inferring gene regulatory networks from gene expression data by path consistency algorithm based on conditional mutual information.Bioinformatics, 28(1):98–104, 2012

    Xiujun Zhang, Xing-Ming Zhao, Kun He, Le Lu, Yongwei Cao, Jingdong Liu, Jin-Kao Hao, Zhi-Ping Liu, and Luonan Chen. Inferring gene regulatory networks from gene expression data by path consistency algorithm based on conditional mutual information.Bioinformatics, 28(1):98–104, 2012

  6. [6]

    A new dynamic bayesian network (dbn) approach for identi- fying gene regulatory networks from time course microarray data.Bioinformatics, 21(1):71–79, 2005

    Min Zou and Suzanne D Conzen. A new dynamic bayesian network (dbn) approach for identi- fying gene regulatory networks from time course microarray data.Bioinformatics, 21(1):71–79, 2005

  7. [7]

    Inference of gene regulatory network based on local bayesian networks.PLoS computational biology, 12(8):e1005024, 2016

    Fei Liu, Shao-Wu Zhang, Wei-Feng Guo, Ze-Gang Wei, and Luonan Chen. Inference of gene regulatory network based on local bayesian networks.PLoS computational biology, 12(8):e1005024, 2016

  8. [8]

    Spred: a simulation-supervised neural network tool for gene regulatory network reconstruction.Bioinformatics Advances, 4(1):vbae011, 2024

    Zijun Wu and Saurabh Sinha. Spred: a simulation-supervised neural network tool for gene regulatory network reconstruction.Bioinformatics Advances, 4(1):vbae011, 2024

  9. [9]

    A hybrid constrained continuous optimization approach for optimal causal discovery from biological data.Bioinformatics, 40(Supplement_2):ii87–ii97, 2024

    Yuehua Zhu, Panayiotis V Benos, and Maria Chikina. A hybrid constrained continuous optimization approach for optimal causal discovery from biological data.Bioinformatics, 40(Supplement_2):ii87–ii97, 2024

  10. [10]

    Computational methods for gene regulatory networks reconstruction and analysis: a review.Artificial intelligence in medicine, 95:133–145, 2019

    Fernando M Delgado and Francisco Gómez-Vela. Computational methods for gene regulatory networks reconstruction and analysis: a review.Artificial intelligence in medicine, 95:133–145, 2019

  11. [11]

    Review of causal discovery methods based on graphical models.Frontiers in genetics, 10:524, 2019

    Clark Glymour, Kun Zhang, and Peter Spirtes. Review of causal discovery methods based on graphical models.Frontiers in genetics, 10:524, 2019

  12. [12]

    A fast pc algo- rithm for high dimensional causal discovery with multi-core pcs.IEEE/ACM transactions on computational biology and bioinformatics, 16(5):1483–1495, 2016

    Thuc Duy Le, Tao Hoang, Jiuyong Li, Lin Liu, Huawen Liu, and Shu Hu. A fast pc algo- rithm for high dimensional causal discovery with multi-core pcs.IEEE/ACM transactions on computational biology and bioinformatics, 16(5):1483–1495, 2016

  13. [13]

    Construction of pan-cancer regulatory networks based on causal inference.Biosystems, 243:105279, 2024

    Ruirui Ji, Mengfei Yan, Meng Zhao, and Yi Geng. Construction of pan-cancer regulatory networks based on causal inference.Biosystems, 243:105279, 2024

  14. [14]

    Gene network inference using con- tinuous time bayesian networks: a comparative study and application to th17 cell differentiation

    Enzo Acerbi, Teresa Zelante, Vipin Narang, and Fabio Stella. Gene network inference using con- tinuous time bayesian networks: a comparative study and application to th17 cell differentiation. BMC bioinformatics, 15(1):387, 2014

  15. [15]

    Hamda B Ajmal and Michael G Madden. Inferring dynamic gene regulatory networks with low-order conditional independencies–an evaluation of the method.Statistical Applications in Genetics and Molecular Biology, 19(4-6):20200051, 2020

  16. [16]

    Co-expressed gene assessment based on the path consistency algorithm: operon detention in escherichia coli

    Shigeru Saito and Katsuhisa Horimoto. Co-expressed gene assessment based on the path consistency algorithm: operon detention in escherichia coli. In2009 IEEE International Conference on Systems, Man and Cybernetics, pages 4280–4286. IEEE, 2009. 10

  17. [17]

    Mrpc: An r package for inference of causal graphs.Frontiers in Genetics, 12:651812, 2021

    Md Bahadur Badsha, Evan A Martin, and Audrey Qiuyan Fu. Mrpc: An r package for inference of causal graphs.Frontiers in Genetics, 12:651812, 2021

  18. [18]

    Using bayesian networks to analyze expression data

    Nir Friedman, Michal Linial, Iftach Nachman, and Dana Pe’er. Using bayesian networks to analyze expression data. InProceedings of the fourth annual international conference on Computational molecular biology, pages 127–135, 2000

  19. [19]

    An improved bayesian network method for reconstructing gene regulatory network based on candidate auto selection.BMC genomics, 18(Suppl 9):844, 2017

    Linlin Xing, Maozu Guo, Xiaoyan Liu, Chunyu Wang, Lei Wang, and Yin Zhang. An improved bayesian network method for reconstructing gene regulatory network based on candidate auto selection.BMC genomics, 18(Suppl 9):844, 2017

  20. [20]

    High-dimensional bayesian network inference from systems genetics data using genetic node ordering.Frontiers in genetics, 10:1196, 2019

    Lingfei Wang, Pieter Audenaert, and Tom Michoel. High-dimensional bayesian network inference from systems genetics data using genetic node ordering.Frontiers in genetics, 10:1196, 2019

  21. [21]

    Gene network inference via structural equation modeling in genetical genomics experiments.Genetics, 178(3):1763–1776, 2008

    Bing Liu, Alberto de La Fuente, and Ina Hoeschele. Gene network inference via structural equation modeling in genetical genomics experiments.Genetics, 178(3):1763–1776, 2008

  22. [22]

    Inference of gene regula- tory networks with sparse structural equation models exploiting genetic perturbations.PLoS computational biology, 9(5):e1003068, 2013

    Xiaodong Cai, Juan Andrés Bazerque, and Georgios B Giannakis. Inference of gene regula- tory networks with sparse structural equation models exploiting genetic perturbations.PLoS computational biology, 9(5):e1003068, 2013

  23. [23]

    Bivariate causal discovery and its applications to gene expression and imaging data analysis.Frontiers in genetics, 9:347, 2018

    Rong Jiao, Nan Lin, Zixin Hu, David A Bennett, Li Jin, and Momiao Xiong. Bivariate causal discovery and its applications to gene expression and imaging data analysis.Frontiers in genetics, 9:347, 2018

  24. [24]

    Differential analysis of gene regulatory networks modeled with structural equation models.Journal of Ambient Intelligence and Humanized Computing, 12(10):9181–9192, 2021

    Yan Li, Dayou Liu, Yungang Zhu, and Jie Liu. Differential analysis of gene regulatory networks modeled with structural equation models.Journal of Ambient Intelligence and Humanized Computing, 12(10):9181–9192, 2021

  25. [25]

    Expression profiling using cdna microarrays.Nature genetics, 21(1):10–14, 1999

    David J Duggan, Michael Bittner, Yidong Chen, Paul Meltzer, and Jeffrey M Trent. Expression profiling using cdna microarrays.Nature genetics, 21(1):10–14, 1999

  26. [26]

    Affymetrix genechip microarray preprocessing for multivariate analyses.Briefings in bioinformatics, 13(5):536–546, 2012

    Matthew N McCall and Anthony Almudevar. Affymetrix genechip microarray preprocessing for multivariate analyses.Briefings in bioinformatics, 13(5):536–546, 2012

  27. [27]

    Microarray data normalization and transformation.Nature genetics, 32(4):496–501, 2002

    John Quackenbush. Microarray data normalization and transformation.Nature genetics, 32(4):496–501, 2002

  28. [28]

    Normalization of microarray data: single-labeled and dual-labeled arrays.Molecules and cells, 22(3):254–261, 2006

    Jin Hwan Do and Dong-Kug Choi. Normalization of microarray data: single-labeled and dual-labeled arrays.Molecules and cells, 22(3):254–261, 2006

  29. [29]

    Mapping and quantifying mammalian transcriptomes by rna-seq.Nature methods, 5(7):621– 628, 2008

    Ali Mortazavi, Brian A Williams, Kenneth McCue, Lorian Schaeffer, and Barbara Wold. Mapping and quantifying mammalian transcriptomes by rna-seq.Nature methods, 5(7):621– 628, 2008

  30. [30]

    Rna-seq: a revolutionary tool for transcrip- tomics.Nature reviews genetics, 10(1):57–63, 2009

    Zhong Wang, Mark Gerstein, and Michael Snyder. Rna-seq: a revolutionary tool for transcrip- tomics.Nature reviews genetics, 10(1):57–63, 2009

  31. [31]

    A scaling normalization method for differential expres- sion analysis of rna-seq data.Genome biology, 11(3):R25, 2010

    Mark D Robinson and Alicia Oshlack. A scaling normalization method for differential expres- sion analysis of rna-seq data.Genome biology, 11(3):R25, 2010

  32. [32]

    Revisiting global gene expression analysis.Cell, 151(3):476–482, 2012

    Jakob Lovén, David A Orlando, Alla A Sigova, Charles Y Lin, Peter B Rahl, Christopher B Burge, David L Levens, Tong Ihn Lee, and Richard A Young. Revisiting global gene expression analysis.Cell, 151(3):476–482, 2012

  33. [33]

    Moderated estimation of fold change and dispersion for rna-seq data with deseq2.Genome biology, 15(12):550, 2014

    Michael I Love, Wolfgang Huber, and Simon Anders. Moderated estimation of fold change and dispersion for rna-seq data with deseq2.Genome biology, 15(12):550, 2014

  34. [34]

    Gene expression classification of colon cancer into molecular subtypes: characterization, validation, and prognostic value.PLoS medicine, 10(5):e1001453, 2013

    Laetitia Marisa, Aurélien de Reyniès, Alex Duval, Janick Selves, Marie Pierre Gaub, Laure Vescovo, Marie-Christine Etienne-Grimaldi, Renaud Schiappa, Dominique Guenot, Mira Ayadi, et al. Gene expression classification of colon cancer into molecular subtypes: characterization, validation, and prognostic value.PLoS medicine, 10(5):e1001453, 2013. 11

  35. [35]

    Multiple common variants for celiac disease influencing immune gene expression.Nature genetics, 42(4):295–302, 2010

    Patrick CA Dubois, Gosia Trynka, Lude Franke, Karen A Hunt, Jihane Romanos, Alessandra Curtotti, Alexandra Zhernakova, Graham AR Heap, Róza Ádány, Arpo Aromaa, et al. Multiple common variants for celiac disease influencing immune gene expression.Nature genetics, 42(4):295–302, 2010

  36. [36]

    Unraveling the role of the mitochondrial one-carbon pathway in undifferentiated thyroid cancer by multi-omics analyses

    Seong Eun Lee, Seongyeol Park, Shinae Yi, Na Rae Choi, Mi Ae Lim, Jae Won Chang, Ho- Ryun Won, Je Ryong Kim, Hye Mi Ko, Eun-Jae Chung, et al. Unraveling the role of the mitochondrial one-carbon pathway in undifferentiated thyroid cancer by multi-omics analyses. Nature Communications, 15(1):1163, 2024

  37. [37]

    Subtype- selective prenylated isoflavonoids disrupt regulatory drivers of mycn-amplified cancers.Cell chemical biology, 31(4):805–819, 2024

    Michael E Stokes, Alessandro Vasciaveo, Jonnell Candice Small, Arie Zask, Eduard Reznik, Nailah Smith, Qian Wang, Jacob Daniels, Farhad Forouhar, Presha Rajbhandari, et al. Subtype- selective prenylated isoflavonoids disrupt regulatory drivers of mycn-amplified cancers.Cell chemical biology, 31(4):805–819, 2024

  38. [38]

    A multiplexed single-cell crispr screening platform enables systematic dissection of the unfolded protein response.Cell, 167(7):1867–1882, 2016

    Britt Adamson, Thomas M Norman, Marco Jost, Min Y Cho, James K Nuñez, Yuwen Chen, Jacqueline E Villalta, Luke A Gilbert, Max A Horlbeck, Marco Y Hein, et al. A multiplexed single-cell crispr screening platform enables systematic dissection of the unfolded protein response.Cell, 167(7):1867–1882, 2016

  39. [39]

    Exploring genetic interaction manifolds constructed from rich single-cell phenotypes.Science, 365(6455):786–793, 2019

    Thomas M Norman, Max A Horlbeck, Joseph M Replogle, Alex Y Ge, Albert Xu, Marco Jost, Luke A Gilbert, and Jonathan S Weissman. Exploring genetic interaction manifolds constructed from rich single-cell phenotypes.Science, 365(6455):786–793, 2019

  40. [40]

    Massively multiplex chemical transcriptomics at single-cell resolution.Science, 367(6473):45– 51, 2020

    Sanjay R Srivatsan, José L McFaline-Figueroa, Vijay Ramani, Lauren Saunders, Junyue Cao, Jonathan Packer, Hannah A Pliner, Dana L Jackson, Riza M Daza, Lena Christiansen, et al. Massively multiplex chemical transcriptomics at single-cell resolution.Science, 367(6473):45– 51, 2020

  41. [41]

    Genome-wide crispri/a screens in human neurons link lysosomal failure to ferroptosis.Nature neuroscience, 24(7):1020–1034, 2021

    Ruilin Tian, Anthony Abarientos, Jason Hong, Sayed Hadi Hashemi, Rui Yan, Nina Dräger, Kun Leng, Mike A Nalls, Andrew B Singleton, Ke Xu, et al. Genome-wide crispri/a screens in human neurons link lysosomal failure to ferroptosis.Nature neuroscience, 24(7):1020–1034, 2021

  42. [42]

    Curation of over 10 000 transcriptomic studies to enable data reuse.Database, 2021:baab006, 2021

    Nathaniel Lim, Stepan Tesar, Manuel Belmadani, Guillaume Poirier-Morency, Burak Ogan Mancarci, Jordan Sicherman, Matthew Jacobson, Justin Leong, Patrick Tan, and Paul Pavlidis. Curation of over 10 000 transcriptomic studies to enable data reuse.Database, 2021:baab006, 2021. 12 A Proof for Functional-form Consistency A.1 Proof for Theorem 4.3 Theorem 4.3(R...