Alternative RNA Splicing as a Potential Major Source of Untapped Molecular Targets in Precision Oncology and Cancer Disparities

Genetic ancestry concordant RNA splicing in prostate cancer involves oncogenic genes and associates with recurrence

Black men suffer disproportionately from prostate cancer (PCa) compared to men of other races and ethnicities. Comparing the molecular landscape of PCa among Black and White patients has the potential to identify targets for development of new precision medicine interventions. Herein, we conducted transcriptomic analysis of prostate tumors and paired tumor-adjacent normals from self-reported Black and White PCa patients and estimated patient genetic ancestry. Clinical follow-up revealed increased biochemical recurrence (BCR) among Black patients compared to White patients with high-grade PCa. Transcriptomic analysis identified differential alternative RNA splicing events (ARSs) between Black and White PCa patients. Genes undergoing genetic ancestry-concordant ARSs in high-grade or low-grade tumors involved cancer promoting genes. Most genes undergoing genetic ancestry-concordant ARSs did not exhibit differential aggregate gene expression or alternative polyadenylation. A number of the genetic ancestry-concordant ARSs associated with BCR; thus, genetic ancestry-concordant RNA splice variants may represent unique targets for PCa precision oncology.

Global profiling of alternative splicing in non-small cell lung cancer reveals novel histological and population differences

Lung cancer is one of the most frequently diagnosed cancers in the US. African-American (AA) men are more likely to develop lung cancer with higher incidence and mortality rates than European-American (EA) men. Herein, we report high-confidence alternative splicing (AS) events from high-throughput, high-depth total RNA sequencing of lung tumors and non-tumor adjacent tissues (NATs) in two independent cohorts of patients with adenocarcinoma (LUAD) and squamous cell carcinoma (LUSC). We identified novel AS biomarkers with notable differential percent spliced in (PSI) values between lung tumors and NATs enriched in the AA and EA populations, which were associated with oncogenic signaling pathways. We also uncovered tumor subtype- and population-specific AS events associated with cell surface proteins and cancer driver genes. We highlighted significant AS events in SYNE2 specific to LUAD in both populations, as well as those in CD44 from EAs and TMBIM6 from AAs specific to LUAD. Here, we also present the validation of cancer signatures based on direct high-throughput reverse transcription-PCR. Our large survey of lung tumors presents a rich data resource that may help to understand molecular subtypes of lung tumor between AAs and EAs and reveal new therapeutic vulnerabilities that potentially advance health equity.

Integrative analysis based on ATAC-seq and RNA-seq reveals a novel oncogene PRPF3 in hepatocellular carcinoma

BACKGROUND

Assay of Transposase Accessible Chromatin Sequencing (ATAC-seq) is a high-throughput sequencing technique that detects open chromatin regions across the genome. These regions are critical in facilitating transcription factor binding and subsequent gene expression. Herein, we utilized ATAC-seq to identify key molecular targets regulating the development and progression of hepatocellular carcinoma (HCC) and elucidate the underlying mechanisms.

METHODS

We first compared chromatin accessibility profiles between HCC and normal tissues. Subsequently, RNA-seq data was employed to identify differentially expressed genes (DEGs). Integrating ATAC-seq and RNA-seq data allowed the identification of transcription factors and their putative target genes associated with differentially accessible regions (DARs). Finally, functional experiments were conducted to investigate the effects of the identified regulatory factors and corresponding targets on HCC cell proliferation and migration.

RESULTS

Enrichment analysis of DARs between HCC and adjacent normal tissues revealed distinct signaling pathways and regulatory factors. Upregulated DARs in HCC were enriched in genes related to the MAPK and FoxO signaling pathways and associated with transcription factor families like ETS and AP-1. Conversely, downregulated DARs were associated with the TGF-β, cAMP, and p53 signaling pathways and the CTCF family. Integration of the datasets revealed a positive correlation between specific DARs and DEGs. Notably, PRPF3 emerged as a gene associated with DARs in HCC, and functional assays demonstrated its ability to promote HCC cell proliferation and migration. To the best of our knowledge, this is the first report highlighting the oncogenic role of PRPF3 in HCC. Furthermore, ZNF93 expression positively correlated with PRPF3, and ChIP-seq data indicated its potential role as a transcription factor regulating PRPF3 by binding to its promoter region.

CONCLUSION

This study provides a comprehensive analysis of the epigenetic landscape in HCC, encompassing both chromatin accessibility and the transcriptome. Our findings reveal that ZNF93 promotes the proliferation and motility of HCC cells through transcriptional regulation of a novel oncogene, PRPF3.

Isoform Switching Regulates the Response to Ionizing Radiation Through SRSF1

PURPOSE

This study investigated how isoform switching affects the cellular response to ionizing radiation (IR), an understudied area despite its relevance to radiation therapy in cancer treatment. We aimed to identify changes in transcript isoform expression post-IR exposure and the proteins mediating these changes, with a focus on their potential to modulate radiosensitivity.

