Differential Expression analysis workflow

Anna
2010-04-26
2013-06-05
  • Anna
    Anna
    2010-04-26

    Hi,

    I'm hoping that someone can help me, as I couldn't work out how to do
    this from the manual.  Would someone be able to give me the steps in a
    differential expression analysis?

    I have run tophat with the following command for each of my two solexa
    sequence.txt RNA-seq files seperately:

    tophat -solexa1.3-quals -p 2 -o 101/100315/tophat/
    ~/software/bowtie-0.12.2/indexes/m_musculus_ncbi37
    101/100315/s_2_sequence.txt &> 101/100315/tophat/tophat.out &

    I would like to get the expression levels for all Ensembl transcripts. 
    I have downloaded this gtf file from Ensembl,
    ftp://ftp.ensembl.org/pub/current_gtf/mus_musculus,

    However when I run
    cuffdiff ~/data/gtf/Mus_musculus.NCBIM37.57.gtf
    101/100315/tophat/accepted_hits.sam 95/100315/tophat/accepted_hits.sam
    &> cuffdiff.out &

    or
    cufflinks -G ~/data/gtf/Mus_musculus.NCBIM37.57.gtf 101/100315/tophat/accepted_hits.sam &>
    101/100315/tophat/cufflinks.out &

    I get the following error

    Error: duplicate GFF ID 'ENSMUST00000127664' (or exons too far apart)!

    I'm pretty sure I've misunderstood  the workflow, if someone could give me
    an overview of the steps and what gtf file I should be using that would
    be great.

    Many Thanks
    Anna

    (Cross posted to SeqAnswers)

     
  • mnqbvk I really liked your post.Really looking forward to read more. Will read on…

     
  • t62eEk I really liked your article.Much thanks again. Fantastic.

     
  • RoFe75 I value the blog. Keep writing.

     
  • VEE45H I really liked your blog article.

     
  • NaQT27 I cannot thank you enough for the article.Really thank you! Cool.

     
  • Hello world! my name is Tom,I come from the earth! insanity http://insanity05081.webs.com