Comprehensive Analysis of Alternative Splicing Across Tumors from 8,705 Patients

KIAA1429 regulates alternative splicing events of cancer-related genes in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) remains one of the most fatal malignancies with high morbidity and mortality rates in the world, whose molecular pathogenesis is incompletely understood. As an RNA-binding protein participating in the processing and modification of RNA, KIAA1429 has been proved to be implicated in the pathogenesis of multiple cancers. However, how KIAA1429 functions in alternative splicing is not fully reported. In the current study, multi-omics sequencing data were used to analyze and decipher the molecular functions and the underlying mechanisms of KIAA1429 in HCC samples. RNA sequencing data (RNA-seq) analysis demonstrated that in HCCLM3 cells, alternative splicing (AS) profiles were mediated by KIAA1429. Regulated AS genes (RASGs) by KIAA1429 were enriched in cell cycle and apoptosis-associated pathways. Furthermore, by integrating the RNA immunoprecipitation and sequencing data (RIP-seq) of KIAA1429, we found that KIAA1429-bound transcripts were highly overlapping with RASGs, indicating that KIAA1429 could globally regulate the alternative splicing perhaps by binding to their transcripts in HCCLM3 cells. The overlapping RASGs were also clustered in cell cycle and apoptosis-associated pathways. In particular, we validated the regulated AS events of three genes using clinical specimens from HCC patients, including the exon 6 of BPTF gene and a marker gene of HCC. In summary, our results shed light on the regulatory functions of KIAA1429 in the splicing process of pre-mRNA and provide theoretical basis for the targeted therapy of HCC.

Multi-omics approach to identifying isoform variants as therapeutic targets in cancer patients

Cancer-specific alternatively spliced events (ASE) play a role in cancer pathogenesis and can be targeted by immunotherapy, oligonucleotide therapy, and small molecule inhibition. However, identifying actionable ASE targets remains challenging due to the uncertainty of its protein product, structure impact, and proteoform (protein isoform) function. Here we argue that an integrated multi-omics profiling strategy can overcome these challenges, allowing us to mine this untapped source of targets for therapeutic development. In this review, we will provide an overview of current multi-omics strategies in characterizing ASEs by utilizing the transcriptome, proteome, and state-of-art algorithms for protein structure prediction. We will discuss limitations and knowledge gaps associated with each technology and informatics analytics. Finally, we will discuss future directions that will enable the full integration of multi-omics data for ASE target discovery.

Retained introns in long RNA-seq reads are not reliably detected in sample-matched short reads

BACKGROUND

There is growing interest in retained introns in a variety of disease contexts including cancer and aging. Many software tools have been developed to detect retained introns from short RNA-seq reads, but reliable detection is complicated by overlapping genes and transcripts as well as the presence of unprocessed or partially processed RNAs.

RESULTS