METHODS AND MATERIALS

Using RNA sequencing, we analyzed the B-cell lines derived from 10 healthy individuals at 3 timepoints, applying the mixture of isoforms algorithm to quantify alternative splicing. We examined RNA binding protein motifs within the sequences of IR-responsive isoforms and validated the serine/arginine-rich splicing factor 1 (SRSF1) as a predominant mediator through RNA immunoprecipitation. We further investigated the effects of SRSF1 on radiosensitivity by RNA interference and by analyzing publicly available data on patients with cancer.

RESULTS

We identified ∼1900 radiation-responsive alternatively spliced isoforms. Many isoforms were differentially expressed without changes in their overall gene expression. Over a third of these transcripts underwent exon skipping, while others used proximal last exons. These IR-responsive isoforms tended to be shorter transcripts missing vital domains for preventing apoptosis and promoting cell division but retaining those necessary for DNA repair. Our combined computational, genetic, and molecular analyses identified the proto-oncogene SRSF1 as a mediator of these radiation-induced isoform-switching events that promote apoptosis. After exposure to DNA double-strand break-inducing agents, SRSF1 expression decreased. A reduction in SRSF1 increased radiosensitivity in vitro and among patients with cancer.

CONCLUSIONS

We establish a pivotal role for isoform switching in the cellular response to IR and propose SRSF1 as a promising biomarker for assessing radiation therapy effectiveness.

The 29th Annual Prostate Cancer Foundation Scientific Retreat Report

BACKGROUND

The 29th Annual Prostate Cancer Foundation (PCF) Scientific Retreat was held from October 27 to 29, 2022, at the Omni La Costa Resort in Carlsbad, CA. This was the first-ever hybrid PCF Retreat.

METHODS

The Annual PCF Scientific Retreat is a prominent international scientific gathering centered on groundbreaking, unpublished, and influential studies in basic, translational, and clinical prostate cancer research. It also covers research from related fields with a strong potential for influencing prostate cancer research and patient care.

RESULTS

Key areas of research that were focused on at the 2022 PCF Retreat included: (i) the contributions of molecular and genomic factors to prostate cancer disparities; (ii) novel clinical trial updates; (iii) lessons from primary prostate cancer; (iv) lessons from single-cell studies; (v) genetic, epigenetic, epitranscriptomic and posttranslational mechanisms and clinical heterogeneity in prostate cancer; (vi) biology of neuroendocrine and lineage-plastic prostate cancer; (vii) next generation prostate cancer theranostics and combination therapies; (viii) the biology and therapeutic potential of targeting phosphoinositide 3-kinases pathways; (ix) combining immunomodulatory treatments for prostate cancer; (x) novel gamma delta (γδ) T-cell therapy platforms for oncology; and (xi) lessons from other cancers.

CONCLUSIONS

This article provides a summary of the presentations from the 2022 PCF Scientific Retreat. By disseminating this knowledge, we hope to enhance our understanding of the present research landscape and guide future strides in both prostate cancer research and patient care.

Long-Read Sequencing Reveals Alternative Splicing-Driven, Shared Immunogenic Neoepitopes Regardless of SF3B1 Status in Uveal Melanoma

Tumor-specific neoepitopes are promising targets in cancer immunotherapy. However, the identification of functional tumor-specific neoepitopes remains challenging. In addition to the most common source, single-nucleotide variants (SNV), alternative splicing (AS) represents another rich source of neoepitopes and can be utilized in cancers with low SNVs such as uveal melanoma (UM). UM, the most prevalent adult ocular malignancy, has poor clinical outcomes due to a lack of effective therapies. Recent studies have revealed the promise of harnessing tumor neoepitopes to treat UM. Previous studies have focused on neoepitope targets associated with mutations in splicing factor 3b subunit 1 (SF3B1), a key splicing factor; however, little is known about the neoepitopes that are commonly shared by patients independent of SF3B1 status. To identify the AS-derived neoepitopes regardless of SF3B1 status, we herein used a comprehensive nanopore long-read-sequencing approach to elucidate the landscape of AS and novel isoforms in UM. We also performed high-resolution mass spectrometry to further validate the presence of neoepitope candidates and analyzed their structures using the AlphaFold2 algorithm. We experimentally evaluated the antitumor effects of these neoepitopes and found they induced robust immune responses by stimulating interferon (IFN)γ production and activating T cell-based UM tumor killing. These results provide novel insights into UM-specific neoepitopes independent of SF3B1 and lay the foundation for developing therapies by targeting these actionable neoepitopes.

Mitochondrial DNA methylation is a predictor of immunotherapy response and prognosis in breast cancer: scRNA-seq and bulk-seq data insights

Background

Alterations in Mitochondrial DNA methylation (MTDM) exist in many tumors, but their role in breast cancer (BC) development remains unclear.

Methods

