Question: Treatment VS Control in Single Cell RNAseq analysis
0
gravatar for grayapply2009
7 months ago by
grayapply2009240
United States
grayapply2009240 wrote:

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?

ADD COMMENTlink modified 7 months ago by ATpoint44k • written 7 months ago by grayapply2009240
2
gravatar for ATpoint
7 months ago by
ATpoint44k
ATpoint44k wrote:

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

ADD COMMENTlink modified 7 months ago • written 7 months ago by ATpoint44k

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.

ADD REPLYlink written 7 months ago by grayapply2009240

Sorry, misread the question. The point is the same though. You have to use the unintegrated values for DE analysis, not the ones that the integration produces, see https://osca.bioconductor.org/multi-sample-comparisons.html#sacrificing-differences

ADD REPLYlink written 7 months ago by ATpoint44k

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

ADD REPLYlink written 7 months ago by grayapply2009240
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: 1209 users visited in the last hour
_