Aberrant and Alternative Splicing in Cancer

Molecular targets and mechanisms of different aberrant alternative splicing in metastatic liver cancer

Metastasis remains a major challenge in the successful management of malignant diseases. The liver is a major site of metastatic disease and a leading cause of death from gastrointestinal malignancies such as colon, stomach, and pancreatic cancers, as well as melanoma, breast cancer, and sarcoma. As an important factor that influences the development of metastatic liver cancer, alternative splicing drives the diversity of RNA transcripts and protein subtypes, which may provide potential to broaden the target space. In particular, the dysfunction of splicing factors and abnormal expression of splicing variants are associated with the occurrence, progression, aggressiveness, and drug resistance of cancers caused by the selective splicing of specific genes. This review is the first to provide a detailed summary of the normal splicing process and alterations that occur during metastatic liver cancer. It will cover the role of alternative splicing in the mechanisms of metastatic liver cancer by examining splicing factor changes, abnormal splicing, and the contribution of hypoxia to these changes during metastasis.

RNA-binding proteins potentially regulate the alternative splicing of apoptotic genes during knee osteoarthritis progression

BACKGROUND

Alternative splicing (AS) is a principal mode of genetic regulation and one of the most widely used mechanisms to generate structurally and functionally distinct mRNA and protein variants. Dysregulation of AS may result in aberrant transcription and protein products, leading to the emergence of human diseases. Although considered important for regulating gene expression, genome-wide AS dysregulation, underlying mechanisms, and clinical relevance in knee osteoarthritis (OA) remain unelucidated. Therefore, in this study, we elucidated and validated AS events and their regulatory mechanisms during OA progression.

RESULTS

In this study, we identified differentially expressed genes between human OA and healthy meniscus samples. Among them, the OA-associated genes were primarily enriched in biological pathways such as extracellular matrix organization and ossification. The predominant OA-associated regulated AS (RAS) events were found to be involved in apoptosis during OA development. The expression of the apoptosis-related gene BCL2L13, XAF1, and NF2 were significantly different between OA and healthy meniscus samples. The construction of a covariation network of RNA-binding proteins (RBPs) and RAS genes revealed that differentially expressed RBP genes LAMA2 and CUL4B may regulate the apoptotic genes XAF1 and BCL2L13 to undergo AS events during OA progression. Finally, RT-qPCR revealed that CUL4B expression was significantly higher in OA meniscus samples than in normal controls and that the AS ratio of XAF1 was significantly different between control and OA samples; these findings were consistent with their expected expression and regulatory relationships.

CONCLUSIONS

Differentially expressed RBPs may regulate the AS of apoptotic genes during knee OA progression. XAF1 and its regulator, CUL4B, may serve as novel biomarkers and potential therapeutic targets for this disease.

The nexus of long noncoding RNAs, splicing factors, alternative splicing and their modulations

The process of alternative splicing (AS) is widely deregulated in a variety of cancers. Splicing is dependent upon splicing factors. Recently, several long noncoding RNAs (lncRNAs) have been shown to regulate AS by directly/indirectly interacting with splicing factors. This review focuses on the regulation of AS by lncRNAs through their interaction with splicing factors. AS mis-regulation caused by either mutation in splicing factors or deregulated expression of splicing factors and lncRNAs has been shown to be involved in cancer development and progression, making aberrant splicing, splicing factors and lncRNA suitable targets for cancer therapy. This review also addresses some of the current approaches used to target AS, splicing factors and lncRNAs. Finally, we discuss research challenges, some of the unanswered questions in the field and provide recommendations to advance understanding of the nexus of lncRNAs, AS and splicing factors in cancer.

Human histone H1 variants impact splicing outcome by controlling RNA polymerase II elongation

Histones shape chromatin structure and the epigenetic landscape. H1, the most diverse histone in the human genome, has 11 variants. Due to the high structural similarity between the H1s, their unique functions in transferring information from the chromatin to mRNA-processing machineries have remained elusive. Here, we generated human cell lines lacking up to five H1 subtypes, allowing us to characterize the genomic binding profiles of six H1 variants. Most H1s bind to specific sites, and binding depends on multiple factors, including GC content. The highly expressed H1.2 has a high affinity for exons, whereas H1.3 binds intronic sequences. H1s are major splicing regulators, especially of exon skipping and intron retention events, through their effects on the elongation of RNA polymerase II (RNAPII). Thus, H1 variants determine splicing fate by modulating RNAPII elongation.

RNA splicing alterations in lung cancer pathogenesis and therapy

RNA splicing alterations are widespread and play critical roles in cancer pathogenesis and therapy. Lung cancer is highly heterogeneous and causes the most cancer-related deaths worldwide. Large-scale multi-omics studies have not only characterized the mutational landscapes but also discovered a plethora of transcriptional and post-transcriptional changes in lung cancer. Such resources have greatly facilitated the development of new diagnostic markers and therapeutic options over the past two decades. Intriguingly, altered RNA splicing has emerged as an important molecular feature and therapeutic target of lung cancer. In this review, we provide a brief overview of splicing dysregulation in lung cancer and summarize the recent progress on key splicing events and splicing factors that contribute to lung cancer pathogenesis. Moreover, we describe the general strategies targeting splicing alterations in lung cancer and highlight the potential of combining splicing modulation with currently approved therapies to combat this deadly disease. This review provides new mechanistic and therapeutic insights into splicing dysregulation in cancer.

Introns and Their Therapeutic Applications in Biomedical Researches

Context

Although for a long time, it was thought that intervening sequences (introns) were junk DNA without any function, their critical roles and the underlying molecular mechanisms in genome regulation have only recently come to light. Introns not only carry information for splicing, but they also play many supportive roles in gene regulation at different levels. They are supposed to function as useful tools in various biological processes, particularly in the diagnosis and treatment of diseases. Introns can contribute to numerous biological processes, including gene silencing, gene imprinting, transcription, mRNA metabolism, mRNA nuclear export, mRNA localization, mRNA surveillance, RNA editing, NMD, translation, protein stability, ribosome biogenesis, cell growth, embryonic development, apoptosis, molecular evolution, genome expansion, and proteome diversity through various mechanisms.

Evidence Acquisition

In order to fulfill the objectives of this study, the following databases were searched: Medline, Scopus, Web of Science, EBSCO, Open Access Journals, and Google Scholar. Only articles published in English were included.

Results & Conclusions

The intervening sequences of eukaryotic genes have critical functions in genome regulation, as well as in molecular evolution. Here, we summarize recent advances in our understanding of how introns influence genome regulation, as well as their effects on molecular evolution. Moreover, therapeutic strategies based on intron sequences are discussed. According to the obtained results, a thorough understanding of intron functional mechanisms could lead to new opportunities in disease diagnosis and therapies, as well as in biotechnology applications.

Target-capture full-length double-stranded cDNA long-read sequencing through Nanopore revealed novel intron retention in patient with tuberous sclerosis complex

Tuberous sclerosis complex (TSC) is a relatively common autosomal dominant disorder characterized by multiple dysplastic organ lesions and neuropsychiatric symptoms caused by loss-of-function mutation of either TSC1 or TSC2. The genetic diagnosis of inherited diseases, including TSC, in the clinical field is widespread using next-generation sequencing. The mutations in protein-coding exon tend to be verified because mutations directly cause abnormal protein. However, it is relatively difficult to verify mutations in the intron region because it is required to investigate whether the intron mutations affect the abnormal splicing of transcripts. In this study, we developed a target-capture full-length double-stranded cDNA sequencing method using Nanopore long-read sequencer (Nanopore long-read target sequencing). This method revealed the occurrence of intron mutation in the TSC2 gene and found that the intron mutation produces novel intron retention splicing transcripts that generate truncated proteins. The protein-coding transcripts were decreased due to the expression of the novel intron retention transcripts, which caused TSC in patients with the intron mutation. Our results indicate that Nanopore long-read target sequencing is useful for the detection of mutations and confers information on the full-length alternative splicing of transcripts for genetic diagnosis.