We analyzed BC patient data by combining scRNA-seq and bulk sequencing. Weighted co-expression network analysis (WGCNA) of TCGA data identified mitochondrial DNA methylation (MTDM)-associated genes in BC. COX regression and LASSO regression were used to build prognostic models. The biological function of MTDM was assessed using various methods, such as signaling pathway enrichment analysis, copynumber karyotyping analysis, and quantitative analysis of the cell proliferation rate. We also evaluated MTDM-mediated alterations in the immune microenvironment using immune microenvironment, microsatellite instability, mutation, unsupervised clustering, malignant cell subtype differentiation, immune cell subtype differentiation, and cell-communication signature analyses. Finally, we performed cellular experiments to validate the role of the MTDM-associated prognostic gene NCAPD3 in BC.

Results

In this study, MTDM-associated prognostic models divided BC patients into high/low MTDM groups in TCGA/GEO datasets. The difference in survival time between the two groups was statistically significant (P<0.001). We found that high MTDM status was positively correlated with tumor cell proliferation. We analyzed the immune microenvironment and found that low-MTDM group had higher immune checkpoint gene expression/immune cell infiltration, which could lead to potential benefits from immunotherapy. In contrast, the high MTDM group had higher proliferation rates and levels of CD8+T cell exhaustion, which may be related to the secretion of GDF15 by malignant breast epithelial cells with a high MTDM status. Cellular experiments validated the role of the MTDM-associated prognostic gene NCAPD3 (the gene most positively correlated with epithelial malignant cell proliferation in the model) in BC. Knockdown of NCAPD3 significantly reduced the activity and proliferation of MDA-MB-231 and BCAP-37 cells, and significantly reduced their migration ability of BCAP-37 cell line.

Conclusion

This study presented a holistic evaluation of the multifaceted roles of MTDM in BC. The analysis of MTDM levels not only enables the prediction of response to immunotherapy but also serves as an accurate prognostic indicator for patients with BC. These insightful discoveries provide novel perspectives on tumor immunity and have the potentially to revolutionize the diagnosis and treatment of BC.

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.

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.

Phosphatidylinositol glycan anchor biosynthesis, class C is a prognostic biomarker and correlates with immune infiltrates in hepatocellular carcinoma

Background and aims: The exact function of Phosphatidylinositol Glycan Anchor Biosynthesis, Class C (PIGC) gene has yet to be elucidated. In the study, we attempted to clarify the correlations of PIGC to prognosis and tumor-infiltrating lymphocytes in hepatocellular carcinoma (HCC).

Methods:PIGC expression was analyzed via the Oncomine database, Gene Expression Profiling Interactive Analysis, Hepatocellular carcinoma data base, Human Protein Atlas database and Tumor Immune Estimation Resource (TIMER). We showed the correlation of PIGC with the clinical characteristics using UALCAN. We evaluated the influence of PIGC on clinical prognosis using Kaplan-Meier plotter databases. And co-expressed genes with PIGC and its regulators were identified using LinkedOmics. The correlations between PIGC and cancer immune infiltrates were investigated via TIMER. We analyzed the drug sensitivity and immunotherapy response via R package.

Results:PIGC was found up-regulated in tumor tissues in multiple HCC cohorts, also increased in HCC patient with different clinical characteristics. High PIGC expression was associated with poorer overall survival. PIGC expression showed a strong positive association with the expression of ACBD6, a strong negative association with AGXT212. The cell components and distribution in treatment and non-treatment of HCC patients were quite distinct, which may reveal the relationship between the immunotherapy with tumor microenvironment. Notably, PIGC expression was positively correlated with infiltrating levels of immune cells.

Conclusion: These findings suggest that PIGC is correlated with prognosis and immune infiltrating in HCC, which can be used as a prognostic biomarker for determining prognosis, laying a foundation for further study of the immune regulatory role of PIGC in HCC.

An Integrated Analysis of the Identified PRPF19 as an Onco-immunological Biomarker Encompassing the Tumor Microenvironment, Disease Progression, and Prognoses in Hepatocellular Carcinoma

Background: Targeting the mRNA splicing process has been identified as a therapeutic strategy for human cancer. PRPF19 is an RNA binding protein that is involved in pre-mRNA processing and repairing DNA damage; the aberrant expression of PRPF19 is potentially associated with carcinogenesis. However, the biological role of PRPF19 in hepatocellular carcinoma (HCC) is still elusive.

Methods: Data obtained from TCGA, Oncomine, and GEO were used to investigate the PRPF19 expression level and its role in tumor immune infiltration, prognosis, and the tumor progression of cohorts from HCC. Using various databases and tools (UALCAN, TIMER, TISMO, and PathCards), we presented the potential mechanisms of PFPF19 upregulation, PRPF19-related pathways, and its biological functions in liver cancer.