We compared introns detected by 8 tools using short RNA-seq reads with introns observed in long RNA-seq reads from the same biological specimens. We found significant disagreement among tools (Fleiss' [Formula: see text]) such that 47.7% of all detected intron retentions were not called by more than one tool. We also observed poor performance of all tools, with none achieving an F1-score greater than 0.26, and qualitatively different behaviors between general-purpose alternative splicing detection tools and tools confined to retained intron detection.

CONCLUSIONS

Short-read tools detect intron retention with poor recall and precision, calling into question the completeness and validity of a large percentage of putatively retained introns called by commonly used methods.

Systematic characterization of cancer transcriptome at transcript resolution

Transcribed RNAs undergo various regulation and modification to become functional transcripts. Notably, cancer transcriptome has not been fully characterized at transcript resolution. Herein, we carry out a reference-based transcript assembly across >1000 cancer cell lines. We identify 498,255 transcripts, approximately half of which are unannotated. Unannotated transcripts are closely associated with cancer-related hallmarks and show clinical significance. We build a high-confidence RNA binding protein (RBP)-transcript regulatory network, wherein most RBPs tend to regulate transcripts involved in cell proliferation. We identify numerous transcripts that are highly associated with anti-cancer drug sensitivity. Furthermore, we establish RBP-transcript-drug axes, wherein PTBP1 is experimentally validated to affect the sensitivity to decitabine by regulating KIAA1522-a6 transcript. Finally, we establish a user-friendly data portal to serve as a valuable resource for understanding cancer transcriptome diversity and its potential clinical utility at transcript level. Our study substantially extends cancer RNA repository and will facilitate anti-cancer drug discovery.

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.

ESRP1-regulated isoform switching of LRRFIP2 determines metastasis of gastric cancer

Although accumulating evidence indicates that alternative splicing is aberrantly altered in many cancers, the functional mechanism remains to be elucidated. Here, we show that epithelial and mesenchymal isoform switches of leucine-rich repeat Fli-I-interacting protein 2 (LRRFIP2) regulated by epithelial splicing regulatory protein 1 (ESRP1) correlate with metastatic potential of gastric cancer cells. We found that expression of the splicing variants of LRRFIP2 was closely correlated with that of ESRP1. Surprisingly, ectopic expression of the mesenchymal isoform of LRRFIP2 (variant 3) dramatically increased liver metastasis of gastric cancer cells, whereas deletion of exon 7 of LRRFIP2 by the CRISPR/Cas9 system caused an isoform switch, leading to marked suppression of liver metastasis. Mechanistically, the epithelial LRRFIP2 isoform (variant 2) inhibited the oncogenic function of coactivator-associated arginine methyltransferase 1 (CARM1) through interaction. Taken together, our data reveals a mechanism of LRRFIP2 isoform switches in gastric cancer with important implication for cancer metastasis.

Analyses of Transcriptomics Cell Signalling for Pre-Screening Applications in the Integrated Approach for Testing and Assessment of Non-Genotoxic Carcinogens

With recent rapid advancement of methodological tools, mechanistic understanding of biological processes leading to carcinogenesis is expanding. New approach methodologies such as transcriptomics can inform on non-genotoxic mechanisms of chemical carcinogens and can be developed for regulatory applications. The Organisation for the Economic Cooperation and Development (OECD) expert group developing an Integrated Approach to the Testing and Assessment (IATA) of Non-Genotoxic Carcinogens (NGTxC) is reviewing the possible assays to be integrated therein. In this context, we review the application of transcriptomics approaches suitable for pre-screening gene expression changes associated with phenotypic alterations that underlie the carcinogenic processes for subsequent prioritisation of downstream test methods appropriate to specific key events of non-genotoxic carcinogenesis. Using case studies, we evaluate the potential of gene expression analyses especially in relation to breast cancer, to identify the most relevant approaches that could be utilised as (pre-) screening tools, for example Gene Set Enrichment Analysis (GSEA). We also consider how to address the challenges to integrate gene panels and transcriptomic assays into the IATA, highlighting the pivotal omics markers identified for assay measurement in the IATA key events of inflammation, immune response, mitogenic signalling and cell injury.

Boosting Antitumor Immunity with an Expanded Neoepitope Landscape

Immune-checkpoint blockade therapy has been successfully applied to many cancers, particularly tumors that harbor a high mutational burden and consequently express a high abundance of neoantigens. However, novel approaches are needed to improve the efficacy of immunotherapy for treating tumors that lack a high load of classic genetically derived neoantigens. Recent discoveries of broad classes of nongenetically encoded and inducible neoepitopes open up new avenues for therapeutic development to enhance sensitivity to immunotherapies. In this review, we discuss recent work on neoantigen discovery, with an emphasis on novel classes of noncanonical neoepitopes.

Roles of Krüppel-like factor 6 splice variant 1 in the development, diagnosis, and possible treatment strategies for non-small cell lung cancer

Krüppel-like factor 6 (KLF6) is a nuclear transcriptional regulator found in mammalian tissue that has been identified as a tumor suppressor gene in several malignancies. As a result of loss of heterozygosity, DNA methylation, and alternative splicing, it is frequently inactivated in various malignancies. Krüppel-like factor 6 splice variant 1 (KLF6-SV1), Krüppel-like factor 6 splice variant 2, and Krüppel-like factor 6 splice variant 3 alternatively spliced isoforms that emerge from a single nucleotide polymorphism in the KLF6 gene. KLF6-SV1 is generally upregulated in multiple cancers, and its biological function is well understood. Overexpression of KLF6-SV1 inhibits the KLF6 gene function while promoting tumor progression, which is associated with a poor prognosis in patients with various malignancies. We reviewed the progress of KLF6-SV1 research in NSCLC over the last several years to understand the molecular mechanisms of tumorigenesis, tumor development, and therapy resistance. Finally, this review emphasizes the therapeutic potential of small interfering RNA targeted silencing of KLF6-SV1 as a novel strategy for managing chemotherapy resistance in NSCLC patients.

Landscape of exitrons in gastric cancer

BACKGROUND

Exitron is a new type of non-canonical alternative splicing. Accumulating evidence implies exitron may have pathological function and contribute to another source of anti-tumor immunogenicity in various cancers. Its role in gastric cancer remains poorly understood. Large-scale, multi-omics analysis could comprehensively characterize the landscape of exitrons in gastric cancer, reveal undiscovered mechanism and hopefully identify molecular biomarkers for predicting immunotherapy response.

METHODS

We collected datasets from five studies for analysis. RNA sequencing was used for exitron identification. Somatic mutations were detected by whole exome sequencing. Neopeptides were confirmed by proteome mass spectrometry.

FINDINGS

42174 gastric cancer-specific exitrons (GCSEs) were identified in 632 patients. GCSEs were clinically relevant to gender, age, Lauren type, tumor stage and prognosis. Tissue specificity test and pathogenic exitron prediction revealed their unique functional impact. GCSEs were mutually exclusive with mutations and demonstrated both unique and complementary function against TP53 mutation in gastric cancer. We further established splicing regulatory network to reveal upstream regulation of exitron splicing. We also evaluated the immunogenicity and diagnostic potential of GCSEs. Evidence of GCSEs-derived neopeptide expression was validated by whole proteome mass spectrometry. PD-1 and Siglecs were significantly increased in high neoantigen load patients. But exitron-related biomarkers failed to predict immunotherapy response, possibly due to small sample size and insufficient sequencing depth.

INTERPRETATION

The present study provided a comprehensive multidimensional landscape of gastric cancer exitrons and underscores insights into underexplored mechanism in gastric cancer pathology.

FUNDING

The Guangdong Provincial Key Laboratory of Precision Medicine for Gastroinstestinal Cancer (2020B121201004), the Guangdong Provincial Major Talents Project (No. 2019JC05Y361) and National Natural Science Foundation of China (grant number:82172960 and 81872013).

FOXA1 regulates alternative splicing in prostate cancer

Dysregulation of alternative splicing in prostate cancer is linked to transcriptional programs activated by AR, ERG, FOXA1, and MYC. Here, we show that FOXA1 functions as the primary orchestrator of alternative splicing dysregulation across 500 primary and metastatic prostate cancer transcriptomes. We demonstrate that FOXA1 binds to the regulatory regions of splicing-related genes, including HNRNPK and SRSF1. By controlling trans-acting factor expression, FOXA1 exploits an "exon definition" mechanism calibrating alternative splicing toward dominant isoform production. This regulation especially impacts splicing factors themselves and leads to a reduction of nonsense-mediated decay (NMD)-targeted isoforms. Inclusion of the NMD-determinant FLNA exon 30 by FOXA1-controlled oncogene SRSF1 promotes cell growth in vitro and predicts disease recurrence. Overall, we report a role for FOXA1 in rewiring the alternative splicing landscape in prostate cancer through a cascade of events from chromatin access, to splicing factor regulation, and, finally, to alternative splicing of exons influencing patient survival.

Proteotranscriptomic Integration analyses reveals new mechanistic insights regarding Bombyx mori fluorosis

The commercial value of silkworms has been widely explored and the effects of fluoride exposure on silkworms' breeding and silk production cannot be ignored. Bombyx mori is a commonly used model to explore the mechanisms of fluorosis. In the present study, we analyzed the differences in physiological and biochemical indicators after exposing larva to NaF, then evaluated differential genes and proteins. Compared to control, larvae exposed to 600 mg L^-1 NaF presented decreased bodyweight, damaged midgut tissue, and were accompanied by oxidative stress. The RNA-seq showed 1493 differentially expressed genes (574 upregulated and 919 downregulated). Meanwhile, the TMT detected 189 differentially expressed proteins (133 upregulated and 56 downregulated). The integrative analysis led to 4 upregulated and 9 downregulated genes and proteins. Finally, we hypothesized that fluoride exposure might affect the intestinal digestion of silkworms, inhibit the gene expression of detoxification enzymes and stimulate cellular immune responses. Our current findings provided new insights into insect fluorosis.

DDX17 modulates the expression and alternative splicing of genes involved in apoptosis and proliferation in lung adenocarcinoma cells

Background

The DEAD-box RNA-binding protein (RBP) DDX17 has been found to be involved in the tumorigenesis of many types of cancers. However, the role of DDX17 in lung adenocarcinoma (LUAD) remains unclear.

Methods

We silenced DDX17 expression in A549 LUAD cells by small interfering RNA (siRNA). Cell proliferation and apoptosis assays were performed to explore the functions of DDX17. Knockdown of DDX17 by siRNA significantly inhibited proliferation and induced apoptosis in A549 cells. We used high-throughput RNA sequencing (RNA-seq) to identify differentially expressed genes (DEGs) and alternative splicing (AS) events in DDX17 knockdown LUAD cells.

Results

DDX17 knockdown increased the expression levels of proapoptotic genes and decreased those of proproliferative genes. Moreover, the DDX17-regulated AS events in A549 cells revealed by computational analysis using ABLas software were strongly validated by quantitative reverse transcription-polymerase chain reaction (RT-qPCR) and were also validated by analysis of The Cancer Genome Atlas (TCGA)-LUAD dataset. These findings suggest that DDX17 may function as an oncogene by regulating both the expression and AS of proliferation- and apoptosis-associated genes in LUAD cells. Our findings may offer new insights into understanding the molecular mechanisms of LUAD and provide a new therapeutic direction for LUAD.

Integration of TE Induces Cancer Specific Alternative Splicing Events

Alternative splicing of messenger RNA (mRNA) precursors contributes to genetic diversity by generating structurally and functionally distinct transcripts. In a disease state, alternative splicing promotes incidence and development of several cancer types through regulation of cancer-related biological processes. Transposable elements (TEs), having the genetic ability to jump to other regions of the genome, can bring about alternative splicing events in cancer. TEs can integrate into the genome, mostly in the intronic regions, and induce cancer-specific alternative splicing by adjusting various mechanisms, such as exonization, providing splicing donor/acceptor sites, alternative regulatory sequences or stop codons, and driving exon disruption or epigenetic regulation. Moreover, TEs can produce microRNAs (miRNAs) that control the proportion of transcripts by repressing translation or stimulating the degradation of transcripts at the post-transcriptional level. Notably, TE insertion creates a cancer-friendly environment by controlling the overall process of gene expression before and after transcription in cancer cells. This review emphasizes the correlative interaction between alternative splicing by TE integration and cancer-associated biological processes, suggesting a macroscopic mechanism controlling alternative splicing by TE insertion in cancer.

Antigens Expressed by Breast Cancer Cells Undergoing EMT Stimulate Cytotoxic CD8+ T Cell Immunity

Antigenic differences formed by alterations in gene expression and alternative splicing are predicted in breast cancer cells undergoing epithelial to mesenchymal transition (EMT) and the reverse plasticity known as MET. How these antigenic differences impact immune interactions and the degree to which they can be exploited to enhance immune responses against mesenchymal cells is not fully understood. We utilized a master microRNA regulator of EMT to alter mesenchymal-like EO771 mammary carcinoma cells to a more epithelial phenotype. A computational approach was used to identify neoantigens derived from the resultant differentially expressed somatic variants (SNV) and alternative splicing events (neojunctions). Using whole cell vaccines and peptide-based vaccines, we find superior cytotoxicity against the more-epithelial cells and explore the potential of neojunction-derived antigens to elicit T cell responses through experiments designed to validate the computationally predicted neoantigens. Overall, results identify EMT-associated splicing factors common to both mouse and human breast cancer cells as well as immunogenic SNV- and neojunction-derived neoantigens in mammary carcinoma cells.

Sources of Cancer Neoantigens beyond Single-Nucleotide Variants

The success of checkpoint blockade therapy against cancer has unequivocally shown that cancer cells can be effectively recognized by the immune system and eliminated. However, the identity of the cancer antigens that elicit protective immunity remains to be fully explored. Over the last decade, most of the focus has been on somatic mutations derived from non-synonymous single-nucleotide variants (SNVs) and small insertion/deletion mutations (indels) that accumulate during cancer progression. Mutated peptides can be presented on MHC molecules and give rise to novel antigens or neoantigens, which have been shown to induce potent anti-tumor immune responses. A limitation with SNV-neoantigens is that they are patient-specific and their accurate prediction is critical for the development of effective immunotherapies. In addition, cancer types with low mutation burden may not display sufficient high-quality [SNV/small indels] neoantigens to alone stimulate effective T cell responses. Accumulating evidence suggests the existence of alternative sources of cancer neoantigens, such as gene fusions, alternative splicing variants, post-translational modifications, and transposable elements, which may be attractive novel targets for immunotherapy. In this review, we describe the recent technological advances in the identification of these novel sources of neoantigens, the experimental evidence for their presentation on MHC molecules and their immunogenicity, as well as the current clinical development stage of immunotherapy targeting these neoantigens.

Dysregulation of splicing variants and spliceosome components in breast cancer

The dysregulation of the splicing process has emerged as a novel hallmark of metabolic and tumor pathologies. In breast cancer (BCa), which represents the most diagnosed cancer type among women worldwide, the generation and/or dysregulation of several oncogenic splicing variants have been described. This is the case of the splicing variants of HER2, ER, BRCA1, or the recently identified by our group, In1-ghrelin and SST5TMD4, which exhibit oncogenic roles, increasing the malignancy, poor prognosis, and resistance to treatment of BCa. This altered expression of oncogenic splicing variants has been closely linked with the dysregulation of the elements belonging to the macromolecular machinery that controls the splicing process (spliceosome components and the associated splicing factors). In this review, we compile the current knowledge demonstrating the altered expression of splicing variants and spliceosomal components in BCa, showing the existence of a growing body of evidence supporting the close implication of the alteration in the splicing process in mammary tumorigenesis.

Quantitative immunopeptidomics reveals a tumor stroma-specific target for T cell therapy

T cell receptor (TCR)-based immunotherapy has emerged as a promising therapeutic approach for the treatment of patients with solid cancers. Identifying peptide-human leukocyte antigen (pHLA) complexes highly presented on tumors and rarely expressed on healthy tissue in combination with high-affinity TCRs that when introduced into T cells can redirect T cells to eliminate tumor but not healthy tissue is a key requirement for safe and efficacious TCR-based therapies. To discover promising shared tumor antigens that could be targeted via TCR-based adoptive T cell therapy, we employed population-scale immunopeptidomics using quantitative mass spectrometry across ~1500 tumor and normal tissue samples. We identified an HLA-A*02:01-restricted pan-cancer epitope within the collagen type VI α-3 (COL6A3) gene that is highly presented on tumor stroma across multiple solid cancers due to a tumor-specific alternative splicing event that rarely occurs outside the tumor microenvironment. T cells expressing natural COL6A3-specific TCRs demonstrated only modest activity against cells presenting high copy numbers of COL6A3 pHLAs. One of these TCRs was affinity-enhanced, enabling transduced T cells to specifically eliminate tumors in vivo that expressed similar copy numbers of pHLAs as primary tumor specimens. The enhanced TCR variants exhibited a favorable safety profile with no detectable off-target reactivity, paving the way to initiate clinical trials using COL6A3-specific TCRs to target an array of solid tumors.

Harnessing gene fusion-derived neoantigens for 'cold' breast and prostate tumor immunotherapy

Breast and prostate cancers are generally considered immunologically 'cold' tumors due to multiple mechanisms rendering them unresponsive to immune checkpoint blockade therapies. With little success in garnering positive outcomes in modern immunotherapeutic clinical trials, it is prudent to re-examine the role of immunogenic neoantigens in these cold tumors. Gene fusions are driver mutations in hormone-driven cancers that can result in alternative mutation-specific neoantigens to promote immunotherapy sensitivity. This review focuses on 1) gene fusion formation mechanisms in neoantigen generation; 2) gene fusion neoantigens in cancer immunotherapeutic strategies and associated clinical trials; and 3) challenges and opportunities in computational and liquid biopsy technologies. This review is anticipated to initiate further research into gene fusion neoantigens of cold tumors for further experimental validation.

