On the Recoverability of Causal Relations from Bulk Gene Expression Data
Pith reviewed 2026-06-28 19:15 UTC · model grok-4.3
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.
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
- 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
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.
Referee Report
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)
- [§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.
- [§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)
- [§2] Notation for the two consistency notions is introduced without an explicit comparison table; a side-by-side definition would improve readability.
- [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
We thank the referee for their thoughtful review and constructive comments. We address each major comment below.
read point-by-point responses
-
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
-
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
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
axioms (2)
- domain assumption Gene expression data arises from a structural causal model
- domain assumption Bulk data results from a non-invertible aggregation of single-cell expressions
Reference graph
Works this paper leans on
-
[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
2008
-
[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
2021
-
[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
2025
-
[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
2006
-
[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
2012
-
[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
2005
-
[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
2016
-
[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
2024
-
[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
2024
-
[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
2019
-
[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
2019
-
[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
2016
-
[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
2024
-
[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
2014
-
[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
2020
-
[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
2009
-
[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
2021
-
[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
2000
-
[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
2017
-
[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
2019
-
[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
2008
-
[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
2013
-
[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
2018
-
[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
2021
-
[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
1999
-
[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
2012
-
[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
2002
-
[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
2006
-
[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
2008
-
[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
2009
-
[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
2010
-
[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
2012
-
[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
2014
-
[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
2013
-
[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
2010
-
[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
2024
-
[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
2024
-
[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
2016
-
[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
2019
-
[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
2020
-
[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
2021
-
[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...
2021
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.