An integrative framework identifies alternative splicing events in colorectal cancer development

Long-read sequencing reveals the landscape of aberrant alternative splicing and novel therapeutic target in colorectal cancer

BACKGROUND

Alternative splicing complexity plays a vital role in carcinogenesis and cancer progression. Improved understanding of novel splicing events and the underlying regulatory mechanisms may contribute new insights into developing new therapeutic strategies for colorectal cancer (CRC).

METHODS

Here, we combined long-read sequencing technology with short-read RNA-seq methods to investigate the transcriptome complexity in CRC. By using experiment assays, we explored the function of newly identified splicing isoform TIMP1 Δ4-5. Moreover, a CRISPR/dCasRx-based strategy to induce the TIMP1 exon 4-5 exclusion was introduced to inhibit neoplasm growth.

RESULTS

A total of 90,703 transcripts were identified, of which > 62% were novel compared with current transcriptome annotations. These novel transcripts were more likely to be sample specific, expressed at relatively lower levels with more exons, and oncogenes displayed a characteristic to generate more transcripts in CRC. Clinical outcome data analysis showed that 1472 differentially expressed alternative splicing events (DEAS) were tightly associated with CRC patients' prognosis, and many novel isoforms were likely to be important determinants for patient survival. Among these, newly identified splicing isoform TIMP1 Δ4-5 was significantly downregulated in CRC. Further in vitro and in vivo assays demonstrated that ectopic expression of TIMP1 Δ4-5 significantly suppresses tumor cell growth and metastasis. Serine/arginine-rich splicing factor 1 (SRSF1) acts as a onco-splicing regulator through sustaining the inclusion of TIMP1 exon 4-5. Furthermore, CRISPR/dCasRx-based strategies designed to induce TIMP1 exon 4-5 exclusion have the potential to restrain the CRC growth.

CONCLUSIONS

This data provides a rich resource for deeper studies of gastrointestinal malignancies. Newly identified splicing isoform TIMP1 Δ4-5 plays an important role in mediating CRC progression and may be a potential therapy target in CRC.

Comprehensive Analysis of Alternative Polyadenylation Events Associated with the Tumor Immune Microenvironment in Colon Adenocarcinoma

Objective

Colon adenocarcinoma (COAD) is one of the leading causes of cancer death worldwide. Alternative polyadenylation (APA) is relevant to the variability of the 3'-UTR of mRNA. However, the posttranscriptional dysregulation of APA in COAD is poorly understood.

Methods

We collected APA data from The Cancer Genome Atlas (TCGA) COAD (n =7692). APA events were evaluated using PDUI values, and the prognostically significant APA events were screened by LASSO Cox regression to construct a prognostic model. Then, prognostic model functions and possible regulatory genes of characteristic APA events were analyzed. Finally, the immune regulatory network based on APA regulatory genes was analyzed and established.

Results

A total of 95 APA events were found to influence the COAD outcomes. Among them, 39 genes were screened as characteristic prognostic APA events by LASSO Cox regression to construct a COAD prognostic signature. The analysis results suggested that a high signature score was associated with poor prognosis and was significantly correlated with a variety of immune cells, including NK and Th1, 2 and 17 cells. Further analysis showed that APA regulators mainly served roles in the prognosis of COAD. Based on the above results, we constructed an immunoregulatory network for APA regulatory genes-APA genes-immune cells.

Conclusion

Our study revealed that APA events in COAD may regulate tumor progression by influencing immune cells, which provides a new direction for exploring the influencing mechanism of the tumor immune microenvironment and is expected to provide a potential new target for COAD immunotherapy.

Alternative splicing patterns reveal prognostic indicator in muscle-invasive bladder cancer

Background

Bladder cancer is one of the most lethal malignancy in urological system, and 20–25% of bladder cancer patients are muscle invasive with unfavorable prognosis. However, the role of alternative splicing (AS) in muscle-invasive bladder cancer (MIBC) remains to be elucidated.

Methods

Percent spliced in (PSI) data obtained from the Cancer Genome Atlas (TCGA) SpliceSeq database (n = 394) were utilized to evaluate the AS events in MIBC. Prognosis-associated AS events were screened out by univariate Cox regression. LASSO Cox regression was used to identify reliable prognostic patterns in a training set and further validated in a test set. Splicing regulatory networks were constructed by correlations between PSI of AS events and RNA expression of splicing factors.

Results

As a result, a total of 2589 prognosis-related AS events in MIBC were identified. Pathways of spliceosomal complex (FDR = 0.017), DNA-directed RNA polymerase II, core complex (FDR = 0.032), and base excision repair (FDR = 0.038) were observed to be significantly enriched. Additionally, we noticed that most of the prognosis-related AS events were favorable factors. According to the LASSO and multivariate Cox regression analyses, 15-AS-based signature was established with the area under curve (AUC) of 0.709, 0.823, and 0.857 at 1-, 3-, and 5- years, respectively. The MIBC patients were further divided into high- and low-risk groups based on median risk sores. Interestingly, we observed that the prevalence of FGFR3 with mutations and focal amplification was significantly higher in low-risk group. Functional and immune infiltration analysis suggested potential signaling pathways and distinct immune states between these two groups. Moreover, splicing correlation network displayed a regulatory mode of prognostic splicing factors (SF) in MIBC patients.

Conclusions

This study not only provided novel insights into deciphering the possible mechanism of tumorgenesis and pathogenesis but also help refine risk stratification systems and potential treatment of decision-making for MIBC.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12957-022-02685-0.

Landscape of Alternative Splicing Events Related to Prognosis and Immune Infiltration in Glioma: A Data Analysis and Basic Verification

Background

Glioma is a prevalent primary brain cancer with high invasiveness and typical local diffuse infiltration. Alternative splicing (AS), as a pervasive transcriptional regulatory mechanism, amplifies the coding capacity of the genome and promotes the progression of malignancies. This study was aimed at identifying AS events and novel biomarkers associated with survival for glioma.

