Estimating slide- and patch-level gene expression profiles from pathology images enables rapid and low-cost molecular analysis with broad clinical impact. Despite strong results, existing approaches treat gene expression as a mere slide- or spot-level signal and do not incorporate the fact that the measured expression arises from the aggregation of underlying cell-level expression. To explicitly introduce this missing cell-resolved guidance, we propose a Cell-type Prototype-informed Neural Network (CPNN) that leverages publicly available single-cell RNA-sequencing datasets. Since single-cell measurements are noisy and not paired with histology images, we first estimate cell-type prototypes-mean expression profiles that reflect stable gene-gene co-variation patterns.CPNN then learns cell-type compositional weights directly from images and models the relationship between prototypes and observed bulk or spatial expression, providing a biologically grounded and structurally regularized prediction framework. We evaluate CPNN on three slide-level datasets and three patch-level spatial transcriptomics datasets. Across all settings, CPNN achieves the highest performance in terms of Spearman correlation. Moreover, by visualizing the inferred compositional weights, our framework provides interpretable insights into which cell types drive the predicted expression. Code is publicly available at https://github.com/naivete5656/CPNN.
Figures & tables
Figure 1: (a) Gene expression is observed together with histology image as the sum of unobserved cell-level expressions. Expression of individual cells is not directly measurable. (b) Cell-level gene expression profiles are obtained from a public single-cell RNA sequencing (scRNA-seq) database, which contains noisy measurements without corresponding image information.
Figure 2: Illustration of concept of Cell-type Prototype-informed Neural Network (CPNN). We estimate gene expression based on cell-type prototypes by modeling the relationship between observed gene expression and cell-level gene expression.
Techniques
Bulk
Single-cell
Spatial
With image
✓
✗
✓
Exp. resolution
Tissue
Single-cell
Spot
Number of samples
High
Moderate
Low
Cost per sample
Low
Medium
High
Table 1: Comparison of bulk, single-cell, and spatial transcriptomics techniques.
Figure 3: Overview of our method. Gene expression is estimated based on cell-type prototype Tˉ and weight w(xi(n)) . The prototype supports the integration of covariance of gene expression among cells.
Method
References
BRCA
KIRC
LUAD
PCC
SCC
PCC
SCC
PCC
SCC
MIL
Max [ 39 ]
Pattern Recog. 2018
0.219
0.225
0.162
0.161
0.170
0.173
Mean [ 39 ]
Pattern Recog. 2018
0.247
0.263
0.211
0.217
0.203
0.213
AbMIL [ 14 ]
ICML 2018
0.280
0.310
0.258
0.260
0.261
0.270
ILRA [ 42 ]
ICLR 2023
0.291
0.310
0.261
0.276
0.245
0.283
S4MIL [ 8 ]
MICCAI2023
0.292
0.301
0.225
0.231
0.243
0.252
Table 2: Comparison performance of slide-level gene expression estimation on TCGA-BRCA, KIRC, and LUAD. MIL: multiple instance learning, MIR: multiple instance regression, PCC: Pearson correlation coefficient, and SCC: spearman correlation coefficient
Method
Ref.
CSCC [ 15 ]
Her2st [ 3 ]
STNet [ 13 ]
PCC
SCC
PCC
SCC
PCC
SCC
BLEEP [ 43 ]
NeurIPS 2023
-0.0050
-0.0024
0.0212
0.0321
-0.0051
-0.0062
EGN [ 45 ]
WACV 2023
0.0552
0.1340
0.0941
0.0479
0.0570
0.0272
STNet [ 13 ]
Nat. Biomed.
0.0856
0.1156
0.0849
0.0772
0.0514
0.0536
w/ Ours
Eng. 2020
0.1008 (+0.0152)
0.1662 (+0.0506)
0.0825 (-0.0024)
0.1111 (+0.0339)
0.0532 (+0.0018)
0.0623 (+0.0087)
TRIPLEX [ 6 ]
CVPR 2024
0.0895
0.1239
0.0900
0.0861
0.0606
0.0546
Table 3: Comparison of methods across CSCC, Her2st, and STNet datasets. Each ( +digit or -digit ) shows the difference from the baseline
Settings
PI
MC
U
R
SCC
w/o PI, MC, R
✗
✗
✓
✗
0.305
w/o MC, U, R
✓
✗
✗
✗
0.174
w/o U, R
✓
✓
✗
✗
0.248
w/o R
✓
✓
✓
✗
0.336
Ours
✓
✓
✓
✓
0.338
Table 4: Ablation study results. PI: prototype initialization of Tˉ , MC: modality correction using α and β , U: prototype update of Tˉ , R: regularization loss.
Coarse (8)
medium (29)
fine (49)
0.317
0.336
0.338
Table 5: Effect of cell-type label granularity on model performance (SCC).
Figure 4: Average weight for each subtype of BRCA slides. We visualized the estimated weights for the cancer-related prototypes.
Method
Setting
SCC
Baseline
10 clusters
0.297 (-0.041)
Baseline
40 clusters
0.300 (-0.038)
Baseline
70 clusters
0.303 (-0.035)
Average prototype
–
0.329 (-0.009)
w/o α and β
–
0.331 (-0.007)
Ours
cell-type informed
0.338
Table 6: Comparison of prototype construction strategies. SCC values are reported, with performance differences from the proposed method shown in parentheses.
Figure 5: Example of estimated proportion.
BRCA
KIRC
LUAD
0.338
0.370
0.253
Table 7: Deconvolution result of cell2location.
Figure 6: The parameter distribution for cell-type prototype on BRCA before and after training. Each scatter plot corresponds to a specific prototype, with individual points representing genes. The horizontal axis denotes parameter values before training, while the vertical axis represents values after training.
Histology-based single-cell spatial transcriptomics (ST) estimation aims to predict gene expression for individual cells from histopathological images and cell locations, reducing the need for costly single-cell ST measurements. Unlike existing histology-to-ST methods that mainly predict spot-level profiles for local regions containing multiple cells, this task requires modeling cell-to-cell expression variability, which is strongly structured by cell type. We propose Genomics-Guided Cell-Type-Specific Mixture-of-Experts (GC-MoE), which estimates cell-type probabilities with a routing network and softly combines cell-type-specific experts for gene expression prediction. To further encode cell-type-dependent gene programs, we introduce the Cell-Type-Specific Co-Expression-Aware Predictor (CAP), together with a lightweight Cell-to-Cell Interaction Attention (C2CA) module for neighboring-cell context. Experiments and ablations on public single-cell ST datasets show consistent improvements over existing single-cell and adapted spot-level baselines.
Predicting gene expression from H&E-stained histology images offers a scalable alternative to costly spatial transcriptomics, yet most existing methods operate at the spot level, where signals from multiple cells are aggregated and critical cellular heterogeneity is obscured. Extending this paradigm to single-cell resolution is non-trivial. Naively applying pathology foundation models faces a scale mismatch: their patch-level representations mix multiple cells, whereas per-cell cropping or resizing distorts morphology and removes local context. Conversely, segmentation-based models without strong pretrained visual encoders often lack the morphological representation capacity needed for accurate molecular prediction and inherit errors from imperfect cell boundary masks. Here, we present CELLO, an efficient end-to-end framework that performs a single pathology foundation model forward pass per image and uses grid sampling to extract location-specific features for all cells simultaneously. We further introduce a distance-decay cross-attention module that refines each cell representation using spatially biased local morphological context. Using 52 public Xenium-H&E pairs from HEST-1k that span 12 organs and approximately 10 million cells, CELLO improves the average predictive accuracy over the evaluated baselines while reducing the mean whole-slide inference time compared to DeepSpot2Cell, a 14.0x speed-up on average that excludes upstream cell segmentation. Our work establishes a scalable foundation for single-cell gene expression prediction from H&E images.
Zijun Gao, Chunbin Gu, Jinxi Xiang +2
The Chinese University of Hong Kong · Stanford University School of Medicine
Spatial transcriptomics can resolve gene expression at single-cell resolution, but it is costly, limited to targeted panels of a few hundred to a few thousand genes, and applicable to only a small number of samples. H&E imaging, by contrast, is cheap and collected routinely at scale. This makes predicting single-cell expression directly from morphology a practical way to bring molecular analysis to large tissue archives. We therefore present VOICE, a multimodal foundation model that predicts single-cell gene expression from H&E images using paired Xenium data. VOICE first aligns cell centered H&E morphology from a pathology foundation model with single-cell expression embeddings from a transcriptome foundation model, trained using contrastive learning over 23 million cells. Next it predicts expression through two branches. One branch directly regresses expression from morphology. The other branch retrieves measured expression from similar reference cells, recovering genes that do not have morphological signal. Because genes vary in morphological predictability, VOICE fuses the two branches with a per-gene weight. After training, VOICE generalizes to heldout patients, slides, and partially overlapping gene panels from Xenium, and it consistently outperforms prior single-cell expression prediction methods on seven metrics.
Xin Luo, Yicheng Tao, Haoxuan Zeng +6
University of Michigan Ann Arbor, Michigan, USA · Weill Cornell Medicine New York, New York, USA