I'm using Homer to perform a differential analysis between KO and WT with ChiP Seq data.
I follow this protocol:
- makeTagDirectory with the different bam files
- Run annotatepeaks.pl in this order: KO KO WT WT
- Run getDiffExpression.pl tablecounts-txt KO KO WT WT (following the same order than I have in my table obtained in the step before)
At the end I want to plot the Fold Change of KO over WT (ratio KO/WT) but I don't fully understand how this is calculated by HOMER with getDiffExpression because the output is :
Output Stats KO vs. WT:
Total Genes: X
Total Up-regulated in WT vs. KO: X ( X %) [log2fold>1, FDR<0.05]
Total Dn-regulated in WT vs. KO: X ( X %) [log2fold<-1, FDR<0.05]
This make me feel confuse about the sense/direction of the ratio (Fold change). Is it doing the KO / WT or WT / KO?
I think understand this is very important to be able to interpret correctly my plots and I don't find the documentation of Homer clear in this regard. Also, I took a look to different papers that use this methods, but didn't explain in methods how they exactly apply getDiffExpression.pl, if they used the order mutant mutant wildtype wildtype or the other way around.
Can someone help me?
Thank you very much!