A patient-derived organoid-based study identified an ASO targeting SNORD14E for endometrial cancer through reducing aberrant FOXM1 Expression and β-catenin nuclear accumulation

BACKGROUND

Most of the endometrial cancer (EC) patients are diagnosis in early stage with a good prognosis while the patients with locally advanced recurrent or metastatic result in a poor prognosis. Adjuvant therapy could benefit the prognosis of patients with high-risk factors. Unfortunately, the molecular classification of great prognostic value has not yet reached an agreement and need to be further refined. The present study aims to identify new targets that have prognostic value in EC based on the method of EC patient-derived organ-like organs (PDOs), and further investigate their efficacy and mechanism.

METHODS

The Cancer Genome Atlas (TCGA) database was used to determine SNORD14E expression. The effects of SNORD14E were investigated using CCK8, Transwell, wound-healing assays, and a xenograft model experiment; apoptosis was measured by flow cytometry. Antisense oligonucleotide (ASO) targeting SNORD14E was designed and patient-derived organoids (PDO) models in EC patients was established. A xenograft mouse and PDO model were employed to evaluate the effects of ASO targeting SNORD14E. RNA-seq, Nm-seq, and RNA immunoprecipitation (RIP) experiments were employed to confirm the alternative splicing (AS) and modification induced by SNORD14E. A minigene reporter gene assay was conducted to confirm AS and splicing factors on a variable exon. Actinomycin-d (Act-D) and Reverse Transcription at Low deoxy-ribonucleoside triphosphate concentrations followed by PCR (RTL-P) were utilized to confirm the effects of 2'-O methylation modification on FOXM1.

RESULTS