Results: For HCC, PRPF19 expression was found upregulated both in single tumor cells and tissues. Furthermore, the increased expression of PRPF19 was significantly correlated to clinical characteristics: advanced stage, vascular invasion, high AFP, and poor prognosis of HCC. According to the tumor-immunological analysis, we found that PRPF19 is positively correlated with infiltrating myeloid-derived suppressor cells (MDSCs). Moreover, the microenvironment of HCC tissues with high expression of PRPF19 is highly immunosuppressive (lower T-lymphocytes, multiple immune checkpoints upregulated). Patients with high expression of PRPF19 and high MDSCs had a worse survival prognosis as well. TP53 mutation may have a positive effect on PRPF19 expression via decreased promoter methylation of PRPF19. By TF-mRNA network analysis, key transcription factors (TFs) in TC-NER and PCS pathways (PRPF19 involved) were identified.

Conclusion: This work implied that PRPF19 is associated with tumor immune evasion and progression, and serves as a prognostic marker for worse clinical outcomes with HCC. Thus, this critical regulator could serve as a potential therapeutic target of HCC.

RBPMS2, as a novel biomarker for predicting lymph node metastasis, guides therapeutic regimens in gastric cancer

RNA-binding protein with multiple splicing 2 (RBPMS2) is a critical gene that encodes a member of the RNA-recognition-motif-containing protein family and is involved in the development and dedifferentiation of digestive smooth muscle cells. However, whether RBPMS2 has an effect on the prognosis of gastric cancer (GC) and its possible mechanism during GC progression remain unknown. Here, we collected 596 GC patients who underwent curative surgical resection to evaluate expression of RBPMS2. RBPMS2 was upregulated in GC tissues, and was associated with lymph node metastasis. We analysis the KEGG pathway and GO biological processes, cellular component and molecular function of RBPMS2 interactive genes, as well as RBPMS2-binding partner ESRP1 interactive genes. We also focus on the correlation analysis between EMT-related genes and RBPMS2/ESRP1. Finally, we analysed the correlation between RBPMS2 expression and chemotherapeutic drugs which may assist in GC therapy and demonstrated that RBPMS2 expression was associated with tumour mutation burden (TMB) and Microsatellite Instability (MSI), as well as immune infiltration level in GC.

Genome-Wide DNA Methylation Signatures Predict the Early Asymptomatic Doxorubicin-Induced Cardiotoxicity in Breast Cancer

Chemotherapy with doxorubicin (DOX) may cause unpredictable cardiotoxicity. This study aimed to determine whether the methylation signature of peripheral blood mononuclear cells (PBMCs) prior to and after the first cycle of DOX-based chemotherapy could predict the risk of cardiotoxicity in breast cancer patients. Cardiotoxicity was defined as a decrease in left ventricular ejection fraction (LVEF) by >10%. DNA methylation of PBMCs from 9 patients with abnormal LVEF and 10 patients with normal LVEF were examined using Infinium HumanMethylation450 BeadChip. We have identified 14,883 differentially methylated CpGs at baseline and 18,718 CpGs after the first cycle of chemotherapy, which significantly correlated with LVEF status. Significant differentially methylated regions (DMRs) were found in the promoter and the gene body of SLFN12, IRF6 and RNF39 in patients with abnormal LVEF. The pathway analysis found enrichment for regulation of transcription, mRNA splicing, pathways in cancer and ErbB2/4 signaling. The preliminary results from this study showed that the DNA methylation profile of PBMCs may predict the risk of DOX-induced cardiotoxicity prior to chemotherapy. Further studies with larger cohorts of patients are needed to confirm these findings.

HnRNP-L-regulated circCSPP1/miR-520h/EGR1 axis modulates autophagy and promotes progression in prostate cancer

The circRNAs, a new subclass of non-coding RNAs that are catalyzed by RNA-binding proteins (RBPs), have been reported to be associated with the progression of multiple types of cancer. We previously discovered that heterogeneous nuclear ribonucleoprotein L (HnRNP-L), a multi-functional RBP, is associated with pro-proliferation and anti-apoptosis activities in prostate tumor cells. In this study, we aim to establish the biological relevance of circCSPP1 (a newly discovered signature circRNA in prostate cancer [PCa]) and HnRNP-L to prostate cancer progression. First, we demonstrated that circCSPP1 expression was higher in prostate cancer tissues than in benign tissues and higher in prostate cancer cells than in benign cells. Then, the in vitro gain- and loss-of-function experiments showed that the circCSPP1 expression in prostate cancer cells was regulated by HnRNP-L, and the increased circCSPP1 significantly induced autophagy, which led to an enhanced potential in proliferation, migration, and invasion of prostate cancer cells. These results were consistent with the in vivo experiment where increased or decreased circCSPP1 was associated with higher or slower growth rate in grafted tumors. Finally, we demonstrated the potential competing endogenous RNA network, involving circCSPP1, miR-520h, and early growth response factor 1 (EGR1), in prostate cancer cells, which may play an important role in prostate cancer progression. Our study indicated that the increase in circCSPP1 in prostate cancer, which may be catalyzed by HnRNP-L, can induce cellular autophagy through the circCSPP1-miR-520h-EGR1 axis, leading to the progression of prostate tumor. This newly discovered circRNA biomarker may be used for clinical prognosis of prostate cancer as well as for development of novel therapy plans.