SF3B1 facilitates HIF1-signaling and promotes malignancy in pancreatic cancer

Mutations in the splicing factor SF3B1 are frequently occurring in various cancers and drive tumor progression through the activation of cryptic splice sites in multiple genes. Recent studies also demonstrate a positive correlation between the expression levels of wild-type SF3B1 and tumor malignancy. Here, we demonstrate that SF3B1 is a hypoxia-inducible factor (HIF)-1 target gene that positively regulates HIF1 pathway activity. By physically interacting with HIF1α, SF3B1 facilitates binding of the HIF1 complex to hypoxia response elements (HREs) to activate target gene expression. To further validate the relevance of this mechanism for tumor progression, we show that a reduction in SF3B1 levels via monoallelic deletion of Sf3b1 impedes tumor formation and progression via impaired HIF signaling in a mouse model for pancreatic cancer. Our work uncovers an essential role of SF3B1 in HIF1 signaling, thereby providing a potential explanation for the link between high SF3B1 expression and aggressiveness of solid tumors.

Multi-omics analysis reveals RNA splicing alterations and their biological and clinical implications in lung adenocarcinoma

Alternative RNA splicing is one of the most important mechanisms of posttranscriptional gene regulation, which contributes to protein diversity in eukaryotes. It is well known that RNA splicing dysregulation is a critical mechanism in tumor pathogenesis and the rationale for the promising splice-switching therapeutics for cancer treatment. Although we have a comprehensive understanding of DNA mutations, abnormal gene expression profiles, epigenomics, and proteomics in lung adenocarcinoma (LUAD), little is known about its aberrant alternative splicing profiles. Here, based on the multi-omics data generated from over 1000 samples, we systematically studied the RNA splicing alterations in LUAD and revealed their biological and clinical implications. We identified 3688 aberrant alternative splicing events (AASEs) in LUAD, most of which were alternative promoter and exon skip. The specific regulatory roles of RNA binding proteins, somatic mutations, and DNA methylations on AASEs were comprehensively interrogated. We dissected the functional implications of AASEs and concluded that AASEs mainly affected biological processes related to tumor proliferation and metastasis. We also found that one subtype of LUAD with a particular AASEs pattern was immunogenic and had a better prognosis and response rate to immunotherapy. These findings revealed novel events related to tumorigenesis and tumor immune microenvironment and laid the foundation for the development of splice-switching therapies for LUAD.

