SON drives oncogenic RNA splicing in glioblastoma by regulating PTBP1/PTBP2 switching and RBFOX2 activity

Comprehensive machine learning-based integration develops a novel prognostic model for glioblastoma

In this study, we developed a new prognostic model for glioblastoma (GBM) based on an integrated machine learning algorithm. We used univariate Cox regression analysis to identify prognostic genes by combining six GBM cohorts. Based on the prognostic genes, 10 machine learning algorithms were integrated into 117 algorithm combinations, and the artificial intelligence prognostic signature (AIPS) with the greatest average C-index was chosen. The AIPS was compared with 10 previously published models by univariate Cox analysis and the C-index. We compared the differences in prognosis, tumor immune microenvironment (TIME), and immunotherapy sensitivity between the high and low AIPS score groups. The AIPS based on the random survival forest algorithm with the highest average C-index (0.868) was selected. Compared with the previous 10 prognostic models, our AIPS has the highest C-index. The AIPS was closely linked to the clinical features of GBM. We discovered that patients in the low score group had improved prognoses, a more active TIME, and were more sensitive to immunotherapy. Finally, we verified the expression of several key genes by western blotting and immunohistochemistry. We identified an ideal prognostic signature for GBM, which might provide new insights into stratified treatment approaches for GBM patients.

Long-read sequencing for brain tumors

Brain tumors and genomics have a long-standing history given that glioblastoma was the first cancer studied by the cancer genome atlas. The numerous and continuous advances through the decades in sequencing technologies have aided in the advanced molecular characterization of brain tumors for diagnosis, prognosis, and treatment. Since the implementation of molecular biomarkers by the WHO CNS in 2016, the genomics of brain tumors has been integrated into diagnostic criteria. Long-read sequencing, also known as third generation sequencing, is an emerging technique that allows for the sequencing of longer DNA segments leading to improved detection of structural variants and epigenetics. These capabilities are opening a way for better characterization of brain tumors. Here, we present a comprehensive summary of the state of the art of third-generation sequencing in the application for brain tumor diagnosis, prognosis, and treatment. We discuss the advantages and potential new implementations of long-read sequencing into clinical paradigms for neuro-oncology patients.

Tumour microenvironment programming by an RNA-RNA-binding protein complex creates a druggable vulnerability in IDH-wild-type glioblastoma

Patients with IDH-wild-type glioblastomas have a poor five-year survival rate along with limited treatment efficacy due to immune cell (glioma-associated microglia and macrophages) infiltration promoting tumour growth and resistance. To enhance therapeutic options, our study investigated the unique RNA-RNA-binding protein complex LOC-DHX15. This complex plays a crucial role in driving immune cell infiltration and tumour growth by establishing a feedback loop between cancer and immune cells, intensifying cancer aggressiveness. Targeting this complex with blood-brain barrier-permeable small molecules improved treatment efficacy, disrupting cell communication and impeding cancer cell survival and stem-like properties. Focusing on RNA-RNA-binding protein interactions emerges as a promising approach not only for glioblastomas without the IDH mutation but also for potential applications beyond cancer, offering new avenues for developing therapies that address intricate cellular relationships in the body.

Integrated pan-cancer analysis reveals the immunological and prognostic potential of RBFOX2 in human tumors

Background

The role of RNA-binding fox one homolog 2 (RBFOX2) in the progression of multiple tumors is increasingly supported by evidence. However, the unclearness pertaining to the expression of RBFOX2, its prognostic potential, and its correlation with the tumor microenvironment (TME) in pan-cancer persists. This study aims to comprehensively investigate the immunological prognostic value of RBFOX2.

Methods

