QIAGEN powered by

Bonferroni and FDR multiple testing corrections too strict for differential expression analyses

Issue description

Calculations for the Bonferroni and FDR multiple testing corrections in Differential Expression for RNA-Seq and Differential Expression in Two Groups used an inflated value for the number of tests performed. In affected software versions, this number included both the number of tests performed and the number of untestable genes/transcripts, i.e. those with NaN as an expression value.

The impact of this is that fewer differentially expressed transcripts/genes are reported after applying a p-value cut-off based on these corrections than should have been the case. The missing transcripts/genes will be those nearest to the p-value cut-off, i.e. they will not be the most significantly differentially expressed transcripts/genes.

Even when using affected software:

Recommendations

For results generated using affected software versions, upgrading to an unaffected version is generally recommended. If analyses are then re-run using the tools Differential Expression for RNA-Seq and Differential Expression in Two Groups, an increased number of differentially expressed transcripts/genes may be reported due to the changes made to the multiple testing correction methods.

Affected versions

A fix was implemented in CLC Genomics Workbench 12.0.3 and CLC Genomics Server 11.0.2.