OncoSplicing: an updated database for clinically relevant alternative splicing in 33 human cancers

The Macrophage migration inhibitory factor is a vital player in Pan-Cancer by functioning as a M0 Macrophage biomarker

BACKGROUND

The role of the macrophage migration inhibitory factor (MIF) has recently attracted considerable attention in cancer research; nonetheless, the insights provided by current investigations remain constrained. Our main objective was to investigate its role and the latent mechanisms within the pan-cancer realm.

METHODS

We used comprehensive pan-cancer bulk sequencing data and online network tools to investigate the association between MIF expression and patient prognosis, genomic instability, cancer cell stemness, DNA damage repair, and immune infiltration. Furthermore, we validated the relationship between MIF expression and M0 macrophages using single-cell datasets, the SpatialDB database, and fluorescence staining. Additionally, we assessed the therapeutic response using the ROC plotter tool.

RESULTS

We observed the upregulation of MIF expression across numerous cancer types. Notably, elevated MIF levels were associated with a decline in genomic stability. We found a significant correlation between increased MIF expression and increased expression of mismatch repair genes, stemness features, and homologous recombination genes across diverse malignancies. Subsequently, through an analysis using ESTIMATE and cytokine results, we revealed the involvement of MIF in immune suppression. Then, we validated MIF as a hallmark of the M0 macrophages involved in tumor immunity. Our study suggests an association with other immune-inhibitory cellular populations and restraint of CD8 + T cells. In addition, we conducted a comparative analysis of MIF expression before and after treatment in three distinct sets of therapy responders and non-responders. Intriguingly, we identified notable disparities in MIF expression patterns in bladder urothelial carcinoma and ovarian cancer following particular therapeutic interventions.

CONCLUSION

Comprehensive pan-cancer analysis revealed notable enrichment of MIF within M0 macrophages, exerting a profound influence on tumor-associated immunosuppression and the intricate machinery of DNA repair.

Exploring Splicing Modulation as an Innovative Approach to Combat Pancreatic Cancer: SF3B1 Emerges as a Prognostic Indicator and Therapeutic Target

Pancreatic ductal adenocarcinoma (PDAC) poses significant challenges in terms of prognosis and treatment. Recent research has identified splicing deregulation as a new cancer hallmark. Herein, we investigated the largely uncharacterized alternative splicing profile and the key splicing factor SF3B1 in PDAC pancreatic cells and tissues as a potential discovery source of plausible drug targets and new predictive biomarkers of clinical outcome. The research involved a transcriptome-wide analysis, comparing profiles of splicing profiles in PDAC primary cells with normal ductal cells. This revealed more than 400 significant differential splicing events in genes involved in regulation of gene expression, primarily related to mRNA splicing, and metabolism of nucleic acids. PDAC cultures were highly sensitive to the SF3B1 modulators, E7107 and Pladienolide-B, showing IC50s in the low nanomolar range. These compounds induced apoptosis, associated to induction of the MCL-1/S splice variant. and reduced cell migration, associated to RON mis-splicing. In an orthotopic mouse model, E7107 showed promising results. Furthermore, we evaluated SF3B1 expression in specimens from 87 patients and found a significant association of SF3B1 expression with progression-free and overall survival. In conclusion, SF3B1 emerges as both a potential prognostic factor and therapeutic target in PDAC, impacting cell proliferation, migration, and apoptosis. These findings warrant future studies on this new therapeutic strategy against PDAC.

Expression landscape of cancer-FOXP3 and its prognostic value in pancreatic adenocarcinoma

FOXP3, a key identifier of Treg, has also been identified in tumor cells, which is referred to as cancer-FOXP3 (c-FOXP3). Human c-FOXP3 undergoes multiple alternative splicing events, generating several isoforms, like c-FOXP3FL and c-FOXP3Δ3. Previous research on c-FOXP3 often ignore its cellular source (immune or tumor cells) and isoform expression patterns, which may obscure our understanding of its clinical significance. Our immunohistochemistry investigations which conducted across 18 tumors using validated c-FOXP3 antibodies revealed distinct expression landscapes for c-FOXP3 and its variants, with the majority of tumors exhibited a predominantly expression of c-FOXP3Δ3. In pancreatic ductal adenocarcinoma (PDAC), we further discovered a potential link between nuclear c-FOXP3Δ3 in tumor cells and poor prognosis. Overexpression of c-FOXP3Δ3 in tumor cells was associated with metastasis. This work elucidates the expression pattern of c-FOXP3 in pan-cancer and indicates its potential as a prognostic biomarker in clinical settings, offering new perspectives for its clinical application.

A 3'-pre-tRNA-derived small RNA tRF-1-Ser regulated by 25(OH)D promotes proliferation and stemness by inhibiting the function of MBNL1 in breast cancer

BACKGROUND

We explored the potential novel anticancer mechanisms of 25-hydroxyvitamin D (25(OH)D), a vitamin D metabolite with antitumour effects in breast cancer. It is stable in serum and is used to assess vitamin D levels in clinical practice. Transfer RNA-derived small RNAs are small noncoding RNAs that generate various distinct biological functions, but more research is needed on their role in breast cancer.

METHODS

Small RNA microarrays were used to explore the novel regulatory mechanism of 25(OH)D. High-throughput RNA-sequencing technology was used to detect transcriptome changes after 25(OH)D treatment and tRF-1-Ser knockdown. RNA pull-down and high-performance liquid chromatography-mass spectrometry/mass spectrometry were used to explore the proteins bound to tRF-1-Ser. In vitro and in vivo functional experiments were conducted to assess the influence of 25(OH)D and tRF-1-Ser on breast cancer. Semi-quantitative PCR was performed to detect alternative splicing events. Western blot assay and qPCR were used to assess protein and mRNA expression.

