Intronic polyadenylation (IpA) usually leads to changes in coding region of an mRNA, and its implication in diseases has been recognized, though at its very beginning status. Conveniently and accurately identifying IpA is of great importance for further evaluating its biological significance. Here, we developed IPAFinder, a bioinformatic method for the de novo identification of intronic poly(A) sites and their dynamic changes from standard RNA-seq data. Applying IPAFinder to 256 pan-cancer tumor/normal pairs across six tumor types, we discovered 490 recurrent dynamically changed IpA events, some of which are novel and derived from cancer-associated genes such as TSC1, SPERD2, and CCND2. Furthermore, IPAFinder revealed that IpA could be regulated by factors related to splicing and m6A modification. In summary, IPAFinder enables the global discovery and characterization of biologically regulated IpA with standard RNA-seq data and should reveal the biological significance of IpA in various processes.

Somatic single nucleotide variants (SNVs) in cancer genome affect gene expression through various mechanisms depending on their genomic location. While somatic SNVs near canonical splice sites have been reported to cause abnormal splicing of cancer-related genes, whether these SNVs can affect gene expression through other mechanisms remains an open question. Here, we analyzed RNA sequencing and exome data from 4,998 cancer patients covering ten cancer types and identified 152 somatic SNVs near splice sites that were associated with abnormal intronic polyadenylation (IPA). IPA-associated somatic variants favored the localization near the donor splice sites compared to the acceptor splice sites. A proportion of SNV-associated IPA events overlapped with premature cleavage and polyadenylation events triggered by U1 small nuclear ribonucleoproteins (snRNP) inhibition. GC content, intron length and polyadenylation signal were three genomic features that differentiated between SNV-associated IPA and intron retention. Notably, IPA-associated SNVs were enriched in tumor suppressor genes (TSGs), including the well-known TSGs such as PTEN and CDH1 with recurrent SNV-associated IPA events. Minigene assay confirmed that SNVs from PTEN, CDH1, VEGFA, GRHL2, CUL3 and WWC2 could lead to IPA. This work reveals that IPA acts as a novel mechanism explaining the functional consequence of somatic SNVs in human cancer.

Retinoblastoma (RB) is the most common pediatric intraocular malignancy and is a serious vision- and life-threatening disease. The biallelic mutation of the retinoblastoma gene RB1 is the initial event in the malignant transformation of RB, but the exact molecular mechanism is still unclear. E2F transcription factors can be activated by RB1 loss of function and lead to uncontrolled cell division. Among E2F family numbers, E2F1 has higher expression abundance than E2F2 and E2F3 in RB clinical samples. By integrating E2F1 ChIP-seq data, RNA-seq profiling from RB samples and RNA-seq profiling upon E2F1 knockdown, together with pathway analysis, literature searching and experimental validation, we identified Cyclin-dependent kinases regulatory subunit 2 (CKS2) as a novel regulator in regulating tumor-associated phenotypes in RB. CKS2 exhibited aberrantly higher expression in RB. Depletion of CKS2 in Y79 retinoblastoma cell line led to reduced cell proliferation, delayed DNA replication and decreased clonogenic growth. Downregulation of CKS2 also slowed tumor xenograft growth in nude mice. Importantly, reversed expression of CKS2 rescued cancer-associated phenotypes. Mechanistically, transcription factor E2F1 enhanced CKS2 expression through binding to its promoter and CKS2 regulated the cancer-associated PI3K–AKT pathway. This study discovered E2F1/CKS2/PTEN signaling axis regulates malignant phenotypes in pediatric retinoblastoma, and CKS2 may serve as a potential therapeutic target for this disease.