Hello guys! I would really appreciate if someone could help me with DE analysis. This is my challenge:
I have four conditions and I have five biological replicates for each condition. Performing a PCA (DESeq2) we can observe that most of the replicates don't cluster together. My first question is, can I simply analyze them using DESeq2 or edgeR or this would be wrong because of this replicate scenario? Second, are there any ways to filter out the genes in which the replicates are not good and keep only the ones that are consistent among the replicates of the same condition to run the DE analysis?
Thank for the help!
2 answers
If the 5 replicates are repeated in each condition (e.g. 5 cell lines), you can do a paired t-test, even if the PCA doesn't look nice. The paired t-test will correct for the individual bias of each replicate. Alternatively, if the 5 replicates are repeated in the conditions, you can perform batch correction on the replicates. As last resort, you can you SVA (surrogate variable analysis) to correct any batch effects.
You can do a paired t-test on Voom-limma transformed data. For edgeR and DEseq2 you will need to put the sample pairing in the linear model as covariate.
That's why I started my sentence with "if". Voom-limma is well accepted and shown to be as good as deseq2 or edger, there are numerous papers about that. Regarding the cell lines, I clearly wrote e.g. (=for example). Swbarnes2, if you know a better method, please suggest it.
I would hazard that the ground truth of your experiment is that your conditions do not affect RNA expression much. You can do DESeq2, but I predict it will return very few genes as DE.
You certainly cannot just throw away genes that don't agree with what you think the experiment should look like.
You might want to investigate and see if there is any experimental reason driving PC1.
You have to decide for yourself if all those genes expressed in only a couple of samples are real or not. You can filter away genes which, say, have fewer than 5 samples with > 10 counts, or whatever, but if that's all really biologically real, then I guess you should keep it. I assume you checked and all of these samples have a healthy number of total gene counts?
You will find many posts at http://support.bioconductor.org where the developer recommends some minimum prefiltering, e.g. at least 10 counts in at least 3 samples. True, that is (I think) not explicitely in the manual, but I think it makes good sense as these genes you show are probably not very reliable, at least not at that sample size. Three are zero and two are high, what do you believe => not easy to say with 5 samples in that group, one would need many more to see whether zero or high is the reproducible truth, hence the conservative approach would be the filtering. For example one of the groups must have at least 10 counts in either all of its samples or in the majority, there is no strict rule for that. edgeR::filterByExpr does this type of filtering with well-tested defaults and is group-aware, you could simple filter out from dds those genes flagged by the edgeR filter.+
Something like:
#/ assuming you had your factor of interest as "group":
raw_counts <- dds(counts)
keep <- filterByExpr(raw_counts, group=dds$group)
dds[keep,]
Log in to answer this question.
-
Hugely different results between edgeR and DESeq2
written by Gama313
-
Estimating power in RNA-seq differential expression
written by ams686 •
-
Quality Check of Samples and Biological Replicates
written by doinelpierrot •
-
Calculate p-value for experiment group with replicates, while control group without replicates
written by zhangdengwei
-
PCA on selected genes.
written by a.archana •
-
FeatureCounts from Galaxy to DESeq2 in RStudio
written by Newbieish •
-
experimental designs that allow/not-allow differential expression analyses
written by mattchen84 •
-
edgeR understanding statistics
written by vm.higareda •
-
Seq. depth for DE analysis
written by grant.hovhannisyan
-
Differential Expression Confidence Values Based On Another Sample Replicate Or Re-Sampling
written by Leszek