The Cancer Genome Atlas Gene Expression Omnibus Genotype-Tissue Expression (GTEx), TIMER2.0, Kaplan-Meier (K-M) Plotter, University of Alabama at Birmingham Cancer data analysis Portal (UALCAN), cbioportal, and Gene Expression Profiling Interactive Analysis 2 (GEPIA2) were utilized for a systematic analysis of RBFOX2. This analysis included studying its expression, prognostic value, DNA methylation, enrichment analysis, immune infiltration cells, and immune-related genes. Additionally, qRT-PCR, CCK-8, colony formation, transwell assays, and immunohistochemistry were employed to analyze the expression and biological function of RBFOX2 in liver cancer.

Results

Variations in RBFOX2 expression have been observed across diverse tumors and have been identified as indicators of unfavorable prognosis. It is closely linked to immune infiltration cells, immune checkpoints, chemokines, and chemokine receptors in the TME. Higher levels of RBFOX2 have been significantly associated with low response and poor prognosis in patients with non-small cell lung cancer (NSCLC) and melanoma who receive immunotherapy. Furthermore, the DNA methylation of RBFOX2 varies across different types of cancer and has shown better prognosis in patients with BLCA, BRCA, CESC, COAD, DLBC, HNSC, LAML, LGG, LUAD, PAAD, SKCM and THYM. Interestingly, RBFOX2 expression was found to be lower in hepatocellular carcinoma (HCC) patients' tumor tissues compared to their paired adjacent tissues. In vitro studies have shown that knockdown of RBFOX2 significantly promotes the growth and metastasis of liver cancer cells.

Conclusion

This study investigates the correlation between DNA methylation, prognostic value, and immune cell infiltration with the expression of RBFOX2 in pan-cancer and indicates its potential role to inhibit metastasis of liver cancer.

MAPP unravels frequent co-regulation of splicing and polyadenylation by RNA-binding proteins and their dysregulation in cancer

Maturation of eukaryotic pre-mRNAs via splicing and polyadenylation is modulated across cell types and conditions by a variety of RNA-binding proteins (RBPs). Although there exist over 1,500 RBPs in human cells, their binding motifs and functions still remain to be elucidated, especially in the complex environment of tissues and in the context of diseases. To overcome the lack of methods for the systematic and automated detection of sequence motif-guided pre-mRNA processing regulation from RNA sequencing (RNA-Seq) data we have developed MAPP (Motif Activity on Pre-mRNA Processing). Applying MAPP to RBP knock-down experiments reveals that many RBPs regulate both splicing and polyadenylation of nascent transcripts by acting on similar sequence motifs. MAPP not only infers these sequence motifs, but also unravels the position-dependent impact of the RBPs on pre-mRNA processing. Interestingly, all investigated RBPs that act on both splicing and 3' end processing exhibit a consistently repressive or activating effect on both processes, providing a first glimpse on the underlying mechanism. Applying MAPP to normal and malignant brain tissue samples unveils that the motifs bound by the PTBP1 and RBFOX RBPs coordinately drive the oncogenic splicing program active in glioblastomas demonstrating that MAPP paves the way for characterizing pre-mRNA processing regulators under physiological and pathological conditions.

RNA splicing analysis deciphers developmental hierarchies and reveals therapeutic targets in adult glioma

Widespread alterations in RNA alternative splicing (AS) have been identified in adult gliomas. However, their regulatory mechanism, biological significance, and therapeutic potential remain largely elusive. Here, using a computational approach with both bulk and single-cell RNA-Seq, we uncover a prognostic AS signature linked with neural developmental hierarchies. Using advanced iPSC glioma models driven by glioma driver mutations, we show that this AS signature could be enhanced by EGFRvIII and inhibited by in situ IDH1 mutation. Functional validations of 2 isoform switching events in CERS5 and MPZL1 show regulations of sphingolipid metabolism and SHP2 signaling, respectively. Analysis of upstream RNA binding proteins reveals PTBP1 as a key regulator of the AS signature where targeting of PTBP1 suppresses tumor growth and promotes the expression of a neuron marker TUJ1 in glioma stem-like cells. Overall, our data highlights the role of AS in affecting glioma malignancy and heterogeneity and its potential as a therapeutic vulnerability for treating adult gliomas.