Methods

RNA splicing patterns were collected from The Cancer Genome Atlas SpliceSeq database, followed by calculating the percentage of splicing index. Expression profiles and related clinical information of glioma were integrated based on the UCSC Xena database. The AS events in glioma were further analyzed, and glioma prognosis-related splicing factors were identified with the use of bioinformatics analysis and laboratory techniques. Further immune infiltration analysis was performed.

Results

Altogether, 9028 AS events were discovered. Upon univariate Cox analysis, 425 AS events were found to be related to the survival of patients with glioma, and 42 AS events were further screened to construct the final prognostic model (area under the curve = 0.974). Additionally, decreased expression of the splicing factors including Neuro-Oncological Ventral Antigen 1 (NOVA1), heterogeneous nuclear ribonucleoprotein C (HNRNPC), heterogeneous nuclear ribonucleoprotein L-like protein (HNRNPLL), and RNA-Binding Motif Protein 4 (RBM4) contributed to the poor survival in glioma. The immune infiltration analysis demonstrated that AS events were related to the proportion of immune cells infiltrating in glioma.

Conclusions

It is of great value for comprehensive consideration of AS events, splicing networks, and related molecular subtype clusters in revealing the underlying mechanism and immune microenvironment remodeling for glioma, which provides clues for the further verification of related therapeutic targets.

Multimerin-1 and cancer: a review

Multimerin-1 (MMRN1) is a platelet protein with a role in haemostasis and coagulation. It is also present in endothelial cells (ECs) and the extracellular matrix (ECM), where it may be involved in cell adhesion, but its molecular functions and protein–protein interactions in these cellular locations have not been studied in detail yet. In recent years, MMRN1 has been identified as a differentially expressed gene (DEG) in various cancers and it has been proposed as a possible cancer biomarker. Some evidence suggest that MMRN1 expression is regulated by methylation, protein interactions, and non-coding RNAs (ncRNAs) in different cancers. This raises the questions if a functional role of MMRN1 is being targeted during cancer development, and if MMRN1’s differential expression pattern correlates with cancer progression. As a result, it is timely to review the current state of what is known about MMRN1 to help inform future research into MMRN1’s molecular mechanisms in cancer.

The pattern of alternative splicing in lung adenocarcinoma shows novel events correlated with tumorigenesis and immune microenvironment

Lung adenocarcinoma (LUAD) is the leading cause of cancer deaths worldwide due to the lack of early diagnostic markers and specific drugs. Previous studies have shown the association of LUAD growth with aberrant alternative splicing (AS). Herein, clinical data of 535 tumor tissues and 59 normal tissues were extracted from The Cancer Genome Atlas (TCGA) database. Each sample was analyzed using the ESTIMATE algorithm; a comparison between higher and lower score groups (stromal or immune) was made to determine the overall- and progression-free survival-related differentially expressed AS (DEAS) events. We then performed unsupervised clustering of these DEASs, followed by determining their relationship with survival rate, immune cells, and the tumor microenvironment (TME). Next, two prognostic signatures were developed using bioinformatics tools to explore the prognosis of cases with LUAD. Five OS- and six PFS-associated DEAS events were implemented to establish a prognostic risk score model. When compared to the high-risk group (HRG), the PFS and OS of the low-risk group (LRG) were found to be considerable. Additionally, a better prognosis was found considerably associated with the ESTIMATE score of the patients as well as immune cells infiltration. Our analysis of AS events in LUAD not only helps to clarify the tumorigenesis mechanism of AS but also provides ideas for revealing potential prognostic biomarkers and therapeutic targets.

CALD1 promotes the expression of PD-L1 in bladder cancer via the JAK/STAT signaling pathway

Background

Bladder cancer (BC) is a common malignant neoplasm with a high rate of recurrence and progression, despite optimal treatment. There is a pressing need to identify new effective biomarkers for the targeted treatment of BC.

Methods

The key gene CALD1 was screened via weighed gene co-expression network analysis (WGCNA) from encoding protein genes of BC. Clinical and prognostic significance was explored in The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. Cell Counting Kit-8 (CCK-8), flow cytometry, transwell chamber experiment and nude mouse xenograft assay were performed to test cell growth, apoptosis, migration, invasion and tumorigenesis capacities. Immune correlation was analyzed in The Tumor Immune Estimation Resource (TIMER) database. Relevant signaling pathways were explored using gene set enrichment analysis (GSEA).

Results

Increased expression of CALD1 was significantly correlated with histological grade, clinical stage, T stage, and lymphatic metastasis. Kaplan-Meier survival curves showed that high CALD1 expression was associated with poor overall survival (OS) and disease-free survival (DFS) in TCGA database, and with poor OS in the four GEO databases. CALD1 promotes growth, migration, invasion, and cell cycle of tumor cell, and inhibits tumor cell apoptosis in vitro and in vivo. CADL1 expression was positively correlated with increased CD274 levels (r=0.357, P=9.71e−14). JAK/STAT signaling pathway was significantly enriched in the high CALD1 expression group. CALD1-mediated PD-L1 overexpression (OE) was via the activation of the JAK/STAT signaling pathway; this effect was blocked by the specific JAK inhibitor Ruxolitinib.

Conclusions

CALD1 is a potential molecular marker associated with prognosis. It promotes the malignant progression of BC and upregulates the PD-L1 expression via the JAK/STAT signaling pathway.

NOVA1-Mediated SORBS2 Isoform Promotes Colorectal Cancer Migration by Activating the Notch Pathway

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.

Identification of the role of mono-ADP-ribosylation in colorectal cancer by integrated transcriptome analysis