RESULTS

The expression of tRF-1-Ser is negatively regulated by 25(OH)D. In our breast cancer (BRCA) clinical samples, we found that the expression of tRF-1-Ser was higher in cancer tissues than in paired normal tissues, and was significantly associated with tumour invasion. Moreover, tRF-1-Ser inhibits the function of MBNL1 by hindering its nuclear translocation. Functional experiments and transcriptome data revealed that the downregulation of tRF-1-Ser plays a vital role in the anticancer effect of 25(OH)D.

CONCLUSIONS

In brief, our research revealed a novel anticancer mechanism of 25(OH)D, unveiled the vital function of tRF-1-Ser in BRCA progression, and suggested that tRF-1-Ser could emerge as a new therapeutic target for BRCA.

Autophagy-related CMTM6 promotes glioblastoma progression by activating Wnt/β-catenin pathway and acts as an onco-immunological biomarker

BACKGROUND

Glioblastoma multiforme (GBM) is identified as one of the most prevalent and malignant brain tumors, characterized by poor treatment outcomes and a limited prognosis. CMTM6, a membrane protein, has been found to upregulate the expression of programmed cell death 1 ligand 1 protein (PD-L1) and acts as an immune checkpoint inhibitor by inhibiting the programmed death 1 protein/PD-L1 signaling pathway. Recent research has demonstrated a high expression of CMTM6 in GBM, suggesting its potential role in influencing the pathogenesis and progression of GBM, as well as its association with immune cell infiltration in the tumor microenvironment. However, the underlying mechanism of CMTM6 in GBM requires further investigation.

METHODS

Data from cancer patients in The Cancer Genome Atlas, Gene Expression Omnibus and Chinese Glioma Genome Atlas cohorts were consolidated for the current study. Through multi-omics analysis, the study systematically examined the expression profile of CMTM6, epigenetic modifications, prognostic significance, biological functions, potential mechanisms of action and alterations in the immune microenvironment. Additionally, the study investigated CMTM6 expression in GBM cell lines and normal cells using reverse transcription PCR and western blot analysis. The impact of CMTM6 on GBM cell proliferation, migration and invasion was evaluated using a combination of cell counting kit-8 assay, clone formation assay, 5-ethynyl-2'-deoxyuridine incorporation assay, wound healing assay and Transwell assay. In order to explore the mechanism of CMTM6, the Wnt/β-catenin signaling pathway and autophagy-related genes were further verified through western blot analysis.

RESULTS

CMTM6 is highly expressed in multiple tumors, particularly GBM. CMTM6 has been shown to be a valuable diagnostic and prognostic biomarker by various bioinformatics approaches. Additionally, CMTM6 plays a pivotal role in the pathogenesis of cancer, specifically GBM, by modulating various biological processes such as DNA methyltransferase expression, RNA modification, copy number variation, genomic heterogeneity, tumor stemness and DNA methylation. The findings of the experiment indicate a significant correlation between elevated CMTM6 expression and the proliferation, invasion, migration and autophagy of GBM cells, with potential key roles mediated through the Wnt/β-catenin signaling pathway. Furthermore, CMTM6 is implicated in modulating tumor immune cell infiltration and is closely linked to the expression of various immune checkpoint inhibitors and immune modulators, particularly within the context of GBM. High levels of CMTM6 expression also enhance the responsiveness of GBM patients to radiotherapy and chemotherapy, thereby offering valuable insights for guiding treatment strategies for GBM.

CONCLUSIONS

Autophagy-related CMTM6 is highly expressed in various types of cancer, especially GBM, and it can regulate GBM progression through the Wnt/β-catenin signaling pathway and is capable of being used as an underlying target for the diagnosis, treatment selection and prognosis of patients with GBM.

CD74 as a prognostic and M1 macrophage infiltration marker in a comprehensive pan-cancer analysis

CD74 is a type-II transmembrane glycoprotein that has been linked to tumorigenesis. However, this association was based only on phenotypic studies, and, to date, no in-depth mechanistic studies have been conducted. In this study, combined with a multi-omics study, CD74 levels were significantly upregulated in most cancers relative to normal tissues and were found to be predictive of prognosis. Elevated CD74 expression was associated with reduced levels of mismatch-repair genes and homologous repair gene signatures in over 10 tumor types. Multiple fluorescence staining and bulk, spatial, single-cell transcriptional analyses indicated its potential as a marker for M1 macrophage infiltration in pan-cancer. In addition, CD74 expression was higher in BRCA patients responsive to conventional chemotherapy and was able to predict the prognosis of these patients. Potential CD74-activating drugs (HNHA and BRD-K55186349) were identified through molecular docking to CD74. The findings indicate activation of CD74 may have potential in tumor immunotherapy.

Multi-omics Analyses Reveal Function of Apolipoprotein E in Alternative Splicing and Tumor Immune Microenvironment in Kidney Renal Clear Cell Carcinoma via Pan-cancer Analysis