Nuclear speckle rejuvenation alleviates proteinopathies at the expense of YAP1

Current treatments targeting individual protein quality control have limited efficacy in alleviating proteinopathies, highlighting the prerequisite for a common upstream druggable target capable of global proteostasis modulation. Building on our prior research establishing nuclear speckles as pivotal organelles responsible for global proteostasis transcriptional control, we aim to alleviate proteinopathies through nuclear speckle rejuvenation. We identified pyrvinium pamoate as a small-molecule nuclear speckle rejuvenator that enhances protein quality control while suppressing YAP1 signaling via decreasing the surface tension of nuclear speckle condensates through interaction with the intrinsically disordered region of nuclear speckle scaffold protein SON. In pre-clinical models, pyrvinium pamoate reduced tauopathy and alleviated retina degeneration by promoting autophagy and ubiquitin-proteasome system. Aberrant nuclear speckle morphology, reduced protein quality control and increased YAP1 activity were also observed in human tauopathies. Our study uncovers novel therapeutic targets for tackling protein misfolding disorders within an expanded proteostasis framework encompassing nuclear speckles and YAP1.

Renal Endothelial Single-Cell Transcriptomics Reveals Spatiotemporal Regulation and Divergent Roles of Differential Gene Transcription and Alternative Splicing in Murine Diabetic Nephropathy

Endothelial cell (EC) injury is a crucial contributor to the progression of diabetic kidney disease (DKD), but the specific EC populations and mechanisms involved remain elusive. Kidney ECs (n = 5464) were collected at three timepoints from diabetic BTBRob/ob mice and non-diabetic littermates. Their heterogeneity, transcriptional changes, and alternative splicing during DKD progression were mapped using SmartSeq2 single-cell RNA sequencing (scRNAseq) and elucidated through pathway, network, and gene ontology enrichment analyses. We identified 13 distinct transcriptional EC phenotypes corresponding to different kidney vessel subtypes, confirmed through in situ hybridization and immunofluorescence. EC subtypes along nephrons displayed extensive zonation related to their functions. Differential gene expression analyses in peritubular and glomerular ECs in DKD underlined the regulation of DKD-relevant pathways including EIF2 signaling, oxidative phosphorylation, and IGF1 signaling. Importantly, this revealed the differential alteration of these pathways between the two EC subtypes and changes during disease progression. Furthermore, glomerular and peritubular ECs also displayed aberrant and dynamic alterations in alternative splicing (AS), which is strongly associated with DNA repair. Strikingly, genes displaying differential transcription or alternative splicing participate in divergent biological processes. Our study reveals the spatiotemporal regulation of gene transcription and AS linked to DKD progression, providing insight into pathomechanisms and clues to novel therapeutic targets for DKD treatment.

NSUN5/TET2-directed chromatin-associated RNA modification of 5-methylcytosine to 5-hydroxymethylcytosine governs glioma immune evasion