Our previous study clarified the carcinogenic properties of arginine-specific mono-ADP-ribosyltransferase 1 (ART1), which results in a critical post-translational modification that changes the structure and function of proteins and is widely involved in important processes. This study provides, for the first time, a comprehensive transcriptomic analysis of colorectal cancer cells with ART1 silencing by Illumina RNA-Seq and related verification experiments. Lentiviral infection was used to construct a CT-26 cell line with stable knockdown of the ART1 gene, and a whole transcriptome sequencing technique was performed to identify differentially expressed genes (DEGs). GO and KEGG classification/enrichment analyses and verification experiments were performed to determine the role of ART1 in the progression of colorectal cancer. A total of 5552 DEGs, GO functions and KEGG pathways with the highest enrichment, various SNPs, and diverse splicing patterns were identified. Importantly, knockdown of ART1 affected the splicing of certain key genes related to tumor cell growth and downregulated the expression of the key gene PTBP1 for alternative splicing. The overall attenuation of the endoplasmic reticulum unfolded protein response (UPR) signaling pathway caused by the inhibition of mono-ADP-ribosylation of GRP78 could disrupt UPR signaling, leading to the occurrence of apoptosis to impede tumorigenesis. ART1, which is clustered in organelles, may promote the development of colorectal cancer by participating in a variety of new mechanisms, including endoplasmic reticulum stress regulation and alternative splicing, and may be a good clinical drug target for targeted therapy of CRC.

Alternative splicing of mRNA in colorectal cancer: new strategies for tumor diagnosis and treatment

Alternative splicing (AS) is an important event that contributes to posttranscriptional gene regulation. This process leads to several mature transcript variants with diverse physiological functions. Indeed, disruption of various aspects of this multistep process, such as cis- or trans- factor alteration, promotes the progression of colorectal cancer. Therefore, targeting some specific processes of AS may be an effective therapeutic strategy for treating cancer. Here, we provide an overview of the AS events related to colorectal cancer based on research done in the past 5 years. We focus on the mechanisms and functions of variant products of AS that are relevant to malignant hallmarks, with an emphasis on variants with clinical significance. In addition, novel strategies for exploiting the therapeutic value of AS events are discussed.

Prognostic Alternative Splicing Signatures in Esophageal Carcinoma

Alternative splicing (AS) is a method of increasing the number of proteins that the genome is capable of coding for, by altering the pre-mRNA during its maturation. This process provides the ability of a broad range of proteins to arise from a single gene. AS events are known to occur in up to 94% of human genes. Cumulative data have shown that aberrant AS functionality is a major factor in human diseases. This review focuses on the contribution made by aberrant AS functionality in the development and progression of esophageal cancer. The changes in the pattern of expression of alternately spliced isoforms in esophageal cancer can be used as diagnostic or prognostic biomarkers. Additionally, these can be used as targets for the development of new treatments for esophageal cancer.

The roles and prognostic significance of ABI1-TSV-11 expression in patients with left-sided colorectal cancer

Abnormally expressed and/or phosphorylated Abelson interactor 1 (ABI1) participates in the metastasis and progression of colorectal cancer (CRC). ABI1 presents as at least 12 transcript variants (TSVs) by mRNA alternative splicing, but it is unknown which of them is involved in CRC metastasis and prognosis. Here, we firstly identified ABI1-TSV-11 as a key TSV affecting the metastasis and prognosis of left-sided colorectal cancer (LsCC) and its elevated expression is related to lymph node metastasis and shorter overall survival (OS) in LsCC by analyzing data from The Cancer Genome Atlas and TSVdb. Secondly, ABI1-TSV-11 overexpression promoted LoVo and SW480 cells adhesion and migration in vitro, and accelerated LoVo and SW480 cells lung metastasis in vivo. Finally, mechanism investigations revealed that ABI1-isoform-11 interacted with epidermal growth factor receptor pathway substrate 8 (ESP8) and regulated actin dynamics to affect LoVo and SW480 cells biological behaviors. Taken together, our data demonstrated that ABI1-TSV-11 plays an oncogenic role in LsCC, it is an independent risk factor of prognosis and may be a potential molecular marker and therapeutic target in LsCC.

Alternative splicing: An important regulatory mechanism in colorectal carcinoma

Alternative splicing (AS) is a process that produces various mRNA splicing isoforms via different splicing patterns of mRNA precursors (pre-mRNAs). AS is the primary mechanism for increasing the types and quantities of proteins to improve biodiversity and influence multiple biological processes, including chromatin modification, signal transduction, and protein expression. It has been reported that AS is involved in the tumorigenesis and development of colorectal carcinoma (CRC). In this review, we delineate the concept, types, regulatory processes, and technical advances of AS and focus on the role of AS in CRC initiation, progression, treatment, and prognosis. This summary of the current knowledge about AS will contribute to our understanding of CRC initiation and development. This study will help in the discovery of novel biomarkers and therapeutic targets for CRC prognosis and treatment.

Alterations of RNA splicing patterns in esophagus squamous cell carcinoma

Alternative splicing (AS) is an important biological process for regulating the expression of various isoforms from a single gene and thus to promote proteome diversity. In this study, RNA-seq data from 15 pairs of matched esophageal squamous cell carcinoma (ESCC) and normal tissue samples as well as two cell lines were analyzed. AS events with significant differences were identified between ESCC and matched normal tissues, which were re-annotated to find protein coding genes or non-coding RNAs. A total of 45,439 AS events were found. Of these, 6019 (13.25%) significant differentially AS events were identified. Exon skipping (SE) events occupied the largest proportion of abnormal splicing events. Fifteen differential splicing events with the same trends of ΔΨ values in ESCC tissues, as well in the two cell lines were found. Four pathways and 20 biological processes related to pro-metastasis cell junction and migration were significantly enriched for the differentially spliced genes. The upregulated splicing factor SF3B4, which regulates 92 gene splicing events, could be a potential prognostic factor of ESCC. Differentially spliced genes, including HNRNPC, VCL, ZNF207, KIAA1217, TPM1 and CALD1 are shown with a sashimi plot. These results suggest that cell junction- and migration-related biological processes are influenced by AS abnormalities, and aberrant splicing events can be affected by splicing factor expression changes. The involved splicing factor SF3B4 was found to be a survival-related gene in ESCC and is presumed to regulate AS in multiple cancers. In summary, we identified significant differentially expressed AS events which may be related to the development of ESCC.

