Python deconvolution tools for bulk RNA-seq data

Question

[AlexStar]

I’m looking for a Python package to perform cell-type deconvolution on bulk RNA-seq datasets. I work with a large number of samples from diverse tissues and sources, and my goal is to develop a Python-based pipeline that can handle this analysis on the fly, to make my life easier.

In the past, I’ve used xCell, but since it isn’t available in Python and isn’t a true deconvolution method that provides actual cell-fraction estimates, I’m now looking for alternatives. I’ve recently transitioned to using Python, so I’d really appreciate recommendations for a robust, well maintained deconvolution tool, along with guidance on how to use it effectively.

Thanks!

[Arup Ghosh]

Check the Bulk2Single option from OmicVerse framework.

Answer 1

[ATpoint]

Imo you are approaching this the wrong way. R packages for deconvolution that are heqavily used and cited do exist, for example MuSiC and BisqueRNA exist. Also, interoperatability packages such as r2py and reticulate do exist. Hence, it is easy to run these packages from inside python, without reinventing the wheel. It comes down most of the time to running a single command, so easy to implement.

[AlexStar]

I see, yes that makes sense. Could you recommend any reliable and accurate package I could use ? ATpoint

MuSiC 2.0 might be a good choice ? The thing is I usually use bulk datasets with either Normal samples or Tumor samples, I rarely use datasets with both conditions in the same dataset.