Malignant glioma exhibits immune evasion characterized by highly expressing the immune checkpoint CD47. RNA 5-methylcytosine(m5C) modification plays a pivotal role in tumor pathogenesis. However, the mechanism underlying m5C-modified RNA metabolism remains unclear, as does the contribution of m5C-modified RNA to the glioma immune microenvironment. In this study, we demonstrate that the canonical 28SrRNA methyltransferase NSUN5 down-regulates β-catenin by promoting the degradation of its mRNA, leading to enhanced phagocytosis of tumor-associated macrophages (TAMs). Specifically, the NSUN5-induced suppression of β-catenin relies on its methyltransferase activity mediated by cysteine 359 (C359) and is not influenced by its localization in the nucleolus. Intriguingly, NSUN5 directly interacts with and deposits m5C on CTNNB1 caRNA (chromatin-associated RNA). NSUN5-induced recruitment of TET2 to chromatin is independent of its methyltransferase activity. The m5C modification on caRNA is subsequently oxidized into 5-hydroxymethylcytosine (5hmC) by TET2, which is dependent on its binding affinity for Fe2+ and α-KG. Furthermore, NSUN5 enhances the chromatin recruitment of RBFOX2 which acts as a 5hmC-specific reader to recognize and facilitate the degradation of 5hmC caRNA. Notably, hmeRIP-seq analysis reveals numerous mRNA substrates of NSUN5 that potentially undergo this mode of metabolism. In addition, NSUN5 is epigenetically suppressed by DNA methylation and is negatively correlated with IDH1-R132H mutation in glioma patients. Importantly, pharmacological blockage of DNA methylation or IDH1-R132H mutant and CD47/SIRPα signaling synergistically enhances TAM-based phagocytosis and glioma elimination in vivo. Our findings unveil a general mechanism by which NSUN5/TET2/RBFOX2 signaling regulates RNA metabolism and highlight NSUN5 targeting as a potential strategy for glioma immune therapy.

The leader RNA of SARS-CoV-2 sequesters polypyrimidine tract binding protein (PTBP1) and influences pre-mRNA splicing in infected cells

The large amount of viral RNA produced during infections has the potential to interact with and effectively sequester cellular RNA binding proteins, thereby influencing aspects of post-transcriptional gene regulation in the infected cell. Here we demonstrate that the abundant 5' leader RNA region of SARS-CoV-2 viral RNAs can interact with the cellular polypyrimidine tract binding protein (PTBP1). Interestingly, the effect of a knockdown of PTBP1 protein on cellular gene expression is also mimicked during SARS-CoV-2 infection, suggesting that this protein may be functionally sequestered by viral RNAs. Consistent with this model, the alternative splicing of mRNAs that is normally controlled by PTBP1 is dysregulated during SARS-CoV-2 infection. Collectively, these data suggest that the SARS-CoV-2 leader RNA sequesters the cellular PTBP1 protein during infection, resulting in significant impacts on the RNA biology of the host cell. These alterations in post-transcriptional gene regulation may play a role in SARS-CoV-2 mediated molecular pathogenesis.

Alternative splicing and related RNA binding proteins in human health and disease

Alternative splicing (AS) serves as a pivotal mechanism in transcriptional regulation, engendering transcript diversity, and modifications in protein structure and functionality. Across varying tissues, developmental stages, or under specific conditions, AS gives rise to distinct splice isoforms. This implies that these isoforms possess unique temporal and spatial roles, thereby associating AS with standard biological activities and diseases. Among these, AS-related RNA-binding proteins (RBPs) play an instrumental role in regulating alternative splicing events. Under physiological conditions, the diversity of proteins mediated by AS influences the structure, function, interaction, and localization of proteins, thereby participating in the differentiation and development of an array of tissues and organs. Under pathological conditions, alterations in AS are linked with various diseases, particularly cancer. These changes can lead to modifications in gene splicing patterns, culminating in changes or loss of protein functionality. For instance, in cancer, abnormalities in AS and RBPs may result in aberrant expression of cancer-associated genes, thereby promoting the onset and progression of tumors. AS and RBPs are also associated with numerous neurodegenerative diseases and autoimmune diseases. Consequently, the study of AS across different tissues holds significant value. This review provides a detailed account of the recent advancements in the study of alternative splicing and AS-related RNA-binding proteins in tissue development and diseases, which aids in deepening the understanding of gene expression complexity and offers new insights and methodologies for precision medicine.

A novel EGFR variant EGFRx maintains glioblastoma stem cells through STAT5

BACKGROUND

Glioblastomas are universally lethal brain tumors containing tumor-propagating glioblastoma stem cells (GSCs). EGFR gene amplification or mutation is frequently detected in GBMs and is associated with poor prognosis. However, EGFR variants in GSCs and their role in the maintenance of GSCs and progression of GBM are unclear.