Identification and validation of critical alternative splicing events and splicing factors in gastric cancer progression

Gene expression and alternative splicing (AS) interact in complex ways to regulate biological process which is associated with cancer development. Here, by integrated analysis of gene expression and AS events, we aimed to identify the hub AS events and splicing factors relevant in gastric cancer development (GC). RNA‐seq data, clinical data and AS events of 348 GC 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 53 GCs. By bioinformatics methods, we found that gene AS event‐ and gene expression‐mediated GC progression shared the same mechanisms, such as PI3K/AKT pathway, but the involved genes were different. Though expression of 17 hub AS events were confirmed in 53 GC tissues, only 10 AS events in seven genes were identified as critical candidates related to GC progression, notably the AS events (Exon Skip) in CLSTN1 and SEC16A. Expression of these AS events in GC correlated with activation of the PI3K/AKT pathway. Genes with AS events associated with clinical parameters and prognosis were different from the genes whose mRNA levels were related to clinical parameters and prognosis. Besides, we further revealed that QKI and NOVA1 were the crucial splicing factors regulating expression of AS events in GC, but not spliceosome genes. Our integrated analysis revealed hub AS events in GC development, which might be the potential therapeutic targets for GC.

Cervical cancer in low and middle-income countries

Cervical cancer is a malignant tumour that occurs in the cervix and is classified into two histological types, adenocarcinoma and squamous cell carcinoma (SCC); SCC is more common and accounts for 70% of all cases. In 2018 there were ~569,000 new cases of cervical cancer diagnosed worldwide and ~311,000 deaths were attributed to cervical cancer. Of these, between 84 and 90% occurred in low- and middle-income countries (LMICs) such as South Africa, India, China and Brazil. The most common cause of cervical cancer is persistent infection caused by the sexually transmitted human papilloma virus. Other factors that contribute to the incidence of cervical cancer include geography, traditional practices and beliefs, the screening levels, socioeconomic status, healthcare access, public awareness, use of oral contraceptives, smoking and co-infection with HIV. An estimated 11 million women from LMICs will be diagnosed with cervical cancer in the next 10-20 years. The aim of this review was to explore various types of genetic and epigenetic factors that influence the development, progression or suppression of cervical cancer.

Post-transcriptional regulation activity through alternative splicing involved in the effects of Aloe vera on the Wnt/β-catenin and Notch pathways in colorectal cancer cells

Aloe vera (L.) Burm.f. is widely used as laxative drugs, cosmetics, and functional food due to a variety of therapeutic effects. However, several studies indicated a colonic carcinogenic activity of Aloe vera. But the underline mechanism has not been well clarified. This study aimed to explore the potential mechanism at the post-transcriptional level. Identification of Differential Expressed Alternative Splicing (DEAS) genes and events and the corresponding functional enrichment analyses were conducted on RKO cells after treated with Aloe vera aqueous extract and its two active components, aloin and aloesin. And RT-qPCR was conducted for validation. Results indicated that they induced 2200, 2342 and 2133 DEAS events, respectively. The GO enrichment and the COG classification results of DEAS genes showed that they were associated with transcription, as well as functions like signal transduction mechanisms. Moreover, DEAS genes related to the two colorectal cancerous pathways, Wnt and Notch pathways, were annotated. In conclusion, aloe extract, aloin and aloesin significantly regulated the DEAS profile of RKO cells. The colonic carcinogenicity of Aloe vera may due to its post-transcriptional regulatory activity through Alternative Splicing (AS) on genes, especially on Wnt-related and Notch-related key genes.

Comprehensive characterization of the alternative splicing landscape in head and neck squamous cell carcinoma reveals novel events associated with tumorigenesis and the immune microenvironment

Alternative splicing (AS) has emerged as a key event in tumor development and microenvironment formation. However, comprehensive analysis of AS and its clinical significance in head and neck squamous cell carcinoma (HNSC) is urgently required. Methods: Genome-wide profiling of AS events using RNA-Seq data from The Cancer Genome Atlas (TCGA) program was performed in a cohort of 464 patients with HNSC. Cancer-associated AS events (CASEs) were identified between paired HNSC and adjacent normal tissues and evaluated in functional enrichment analysis. Splicing networks and prognostic models were constructed using bioinformatics tools. Unsupervised clustering of the CASEs identified was conducted and associations with clinical, molecular and immune features were analyzed. Results: We detected a total of 32,309 AS events and identified 473 CASEs in HNSC; among these, 91 were validated in an independent cohort (n = 15). Functional protein domains were frequently altered, especially by CASEs affecting cancer drivers, such as PCSK5. CASE parent genes were significantly enriched in pathways related to HNSC and the tumor immune microenvironment, such as the viral carcinogenesis (FDR < 0.001), Human Papillomavirus infection (FDR < 0.001), chemokine (FDR < 0.001) and T cell receptor (FDR < 0.001) signaling pathways. CASEs enriched in immune-related pathways were closely associated with immune cell infiltration and cytolytic activity. AS regulatory networks suggested a significant association between splicing factor (SF) expression and CASEs and might be regulated by SF methylation. Eighteen CASEs were identified as independent prognostic factors for overall and disease-free survival. Unsupervised clustering analysis revealed distinct correlations between AS-based clusters and prognosis, molecular characteristics and immune features. Immunogenic features and immune subgroups cooperatively depict the immune features of AS-based clusters. Conclusion: This comprehensive genome-wide analysis of the AS landscape in HNSC revealed novel AS events related to carcinogenesis and immune microenvironment, with implications for prognosis and therapeutic responses.

Oncogenic splicing abnormalities induced by DEAD-Box Helicase 56 amplification in colorectal cancer

