Hello there, I'm kinda new in the field and i'm trying for the first time to analyze scRNA-seq data downloaded from GEO database. I downloaded data form GSE273040 which is a count matrix with raw and normalized counts from every sample. Since i want to proceed the analysis with DESeq2, should I use just the normalized counts, or should I instead find I way to use the raw counts? Moreover, I saw many tutorials saying to create a Seurat obj and then use it to perform the analysis, should I do it or not?
Thanks for the help!
2 answers
A "raw matrix" usually contains lots of empty droplets and lots of droplet rows, therefore you should apply an empty droplet detection or ambient rna detection method.
People usually put only count matrix to GEO database. You can identify whether a matrix is normalized or not by checking the values.
Deseq2 requires only "count" matrix with integer values.
You should apply filtering on cells and genes before doing deseq2. And also if you are also interested in cell type differences you should also identify cell types first.
And also ı dont think you can apply deseq2 on single cell data. You should do pseudobulking first.
Generally speaking, if you're new to something and expert software exists, then use it. For single-cell data there are three big open-source frameworks that are established:
Bioconductor (R): https://bioconductor.org/books/release/OSCA/
Seurat (R): https://satijalab.org/seurat/articles/get_started_v5_new
ScanPy (Python): https://scanpy.readthedocs.io/en/stable/tutorials/index.html
...and for a general read: https://www.sc-best-practices.org/preamble.html
For DESeq2 use raw counts, but there is a lot of upstream processing (like clustering, annotation, subsetting as carolofharvest suggested already) before you do DE analysis on single-cell data so read above first.
Log in to answer this question.