METHODS

EGFR variants were detected through bioinformatic HISAT-StringTie-Ballgown pipeline and verified through 5' RACE, RT-PCR, ribonuclease protection, and northern blotting assays. EGFRx function was investigated through neurosphere, cell viability, intracranial xenograft and RNA-seq assays. EGFRx-STAT5 signaling was investigated through western blotting, coimmunoprecipitation, immunofluorescence, luciferase reporter, RT-PCR and CUT&Tag assays.

RESULTS

We identified a novel EGFR variant (EGFRx), that is specifically expressed in GSCs. Unlike the EGFRvIII variant, which lacks exons 2-7, EGFRx is characterized by the absence of exons 2-14, and encodes an EGFR protein that does not possess the entire extracellular ligand-binding domain. We observed that EGFRx exhibits significant glycosylation, is required for GSC self-renewal, proliferation, and tumorigenesis, and highly active in glioblastomas compared to normal brain tissue. Mechanistically, EGFRx constitutively and specifically activates STAT5 in GSCs through spontaneous asymmetric dimerization of the kinase domain.

CONCLUSIONS

EGFRx plays essential roles in the maintenance of the GSC phenotype through constitutive activation of STAT5 and promotes GBM progression, suggesting that EGFRx-STAT5 signaling represents a promising therapeutic target for GBM.

PTBP1 as a potential regulator of disease

Polypyrimidine tract-binding protein 1 (PTBP1) is a member of the heterogeneous nuclear ribonucleoprotein (hnRNP) family, which plays a key role in alternative splicing of precursor mRNA and RNA metabolism. PTBP1 is universally expressed in various tissues and binds to multiple downstream transcripts to interfere with physiological and pathological processes such as the tumor growth, body metabolism, cardiovascular homeostasis, and central nervous system damage, showing great prospects in many fields. The function of PTBP1 involves the regulation and interaction of various upstream molecules, including circular RNAs (circRNAs), microRNAs (miRNAs) and long non-coding RNAs (lncRNAs). These regulatory systems are inseparable from the development and treatment of diseases. Here, we review the latest knowledge regarding the structure and molecular functions of PTBP1 and summarize its functions and mechanisms of PTBP1 in various diseases, including controversial studies. Furthermore, we recommend future studies on PTBP1 and discuss the prospects of targeting PTBP1 in new clinical therapeutic approaches.

scMuffin: an R package to disentangle solid tumor heterogeneity by single-cell gene expression analysis

INTRODUCTION

Single-cell (SC) gene expression analysis is crucial to dissect the complex cellular heterogeneity of solid tumors, which is one of the main obstacles for the development of effective cancer treatments. Such tumors typically contain a mixture of cells with aberrant genomic and transcriptomic profiles affecting specific sub-populations that might have a pivotal role in cancer progression, whose identification eludes bulk RNA-sequencing approaches. We present scMuffin, an R package that enables the characterization of cell identity in solid tumors on the basis of a various and complementary analyses on SC gene expression data.

RESULTS

scMuffin provides a series of functions to calculate qualitative and quantitative scores, such as: expression of marker sets for normal and tumor conditions, pathway activity, cell state trajectories, Copy Number Variations, transcriptional complexity and proliferation state. Thus, scMuffin facilitates the combination of various evidences that can be used to distinguish normal and tumoral cells, define cell identities, cluster cells in different ways, link genomic aberrations to phenotypes and identify subtle differences between cell subtypes or cell states. We analysed public SC expression datasets of human high-grade gliomas as a proof-of-concept to show the value of scMuffin and illustrate its user interface. Nevertheless, these analyses lead to interesting findings, which suggest that some chromosomal amplifications might underlie the invasive tumor phenotype and the presence of cells that possess tumor initiating cells characteristics.

CONCLUSIONS

The analyses offered by scMuffin and the results achieved in the case study show that our tool helps addressing the main challenges in the bioinformatics analysis of SC expression data from solid tumors.