We found that SNORD14E was overexpressed in EC tissues and patients with high expressed SNORD14E were distributed in the TCGA biomolecular classification subgroups without difference. Further, SNORD14E could reduce disease-free survival (DFS) and recurrence free survival (RFS) of EC patients. SNORD14E promoted proliferation, migration, and invasion and inhibited the apoptosis of EC cells in vitro. ASOs targeting SNORD14E inhibited cell proliferation, migration, invasion while promoted cell apoptosis. ASOs targeting SNORD14E inhibited tumor growth in the xenograft mouse model. TCGA-UCEC database showed that the proportion of patients with high expression of SNORD14E in middle-high risk and high-risk patients recommended by EMSO-ESGO-ESTRO guidelines for adjuvant therapy is more than 50%. Next, we enrolled 8 cases of high-risk and high-risk EC patients according to EMSO-ESGO-ESTRO guidelines and successfully constructed EC-PDOs. ASOs targeting SNORD14E inhibited the EC-PDO growth. Mechanistically, SNORD14E could recognize the mRNA of FOXM1 and recruit SRSF1 to promote the shearing of the variable exon VIIa of FOXM1, resulting in the overexpression of the FOXM1 malignant subtypes FOXM1b and FOXM1c. In addition, SNORD14E modified FOXM1 mRNA with 2`-O-methylation, which prolonged the half-life of FOXM1 mRNA. The nucleus accumulation of β-catenin caused by aberrant FOXM1 expression led to EC progression.

CONCLUSIONS

ASO targeting SNORD14E can be an effective treatment for EC.

Moonlighting functions of the ubiquitin-like protein, Hub1/UBL-5

The faithful splicing of pre-mRNA is critical for accurate gene expression. Dysregulation of pre-mRNA splicing has been associated with several human diseases including cancer. The ubiquitin-like protein Hub1/UBL5 binds to the substrates non-covalently and promotes pre-mRNA splicing. Additionally, UBL5 promotes the common fragile sites stability and the Fanconi anemia pathway of DNA damage repair. These functions strongly suggests that UBL5 could potentially be implicated in cancer. Therefore, we analyzed the UBL5 expression in TCGA tumor sample datasets and observed the differences between tumor and normal tissues among different tumor subtypes. We have noticed the alteration frequency of UBL5 in TCGA tumor samples. Altogether, this review summarizes the UBL5 functions and discusses its putative role in tumorigenesis.

Proteogenomic data and resources for pan-cancer analysis

The National Cancer Institute's Clinical Proteomic Tumor Analysis Consortium (CPTAC) investigates tumors from a proteogenomic perspective, creating rich multi-omics datasets connecting genomic aberrations to cancer phenotypes. To facilitate pan-cancer investigations, we have generated harmonized genomic, transcriptomic, proteomic, and clinical data for >1000 tumors in 10 cohorts to create a cohesive and powerful dataset for scientific discovery. We outline efforts by the CPTAC pan-cancer working group in data harmonization, data dissemination, and computational resources for aiding biological discoveries. We also discuss challenges for multi-omics data integration and analysis, specifically the unique challenges of working with both nucleotide sequencing and mass spectrometry proteomics data.

Granulocyte-Macrophage Colony-Stimulating Factor Influence on Soluble and Membrane-Bound ICOS in Combination with Immune Checkpoint Blockade

With the successful development of immune checkpoint blockade, there remains the continued need to improve efficacy and decrease toxicities. The addition of granulocyte-macrophage colony-stimulating factor (GM-CSF) to ipilimumab has previously demonstrated both an improvement in efficacy and decrease in the incidence of high-grade adverse events. ICOS+CD4+ or ICOS+CD8+ peripheral blood T cells are significantly greater in the patients treated with ipilimumab plus GM-CSF than in the patients treated with ipilimumab alone. To better understand the effects of GM-CSF on inducible T-cell costimulator (ICOS) and clinical outcomes, the relative roles of identified soluble ICOS and membrane-bound ICOS were evaluated. The ICOS splice variant was secreted and found to have immunologic suppressive effects. Changes in soluble ICOS splice variant levels in treated patients correlated with clinical outcomes. GM-CSF enhanced membrane-bound ICOS in an IL12-dependent manner but did not increase soluble ICOS levels. Whereas soluble ICOS plays a role in immune suppression, GM-CSF efficacy involves increasing membrane-bound ICOS and induction of dendritic cell development. Thus, soluble ICOS splice variants may be used as a biomarker for GM-CSF and immune checkpoint blockade-based therapies.

Tandem enrichment of serum exosomes and exosomal RNA with titanium dioxide

Exosomes have great potential as biomarker carriers for disease diagnosis and prognosis. In recent years, exosomal RNA (exoRNA) has become a promising candidate for the early diagnosis and prognosis of cancers, and its pathophysiological roles in various diseases have been revealed. For example, exosome-derived mRNAs, miRNAs, circRNAs, and lncRNAs function as signalling molecules to regulate tumour growth, angiogenesis, invasion, metastasis, and the response to chemotherapy. However, the isolation of exosomes and exoRNA with high quality and purity remains challenging due to the relatively small size of exosomes and the limited amount of RNA in exosomes. In this work, we developed a novel tandem enrichment method to isolate exoRNA from serum based on the specific interaction between titanium dioxide (TiO₂) and the phosphate groups on the lipid bilayer of exosomes and of the exoRNA. TiO₂-based RNA isolation was first demonstrated and optimized in HeLa cells. A total of 130.9 ± 8.34 µg of RNA was rapidly enriched from approximately 5 × 10⁶ HeLa cells within 10 min. This was a 41.5% higher yield than that using a commercial Ultrapure RNA Kit. TiO₂-based tandem enrichment of exoRNA was then performed using human serum, obtaining 64.53±3.41 ng of exoRNA from 500 µL of human serum within 30 min. A total of 2,137,902 reads, including seven types of exoRNAs, were identified from the exosomes. This method is compatible with various downstream RNA processing techniques and does not use toxic or irritating reagents, such as phenol or chloroform, providing a simple, economical, rapid, and safe approach for exoRNA extraction from biological samples.

Integrated analysis of genomic and transcriptomic data for the discovery of splice-associated variants in cancer

Somatic mutations within non-coding regions and even exons may have unidentified regulatory consequences that are often overlooked in analysis workflows. Here we present RegTools ( www.regtools.org ), a computationally efficient, free, and open-source software package designed to integrate somatic variants from genomic data with splice junctions from bulk or single cell transcriptomic data to identify variants that may cause aberrant splicing. We apply RegTools to over 9000 tumor samples with both tumor DNA and RNA sequence data. RegTools discovers 235,778 events where a splice-associated variant significantly increases the splicing of a particular junction, across 158,200 unique variants and 131,212 unique junctions. To characterize these somatic variants and their associated splice isoforms, we annotate them with the Variant Effect Predictor, SpliceAI, and Genotype-Tissue Expression junction counts and compare our results to other tools that integrate genomic and transcriptomic data. While many events are corroborated by the aforementioned tools, the flexibility of RegTools also allows us to identify splice-associated variants in known cancer drivers, such as TP53, CDKN2A, and B2M, and other genes.

CD22 Exon 12 Deletion as an Independent Predictor of Poor Treatment Outcomes in B-ALL

We previously reported a splicing defect (CD22ΔE12) associated with the deletion of exon 12 of the inhibitory co-receptor CD22 (Siglec-2) in leukemia cells from patients with CD19+ B-precursor acute lymphoblastic leukemia (B-ALL). CD22ΔE12 causes a truncating frameshift mutation and yields a dysfunctional CD22 protein that lacks most of the cytoplasmic domain required for its inhibitory function, and it is associated with aggressive in vivo growth of human B-ALL cells in mouse xenograft models. Although CD22ΔE12 with selective reduction of CD22 exon 12 (CD22E12) levels was detected in a high percentage of newly diagnosed as well as relapsed B-ALL patients, its clinical significance remains unknown. We hypothesized that B-ALL patients with very low levels of wildtype CD22 would exhibit a more aggressive disease with a worse prognosis because the missing inhibitory function of the truncated CD22 molecules could not be adequately compensated by competing wildtype CD22. Here, we demonstrate that newly diagnosed B-ALL patients with very low levels of residual wildtype CD22 ("CD22E12low"), as measured by RNAseq-based CD22E12 mRNA levels, have significantly worse leukemia-free survival (LFS) as well as overall survival (OS) than other B-ALL patients. CD22E12low status was identified as a poor prognostic indicator in both univariate and multivariate Cox proportional hazards models. CD22E12low status at presentation shows clinical potential as a poor prognostic biomarker that may guide the early allocation of risk-adjusted, patient-tailored treatment regimens and refine risk classification in high-risk B-ALL.

Elevated expression of ISY1, APOA-1, SYNE1, MTG1, and MMP10 at HCC initiation: HCC specific protein network involving interactions of key regulators of lipid metabolism, EGFR signaling, MAPK, and splicing pathways

Identification of molecular regulators of hepatocellular carcinoma (HCC) initiation and progression is not well understood. We chemically induced HCC in male Wistar rats by administration of diethyl nitrosamine (DEN) and 2-acetylaminofluorene (2-AFF). Using 2D-electrophoresis and MALDI-TOF-MS/MS analyses, we characterized differentially expressed proteins in liver tissues at early stage of HCC progression. Using RT-PCR analysis, we quantified the mRNA expression of the characterized proteins and validated the transcript expression with tumor tissues of clinically confirmed HCC patients. Using bioinformatic tools, we analyzed a network among the introduced proteins that identified their interacting partners and analyzed the molecular mechanisms associated with signaling pathways during HCC progression. We characterized a protein, namely, pre-mRNA splicing factor 1 homolog (ISY1), which is upregulated at both transcriptome and proteome levels at HCC initiation, progression, and tumor stages. We analyzed the interacting partners of ISY1, namely, APOA-1, SYNE1, MMP10, and MTG1. Real-time PCR analysis confirmed elevated expression of APOA-1 mRNA at HCC initiation, progression, and tumor stages in animals undergoing tumorigenesis. The mRNA expression of the interacting partners was validated with tumor tissues of clinically confirmed liver cancer patients; the analysis revealed significant elevation in expression of transcripts. The transcriptome and proteome analyses complement each other and dysregulation in mRNA and protein expression of these regulators may play critical role in HCC initiation and progression.

RNA-Seq Analysis of Extradomain A and Extradomain B Fibronectin as Extracellular Matrix Markers for Cancer

Alternatively spliced forms of fibronectin, called oncofetal fibronectin, are aberrantly expressed in cancer, with little to no expression in normal tissue, making them attractive biomarkers to exploit for tumor-targeted therapeutics and diagnostics. While prior studies have explored oncofetal fibronectin expression in limited cancer types and limited sample sizes, no studies have performed a large-scale pan-cancer analysis in the context of clinical diagnostics and prognostics to posit the utility of these biomarkers across multiple cancer types. In this study, RNA-Seq data sourced from the UCSC Toil Recompute project were extracted and analyzed to determine the correlation between the expression of oncofetal fibronectin, including extradomain A and extradomain B fibronectin, and patient diagnosis and prognosis. We determined that oncofetal fibronectin is significantly overexpressed in most cancer types relative to corresponding normal tissues. In addition, strong correlations exist between increasing oncofetal fibronectin expression levels and tumor stage, lymph node activity, and histological grade at the time of diagnosis. Furthermore, oncofetal fibronectin expression is shown to be significantly associated with overall patient survival within a 10-year window. Thus, the results presented in this study suggest oncofetal fibronectin as a commonly upregulated biomarker in cancer with the potential to be used for tumor-selective diagnosis and treatment applications.

Alternative splicing in multiple myeloma is associated with the non-homologous end joining pathway

Alternative splicing plays a pivotal role in tumorigenesis and proliferation. However, its pattern and pathogenic role has not been systematically analyzed in multiple myeloma or its subtypes. Alternative splicing profiles for 598 newly diagnosed myeloma patients with comprehensive genomic annotation identified primary translocations, 1q amplification, and DIS3 events to have more differentially spliced events than those without. Splicing levels were correlated with expression of splicing factors. Moreover, the non-homologous end joining pathway was an independent factor that was highly associated with splicing frequency as well as an increased number of structural variants. We therefore identify an axis of high-risk disease encompassing expression of the non-homologous end joining pathway, increase structural variants, and increased alternative splicing that are linked together. This indicates a joint pathogenic role for DNA damage response and alternative RNA processing in myeloma.

Genome-wide identification of alternative splicing and splicing regulated in immune infiltration in osteosarcoma patients

Background: Osteosarcoma typically occurs in adolescents, and the survival rate of patients with metastatic and recurrent osteosarcoma remains low. Abnormal regulation of alternative splicing is associated with the development of osteosarcoma. However, there is no genome-wide analysis of the function and regulatory mechanisms of aberrant alternative splicing associated with osteosarcoma. Methods: Published transcriptome data on osteosarcoma (GSE126209) derived from osteosarcoma patient tissue were downloaded. Gene expression profiling by high-throughput sequencing was performed on 9 normal samples and 10 tumor samples for genome-wide identification of osteosarcoma-related alternative splicing events. The potential function of osteosarcoma-associated alternative splicing events was examined by immune infiltration and correlation analysis. Regulation of aberrantly expressed RNA-binding proteins (RBPs) related to alternative splicing in osteosarcoma was clarified by co-expression analysis. Results: A total of 63 alternative splicing events, which are highly credible and dominant, were identified. GO enrichment analysis indicated that alternative splicing may be closely related to the immune response process. Immune infiltration analysis showed significant changes in the percentages of CD8 T cells, resting memory CD4 T cells, activated memory CD4 T cells, monocytes, resting dendritic cells, and activated mast cells in tumors compared to normal tissues, indicating the involvement of these immune cell types in the occurrence of osteosarcoma. Moreover, the analysis identified alternative splicing events that were co-altered with resting memory CD4 T cells, resting dendritic cells, and activated mast cells, events that may be associated with regulation of the osteosarcoma immune microenvironment. In addition, a co-regulatory network (RBP-RAS-immune) of osteosarcoma-associated RBPs with aberrant alternative splicing and altered immune cells was established. These RBPs include NOP58, FAM120C, DYNC1H1, TRAP1, and LMNA, which may serve as molecular targets for osteosarcoma immune regulation. Conclusion: These findings allow us to further understand the causes of osteosarcoma development and provide a new research direction for osteosarcoma immunotherapy or targeted therapy.

The key role of microRNA-766 in the cancer development

Cancer is caused by defects in coding and non-coding RNAs. In addition, duplicated biological pathways diminish the efficacy of mono target cancer drugs. MicroRNAs (miRNAs) are short, endogenous, non-coding RNAs that regulate many target genes and play a crucial role in physiological processes such as cell division, differentiation, cell cycle, proliferation, and apoptosis, which are frequently disrupted in diseases such as cancer. MiR-766, one of the most adaptable and highly conserved microRNAs, is notably overexpressed in several diseases, including malignant tumors. Variations in miR-766 expression are linked to various pathological and physiological processes. Additionally, miR-766 promotes therapeutic resistance pathways in various types of tumors. Here, we present and discuss evidence implicating miR-766 in the development of cancer and treatment resistance. In addition, we discuss the potential applications of miR-766 as a therapeutic cancer target, diagnostic biomarker, and prognostic indicator. This may shed light on the development of novel therapeutic strategies for cancer therapy.

PSSNet-An Accurate Super-Secondary Structure for Protein Segmentation

A super-secondary structure (SSS) is a spatially unique ensemble of secondary structural elements that determine the three-dimensional shape of a protein and its function, rendering SSSs attractive as folding cores. Understanding known types of SSSs is important for developing a deeper understanding of the mechanisms of protein folding. Here, we propose a universal PSSNet machine-learning method for SSS recognition and segmentation. For various types of SSS segmentation, this method uses key characteristics of SSS geometry, including the lengths of secondary structural elements and the distances between them, torsion angles, spatial positions of Cα atoms, and primary sequences. Using four types of SSSs (βαβ-unit, α-hairpin, β-hairpin, αα-corner), we showed that extensive SSS sets could be reliably selected from the Protein Data Bank and AlphaFold 2.0 database of protein structures.

Identification of tumor antigens and immune subtypes in head and neck squamous cell carcinoma for mRNA vaccine development

The mRNA vaccines have been considered effective for combating cancer. However, the core components of the mRNA vaccines against head and neck squamous cell carcinoma (HNSCC) and the effects remain unclear. Our study aims to identify effective antigens in HNSCC to develop mRNA vaccines for corresponding potential patients. Here, we analyzed alternative splicing and mutation of genes in TCGA-HNSCC samples and identified seven potential tumor antigens, including SREBF1, LUC7L3, LAMA5, PCGF3, HNRNPH1, KLC4, and OFD1, which were associated with nonsense-mediated mRNA decay factor expression, overall survival prognosis and the infiltration of antigen-presenting cells. Furthermore, to select suitable patients for vaccination, immune subtypes related to HNSCC were identified by consensus clustering analysis, and visualization of the HNSCC immune landscape was performed by graph-learning-based dimensionality reduction. To address the heterogeneity of the population that is suitable for vaccination, plot cell trajectory and WGCNA were also utilized. HNSCC patients were classified into three prognostically relevant immune subtypes (Cluster 1, Cluster 2, and Cluster 3) possessing different molecular and cellular characteristics, immune modulators, and mutation statuses. Cluster 1 had an immune-activated phenotype and was associated with better survival, while Cluster 2 and Cluster 3 were immunologically cold and linked to increased tumor mutation burden. Therefore, HNSCC patients with immune subtypes Cluster 2 and Cluster 3 are potentially suitable for mRNA vaccination. Moreover, the prognostic module hub genes screened seven genes, including IGKC, IGHV3-15, IGLV1-40, IGLV1-51, IGLC3, IGLC2, and CD79A, which could be potential biomarkers to predict prognosis and identify suitable patients for mRNA vaccines. Our findings provide a theoretical basis for further research and the development of anti-HNSCC mRNA vaccines and the selection of suitable patients for vaccination.

Intravascular NK/T-Cell Lymphoma: What We Know about This Diagnostically Challenging, Aggressive Disease

Intravascular lymphoma is a form of lymphoid malignancy characterized by neoplastic cells growing almost exclusively within the lumina of small- to medium-sized blood vessels. Most cases are of B-cell origin with rare cases of natural killer or T-cell lineage. Extranodal sites are affected, mainly the skin and central nervous system, although any organ may be involved. Intravascular NK/T-cell lymphoma deserves special attention because of its clinicopathologic features and the need for adequate immunophenotyping combined with clonality test for a proper diagnosis. Moreover, intravascular NK/T-cell lymphoma is strongly linked to Epstein-Barr virus (EBV), which is considered to play a role in tumorigenesis and to be responsible for the aggressive behavior of the disease. In this paper, we review the current knowledge on this rare lymphoma and, in particular, the most recent advances about its molecular landscape. The main distinguishing features with other EBV-related entities, such as extranodal NK/T-cell lymphoma, EBV-positive primary nodal T/NK-cell lymphoma, and aggressive NK-cell leukemia, are discussed to help pathologists obtain the correct diagnosis and consequently develop an adequate and prompt therapy response.

SERBP1 affects the apoptotic level by regulating the expression and alternative splicing of cellular and metabolic process genes in HeLa cells

Background

RNA-binding proteins (RBPs) have important roles in orchestrating posttranscriptional regulation and modulating many tumorigenesis events. SERBP1 has been recognized as an important regulator in multiple cancers, while it remains unclear whether SERBP1-regulated gene expression at the transcriptome-wide level is significantly correlated with tumorigenesis.

Methods

We overexpressed SERBP1 in HeLa cells and explored whether SERBP1 overexpression (SERBP1-OE) affects the proliferation and apoptosis of HeLa cells. We analyzed the transcriptome-wide gene expression changes and alternative splicing changes mediated by SERBP1-OE using the transcriptome sequencing method (RNA-seq). RT-qPCR was conducted to assay SERBP1-regulated alternative splicing.

Results

SERBP1-OE induced the apoptosis of HeLa cells. The downregulated genes were strongly enriched in the cell proliferation and apoptosis pathways according to the GO analysis, including FOS, FOSB, PAK6 and RAB26. The genes undergoing at least one SERBP1-regulated alternative splicing event were enriched in transcriptional regulation, suggesting a mechanism of the regulation of gene expression, and in pyruvate and fatty acid metabolic processes critical for tumorigenesis events. The SERBP1-regulated alternative splicing of ME3, LPIN3, CROT, PDP1, SLC27A1 and ALKBH7 was validated by RT-qPCR analysis.

Conclusions

We for the first time demonstrated the cellular function and molecular targets of SERBP1 in HeLa cells at transcriptional and post-transcriptional levels. The SERBP1-regulated gene expression and alternative splicing networks revealed by this study provide important information for exploring the functional roles and regulatory mechanisms of SERBP1 in cancer development and progression.

Comprehensive analysis of full-length transcripts reveals novel splicing abnormalities and oncogenic transcripts in liver cancer

Genes generate transcripts of various functions by alternative splicing. However, in most transcriptome studies, short-reads sequencing technologies (next-generation sequencers) have been used, leaving full-length transcripts unobserved directly. Although long-reads sequencing technologies would enable the sequencing of full-length transcripts, the data analysis is difficult. In this study, we developed an analysis pipeline named SPLICE and analyzed cDNA sequences from 42 pairs of hepatocellular carcinoma (HCC) and matched non-cancerous livers with an Oxford Nanopore sequencer. Our analysis detected 46,663 transcripts from the protein-coding genes in the HCCs and the matched non-cancerous livers, of which 5,366 (11.5%) were novel. A comparison of expression levels identified 9,933 differentially expressed transcripts (DETs) in 4,744 genes. Interestingly, 746 genes with DETs, including the LINE1-MET transcript, were not found by a gene-level analysis. We also found that fusion transcripts of transposable elements and hepatitis B virus (HBV) were overexpressed in HCCs. In vitro experiments on DETs showed that LINE1-MET and HBV-human transposable elements promoted cell growth. Furthermore, fusion gene detection showed novel recurrent fusion events that were not detected in the short-reads. These results suggest the efficiency of full-length transcriptome studies and the importance of splicing variants in carcinogenesis.

Genes generate transcripts of various functions by alternative splicing. However, in most transcriptome studies, short-reads sequencing technologies (next-generation sequencers) have been used, leaving full-length transcripts unobserved directly. In this study, we developed an analysis pipeline named SPLICE for long-read transcriptome sequencing and analyzed cDNA sequences from 42 pairs of hepatocellular carcinoma (HCC), and matched non-cancerous livers with an Oxford Nanopore sequencer. Our analysis detected 5,366 novel transcripts and 9,933 differentially expressed transcripts in 4,744 genes between HCCs and non-cancerous livers. An analysis of hepatitis B virus (HBV) transcripts showed that fusion transcripts of the HBV gene and human transposable elements were overexpressed in HBV-infected HCCs. We also identified fusion genes that were not found in the short-reads. These results suggest that long-reads sequencing technologies provide a fuller understanding of cancer transcripts and that our method contributes to the analysis of transcriptome sequences by such technologies.

The Role of Osteopontin in Tumor Progression Through Tumor-Associated Macrophages

Osteopontin (OPN) is a multifunctional phosphorylated protein. It is widely involved in solid tumor progression, such as intensification of macrophage recruitment, inhibition of T-cell activity, aggravation of tumor interstitial fibrosis, promotion of tumor metastasis, chemotherapy resistance, and angiogenesis. Most of these pathologies are affected by tumor-associated macrophages (TAMs), an important component of the tumor microenvironment (TME). TAMs have been extensively characterized, including their subsets, phenotypes, activation status, and functions, and are considered a promising therapeutic target for cancer treatment. This review focuses on the interaction between OPN and TAMs in mediating tumor progression. We discuss the strategies for targeting OPN and TAMs to treat cancer and factors that may affect the therapeutic outcomes of blocking OPN or depleting TAMs. We also discuss the role of cancer cell- vs. TAM-derived OPN in tumorigenesis, the mechanisms of how OPN affects TAM recruitment and polarization, and why OPN could mediate anti-tumor and pro-tumor effects, as well as previously reported discrepancies.

Identification of Alternative Splicing-Related Genes CYB561 and FOLH1 in the Tumor-Immune Microenvironment for Endometrial Cancer Based on TCGA Data Analysis

Background: Advanced and recurrent endometrial cancer EC remains controversial. Immunotherapy will play a landmark role in cancer treatment, and alternative splicing (AS) of messenger RNA (mRNA) may offer the potential of a broadened target space.

Methods: We downloaded the clinical information and mRNA expression profiles from The Cancer Genome Atlas (TCGA) database. Hub genes were extracted from 11 AS-related genes to analyze the correlation between clinical parameters and the tumor-immune microenvironment. We also analyzed the correlations between the copy numbers, gene expressions of hub genes, and immune cells. The correlation between the risk score and the six most important checkpoint genes was also investigated. The ESTIMATE algorithm was finally performed on each EC sample based on the high- and low-risk groups.

Results: The risk score was a reliable and stable independent risk predictor in the Uterine Corpus Endometrial Carcinoma (UCEC) cohort. CYB561|42921|AP and FOLH1|15817|ES were extracted. The expression of CYB561 and FOLH1 decreased gradually with the increased grade and International Federation of Gynecology and Obstetrics (FIGO) stage (p < 0.05). Gene copy number changes in CYB561 and FOLH1 led to the deletion number of myeloid DC cells and T cell CD8+. Low expression of both CYB561 and FOLH1 was associated with poor prognosis (p < 0.001). The checkpoint genes, CTLA-4 and PDCD1, exhibited a negative correlation with the risk score of AS in UCEC.

Conclusion: AS-related gene signatures were related to the immune-tumor microenvironment and prognosis. These outcomes were significant for studying EC’s immune-related mechanisms and exploring novel prognostic predictors and precise therapy methods.

Poly(A) capture full length cDNA sequencing improves the accuracy and detection ability of transcript quantification and alternative splicing events

The full-length double-strand cDNA sequencing, one of the RNA-Seq methods, is a powerful method used to investigate the transcriptome status of a gene of interest, such as its transcription level and alternative splicing variants. Furthermore, full-length double-strand cDNA sequencing has the advantage that it can create a library from a small amount of sample and the library can be applied to long-read sequencers in addition to short-read sequencers. Nevertheless, one of our previous studies indicated that the full-length double-strand cDNA sequencing yields non-specific genomic DNA amplification, affecting transcriptome analysis, such as transcript quantification and alternative splicing analysis. In this study, it was confirmed that it is possible to produce the RNA-Seq library from only genomic DNA and that the full-length double-strand cDNA sequencing of genomic DNA yielded non-specific genomic DNA amplification. To avoid non-specific genomic DNA amplification, two methods were examined, which are the DNase I-treated full-length double-strand cDNA sequencing and poly(A) capture full-length double-strand cDNA sequencing. Contrary to expectations, the non-specific genomic DNA amplification was increased and the number of the detected expressing genes was reduced in DNase I-treated full-length double-strand cDNA sequencing. On the other hand, in the poly(A) capture full-length double-strand cDNA sequencing, the non-specific genomic DNA amplification was significantly reduced, accordingly the accuracy and the number of detected expressing genes and splicing events were increased. The expression pattern and percentage spliced in index of splicing events were highly correlated. Our results indicate that the poly(A) capture full-length double-strand cDNA sequencing improves transcript quantification accuracy and the detection ability of alternative splicing events. It is also expected to contribute to the determination of the significance of DNA variants to splicing events.

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.

Transcriptome analysis of clock disrupted cancer cells reveals differential alternative splicing of cancer hallmarks genes

Emerging evidence points towards a regulatory role of the circadian clock in alternative splicing (AS). Whether alterations in core-clock components may contribute to differential AS events is largely unknown. To address this, we carried out a computational analysis on recently generated time-series RNA-seq datasets from three core-clock knockout (KO) genes (ARNTL, NR1D1, PER2) and WT of a colorectal cancer (CRC) cell line, and time-series RNA-seq datasets for additional CRC and Hodgkin’s lymphoma (HL) cells, murine WT, Arntl KO, and Nr1d1/2 KO, and murine SCN WT tissue. The deletion of individual core-clock genes resulted in the loss of circadian expression in crucial spliceosome components such as SF3A1 (in ARNTL^KO), SNW1 (in NR1D1^KO), and HNRNPC (in PER2^KO), which led to a differential pattern of KO-specific AS events. All HCT116^KO cells showed a rhythmicity loss of a crucial spliceosome gene U2AF1, which was also not rhythmic in higher progression stage CRC and HL cancer cells. AS analysis revealed an increase in alternative first exon events specific to PER2 and NR1D1 KO in HCT116 cells, and a KO-specific change in expression and rhythmicity pattern of AS transcripts related to cancer hallmarks genes including FGFR2 in HCT116_ARNTL^KO, CD44 in HCT116_NR1D1^KO, and MET in HCT116_PER2^KO. KO-specific changes in rhythmic properties of known spliced variants of these genes (e.g. FGFR2 IIIb/FGFR2 IIIc) correlated with epithelial-mesenchymal-transition signalling. Altogether, our bioinformatic analysis highlights a role for the circadian clock in the regulation of AS, and reveals a potential impact of clock disruption in aberrant splicing in cancer hallmark genes.

A comparison of mRNA sequencing (RNA-Seq) library preparation methods for transcriptome analysis

BACKGROUND

mRNA sequencing is a powerful technique, which is used to investigate the transcriptome status of a gene of interest, such as its transcription level and splicing variants. Presently, several RNA sequencing (RNA-Seq) methods have been developed; however, the relative advantage of each method has remained unknown. Here we used three commercially available RNA-Seq library preparation kits; the traditional method (TruSeq), in addition to full-length double-stranded cDNA methods (SMARTer and TeloPrime) to investigate the advantages and disadvantages of these three approaches in transcriptome analysis.

RESULTS

We observed that the number of expressed genes detected from the TeloPrime sequencing method was fewer than that obtained using the TruSeq and SMARTer. We also observed that the expression patterns between TruSeq and SMARTer correlated strongly. Alternatively, SMARTer and TeloPrime methods underestimated the expression of relatively long transcripts. Moreover, genes having low expression levels were undetected stochastically regardless of any three methods used. Furthermore, although TeloPrime detected a significantly higher proportion at the transcription start site (TSS), its coverage of the gene body was not uniform. SMARTer is proposed to be yielded for nonspecific genomic DNA amplification. In contrast, the detected splicing event number was highest in the TruSeq. The percent spliced in index (PSI) of the three methods was highly correlated.

CONCLUSIONS

TruSeq detected transcripts and splicing events better than the other methods and measured expression levels of genes, in addition to splicing events accurately. However, although detected transcripts and splicing events in TeloPrime were fewer, the coverage at TSS was highest. Additionally, SMARTer was better than TeloPrime with regards to the detected number of transcripts and splicing events among the understudied full-length double-stranded cDNA methods. In conclusion, for short-read sequencing, TruSeq has relative advantages for use in transcriptome analysis.

The aberrant upregulation of exon 10-inclusive SREK1 through SRSF10 acts as an oncogenic driver in human hepatocellular carcinoma

Deregulation of alternative splicing is implicated as a relevant source of molecular heterogeneity in cancer. However, the targets and intrinsic mechanisms of splicing in hepatocarcinogenesis are largely unknown. Here, we report a functional impact of a Splicing Regulatory Glutamine/Lysine-Rich Protein 1 (SREK1) variant and its regulator, Serine/arginine-rich splicing factor 10 (SRSF10). HCC patients with poor prognosis express higher levels of exon 10-inclusive SREK1 (SREK1^L). SREK1^L can sustain BLOC1S5-TXNDC5 (B-T) expression, a targeted gene of nonsense-mediated mRNA decay through inhibiting exon-exon junction complex binding with B-T to exert its oncogenic role. B-T plays its competing endogenous RNA role by inhibiting miR-30c-5p and miR-30e-5p, and further promoting the expression of downstream oncogenic targets SRSF10 and TXNDC5. Interestingly, SRSF10 can act as a splicing regulator for SREK1^L to promote hepatocarcinogenesis via the formation of a SRSF10-associated complex. In summary, we demonstrate a SRSF10/SREK1^L/B-T signalling loop to accelerate the hepatocarcinogenesis.

Alternative splicing is dysregulated in hepatocellular carcinoma. Here, the authors investigate the role of the splice variant of Splicing Regulatory Glutamic Acid and Lysine Rich Protein 1 (SREK1) and its upstream regulator, Serine/arginine-rich splicing factor 10 (SRSF10) in sustaining the oncogenic signal.

The Role of Extracellular Matrix Proteins in Breast Cancer

The extracellular matrix is a structure composed of many molecules, including fibrillar (types I, II, III, V, XI, XXIV, XXVII) and non-fibrillar collagens (mainly basement membrane collagens: types IV, VIII, X), non-collagenous glycoproteins (elastin, laminin, fibronectin, thrombospondin, tenascin, osteopontin, osteonectin, entactin, periostin) embedded in a gel of negatively charged water-retaining glycosaminoglycans (GAGs) such as non-sulfated hyaluronic acid (HA) and sulfated GAGs which are linked to a core protein to form proteoglycans (PGs). This highly dynamic molecular network provides critical biochemical and biomechanical cues that mediate the cell–cell and cell–matrix interactions, influence cell growth, migration and differentiation and serve as a reservoir of cytokines and growth factors’ action. The breakdown of normal ECM and its replacement with tumor ECM modulate the tumor microenvironment (TME) composition and is an essential part of tumorigenesis and metastasis, acting as key driver for malignant progression. Abnormal ECM also deregulate behavior of stromal cells as well as facilitating tumor-associated angiogenesis and inflammation. Thus, the tumor matrix modulates each of the classically defined hallmarks of cancer promoting the growth, survival and invasion of the cancer. Moreover, various ECM-derived components modulate the immune response affecting T cells, tumor-associated macrophages (TAM), dendritic cells and cancer-associated fibroblasts (CAF). This review article considers the role that extracellular matrix play in breast cancer. Determining the detailed connections between the ECM and cellular processes has helped to identify novel disease markers and therapeutic targets.

USP15 and USP4 facilitate lung cancer cell proliferation by regulating the alternative splicing of SRSF1

The deubiquitinating enzyme USP15 is implicated in several human cancers by regulating different cellular processes, including splicing regulation. However, the underlying molecular mechanisms of its functional relevance and the successive roles in enhanced tumorigenesis remain ambiguous. Here, we found that USP15 and its close paralog USP4 are overexpressed and facilitate lung cancer cell proliferation by regulating the alternative splicing of SRSF1. Depletion of USP15 and USP4 impair SRSF1 splicing characterized by the replacement of exon 4 with non-coding intron sequences retained at its C-terminus, resulting in an alternative isoform SRSF1-3. We observed an increased endogenous expression of SRSF1 in lung cancer cells as well, and its overexpression significantly enhanced cancer cell phenotype and rescued the depletion effect of USP15 and USP4. However, the alternatively spliced isoform SRSF1-3 was deficient in such aspects for its premature degradation through nonsense-mediated mRNA decay. The increased USP15 expression contributes to the lung adenocarcinoma (LUAD) development and shows significantly lower disease-specific survival of patients with USP15 alteration. In short, we identified USP15 and USP4 as key regulators of SRSF1 alternative splicing with altered functions, which may represent the novel prognostic biomarker as well as a potential target for LUAD.

Short-read and long-read RNA sequencing of mouse hematopoietic stem cells at bulk and single-cell levels

Hematopoietic stem cells (HSCs) lie at the top of the differentiation hierarchy. Although HSC and their immediate downstream, multipotent progenitors (MPP) have full multilineage differentiation capacity, only long-term (LT-) HSC has the capacity of long-term self-renewal. The heterogeneity within the HSC population is gradually acknowledged with the development of single-cell RNA sequencing and lineage tracing technologies. Transcriptional and post-transcriptional regulations play important roles in controlling the differentiation and self-renewal capacity within HSC population. Here we report a dataset comprising short- and long-read RNA sequencing for mouse long- and short-term HSC and MPP at bulk and single-cell levels. We demonstrate that integrating short- and long-read sequencing can facilitate the identification and quantification of known and unannotated isoforms. Thus, this dataset provides a groundwork for comprehensive and comparative studies on transcriptional diversity and heterogeneity within different HSC cell types.

Down Syndrome Cell Adhesion Molecule Triggers Membrane-to-Nucleus Signaling-Regulated Hemocyte Proliferation against Bacterial Infection in Invertebrates

Down syndrome cell adhesion molecule (Dscam) generates tens of thousands of isoforms by alternative splicing, thereby providing crucial functions during immune responses. In this study, a novel Dscam signaling pathway was investigated in crab, which remains poorly characterized in invertebrates. Bacterial infection induced the cytoplasmic cleavage of Dscam intracellular domains (ICDs) by γ-secretase, and then the released ICDs carrying specific alternatively spliced exons could directly interact with IPO5 to facilitate nuclear translocation. Nuclear imported ICDs thus promoted hemocyte proliferation and protect the host from bacterial infection. Protein-interaction studies revealed that the ectodomain of Dscam bound to a disintegrin and metalloprotease domain 10 (ADAM10) rather than ADAM17. Inhibition or overexpression of ADAM10 impaired or accelerated Dscam shedding activity post-bacterial stimulation, respectively. Moreover, the shedding signal then mediated Dscam with an intact cytoplasmic domain to promote the cleavage of ICDs by γ-secretase. Furthermore, the transcription of ADAM10 was regulated by Dscam-induced canonical signaling, but not nuclear imported ICDs, to serve as a feedback regulation between two different Dscam pathways. Thus, membrane-to-nuclear signaling of Dscam regulated hemocyte proliferation in response to bacterial infection.

Comprehensive characterization of somatic variants associated with intronic polyadenylation in human cancers

Somatic single nucleotide variants (SNVs) in cancer genome affect gene expression through various mechanisms depending on their genomic location. While somatic SNVs near canonical splice sites have been reported to cause abnormal splicing of cancer-related genes, whether these SNVs can affect gene expression through other mechanisms remains an open question. Here, we analyzed RNA sequencing and exome data from 4,998 cancer patients covering ten cancer types and identified 152 somatic SNVs near splice sites that were associated with abnormal intronic polyadenylation (IPA). IPA-associated somatic variants favored the localization near the donor splice sites compared to the acceptor splice sites. A proportion of SNV-associated IPA events overlapped with premature cleavage and polyadenylation events triggered by U1 small nuclear ribonucleoproteins (snRNP) inhibition. GC content, intron length and polyadenylation signal were three genomic features that differentiated between SNV-associated IPA and intron retention. Notably, IPA-associated SNVs were enriched in tumor suppressor genes (TSGs), including the well-known TSGs such as PTEN and CDH1 with recurrent SNV-associated IPA events. Minigene assay confirmed that SNVs from PTEN, CDH1, VEGFA, GRHL2, CUL3 and WWC2 could lead to IPA. This work reveals that IPA acts as a novel mechanism explaining the functional consequence of somatic SNVs in human cancer.

Cancer DEIso: An integrative analysis platform for investigating differentially expressed gene-level and isoform-level human cancer markers

Transcript isoforms regulated by alternative splicing can substantially impact carcinogenesis, leading to a need to obtain clues for both gene differential expression and malfunctions of isoform distributions in cancer studies. The Cancer Genome Atlas (TCGA) project was launched in 2008 to collect cancer-related genome mutation raw data from the population. While many repositories tried to add insights into the raw data in TCGA, no existing database provides both comprehensive gene-level and isoform-level cancer stage marker investigation and survival analysis. We constructed Cancer DEIso to facilitate in-depth analyses for both gene-level and isoform-level human cancer studies. Patient RNA-seq data, sample sheets, patient clinical data, and human genome datasets were collected and processed in Cancer DEIso. And four functions to search differentially expressed genes/isoforms between cancer stages were implemented: (i) Search potential gene/isoform markers for a specified cancer type and its two stages; (ii) Search potentially induced cancer types and stages for a gene/isoform; (iii) Expression survival analysis on a given gene/isoform for some cancer; (iv) Gene/isoform stage expression comparison visualization. As an example, we demonstrate that Cancer DEIso can indicate potential colorectal cancer isoform diagnostic markers that are not easily detected when only gene-level expressions are considered. Cancer DEIso is available at http://cosbi4.ee.ncku.edu.tw/DEIso/.

A pan-cancer analysis of alternative splicing of splicing factors in 6904 patients

Great progress has been made in the investigation on mutation and expression of splicing factor. However, little is known on the role of alternative splicing of splicing factors across cancers. Here, we reported a pan-cancer analysis of alternative splicing of splicing factors spanning 6904 patients across 16 cancer types, and identified 167 splicing factors with implications regulating cancer-specific splicing patterns through alternative splicing. Furthermore, we found that abnormal splicing events of splicing factors could serve as potential common regulators for alternative splicing in different cancers. In addition, we developed a splicing-derived neoepitopes database (ASPNs), which provided the corresponding putative alternative splicing-derived neoepitopes of 16 cancer types. Our results suggested that alternative splicing of splicing factors involved in the pre-RNA splicing process was common across cancer types and may represent an underestimated hallmark of tumorigenesis.

Anastasis: Return Journey from Cell Death

For over 20 years, it has been a dogma that once the integrity of mitochondria is disrupted and proapoptotic proteins that are normally located in the intermembrane space of mitochondria appeared in the cytoplasm, the process of cell death becomes inevitable. However, it has been recently shown that upon removal of the death signal, even at the stage of disturbance in the mitochondria, cells can recover and continue to grow. This phenomenon was named anastasis. Here, we will critically discuss the present knowledge concerning the mechanisms of cell death reversal, or development of anastasis, methods for its detection, and what role signaling from different intracellular compartments plays in anastasis stimulation.

SUVA: splicing site usage variation analysis from RNA-seq data reveals highly conserved complex splicing biomarkers in liver cancer

Most of the current alternative splicing (AS) analysis tools are powerless to analyse complex splicing. To address this, we developed SUVA (Splice sites Usage Variation Analysis) that decomposes complex splicing events into five types of splice junction pairs. By analysing real and simulated data, SUVA showed higher sensitivity and accuracy in detecting AS events than the compared methods. Notably, SUVA detected extensive complex AS events and screened out 69 highly conserved and dominant AS events associated with cancer. The cancer-associated complex AS events in FN1 and the co-regulated RNA-binding proteins were significantly correlated with patient survival.

The Exon Junction Complex Core Represses Cancer-Specific Mature mRNA Re-splicing: A Potential Key Role in Terminating Splicing

Using TSG101 pre-mRNA, we previously discovered cancer-specific re-splicing of mature mRNA that generates aberrant transcripts/proteins. The fact that mRNA is aberrantly re-spliced in various cancer cells implies there must be an important mechanism to prevent deleterious re-splicing on the spliced mRNA in normal cells. We thus postulated that mRNA re-splicing is controlled by specific repressors, and we searched for repressor candidates by siRNA-based screening for mRNA re-splicing activity. We found that knock-down of EIF4A3, which is a core component of the exon junction complex (EJC), significantly promoted mRNA re-splicing. Remarkably, we could recapitulate cancer-specific mRNA re-splicing in normal cells by knock-down of any of the core EJC proteins, EIF4A3, MAGOH, or RBM8A (Y14), implicating the EJC core as the repressor of mRNA re-splicing often observed in cancer cells. We propose that the EJC core is a critical mRNA quality control factor to prevent over-splicing of mature mRNA.

Identification of Novel Alternative Splicing Events Associated With Tumorigenesis, Protein Modification, and Immune Microenvironment in Early-Onset Gastric Cancer

Background

Alternative splicing (AS), e.g. the tandem alternative polyadenylation (TAPA), has emerged as major post-transcriptional modification events in human disease. However, the roles of the AS and TAPA in early-onset gastric cancer (EOGC) have not been revealed.

Methods

The global AS profiles of 80 EOGC patients were analyzed. The EOGC-specific AS events (ESASs) were identified in both the EOGC and adjacent non-tumor tissues. The functional enrichment analysis, Splicing network, Alternative Polyadenylation (APA) core factor network, and cell abundancy analysis were performed. Furthermore, the landscapes of the AS events in the varied subtypes of the EOGC patients were evaluated.

Results

Overall, 66,075 AS events and 267 ESASs were identified in the EOGC. Furthermore, 4809 genes and 6152 gene isoforms were found to be aberrantly expressed in the EOGC. The Gene Ontology (GO) and Kyoto Encyclopedia of Gene and Genome (KEGG) pathway analyses showed that the significant pathway alterations might exist in these AS events, genes, and gene isoforms. Moreover, the Protein-protein interaction (PPI) network analysis revealed that the UBC, NEK2, EPHB2, and DCTN1 genes were the hub genes in the AS events in the EOGC. The immune cell infiltration analysis indicated a correlation between the AS events and the cancer immune microenvironment. The distribution of the AS events in varied EOGC subtypes, protein phosphorylation and glycosylation was uneven.

Conclusion

The study highlighted the vital roles of the AS in the EOGC, including modulating the specific protein modification and reshaping the cancer immune microenvironment, and yielded new insights into the diagnosis of the EOGC as well as cancer treatment.

The landscape and driver potential of site-specific hotspots across cancer genomes

Large sets of whole cancer genomes make it possible to study mutation hotspots genome-wide. Here we detect, categorize, and characterize site-specific hotspots using 2279 whole cancer genomes from the Pan-Cancer Analysis of Whole Genomes project and provide a resource of annotated hotspots genome-wide. We investigate the excess of hotspots in both protein-coding and gene regulatory regions and develop measures of positive selection and functional impact for individual hotspots. Using cancer allele fractions, expression aberrations, mutational signatures, and a variety of genomic features, such as potential gain or loss of transcription factor binding sites, we annotate and prioritize all highly mutated hotspots. Genome-wide we find more high-frequency SNV and indel hotspots than expected given mutational background models. Protein-coding regions are generally enriched for SNV hotspots compared to other regions. Gene regulatory hotspots show enrichment of potential same-patient second-hit missense mutations, consistent with enrichment of hotspot driver mutations compared to singletons. For protein-coding regions, splice-sites, promoters, and enhancers, we see an excess of hotspots associated with cancer genes. Interestingly, missense hotspot mutations in tumor suppressors are associated with elevated expression, suggesting localized amino-acid changes with functional impact. For individual non-coding hotspots, only a small number show clear signs of positive selection, including known sites in the TERT promoter and the 5' UTR of TP53. Most of the new candidates have few mutations and limited driver evidence. However, a hotspot in an enhancer of the oncogene POU2AF1, which may create a transcription factor binding site, presents multiple lines of driver-consistent evidence.

Tumor reversion: a dream or a reality

Reversion of tumor to a normal differentiated cell once considered a dream is now at the brink of becoming a reality. Different layers of molecules/events such as microRNAs, transcription factors, alternative RNA splicing, post-transcriptional, post-translational modifications, availability of proteomics, genomics editing tools, and chemical biology approaches gave hope to manipulation of cancer cells reversion to a normal cell phenotype as evidences are subtle but definitive. Regardless of the advancement, there is a long way to go, as customized techniques are required to be fine-tuned with precision to attain more insights into tumor reversion. Tumor regression models using available genome-editing methods, followed by in vitro and in vivo proteomics profiling techniques show early evidence. This review summarizes tumor reversion developments, present issues, and unaddressed challenges that remained in the uncharted territory to modulate cellular machinery for tumor reversion towards therapeutic purposes successfully. Ongoing research reaffirms the potential promises of understanding the mechanism of tumor reversion and required refinement that is warranted in vitro and in vivo models of tumor reversion, and the potential translation of these into cancer therapy. Furthermore, therapeutic compounds were reported to induce phenotypic changes in cancer cells into normal cells, which will contribute in understanding the mechanism of tumor reversion. Altogether, the efforts collectively suggest that tumor reversion will likely reveal a new wave of therapeutic discoveries that will significantly impact clinical practice in cancer therapy.

Development and structural characterisation of human scFv targeting MDM2 spliced variant MDM215kDa

The murine double minute 2 (MDM2) protein is a major negative regulator of the tumour suppressor protein p53. Under normal conditions, MDM2 constantly binds to p53 transactivation domain and/or ubiquinates p53 via its role as E3 ubiquitin ligase to promote p53 degradation as well as nuclear export to maintain p53 levels in cells. Meanwhile, amplification of MDM2 and appearance of MDM2 spliced variants occur in many tumours and normal tissues making it a prognostic indicator for human cancers. The mutation or deletion of p53 protein in half of human cancers inactivates its tumour suppressor activity. However, cancers with wild type p53 have its function effectively inhibited through direct interaction with MDM2 oncoprotein. Here, we described the construction of a MDM2 spliced variant (rMDM2^15kDa) consisting of SWIB/MDM2 domain and its central region for antibody generation. Biopanning with a human naïve scFv library generated four scFv clones specific to rMDM2^15kDa. Additionally, the selected scFv clones were able to bind to the recombinant full length MDM2 (rMDM2-FL). Computational prediction showed that the selected scFv clones potentially bind to exon 7-8 of MDM2 while leaving the MDM2/SWIB domain free for p53 interaction. The developed antibodies exhibit good specificity can be further investigated for downstream biomedical and research applications.

lncRNA TRMP-S directs dual mechanisms to regulate p27-mediated cellular senescence

Long noncoding RNAs (lncRNAs) undergo extensive alternative splicing, but little is known about isoform functions. A prior investigation of lncRNA RP11-369C8.1 reported that its splice variant TRMP suppressed p27 translation through PTBP1. Here we characterize a second major splice variant, TRMP-S (short variant), whose enforced loss promotes cancer cell-cycle arrest and p27-dependent entry into cellular senescence. Remarkably, despite sharing a single common exon with TRMP, TRMP-S restrains p27 expression through distinct mechanisms. First, TRMP-S stabilizes UHRF1 protein levels, an epigenetic inhibitor of p27, by promoting interactions between UHRF1 and its deubiquitinating enzyme USP7. Alternatively, binding interactions between TRMP-S and FUBP3 prevent p53 mRNA interactions with RPL26 ribosomal protein, the latter essential for promoting p53 translation with ensuing suppression of p53 translation limiting p27 expression. Significantly, as TRMP-S is itself transactivated by p53, this identifies negative feedback regulation between p53 and TRMP-S. Different splicing variants of the RP11-369C8.1 gene thereby exert distinct roles that converge on the homeostatic control of p27 expression, providing an important precedent for understanding the actions of alternatively spliced lncRNAs.

Development of alternative splicing signature in lung squamous cell carcinoma

Increasing evidence demonstrated that alternative splicing (AS) plays a vital role in tumorigenesis and clinical outcome of patient. However, systematical analysis of AS in lung squamous cell carcinoma (LUSC) is lacking and greatly necessary. Thus, this study was to systematically estimate the function of AS events served as prognostic indicators in LUSC. Among 31,345 mRNA AS events in 9633 genes, we detected 1996 AS in 1409 genes which have significant connection with overall survival (OS) of LUSC patients. Then, prognostic model based on seven types of AS events was established and we further constructed a combined prognostic model. The Kaplan–Meier curve results suggested that seven types of AS signatures and the combined prognostic model could exhibit robust performance in predicting prognosis. Patients in the high-risk group had significantly shorter OS than those in the low-risk group. The ROC showed all prognostic models had high accuracy and powerful predictive performance with different AUC ranging from 0.837 to 0.978. Moreover, the combined prognostic model had highest performance in risk stratification and predictive accuracy than single prognostic models and had higher accuracy than other mRNA model. Finally, a significant correlation network between survival-related AS genes and prognostic splicing factors (SFs) was established. In conclusion, our study provided several potential prognostic AS models and constructed splicing network between AS and SFs in LUSC, which could be used as potential indicators and treatment targets for LUSC patients.