Construction of a breast cancer prognosis model based on alternative splicing and immune infiltration

BACKGROUND

Breast cancer (BC) is the most common malignancy among women in the world. Alternative splicing (AS) is an important mechanism for regulating gene expression and producing proteome diversity, which is closely related to tumorigenesis. Understanding the role of AS in BC may be helpful to reveal new therapeutic targets for clinical interventions.

METHODS

RNA-seq, clinical and AS data of TCGA-BRCA were downloaded from TCGA and TCGA SpliceSeq databases. AS events associated with prognosis were filtered by univariate Cox regression. The AS risk model of BC was built by Lasso regression, random forest and multivariate Cox regression. The accuracy of the AS risk model and clinicopathological factors were evaluated by time-dependent receiver operating characteristic (ROC) curves. The significant factors were used to construct the nomogram model. Tumor microenvironment analysis, immune infiltration and immune checkpoint analysis were performed to show the differences between the high and low AS risk groups. The expression differences of genes of AS events constituting the risk model in tumor tissues and normal tissues were analyzed, the genes with significant differences were screened, and their relationship with prognosis, tumor microenvironment, immune infiltration and immune checkpoint were analyzed. Finally, Pearson correlation analysis was used to calculate the correlation coefficient between splicing factors (SF) and prognostic AS events in TCGA-BRCA. The results were imported into Cytoscape, and the associated network was constructed.

RESULTS

A total of 21,232 genes had 45,421 AS events occurring in TCGA-BRCA, while 1604 AS events were found to be significantly correlated with survival. The BRCA risk model consisted of 5 AS events, (TTC39C|44853|AT*- 2.67) + (HSPBP1|52052|AP*- 4.28) + (MAZ|35942|ES*2.34) + (ANK3|11845|AP*1.18) + (ZC3HAV1|81940|AT*1.59), which were confirmed to be valuable for predicting BRCA prognosis to a certain degree, including ROC curve, survival analysis, tumor microenvironment analysis, immune infiltration and immune checkpoint analysis. Based on this, we constructed a nomogram prediction model composed of clinicopathological features and the AS risk signature. Furthermore, we found that MAZ was a core gene indicating the connection of tumor prognosis and AS events. Ultimately, a network of SF-AS regulation was established to reveal the relationship between them.

CONCLUSIONS

We constructed a nomogram model combined with clinicopathological features and AS risk score to predict the prognosis of BC. The detailed analysis of tumor microenvironment and immune infiltration in the AS risk model may further reveal the potential mechanisms of BC recurrence and development.

XBP1 variant 1 promotes mitosis of cancer cells involving upregulation of the polyglutamylase TTLL6

XBP1 variant 1 (Xv1) is the most abundant XBP1 variant and is highly enriched across cancer types but nearly none in normal tissues. Its expression is associated with poor patients' survival and is specifically required for survival of malignant cells, but the underlying mechanism is not known. Here we report that Xv1 upregulates the polyglutamylase tubulin tyrosine ligase-like 6 (TTLL6) and promotes mitosis of cancer cells. Like the canonical XBP1, Xv1 mRNA undergoes unconventional splicing by IRE1α under endoplasmic reticulum stress, but it is also constitutively spliced by IRE1β. The spliced Xv1 mRNA encodes the active form of Xv1 protein (Xv1s). RNA sequencing in HeLa cells revealed that Xv1s overexpression regulates expression of genes that are not involved in the canonical unfolded protein response, including TTLL6 as a highly upregulated gene. Gel shift assay and chromatin immunoprecipitation revealed that Xv1s bind to the TTLL6 promoter region. Knockdown of TTLL6 caused death of cancer cells but not benign and normal cells, similar to the effects of knocking down Xv1. Moreover, overexpression of TTLL6 partially rescued BT474 cells from apoptosis induced by either TTLL6 or Xv1 knockdown, supporting TTLL6 as an essential downstream effector of Xv1 in regulating cancer cell survival. TTLL6 is localized in the mitotic spindle of cancer cells. Xv1 or TTLL6 knockdown resulted in decreased spindle polyglutamylation and interpolar spindle, as well as congression failure, mitotic arrest and cell death. These findings suggest that Xv1 is essential for cancer cell mitosis, which is mediated, at least in part, by increasing TTLL6 expression.

Identification of personalized neoantigen-based vaccines and immune subtype characteristic analysis of glioblastoma based on abnormal alternative splicing

The development of personalized neoantigen-based vaccines in cancer immunotherapy has shown promise. In this study, a large-scale bioinformatics analysis was performed to identify potential GBM-associated neoantigens based on abnormal alternative splicing, and then screen suitable patients for vaccination. Gene expression profiles and clinical information were collected from TCGA. We filtered the percent-spliced-in (PSI) spectrum of alternative splicing events in the dataset to identify abnormal alternative splicing events. MAF package was used to identify and analyse tumour mutation burden (TMB) in cancer samples. Tumour Immune Estimation Resource (TIMER) was used to calculate and visualize the infiltration of antigen presenting cells (APCs). In addition, consistent clustering algorithm utilized to identify immune subtypes of GBM. Five potential tumour neoantigens (LRP1, TCF12, DERL3, WIPI2, and TSHZ3) were identified in GBM by selecting genes both with abnormal alternative splicing (upregulated) and gene frameshift mutations, in which LRP1 was significantly associated with APCs. According to the expressions of five potential tumour neoantigens, 160 patients with GBM were divided into three immune subtypes. Patients in cluster3 exhibited good prognoses. Furthermore, the characteristics, including TMB, abnormal alternative splicing events, immune activity, immune cells proportion, and association with tumour biomarkers, were unique in each immune subtypes. The characteristics of cluster3 illustrated that cluster3 participants were more suitable candidates for vaccination. LRP1 was identified as a potential neoantigen for immunotherapy against GBM, and patients in cluster3 were more suitable for vaccination. Our findings provide important guidance for the development of novel neoantigens and therapeutic targets in patients with GBM.

Detection of orthologous exons and isoforms using EGIO

MOTIVATION

Alternative splicing is an important mechanism to generate transcriptomic and phenotypic diversity. Existing methods have limited power to detect orthologous isoforms.

RESULTS

We develop a new method, EGIO, to detect orthologous exons and orthologous isoforms from two species. EGIO uses unique exonic regions to construct exon groups, in which process dynamic programming strategy is used to do exon alignment. EGIO could cover all the coding exons within orthologous genes. A comparison between EGIO and ExTraMapper shows that EGIO could detect more orthologous isoforms with conserved sequence and exon structures. We apply EGIO to compare human and chimpanzee protein-coding isoforms expressed in the frontal cortex and identify 6912 genes that express human unique isoforms. Unexpectedly, more human unique isoforms are detected than those conserved between humans and chimpanzees.

AVAILABILITY AND IMPLEMENTATION

