Question: (Closed) confusing cuffdifff results
0
gravatar for F
3.2 years ago by
F3.4k
Iran
F3.4k wrote:

Hi,

I got fully confused and frustrated by watching the result

I had 6 samples paired-end reads by tophat2

I used ensemble gtf and pre-index bowtie2

Left reads:
     Input   : 19667353
     Mapped  : 18672982 (94.9% of input)
      of these:  1127159 ( 6.0%) have multiple alignments (939 have >20)
Right reads:
     Input   : 19667353
     Mapped  : 17146234 (87.2% of input)
      of these:  1019710 ( 5.9%) have multiple alignments (848 have >20)
91.1% overall read mapping rate.

Aligned pairs: 16994149
  of these:  1010648 ( 5.9%) have multiple alignments
         25673 ( 0.2%) are discordant alignments
**86.3% concordant pair alignment rate.**

then I used cuffdiff to get differential expressed genes

This is my experiment details

In my experiment, before high temperature treatment (48 c), Arabidopsis cotyledons were pre-heated for 2, 8 and 24 hours. In the other hand I have cotyledons keeping in room temperature for 2, 8 and 24 hours, I mean without pre-heating these cotyledons directly face the high temperature. I don have any replication

Something like below

  • pre-heated-2h
  • pre-heated-8h
  • pre-heated-24h
  • unpre-heated-2h
  • unpre-heated-8h
  • unpre-heated-24h

I tried 3 cuffdiff syntaxs,but in ALL, in the gene_exp.diff file the significant column was totally NO, not only this file, the significant column wherever is NO , even when I labeled 6 samples, in the output I have only to columns for samples

cuffdiff -o diff_out -b genome.fa -p 8 \
    -L Pre-2h,unPre-2h,Pre-8h,unPre-8h,Pre-24h,unPre-24h \
    -u merged_asm/merged.gtf \
    -T pre-heated-2h.bam \
    unpre-heated-2h.bam pre-heated-8h.bam \
    unpre-heated-8h.bam pre-heated-24h.bam \
    unpre-heated-24h.bam #syntax1

cuffdiff -o diff_out -b genome.fa -p 8 \
    -L Pre-2h,unPre-2h,Pre-8h,unPre-8h,Pre-24h,unPre-24h \
    -u merged_asm/merged.gtf \
    pre-heated-2h.bam unpre-heated-2h.bam \
    pre-heated-8h.bam unpre-heated-8h.bam \
    pre-heated-24h.bam unpre-heated-24h.bam #syntax2

cuffdiff -o diff_out -b genome.fa -p 8 \
    -L pre,unpre \
    -u merged_asm/merged.gtf \
    pre-heated-2h.bam, pre-heated-8h.bam, \
    pre-heated-24h.bam unpre-heated-2h.bam \
    unpre-heated-8h.bam, unpre-heated-24h.bam #syntax3

How possible that there is no any differential expressed genes in my experiment??? :(

Thank you

EDIT (@RamRS): Code folder for better display

myposts rna-seq next-gen • 861 views
ADD COMMENTlink modified 7 months ago by RamRS21k • written 3.2 years ago by F3.4k
2

Hello Fereshteh!

Questions similar to yours can already be found at:

We have closed your question to allow us to keep similar content in the same thread.

If you disagree with this please tell us why in a reply below. We'll be happy to talk about it.

Cheers!

PS: Why open another question with the same content?

ADD REPLYlink modified 7 months ago • written 3.2 years ago by RamRS21k
1

Its very funny to try different commands to check differentially expressed genes. Lets say the command is wrong but you get differentially expressed genes ( some junk ), what do you do ? 

ADD REPLYlink modified 3.2 years ago • written 3.2 years ago by geek_y9.4k

you mean i did wrong commands? if so please tell me the right one

i should find some sense result to give to my supervisor...you can't imagine how much yesterday i cried, yesterday was one of the worst days in my life. she told she won't give me the data sets because i could not infer any results then without the data sets how i could defend my dissertation?

i know this is not sense to express such topic in this thread but you may know why i am asking silly questions

ADD REPLYlink written 3.2 years ago by F3.4k

What Goutham means is that you should judge your results by your approach, not the other way around. Data science will yield results if you keep tweaking your question, but along the way, your question itself will lose meaning. You'll get a flawed answer with no biological significance.

Your approach must be to ensure your command aligns with what you know of the data. And remember, no result is not research failing - it's research yielding a different yet important result - that there is no significant answer to your question. That does not make the question invalid, it just means that other questions need to be explored.

ADD REPLYlink written 3.2 years ago by RamRS21k

thank you Ram,

actually when I saved the gene_exp.diff file as txt and inspected for any significancy, I found more than 700 yes. now another step started...anyway many thanks

ADD REPLYlink modified 3.2 years ago • written 3.2 years ago by F3.4k
Please log in to add an answer.
The thread is closed. No new answers may be added.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2112 users visited in the last hour