Hi all,

I have a scRNA-seq dataset, which has 6 patients and each patient has 2 sample types (normal, tumor). So I have **12 folders** and each folder contain its scRNA data: `barcodes.tsv.zip`, `features.tsv.zip`, `matirx.mtx.tsv.zip`. I would like to use this data to practice and learn the `Seurat (version 4.0)`  r package workflow. 

I would like to process all the 12 samples into cluster and downstream analysis. For the preprocessing, I now have no idea which point I should merge all the datasets. 

Should I **read and QC each folder** (e.g patient1-normal -> patient1-tumor -> patient2-normal .... so on), and then merge all data? I am still learning the scRNA seq, I believe each data will have different number of genes and cells, I need to normalize it after merging all the data.

Thank you for your help and advice.


Yes, you should QC each sample individually (you can automate this of course but be sure to manually inspect each QC plot). This is at least what I do. There are approaches like 3xMAF to automatically filter data, and I use them in support of what my eyes see in the QC plots, but I personally think that the human brain is an awesome diagnostic tool so just looking at the plots makes a lot of sense to see whether the automated cutoffs are actually meaningful. With skewed data they might not be. And then for clustering you should integrate your data. Please read the Seurat integration vignette. I suggest though that you download some of the Seurat example data to practice with because with 12 samples it will take both time and memory. Better to go easy with a smaller dataset or just two patients to start with. As an alternative to Seurat with lots of code to follow is https://osca.bioconductor.org/, I recommend it as documentation is (for me) much clearer. It has a section on multi-sample integration as well.