Neoantigen-Reactive T Cells: The Driving Force behind Successful Melanoma Immunotherapy

Simple Summary

Cancer immunotherapy is a revolutionary type of cancer therapy. It uses the patient’s own immune system to fight and potentially cure cancer. The first major breakthrough of immunotherapy came from successful clinical trials for melanoma treatments. Since then, researchers have focused on understanding the science behind immunotherapy, so that patients with other types of cancer may also benefit. One of the major findings is that the T cells in melanoma patients may recognize a specific type of tumor antigen, called neoantigens, and then kill tumor cells that present these neoantigens. The neoantigens mainly arise from the DNA mutations found in tumor cells. These mutations are translated into mutated proteins that are then distinguished by T cells. In this article, we discuss the critical role of T cells in immunotherapy, as well as the clinical trials that shaped the treatments for melanoma.

Abstract

Patients with metastatic cutaneous melanoma have experienced significant clinical responses after checkpoint blockade immunotherapy or adoptive cell therapy. Neoantigens are mutated proteins that arise from tumor-specific mutations. It is hypothesized that the neoantigen recognition by T cells is the critical step for T-cell-mediated anti-tumor responses and subsequent tumor regressions. In addition to describing neoantigens, we review the sentinel and ongoing clinical trials that are helping to shape the current treatments for patients with cutaneous melanoma. We also present the existing evidence that establishes the correlations between neoantigen-reactive T cells and clinical responses in melanoma immunotherapy.

Neoantigen vaccine: An emerging immunotherapy for hepatocellular carcinoma

Tumor-specific neoantigens, which are expressed on tumor cells, can induce an effective antitumor cytotoxic T-cell response and mediate tumor regression. Among tumor immunotherapies, neoantigen vaccines are in early human clinical trials and have demonstrated substantial efficiency. Compared with more neoantigens in melanoma, the paucity and inefficient identification of effective neoantigens in hepatocellular carcinoma (HCC) remain enormous challenges in effectively treating this malignancy. In this review, we highlight the current development of HCC neoantigens in its generation, screening, and identification. We also discuss the possibility that there are more effective neoantigens in hepatitis B virus (HBV)-related HCC than in non-HBV-related HCC. In addition, since HCC is an immunosuppressive tumor, strategies that reverse immunosuppression and enhance the immune response should be considered for the practical exploitation of HCC neoantigens. In summary, this review offers some strategies to solve existing problems in HCC neoantigen research and provide further insights for immunotherapy.

Tumor Fibroblast Growth Factor Receptor 4 Level Predicts the Efficacy of Lenvatinib in Patients With Advanced Hepatocellular Carcinoma

OBJECTIVES

Biomarkers for optimizing the outcome of treatment with lenvatinib in patients with advanced hepatocellular carcinoma remain to be established despite intensive and comprehensive genomic research. Lenvatinib is characterized by its prominent inhibitory potency for fibroblast growth factor receptor (FGFR) 4 compared with earlier tyrosine kinase inhibitors. Thus, in this study, we focused on simplified quantification of FGFR4 in tumors as a potential predictive indicator.

METHODS

According to The Cancer Genome Atlas data set curation, FGFR4 messenger RNA is broadly overexpressed in hepatocellular carcinoma in the absence of gene alteration. Gene set enrichment analysis revealed that the aggressiveness of the tumor was closely related to the FGFR4 level. To confirm the relationship between the benefits of lenvatinib and tumor addiction to the FGFR4 pathway, we analyzed protein levels in tumors and peripheral blood obtained from 57 prospectively registered patients treated with lenvatinib.

RESULTS

Positive immunohistochemistry (>10% of tumor cells) for FGFR4 in biopsy samples before treatment was associated with a longer progression-free survival (2.5 vs 5.5 months, P = 0.01) and a favorable objective response rate (31% vs 81%, P = 0.006). By contrast, the concentration of soluble FGFR4 in peripheral blood as measured by an enzyme-linked immunosorbent assay was not associated with survival outcomes, because its fluctuations reflect hepatic fibrosis. Additional RNA sequencing analysis using archival surgical specimens (n = 90) suggested that alternative RNA splicing of FGFR4 in cancer may also explain this discrepancy.

DISCUSSION

The tumor FGFR4 level was an independent predictor of response to lenvatinib.

SF3B1 mutation in pancreatic cancer contributes to aerobic glycolysis and tumor growth through a PP2A-c-Myc axis

