Alternative splicing downstream of EMT enhances phenotypic plasticity and malignant behaviour in colon cancer
Abstract
Phenotypic plasticity allows carcinoma cells to transiently acquire the quasi-mesenchymal features necessary to detach from the primary mass and proceed along the invasion-metastasis cascade. A broad spectrum of epigenetic mechanisms is likely to cause the epithelial-to-mesenchymal (EMT) and mesenchymal-to-epithelial (MET) transitions necessary to allow local dissemination and distant metastasis. Here, we report on the role played by alternative splicing (AS) in eliciting phenotypic plasticity in epithelial malignancies with focus on colon cancer. By taking advantage of the coexistence of subpopulations of fully epithelial (EpCAMhi) and quasi-mesenchymal and highly metastatic (EpCAMlo) cells in conventional human cancer cell lines, we here show that the differential expression of ESRP1 and other RNA-binding proteins (RBPs) downstream of the EMT master regulator ZEB1, alters the AS pattern of a broad spectrum of targets including CD44 and NUMB, thus resulting in the generation of specific isoforms functionally associated with increased invasion and metastasis. Additional functional and clinical validation studies indicate that both the newly identified RBPs and the CD44s and NUMB2/4 splicing isoforms promote local invasion and distant metastasis and are associated with poor survival in colon cancer. The systematic elucidation of the spectrum of EMT-related RBPs and AS targets in epithelial cancers, apart from the insights in the mechanisms underlying phenotypic plasticity, will lead to the identification of novel and tumor-specific therapeutic targets.
Data availability
The RNA-sequencing data from this study have been submitted to the Gene Expression Omnibus (GEO) database under the accession number GSE192877. Other data referenced in this study are publicly available and can be accessed from the GEO using GSE154927, GSE154730 and Synapse using identifier syn2623706 .
-
CD44s and CD44v6 overexpressed RNAseq profiles of colon cancer cell lines HCT116 and SW480NCBI Gene Expression Omnibus, GSE192877.
-
Transcriptome datasets of the Quaking (QKI) gene knock-down human oral squamous cell carcinoma (OSCC) cellsNCBI Gene Expression Omnibus, GSE154730.
Article and author information
Author details
Funding
China Scholarship Council (201806300047)
- Tong Xu
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Ethics
Animal experimentation: All mice experiments were implemented according to the Code of Practice for Animal Experiment in Cancer Research from the Netherlands Inspectorate for Health Protections, Commodities and Veterinary Public Health. Permit number AVD1010020171344.
Reviewing Editor
- Lynne-Marie Postovit, University of Alberta, Canada
Publication history
- Preprint posted: March 1, 2022 (view preprint)
- Received: July 19, 2022
- Accepted: November 7, 2022
- Accepted Manuscript published: November 8, 2022 (version 1)
- Version of Record published: November 18, 2022 (version 2)
Copyright
© 2022, Xu et al.
This article is distributed under the terms of the Creative Commons Attribution License permitting unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 948
- Page views
-
- 160
- Downloads
-
- 0
- Citations
Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Cancer Biology
In embryonal rhabdomyosarcoma (ERMS) and generally in sarcomas, the role of wild-type and loss or gain-of-function TP53 mutations remains largely undefined. Eliminating mutant or restoring wild-type p53 is challenging; nevertheless, understanding p53 variant effects on tumorigenesis remains central to realizing better treatment outcomes. In ERMS, >70% of patients retain wild-type TP53, yet mutations when present are associated with worse prognosis. Employing a kRASG12D-driven ERMS tumor model and tp53 null (tp53-/-) zebrafish, we define wild-type and patient-specific TP53 mutant effects on tumorigenesis. We demonstrate that tp53 is a major suppressor of tumorigenesis, where tp53 loss expands tumor initiation from <35% to >97% of animals. Characterizing three patient-specific alleles reveals that TP53C176F partially retains wild-type p53 apoptotic activity that can be exploited, whereas TP53P153D and TP53Y220C encode two structurally related proteins with gain-of-function effects that predispose to head musculature ERMS. TP53P153D unexpectedly also predisposes to hedgehog expressing medulloblastomas in the kRASG12D-driven ERMS-model.
-
- Cancer Biology
The presence of lymph node metastasis (LNM) affects treatment strategy decisions in T1NxM0 colorectal cancer (CRC), but the currently used clinicopathological-based risk stratification cannot predict LNM accurately. In this study, we detected proteins in formalin-fixed paraffin-embedded (FFPE) tumor samples from 143 LNM-negative and 78 LNM-positive patients with T1 CRC and revealed changes in molecular and biological pathways by label-free liquid chromatography tandem mass spectrometry (LC-MS/MS) and established classifiers for predicting LNM in T1 CRC. An effective 55-proteins prediction model was built by machine learning and validated in a training cohort (N=132) and two validation cohorts (VC1, N=42; VC2, N=47), achieved an impressive AUC of 1.00 in the training cohort, 0.96 in VC1 and 0.93 in VC2, respectively. We further built a simplified classifier with nine proteins, and achieved an AUC of 0.824. The simplified classifier was performed excellently in two external validation cohorts. The expression patterns of 13 proteins were confirmed by immunohistochemistry, and the IHC score of five proteins was used to build an IHC predict model with an AUC of 0.825. RHOT2 silence significantly enhanced migration and invasion of colon cancer cells. Our study explored the mechanism of metastasis in T1 CRC and can be used to facilitate the individualized prediction of LNM in patients with T1 CRC, which may provide a guidance for clinical practice in T1 CRC.