Gene set analysis for time-to-event outcome: comparison of a new approach based on the generalized Berk–Jones statistic with existing methods in presence of intra gene-set correlation

Source: PubMed Central Open Access, NCBI / U.S. National Library of Medicine

Abstract Gene set analysis evaluates the collective impact of groups of genes on an outcome of interest, such as disease occurrence. By incorporating biological knowledge through predefined gene sets, this approach enhances the interpretability of results and improves statistical power compared with gene-wise analyses. In the context of time-to-event data, existing methods are limited and fail to account for potentially strong correlations within gene sets. Given the strong performance of the Generalized Berk-Jones () statistic, which effectively incorporates correlation within the test statistic, we adapted this method to the time-to-event framework using a Cox model. We then compared its performance with established methods, including the,,,, and. We further benchmarked these methods in two different real-world datasets: gliomas and breast cancer. Our proposed method,, shows an overcontrol of Type I error, leading to reduced statistical power compared with other methods in numerical studies particularly when the number of genes is greater than or equal to the number of observations. Theandgenerally exhibited the highest power, except in very high-correlation settings for the, while thecould not adjust for confounders in current implementations.