Alternative splicing, regulated by DEAD‐Box Helicase (DDX) families, plays an important role in cancer. However, the relationship between the DDX family and cancer has not been fully elucidated. In the present study, we identified a candidate oncogene DDX56 on Ch.7p by a bioinformatics approach using The Cancer Genome Atlas (TCGA) dataset of colorectal cancer (CRC). DDX56 expression was measured by RT‐qPCR and immunochemical staining in 108 CRC patients. Clinicopathological and survival analyses were carried out using three CRC datasets. Biological roles of DDX56 were explored by gene set enrichment analysis (GSEA), and cell proliferation in vitro and in vivo, cell cycle assays, and using DDX56‐knockdown or overexpressed CRC cells. RNA sequencing was carried out to elucidate the effect of DDX56 on mRNA splicing. We found that DDX56 expression was positively correlated with the amplification of DDX56 and was upregulated in CRC cells. High DDX56 expression was associated with lymphatic invasion and distant metastasis and was an independent poor prognostic factor. In vitro analysis, in vivo analysis and GSEA showed that DDX56 promoted proliferation ability through regulating the cell cycle. DDX56 knockdown reduced intron retention and tumor suppressor WEE1 expression, which functions as a G2‐M DNA damage checkpoint. We have identified DDX56 as a novel oncogene and prognostic biomarker of CRC that promotes alternative splicing of WEE1.

Determining the prognostic significance of alternative splicing events in soft tissue sarcoma using data from The Cancer Genome Atlas

BACKGROUND

Surgery, adjuvant chemotherapy, and radiotherapy are the primary treatment options for soft tissue sarcomas (STSs). However, identifying ways to improve the prognosis of patients with STS remains a considerable challenge. Evidence shows that the dysregulation of alternative splicing (AS) events is involved in tumor pathogenesis and progression. The present study objective was to identify survival-associated AS events that could serve as prognostic biomarkers and potentially serve as tumor-selective STS drug targets.

METHODS

STS-specific 'percent spliced in' (PSI) values for splicing events in 206 STS samples were downloaded from The Cancer Genome Atlas SpliceSeq^® database. Prognostic analyses were performed on seven types of AS events to determine their prognostic value in STS patients, for which prediction models were constructed with the risk score formula [Formula: see text]. Prediction models were also constructed to determine the prognostic value of AS events, and Spearman's rank correlation coefficients were calculated to determine the degree of correlation between splicing factor expression and the PSI values.

RESULTS

A total 10,439 events were found to significantly correlate with patient survival rates. The area under the time-dependent receiver operating characteristic curve for the prognostic predictor of STS overall survival was 0.826. Notably, the splicing events of certain STS key genes were significantly associated with STS 2-year overall survival in the present study, including exon skip (ES) events in MDM2 and EWSR1, alternate terminator events in CDKN2A and HMGA2 for dedifferentiated liposarcoma, ES in MDM2 and alternate promoter events in CDKN2A for leiomyosarcoma, and ES in EWSR1 for undifferentiated pleomorphic sarcoma. Moreover, splicing correlation networks between AS events and splicing factors revealed that almost all of the AS events showed negatively correlations with the expression of splicing factors.

CONCLUSION

An in-depth analysis of alternative RNA splicing could provide new insights into the mechanisms of STS oncogenesis and the potential for novel approaches to this type of cancer therapy.

The Intricate Interplay between Epigenetic Events, Alternative Splicing and Noncoding RNA Deregulation in Colorectal Cancer

Colorectal cancer (CRC) results from a transformation of colonic epithelial cells into adenocarcinoma cells due to genetic and epigenetic instabilities, alongside remodelling of the surrounding stromal tumour microenvironment. Epithelial-specific epigenetic variations escorting this process include chromatin remodelling, histone modifications and aberrant DNA methylation, which influence gene expression, alternative splicing and function of non-coding RNA. In this review, we first highlight epigenetic modulators, modifiers and mediators in CRC, then we elaborate on causes and consequences of epigenetic alterations in CRC pathogenesis alongside an appraisal of the complex feedback mechanisms realized through alternative splicing and non-coding RNA regulation. An emphasis in our review is put on how this intricate network of epigenetic and post-transcriptional gene regulation evolves during the initiation, progression and metastasis formation in CRC.

Genome-Wide Profiling of Prognostic Alternative Splicing Signature in Colorectal Cancer

Background: This study was to explore differential RNA splicing patterns and elucidate the function of the splice variants served as prognostic biomarkers in colorectal cancer (CRC).

Methods: Genome-wide profiling of prognostic alternative splicing (AS) events using RNA-seq data from The Cancer Genome Atlas (TCGA) program was conducted to evaluate the roles of seven AS patterns in 330 colorectal cancer cohort. The prognostic predictors models were assessed by integrated Cox proportional hazards regression. Based on the correlations between survival associated AS events and splicing factors, splicing networks were built.

Results: A total of 2,158 survival associated AS events in CRC were identified. Interestingly, most of these top 20 survival associated AS events were adverse prognostic factors. The prognostic models were built by each type of splicing patterns, performing well for risk stratification in CRC patients. The area under curve (AUC) of receiver operating characteristic (ROC) for the combined prognostic predictors model could reach 0.963. Splicing network also suggested distinguished correlation between the expression of splicing factors and AS events in CRC patients.

Conclusion: The ideal prognostic predictors model for risk stratification in CRC patients was constructed by differential splicing patterns of 13 genes. Our findings enriched knowledge about differential RNA splicing patterns and the regulation of splicing, providing generous biomarker candidates and potential targets for the treatment of CRC.

Profiles of alternative splicing in colorectal cancer and their clinical significance: A study based on large-scale sequencing data

Background

Alternative splicing (AS), as a potent and pervasive mechanism of transcriptional regulatory, expands the genome's coding capacity and involves in the initiation and progression of cancer. Systematic analysis of alternative splicing in colorectal cancer (CRC) is lacking and greatly needed.

