Abstract Motivation O-GlcNAcylation, a dynamic post-translational modification regulated by O-GlcNAc transferase (OGT) and O-GlcNAcase (OGA), influences critical biological processes and is dysregulated in cancers. Direct measurement of O-GlcNAcylation dysregulation is challenging due to its instability and low-throughput nature, limiting large-scale studies. However, the regulatory simplicity of this system and the availability of transcriptomic data enable inference of dysregulation from OGT and OGA expression. Results We introduce a nonparametric kernel density estimation-based approach to quantify O-GlcNAcylation dysregulation using joint OGT and OGA expression. In simulated datasets with varied expression patterns and controlled dysregulation levels, our method consistently outperformed canonical metrics in quantifying dysregulation. In TCGA data from six cancer types, inferred regulation scores were significantly lower in cancer samples (0. 25–0. 30 vs. 0. 49–0. 51) and showed strong distributional differences (Kolmogorov–Smirnov p-values 5. 95e-11; D-statistics 0. 31) compared to those from healthy samples. The scores also allow for accurate classification of cancer status (AUROC: 0. 71–0. 75) and generalized well to external datasets without retraining. This transcriptomics-based framework offers a scalable approach for interpretable quantification of O-GlcNAcylation dysregulation in cancer. Availability and implementation The code and datasets used in this study are freely available at https: //github. com/wonder-ai/O-GlcNAcylationProject under an open-source license.
Stojšin et al. (Fri,) studied this question.