RNA-seq: differentially expressed genes without genome reference
3
2
Entering edit mode
5.8 years ago
guiu ▴ 60

Hello,

I have RNA-seq data for an organism for which the transcriptome is available, but not the genome (it is not assembled yet). I am interested in differential expression analysis at gene level.

I have been reading about RSEM, CuffDiff, Sailfish, Kallisto, HTSeq, but I didn't get a clear idea of which is the best approach to get differentially expressed genes.

What is your recommendation? Thanks

 

RNA-Seq • 2.9k views
ADD COMMENT
6
Entering edit mode
5.8 years ago
Rob 5.0k

I'd recommend using Sailfish or Salmon with tx-import to get fast accurate gene-level results (example). tx-import will allow you to import these results directly into DESeq2 (or prepare them easily for most other downstream DE packages).

ADD COMMENT
3
Entering edit mode
5.8 years ago

You have two choices in the absence of a reference genome 1. Annotate you transcriptome with tools like Trinotate, do the mapping taking into account this information (Bowtie2, but, kallisto), and then the DE analysis (limma, DESeq2, edgeR, NOISeq) 2. Do the mapping and DE analysis and annotate afterwards

ADD COMMENT
1
Entering edit mode
5.8 years ago
Kamil ★ 2.1k

If you know which transcripts belong to a gene, then you can sum their TPM values to get the gene's TPM.

Consider using Kallisto to get quick and accurate results.

Next, use your preferred method for differential expression analysis: limma, DESeq, edgeR, etc.

ADD COMMENT
0
Entering edit mode

Do you I need annotation for that ? or only transcriptome is enough ?

ADD REPLY
0
Entering edit mode

Kallisto requires sequencing reads and a file with target sequences.

Please see the kallisto manual for more information.

ADD REPLY

Login before adding your answer.

Traffic: 1008 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6