Methods

RNA-Seq data and corresponding clinical information of CRC cohort were downloaded from the TCGA data portal. Then, a java application, known as SpliceSeq, was used to evaluate the RNA splicing patterns and calculate the Percent Spliced In (PSI) value. Differently expressed AS events (DEAS) were identified based on PSI value between paired CRC and adjacent tissues. DEAS and its splicing networks were further analyzed by bioinformatics methods. Kaplan-Meier, Cox proportional regression and unsupervised clustering analysis were used to evaluate the association between DEAS and patients' clinical features.

Results

After strict filtering, a total of 34,334 AS events were identified, among which 421 AS events were found expressed differently. Parent genes of these DEAS play a important role in regulating CRC-related processes such as protein kinase activity (FDR<0.0001), PI3K-Akt signaling pathway (FDR = 0.0024) and p53 signaling pathway (FDR = 0.0143). 37 DEAS events were found to be associated with OS, and 68 DEAS events were found to be associated with DFS. Stratifying patients according to the PSI value of AT in CXCL12 and RI in CSTF3 formed significant Kaplan-Meier curves in both OS and DFS survival analysis. Unsupervised clustering analysis using DEAS revealed four clusters with distinct survival patterns, and associated with consensus molecular subtypes.

Conclusions

Large differences of AS events in CRC appear to exist, and these differences are likely to be important determinants of both prognosis and biological regulation. Our identified CRC-related AS events and uncovered splicing networks are valuable in deciphering the underlying mechanisms of AS in CRC, and provide clues of therapeutic targets to further validations.

A Systems-Level Analysis Reveals Circadian Regulation of Splicing in Colorectal Cancer

Accumulating evidence points to a significant role of the circadian clock in the regulation of splicing in various organisms, including mammals. Both dysregulated circadian rhythms and aberrant pre-mRNA splicing are frequently implicated in human disease, in particular in cancer. To investigate the role of the circadian clock in the regulation of splicing in a cancer progression context at the systems-level, we conducted a genome-wide analysis and compared the rhythmic transcriptional profiles of colon carcinoma cell lines SW480 and SW620, derived from primary and metastatic sites of the same patient, respectively. We identified spliceosome components and splicing factors with cell-specific circadian expression patterns including SRSF1, HNRNPLL, ESRP1, and RBM 8A, as well as altered alternative splicing events and circadian alternative splicing patterns of output genes (e.g., VEGFA, NCAM1, FGFR2, CD44) in our cellular model. Our data reveals a remarkable interplay between the circadian clock and pre-mRNA splicing with putative consequences in tumor progression and metastasis.

The Usage of Exon-Exon Splice Junctions for the Detection of Alternative Splicing using the REIDS model

Alternative gene splicing is a common phenomenon in which a single gene gives rise to multiple transcript isoforms. The process is strictly guided and involves a multitude of proteins and regulatory complexes. Unfortunately, aberrant splicing events have been linked to genetic disorders. Therefore, understanding mechanisms of alternative splicing regulation and differences in splicing events between diseased and healthy tissues is crucial in advancing personalized medicine and drug developments. We propose a linear mixed model, Random Effects for the Identification of Differential Splicing (REIDS), for the identification of alternative splicing events using Human Transcriptome Arrays (HTA). For each exon, a splicing score is calculated based on two scores, an exon score and an array score. The junction information is used to rank the identified exons from strongly confident to less confident candidates for alternative splicing. The design of junctions was also discussed to highlight the complexity of exon-exon and exon-junction interactions. Based on a list of Rt-PCR validated probe sets, REIDS outperforms AltAnalyze and iGems in the % recall rate.

RNA analysis of cancer predisposing genes in formalin-fixed paraffin-embedded tissue determines aberrant splicing

High-throughput sequencing efforts in molecular tumour diagnostics detect increasing numbers of novel variants, including variants predicted to affect splicing. In silico prediction tools can reliably predict the effect of variant disrupting canonical splice sites; however, experimental validation is required to confirm aberrant splicing. Here, we present RNA analysis performed for 13 canonical splice site variants predicted or known to result in splicing in the cancer predisposition genes MLH1, MSH2, MSH6, APC and BRCA1. Total nucleic acid was successfully isolated for 10 variants from eight formalin-fixed paraffin-embedded (FFPE) tumour tissues and two B-cell lines. Aberrant splicing was confirmed in all six variants known to result in splicing. Of one known variant in the B-cell line, aberrant splicing could only be detected after formalin fixation, which indicated that formalin fixation could possibly inhibit RNA degradation. Aberrant splicing was concluded in three of four predicted splice variants of uncertain significance, supporting their pathogenic effect. With this assay, somatic splice variants can be easily and rapidly analysed, enabling retrospective analysis to support the pathogenicity of variants predicted to result in splicing when only FFPE material is available.

Proteome Heterogeneity in Colorectal Cancer

Tumor heterogeneity is an important feature of colorectal cancer (CRC) manifested by dynamic changes in gene expression, protein expression, and availability of different tumor subtypes. Recent publications in the past 10 years have revealed proteome heterogeneity between different colorectal tumors and within the same tumor site. This paper reviews recent research works on the proteome heterogeneity in CRC, which includes the heterogeneity within a single tumor (intratumor heterogeneity), between different anatomical sites at the same organ, and between primary and metastatic sites (intertumor heterogeneity). The potential use of proteome heterogeneity in precision medicine and its implications in biomarker discovery and therapeutic outcomes will be discussed. Identification of the unique proteome landscape between and within individual tumors is imperative for understanding cancer biology and the management of CRC patients.

Integrated analysis of differential expression and alternative splicing of non-small cell lung cancer based on RNA sequencing

