DECONVersation leverages embedding representations from large-scale, LLM-based foundation models to perform deconvolution of bulk RNA-seq data. This takes advantage of the strengths of scFMs in faithfully representing transcriptomes, learning meaningful biological networks, and minimizing batch effect and noise. Currently, cell embeddings from Geneformer, Cell2Sentence, CellHermes, and scGPT are supported (+PCA and scVI for comparison).
DECONVersation enables end-to-end deconvolution through a set of easy-to-use functions. Embeddings can be extracted from both bulk and single-cell datasets, with single-cell embeddings used to construct robust signature matrices from .h5ad references. Cell type proportions are then estimated via NNLS directly in embedding space. Built-in benchmarking tools evaluate predictions against ground truth using RMSE and Pearson correlation, complemented by visualization utilities for assessing method performance. DECONVersation also supports testing and validation with in-built pseudobulk functions, model fine-tuning and attention extraction.
DECONVersation was benchmarked across 6 real bulk RNA-seq datasets with ground truths and 2 pseudobulk dataset, spanning diverse tissue types and experimental conditions, to evaluate deconvolution performance and generalizability.
Summary
Across 6 benchmarked real bulk and 2 pseudobulk datasets, we calculate overall RMSE and correlation coefficient alongside mean RMSE and correlation averaged across cell types. Fine-tuned Cell2Sentence and Geneformer-based embeddings both demonstrate consistent deconvolution performance across all datasets, with fine-tuned models outperforming their zero-shot counterparts in each case. Though zero-shot performance is already comparable to some common tools in the field, this highlights the benefit of fine-tuning (training models to predict cell type annotations from a single-cell reference). Among the tested tools, only DWLS R package achieves comparable performance to the fine-tuned embedding-based approaches available in DECONVersation.
| # | Dataset | Source | Ground Truth | Cell Type # |
|---|---|---|---|---|
| 1 | PBMC (Hoek) | PBMC | FACS | 5 |
| 2 | PBMC (Finotello) | PBMC | FACS | 5 |
| 3 | PBMC (Morandini) | PBMC | FACS | 5 |
| 4 | Cell Line Mixture (Cobos) | Cell Line Mixture | Mixture Count | 6 |
| 5 | Pre-Frontal Cortex (Huuki-Myers) | DLPFC | RNAScope/IF | 6 |
| 6 | Retina (Guo) | Retina | snRNA | 6 |
| 7 | HECA (Marečková) | Endometrium | pseudobulk | 9 |
| 8 | Lung tumor (Guimarães) | Lung | pseudobulk | 10 |
While DECONVersation itself is lightweight and easy to install with pip install DECONVersation, the various single cell foundation models themselves are not. In fact, due to dependency restrictions, they will never be compatible in the same python environment. DECONVersation works around this by detecting and only loading the available model(s). For each scFM model and package, users should consult the corresponding official installation guides.
We also provide conda env yaml files in the envs directory that are reproducibly operational on our hardware (NVIDIA L40S), each corresponding to one of the scFMs + DECONVersation. They can be installed with e.g. conda env create -f deconv_gf.yml. Apptainer def files are also included. Alternatively, we provide prebuilt CUDA12.9 singularity images for each scFM, with DECONVersation v0.1.0 installed, on Zenodo.
# python
import deconversation
res = deconversation.deconverse(
bulk_df = "bulk_rna.csv",
sig_df = "signature_matrix.csv", # or adata = "reference.h5ad", needs one or the other
model = "path_to/Geneformer-V2-316M",
temp_output_dir = "temp"
)
# command line
deconverse --help
deconverse --demo --model path_to_model
deconverse --bulk bulk_rna.csv --sig signature_matrix.csv --model path_to_model -o deconv_results.csv
deconverse --bulk bulk_rna.csv --adata reference.h5ad --cell_type_col celltype --model path_to_model
# with SIF mounting
singularity exec --nv \
--bind path_to_input_dir:/data:ro \
--bind path_to_output_dir:/results \
--bind path_to_model:/model:ro \
deconversation-geneformer.sif \
deconverse \
--bulk /data/bulk_rna.csv \
--sig /data/signature_matrix.csv \
--mode geneformer \
--model /model \
--temp-output-dir /results/temp \
--output /results/proportions.csv
- DECONVersation on bulk RNA-seq using Geneformer: Extract embeddings and deconvolute on bulk against a single cell reference.
- DECONVersation on pseudobulk using Geneformer: Validate deconvolution using pseudobulk data.
- Fine-tuning geneformer for cell type classification: Fine-tune geneformer for cell type classification
- Extracting attention weights using DECONVersation: Extract scFM attention weights (Geneformer, Cell2Sentence & cellHermes)
- Geneformer Transfer learning enables predictions in network biology
- Cell2Sentence Cell2Sentence: Teaching Large Language Models the Language of Biology
- CellHermes Language may be all omics needs: Harmonizing multimodal data for omics understanding with CellHermes
- scGPT scGPT: toward building a foundation model for single-cell multi-omics using generative AI