Apolipoprotein E (APOE) regulates lipid metabolism, associated with the development of various cancers. However, its precise prognostic significance and functions in alternative splicing and the tumor immune microenvironment remain unclear. In this study, we extracted APOE expression in pan-cancer from TCGA and analyzed mRNA transcriptome, cell lines, and protein levels. Furthermore, we analyzed the alternative splicing expression of the APOE gene transcript with prognostic profiles using the OncoSplicing database. We obtained 73 common APOE genes to perform functional enrichment analysis, assess the correlation between genes and immune cells using TIMER, EPIC, and ssGSEA methods, and examine the prognostic significance using the UALCAN database. Finally, single-cell data was employed to assess the correlation between APOE genes and cell functions. Our findings revealed that APOE expression varies across different tumor types and cancer cell lines. The alternative splicing analysis demonstrated that APOE transcript expression levels have prognostic value in cancers such as LGG, KIRC, and KIRP. Functional enrichment analysis indicated significant associations between APOE and various immune cells, such as macrophages, CD8 T cells, and NK cells, with significant implications for prognosis. Moreover, single-cell data indicated that APOE was primarily expressed in renal epithelial cells among stromal cells and in macrophages among immune cells, significantly negatively correlated with five functional states. Our study represents the first comprehensive exploration of APOE's function in pan-cancers and identifies APOE as a potential biomarker in cancer pathogenesis, prognosis, and immune therapeutic target.

Bioinformatics analysis of the prognostic role of alternative splicing data in lung adenocarcinoma

Background

As a post-transcriptional regulatory mechanism, alternative splicing (AS) is engaged in a variety of pathophysiological processes, and it has been widely reported in connection with the occurrence, progression, metastasis, and drug resistance of cancer. However, the research on AS in lung adenocarcinoma (LUAD) is very limited. In addition, the prognostic effect of AS event (ASE) on LUAD and its related mechanism are not clear. This study aimed to explore the role and potential prognostic value of ASE in LUAD.

Methods

Relevant data and ASE datasets of the sample were acquired from The Cancer Genome Atlas (TCGA) and TCGASpliceSeq databases. We constructed a new prognostic criterion based on ASEs. Then, Cox regression and least absolute shrinkage and selection operator (LASSO) regression analysis were used to construct the model. Based on this model, the risk score of each ASE was calculated, and the reliability of this model was evaluated by Kaplan-Meier survival and receiver operating characteristic (ROC) curve analyses. Finally, these results were verified on different network platforms.

Results

We identified seven types of ASEs related to survival. The prognostic risk model for ASEs was established. The Kaplan-Meier curve showed that compared to the low-risk group, the overall survival (OS) rate of LUAD patients in the high-risk group was lower. ROC curve analysis showed that the prognostic risk model of LUAD patients was well predicted, and the area under the curve (AUC) also confirmed this.

Conclusions

This study screened the ASE related to the prognosis of LUAD patients, and provided a theoretical basis for further study of the correlation between ASE and the prognosis of LUAD patients. It has provided new ideas for developing new biomarkers and therapeutic targets for LUAD patients.

PSME2 offers value as a biomarker of M1 macrophage infiltration in pan-cancer and inhibits osteosarcoma malignant phenotypes

A growing number of studies have revealed an association between proteasome activator complex subunit 2 (PSME2) and the progression of various forms of cancer. However, the effect of PSME2 on osteosarcoma progression is unknown. Pan-cancer analyses focused on the immunological activity and prognostic relevance of PSME2 have yet to be conducted. The Cancer Genome Atlas and Genome-Tissue Expression databases were leveraged to evaluate PSME2 expression and activity across 33 cancer types. Significant PSME2 dysregulation was noted in a wide range of cancer types and this gene was found to offer significant diagnostic and prognostic utility in most analyzed cancers. From a mechanistic perspective, PSME2 expression levels were correlated with DNA methylation, DNA repair, genomic instability, and TME scores in multiple cancer types. PSME2 was subsequently established as a pan-cancer biomarker of M1 macrophage infiltration based on a combination of bulk, single-cell, and spatial transcriptomic data and confirmatory fluorescent staining results. In osteosarcoma cells, overexpressing PSME2 significantly suppressed tumor proliferative, migratory, and invasive activity. Screening efforts also successfully identified the PSME2-activating drug irinotecan, which can synergistically promote the death of osteosarcoma cells when combined with the chemotherapeutic drug paclitaxel. As a biomarker of M1 macrophage infiltration, PSME2 expression levels may offer insight into tumor development and progression for a wide range of cancers including osteosarcoma, emphasizing its potential utility as a prognostic and therapeutic target worthy of further study.

MAJIQlopedia: an encyclopedia of RNA splicing variations in human tissues and cancer

Quantification of RNA splicing variations based on RNA-Sequencing can reveal tissue- and disease-specific splicing patterns. To study such splicing variations, we introduce MAJIQlopedia, an encyclopedia of splicing variations that encompasses 86 human tissues and 41 cancer datasets. MAJIQlopedia reports annotated and unannotated splicing events for a total of 486 175 alternative splice junctions in normal tissues and 338 317 alternative splice junctions in cancer. This database, available at https://majiq.biociphers.org/majiqlopedia/, includes a user-friendly interface that provides graphical representations of junction usage quantification for each junction across all tissue or cancer types. To demonstrate case usage of MAJIQlopedia, we review splicing variations in genes WT1, MAPT and BIN1, which all have known tissue or cancer-specific splicing variations. We also use MAJIQlopedia to highlight novel splicing variations in FDX1 and MEGF9 in normal tissues, and we uncover a novel exon inclusion event in RPS6KA6 that only occurs in two cancer types. Users can download the database, request the addition of data to the webtool, or install a MAJIQlopedia server to integrate proprietary data. MAJIQlopedia can serve as a reference database for researchers seeking to understand what splicing variations exist in genes of interest, and those looking to understand tissue- or cancer-specific splice isoform usage.