Non-small cell lung cancer (NSCLC) is the most common type of lung cancer, with high morbidity and mortality rates. Numerous diagnosis and treatment methods have been proposed, and the prognosis of NSCLC has improved to a certain extent. However, the mechanisms of NSCLC remain largely unknown, and additional studies are required. In the present study, the RNA sequencing dataset of NSCLC was downloaded from the Gene Expression Omnibus (http://www.ncbi.nlm.nih.gov/geo/). The clean reads obtained from the raw data were mapped to the University of California Santa Cruz human genome (hg19), based on TopHat, and were assembled into transcripts via Cufflink. The differential expression (DE) and differential alternative splicing (DAS) genes were screened out through Cuffdiff and rMATS, respectively. The significantly enriched gene ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes pathways were obtained through the Database of Annotation, Visualization and Integrated Discovery (DAVID). Different numbers of DE and DAS genes were identified in different types of NSCLC samples, but a number of common functions and pathways were obtained, including biological processes associated with abnormal immune and cell activity. GO terms and pathways associated with substance metabolism, including the insulin signaling pathway and oxidative phosphorylation, were enriched in DAS genes rather than DE genes. Integrated analysis of differential expression and alternative splicing may be helpful in understanding the mechanisms of NSCLC, in addition to its early diagnosis and treatment.

A random effects model for the identification of differential splicing (REIDS) using exon and HTA arrays

Background

Alternative gene splicing is a common phenomenon in which a single gene gives rise to multiple transcript isoforms. The process is strictly guided and involves a multitude of proteins and regulatory complexes. Unfortunately, aberrant splicing events do occur which have been linked to genetic disorders, such as several types of cancer and neurodegenerative diseases (Fan et al., Theor Biol Med Model 3:19, 2006). Therefore, understanding the mechanism of alternative splicing and identifying the difference in splicing events between diseased and healthy tissue is crucial in biomedical research with the potential of applications in personalized medicine as well as in drug development.

Results

We propose a linear mixed model, Random Effects for the Identification of Differential Splicing (REIDS), for the identification of alternative splicing events. Based on a set of scores, an exon score and an array score, a decision regarding alternative splicing can be made. The model enables the ability to distinguish a differential expressed gene from a differential spliced exon. The proposed model was applied to three case studies concerning both exon and HTA arrays.

Conclusion

The REIDS model provides a work flow for the identification of alternative splicing events relying on the established linear mixed model. The model can be applied to different types of arrays.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-017-1687-8) contains supplementary material, which is available to authorized users.

A saga of cancer epigenetics: linking epigenetics to alternative splicing

The discovery of an increasing number of alternative splicing events in the human genome highlighted that ∼94% of genes generate alternatively spliced transcripts that may produce different protein isoforms with diverse functions. It is now well known that several diseases are a direct and indirect consequence of aberrant splicing events in humans. In addition to the conventional mode of alternative splicing regulation by 'cis' RNA-binding sites and 'trans' RNA-binding proteins, recent literature provides enormous evidence for epigenetic regulation of alternative splicing. The epigenetic modifications may regulate alternative splicing by either influencing the transcription elongation rate of RNA polymerase II or by recruiting a specific splicing regulator via different chromatin adaptors. The epigenetic alterations and aberrant alternative splicing are known to be associated with various diseases individually, but this review discusses/highlights the latest literature on the role of epigenetic alterations in the regulation of alternative splicing and thereby cancer progression. This review also points out the need for further studies to understand the interplay between epigenetic modifications and aberrant alternative splicing in cancer progression.

CircRNAs in hematopoiesis and hematological malignancies

Cell states in hematopoiesis are controlled by master regulators and by complex circuits of a growing family of RNA species impacting cell phenotype maintenance and plasticity. Circular RNAs (circRNAs) are rapidly gaining the status of particularly stable transcriptome members with distinctive qualities. RNA-seq identified thousands of circRNAs with developmental stage- and tissue-specific expression corroborating earlier suggestions that circular isoforms are a natural feature of the cell expression program. CircRNAs are abundantly expressed also in the hematopoietic compartment. There are a number of studies on circRNAs in blood cells, a specific overview is however lacking. In this review we first present current insight in circRNA biogenesis discussing the relevance for hematopoiesis of the highly interleaved processes of splicing and circRNA biogenesis. Regarding molecular functions circRNAs modulate host gene expression, but also compete for binding of microRNAs, RNA-binding proteins or translation initiation and participate in regulatory circuits. We examine circRNA expression in the hematopoietic compartment and in hematologic malignancies and review the recent breakthrough study that identified pathogenic circRNAs derived from leukemia fusion genes. CircRNA high and regulated expression in blood cell types indicate that further studies are warranted to inform the position of these regulators in normal and malignant hematopoiesis.

Genome-wide profiling of chemoradiation‑induced changes in alternative splicing in colon cancer cells

Alternative splicing is a key mechanism that regulates protein diversity and has been found to be associated with colon cancer progression and metastasis. However, the function of alternative splicing in chemoradiation‑resistant colon cancer remains elusive. In this study, we constructed a chemoradiation‑resistant colon cancer cell line. Through RNA-sequencing of normal and chemoradiation‑resistant colon cancer cells (HCT116), we found 818 genes that were highly expressed in the normal HCT116 cells, whereas 285 genes were highly expressed in the chemoradiation-resistant HCT116 (RCR-HCT116) cells. Gene ontology (GO) analysis showed that genes that were highly expressed in the HCT116 cells were enriched in GO categories related to cell cycle and cell division, whereas genes that were highly expressed in the RCR-HCT116 cells were associated with regulation of system processes and response to wounding. Analysis of alternative splicing events revealed that exon skipping was significantly increased in the chemoradiation‑resistant colon cancer cells. Moreover, we identified 323 alternative splicing events in 293 genes that were significantly different between the two different HCT116 cell types. These alternative splicing‑related genes were clustered functionally into several groups related with DNA replication, such as deoxyribonucleotide metabolic/catabolic processes, response to DNA damage stimulus and helicase activity. These findings enriched our knowledge by elucidating the function of alternative splicing in chemoradiation-resistant colon cancer.

