Question: Hill shape p-value histogram after differential miRNA expression analysis with DESeq2
0
gravatar for L_K
3 months ago by
L_K0
L_K0 wrote:

Hi everyone

I have an issue with my differential miRNA expression analysis project. I have two different comparisons treatment1 (T1) vs control and treatment2 (T2) vs control and I performed standard differential expression analysis with DESeq2 and the analysis resulted in 5 significantly differentialy expressed miRNA between T1 and Control and 0 significant miRNA between T2 and controls (padj<0.05). Thus I started to investage what could have gone wrong and thus I plotted the p-value histograms which looked as follows:

histogramT1vC

histogramT2vC

histogram1T1vC

histogramT2vC

Alright, thus i followed the instructions on following website https://www.huber.embl.de/users/klaus/Teaching/DESeq2Predoc2014.html where they use fdrtool to estimate the true null model variance. I did this for both datasets T1vC and T2v C according to the link above. output of the result for one dataset (the second one looks similar): fdrtools null model sd https://paste.pics/4SXS5

So I redid the testing etc acording to the embl protocol from huber lab and now to my main question: new p-value histogram is more what one would expect new pvalue hist https://paste.pics/4SXU4 However out of about 550 tested miRNAs I get 68 significantly differentially expressed miRNAs for 1 dataset T1vC and about 120 significant miRNAs for T2vC. And this number just seems way too high, no? I would never expect 1/5 of tested miRNAs to be significantly differentially expressed, or would you?

Which results should & can I trust more?

I appreciate every help, as I am not a statistican...

Best regards!

edit: Thanks for the input. Added 'fdrtool' as a tag

mirna deseq2 fdrtool • 114 views
ADD COMMENTlink modified 3 months ago by Robert Castelo2.3k • written 3 months ago by L_K0
Answer: Hill shape p-value histogram after differential miRNA expression analysis with D
0
gravatar for Michael Love
3 months ago by
Michael Love23k
United States
Michael Love23k wrote:

This is a good question for the fdrtool author Bernd Klaus. Can you add this package as a tag?

ADD COMMENTlink modified 3 months ago • written 3 months ago by Michael Love23k

True, thanks for the input. I added 'fdrtool' as a tag!

ADD REPLYlink written 3 months ago by L_K0

True, thanks for the input. I added 'fdrtool' as a tag!

ADD REPLYlink written 3 months ago by L_K0
Answer: Hill shape p-value histogram after differential miRNA expression analysis with D
0
gravatar for Lorena Pantano
3 months ago by
Boston
Lorena Pantano90 wrote:

Hi,

Although I am not an expert in the use of fdrtool, normally I check I few things to understand the results of differential expression analysis.

Just to add some information to your question, How many replicates do you have in each group?

Did you plot the new genes are significant to know whether there is an outlier in your groups (an example explained in your link as the cause of the uphill pvalue distribution)? I would inspect a bunch of them, and see if you catch up some patterns.

Did you see whether the fold change are big or small for the significant genes?

As for the number of miRNA deregulated, it really depends on the biology involved. I had similar results for fibrosis study where the cells change a lot. So, I wouldn't discard that number because it is big. Many miRNA are involved in similar biological processes, and they are regulated in clusters. This doesn't mean that suddenly things make sense, but I think looking at the specific miRNA expression to catch up some obvious pattern, the fold change and the biological context should help to decide whether something is off or not.

Thanks for the question, it is very interesting for me.

Cheers

ADD COMMENTlink written 3 months ago by Lorena Pantano90
Answer: Hill shape p-value histogram after differential miRNA expression analysis with D
0
gravatar for Robert Castelo
3 months ago by
Robert Castelo2.3k
Spain/Barcelona/Universitat Pompeu Fabra
Robert Castelo2.3k wrote:

hi, i'd say that such a right-skewed distribution is often a product of some confounding factor that is not being adjusted. You can find a discussion on this in this old thread from the former Bioconductor mailing list.

there has been also a more recent discussion about hill-shaped histograms and the hint that Mike was referring to about fdrtool, in this other thread.

cheers, robert.

ADD COMMENTlink modified 3 months ago • written 3 months ago by Robert Castelo2.3k
Please log in to add an answer.

Help
Access

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