Aurora kinase A regulates cancer-associated RNA aberrant splicing in breast cancer

The contribution of oncogenes to tumor-associated RNA splicing and the relevant molecular mechanisms therein require further elaboration. Here, we show that oncogenic Aurora kinase A (AURKA) promotes breast cancer-related RNA aberrant splicing in a context-dependent manner. AURKA regulated pan-breast cancer-associated RNA splicing events including GOLGA4, RBM4 and UBQLN1. Aberrant splicing of GOLGA4 and RBM4 was closely related to breast cancer development. Mechanistically, AURKA interacted with the splicing factor YBX1 and promoted AURKA-YBX1 complex-mediated GOLGA4 exon inclusion. AURKA binding to the splicing factor hnRNPK promoted AURKA-hnRNPK complex-mediated RBM4 exon skipping. Analysis of clinical data identified an association between the AURKA-YBX1/hnRNPK complex and poor prognosis in breast cancer. Blocking AURKA nuclear translocation with small molecule drugs partially reversed the oncogenic splicing of RBM4 and GOLGA4 in breast cancer cells. In summary, oncogenic AURKA executes its function on modulating breast cancer-related RNA splicing, and nuclear AURKA is distinguished as a hopeful target in the case of treating breast cancer.

ZBTB7A suppresses glioblastoma tumorigenesis through the transcriptional repression of EPB41L5

Glioblastoma multiforme (GBM), the most aggressive and malignant glioma, has a poor prognosis. Although patients with GBM are treated with surgery, chemotherapy, and radiation therapy, GBM is highly resistant to treatment, making it difficult and expensive to treat. In this study, we analyzed the Gene Expression Profiling Interactive Analysis dataset, the Cancer Genome Atlas dataset, and Gene Expression Omnibus array data. ZBTB7A (also called FBI1/POKEMON/LRF) was found to be highly expressed in low-grade glioma but significantly downregulated in patients with GBM. ZBTB7A is a transcription factor that plays an important role in many developmental stages, including cell proliferation. The activation of epithelial-mesenchymal transition (EMT) is a key process in cancer progression and metastasis. Erythrocyte membrane protein band 4.1 like 5 (EPB41L5) is an essential protein for EMT progression and metastasis in various types of cancer. We found that ZBTB7A depletion in U87 cells induced GBM progression and metastasis. Based on RNA sequencing data, ZBTB7A directly binds to the promoter of the EPB41L5 gene, reducing its expression and inhibiting GBM progression. We demonstrated that ZBTB7A dramatically inhibits GBM tumor growth through transcriptional repression of EPB41L5. Thus, both ZBTB7A and EPB41L5 may be potential biomarkers and novel therapeutic targets for GBM treatment. Overall, we discovered the role of a novel tumor suppressor that directly inhibits GBM progression (ZBTB7A) and identified EPB41L5 as a therapeutic target protein for patients with GBM.

The Chromatin-Oxygen Sensor Gene KDM5C Associates with Novel Hypoxia-Related Signatures in Glioblastoma Multiforme

