catavallejos / BASiCS

BASiCS: Bayesian Analysis of Single-Cell Sequencing Data. This is an unstable experimental version. Please see http://bioconductor.org/packages/BASiCS/ for the official release version
84 stars 17 forks source link

Processing raw data with CellBender prior to BASiCS analysis #233

Closed binyaminZ closed 2 years ago

binyaminZ commented 2 years ago

Hi, (this is not a technical issue rather a usage consultation...) I saw a few recent studies (like this) that perform ambient RNA background removal with CellBender prior to their analysis and see substantial improvement in quantification quality. If I understand correctly, your model doesn't account for this type of noise, which is assumed to be global (not cell-specific) and uniform (unlike the "unexplained noise" component in your model). CellBender also deals with barcode swapping. This is particularly important when dealing with single-nucleus RNAseq data, as in the study mentioned above and in my case. Do you think it would be beneficial to run CellBender prior to BASiCS? Or it might introduce some bias that could interfere with correct BASiCS performance? Thanks, Binyamin

alanocallaghan commented 2 years ago

Hi Binyamin,

I haven't read the CellBender paper in detail, but I don't think there would be any issues using BASiCS on the corrected count matrix that it creates. As you say it should handle extra noise sources that BASiCS doesn't account for, so it seems like a good idea in principle. If you do use it feel free to post any issues or feedback that result from it Cheers Alan