Endoluminal high-resolution MR imaging protocol for colon walls analysis in a mouse model of colitis

OBJECTIVE

An endoluminal magnetic resonance (MR) imaging protocol including the design of an endoluminal coil (EC) was defined for high-spatial-resolution MR imaging of mice gastrointestinal walls at 4.7 T.

MATERIALS AND METHODS

A receive-only radiofrequency single-loop coil was developed for mice colon wall imaging. Combined with a specific protocol, the prototype was first characterized in vitro on phantoms and on vegetables. Signal-to-noise ratio (SNR) profiles were compared with a quadrature volume birdcage coil (QVBC). Endoluminal MR imaging protocol combined with the EC was assessed in vivo on mice.

RESULTS

The SNR measured close to the coil is significantly higher (10 times and up to 3 mm of the EC center) than the SNR measured with the QVBC. The gain in SNR can be used to reduce the in-plane pixel size up to 39 × 39 µm(2) (234 µm slice thickness) without time penalty. The different colon wall layers can only be distinguished on images acquired with the EC.

CONCLUSION

Dedicated EC provides suitable images for the assessment of mice colon wall layers. This proof of concept provides gains in spatial resolution and leads to adequate protocols for the assessment of human colorectal cancer, and can now be used as a new imaging tool for a better understanding of the pathology.

The Associations between RNA Splicing Complex Gene SF3A1 Polymorphisms and Colorectal Cancer Risk in a Chinese Population

BACKGROUND

Aberrant alternative splicing included alterations in components of the mRNA splicing machinery often occurred in colon cancer. However, the role of SF3A1, one key component of the mRNA splicing machinery, on colorectal cancer (CRC) risk was still not elucidated.

METHOD AND FINDINGS

We performed a hospital-based case-control study containing 801 CRC patients and 817 cancer-free controls to examine the association between SF3A1 polymorphisms and CRC risk in a Chinese population. Four candidate SNPs (rs10376, rs5753073, rs2839998 and rs2074733) were selected based on bioinformatics analysis and previous findings. The results showed no significant associations between these SNPs and CRC risk (P > 0.05). Besides, the stratified analysis based on the smoking and alcohol use status obtained no statistically significant results.

CONCLUSION

Our study was the first one to investigate the association between SF3A1 polymorphisms and CRC risk. The results suggested these four SNPs in SF3A1 were not associated with CRC risk in a Chinese population, however, further more studies are needed to confirm our findings.

Integrated exon level expression analysis of driver genes explain their role in colorectal cancer

Integrated analysis of genomic and transcriptomic level changes holds promise for a better understanding of colorectal cancer (CRC) biology. There is a pertinent need to explain the functional effect of genome level changes by integrating the information at the transcript level. Using high resolution cytogenetics array, we had earlier identified driver genes by ‘Genomic Identification of Significant Targets In Cancer (GISTIC)’ analysis of paired tumour-normal samples from colorectal cancer patients. In this study, we analyze these driver genes at three levels using exon array data – gene, exon and network. Gene level analysis revealed a small subset to experience differential expression. These results were reinforced by carrying out separate differential expression analyses (SAM and LIMMA). ATP8B1 was found to be the novel gene associated with CRC that shows changes at cytogenetic, gene and exon levels. Splice index of 29 exons corresponding to 13 genes was found to be significantly altered in tumour samples. Driver genes were used to construct regulatory networks for tumour and normal groups. There were rearrangements in transcription factor genes suggesting the presence of regulatory switching. The regulatory pattern of AHR gene was found to have the most significant alteration. Our results integrate data with focus on driver genes resulting in highly enriched novel molecules that need further studies to establish their role in CRC.

BCLAF1 and its splicing regulator SRSF10 regulate the tumorigenic potential of colon cancer cells

Bcl-2-associated transcription factor 1 (BCLAF1) is known to be involved in multiple biological processes. Although several splice variants of BCLAF1 have been identified, little is known about how BCLAF1 splicing is regulated or the contribution of alternative splicing to its developmental functions. Here we find that inclusion of alternative exon5a was significantly increased in colorectal cancer (CRC) samples. Knockdown of the BCLAF1 protein isoform resulting from exon5a inclusion inhibited growth and that its overexpression increased tumorigenic potential. We also found that the splicing factor SRSF10 stimulates inclusion of exon5a and has growth-inducing activity. Importantly, the upregulation of SRSF10 expression observed in clinical CRC samples parallels the increased inclusion of BCLAF1 exon5a, both of which are associated with higher tumour grade. These findings identify SRSF10 as a key regulator of BCLAF1 pre-mRNA splicing and the maintenance of oncogenic features in human colon cancer cells.

An integrative framework identifies alternative splicing events in colorectal cancer development

Alternative splicing (AS) is a common mechanism which creates diverse RNA isoforms from a single gene, potentially increasing protein variety. Growing evidence suggests that this mechanism is closely related to cancer progression. In this study, whole transcriptome analysis was performed with GeneChip Human exon 1.0 ST Array from 80 samples comprising 23 normal colon mucosa, 30 primary colorectal cancer and 27 liver metastatic specimens from 46 patients, to identify AS events in colorectal cancer progression. Differentially expressed genes and exons were estimated and AS events were reconstructed by combining exon-level analyses with AltAnalyze algorithms and transcript-level estimations (MMBGX probabilistic method). The number of AS genes in the transition from normal colon mucosa to primary tumor was the most abundant, but fell considerably in the next transition to liver metastasis. 206 genes with probable AS events in colon cancer development and progression were identified, that are involved in processes and pathways relevant to tumor biology, as cell-cell and cell-matrix interactions. Several AS events in VCL, CALD1, B3GNT6 and CTHRC1 genes, differentially expressed during tumor development were validated, at RNA and at protein level. Taken together, these results demonstrate that cancer-specific AS is common in early phases of colorectal cancer natural history.