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
3 answers
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
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.
Kallisto requires sequencing reads and a file with target sequences.
Please see the kallisto manual for more information.
Log in to answer this question.
-
Differential expression genes
written by m90 •
-
Summarizing transcript level annotation from Trinity de novo assembly to gene level annotation
written by MB •
-
RNAseq differential expression analysis results: Kallisto vs STAR
written by marcos.georgiades.18 •
-
Differential expression analysis suggestions needed
written by sierraallinone •
-
Best practises RNAseq - Sleuth vs. NA's in annotation
written by BioBing
-
RNASeq reads from two bacterial species
written by alyamahmoud •
-
how does cuffdiff treats replicates with high std?
written by galozs •
-
Why does cuffdiff test differential isoform expression for genes with only one annotated transcript?
written by unsupervised_learner •
-
Forum: Need suggestions for rna-seq pipeline.
written by GouthamAtla
-
Rna Seq Differential Expression
written by Empyrean •