Source code and test data of EGIO are available at https://github.com/wu-lab-egio/EGIO.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Isolation, Characterization, and Autophagy Function of BECN1-Splicing Isoforms in Cancer Cells

Alternative splicing allows the synthesis of different protein variants starting from a single gene. Human Beclin 1 (BECN1) is a key autophagy regulator that acts as haploinsufficient tumor suppressor since its decreased expression correlates with tumorigenesis and poor prognosis in cancer patients. Recent studies show that BECN1 mRNA undergoes alternative splicing. Here, we report on the isolation and molecular and functional characterization of three BECN1 transcript variants (named BECN1-α, -β and -γ) in human cancer cells. In ovarian cancer NIHOVCAR3, these splicing variants were found along with the canonical wild-type. BECN1-α lacks 143 nucleotides at its C-terminus and corresponds to a variant previously described. BECN1-β and -γ lack the BCL2 homology 3 domain and other regions at their C-termini. Following overexpression in breast cancer cells MDA-MB231, we found that BECN1-α stimulates autophagy. Specifically, BECN1-α binds to Parkin and stimulates mitophagy. On the contrary, BECN1-β reduces autophagy with a dominant negative effect over the endogenous wild-type isoform. BECN1-γ maintains its ability to interact with the vacuolar protein sorting 34 and only has a slight effect on autophagy. It is possible that cancer cells utilize the alternative splicing of BECN1 for modulating autophagy and mitophagy in response to environmental stresses.

Nano-omics: nanotechnology-based multidimensional harvesting of the blood-circulating cancerome

Over the past decade, the development of 'simple' blood tests that enable cancer screening, diagnosis or monitoring and facilitate the design of personalized therapies without the need for invasive tumour biopsy sampling has been a core ambition in cancer research. Data emerging from ongoing biomarker development efforts indicate that multiple markers, used individually or as part of a multimodal panel, are required to enhance the sensitivity and specificity of assays for early stage cancer detection. The discovery of cancer-associated molecular alterations that are reflected in blood at multiple dimensions (genome, epigenome, transcriptome, proteome and metabolome) and integration of the resultant multi-omics data have the potential to uncover novel biomarkers as well as to further elucidate the underlying molecular pathways. Herein, we review key advances in multi-omics liquid biopsy approaches and introduce the 'nano-omics' paradigm: the development and utilization of nanotechnology tools for the enrichment and subsequent omics analysis of the blood-circulating cancerome.

Splicing dysregulation in human hematologic malignancies: beyond splicing mutations

Splicing is a fundamental process in pre-mRNA maturation. Whereas alternative splicing (AS) enriches the diversity of the proteome, its aberrant regulation can drive oncogenesis. So far, most attention has been given to spliceosome mutations (SMs) in the context of splicing dysregulation in hematologic diseases. However, in recent years, post-translational modifications (PTMs) and transcriptional alterations of splicing factors (SFs), just as epigenetic signatures, have all been shown to contribute to global splicing dysregulation as well. In addition, the contribution of aberrant splicing to the neoantigen repertoire of cancers has been recognized. With the pressing need for novel therapeutics to combat blood cancers, this article provides an overview of emerging mechanisms that contribute to aberrant splicing, as well as their clinical potential.

Analysis of Alternative mRNA Splicing in Vemurafenib-Resistant Melanoma Cells

Alternative mRNA splicing is common in cancers. In BRAF V600E-mutated malignant melanoma, a frequent mechanism of acquired resistance to BRAF inhibitors involves alternative splicing (AS) of BRAF. The resulting shortened BRAF protein constitutively dimerizes and conveys drug resistance. Here, we have analysed AS in SK-MEL-239 melanoma cells and a BRAF inhibitor (vemurafenib)-resistant derivative that expresses an AS, shortened BRAF V600E transcript. Transcriptome analysis showed differential expression of spliceosome components between the two cell lines. As there is no consensus approach to analysing AS events, we used and compared four common AS softwares based on different principles, DEXSeq, rMATS, ASpli, and LeafCutter. Two of them correctly identified the BRAF V600E AS in the vemurafenib-resistant cells. Only 12 AS events were identified by all four softwares. Testing the AS predictions experimentally showed that these overlapping predictions are highly accurate. Interestingly, they identified AS caused alterations in the expression of melanin synthesis and cell migration genes in the vemurafenib-resistant cells. This analysis shows that combining different AS analysis approaches produces reliable results and meaningful, biologically testable hypotheses.

A Regulatory Axis between Epithelial Splicing Regulatory Proteins and Estrogen Receptor α Modulates the Alternative Transcriptome of Luminal Breast Cancer

Epithelial splicing regulatory proteins 1 and 2 (ESRP1/2) control the splicing pattern during epithelial to mesenchymal transition (EMT) in a physiological context and in cancer, including breast cancer (BC). Here, we report that ESRP1, but not ESRP2, is overexpressed in luminal BCs of patients with poor prognosis and correlates with estrogen receptor α (ERα) levels. Analysis of ERα genome-binding profiles in cell lines and primary breast tumors showed its binding in the proximity of ESRP1 and ESRP2 genes, whose expression is strongly decreased by ERα silencing in hormone-deprived conditions. The combined knock-down of ESRP1/2 in MCF-7 cells followed by RNA-Seq, revealed the dysregulation of 754 genes, with a widespread alteration of alternative splicing events (ASEs) of genes involved in cell signaling, metabolism, cell growth, and EMT. Functional network analysis of ASEs correlated with ESRP1/2 expression in ERα+ BCs showed RAC1 as the hub node in the protein-protein interactions altered by ESRP1/2 silencing. The comparison of ERα- and ESRP-modulated ASEs revealed 63 commonly regulated events, including 27 detected in primary BCs and endocrine-resistant cell lines. Our data support a functional implication of the ERα-ESRP1/2 axis in the onset and progression of BC by controlling the splicing patterns of related genes.

RNA Instant Quality Check: Alignment-Free RNA-Degradation Detection

With the constant increase of large-scale genomic data projects, automated and high-throughput quality assessment becomes a crucial component of any analysis. Whereas small projects often have a more homogeneous design and a manageable structure allowing for a manual per-sample analysis of quality, large-scale studies tend to be much more heterogeneous and complex. Many quality metrics have been developed to assess the quality of an individual sample on the raw read level. Degradation effects are typically assessed based on the RNA integrity (RIN) score, or on postalignment data. In this study, we show that single commonly used quality criteria such as the RIN score alone are not sufficient to ensure RNA sample quality. We developed a new approach and provide an efficient tool that estimates RNA sample degradation by computing the 5'/3' bias based on all genes in an alignment-free manner. That enables degradation assessment right after data generation and not during the analysis procedure allowing for early intervention in the sample handling process. Our analysis shows that this strategy is fast, robust to annotation and differences in library size, and provides complementary quality information to RIN scores enabling the accurate identification of degraded samples.

The role of neoantigens in tumor immunotherapy

Tumor neoantigens are aberrant polypeptides produced by tumor cells as a result of genomic mutations. They are also tumor-specific antigens (TSA). Neoantigens are more immunogenic than tumor-related antigens and do not induce autoimmunity. Based on the rapid development of bioinformatics and the continuous update of sequencing technology, cancer immunotherapy with tumor neoantigens has made promising breakthroughs and progress. In this review, the generation, prediction, and identification of novel antigens, as well as the individualized treatments of neoantigens, were first introduced. Secondly, the mechanism of Chimeric Antigen Receptor T-Cell Immunotherapy (CAR-T) therapy and immune checkpoint blockade therapy in the treatment of tumors were outlined, and the three treatment methods were compared. Thirdly, the application of neoantigens in CAR-T therapy and PD-1/PD-L1 blockade therapy was briefly described. The benefits of the neoantigen vaccines over common vaccines were summarized as well. Finally, the prospect of neoantigen therapy was presented.