Integrated Multi-omics Analyses Identify CDCA5 as a Novel Biomarker Associated with Alternative Splicing, Tumor Microenvironment, and Cell Proliferation in Colon Cancer Via Pan-cancer Analysis

Background: CDCA5 has been reported as a gene involved in the cell cycle, however current research provides little details. Our goal was to figure out its functions and probable mechanisms in pan-cancer. Methods: Pan-cancer bulk sequencing data and web-based analysis tools were applied to analyze CDCA5's correlations with the gene expression, clinical prognosis, genetic alterations, promoter methylation, alternative splicing, immune checkpoints, tumor microenvironment and enrichment. Real‑time PCR, cell clone formation assay, CCK-8 assay, cell proliferation assay, migration assay, invasion assay and apoptosis assay were used to evaluate the effect of CDCA5 silencing on colon cancer cell lines. Results: CDCA5 is highly expressed in most tumors, which has been linked to a poor prognosis. Immune checkpoints analysis revealed that CDCA5 was associated with the immune gene CD276 in various tumors. Single-cell analysis showed that CDCA5 correlated with proliferating T cell infiltration in COAD. Enrichment analysis demonstrated that CDCA5 may modify cell cycle genes to influence p53 signaling. The examination of DLD1 cells revealed that CDCA5 increased the proliferation and blocked cell apoptosis. Conclusion: This study contributes to the knowledge of the role of CDCA5 in carcinogenesis, highlighting the prognostic potential and carcinogenic involvement of CDCA5 in pan-cancer.

ESRP1-driven alternative splicing of CLSTN1 inhibits the metastasis of gastric cancer

Tumor metastasis severely limits the prognosis of gastric cancer patients. RNA-binding proteins (RBPs) are crucial in tumor metastasis, yet there is limited research into their involvement in gastric cancer. Here, we found that ESRP1, a RBP specific in epithelial cells, is important in regulating the metastasis of gastric cancer cells. ESRP1 is negatively correlated with distant metastasis and lymph node metastasis in gastric cancer patients. And we demonstrated that ESRP1 inhibit migration and invasion of gastric cancer in vitro and in vivo. Mechanistically, ESRP1 promotes exon 11 alternative splicing of CLSTN1 pre-mRNA. The post-splicing short CLSTN1 stabilizes the Ecadherin/β-catenin binding structure, and promotes β-catenin protein ubiquitination and degradation, thereby inhibiting the migration and invasion of gastric cancer cells. Our study highlights the role of ESRP1 in regulating metastasis of gastric cancer and extends its mechanism. These results provide a possibility for ESRP1 and CLSTN1 to become therapeutic targets for metastasis of gastric cancer.

Comparative RNA-Seq Analysis Revealed Tissue-Specific Splicing Variations during the Generation of the PDX Model

Tissue-specific gene expression generates fundamental differences in the function of each tissue and affects the characteristics of the tumors that are created as a result. However, it is unclear how much the tissue specificity is conserved during grafting of the primary tumor into an immune-compromised mouse model. Here, we performed a comparative RNA-seq analysis of four different primary-patient derived xenograft (PDX) tumors. The analysis revealed a conserved RNA biotype distribution of primary-PDX pairs, as revealed by previous works. Interestingly, we detected significant changes in the splicing pattern of PDX, which was mainly comprised of skipped exons. This was confirmed by splicing variant-specific RT-PCR analysis. On the other hand, the correlation analysis for the tissue-specific genes indicated overall strong positive correlations between the primary and PDX tumor pairs, with the exception of gastric cancer cases, which showed an inverse correlation. These data propose a tissue-specific change in splicing events during PDX formation as a variable factor that affects primary-PDX integrity.

Targeting PRMT1-mediated SRSF1 methylation to suppress oncogenic exon inclusion events and breast tumorigenesis

PRMT1 plays a vital role in breast tumorigenesis; however, the underlying molecular mechanisms remain incompletely understood. Herein, we show that PRMT1 plays a critical role in RNA alternative splicing, with a preference for exon inclusion. PRMT1 methylome profiling identifies that PRMT1 methylates the splicing factor SRSF1, which is critical for SRSF1 phosphorylation, SRSF1 binding with RNA, and exon inclusion. In breast tumors, PRMT1 overexpression is associated with increased SRSF1 arginine methylation and aberrant exon inclusion, which are critical for breast cancer cell growth. In addition, we identify a selective PRMT1 inhibitor, iPRMT1, which potently inhibits PRMT1-mediated SRSF1 methylation, exon inclusion, and breast cancer cell growth. Combination treatment with iPRMT1 and inhibitors targeting SRSF1 phosphorylation exhibits an additive effect of suppressing breast cancer cell growth. In conclusion, our study dissects a mechanism underlying PRMT1-mediated RNA alternative splicing. Thus, PRMT1 has great potential as a therapeutic target in breast cancer treatment.

SpliceWiz: interactive analysis and visualization of alternative splicing in R

