This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Treatment VS Control in Single Cell RNAseq analysis

I have scRNA-seq data for 4 groups Treatment 1, Control 1, Treatment 2, Control 2.

My goal is to compare gene expression in each cell type between Treatment 1 and Control 1, between Treatment 2 and Control 2, and between Treatment 1 and Treatment 2.

Should I perform QC, Normalization, and clustering on the combined 4 scRNA-seq datasets and then do DE analysis in each identified cluster?

Any suggestions?

single cell rnaseq de analysis treatment control
• 5,701 views • link
updated by ATpoint ATpoint • written by ngs101.com ngs101.com

Hello, I have more or less same experimental setting and was wondering how you proceeded with your analysis? Can I reach its GitHub repository to see codes or paper to see comparisons on plots?

• 0 views • link
written by kilcdincer • kilcdincer

Open a new question, with details.

• 0 views • link
written by ATpoint ATpoint

1 answer

I suggest you go through https://osca.bioconductor.org/ extensively before diving into analysis. Read all paragraphs, it answers 99.9% of all standard questions.

  • QC, initial normalization and selection of informative genes typically is done per dataset separately.
  • Then an integration step to merge all datasets, e.g. with fastMNN or the Seurat anchoring framework.
  • clustering based on the integrated dataset
  • DE analysis between clusters based on the unintegrated dataset with the clusters defined with the integrated dataset

For an extensive discussion see the linked workflow. See also C: Batch correction in scRNA-seq data via fastMNN | pro/contra

• 0 views • link
written by ATpoint ATpoint

Thank you for your prompt reply, ATpoint. What do you mean by "DE analysis between clusters based on the unintegrated dataset with the clusters defined with the integrated dataset"? I'm interested in the DE analysis between treatment groups and control groups in the same cluster. I don't care about the comparison between different clusters except when I'm identifying cluster marker genes.

• 0 views • link
written by ngs101.com ngs101.com

Thank you, my friend. That was helpful. I'm reading it now.

• 0 views • link
written by ngs101.com ngs101.com

Log in to answer this question.

More posts like this