Hot spot gene mutations in splicing factor 3b subunit 1 (SF3B1) are observed in many types of cancer and create abundant aberrant mRNA splicing, which is profoundly implicated in tumorigenesis. Here, we identified that the SF3B1 K700E (SF3B1^K700E) mutation is strongly associated with tumor growth in pancreatic ductal adenocarcinoma (PDAC). Knockdown of SF3B1 significantly retarded cell proliferation and tumor growth in a cell line (Panc05.04) with the SF3B1^K700E mutation. However, SF3B1 knockdown had no notable effect on cell proliferation in two cell lines (BxPC3 and AsPC1) carrying wild‐type SF3B1. Ectopic expression of SF3B1^K700E but not SF3B1^WT in SF3B1‐knockout Panc05.04 cells largely restored the inhibitory role induced by SF3B1 knockdown. Introduction of the SF3B1^K700E mutation in BxPC3 and AsPC1 cells also boosted cell proliferation. Gene set enrichment analysis demonstrated a close correlation between SF3B1 mutation and aerobic glycolysis. Functional analyses showed that the SF3B1^K700E mutation promoted tumor glycolysis, as evidenced by glucose consumption, lactate release, and extracellular acidification rate. Mechanistically, the SF3B1 mutation promoted the aberrant splicing of PPP2R5A and led to the activation of the glycolytic regulator c‐Myc via post‐translational regulation. Pharmacological activation of PP2A with FTY‐720 markedly compromised the growth advantage induced by the SF3B1^K700E mutation in vitro and in vivo. Taken together, our data suggest a novel function for SF3B1 mutation in the Warburg effect, and this finding may offer a potential therapeutic strategy against PDAC with the SF3B1^K700E mutation.

Differential alternative RNA splicing and transcription events between tumors from African American and White patients in The Cancer Genome Atlas

Individuals of African ancestry suffer disproportionally from higher incidence, aggressiveness, and mortality for particular cancers. This disparity likely results from an interplay among differences in multiple determinants of health, including differences in tumor biology. We used The Cancer Genome Atlas (TCGA) SpliceSeq and TCGA aggregate expression datasets and identified differential alternative RNA splicing and transcription events (ARS/T) in cancers between self-identified African American (AA) and White (W) patients. We found that retained intron events were enriched among race-related ARS/T. In addition, on average, 12% of the most highly ranked race-related ARS/T overlapped between any two analyzed cancers. Moreover, the genes undergoing race-related ARS/T functioned in cancer-promoting pathways, and a number of race-related ARS/T were associated with patient survival. We built a web-application, CanSplice, to mine genomic datasets by self-identified race. The race-related targets have the potential to aid in the development of new biomarkers and therapeutics to mitigate cancer disparity.

HIV-1: To Splice or Not to Splice, That Is the Question

The transcription of the HIV-1 provirus results in only one type of transcript-full length genomic RNA. To make the mRNA transcripts for the accessory proteins Tat and Rev, the genomic RNA must completely splice. The mRNA transcripts for Vif, Vpr, and Env must undergo splicing but not completely. Genomic RNA (which also functions as mRNA for the Gag and Gag/Pro/Pol precursor polyproteins) must not splice at all. HIV-1 can tolerate a surprising range in the relative abundance of individual transcript types, and a surprising amount of aberrant and even odd splicing; however, it must not over-splice, which results in the loss of full-length genomic RNA and has a dramatic fitness cost. Cells typically do not tolerate unspliced/incompletely spliced transcripts, so HIV-1 must circumvent this cell policing mechanism to allow some splicing while suppressing most. Splicing is controlled by RNA secondary structure, cis-acting regulatory sequences which bind splicing factors, and the viral protein Rev. There is still much work to be done to clarify the combinatorial effects of these splicing regulators. These control mechanisms represent attractive targets to induce over-splicing as an antiviral strategy. Finally, splicing has been implicated in latency, but to date there is little supporting evidence for such a mechanism. In this review we apply what is known of cellular splicing to understand splicing in HIV-1, and present data from our newer and more sensitive deep sequencing assays quantifying the different HIV-1 transcript types.

RNA splicing and aggregate gene expression differences in lung squamous cell carcinoma between patients of West African and European ancestry

OBJECTIVES

Despite disparities in lung cancer incidence and mortality, the molecular landscape of lung cancer in patients of African ancestry remains underexplored, and race-related differences in RNA splicing remain unexplored.

MATERIALS AND METHODS

We identified differentially spliced genes (DSGs) and differentially expressed genes (DEGs) in biobanked lung squamous cell carcinoma (LUSC) between patients of West African and European ancestry, using ancestral genotyping and Affymetrix Clariom D array. DSGs and DEGs were validated independently using the National Cancer Institute Genomic Data Commons. Associated biological processes, overlapping canonical pathways, enriched gene sets, and cancer relevance were identified using Gene Ontology Consortium, Ingenuity Pathway Analysis, Gene Set Enrichment Analysis, and CancerMine, respectively. Association with LUSC survival was conducted using The Cancer Genome Atlas.

RESULTS

