Increasing evidence shows that epistasis, defined as interactive effects between genetic loci, may contribute to the missing heritability of cancer. However, systematic genome-wide epistasis identification in cancer remains challenging. Here, by leveraging genotype and clinical data from 380,983 samples in the UK Biobank, we identified 202,032 candidate epistatic single nucleotide polymorphism (epiSNP) pairs associated with cancer risk across 16 cancer types. Notably, multivariable Cox regression identified 123 epiSNP pairs with significant interaction effects on overall survival, suggesting that interaction-level genetic signals can provide prognostic information beyond individual SNP effects. Through functional analysis of the 202,032 candidate epiSNP pairs, we identified 7152 pairs supported by gene co-expression data and 12,326 pairs with protein–protein interaction (PPI) evidence. By mapping epiSNP pairs to corresponding gene pairs and then linking these gene pairs to drug–target databases, we identified 1040 epistatic gene pairs with FDA-approved drug–target records. Additionally, through KM survival analysis of the candidate epiSNP pairs, we detected 7068 pairs significantly associated with patient overall survival. Finally, we constructed an open-access database, EpiSNPdb, to facilitate cancer epistasis research.

The identification of cancer prognostic biomarkers is crucial for predicting disease progression, optimizing personalized therapies, and improving patient survival. Molecular biomarkers are increasingly being identified for cancer prognosis estimation. However, existing studies and databases often focus on single-type molecular biomarkers, deficient in comprehensive multi-omics data integration, which constrains the comprehensive exploration of biomarkers and underlying mechanisms. To fill this gap, we conducted a systematic prognostic analysis using over 10,000 samples across 33 cancer types from The Cancer Genome Atlas (TCGA). Our study integrated nine types of molecular biomarker-related data: single-nucleotide polymorphism (SNP), copy number variation (CNV), alternative splicing (AS), alternative polyadenylation (APA), coding gene expression, DNA methylation, lncRNA expression, miRNA expression, and protein expression. Using log-rank tests, univariate Cox regression (uni-Cox), and multivariate Cox regression (multi-Cox), we evaluated potential biomarkers associated with four clinical outcome endpoints: overall survival (OS), disease-specific survival (DSS), disease-free interval (DFI), and progression-free interval (PFI). As a result, we identified 4,498,523 molecular biomarkers significantly associated with cancer prognosis. Finally, we developed SurvDB, an interactive online database for data retrieval, visualization, and download, providing a comprehensive resource for biomarker discovery and precision oncology research.

Multi-nucleotide variants (MNVs) represent a crucial yet underexplored category of genetic variation. Despite previous studies highlighting the prevalence and potential biological impact of MNVs in populations, comprehensive identification and detailed functional annotation of MNVs remain challenging. Here, we develop MNVAnno, a toolbox for rapid identification and annotation of complex MNVs, and utilize it to identify 3,984,258 MNVs from 700,134 human samples, expanding the human MNV list to 8,199,654. Our analysis reveals that MNVs can not only lead to distinct amino acid changes from their constituent single-nucleotide variants, but also significantly impact the function of non-coding regions. Furthermore, through genome-wide association studies, we identify some MNVs associated with multiple cancers, and establish the Human MNV Database to facilitate MNV research. Our study emphasizes the importance of MNV annotation, broadens the human MNV landscape, and opens avenues for exploring genetic variation in phenotypes and diseases.