Glioblastoma multiforme (GBM) is a fatal brain tumor without effective drug treatment. In this study, we highlight, for the first time, the contribution of chromatin remodeling gene Lysine (K)-specific demethylase 5C (KDM5C) in GBM via an extensive analysis of clinical, expression, and functional data, integrated with publicly available omic datasets. The expression analysis on GBM samples (N = 37) revealed two informative subtypes, namely KDM5C^High and KDM5C^Low, displaying higher/lower KDM5C levels compared to the controls. The former subtype displays a strong downregulation of brain-derived neurotrophic factor (BDNF)—a negative KDM5C target—and a robust overexpression of hypoxia-inducible transcription factor-1A (HIF1A) gene, a KDM5C modulator. Additionally, a significant co-expression among the prognostic markers HIF1A, Survivin, and p75 was observed. These results, corroborated by KDM5C overexpression and hypoxia-related functional assays in T98G cells, suggest a role for the HIF1A-KDM5C axis in the hypoxic response in this tumor. Interestingly, fluorescence-guided surgery on GBM sections further revealed higher KDM5C and HIF1A levels in the tumor rim niche compared to the adjacent tumor margin, indicating a regionally restricted hyperactivity of this regulatory axis. Analyzing the TCGA expression and methylation data, we found methylation changes between the subtypes in the genes, accounting for the hypoxia response, stem cell differentiation, and inflammation. High NANOG and IL6 levels highlight a distinctive stem cell-like and proinflammatory signature in the KDM5C^High subgroup and GBM niches. Taken together, our results indicate HIF1A-KDM5C as a new, relevant cancer axis in GBM, opening a new, interesting field of investigation based on KDM5C as a potential therapeutic target of the hypoxic microenvironment in GBM.

PTBPs: An immunomodulatory-related prognostic biomarker in pan-cancer

Background: The polypyrimidine tract-binding protein (PTBP) nuclear ribonucleoprotein family of proteins, including PTBP1, PTBP2 and PTBP3, regulate the process of cell proliferation, differentiation, apoptosis and carcinogenesis. PTBPs exhibit oncogenic effects in certain tumors. However, the role of PTBPs in pan-cancer remains unclear. Our study examined the clinical significance and mechanism of PTBPs in pan-cancer.

Methods: We compared the expression of PTBPs in paired and unpaired tissue samples from the Cancer Genome Atlas (TCGA) database. Univariate and multivariate Cox regression, Kaplan–Meier curves, and time-dependent receiver operating characteristic (ROC) curves were used to assess the prognostic significance of PTBPs in pan-cancer. The cBioPortal database also identified genomic abnormalities in PTBPs. TISIDB, TCGA, and Cellminer were used to investigate the relationship between PTBP expression and immune subtypes, immune checkpoint (ICP) genes, tumor mutational burden (TMB), microsatellite instability (MSI), tumor-infiltrating immune cells, and chemosensitivity. cBioPortal was used to search for PTBP co-expressing genes in pan-cancer, and GO and KEGG enrichment analyses were performed to search for PTBP-related signaling pathways.

Results:PTBPs were shown to be widely upregulated in human tumor tissues. PTBP1 showed good prognostic value in ACC, KIRP, and LGG; PTBP2 in ACC and KICH; and PTBP3 in ACC, LGG, and PAAD, with AUC >0.7. PTBPs were differentially expressed in tumor immune subtypes and had a strong correlation with tumor-infiltrating lymphocytes (TILs) in the tumor microenvironment (TME). In addition, PTBP expressions were related to ICP, TMB, and MSI, suggesting that these three PTBPs may be potential tumor immunotherapeutic targets and predict the efficacy of immunotherapy. Enrichment analysis of co-expressed genes of PTBPs showed that they may be involved in alternative splicing, cell cycle, cellular senescence, and protein modification.

Conclusion: PTBPs are involved in the malignant progression of tumors. PTBP1, PTBP2 and PTBP3 may be potential biomarkers for prognosis and immunotherapy in pan-cancer and may be novel immunotherapeutic targets.

A Novel Splice Variant of BCAS1 Inhibits β-Arrestin 2 to Promote the Proliferation and Migration of Glioblastoma Cells, and This Effect Was Blocked by Maackiain