4,829 DSGs and 267 DEGs were identified, including novel targets in NSCLC as well as genes identified previously to have relevance to NSCLC. RNA splicing events within 3 DSGs as well as 1 DEG were validated in the independent cohort. 853 DSGs and 29 DEGs have been implicated as potential drivers, oncogenes and/or tumor suppressor genes. Biological processes enriched among DSGs and DEGs included metabolic process, biological regulation, and multicellular organismal process and, among DSGs, ion transport. Overlapping canonical pathways among DSGs included neuronal signaling pathways and, among DEGs, cell metabolism involving biosynthesis. Gene sets enriched among DSGs included KRAS Signaling, UV Response, E2 F Targets, Glycolysis, and Coagulation. 355 RNA splicing events within DSGs and 18 DEGs show potential association with LUSC patient survival.

CONCLUSION

These DSGs and DEGs, which show potential biological and clinical relevance, could have the ability to drive novel biomarker and therapeutic development to mitigate LUSC disparities.

Systematic Analyses of the Differentially Expressed Alternative Splicing Events in Gastric Cancer and Its Clinical Significance

Accumulation of evidence has indicated a close relationship between alternative splicing (AS) and gastric cancer (GC), whereas systematic analyses of the differentially expressed AS events (DEAS) between GC and normal tissues are lacking. RNA-Seq data and the corresponding clinical information were downloaded from TCGA GC cohort. The percent spliced-in (PSI) value calculated in the GC tissues and normal tissues was employed to quantify the DEAS. Further, survival-associated DEAS and DEAS signatures were identified by univariate and multivariate cox regression analyses. To evaluate the association between DEAS and patients' clinical features, Kaplan-Meier analysis, receiver operator characteristic (ROC) curve, Cox proportional regression and nomograms incorporating the DEAS signatures were performed. DEAS and their splicing networks were finally analyzed by bioinformatics methods. In addition, we use the method of random grouping to divide the samples into the training group and the test group. The final results of the two groups are consistent. After strict filtering, a total of 44,935 AS events were identified, among which 11,141 DEAS were preliminarily screened from 5032 genes. A total of 454 DEAS was associated with OS, and 872 DEAS were associated with DFS. The final prognostic signatures were constructed from the survival-associated DEAS with an area under the receiver operating characteristic (ROC) curve (AUC) greater than 0.6. Only ES in ABI1 was simultaneously associated with OS and DFS. Finally, we identified the splicing correlation network between the prognostic splicing factors (SF) and DEAS in GC. Our study provided a systematic portrait of survival-associated DEAS in GC and uncovered splicing networks that are valuable in deciphering the underlying mechanisms of AS in GC.

RNA splicing alteration in the response to platinum chemotherapy in ovarian cancer: A possible biomarker and therapeutic target

Since its discovery, alternative splicing has been recognized as a powerful way for a cell to amplify the genetic information and for a living organism to adapt, evolve, and survive. We now know that a very high number of genes are regulated by alternative splicing and that alterations of splicing have been observed in different types of human diseases, including cancer. Here, we review the accumulating knowledge that links the regulation of alternative splicing to the response to chemotherapy, focusing our attention on ovarian cancer and platinum-based treatments. Moreover, we discuss how expanding information could be exploited to identify new possible biomarkers of platinum response, to better select patients, and/or to design new therapies able to overcome platinum resistance.

Key Parameters of Tumor Epitope Immunogenicity Revealed Through a Consortium Approach Improve Neoantigen Prediction

Many approaches to identify therapeutically relevant neoantigens couple tumor sequencing with bioinformatic algorithms and inferred rules of tumor epitope immunogenicity. However, there are no reference data to compare these approaches, and the parameters governing tumor epitope immunogenicity remain unclear. Here, we assembled a global consortium wherein each participant predicted immunogenic epitopes from shared tumor sequencing data. 608 epitopes were subsequently assessed for T cell binding in patient-matched samples. By integrating peptide features associated with presentation and recognition, we developed a model of tumor epitope immunogenicity that filtered out 98% of non-immunogenic peptides with a precision above 0.70. Pipelines prioritizing model features had superior performance, and pipeline alterations leveraging them improved prediction performance. These findings were validated in an independent cohort of 310 epitopes prioritized from tumor sequencing data and assessed for T cell binding. This data resource enables identification of parameters underlying effective anti-tumor immunity and is available to the research community.

Splice variants of RAS-translational significance

One of the mechanisms potentially explaining the discrepancy between the number of human genes and the functional complexity of organisms is generating alternative splice variants, an attribute of the vast majority of multi-exon genes. Members of the RAS family, such as NRAS, KRAS and HRAS, all of which are of significant importance in cancer biology, are no exception. The structural and functional differences of these splice variants, particularly if they contain the canonical (and therefore routinely targeted for diagnostic purposes) hot spot mutations, pose a significant challenge for targeted therapies. We must therefore consider whether these alternative splice variants constitute a minor component as originally thought and how therapies targeting the canonical isoforms affect these alternative splice variants and their overall functions.

Mutations in transmembrane proteins: diseases, evolutionary insights, prediction and comparison with globular proteins

