Question: Problem in analyzing Rice RNA-Seq data
gravatar for Whoknows
6.1 years ago by
Whoknows840 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:


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 :

log2 Ratio(B/A)

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)


deseq deseq2 tophat cufflinks htseq • 2.3k views
ADD COMMENTlink modified 6.1 years ago • written 6.1 years ago by Whoknows840

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 6.1 years ago by Devon Ryan97k

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?




ADD REPLYlink written 6.1 years ago by Whoknows840

This should probably be its own post.

Having said that, that paper should have been rejected.

ADD REPLYlink written 6.1 years ago by Devon Ryan97k
gravatar for Istvan Albert
6.1 years ago by
Istvan Albert ♦♦ 85k
University Park, USA
Istvan Albert ♦♦ 85k 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 6.1 years ago by Istvan Albert ♦♦ 85k

thanks Istvan 

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


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