Knockdown of hnRNPA1 Promotes NSCLC Metastasis and EMT by Regulating Alternative Splicing of LAS1L exon 9

Tumor metastasis is still an insurmountable obstacle in tumor treatment. Lung cancer represents one of the most common malignancies with high morbidity worldwide. hnRNPA1 has been reported to be involved in the regulation of tumor metastasis, while its specific role in tumor metastasis seems to be controversial and its molecular mechanism in lung cancer metastasis remains to be further elucidated. In this study, we confirmed that knockdown of the hnRNPA1 led to enhanced migration, invasion and EMT transition in lung cancer cells. Bioinformatics analysis of the GSE34992 dataset revealed that hnRNPA1 may regulate the alternative splicing (AS) of LAS1L exon 9. Further AGE assays and RIP assays revealed that hnRNPA1 can directly bind to the LAS1L pre-mRNA to inhibit the splicing of LAS1L exon 9. The RNA pull-down assays showed that hnRNPA1 can specifically bind to the two sites (UAGGGU(WT1) and UGGGGU(WT3)) of LAS1L Intron 9. Further Transwell assays indicated that the expression ratio of LAS1L-L/LAS1L-S regulated by hnRNPA1 can further promote the migration, invasion and EMT transition in lung cancer cells. Moreover, hnRNPA1 expression showed significant heterogeneity in lung cancer tissues, which may contain new research directions and potential therapeutic targets. Our results indicate that hnRNPA1 can affect the metastasis of lung cancer cells by modulating the AS of LAS1L exon 9, highlighting the potential significance of hnRNPA1 in lung cancer metastasis.

Derivation of splice junction-specific antibodies using a unique hapten targeting strategy and directed evolution

Alternative splicing of RNA occurs frequently in eukaryotic cells and can result in multiple protein isoforms that are nearly identical in amino acid sequence, but have unique biological roles. Moreover, the relative abundance of these unique isoforms can be correlative with diseased states and potentially used as biomarkers or therapeutic targets. However, due to high sequence similarities among isoforms, current proteomic methods are incapable of differentiating native protein isoforms derived from most alternative splicing events. Herein, a strategy employing a nonsynonymous, non-native amino acid (nnAA) pseudo-hapten (i.e. an amino acid or amino acid derivative that is different from the native amino acid at a particular position) as a targeting epitope in splice junction-spanning peptides was successful in directed antibody derivation. After isolating nnAA-specific antibodies, directed evolution reduced the antibody's binding dependence on the nnAA pseudo-hapten and improved binding to the native splice junction epitope. The resulting antibodies demonstrated codependent binding affinity to each exon of the splice junction and thus are splice junction- and isoform-specific. Furthermore, epitope scanning demonstrated that positioning of the nnAA pseudo-hapten within a peptide antigen can be exploited to predetermine the isolated antibody's specificity at, or near, amino acid resolution. Thus, this nnAA targeting strategy has the potential to robustly derive splice junction- and site-specific antibodies that can be used in a wide variety of research endeavors to unambiguously differentiate native protein isoforms.

A Novel Proteogenomic Integration Strategy Expands the Breadth of Neo-Epitope Sources

Tumor-specific antigens can activate T cell-based antitumor immune responses and are ideal targets for cancer immunotherapy. However, their identification is still challenging. Although mass spectrometry can directly identify human leukocyte antigen (HLA) binding peptides in tumor cells, it focuses on tumor-specific antigens derived from annotated protein-coding regions constituting only 1.5% of the genome. We developed a novel proteogenomic integration strategy to expand the breadth of tumor-specific epitopes derived from all genomic regions. Using the colorectal cancer cell line HCT116 as a model, we accurately identified 10,737 HLA-presented peptides, 1293 of which were non-canonical peptides that traditional database searches could not identify. Moreover, we found eight tumor neo-epitopes derived from somatic mutations, four of which were not previously reported. Our findings suggest that this new proteogenomic approach holds great promise for increasing the number of tumor-specific antigen candidates, potentially enlarging the tumor target pool and improving cancer immunotherapy.

HLA-I-restricted CD8+ T cell immunity may accelerate tumorigenesis in conjunction with VHL inactivation

CD8⁺ T cells recognize and kill tumor cells with HLA-I tumor antigens in early tumorigenesis, the efficiency of which differs according to antigen-recognition coverage, as shown in earlier tumor onset in HLA-I homozygosity. However, the universality of these associations remains unknown. Here, we assessed the tumor type and driver mutation specificity in the association between tumor onset age and HLA-I zygosity. Statistical analyses identified an unexpected negative relationship in tumors with VHL biallelic loss, wherein HLA-I heterozygosity was associated with earlier tumor onset, while all others showed either no or a positive association. Testing on an independent dataset reproduced the VHL-dependent acceleration of tumor onset in the HLA-I heterozygous group, confirming the association. Further speculation proposed VEGF-A-mediated T cell exhaustion under VHL inactivation as a potential mechanism. Our findings suggest that CD8⁺ T cell immunity in early tumor suppression can be conditional to the genetic status of tumors and may even lead to adverse consequences.

•

HLA homozygosity reduces antigen coverage and is associated with earlier tumor onset

•

Tumors with VHL^−/−, such as ccRCC, have the opposite association

•

In VHL^−/− tumors, CD8⁺ T cell immunity may have adverse effects in imunosurveillance

The Oncogenic PI3K-Induced Transcriptomic Landscape Reveals Key Functions in Splicing and Gene Expression Regulation

The phosphoinositide 3-kinase (PI3K) pathway regulates proliferation, survival, and metabolism and is frequently activated across human cancers. A comprehensive elucidation of how this signaling pathway controls transcriptional and cotranscriptional processes could provide new insights into the key functions of PI3K signaling in cancer. Here, we undertook a transcriptomic approach to investigate genome-wide gene expression and transcription factor activity changes, as well as splicing and isoform usage dynamics, downstream of PI3K. These analyses uncovered widespread alternatively spliced isoforms linked to proliferation, metabolism, and splicing in PIK3CA-mutant cells, which were reversed by inhibition of PI3Kα. Analysis of paired tumor biopsies from patients with PIK3CA-mutated breast cancer undergoing treatment with PI3Kα inhibitors identified widespread splicing alterations that affect specific isoforms in common with the preclinical models, and these alterations, namely PTK2/FRNK and AFMID isoforms, were validated as functional drivers of cancer cell growth or migration. Mechanistically, isoform-specific splicing factors mediated PI3K-dependent RNA splicing. Treatment with splicing inhibitors rendered breast cancer cells more sensitive to the PI3Kα inhibitor alpelisib, resulting in greater growth inhibition than alpelisib alone. This study provides the first comprehensive analysis of widespread splicing alterations driven by oncogenic PI3K in breast cancer. The atlas of PI3K-mediated splicing programs establishes a key role for the PI3K pathway in regulating splicing, opening new avenues for exploiting PI3K signaling as a therapeutic vulnerability in breast cancer.

SIGNIFICANCE

Transcriptomic analysis reveals a key role for the PI3K pathway in regulating RNA splicing, uncovering new mechanisms by which PI3K regulates proliferation and metabolism in breast cancer. See related commentary by Claridge and Hopkins, p. 2216.

