Question: analysis of multiple batches and multiple treatment conditions in scRNA-seq
0
gravatar for Bogdan
18 days ago by
Bogdan850
Palo Alto, CA, USA
Bogdan850 wrote:

Dear all, please would you advise us on the following :

for example, shall we have 2 batches of scRNA-seq data of 2 conditions (WT_batch1, WT_batch2, DISEASE_batch1, DISEASE_batch2), would the following approach be statistically legitimate in order to account/correct for the batch effect :

1 -- use CCA or MNNcorrect to account for the batch effects

2 -- followed by TSNE and network_based_clustering, in order to place correctly the cells in CORRECT CLUSTERS

3 -- and perform differential expression (with wilcoxon test, limma, edgeR, etc) between the CLUSTERS

We know that CCA or MNNcorrect only place the cells in more "correct" clusters after batch correction, and do NOT provide a batch - corrected expression value.

In this case, considering for instance cluster_0, could we combine :

a -- the matrix of cells : normalized_expression in cluster-0 in WT_batch1

with the matrix of cells : normalized_expresion in cluster-0 in WT_batch2

(let's call this matrix WT_batch1_batch2)

b -- the matrix of cells : normalized_expression in cluster-0 in DISEASE_batch1

with the matrix of cells : normalized_expresion in cluster-0 in DISEASE_batch2

(let's call this matrix DISEASE_batch1_batch2)

c -- and use limma or edgeR or DESeq2 on WT_batch1_batch2 versus DISEASE_batch1_batch2 in order to get the differential expression

we would prefer to combine the batches into WT_batch1_batch2 and respectively, DISEASE_batch1_batch2, as, sometimes, the number of cells in a cluster may be small (ie less than 200 cells)

or if there is any other approach that you'd recommend ..

thank you,

bogdan

ADD COMMENTlink modified 16 days ago by sofiagreen7221130 • written 18 days ago by Bogdan850
3
gravatar for sofiagreen72211
16 days ago by
sofiagreen7221130 wrote:

Seurat has recently introduced a new workflow to deal with batch correction:

https://satijalab.org/seurat/v3.0/integration.html https://github.com/satijalab/seurat/issues/1224

Hope this helps.

ADD COMMENTlink written 16 days ago by sofiagreen7221130

thank you Sofia ... I m glad that the new Seurat pipelines address the question that we had.

ADD REPLYlink written 16 days ago by Bogdan850
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1118 users visited in the last hour