Hi. I have two libraries for treatment and control. I have no replicate. Can I use TPM for each library as input for DESeq2? How can I calculate log 2 fold change by GFOLD? Is there any software for calculate it? How can I calculate reliable fold change and obtain DEGs when I have not any replicates?Is it correct that I use of TPM as input for DESeq2? What is your suggestion?
thanks.
So DESeq2/edgeR and limma works on count data and not TPM for any of the statistics that will build the model for DE analysis. So do not use TPM. Having said that it does not work for one sample per condition for an estimation of DESeq2. The have always maintained to not trust those answers. TPM is mostly used for visualization and not for estimating DEGs. GFOLD works but as far as I know it is still using count data. Take a look here. Other DE tools will not be advisable for usage for DE testing.
Thanks vchris. Its meaning that i use count of read that produced by miRDeep2 as input for DESeq2? What is seq(norm) that produced by miRDeep2?When I dont have any replicate, what can I do for calculating differentially expressed miRNAs?
Do you have bam file? you can use featurecounts of HT-Seq to pull the count data from your bam file. I have never used miRDeep2 so cannot tell you about it. What I understand from seq(norm) as you state is probably normalized by library size produced by mirDeep2. You need to check that in the manual. However, the link I sent you earlier shows how you can pull out counts from bam file with GFOLD as well and feed it to GFOLD for DE analysis.
Bam file produced by????I have sam file from bowtie 2. You are right about seq(norm). Now can I use seq(norm) for DESeq2 input? Its very confusing. Thanks a lot
you are missing the point, please read the manuals of GFOLD, the link I posted, it is your work to understand how and why you need to use tools. The link I have given for GFOLD explicitly comes with a tutorial how to use from sam files count data and use GFOLD for DE analysis without replicates.
For DESeq2 I will never use it for 1 vs 1 and to top it no you cannot use seq(norm) data in DESeq2 since DESeq normalization on its own relies on geometric norm which in your case will not work. So better to avoid it. However, if you are hell-bent on using DESeq2, please read the link below what ATPoint has provided. Without reading and understanding, it does not help.
Ok. Thank you very much.