Question: Problem in analyzing Rice RNA-Seq data
0
gravatar for Whoknows
4.9 years ago by
Whoknows750
Tehran,Iran
Whoknows750 wrote:

Hi friends 
I have project on Rice RNA-Seq, in this project i have only 1 replicate for each of 3 conditions.(N,M,S)
I should to say, Number of mapped reads (using tophat2) for conditions are: 

N= 14,640,080 
M= 14,764,894 
S= 14,935,162
And you can see "S" has more mapped reads( in contrast with other ~17000 and ~29000).

I analyzed with cuffdiff (with blind mode) which generates:

N_vs_M=30
N_vs_S=350
M_vs_S=319

i know N_vs_M=30 has a problem in gene count and not fine.

I can not find out the problem , i also used DESeq,DESeq2 and edgeR but they could not show any significant DE genes.

I literature review and found same project where done by BGI company, i could not gain any information about the analysis process. They produced these columns :

geneID
geneLength
A-Expression
B-Expression
A-RPKM
B-RPKM
log2 Ratio(B/A)
Up-Down-Regulation(B/A)
P-value
FDR


I have 2 questions:
1- Which software can produce these output
2- How can i generate P-value and FDR without replicate? (i mean software)

Thanks.

deseq deseq2 tophat cufflinks htseq • 2.1k views
ADD COMMENTlink modified 4.9 years ago • written 4.9 years ago by Whoknows750

Hello pcsam.2008!

It appears that your post has been cross-posted to another site: SEQanswers

This is typically not recommended as it runs the risk of annoying people in both communities.

ADD REPLYlink written 4.9 years ago by Devon Ryan91k

Do you have another idea about this paper "Transcriptomic Analysis of Rice (Oryza sativa) Developing Embryos Using the RNA-Seq Technique" , they have only 1 sample per condition and also said:

We applied the R package DEGseq to identify DEGs with the
random sampling model based on the read count for each gene at
different developmental stages [33]. We used ‘‘FDR <=0.001 and
the absolute value of log2Ratio >=1’’ as the threshold for judging
the significance of each gene expression difference. More stringent
criteria with a smaller FDR and bigger fold-change value can be
used to identify DEGs

What do you think about this issue? - I mean 1 sample and FDR 0.001 and FC=2

How is it possible?

thanks..,,

 

 

ADD REPLYlink written 4.9 years ago by Whoknows750

This should probably be its own post.

Having said that, that paper should have been rejected.

ADD REPLYlink written 4.9 years ago by Devon Ryan91k
0
gravatar for Istvan Albert
4.9 years ago by
Istvan Albert ♦♦ 81k
University Park, USA
Istvan Albert ♦♦ 81k wrote:

You need to contact the institution/company/individual that gave you the analysis for details. 

That being said if you have no replication across conditions then there is very little chance of being able to meaningfully compare these conditions. 

ADD COMMENTlink written 4.9 years ago by Istvan Albert ♦♦ 81k

thanks Istvan 

ADD REPLYlink modified 4.9 years ago • written 4.9 years ago by Whoknows750
Please log in to add an answer.

Help
Access

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