Brain-enriched myelin-associated protein 1 (BCAS1) is frequently highly expressed in human cancer, but its detailed function is unclear. Here, we identified a novel splice variant of the BCAS1 gene in glioblastoma multiforme (GBM) named BCAS1-SV1. The expression of BCAS1-SV1 was weak in heathy brain cells but high in GBM cell lines. The overexpression of BCAS1-SV1 significantly increased the proliferation and migration of GBM cells, whereas the RNA-interference-mediated knockdown of BCAS1-SV1 reduced proliferation and migration. Moreover, using a yeast-two hybrid assay, immunoprecipitation, and immunofluorescence staining, we confirmed that β-arrestin 2 is an interaction partner of BCAS1-SV1 but not BCAS1. The downregulation of β-arrestin 2 directly enhanced the malignancy of GBM and abrogated the effects of BCAS1-SV1 on GBM cells. Finally, we used a yeast two-hybrid-based growth assay to identify that maackiain (MK) is a potential inhibitor of the interaction between BCAS1-SV1 and β-arrestin 2. MK treatment lessened the proliferation and migration of GBM cells and prolonged the lifespan of tumor-bearing mice in subcutaneous xenograft and intracranial U87-luc xenograft models. This study provides the first evidence that the gain-of-function BCAS1-SV1 splice variant promotes the development of GBM by suppressing the β-arrestin 2 pathway and opens up a new therapeutic perspective in GBM.

FOXD1 Is a Transcription Factor Important for Uveal Melanocyte Development and Associated with High-Risk Uveal Melanoma

Simple Summary

Despite successful treatment of primary uveal melanoma (UM), metastases still occur in approximately 50% of the patients. Unfortunately, little is known about the mechanism behind metastasized UM. By reanalyzing publicly available single-cell RNA sequencing data of embryonic zebrafish larvae and validating the results with UM data, we have identified five transcription regulators of interest: ELL2, KDM5B, REXO4, RBFOX2 and FOXD1. The most significant finding is FOXD1, which is nearly exclusively expressed in high-risk UM and is associated with poor survival. FOXD1 is a novel gene which could be involved in the metastatic capability of UM. Elucidating its function and role in metastatic UM could help to understand and develop treatment for UM.

Abstract

Uveal melanoma (UM) is a deadly ocular malignancy, originating from uveal melanocytes. Although much is known regarding prognostication in UM, the exact mechanism of metastasis is mostly unknown. Metastatic tumor cells are known to express a more stem-like RNA profile which is seen often in cell-specific embryonic development to induce tumor progression. Here, we identified novel transcription regulators by reanalyzing publicly available single cell RNA sequencing experiments. We identified five transcription regulators of interest: ELL2, KDM5B, REXO4, RBFOX2 and FOXD1. Our most significant finding is FOXD1, as this gene is nearly exclusively expressed in high-risk UM and its expression is associated with a poor prognosis. Even within the BAP1-mutated UM, the expression of FOXD1 is correlated with poor survival. FOXD1 is a novel factor which could potentially be involved in the metastatic capacity of high-risk UM. Elucidating the function of FOXD1 in UM could provide insight into the malignant transformation of uveal melanocytes, especially in high-risk UM.

LPS Administration Impacts Glial Immune Programs by Alternative Splicing

We performed transcriptome analysis in the hippocampus 24 h after lipopolysaccharide (LPS) administration. We observed glial-specific genes, comprised of two-thirds of all differentially expressed genes (DEGs). We found microglial DEGs that were the most numerous in LPS group. On the contrary, differential alternative splicing (DAS) analysis revealed the most numerous DAS events in astrocytes. Besides, we observed distinct major isoform switching in the Ptbp1 gene, with skipping of exon 8 in LPS group. Ptbp1 usually considered a pluripotency sustaining agent in brain embryonic development, according to the previous studies. Analyzing the splicing tune-up upon LPS exposure, we came to a supposition that the short Ptbp1 isoform de-represses immune-specific response by Ptbp1 adjusted splicing architecture. Additionally, the Ptbp3 (NOD1) immune-specific splicing factor has apparently been de-repressed by the Ptbp1 short isoform in glial cells. Notably, both the Ptbp1 and Ptbp3 genes express primarily in microglial/endothelial brain cells. We also report immune-related genes, altering their major isoforms upon LPS exposure. The results revealed immune modulating role of alternative splicing in brain.