Causal Cellular Context Transfer Learning (C3TL): An Efficient Architecture for Prediction of Unseen Perturbation Effects
Authors: Michael Scholkemper, Sach Mukherjee
Organizations: Statistics and Machine Learning, German Center for Neurodegenerative Diseases (DZNE), Bonn, Germany. · Statistics and Machine Learning, German Center for Neurodegenerative Diseases (DZNE) Bonn, Germany, University of Bonn, Germany & University of Cambridge, United Kingdom.
Predicting the effects of chemical and genetic perturbations on quantitative cell states is a central challenge in computational biology, molecular medicine and drug discovery. Recent work has leveraged large-scale single-cell data and massive foundation models to address this task. However, such computational resources and extensive datasets are not always accessible in academic or clinical settings, hence limiting utility. Here we propose a lightweight framework for perturbation effect prediction that exploits the structured nature of biological interventions and specific inductive biases/invariances. Our approach leverages available information concerning perturbation effects to allow generalization to novel contexts and requires only widely-available bulk molecular data. Extensive testing, comparing predictions of context-specific perturbation effects against real, large-scale interventional experiments, demonstrates accurate prediction in new contexts. The proposed approach is competitive with SOTA foundation models but requires simpler data, much smaller model sizes and less time. Focusing on robust bulk signals and efficient architectures, we show that accurate prediction of perturbation effects is possible without proprietary hardware or very large models, hence opening up ways to leverage causal learning approaches in biomedicine generally.
Figures & tables
Figure 1: Proposed architecture. The model compresses the high-dimensional gene expression input {xpγ}γ from multiple contexts (indexed by γ ) but under the same perturbation p into a lower-dimensional latent representation z^p of the perturbation itself. Similarly, it forms a representation ψ^c of the context from all the different perturbations (indexed by π ) observed in the specific context c . From these latent representations, it then reconstructs the original gene expression x^pc using a decoder Γ .
Replogle
Parse
Tahoe
Method
Correlation ( ↑ )
MSE ( ↓ )
Epoch
Correlation ( ↑ )
MSE ( ↓ )
Epoch
Correlation ( ↑ )
MSE ( ↓ )
Epoch
C3TL (Ours)
0.491 ± 0.031
0.016 ± 0.005
0.25s
0.670 ± 0.202
0.019 ± 0.018
0.11s
0.777 ± 0.074
0.002 ± 0.000
1.7s
State
0.474 ± 0.047
0.018 ± 0.006
10s
0.608 ± 0.099
0.104 ± 0.028
5.4s
0.778 ± 0.069
0.004 ± 0.003
59s
CPA
0.008 ± 0.026
0.028 ± 0.011
8.9s
0.000 ± 0.000
0.589 ± 0.102
3.8s
0.13 ± 0.015
0.010 ± 0.003
87s
Mean Baseline
0.370 ± 0.041
62.453 ± 24.680
-
0.627 ± 0.151
4.602 ± 1.400
-
0.725 ± 0.083
1.684 ± 0.345
-
Closest Cell Baseline
0.319 ± 0.059
90.141 ± 30.987
-
0.523 ± 0.147
0.926 ± 0.641
-
0.630 ± 0.067
0.682 ± 0.300
-
Table 1: Performance of the models on perturbation datasets. The table shows the mean of the Pearson Correlation of the prediction of the individual algorithms with the true target for each gene as well as the Mean Squared Error. Metrics are reported as the mean ( ± standard deviation) over a 5-fold cross validation (4-fold for the Replogle dataset). Best and second best performing methods are colored in blue and gray respectively. In all experiments, metrics are only computed on interventions that have not been seen by the models in the target context.
Figure 2: Scatter plot of model predictions. Comparison of model outputs versus true target values for the Tahoe-100 dataset. The dashed line represents perfect prediction.
Figure 3: Sensitivity to data scarcity. Performance evaluation under varying training data availability. Panels (left to right) display regimes with 43, 20, 10, and 5 randomly sampled training contexts. The x-axis represents the number of interventions in the test context available for adaptation (implying the complement was used for testing). This ticks correspond to (from left to right) 80% , 50% , 30% , 20% , 10% , 5% , 1% of the data available in the test contexts in total. The y-axis displays the Pearson Correlation between prediction and target ( ± standard deviation across 5 test contexts). In all experiments, metrics are only computed on interventions that have not been seen by the models in the target context.
Appendix figures & tables1 asset
Supplementary material from the paper’s appendix.
Appendix
Tahoe-100
Replogle
Parse
Type
Compound
Genetic
Signaling
# Contexts
48
4
24
# Pert. / Context
1138
1677
90
Reference
[ 20 ]
[ 2 , 3 ]
[ 21 ]
Appendix
Table 2: Overview of datasets used for evaluation. The datasets cover different perturbation modalities including compound, genetic, and signaling perturbations.
Planning effective interventions in biological systems requires treatment-effect models that adapt to unseen biological contexts by identifying their specific underlying mechanisms. Yet single-cell perturbation datasets span only a handful of biological contexts, and existing methods cannot leverage new interventional evidence at inference time to adapt beyond their training data. To meta-learn a perturbation effect estimator, we present MapPFN, a prior-data fitted network (PFN) pre-trained on a synthetic biological prior with causal interventions, decoupling pre-training from limited wet-lab data. Unlike existing methods, MapPFN uses in-context learning to map a sequence of experiments to a post-perturbation distribution, enabling a single pre-trained model to adapt to new datasets and gene sets. Our results show that combining synthetic pre-training with fine-tuning on real single-cell data improves perturbation predictions across biological contexts. The code, model and data are available at https://marvinsxtr.github.io/MapPFN.
Marvin Sextro, Weronika Kłos, Gabriel Dernbach
Machine Learning Group, Technische Universität Berlin, Berlin, Germany · Berlin Institute for the Foundations of Learning and Data (BIFOLD) · Aignostics, Berlin, Germany +1
Single-cell perturbation prediction aims to infer how cells respond to unseen interventions and to achieve out-of-distribution (OOD) generalization, providing a computational route to understanding how perturbations reshape cellular programs over time. Existing machine learning methods have made important progress, but typically capture only one side of the response. Latent causal approaches seek mechanisms that support generalization and interpretation, yet often treat perturbation effects as static outcomes. Temporal models describe how gene expression changes across time, but usually do not explicitly recover the latent causal generative mechanisms driving these changes. In practice, perturbation effects are both latent and dynamical: interventions act through unobserved cellular programs, whose states evolve over time and give rise to observed expression profiles. Motivated by this view, we propose a latent dynamical causal generative model for single-cell perturbation data that jointly captures latent cellular programs, perturbation-conditioned mechanisms, and temporal evolution. We further provide an identifiability analysis showing that, under suitable conditions, the latent causal variables are recoverable up to standard equivalence classes. Guided by this analysis, we develop CITE-VAE, a learning framework for recovering latent cellular programs and their perturbation-driven dynamics from single-cell sequencing data. Experiments on Causal-3DIdent validate the theoretical results and the effectiveness of the proposed method in controlled settings. Additional experiments on real-world CRISPR-based single-cell perturbation data show improved generalization to unseen perturbations compared with state-of-the-art baselines, highlighting the practical robustness of our approach.
Wenkang Jiang, Yuhang Liu, Erdun Gao +3
AIML, Adelaide University Adelaide, Australia · Responsible AI Research Centre AIML, Adelaide University Adelaide, Australia · Monash University Melbourne, Australia +1
Predicting transcriptional responses to genetic perturbations is fundamental to functional genomics and therapeutic discovery. Recent deep learning models have shown promise in single-cell perturbation response prediction, but they typically generate each response in isolation, without explicitly leveraging related perturbations. We introduce PT-RAG (Perturbation-aware Two-stage Retrieval-Augmented Generation), a plug-in retrieval-and-conditioning module for generative cellular perturbation response. PT-RAG augments an existing perturbation-response backbone with learned access to related perturbation contexts. The key challenge is that relevance is not fixed in this setting: functionally related genes may elicit different effects across cell types. PT-RAG addresses this with a two-stage retrieval mechanism: GenePT-based semantic retrieval first identifies K candidate perturbations, after which a differentiable Gumbel-Softmax selector adaptively selects retrieved contexts conditioned on the control cell state, the query perturbation, and each candidate perturbation. We evaluate PT-RAG on two backbones, a STATE-style generator used as a frozen random reservoir and a fully trained scGPT, across cross-cell-type and cross-perturbation generalization tasks. PT-RAG consistently improves distributional similarity and often overall predictive quality; for example, on scGPT cross-cell-type results, the 2-Wasserstein distance drops by 5.9% relative to scGPT alone. The code to reproduce our experiments is available at https://github.com/difra100/PT-RAG_NIPS.
Andrea Giuseppe Di Francesco, Andrea Rubbi, Rishabh Jain +1