Proteomic Profiling Identifies Co-Regulated Expression of Splicing Factors as a Characteristic Feature of Intravenous Leiomyomatosis

Intravenous leiomyomatosis (IVLM) is a rare benign smooth muscle tumour that is characterised by intravenous growth in the uterine and pelvic veins. Previous DNA copy number and transcriptomic studies have shown that IVLM harbors unique genomic and transcriptomic alterations when compared to uterine leiomyoma (uLM), which may account for their distinct clinical behaviour. Here we undertake the first comparative proteomic analysis of IVLM and other smooth muscle tumours (comprising uLM, soft tissue leiomyoma and benign metastasizing leiomyoma) utilising data-independent acquisition mass spectrometry. We show that, at the protein level, IVLM is defined by the unique co-regulated expression of splicing factors. In particular, IVLM is enriched in two clusters composed of co-regulated proteins from the hnRNP, LSm, SR and Sm classes of the spliceosome complex. One of these clusters (Cluster 3) is associated with key biological processes including nascent protein translocation and cell signalling by small GTPases. Taken together, our study provides evidence of co-regulated expression of splicing factors in IVLM compared to other smooth muscle tumours, which suggests a possible role for alternative splicing in the pathogenesis of IVLM.

GLTSCR1 coordinates alternative splicing and transcription elongation of ZO1 to regulate colorectal cancer progression

Alternative splicing (AS) and transcription elongation are vital biological processes, and their dysregulation causes multiple diseases, including tumors. However, the coregulatory mechanism of AS and transcription elongation in tumors remains unclear. This study demonstrates a novel AS pattern of tight junction protein 1 (ZO1) regulated by the RNA polymerase II elongation rate in colorectal cancer (CRC). Glioma tumor suppressor candidate region gene 1 (GLTSCR1) decreases the transcription elongation rate of ZO1 to provide a time window for binding of the splicing factor HuR to the specific motif in intron 22 of ZO1 and spliceosome recognition of the weak 3' and 5' splice sites in exon 23 to promote exon 23 inclusion. Since exon 23 inclusion in ZO1 suppresses migration and invasion of CRC cells, our findings suggest a novel potential therapeutic target for CRC.

Novel γδ T cell-based prognostic signature to estimate risk and aid therapy in hepatocellular carcinoma

BACKGROUND

Numerous studies have revealed that gamma delta (γδ) T cell infiltration plays a crucial regulatory role in hepatocellular carcinoma (HCC) development. Nonetheless, a comprehensive analysis of γδ T cell infiltration in prognosis evaluation and therapeutic prediction remains unclear.

METHODS

Multi-omic data on HCC patients were obtained from public databases. The CIBERSORT algorithm was applied to decipher the tumor immune microenvironment (TIME) of HCC. Weighted gene co-expression network analysis (WGCNA) was performed to determine significant modules with γδ T cell-specific genes. Kaplan-Meier survival curves and receiver operating characteristic analyses were used to validate prognostic capability. Additionally, the potential role of RFESD inhibition by si-RFESD in vitro was investigated using EdU and CCK-8 assays.

RESULTS

A total of 16,421 genes from 746 HCC samples (616 cancer and 130 normal) were identified based on three distinct cohorts. Using WGCNA, candidate modules (brown) with 1755 significant corresponding genes were extracted as γδ T cell-specific genes. Next, a novel risk signature consisting of 11 hub genes was constructed using multiple bioinformatic analyses, which presented great prognosis prediction reliability. The risk score exhibited a significant correlation with ICI and chemotherapeutic targets. HCC samples with different risks experienced diverse signalling pathway activities. The possible interaction of risk score with tumor mutation burden (TMB) was further analyzed. Subsequently, the potential functions of the RFESD gene were explored in HCC, and knockdown of RFESD inhibited cell proliferation in HCC cells. Finally, a robust prognostic risk-clinical nomogram was developed and validated to quantify clinical outcomes.

CONCLUSIONS

Collectively, comprehensive analyses focusing on γδ T cell patterns will provide insights into prognosis prediction, the mechanisms of immune infiltration, and advanced therapy strategies in HCC.

Clonally expanded CD8 T cells characterize amyotrophic lateral sclerosis-4

Amyotrophic lateral sclerosis (ALS) is a heterogenous neurodegenerative disorder that affects motor neurons and voluntary muscle control1. ALS heterogeneity includes the age of manifestation, the rate of progression and the anatomical sites of symptom onset. Disease-causing mutations in specific genes have been identified and define different subtypes of ALS1. Although several ALS-associated genes have been shown to affect immune functions2, whether specific immune features account for ALS heterogeneity is poorly understood. Amyotrophic lateral sclerosis-4 (ALS4) is characterized by juvenile onset and slow progression3. Patients with ALS4 show motor difficulties by the time that they are in their thirties, and most of them require devices to assist with walking by their fifties. ALS4 is caused by mutations in the senataxin gene (SETX). Here, using Setx knock-in mice that carry the ALS4-causative L389S mutation, we describe an immunological signature that consists of clonally expanded, terminally differentiated effector memory (TEMRA) CD8 T cells in the central nervous system and the blood of knock-in mice. Increased frequencies of antigen-specific CD8 T cells in knock-in mice mirror the progression of motor neuron disease and correlate with anti-glioma immunity. Furthermore, bone marrow transplantation experiments indicate that the immune system has a key role in ALS4 neurodegeneration. In patients with ALS4, clonally expanded TEMRA CD8 T cells circulate in the peripheral blood. Our results provide evidence of an antigen-specific CD8 T cell response in ALS4, which could be used to unravel disease mechanisms and as a potential biomarker of disease state.

NSrp70 suppresses metastasis in triple-negative breast cancer by modulating Numb/TβR1/EMT axis

Alternative splicing of mRNA precursors allows cancer cells to create different protein isoforms that promote growth and survival. Compared to normal cells, cancer cells frequently exhibit a higher diversity of their transcriptomes. A comprehensive understanding of splicing regulation is required to correct the splicing alterations for the future precision oncology. A quantitative proteomic screen was performed to identify the regulators associated the metastasis in triple-negative breast cancer. Multiple in vitro and in vivo functional analyses were used to study the effects of NSrp70 on breast cancer metastasis. Next, transcriptomic sequencing (RNA-seq) and alternative splicing bioinformatics analysis was applied to screen the potential targets of NSrp70. Moreover, in vitro splicing assays, RNA pull-down, and RNA immunoprecipitation assay were used to confirm the specific binding between NSrp70 and downstream target genes. Furthermore, the prognostic value of NSrp70 was analyzed in a cohort of patients by performing IHC. We uncovered NSrp70 as a novel suppressor of breast cancer metastasis. We discovered that NSrp70 inhibited the skipped exon alternative splicing of NUMB, promoted the degradation of transforming growth factor receptor 1 through lysosome pathway, and regulated TGFβ/SMAD-mediated epithelial-mesenchymal transition phenotype in breast cancer cells. Furthermore, high NSrp70 expression correlated with a better prognosis in breast cancer patients. Our findings revealed that splicing regulator NSrp70 serves as a metastasis suppressor.

Editorial overview: Vaccines: Reinvigorating therapeutic cancer vaccines

Lessons learned from the rapid deployment of vaccines during the COVID-19 pandemic are reinvigorating the cancer vaccine field. Using delivery platforms including mRNA and synthetic long peptides, recent clinical trials have demonstrated that cancer vaccines are safe, feasible, and can be associated with the generation of antigen-specific memory T cells and, in some cases, durable clinical responses. Despite these advances, fundamental questions remain regarding the optimal delivery platforms and antigen targets to use in cancer vaccines. Ongoing and future studies that harness advances in the identification of novel sources of antigens, the prediction of immunogenic antigens, and the use of single-cell technologies to profile antigen-specific T cells will hopefully reveal correlates with clinical outcomes and provide a mechanistic basis for future progress.