Membrane proteins are unique in that they interact with lipid bilayers, making them indispensable for transporting molecules and relaying signals between and across cells. Due to the significance of the protein's functions, mutations often have profound effects on the fitness of the host. This is apparent both from experimental studies, which implicated numerous missense variants in diseases, as well as from evolutionary signals that allow elucidating the physicochemical constraints that intermembrane and aqueous environments bring. In this review, we report on the current state of knowledge acquired on missense variants (referred to as to single amino acid variants) affecting membrane proteins as well as the insights that can be extrapolated from data already available. This includes an overview of the annotations for membrane protein variants that have been collated within databases dedicated to the topic, bioinformatics approaches that leverage evolutionary information in order to shed light on previously uncharacterized membrane protein structures or interaction interfaces, tools for predicting the effects of mutations tailored specifically towards the characteristics of membrane proteins as well as two clinically relevant case studies explaining the implications of mutated membrane proteins in cancer and cardiomyopathy.

Hypoxia-induced alternative splicing: the 11th Hallmark of Cancer

Hypoxia-induced alternative splicing is a potent driving force in tumour pathogenesis and progression. In this review, we update currents concepts of hypoxia-induced alternative splicing and how it influences tumour biology. Following brief descriptions of tumour-associated hypoxia and the pre-mRNA splicing process, we review the many ways hypoxia regulates alternative splicing and how hypoxia-induced alternative splicing impacts each individual hallmark of cancer. Hypoxia-induced alternative splicing integrates chemical and cellular tumour microenvironments, underpins continuous adaptation of the tumour cellular microenvironment responsible for metastatic progression and plays clear roles in oncogene activation and autonomous tumour growth, tumor suppressor inactivation, tumour cell immortalization, angiogenesis, tumour cell evasion of programmed cell death and the anti-tumour immune response, a tumour-promoting inflammatory response, adaptive metabolic re-programming, epithelial to mesenchymal transition, invasion and genetic instability, all of which combine to promote metastatic disease. The impressive number of hypoxia-induced alternative spliced protein isoforms that characterize tumour progression, classifies hypoxia-induced alternative splicing as the 11th hallmark of cancer, and offers a fertile source of potential diagnostic/prognostic markers and therapeutic targets.

Accessing a New Dimension in TP53 Biology: Multiplex Long Amplicon Digital PCR to Specifically Detect and Quantitate Individual TP53 Transcripts

TP53, the most commonly-mutated gene in cancer, undergoes complex alternative splicing. Different TP53 transcripts play different biological roles, both in normal function and in the progression of diseases such as cancer. The study of TP53’s alternative RNA splice forms and their use as clinical biomarkers has been hampered by limited specificity and quantitative accuracy of current methods. TP53 RNA splice variants differ at both 5’ and 3’ ends, but because they have a common central region of 618 bp, the individual TP53 transcripts are impossible to specifically detect and precisely quantitate using standard PCR-based methods or short-read RNA sequencing. Therefore, we devised multiplex probe-based long amplicon droplet digital PCR (ddPCR) assays, which for the first time allow precise end-to-end quantitation of the seven major TP53 transcripts, with amplicons ranging from 0.85 to 1.85 kb. Multiple modifications to standard ddPCR assay procedures were required to enable specific co-amplification of these long transcripts and to overcome issues with secondary structure. Using these assays, we show that several TP53 transcripts are co-expressed in breast cancers, and illustrate the potential for this method to identify novel TP53 transcripts in tumour cells. This capability will facilitate a new level of biological and clinical understanding of the alternatively-spliced TP53 isoforms.

Prognostic potential of PRPF3 in hepatocellular carcinoma

pre-mRNA processing factor 3 (PRPF3) is an RNA binding protein in a core component of the exon junction complex. Abnormal PRPF3 expression is potentially associated with carcinogenesis. However, the biological role of PRPF3 in hepatocellular carcinoma (HCC) remains to be determined. We analyzed PRPF3 expression via multiple gene expression databases and identified its genetic alterations and functional networks using cBioPortal. Co-expressed genes with PRPF3 and its regulators were identified using LinkedOmics. The correlations between PRPF3 and cancer immune infiltrates were investigated via Tumor Immune Estimation Resource (TIMER). PRPF3 was found up-regulated with amplification in tumor tissues in multiple HCC cohorts. High PRPF3 expression was associated with poorer overall survival (OS) and disease-free survival (DFS). Functional network analysis suggested that PRPF3 regulates spliceosome, DNA replication, and cell cycle signaling via pathways involving several cancer-related kinases and E2F family. Notably, PRPF3 expression was positively correlated with infiltrating levels of CD4+ T and CD8+ T cells, macrophages, neutrophils, and dendritic cells. PRPF3 expression showed strong correlations with diverse immune marker sets in HCC. These findings suggest that PRPF3 is correlated with prognosis and immune infiltrating in HCC, laying a foundation for further study of the immune regulatory role of PRPF3 in HCC.