Dear all,
referring to scRNA-seq analysis, many thanks to Aaron Lun and Helena Crowell, for writing the following tutorials on OSCA and MUSCAT :
https://osca.bioconductor.org/multi-sample-comparisons.html
'd appreciate also having your suggestions on the following case of scRNAseq analysis (that includes Seurat 3.1, beside SimpleSingleCell/BioC packages) ; the question is : what analysis strategy would you recommend ?
Shall we have 4 batches of scRNA-seq data of these experiments :
WTbatch1, WTbatch2, Abatch1, Abatch2
WTbatch3, WTbatch4, Bbatch3, Bbatch4
What is the optimal way to analyze the data-sets ? Several analysis strategies are possible :
STRATEGY A : using MUSCAT (as outlined above)
STRATEGY B :
- to use CELLRANGER AGGR (with NORMALIZATION = TRUE) on :
WTbatch1, WTbatch2 : to produce WTbatch1_2
Abatch1, Abatch2 : to produce Abatch1_2
WTbatch3, WTbatch4 : to produce WTbatch3_4
Bbatch3, Bbatch4 : to produce Bbatch3_4
- and, to follow the descriptions of SimpleSingleCell/fastMNN or SEURAT pipelines with fastMNN (or with LIGER, HARMONY, CONOS) on WTbatch12, WTbatch34, Abatch12, Bbatch34 :
or,
STRATEGY C.
to join all the data (WTbatch1, WTbatch2, Abatch1, Abatch2, WTbatch3, WTbatch4, Bbatch3, Bbatch4) in a large MATRIX
and to follow the "pseudo-bulk" approach described in https://osca.bioconductor.org/multi-sample-comparisons.html
or,
Any other analysis strategy ? Any suggestions, comments would be very welcome ! Thanks a lot !
-- bogdan
Cross-posted: https://www.biostars.org/p/418552/
Also to bioinfo SE: https://bioinformatics.stackexchange.com/questions/11248/analysis-of-multiple-time-points-2-replicates-of-scrna-seq