Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • ChrisL
    replied
    No need for Cufflinks or Cuffcompare if you are happy to stick within an existing genome annotation like RefSeq, UCSC or Ensembl. Just take the sam files from Tophat and go straight to Cuffdiff. There is a little trick to get a gtf file that is suitable to use with Cuffdiff with tss id, etc. Just use Cuffcompare and feed it the reference annotation (RefSeq, Ensembl, etc) twice and it will give you a gtf file that is suitable.

    Hope this helps.

    Leave a comment:


  • poisson200
    replied
    Thanks Olgy,
    I am sure your intelligence and knowledge of other things is vast.

    ===============================

    Dear Enrico,
    thanks for your answer. I since found out this is the answer to my questions, which are (very similar to yours with the novel and known case)

    1) Known transcripts only; Find differentially expressed known Refseq transcripts between case and control data (no novel transcripts wanted);
    a. Run TopHat on both fastq read files using a GFF file of Refseq
    b. Run cuffdiff on both sam outputs, using a Refseq GTF file

    2) Known and novel transcripts; Find differentially expressed transcripts, both novel and known Refseq
    a. Run TopHat on both fastq (without GFF)
    b. Run cufflinks without GTF file on each SAM file
    c. Run cuffcompare with Refseq.GTF sam1.gtf and sam2.gtf
    d. Run cuffdiff on the combined.GTF and sam files

    Thanks for the link too.

    Kind regards,

    J

    Leave a comment:


  • Enrico Palazzo
    replied
    Originally posted by poisson200 View Post
    Dear Bioinformaticians,
    The experiment; RNA-seq, transcriptomics; we simply want to find differentially expressed genes between 2 Illumina samples, case vs. control. No novel transcripts needed, just want to find differentially expressed known genes and/or known transcript isoforms.

    Do I make the right assumption that there is no need to use cufflinks or cuffcompare?
    E.g. Cuffcompare does;
    1) Compare your assembled transcripts to a reference annotation
    2) Track Cufflinks transcripts across multiple experiments
    ????

    Our proposed approach to finding known differentially expressed genes/isoforms:
    1) Run Tophat with a Refseq GTF file on both case and control
    2) Run cuffdiff on both sam files generated from 1 and the GTF file used in 1 as well.

    Please let me know if this is or is not a valid approach?

    Thank you,

    J.
    Hi J,
    I'm new to this field and this is just a guess:
    I think you have to run cuffcompare on the two GTF files produced in step 2 and your reference annotation file:
    cuffcompare -r Mus_musculus.test.gtf -R -o prefix transcripts1.gtf transcripts2.gtf

    Then you run cuffdiff on the combined output GTF file (from cuffcompare) and the 2 SAM files which had been generated by TopHat in step 1:
    cuffdiff -o diff_out/ combined.gtf accepted_hits1.sam accepted_hits2.sam

    DE genes/isoforms will be in 0_1_gene_exp.diff/0_1_isoform_exp.diff
    but it seems that there is no output for coding sequences as mentioned in an other thread http://seqanswers.com/forums/showthread.php?t=4989. If you have trouble with duplicate entries in your reference annotation filles: check this: http://seqanswers.com/forums/showthread.php?t=3493

    Leave a comment:


  • Olgy
    replied
    I would really want to help you, but you know that my knowlage about programming are very small

    Olgy Wolgy

    Leave a comment:


  • poisson200
    replied
    There are 49 views but nobody knows the answer?


    Or is this is a stupid question?

    Leave a comment:


  • poisson200
    started a topic To cuffcompare or not to cuffcompare

    To cuffcompare or not to cuffcompare

    Dear Bioinformaticians,
    The experiment; RNA-seq, transcriptomics; we simply want to find differentially expressed genes between 2 Illumina samples, case vs. control. No novel transcripts needed, just want to find differentially expressed known genes and/or known transcript isoforms.

    Do I make the right assumption that there is no need to use cufflinks or cuffcompare?
    E.g. Cuffcompare does;
    1) Compare your assembled transcripts to a reference annotation
    2) Track Cufflinks transcripts across multiple experiments
    ????

    Our proposed approach to finding known differentially expressed genes/isoforms:
    1) Run Tophat with a Refseq GTF file on both case and control
    2) Run cuffdiff on both sam files generated from 1 and the GTF file used in 1 as well.

    Please let me know if this is or is not a valid approach?

    Thank you,

    J.
    Last edited by poisson200; 07-29-2010, 04:16 AM.

Latest Articles

Collapse

  • seqadmin
    Latest Developments in Precision Medicine
    by seqadmin



    Technological advances have led to drastic improvements in the field of precision medicine, enabling more personalized approaches to treatment. This article explores four leading groups that are overcoming many of the challenges of genomic profiling and precision medicine through their innovative platforms and technologies.

    Somatic Genomics
    “We have such a tremendous amount of genetic diversity that exists within each of us, and not just between us as individuals,”...
    05-24-2024, 01:16 PM
  • seqadmin
    Recent Advances in Sequencing Analysis Tools
    by seqadmin


    The sequencing world is rapidly changing due to declining costs, enhanced accuracies, and the advent of newer, cutting-edge instruments. Equally important to these developments are improvements in sequencing analysis, a process that converts vast amounts of raw data into a comprehensible and meaningful form. This complex task requires expertise and the right analysis tools. In this article, we highlight the progress and innovation in sequencing analysis by reviewing several of the...
    05-06-2024, 07:48 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, 05-24-2024, 07:15 AM
0 responses
13 views
0 likes
Last Post seqadmin  
Started by seqadmin, 05-23-2024, 10:28 AM
0 responses
17 views
0 likes
Last Post seqadmin  
Started by seqadmin, 05-23-2024, 07:35 AM
0 responses
20 views
0 likes
Last Post seqadmin  
Started by seqadmin, 05-22-2024, 02:06 PM
0 responses
10 views
0 likes
Last Post seqadmin  
Working...
X