An extension of the Wilcoxon-Mann-Whitney test for analyzing RT-qPCR data

Classical approaches for analyzing reverse transcription quantitative polymerase chain reaction (RT-qPCR) data commonly require normalization before assessing differential expression (DE). Normalization often has a substantial effect on the interpretation and validity of the subsequent analysis steps, but at the same time it causes a reduction in variance and introduces dependence among the normalized outcomes. These effects can be substantial, however, they are typically ignored. Most normalization techniques and methods for DE focus on mean expression and are sensitive to outliers. Moreover, in cancer studies, for example, oncogenes are often only expressed in a subsample of the populations during sampling. This primarily affects the skewness and the tails of the distribution and the mean is therefore not necessarily the best effect size measure within these experimental setups. In our contribution, we propose an extension of the Wilcoxon-Mann-Whitney test which incorporates a robust normalization, and the uncertainty associated with normalization is propagated into the final statistical summaries for DE. Our method relies on semiparametric regression models that focus on the probability P{Y ≤ Y'}, where Y and Y' denote independent responses for different subject groups. This effect size is robust to outliers, while remaining informative and intuitive when DE affects the shape of the distribution instead of only the mean. We also extend our approach for assessing DE for multiple features simultaneously. Simulation studies show that the test has a good performance, and that it is very competitive with standard methods for this platform. The method is illustrated on two neuroblastoma studies.