Alternative splicing (AS) is a crucial mechanism for regulating gene expression and isoform diversity in eukaryotes. However, the analysis and visualization of AS events from RNA sequencing data remains challenging. Most tools require a certain level of computer literacy and the available means of visualizing AS events, such as coverage and sashimi plots, have limitations and can be misleading. To address these issues, we present SpliceWiz, an R package with an interactive Shiny interface that allows easy and efficient AS analysis and visualization at scale. A novel normalization algorithm is implemented to aggregate splicing levels within sample groups, thereby allowing group differences in splicing levels to be accurately visualized. The tool also offers downstream gene ontology enrichment analysis, highlighting ASEs belonging to functional pathways of interest. SpliceWiz is optimized for speed and efficiency and introduces a new file format for coverage data storage that is more efficient than BigWig. Alignment files are processed orders of magnitude faster than other R-based AS analysis tools and on par with command-line tools. Overall, SpliceWiz streamlines AS analysis, enabling reliable identification of functionally relevant AS events for further characterization. SpliceWiz is a Bioconductor package and is also available on GitHub (https://github.com/alexchwong/SpliceWiz).

DASES: a database of alternative splicing for esophageal squamous cell carcinoma

Esophageal carcinoma ranks as the sixth leading cause of cancer-related mortality globally, with esophageal squamous cell carcinoma (ESCC) being particularly prevalent among Asian populations. Alternative splicing (AS) plays a pivotal role in ESCC development and progression by generating diverse transcript isoforms. However, the current landscape lacks a specialized database focusing on alternative splicing events (ASEs) derived from a large number of ESCC cases. Additionally, most existing AS databases overlook the contribution of long non-coding RNAs (lncRNAs) in ESCC molecular mechanisms, predominantly focusing on mRNA-based ASE identification. To address these limitations, we deployed DASES (http://www.hxdsjzx.cn/DASES). Employing a combination of publicly available and in-house ESCC RNA-seq datasets, our extensive analysis of 346 samples, with 93% being paired tumor and adjacent non-tumor tissues, led to the identification of 257 novel lncRNAs in esophageal squamous cell carcinoma. Leveraging a paired comparison of tumor and adjacent normal tissues, DASES identified 59,094 ASEs that may be associated with ESCC. DASES fills a critical gap by providing comprehensive insights into ASEs in ESCC, encompassing lncRNAs and mRNA, thus facilitating a deeper understanding of ESCC molecular mechanisms and serving as a valuable resource for ESCC research communities.

Multi-omics analysis reveals interferon-stimulated gene OAS1 as a prognostic and immunological biomarker in pan-cancer

Introduction

OAS1(2'-5'-oligoadenylate synthetase 1) is a member of the Interferon-Stimulated Genes which plays an important role in the antiviral process. In recent years, the role of OAS1 in tumors has attracted attention, and it was found to be associated with prognosis in several tumors. However, the mechanism by which OAS1 affects tumors is unclear and pan-cancer study of OAS1 is necessary to better understand its implication in cancers.

Methods

The expression, prognostic value, genetic alteration, alternative splicing events of OAS1 in pan-cancers were analyzed using TCGA, GTEx, HPA, GEPIA and OncoSplicing databases. OAS1 associated immune cell infiltration was evaluated using the ESTIMATE, xCell, CIBERSORT and QUANTISEQ algorithm. Single cell transcriptome data download using TISH database. Finally, the roles of the OAS1 on apoptosis, migration and invasion were investigated in two pancreatic cancer cells.

Results

Our results revealed significant differences in OAS1 expression among various tumors, which had prognostic implications. In addition, we investigated the impact of OAS1 on genomic stability, methylation status, and other factors across different types of cancer, and the effects of these factors on prognosis. Notably, our study also demonstrated that OAS1 overexpression can contribute to CTL dysfunction and macrophage M2 polarization. In addition, cell experiments showed that the knockdown of OAS1 could reduce the invasive ability and increased the apoptosis rate of PAAD cells.

Discussion

These results confirmed that OAS1 could be a prognostic biomarker and therapeutic target for its potential role in CTL dysfunction and macrophage M2 polarization.

An Inflammatory Checkpoint Generated by IL1RN Splicing Offers Therapeutic Opportunity for KRAS-Mutant Intrahepatic Cholangiocarcinoma

KRAS mutations are causally linked to protumor inflammation and are identified as driving factors in tumorigenesis. Here, using multiomics data gathered from a large set of patients, we showed that KRAS mutation was associated with a specific landscape of alternative mRNA splicing that connected to myeloid inflammation in intrahepatic cholangiocarcinoma (iCCA). Then, we identified a negative feedback mechanism in which the upregulation of interleukin 1 receptor antagonist (IL1RN)-201/203 due to alternative splicing confers vital anti-inflammatory effects in KRAS-mutant iCCA. In KRAS-mutant iCCA mice, both IL1RN-201/203 upregulation and anakinra treatment ignited a significant antitumor immune response by altering neutrophil recruitment and phenotypes. Furthermore, anakinra treatment synergistically enhanced anti-PD-1 therapy to activate intratumoral GZMB+ CD8+ T cells in KRAS-mutant iCCA mice. Clinically, we found that high IL1RN-201/203 levels in patients with KRAS-mutant iCCA were significantly associated with superior response to anti-PD-1 immunotherapy.

SIGNIFICANCE

This work describes a novel inflammatory checkpoint mediated by IL1RN alternative splicing variants that may serve as a promising basis to develop therapeutic options for KRAS-mutant iCCA and other cancers. This article is featured in Selected Articles from This Issue, p. 2109.

Identification of alternative splicing associated with clinical features: from pan-cancers to genitourinary tumors

Background

Alternative splicing events (ASEs) are vital causes of tumor heterogeneity in genitourinary tumors and many other cancers. However, the clinicopathological relevance of ASEs in cancers has not yet been comprehensively characterized.

Methods

By analyzing splicing data from the TCGA SpliceSeq database and phenotype data for all TCGA samples from the UCSC Xena database, we identified differential clinical feature-related ASEs in 33 tumors. CIBERSORT immune cell infiltration data from the TIMER2.0 database were used for differential clinical feature-related immune cell infiltration analysis. Gene function enrichment analysis was used to analyze the gene function of ASEs related to different clinical features in tumors. To reveal the regulatory mechanisms of ASEs, we integrated race-related ASEs and splicing quantitative trait loci (sQTLs) data in kidney renal clear cell carcinoma (KIRC) to comprehensively assess the impact of SNPs on ASEs. In addition, we predicted regulatory RNA binding proteins in bladder urothelial carcinoma (BLCA) based on the enrichment of motifs around alternative exons for ASEs.

Results

Alternative splicing differences were systematically analyzed between different groups of 58 clinical features in 33 cancers, and 30 clinical features in 24 cancer types were identified to be associated with more than 50 ASEs individually. The types of immune cell infiltration were found to be significantly different between subgroups of primary diagnosis and disease type. After integrating ASEs with sQTLs data, we found that 63 (58.9%) of the race-related ASEs were significantly SNP-correlated ASEs in KIRC. Gene function enrichment analyses showed that metastasis-related ASEs in KIRC mainly enriched Rho GTPase signaling pathways. Among those ASEs associated with metastasis, alternative splicing of GIT2 and TUBB3 might play key roles in tumor metastasis in KIRC patients. Finally, we identified several RNA binding proteins such as PCBP2, SNRNP70, and HuR, which might contribute to splicing differences between different groups of neoplasm grade in BLCA.

Conclusion

We demonstrated the significant clinical relevance of ASEs in multiple cancer types. Furthermore, we identified and validated alternative splicing of TUBB3 and RNA binding proteins such as PCBP2 as critical regulators in the progression of urogenital cancers.

Identification of the H3K36me3 reader LEDGF/p75 in the pancancer landscape and functional exploration in clear cell renal cell carcinoma

Lens epithelium-derived growth factor (LEDGF/p75) is a reader of epigenetic marks and a potential target for therapeutic intervention. Its involvement in human immunodeficiency virus (HIV) integration and the development of leukemia driven by MLL (also known as KMT2A) gene fusion make it an attractive candidate for drug development. However, exploration of LEDGF/p75 as an epigenetic reader of H3K36me3 in tumors is limited. Here, for the first time, we analyze the role of LEDGF/p75 in multiple cancers via multiple online databases and in vitro experiments. We used pancancer bulk sequencing data and online tools to analyze correlations of LEDGF/p75 with prognosis, genomic instability, DNA damage repair, prognostic alternative splicing, protein interactions, and tumor immunity. In summary, the present study identified that LEDGF/p75 may serve as a prognostic predictor for tumors such as adrenocortical carcinoma, kidney chromophobe, liver hepatocellular carcinoma, pancreatic adenocarcinoma, skin cutaneous melanoma, and clear cell renal cell carcinoma (ccRCC). In addition, in vitro experiments and gene microarray sequencing were performed to explore the function of LEDGF/p75 in ccRCC, providing new insights into the pathogenesis of the nonmutated SETD2 ccRCC subtype.

Protein diversification through post-translational modifications, alternative splicing, and gene duplication

Proteins provide the basis for cellular function. Having multiple versions of the same protein within a single organism provides a way of regulating its activity or developing novel functions. Post-translational modifications of proteins, by means of adding/removing chemical groups to amino acids, allow for a well-regulated and controlled way of generating functionally distinct protein species. Alternative splicing is another method with which organisms possibly generate new isoforms. Additionally, gene duplication events throughout evolution generate multiple paralogs of the same genes, resulting in multiple versions of the same protein within an organism. In this review, we discuss recent advancements in the study of these three methods of protein diversification and provide illustrative examples of how they affect protein structure and function.

Alternative splicing: a new breakthrough for understanding tumorigenesis and potential clinical applications

BACKGROUND

Alternative splicing (AS) is a post-transcriptional process that produces transcript variants, thus leading to transcriptome complexity. Recently, the scope of AS studies has been greatly expanded toward clinical applications owing to the abundance of RNA sequencing data.

OBJECTIVE

This review consists of two parts. We first summarize bioinformatic resources that are useful for large-scale cancer-related AS studies. We then highlight the research efforts to utilize AS events for predicting clinical outcomes and planning therapeutic strategies.

RESULTS

Computational approaches to interrogate AS events have been reviewed under three categories: (1) databases to provide functional and clinical annotation of AS events, (2) analytical tools to identify cancer-associated AS event, and (3) methods to identify splicing-related DNA variants and splicing-derived neoantigens. We also present the recent progress in exploring the clinical utility of AS under four categories: (1) identification of AS events for cancer prognosis, (2) utilization of AS events in molecular classification of various cancers, (3) regulatory mechanisms of AS underlying drug resistance, and (4) potential use of AS in cancer therapy.

CONCLUSION

This review will be helpful for understanding the biological implications of AS in cancer and facilitate the development of AS markers for cancer prognosis and treatment. We anticipate that future studies will lead to the application of genome-wide AS profiles in cancer precision medicine.

The short isoform of MS4A7 is a novel player in glioblastoma microenvironment, M2 macrophage polarization, and tumor progression

BACKGROUND

The unique intracranial tumor microenvironment (TME) contributes to the immunotherapy failure for glioblastoma (GBM), thus new functional protein targets are urgently needed. Alternative splicing is a widespread regulatory mechanism by which individual gene can express variant proteins with distinct functions. Moreover, proteins located in the cell plasma membrane facilitate targeted therapies. This study sought to obtain functional membrane protein isoforms from GBM TME.

METHODS

With combined single-cell RNA-seq and bulk RNA-seq analyses, novel candidate membrane proteins generated by prognostic splicing events were screened within GBM TME. The short isoform of MS4A7 (MS4A7-s) was selected for evaluation by RT-PCR and western blotting in clinical specimens. Its clinical relevance was evaluated in a GBM patient cohort. The function of MS4A7-s was identified by in vitro and in vivo experiments. MS4A7-s overexpression introduced transcriptome changes were analyzed to explore the potential molecular mechanism.

RESULTS

The main expression product, isoform MS4A7-s, generated by exon skipping, is an M2-specific plasma membrane protein playing a pro-oncogenic role in GBM TME. Higher expression of MS4A7-s correlates with poor prognosis in a GBM cohort. In vitro cell co-culture experiments, intracranial co-injection tumorigenesis assay, and RNA-seq suggest MS4A7-s promotes activation of glioma-associated macrophages' (GAMs) PI3K/AKT/GSK3β pathway, leading to M2 polarization, and drives malignant progression of GBM.

CONCLUSIONS

MS4A7-s, a novel splicing isoform of MS4A7 located on the surface of GAMs in GBM TME, is a predictor of patient outcome, which contributes to M2 polarization and the malignant phenotype of GBM. Targeting MS4A7-s may constitute a promising treatment for GBM.

ChromLoops: a comprehensive database for specific protein-mediated chromatin loops in diverse organisms

Chromatin loops (or chromatin interactions) are important elements of chromatin structures. Disruption of chromatin loops is associated with many diseases, such as cancer and polydactyly. A few methods, including ChIA-PET, HiChIP and PLAC-Seq, have been proposed to detect high-resolution, specific protein-mediated chromatin loops. With rapid progress in 3D genomic research, ChIA-PET, HiChIP and PLAC-Seq datasets continue to accumulate, and effective collection and processing for these datasets are urgently needed. Here, we developed a comprehensive, multispecies and specific protein-mediated chromatin loop database (ChromLoops, https://3dgenomics.hzau.edu.cn/chromloops), which integrated 1030 ChIA-PET, HiChIP and PLAC-Seq datasets from 13 species, and documented 1 491 416 813 high-quality chromatin loops. We annotated genes and regions overlapping with chromatin loop anchors with rich functional annotations, such as regulatory elements (enhancers, super-enhancers and silencers), variations (common SNPs, somatic SNPs and eQTLs), and transcription factor binding sites. Moreover, we identified genes with high-frequency chromatin interactions in the collected species. In particular, we identified genes with high-frequency interactions in cancer samples. We hope that ChromLoops will provide a new platform for studying chromatin interaction regulation in relation to biological processes and disease.

AS-CMC: a pan-cancer database of alternative splicing for molecular classification of cancer

Alternative splicing (AS) is a post-transcriptional regulation that leads to the complexity of the transcriptome. Despite the growing importance of AS in cancer research, the role of AS has not been systematically studied, especially in understanding cancer molecular classification. Herein, we analyzed the molecular subtype-specific regulation of AS using The Cancer Genome Atlas data and constructed a web-based database, named Alternative Splicing for Cancer Molecular Classification (AS-CMC). Our system harbors three analysis modules for exploring subtype-specific AS events, evaluating their phenotype association, and performing pan-cancer comparison. The number of subtype-specific AS events was found to be diverse across cancer types, and some differentially regulated AS events were recurrently found in multiple cancer types. We analyzed a subtype-specific AS in exon 11 of mitogen-activated protein kinase kinase 7 (MAP3K7) as an example of a pan-cancer AS biomarker. This AS marker showed significant association with the survival of patients with stomach adenocarcinoma. Our analysis revealed AS as an important determinant for cancer molecular classification. AS-CMC is the first web-based resource that provides a comprehensive tool to explore the biological implications of AS events, facilitating the discovery of novel AS biomarkers.

ASCancer Atlas: a comprehensive knowledgebase of alternative splicing in human cancers

Alternative splicing (AS) is a fundamental process that governs almost all aspects of cellular functions, and dysregulation in this process has been implicated in tumor initiation, progression and treatment resistance. With accumulating studies of carcinogenic mis-splicing in cancers, there is an urgent demand to integrate cancer-associated splicing changes to better understand their internal cross-talks and functional consequences from a global view. However, a resource of key functional AS events in human cancers is still lacking. To fill the gap, we developed ASCancer Atlas (https://ngdc.cncb.ac.cn/ascancer), a comprehensive knowledgebase of aberrant splicing in human cancers. Compared to extant databases, ASCancer Atlas features a high-confidence collection of 2006 cancer-associated splicing events experimentally proved to promote tumorigenesis, a systematic splicing regulatory network, and a suit of multi-scale online analysis tools. For each event, we manually curated the functional axis including upstream splicing regulators, splicing event annotations, downstream oncogenic effects, and possible therapeutic strategies. ASCancer Atlas also houses about 2 million computationally putative splicing events. Additionally, a user-friendly web interface was built to enable users to easily browse, search, visualize, analyze, and download all splicing events. Overall, ASCancer Atlas provides a unique resource to study the functional roles of splicing dysregulation in human cancers.

A Pan-Cancer Analysis Reveals CLEC5A as a Biomarker for Cancer Immunity and Prognosis

Background

CLEC5A is a member of the C-type lectin superfamily. It can activate macrophages and lead to a series of immune-inflammation reactions. Previous studies reveal the role of CLEC5A in infection and inflammation diseases.

Method

We acquire and analyze data from The Cancer Genome Atlas (TCGA) database, Genotype-Tissue Expression (GTEx) database, and other comprehensive databases via GSCALite, cBioPortal, and TIMER 2.0 platforms or software. Single-cell sequencing analysis was performed for quantifying the tumor microenvironment of several types of cancers.

Results

CLEC5A is differentially expressed in a few cancer types, of which overexpression accompanies low overall survival of patients. DNA methylation mainly negatively correlates with CLEC5A expression. Moreover, CLEC5A is positively related to immune infiltration, including macrophages, cancer-associated fibroblasts (CAFs), and regulatory T cells (Tregs). Immune checkpoint genes are significantly associated with CLEC5A expression in diverse cancers. In addition, CLEC5A expression correlates with mismatch repair (MMR) in several cancers. Tumor mutation burden (TMB), microsatellite instability (MSI), and neoantigens show a positive association with CLEC5A expression in several cancers. Furthermore, CLEC5A in cancer correlates with signal transduction, the immune system, EMT, and apoptosis process. The drug sensitivity analysis screens out potential therapeutic agents associated with CLEC5A expression, including FR-180204, Tivozanib, OSI-930, Linifanib, AC220, VNLG/124, Bexarotene, omacetaxine mepesuccinate, narciclasine, leptomycin B, PHA-793887, LRRK2-IN-1, and CR-1-31B.

Conclusion

CLEC5A overexpresses in multiple cancers in contrast to normal tissues, and high CLEC5A expression predicts poor prognosis of patients and immune infiltration. CLEC5A is a potential prognostic biomarker of diverse cancers and a target for anti-tumor therapy.

LncRNA RP3-525N10.2-NFKB1-PROS1 triplet-mediated low PROS1 expression is an onco-immunological biomarker in low-grade gliomas: a pan-cancer analysis with experimental verification

Background

Glioma is the most common cancer in the central nervous system, and low grade gliomas are notorious for many types of tumors and heterogeneity. PROS1 not only plays an important role in the blood coagulation system, and recent studies have found that it was correlated with the development of tumors, especially related to tumor immune infiltration. However, the study of underlying role and mechanism of PROS1 in gliomas, especially in low-grade gliomas, is almost absent.

Methods

We integrated the information of patients with LGG in The Cancer Genome Atlas (TCGA) cohort and Chinese Glioma Genome Atlas (CGGA) cohort. Then, we systematically demonstrated the differences and prognostic prognosis value of PROS1 based on multi-omics analyses. In addition, Cell counting kit-8 (CCK-8) assay, colony formation assay, 5-Ethynyl-2’-deoxyuridine (EdU) incorporation assay, and Transwell assays were performed to evaluate cell proliferation and invasion. qRT-PCR and immunohistochemistry were used to evaluate the expression of PROS1 in LGG.

Results

Various bioinformatics approaches revealed that PROS1 was a valuable prognostic marker and may influence tumour development via distinct mechanisms, including expression of DNA methyltransferase, RNA modification, and DNA mismatch repair system genes, copy number variation, single nucleotide variation frequency, genomic heterogeneity, cancer stemness, DNA methylation, and alternative PROS1 splicing. Our analyses indicated that the long non-coding RNA RP3-525N10.2 may “decoy” or “guide” the transcription factor NFKB1 and prevent its association with PROS1, thereby reducing PROS1 expression and improving poor LGG prognosis. PROS1 expression was also closely associated with tumour infiltration by immune cells, especially tumour-associated macrophages, as well as the expression of various immune checkpoint inhibitors, immunomodulators, and immune cell markers.

Conclusion

long non-coding RNA RP3-525N10.2-NFKB1-PROS1 triplet-mediated PROS1 expression could serve as a biomarker for cancer diagnosis, prognosis, therapy selection, and follow-up in LGG patients.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12967-022-03536-y.

JMJD8 Is an M2 Macrophage Biomarker, and It Associates With DNA Damage Repair to Facilitate Stemness Maintenance, Chemoresistance, and Immunosuppression in Pan-Cancer

Background

JMJD8 has recently been identified as a cancer-related gene, but current studies provide limited information. We aimed to clarify its roles and the potential mechanisms in pan-cancer.

Methods

Pan-cancer bulk sequencing data and online web tools were applied to analyze JMJD8’s correlations with prognosis, genome instability, cancer stemness, DNA repair, and immune infiltration. Moreover, single-cell datasets, SpatialDB database, and multiple fluorescence staining were used to validate the association between JMJD8 expression and M2 macrophages. Further, we utilized ROCplotter and cMap web tool to analyze the therapeutic responses and screened JMJD8-targeted compounds, respectively, and we used AlphaFold2 and Discovery Studio to conduct JMJD8 homology modeling and molecular docking.

Results

We first noticed that JMJD8 was an oncogene in many cancer types. High JMJD8 was associated with lower genome stability. We then found that high JMJD8 correlated with high expression of mismatch repair genes, stemness, homologous repair gene signature in more than 9 cancers. ESTIMATE and cytokine analyses results presented JMJD8’s association with immunosuppression. Also, immune checkpoint CD276 was positively relevant to JMJD8. Subsequently, we validated JMJD8 as the M2 macrophage marker and showed its connection with other immunosuppressive cells and CD8+ T-cell depression. Finally, potential JMJD8-targeted drugs were screened out and docked to JMJD8 protein.

Conclusion

We found that JMJD8 was a novel oncogene, and it correlated with immunosuppression and DNA repair. JMJD8 was highly associated with immune checkpoint CD276 and was an M2 macrophage biomarker in many cancers. This study will reveal JMJD8’s roles in pan-cancer and its potential as a novel therapeutic target.