Differential expressiong with two groups in BASiCS
2
0
Entering edit mode
@muadabdelhay-14075
Last seen 3.8 years ago

I have two groups of cells with two predefined phenotypes. There are 30 in each group (60 in total) and I wanted to use BASiCS to perform differential expression between the two groups. How do I go about it?

 

The vignette suggests providing two chains into the DE-function of BASiCS. Does that mean that I have to run MCMC for each subgroup separately and then do DE providing each group separately?

 

Or do I run MCMC on the whole group of cells (60) and then split them into two chains, one for each group? I haven't found a was to subset a chain yet (maybe it is not possible...?)

BASiCS single-cell bayesian differential expression • 1.4k views
ADD COMMENT
1
Entering edit mode
CataVallejos ▴ 10
@catavallejos-13762
Last seen 4.5 years ago
MRC Human Genetics Unit & The Alan Turi…

Yes, you should run `BASiCS_MCMC` separately for each group of samples and then use those as the input of `BASiCS_TestDE`. Let me know if you have any issues while running the analysis.

 

 

ADD COMMENT
0
Entering edit mode

I ran it as you suggested. The results have similarities to what I get from DESeq2.

 

Is there something I need to take into account if I used the WithSpikes = FALSE option (as my data lacks spike-ins)?

ADD REPLY
0
Entering edit mode
CataVallejos ▴ 10
@catavallejos-13762
Last seen 4.5 years ago
MRC Human Genetics Unit & The Alan Turi…

Indeed, I will expect that DESeq2 and BASiCS report similar results when assessing changes in mean expression between the groups. Said this, the main advantage of using BASiCS is the differential variability testing (through over-dispersion or, more recently, residual over-dispersion https://www.biorxiv.org/content/early/2017/12/21/237214). Moreover, downstream analysis should not be affected by the lack of spikes (as seen in the biorXiv pre-print, the spikes and non-spikes models lead to very similar posterior inference). 

 

 

ADD COMMENT

Login before adding your answer.

Traffic: 972 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6