Turning cold tumors hot: from molecular mechanisms to clinical applications

Immune checkpoint blockade (ICB) therapies have achieved clinical benefit, but most 'immune-cold' solid tumors are not responsive. The diversity of immune evasion mechanisms remains a key obstacle in turning nonresponsive 'cold' tumors into responsive 'hot' ones. Therefore, exploring the mechanisms of such transitions and tumor immunotyping can provide significant insights into designing effective therapeutic strategies against cancer. Here, we focus on the latest advances regarding local and systemic regulatory mechanisms of immune responses in cold and hot tumors. We also highlight the necessity for tumor immunotyping through the assessment of multiple immunological variables using various diagnostic techniques and biomarkers. Finally, we discuss the challenges and potential clinical applications of immunophenotyping to turn cold tumors hot, which may further guide combined immunotherapies.

Identification of gene signatures for COAD using feature selection and Bayesian network approaches

The combination of TCGA and GTEx databases will provide more comprehensive information for characterizing the human genome in health and disease, especially for underlying the cancer genetic alterations. Here we analyzed the gene expression profile of COAD in both tumor samples from TCGA and normal colon tissues from GTEx. Using the SNR-PPFS feature selection algorithms, we discovered a 38 gene signatures that performed well in distinguishing COAD tumors from normal samples. Bayesian network of the 38 genes revealed that DEGs with similar expression patterns or functions interacted more closely. We identified 14 up-DEGs that were significantly correlated with tumor stages. Cox regression analysis demonstrated that tumor stage, STMN4 and FAM135B dysregulation were independent prognostic factors for COAD survival outcomes. Overall, this study indicates that using feature selection approaches to select key gene signatures from high-dimensional datasets can be an effective way for studying cancer genomic characteristics.

Altered splicing of ATG16-L1 mediates acquired resistance to tyrosine kinase inhibitors of EGFR by blocking autophagy in non-small cell lung cancer

Despite the initial efficacy of using tyrosine kinase inhibitors of epidermal growth factor receptor (EGFR-TKIs) for treating patients with non-small cell lung cancer (NSCLC), resistance inevitably develops. Recent studies highlight a link between alternative splicing and cancer drug response. Therefore, we aimed to identify deregulated splicing events that play a role in resistance to EGFR-TKI. By using RNA sequencing, reverse transcription PCR (RT-PCR) and RNA interference, we showed that overexpression of a splice variant of the autophagic gene ATG16-L1 that retains exon 8 and encodes the β-isoform of autophagy-related protein 16-1 (ATG16-L1-β) concurs acquired resistance to EGFR-TKI in NSCLC cells. Using matched biopsies, we found increased levels of ATG16-L1-β at the time of progression in 3 of 11 NSCLC patients treated with EGFR-TKI. Mechanistically, gefitinib-induced autophagy was impaired in resistant cells that accumulated ATG16-L1-β. Neutralization of ATG16-L1-β restored autophagy in response to gefitinib, induced apoptosis and inhibited the growth of in ovo tumor xenografts. Conversely, overexpression of ATG16-L1-β in parental sensitive cells prevented gefitinib-induced autophagy and increased cell survival. These results support a role for defective autophagy in acquired resistance to EGFR-TKIs and identify splicing regulation of ATG16-L1 as a therapeutic vulnerability that could be explored for improving EGFR-targeted cancer therapy.

Alternative Splicing of the Aryl Hydrocarbon Receptor Nuclear Translocator (ARNT) Is Regulated by RBFOX2 in Lymphoid Malignancies

Aberrant alternative splicing (AS) of pre-mRNAs promotes the development and proliferation of cancerous cells. Accordingly, we had previously observed higher levels of the aryl hydrocarbon receptor nuclear translocator (ARNT) spliced variant isoform 1 in human lymphoid malignancies compared to that in normal lymphoid cells, which is a consequence of increased inclusion of alternative exon 5. ARNT is a transcription factor that has been implicated in the survival of various cancers. Notably, we found that ARNT isoform 1 promoted the growth and survival of lymphoid malignancies, but the regulatory mechanism controlling ARNT AS is unclear. Here, we report cis- and trans-regulatory elements which are important for the inclusion of ARNT exon 5. Specifically, we identified recognition motifs for the RNA-binding protein RBFOX2, which are required for RBFOX2-mediated exon 5 inclusion. RBFOX2 upregulation was observed in lymphoid malignancies, correlating with the observed increase in ARNT exon 5 inclusion. Moreover, suppression of RBFOX2 significantly reduced ARNT isoform 1 levels and cell growth. These observations reveal RBFOX2 as a critical regulator of ARNT AS in lymphoid malignancies and suggest that blocking the ARNT-specific RBFOX2 motifs to decrease ARNT isoform 1 levels is a viable option for targeting the growth of lymphoid malignancies.

Immune-Based Therapeutic Interventions for Acute Myeloid Leukemia

Acute myeloid leukemia (AML) is an aggressive, clonally heterogeneous, myeloid malignancy, with a 5-year overall survival of approximately 27%. It constitutes the most common acute leukemia in adults, with an incidence of 3-5 cases per 100,000 in the United States. Despite great advances in understanding the molecular mechanisms underpinning leukemogenesis, the past several decades had seen little change to the backbone of therapy, comprised of an anthracycline-based induction regimen for those who are fit enough to receive it, followed by risk-stratified post-remission therapy with consolidation cytarabine or allogeneic stem cell transplantation (allo-SCT). Allo-SCT is the most fundamental form of immunotherapy in which donor cytotoxic T and NK cells recognize and eradicate residual AML in the graft-versus-leukemia (GvL) effect. Building on that, several alternative or synergistic approaches to exploit both self and foreign immunity against AML have been developed. Checkpoint inhibitors, for example, CTLA-4 inhibitors, PD-1 inhibitors, and PD-L1 inhibitors block proteins found on T cells or cancer cells that stop the immune system from attacking the cancer cells. They have been used with limited success in both the AML relapsed/refractory (R/R) and post SCT settings. AML tumor mutational burden is low compared to solid tumors and thus, it is less likely to generate neoantigens and respond to antibody-mediated checkpoint blockade that has shown unprecedented results in solid tumors. Therefore, alternative therapeutic strategies that work independently of the T cell receptor (TCR) specificity have been developed. They include bispecific antibodies, which recruit T cells through CD3 engagement, and in AML have shown an overall response rate ranging between 14 and 30% in early phase trials. Chimeric Antigen Receptor (CAR) T cell therapy is a type of treatment in which T cells are genetically engineered to produce a recombinant receptor that redirects the specificity and function of T lymphocytes. However, lack of cell surface targets exclusively expressed on AML cells including Leukemic Stem Cells (LSCs) combined with clonal heterogeneity represents the biggest challenge in developing CAR therapy for AML. Antibody-Drug Conjugates (ADC) constitute the only FDA-approved immunotherapy to treat AML with Gemtuzumab Ozogamicin, a CD33-specific ADC used in CEBPα-mutated AML. The identification of additional cell surface targets is critical for the development of other ADC's potentially useful in the induction and maintenance regimens, given the ease at which these reagents can be generated and managed. Here, we will review those immune-based therapeutic interventions and highlight active areas of research investigations toward fulfillment of the great promise of immunotherapy to AML.