Abstract
Background: Gene expression and alternative splicing (AS) can promote cancer development via complex mechanisms. We aimed to identify and verify the hub AS events and splicing factors associated with the progression of colorectal cancer (CRC).
Methods: RNA-Seq data, clinical data, and AS events of 590 CRC samples were obtained from the TCGA and TCGASpliceSeq databases. Cox univariable and multivariable analyses, KEGG, and GO pathway analyses were performed to identify hub AS events and splicing factor/spliceosome genes, which were further validated in five CRCs.
Results: In this study, we first compared differentially expressed genes and gene AS events between normal and tumor tissues. Differentially expressed genes were different from genes with differentially expressed AS events. Prognostic analysis and co-expression network analysis of gene expression and gene AS events were conducted to screen five hub gene AS events involved in CRC progression: EPB41L2, CELF2, TMEM130, VCL, and SORBS2. Using qRT-PCR, we also verified that the gene AS events SORBS2 were downregulated in tumor tissue, and gene AS events EPB41L2, CELF2, TMEM130, and VCL were upregulated in tumor tissue. The genes whose mRNA levels were significantly related to the five hub gene AS events were significantly enriched in the GO term of cell division and Notch signaling pathway. Further coexpression of gene AS events and alternative splicing factor genes revealed NOVA1 as a crucial factor regulating the hub gene AS event expression in CRC. Through in vitro experiments, we found that NOVA1 inhibited gene AS event SORBS2, which induced the migration of CRC cells via the Notch pathway.
Conclusion: Integrated analysis of gene expression and gene AS events and further experiments revealed that NOVA1-mediated SORBS2 promoted the migration of CRC, indicating its potential as a therapeutic target.
Introduction
Colorectal cancer is one of the most frequently diagnosed malignant cancers and the second leading cause of cancer-related deaths worldwide (). Despite advances in diagnosis and treatment, the clinical outcomes for individual cases remain unsatisfactory (). Unfortunately, approximately 50% of CRC patients will develop recurrence after tumor resection, most of them within 2 years (). Accumulated evidence in the past decades has indicated that recurrence and metastasis are complex and tightly regulated processes involved in numerous genetic alternations (). Therefore, exploring the biological mechanisms involved in CRC metastasis is vital in developing targeted therapies and improving prognosis of patients.
Although the total number of protein-coding genes in the human genome is nearly 20,000, the transcriptome is tremendously more complex (). A single RNA precursor can be spliced in different arrangements to synthesize distinct mRNA and protein variants with different structures and functions. This may be one of the main reasons for the increased complexity of biological processes (). The basic splicing patterns include alternate acceptor sites (AA), alternate donor sites (AD), exon skipping (ES), mutually exclusive exons (ME), alternate terminators (ATs), alternate promoters (APs), and retained introns (RI). Although exon array-based studies could help reveal the complexity of AS, the limitations of this technology are also apparent (). With the rapid development of high-throughput technology, the functions of AS events can be widely and largely used in CRC. Using the Cancer Genome Atlas (TCGA) project data, performed a systematic analysis to reveal critical AS events involved in CRC progression. Similar studies have also been performed in other types of solid tumors, such as lung cancer, esophageal carcinoma, and gastric cancer (; ; ). Specifically, the splicing factor SRSF6 directly binds to exon23 of ZO-1 mRNA to produce the aberrant ZO-1 isoform, which functions as a tumor-promoting gene in colorectal cancer (). NOVA1, another splicing factor, promotes the inclusion of exons in the reverse transcriptase domain of hTERT to produce full-length hTERT transcripts, which accelerates non-small cell lung cancer growth ().
Notch signaling is a primordial, evolutionarily conserved cell fate determination pathway involved in cancer biology, such as cancer stem cells, angiogenesis, and tumor immunity (; ). Activation of Notch signaling can induce cancer cell proliferation, invasion, and metastasis, such as colorectal cancer, gastric cancer, tongue cancer, and cervical cancer (; ; ; ; ). Different gene isoforms regulate Notch signaling. presented evidence that the predominant TP63 isoform (ΔNp63α) promotes cellular quiescence via activation of Notch signaling. reported that an isoform of the “a” subunit of V-ATPase (named as a2V-ATPase) regulated processing of Notch receptor and altered Notch signaling in breast cancer. QKI-5, an RNA-binding protein, regulates the alternative splicing of NUMB to suppress cell proliferation and prevent Notch pathway activation ().
Due to the complex mutual regulation mediated by gene expression and AS events (; ), integrated analysis of gene expression and AS helps to better understand the mechanisms involved in cancer progression (; ; ). In this study, we performed an integrated analysis of gene expression and gene AS events and identified critical gene AS events in CRC progression. Using GO analysis and a series of in vitro experiments, we further explored the potential mechanisms regulating gene AS events and AS events that induce CRC progression.
Materials and Methods
Data Extract
The transcriptome, AS events, and the prognostic data of 590 colorectal cancers were downloaded from the TCGA1 and TCGASpliceSeq2 databases, which included 50 colorectal cancers with matched normal tissue. Sixty-six genes of splicing factor genes were obtained from the SpliceAid 2 database3.
Patients
Five consecutive, unselected CRC tissues and paired normal tissues were collected from surgical resections of CRC patients at Ruijin Hospital. None of the patients received adjuvant therapy before the operation. CRC tumors were histologically confirmed by two experienced pathologists independently. Written informed consent was obtained from all the patients with CRC. All aspects of this study were approved by the Research Ethics Committee of Shanghai Jiao Tong University. The UICC TNM classification system was used to determine the tumor stage. All tissue samples were stored in liquid nitrogen. The patient information is illustrated in Supplementary Table 2.
Cell Lines, Cell Culture, and Transient Transfection
The CRC cell lines were purchased from the American Type Culture Collection and cultured in RPMI 1640 supplemented with 10% fetal bovine serum and 1% penicillin–streptomycin at 37°C with 5% CO2. Cancer cells were grown to 60–70% confluence, and the siRNAs were transfected using Lipofectamine 3000 according to the manufacturer’s protocol (Thermo Fisher Scientific, L3000015). After 48 h of transfer, RNA was collected, and then at 72 h, proteins were collected, and functional experiments were performed. The siRNA sequences used in this study were purchased from Bioegene (Shanghai) and are listed in Supplementary Table 3.
Western Blot Analysis
The total protein of cancer cells was extracted using RIPA buffer supplemented with 1 mM phenylmethylsulfonyl fluoride. Protein concentration was measured using the BCA Protein Assay Kit (Cat. No. P0012S. Beyotime, China). Western blotting was performed according to our previous protocols (). The following primary antibodies were used: anti-GAPDH (ABCAM, ab8245, 1:1,000), anti-NOVA1 (SAB, 27594, 1:1,000), anti-NOTCH1 (CST, 3608S, 1:1,000), and anti-SORBS2 (SAB, 24643, 1:1,000). GAPDH was used as the internal control. Goat anti-mouse or goat anti-rabbit horseradish peroxidase-conjugated IgG was used as the secondary antibody (CST, 7076, 7074, 1:10,000). The membranes were incubated with secondary antibody for 2 h at room temperature, and bands were visualized using an enhanced chemiluminescence detection system (Amersham Bioscience, Piscataway, NJ, United States) according to the instructions of the manufacturer.
Transwell Assays
After digestion with trypsin, 5 × 104 cells were seeded in the top chamber containing a polycarbonate membrane with 8-μm pores (Corning, Cambridge, MA, United States). The lower chamber was filled with 10% FBS as a chemoattractant. After 48 h, the cells of the inner membrane were scraped using a cotton swab, and the cells on the lower surface were fixed with 4% paraformaldehyde and stained with 0.1% crystal violet for counting.
Quantitative Real-Time PCR and Percent Spliced-in Analysis
Total RNA from cancer cells and fresh tissues was extracted using TRIzol Reagent (Invitrogen) according to the protocols of the manufacturer. Then, the RNA molecules were dissolved in RNase-free DNase (TaKaRa, Dalian, China) and stored at −80°C. cDNA synthesis was performed using HiScript II Q Select RT SuperMix for qPCR (Vazyme, Nanjing, China). Finally, quantitative real-time PCR (qRT-PCR) was conducted using an Applied Biosystems 7500 Sequence Detection System (Applied Biosystems, Foster City, CA, United States) with AceQ® Universal SYBR qPCR Master Mix (Vazyme, Nanjing, China). The sequences of the sense and antisense primers are listed in Supplementary Table 4. GAPDH was used as an internal control. Relative gene expression was calculated using the 2–ΔΔCT. According to the PSI analysis in the TCGASpliceSeq database, two paired primers were designed for quantitative analysis. One pair of primers was used to amplify the alternative exon (V1) and the other for the common exon (V2). The relative PSI value of every AS event in the CRC tissues was calculated as V1/V2.
Coexpression Network Analysis
For these genes and AS events that are associated with poor or good prognosis, we introduced an algorithm called CoExpNetViz (version 2.0.2) () in the Cytoscape software (Version 3.6.0, citation) to create a coexpression network. First, all the correlation between these genes were calculated, then the lower and upper percentile of the distribution of coexpression values are taken, and these genes were considered as coexpressed. Finally, these coexpressed that passed the cutoff were used to constructed a coexpression network.
Gene Ontology Analysis
After we calculated and selected these differentially expressed genes (DEGs), we imported these DEGs into a database called Metascape4 () for Gene Ontology enrichment analysis. Then we selected the “Express Analysis” button for default analysis mode. Terms with a p-value < 0.01, a minimum count of 3, and an enrichment factor >1.5 are collected and grouped into clusters based on their membership similarities. More specifically, p-values are calculated based on the accumulative hypergeometric distribution 2, and q-values are calculated using the Banjamini–Hochberg procedure to account for multiple testings.
Hub Gene Detection
After the gene coexpression network was created, these gene pairs were then imported into CytoHubba, an algorithm used for identifying hub genes from complex interactomes (). We ranked these genes with MCC score calculated from CytoHubba and selected the top 5 ones as the most important hub genes for downstream validation.
Statistics Analysis
To compare the differentially expressed AS events between normal and tumor tissue, ANOVA analysis was performed using “stats” package in R. The differentially expressed genes between normal and tumor tissue were obtained using the “limma” package in R. Differentially expressed genes and AS events are displayed in a heatmap. Univariate Cox regression analysis was performed to obtain prognostic gene expression and gene AS events using the “survival” package in R and further displayed by forest plot. Other statistical analyses were performed using GraphPad Prism 6.0 (Inc., La Jolla, CA, United States). Comparisons between two groups were performed using a two-tailed unpaired Student’s t-test. Data are displayed as mean ± SD. Survival analysis was performed using the Kaplan–Meier method and compared using the log-rank test. Pearson’s correlation analysis was conducted to explore the correlation between the two variables. A value of p < 0.05 was considered statistically significant.
Results
Gene Expression Analysis and Gene Alternative Splicing Event Analysis Reveal Different Transcriptome Characteristics Involved in Colorectal Cancer Development
By comparing differentially expressed genes between 50 paired normal and tumor tissues, 2,313 genes were obtained and displayed by heatmap (Log2FC >1, adjusted p-value < 0.01, Figure 1A and Supplementary Table 4). Through ANOVA analysis, we found 338 AS events of 319 genes that were differentially expressed between normal and tumor tissues (difference >0.2, FDR <0.01, Figure 1B and Supplementary Table 5). Using the Venn diagram, we found that most of the genes with differentially expressed mRNA levels were different from the genes with differentially expressed AS events (Figure 1C). GO analysis revealed that pathways enriched by genes with differentially expressed mRNA levels (Figure 1D) were also inconsistent with pathways enriched by genes with differentially expressed AS events (Figure 1E). These results indicate that the carcinogenic mechanisms involved in gene expression differed from gene AS events.
FIGURE 1
Integrated Analysis Revealed Hub Gene Alternative Splicing Events in Colorectal Cancer
Of the 338 differentially expressed AS events, 12 upregulated AS events in tumor tissue were associated with poor prognosis, and four downregulated AS events in tumor tissue were associated with good prognosis (Figure 2A). However, of the 16 genes with prognostic AS events, the mRNA level of only one gene was related to poor prognosis (ACOXL, Figure 2B), and the mRNA level of only one gene correlated with good prognosis (EPB41L2, Figure 2B). By constructing a coexpression network of gene expression and gene AS events (R > 0.3, p < 0.001), we identified five hub gene AS events: EPB41L2, CELF2, TMEM130, VCL, and SORBS2 (Figure 2C and Supplementary Table 1). GO analysis revealed that the 300 genes whose mRNA levels were significantly related to the gene AS events were significantly involved in the GO terms (Figure 2D). We further examined the hub gene AS events and hub gene mRNA expression between five paired normal and tumor tissues by qRT-PCR, and found that the PSI of EPB41L2, TMEM130, and SORBS2 were different from tumor to normal (Figure 2E).
FIGURE 2
Identification of Splicing Factor Regulating the Hub Gene Alternative Splicing Events
Univariate Cox regression analysis revealed that three splicing factors were associated with good prognosis, and one splicing factor was related to poor prognosis (Figure 3A). By constructing a coexpression network of splicing factor gene expression and gene AS events (p < 0.05), four of the above five hub gene AS events were significantly related to NOVA1 mRNA expression (Figure 3B). Then we detected the expression of NOVA1 by qRT-PCR, and the results showed that the tumor was higher than normal (Figure 3C). We examined the expression of NOVA1 and selected CRC cell lines with a high expression (RKO) of NOVA1 for further experiments (Figure 3D).
FIGURE 3
NOVA1-Mediated SORBS2 Alter Splicing Promotes Colorectal Cancer Cell Migration via Notch Signaling Pathway
The above results indicated that the alternative splice of SORBS2-exon3 was downregulated in CRC tissue, which was defined as SORBS2-exon3. To explore the function of SORBS2-exon3 in CRC, after knocking down NOVA1 in CRC cells (Figure 4A), we also detected decreased AS events in TMEM130 and increased SORBS2 (Figure 4B). Of the four hub gene AS events, we found that SORBS2 was most apparently changed by NOVA1 knockdown. To verify the relationship between NOVA1, SORBS2-exon3, and NOVA1, RKO cells were transfected with siNOVA1 and/or siSORBS2-exon3. Western blotting showed that after knocking down NOVA1, SORBS2 containing exon3 increased, and the expression of NOTCH1 was suppressed. After knocking down SORBS2 containing exon3, the expression of NOTCH1 was increased. It is also surprising that the expression of NOVA1 was increased, indicating a negative feedback between NOVA1 and SORBS2-exon3, which may be the reason for CRC progression (Figure 4C). The qPCR results were consistent with the results of Western blotting (Figure 4D). Transwell experiments were preformed to detect the migration of CRC, and the result showed that NOVA1 knockdown-mediated CRC cell migration inhibition could be reversed by SORBS2-exon3 downregulation (Figure 4E), which indicated that NOVA1/SORBS2-exon3/NOTCH1 was involved in the migration of CRC (Figure 4F).
FIGURE 4
Discussion
Colorectal cancer is a major threat to the lives of individuals worldwide, and the mechanisms involved in oncogenesis and cancer development are very complex (). Gene alternative splicing makes the mechanisms more implicated because different isoforms of a gene may function in diverse ways. The inconsistency of different gene isoforms reminds us that we need to consider gene isoforms when studying genes. Using bioinformatics analysis, we first identified five hub gene AS events involved in CRC progression. We also verified that the gene AS events SORBS2 were upregulated in normal tissue, and gene AS events EPB41L2, CELF2, TMEM130, and VCL were downregulated in normal tissue by qRT-PCR. NOVA1 was identified as a critical splicing factor that regulates SORBS2 expression in CRC. Finally, in vitro experiments showed that NOVA1/SORBS2-exon3 promoted CRC cell migration through the Notch pathway. Therefore, a therapeutic regimen targeting SORBS2 may be a promising strategy to overcome CRC progression.
A previous analysis was performed to reveal the critical AS events associated with CRC development, such as AT in CXCL12, RI in CSTF3, AP in LBH, and ES in ALDH4A1 (; ). However, these studies did not simultaneously consider gene expression. Integrated analysis of gene expression and AS events can provide additional insights into the mechanisms of cancer progression (; ). By comparing differentially expressed genes and gene AS events, we found that the genes whose mRNA levels were significantly changed were different from those whose AS events were significantly changed. GO analysis revealed that the enriched GO terms using the differentially expressed genes also differed from the enriched GO terms using the genes with differentially expressed AS events. Coexpression network analysis of gene AS events and gene expression revealed five hub gene AS events involved in CRC progression: ES in EPB41L2, AP in CELF2, AP in TMEM130, ES in VCL, and AP in SORBS2. Overexpression of EPB41L2 in ovarian cancer has been reported to be associated with chemotherapy resistance (). CELF2 suppression promotes gastric cancer cell proliferation and migration (). Loss of vinculin can promote CRC metastasis and predict poor prognosis in patients with CRC (). In ovarian cancer, SORBS2 suppresses metastatic colonization of cancer cells by eliciting a tumor-suppressive immune microenvironment ().
GO analysis indicated that the hub gene AS events and NOVA1 were associated with Notch pathway activation. The splice factor and splice factor kinase antiapoptotic isoform CLK1 can be induced by hypoxia to produce the antiapoptotic isoform caspase 9b in prostate cancer cells, which may be involved in antitumor therapy resistance (). Our results showed that NOVA1 was negatively correlated with SORBS2-exon3, and NOVA1 was positively correlated with NOTCH1 expression. After knocking down NOVA1, the migration cell number of RKO cells was less than that of NC. After knocking down SORBS2-exon3, the migration cell number of RKO cells was greater than that of NC. Furthermore, NOVA1 knockdown-mediated CRC cell migration inhibition could be reversed by SORBS2-exon3 downregulation. As we know, including mutations in KRAS, NRAS, PI3K, and BRAF genes are considered to be the key factors in the transformation of CRC (; ). In addition, increasing number of scholars are paying more attention to the impact of epigenetics, including altered splicing, on tumor progression (; ; ; ). Our data showed that after knocking down SORBS2-exon3, even if NOVA1 was knocked down together, the loss of the CRC migration function mediated by NOVA1 would be rescued. Remarkably, SORBS2-exon3 will reverse the expression of NOVA1. More attention should be paid to the altered splicing of SORBS2 in clinical treatment, which may be helpful for the treatment of CRC metastasis. We will further verify this in future experiments.
In conclusion, we revealed that gene expression-mediated CRC development mechanisms differ from the gene AS events induced by CRC progression. Further screening and validation confirmed that NOVA1 mediated SORBS2 alternative splicing induced CRC cell migration via the Notch pathway. Therefore, our results reveal a promising therapeutic strategy targeting SORBS2-exon3 in CRC.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Statements
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Ethics statement
The studies involving human participants were reviewed and approved by Ethics Committee of Shanghai Ruijin Hospital affiliated to Shanghai Jiao Tong University School of Medicine. The patients/participants provided their written informed consent to participate in this study.
Author contributions
RZ and KL were responsible for the conceptualizing and designing of the study. TZ, SC, and YP analyzed and interpreted the patient data and were major contributors in the writing of the manuscript. TZ, CW, XC, and KL collected and assembled the data. YP and SC performed the experiments. All authors read and approved the final manuscript.
Funding
This study was supported by the National Nature Science Foundation of China (NSFC: 82002475), Shanghai Hospital Development Center (Number: 16CR2064B), Shanghai Anti-Cancer Association (Number: SCCRC17002), and Youth Development Program of Ruijin Hospital and Shanghai Jiao Tong University School of Medicine (2019QNPY01024).
Acknowledgments
We would like to thank Editage (www.editage.cn) for English language editing.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2021.673873/full#supplementary-material
Footnotes
1.^https://tcga-data.nci.nih.gov/tcga/
2.^http://projects.insilico.us.com/TCGASpliceSeq/index.jsp
References
1
AmirkhahR.Naderi-MeshkinH.ShahJ. S.DunneP. D.SchmitzU. (2019). The intricate interplay between epigenetic events, alternative splicing and noncoding RNA deregulation in colorectal cancer.Cells8:929. 10.3390/cells8080929
2
BisogninA.PizziniS.PerilliL.EspositoG.MocellinS.NittiD.et al (2014). An integrative framework identifies alternative splicing events in colorectal cancer development.Mol. Oncol.8129–141. 10.1016/j.molonc.2013.10.004
3
BowlerE.PorazinskiS.UzorS.ThibaultP.DurandM.LapointeE.et al (2018). Hypoxia leads to significant changes in alternative splicing and elevated expression of CLK splice factor kinases in PC3 prostate cancer cells.BMC Cancer18:355. 10.1186/s12885-018-4227-7
4
BrayF.FerlayJ.SoerjomataramI.SiegelR. L.TorreL. A.JemalA. (2018). Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries.CA Cancer J. Clin.68394–424. 10.3322/caac.21492
5
ChenY. S.LiuC. W.LinY. C.TsaiC. Y.YangC. H.LinJ. C. (2020). The SRSF3-MBNL1-Acin1 circuit constitutes an emerging axis to lessen DNA fragmentation in colorectal cancer via an alternative splicing mechanism.Neoplasia22702–713. 10.1016/j.neo.2020.10.002
6
ChinC.-H.ChenS.-H.WuH.-H.HoC.-W.KoM.-T.LinC.-Y. (2014). cytoHubba: identifying hub objects and sub-networks from complex interactome.BMC Syst. Biol.8:S11. 10.1186/1752-0509-8-S4-S11
7
Climente-GonzalezH.Porta-PardoE.GodzikA.EyrasE. (2017). The functional impact of alternative splicing in cancer.Cell Rep.202215–2226. 10.1016/j.celrep.2017.08.012
8
CuiQ.XingJ.GuY.NanX.MaW.ChenY.et al (2019). GXYLT2 accelerates cell growth and migration by regulating the Notch pathway in human cancer cells.Exp. Cell Res.3761–10. 10.1016/j.yexcr.2019.01.023
9
FengH.JinZ.LiuK.PengY.JiangS.WangC.et al (2020a). Identification and validation of critical alternative splicing events and splicing factors in gastric cancer progression.J. Cell. Mol. Med.2412667–12680. 10.1111/jcmm.15835
10
FengH.LiuK.ShenX.LiangJ.WangC.QiuW.et al (2020b). Correction: targeting tumor cell-derived CCL2 as a strategy to overcome Bevacizumab resistance in ETV5(+) colorectal cancer.Cell Death Discov.11:1006. 10.1038/s41419-020-03208-z
11
GabutM.Samavarchi-TehraniP.WangX.SlobodeniucV.O’HanlonD.SungH. K.et al (2011). An alternative splicing switch regulates embryonic stem cell pluripotency and reprogramming.Cell147132–146. 10.1016/j.cell.2011.08.023
12
HampelH.PearlmanR.BeightolM.ZhaoW.JonesD.FrankelW. L.et al (2018). Assessment of tumor sequencing as a replacement for lynch syndrome screening and current molecular tests for patients with colorectal cancer.JAMA Oncol.4806–813. 10.1001/jamaoncol.2018.0104
13
HeS.TongX.LuK.LuY.LuoJ.YangW.et al (2016). Comparative analysis of transcriptomes among bombyx mori strains and sexes reveals the genes regulating melanic morph and the related phenotypes.PLoS One11:e0155061. 10.1371/journal.pone.0155061
14
KalariK. R.RossellD.NecelaB. M.AsmannY. W.NairA.BahetiS.et al (2012). Deep sequence analysis of non-small cell lung cancer: integrated analysis of gene expression, alternative splicing, and single nucleotide variations in lung adenocarcinomas with and without oncogenic KRAS mutations.Front. Oncol.2:12. 10.3389/fonc.2012.00012
15
KentS.HutchinsonJ.BalboniA.DecastroA.CherukuriP.DirenzoJ. (2011). DeltaNp63alpha promotes cellular quiescence via induction and activation of Notch3.Cell Cycle103111–3118. 10.4161/cc.10.18.17300
16
KindlerH. L.ShulmanK. L. (2001). Metastatic colorectal cancer.Curr. Treat. Options Oncol.2459–471. 10.1007/s11864-001-0068-7
17
KochL. (2017). Alternative splicing: a thermometer controlling gene expression.Nat. Rev. Genet.18:515. 10.1038/nrg.2017.61
18
LaissueP. (2019). The forkhead-box family of transcription factors: key molecular players in colorectal cancer pathogenesis.Mol. Cancer18:5. 10.1186/s12943-019-0938-x
19
LiW.WangD.SunX.ZhangY.WangL.SuoJ. (2019). ADAM17 promotes lymph node metastasis in gastric cancer via activation of the Notch and Wnt signaling pathways.Int. J. Mol. Med.43914–926. 10.3892/ijmm.2018.4028
20
LiY.SunN.LuZ.SunS.HuangJ.ChenZ.et al (2017). Prognostic alternative mRNA splicing signature in non-small cell lung cancer.Cancer Lett.39340–51. 10.1016/j.canlet.2017.02.016
21
LinP.HeR. Q.MaF. C.LiangL.HeY.YangH.et al (2018). Systematic analysis of survival-associated alternative splicing signatures in gastrointestinal pan-adenocarcinomas.EBioMedicine3446–60. 10.1016/j.ebiom.2018.07.040
22
LinR. Z.Moreno-LunaR.LiD.JaminetS. C.GreeneA. K.Melero-MartinJ. M. (2014). Human endothelial colony-forming cells serve as trophic mediators for mesenchymal stem cell engraftment via paracrine signaling.Proc. Natl. Acad. Sci. U. S. A.11110137–10142. 10.1073/pnas.1405388111
23
LiuJ.LiH.ShenS.SunL.YuanY.XingC. (2018). Alternative splicing events implicated in carcinogenesis and prognosis of colorectal cancer.J. Cancer91754–1764. 10.7150/jca.24569
24
LudlowA. T.WongM. S.RobinJ. D.BattenK.YuanL.LaiT. P.et al (2018). NOVA1 regulates hTERT splicing and cell growth in non-small cell lung cancer.Nat. Commun.9:3112. 10.1038/s41467-018-05582-x
25
MaoS.LiY.LuZ.CheY.SunS.HuangJ.et al (2019). Survival-associated alternative splicing signatures in esophageal carcinoma.Carcinogenesis40121–130. 10.1093/carcin/bgy123
26
MenyhartO.FeketeJ. T.GyorffyB. (2019). Gene expression indicates altered immune modulation and signaling pathway activation in ovarian cancer patients resistant to topotecan.Int. J. Mol. Sci.20:2750. 10.3390/ijms20112750
27
MojicaW.HawthornL. (2010). Normal colon epithelium: a dataset for the analysis of gene expression and alternative splicing events in colon disease.BMC Genomics11:5. 10.1186/1471-2164-11-5
28
MuY.ZouH.ChenB.FanY.LuoS. (2017). FAM83D knockdown regulates proliferation, migration and invasion of colorectal cancer through inhibiting FBXW7/Notch-1 signalling pathway.Biomed. Pharmacother.90548–554. 10.1016/j.biopha.2017.03.073
29
PamarthyS.MaoL.KataraG. K.FleetwoodS.KulshreshtaA.Gilman-SachsA.et al (2016). The V-ATPase a2 isoform controls mammary gland development through Notch and TGF-beta signaling.Cell Death Dis.7:e2443. 10.1038/cddis.2016.347
30
PannutiA.ForemanK.RizzoP.OsipoC.GoldeT.OsborneB.et al (2010). Targeting notch to target cancer stem cells.Clin. Cancer Res.163141–3152. 10.1158/1078-0432.CCR-09-2823
31
PimentelH.ParraM.GeeS.GhanemD.AnX.LiJ.et al (2014). A dynamic alternative splicing program regulates gene expression during terminal erythropoiesis.Nucleic Acids Res.424031–4042. 10.1093/nar/gkt1388
32
SongJ.LiuY. D.SuJ.YuanD.SunF.ZhuJ. (2019). Systematic analysis of alternative splicing signature unveils prognostic predictor for kidney renal clear cell carcinoma.J. Cell. Physiol.23422753–22764. 10.1002/jcp.28840
33
SveenA.AgesenT. H.NesbakkenA.RognumT. O.LotheR. A.SkotheimR. I. (2011). Transcriptome instability in colorectal cancer identified by exon microarray analyses: associations with splicing factor expression levels and patient survival.Genome Med.3:32. 10.1186/gm248
34
TakebeN.MieleL.HarrisP. J.JeongW.BandoH.KahnM.et al (2015). Targeting Notch, Hedgehog, and Wnt pathways in cancer stem cells: clinical update.Nat. Rev. Clin. Oncol.12445–464. 10.1038/nrclinonc.2015.61
35
TangX. J.WangW.HannS. S. (2019). Interactions among lncRNAs, miRNAs and mRNA in colorectal cancer.Biochimie16358–72. 10.1016/j.biochi.2019.05.010
36
TzfadiaO.DielsT.De MeyerS.VandepoeleK.AharoniA.Van de PeerY. (2015). CoExpNetViz: comparative co-expression networks construction and visualization tool.Front. Plant Sci.6:1194. 10.3389/fpls.2015.01194
37
VakianiE.JanakiramanM.ShenR.SinhaR.ZengZ.ShiaJ.et al (2012). Comparative genomic analysis of primary versus metastatic colorectal carcinomas.J. Clin. Oncol.302956–2962. 10.1200/JCO.2011.38.2994
38
WanL.YuW.ShenE.SunW.LiuY.KongJ.et al (2019). SRSF6-regulated alternative splicing that promotes tumour progression offers a therapy target for colorectal cancer.Gut68118–129. 10.1136/gutjnl-2017-314983
39
WangJ.LiuL.SunY.XueY.QuJ.PanS.et al (2018). miR-615-3p promotes proliferation and migration and inhibits apoptosis through its potential target CELF2 in gastric cancer.Biomed. Pharmacother.101406–413. 10.1016/j.biopha.2018.02.104
40
WangX.LiuJ.WangD.FengM.WuX. (2020). Epigenetically regulated gene expression profiles reveal four molecular subtypes with prognostic and therapeutic implications in colorectal cancer.Brief. Bioinform.22:bbaa309. 10.1093/bib/bbaa309
41
XiongY.DengY.WangK.ZhouH.ZhengX.SiL.et al (2018). Profiles of alternative splicing in colorectal cancer and their clinical significance: a study based on large-scale sequencing data.EBioMedicine36183–195. 10.1016/j.ebiom.2018.09.021
42
YueY.ZhouK.LiJ.JiangS.LiC.MenH. (2018). MSX1 induces G0/G1 arrest and apoptosis by suppressing Notch signaling and is frequently methylated in cervical cancer.Onco. Targets Ther.114769–4780. 10.2147/OTT.S165144
43
ZhangT. H.LiangL. Z.LiuX. L.WuJ. N.SuK.ChenJ. Y.et al (2019). LncRNA UCA1/miR-124 axis modulates TGFbeta1-induced epithelial-mesenchymal transition and invasion of tongue cancer cells through JAG1/Notch signaling.J. Cell. Biochem.12010495–10504. 10.1002/jcb.28334
44
ZhaoL.WangW.HuangS.YangZ.XuL.YangQ.et al (2018). The RNA binding protein SORBS2 suppresses metastatic colonization of ovarian cancer by stabilizing tumor-suppressive immunomodulatory transcripts.Genome Biol.19:35. 10.1186/s13059-018-1412-6
45
ZhouY.ZhouB.PacheL.ChangM.KhodabakhshiA. H.TanaseichukO.et al (2019). Metascape provides a biologist-oriented resource for the analysis of systems-level datasets.Nat. Commun.10:1523. 10.1038/s41467-019-09234-6
46
ZongF. Y.FuX.WeiW. J.LuoY. G.HeinerM.CaoL. J.et al (2014). The RNA-binding protein QKI suppresses cancer-associated aberrant splicing.PLoS Genet.10:e1004289. 10.1371/journal.pgen.1004289
Summary
Keywords
colorectal cancer, alternative splice, SORBS2-exon3, NOVA1, migration
Citation
Zhang T, Chen S, Peng Y, Wang C, Cheng X, Zhao R and Liu K (2021) NOVA1-Mediated SORBS2 Isoform Promotes Colorectal Cancer Migration by Activating the Notch Pathway. Front. Cell Dev. Biol. 9:673873. doi: 10.3389/fcell.2021.673873
Received
28 February 2021
Accepted
08 September 2021
Published
08 October 2021
Volume
9 - 2021
Edited by
Lucio Miele, Louisiana State University, United States
Reviewed by
Mary C. Farach-Carson, University of Texas Health Science Center at Houston, United States; Lilian Sales Gomez, Mayo Clinic, United States
Updates
Copyright
© 2021 Zhang, Chen, Peng, Wang, Cheng, Zhao and Liu.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Xi Cheng, drchengxi@126.comRen Zhao, rjzhaoren@139.comKun Liu, lookiere@126.com
†These authors have contributed equally to this work
This article was submitted to Cell Adhesion and Migration, a section of the journal Frontiers in Cell and Developmental Biology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.