Alternative polyadenylation (APA) plays a role in gene expression regulation generally by shortening of 3'UTRs (untranslated regions) upon proliferative signals and relieving microRNA-mediated repression. Owing to high proliferative indices of triple negative breast cancers (TNBCs), we hypothesized APA to cause 3'UTR length changes in this aggressive subgroup of breast cancers. Our probe-based meta-analysis approach identified 3'UTR length alterations where the significant majority was shortening events (∼70%, 113 of 165) of mostly proliferation-related transcripts in 520 TNBC patients compared with controls. Representative shortening events were further investigated for their microRNA binding potentials by computational predictions and dual-luciferase assay. In silico-predicted 3'UTR shortening events were experimentally confirmed in patient and cell line samples. To begin addressing the underlying mechanisms, we found CSTF2 (cleavage stimulation factor 2), a major regulator of 3'UTR shortening to be up-regulated in response to epidermal growth factor (EGF). EGF treatment also resulted with further shortening of the 3'UTRs. To investigate the contribution of CSTF2 and 3'UTR length alterations to the proliferative phenotype, we showed pharmacological inhibition of the EGF pathway to lead to a reduction in CSTF2 levels. Accordingly, RNAi-induced silencing of CSTF2 decreased the proliferative rate of cancer cells. Therefore, our computational and experimental approach revealed a pattern of 3'UTR length changes in TNBC patients and a potential link between APA and EGF signaling. Overall, detection of 3'UTR length alterations of various genes may help the discovery of new cancer-related genes, which may have been overlooked in conventional microarray gene expression analyses.

Certain aspects of diagnosis, prognosis, and treatment of cancer patients are still important challenges to be addressed. Therefore, we propose a pipeline to uncover patterns of alternative polyadenylation (APA), a hidden complexity in cancer transcriptomes, to further accelerate efforts to discover novel cancer genes and pathways. Here, we analyzed expression data for 1045 cancer patients and found a significant shift in usage of poly(A) signals in common tumor types (breast, colon, lung, prostate, gastric, and ovarian) compared to normal tissues. Using machine-learning techniques, we further defined specific subsets of APA events to efficiently classify cancer types. Furthermore, APA patterns were associated with altered protein levels in patients, revealed by antibody-based profiling data, suggesting functional significance. Overall, our study offers a computational approach for use of APA in novel gene discovery and classification in common tumor types, with important implications in basic research, biomarker discovery, and precision medicine approaches.

Epidermal growth factor receptor (EGFR) has critical roles in epithelial cell physiology. Over-expression and over-activation of EGFR have been implicated in diverse cancers, including triple-negative breast cancers (TNBCs), prompting anti-EGFR therapies. Therefore, developing potent therapies and addressing the inevitable drug resistance mechanisms necessitates deciphering of EGFR related networks. Here, we describe Sorting Nexin 3 (SNX3), a member of the recycling retromer complex, as a critical player in the epidermal growth factor (EGF) stimulated EGFR network in TNBCs. We show that SNX3 is an immediate and sustained target of EGF stimulation initially at the protein level and later at the transcriptional level, causing increased SNX3 abundance. Using a proximity labeling approach, we observed increased interaction of SNX3 and EGFR upon EGF stimulation. We also detected colocalization of SNX3 with early endosomes and endocytosed EGF. Moreover, we show that EGFR protein levels are sensitive to SNX3 loss. Transient RNAi models of SNX3 downregulation have a temporary reduction in EGFR levels. In contrast, long-term silencing forces cells to recover and overexpress EGFR mRNA and protein, resulting in increased proliferation, colony formation, migration, invasion in TNBC cells, and increased tumor growth and metastasis in syngeneic models. Consistent with these results, low SNX3 and high EGFR mRNA levels correlate with poor relapse-free survival in breast cancer patients. Overall, our results suggest that SNX3 is a critical player in the EGFR network in TNBCs with implications for other cancers dependent on EGFR activity.