Question: Data with only 1 pair for DEG analysis
gravatar for bharata1803
3.9 years ago by
bharata1803430 wrote:


So, currently, I have 1 pair of normal-cancer data from RNA-seq. I already generated gene level read counts with Salmon. What I want to do is to check for a certain gene whether it is up or down regulated. So, is there any problem with just 1 pair of sample? What tools will be suitable for 1 pair sample like this (DESeq2, limma, edgeR,etc)?

Thank you for your answers and suggestions.

rna-seq deseq2 R • 1.2k views
ADD COMMENTlink modified 3.9 years ago by Chirag Parsania1.6k • written 3.9 years ago by bharata1803430
gravatar for Chirag Parsania
3.9 years ago by
Chirag Parsania1.6k
University of Macau
Chirag Parsania1.6k wrote:

HI Bharata,

It doesn't matter how many pairs/samples you have in differential analysis. What matters is whether you have biological replicates or not for each of your sample. If you have biological replicates then check the correlation of those replicates. If they are correlating well (corr coeff. >= 0.9) great. You can go for diff analysis. I have good exp with DESeq2. It works well for me. You can use other package as well (limma edgeR)and then make decision on the basis of what biology you are getting. whether you are getting relevant genes with significant p val or not in diff exp list. 





ADD COMMENTlink written 3.9 years ago by Chirag Parsania1.6k

Thank you. So, for this data, there is no replicates. I only have 1 fastq file from normal cell and 1 fastq file from cancer cell. That's why I am currently don't know what to do.

ADD REPLYlink written 3.9 years ago by bharata1803430

It's little difficult to convince people on your findings without replicates. Anyway you can run DESeq2 without replicates as well.  But in that case there is no meaning of your p-value and statistics. To overcome this you can filter highly up and down genes , confirm those in genome browser (for eg IGV ) and perform downstream analysis with those confidant diff. exp. genes. 



ADD REPLYlink written 3.9 years ago by Chirag Parsania1.6k

Thank you for your suggestion. I will try that. Also, I want to ask about replicates. So, if there are 3 cancer samples but from different patient/person, can those samples can be consider as replicates or it should come from same patient/person? Until now, the RNA-seq published via NCBI mostly compare several normal-cancer pair but from different patient/person. I know this replicates problem is very important for statistical significance so I want to understand more about this.

ADD REPLYlink written 3.9 years ago by bharata1803430

The answer is yes and no. Yes because if those cancers are same you can consider them as biological replicates. Again I would suggest you to calculate correlation between those patients samples. If samples don't correlate  though cancers are same (in different patients) you will be in trouble. You can not consider them as biological replicates. The reason for not correlating possibly technical or real biological phenomena. We don't know.  So far what I have observed is  people do sequencing of  same sample in duplicates (one patient) to check the technical variability vs biological phenomena. That will give you more confidant but definitely more expensive. 


ADD REPLYlink written 3.9 years ago by Chirag Parsania1.6k
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: 1182 users visited in the last hour