Biostars beta testing.
Question: transcription expression matrix for DESEQ2
0
Entering edit mode

Dear all, I am kind of new to DESEQ2 and I have an expression matrix csv that is consisted of ensemble transcription IDs (not gene ids) and their expression value in each condition. I am wondering if I can use this matrix directly as an input in the deseq2 pipeline.

Thanks Cheers

ADD COMMENTlink 17 months ago Duckula • 20 • updated 17 months ago andrew.j.skelton73 5.7k
Entering edit mode
1

Where did this expression matrix come from? - You can provide matrix input for DESeq2, but this must be integer counts of genes. For transcript level count estimates, you need to follow the tximport protocol, from either Salmon or Kallisto.

ADD REPLYlink 17 months ago
andrew.j.skelton73
5.7k
Entering edit mode
0

Hi, Thanks for your reply. These are actually raw transcript counts that are coming from bitSEQ. The values are integer. so would be possible using tximport package?

ADD REPLYlink 17 months ago
Duckula
• 20
1
Entering edit mode

bitSeq is more an all in one for differential transcript expression. You wouldn't be able to input those estimates into DESeq2, I'd recommend following one of these methodologies if you want to do differential gene expression or differential transcript expression from the same quantifications:

Salmon

  • Quantify your raw fastq with Salmon
  • Read in using tximport
  • Differential gene expression with DESeq2 or Limma
  • Convert Salmon output for compatibility with Wasabi
  • Run Sleuth for Differential Transcript Expression

Kallisto

  • Quantify your raw fastq with Kallisto
  • Read in using tximport
  • Differential gene expression with DESeq2 or Limma
  • Run Sleuth for Differential Transcript Expression
ADD COMMENTlink 17 months ago andrew.j.skelton73 5.7k
Entering edit mode
0

Dear Andrew, Thanks a lot for your reply. The thing is that in the paper that this dataset is originally published they have mentioned they have put the raw transcription data into deseq for further analysis. That is why I was wondering how can I do it and how should I proceed with this matrix. Do you have any idea that might be worth to try ?

ADD REPLYlink 17 months ago
Duckula
• 20
Entering edit mode
0

In this case, it's probably worth contacting the author, or looking at potentially more detailed methods in the supplementary. Keep in mind, "just because you can, doesn't mean you should" - It's worth critically assessing what DESeq2 requires as inputs, against what they've provided from bitSeq.

ADD REPLYlink 17 months ago
andrew.j.skelton73
5.7k

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.0