MERIT: Systematic Analysis and Characterization of Mutational Effect on RNA Interactome Topology

Single-cell RNA binding protein regulatory network analyses reveal oncogenic HNRNPK-MYC signalling pathway in cancer

RNA-binding proteins (RBPs) are key players of gene expression and perturbations of RBP-RNA regulatory network have been observed in various cancer types. Here, we propose a computational method, RBPreg, to identify the RBP regulators by integration of single cell RNA-Seq (N = 233,591) and RBP binding data. Pan-cancer analyses suggest that RBP regulators exhibit cancer and cell specificity and perturbations of RBP regulatory network are involved in cancer hallmark-related functions. We prioritize an oncogenic RBP-HNRNPK, which is highly expressed in tumors and associated with poor prognosis of patients. Functional assays performed in cancer cells reveal that HNRNPK promotes cancer cell proliferation, migration, and invasion in vitro and in vivo. Mechanistic investigations further demonstrate that HNRNPK promotes tumorigenesis and progression by directly binding to MYC and perturbed the MYC targets pathway in lung cancer. Our results provide a valuable resource for characterizing RBP regulatory networks in cancer, yielding potential biomarkers for precision medicine.

Long noncoding RNA LINC01132 enhances immunosuppression and therapy resistance via NRF1/DPP4 axis in hepatocellular carcinoma

Background

Long noncoding RNAs (lncRNAs) are emerging as critical regulators of gene expression and play fundamental roles in various types of cancer. Current developments in transcriptome analyses unveiled the existence of lncRNAs; however, their functional characterization remains a challenge.

Methods

A bioinformatics screen was performed by integration of multiple omics data in hepatocellular carcinoma (HCC) prioritizing a novel oncogenic lncRNA, LINC01132. Expression of LINC01132 in HCC and control tissues was validated by qRT-PCR. Cell viability and migration activity was examined by MTT and transwell assays. Finally, our results were confirmed in vivo mouse model and ex vivo patient derived tumor xenograft experiments to determine the mechanism of action and explore LINC01132-targeted immunotherapy.

Results

Systematic investigation of lncRNAs genome-wide expression patterns revealed LINC01132 as an oncogene in HCC. LINC01132 is significantly overexpressed in tumor and associated with poor overall survival of HCC patients, which is mainly driven by copy number amplification. Functionally, LINC01132 overexpression promoted cell growth, proliferation, invasion and metastasis in vitro and in vivo. Mechanistically, LINC01132 acts as an oncogenic driver by physically interacting with NRF and enhancing the expression of DPP4. Notably, LINC01132 silencing triggers CD8+ T cells infiltration, and LINC01132 knockdown combined with anti-PDL1 treatment improves antitumor immunity, which may prove a new combination therapy in HCC.

Conclusions

LINC01132 functions as an oncogenic driver that induces HCC development via the NRF1/DPP4 axis. Silencing LINC01132 may enhance the efficacy of anti-PDL1 immunotherapy in HCC patients.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13046-022-02478-z.

Robust prognostic model based on immune infiltration-related genes and clinical information in ovarian cancer

Immune infiltration of ovarian cancer (OV) is a critical factor in determining patient's prognosis. Using data from TCGA and GTEx database combined with WGCNA and ESTIMATE methods, 46 genes related to OV occurrence and immune infiltration were identified. Lasso and multivariate Cox regression were applied to define a prognostic score (IGCI score) based on 3 immune genes and 3 types of clinical information. The IGCI score has been verified by K‐M curves, ROC curves and C‐index on test set. In test set, IGCI score (C‐index = 0.630) is significantly better than AJCC stage (C‐index = 0.541, p < 0.05) and CIN25 (C‐index = 0.571, p < 0.05). In addition, we identified key mutations to analyse prognosis of patients and the process related to immunity. Chi‐squared tests revealed that 6 mutations are significantly (p < 0.05) related to immune infiltration: BRCA1, ZNF462, VWF, RBAK, RB1 and ADGRV1. According to mutation survival analysis, we found 5 key mutations significantly related to patient prognosis (p < 0.05): CSMD3, FLG2, HMCN1, TOP2A and TRRAP. RB1 and CSMD3 mutations had small p‐value (p < 0.1) in both chi‐squared tests and survival analysis. The drug sensitivity analysis of key mutation showed when RB1 mutation occurs, the efficacy of six anti‐tumour drugs has changed significantly (p < 0.05).

Global RNA editing identification and characterization during human pluripotent-to-cardiomyocyte differentiation

RNA editing is widely involved in stem cell differentiation and development; however, RNA editing events during human cardiomyocyte differentiation have not yet been characterized and elucidated. Here, we identified genome-wide RNA editing sites and systemically characterized their genomic distribution during four stages of human cardiomyocyte differentiation. It was found that the expression level of ADAR1 affected the global number of adenosine to inosine (A-to-I) editing sites but not the editing degree. Next, we identified 43, 163, 544, and 141 RNA editing sites that contribute to changes in amino acid sequences, variation in alternative splicing, alterations in miRNA-target binding, and changes in gene expression, respectively. Generally, RNA editing showed a stage-specific pattern with 211 stage-shared editing sites. Interestingly, cardiac muscle contraction and heart-disease-related pathways were enriched by cardio-specific editing genes, emphasizing the connection between cardiomyocyte differentiation and heart diseases from the perspective of RNA editing. Finally, it was found that these RNA editing sites are also related to several congenital and noncongenital heart diseases. Together, our study provides a new perspective on cardiomyocyte differentiation and offers more opportunities to understand the mechanisms underlying cell fate determination, which can promote the development of cardiac regenerative medicine and therapies for human heart diseases.

ImmReg: the regulon atlas of immune-related pathways across cancer types

Immune system gene regulation perturbation has been found to be a major cause of the development of various types of cancer. Numbers of mechanisms contribute to gene expression regulation, thus, systematically identification of potential regulons of immune-related pathways is critical to cancer immunotherapy. Here, we comprehensively chart the landscape of transcription factors, microRNAs, RNA binding proteins and long noncoding RNAs regulation in 17 immune-related pathways across 33 cancers. The potential immunology regulons are likely to exhibit higher expressions in immune cells, show expression perturbations in cancer, and are significantly correlated with immune cell infiltrations. We also identify a panel of clinically relevant immunology regulons across cancers. Moreover, the regulon atlas of immune-related pathways helps prioritizing cancer-related genes (i.e. ETV7, miR-146a-5p, ZFP36 and HCP5). We further identified two molecular subtypes of glioma (cold and hot tumour phenotypes), which were characterized by differences in immune cell infiltrations, expression of checkpoints, and prognosis. Finally, we developed a user-friendly resource, ImmReg (http://bio-bigdata.hrbmu.edu.cn/ImmReg/), with multiple modules to visualize, browse, and download immunology regulation. Our study provides a comprehensive landscape of immunology regulons, which will shed light on future development of RNA-based cancer immunotherapies.

Development and Validation of a Prognostic Model of RNA-Binding Proteins in Colon Adenocarcinoma: A Study Based on TCGA and GEO Databases

Background

Previous studies reported that dysregulation of RNA-binding proteins (RBPs) is significantly associated with the development of cancer. However, there are few studies to date on the role of RBPs in colon adenocarcinoma (COAD).

Methods

RNA sequencing and clinical data for COAD patients were downloaded from The Cancer Genome Atlas (TCGA) database to identify differentially expressed (DE) RBPs between COAD tissue and normal colon tissue, and then the expression and prognostic significance of these RBPs were investigated in detail by systematic bioinformatics analysis. qRT-PCR was used to validate the expressions of prognosis-related RBP-encoding genes.

Results

Seven RBPs (RPL10L, ERI1, POP1, CAPRIN2, TDRD7, SNIP1 and PPARGC1A) were identified as hub genes associated with prognosis by a series of regression analyses, and were then used to construct a prognostic model. Further analysis based on this model indicated that the overall survival (OS) of the high-risk groups was lower than that of the low-risk groups. In this prognostic model, the area under the ROC curve (AUC) was 0.694, 0.709 and 0.665 for the TCGA cohort at 1, 3 and 5 years, respectively, while the AUC was 0.671, 0.633 and 0.601 for the GEO combined cohort at 1, 3 and 5 years, respectively, indicating the good predictive ability of the model. We also built a nomogram based on the 7 RBPs in the TCGA cohort, and the model showed good discriminatory ability for COAD.

Conclusion

We screened seven prognosis-related genes in COAD patients based on RBP-related genes, validated the expressions of the seven prognosis-related RBP-encoding genes by qRT-PCR and constructed a prognosis-related nomogram for patients with COAD.

e-MutPath: computational modeling reveals the functional landscape of genetic mutations rewiring interactome networks

Understanding the functional impact of cancer somatic mutations represents a critical knowledge gap for implementing precision oncology. It has been increasingly appreciated that the interaction profile mediated by a genomic mutation provides a fundamental link between genotype and phenotype. However, specific effects on biological signaling networks for the majority of mutations are largely unknown by experimental approaches. To resolve this challenge, we developed e-MutPath (edgetic Mutation-mediated Pathway perturbations), a network-based computational method to identify candidate ‘edgetic’ mutations that perturb functional pathways. e-MutPath identifies informative paths that could be used to distinguish disease risk factors from neutral elements and to stratify disease subtypes with clinical relevance. The predicted targets are enriched in cancer vulnerability genes, known drug targets but depleted for proteins associated with side effects, demonstrating the power of network-based strategies to investigate the functional impact and perturbation profiles of genomic mutations. Together, e-MutPath represents a robust computational tool to systematically assign functions to genetic mutations, especially in the context of their specific pathway perturbation effect.

Molecular biomarkers in early stage lung cancer

Low dose computed tomography (LDCT) screening, together with the recent advances in targeted and immunotherapies, have shown to improve non-small cell lung cancer (NSCLC) survival. Furthermore, screening has increased the number of early stage-detected tumors, allowing for surgical resection and multimodality treatments when needed. The need for improved sensitivity and specificity of NSCLC screening has led to increased interest in combining clinical and radiological data with molecular data. The development of biomarkers is poised to refine inclusion criteria for LDCT screening programs. Biomarkers may also be useful to better characterize the risk of indeterminate nodules found in the course of screening or to refine prognosis and help in the management of screening detected tumors. The clinical implications of these biomarkers are still being investigated and whether or not biomarkers will be included in further decision-making algorithms in the context of screening and early lung cancer management still needs to be determined. However, it seems clear that there is much room for improvement even in early stage lung cancer disease-free survival (DFS) rates; thus, biomarkers may be the key to refine risk-stratification and treatment of these patients. Clinicians’ capacity to register, integrate, and analyze all the available data in both high risk individuals and early stage NSCLC patients will lead to a better understanding of the disease’s mechanisms, and will have a direct impact in diagnosis, treatment, and follow up of these patients. In this review, we aim to summarize all the available data regarding the role of biomarkers in LDCT screening and early stage NSCLC from a multidisciplinary perspective. We have highlighted clinical implications, the need to combine risk stratification, clinical data, radiomics, molecular information and artificial intelligence in order to improve clinical decision-making, especially regarding early diagnostics and adjuvant therapy. We also discuss current and future perspectives for biomarker implementation in routine clinical practice.

Dynamic Expression of m6A Regulators During Multiple Human Tissue Development and Cancers

N6-methyladenosine (m⁶A) plays critical roles in human development and cancer progression. However, our knowledge regarding the dynamic expression of m⁶A regulators during human tissue development is still lacking. Here, we comprehensively analyzed the dynamic expression alterations of m⁶A regulators during seven tissue development and eight cancer types. We found that m⁶A regulators globally exhibited decreased expression during development. In addition, IGF2BP1/2/3 (insulinlike growth factor 2 MRNA-binding protein 1/2/3) exhibited reverse expression pattern in cancer progression, suggesting an oncofetal reprogramming in cancer. The expressions of IGF2BP1/2/3 were regulated by genome alterations, particularly copy number amplification in cancer. Clinical association analysis revealed that higher expressions of IGF2BP1/2/3 were associated with worse survival of cancer patients. Finally, we found that genes significantly correlated with IGF2BP1/2/3 were significantly enriched in cancer hallmark-related pathways. In summary, dynamic expression analysis will guide both mechanistic and therapeutic roles of m⁶A regulators during tissue development and cancer progression.

Development of prognostic signature based on RNA binding proteins related genes analysis in clear cell renal cell carcinoma

RNA binding proteins (RBPs) play significant roles in the development of tumors. However, a comprehensive analysis of the biological functions of RBPs in clear cell renal cell carcinoma (ccRCC) has not been performed. Our study aimed to construct an RBP-related risk model for prognosis prediction in ccRCC patients. First, RNA sequencing data of ccRCC were downloaded from The Cancer Genome Atlas (TCGA) database. Three RBP genes (EIF4A1, CARS, and RPL22L1) were validated as prognosis-related hub genes by univariate and multivariate Cox regression analyses and were integrated into a prognostic model by least absolute shrinkage and selection operator (LASSO) Cox regression analysis. According to this model, patients with high risk scores displayed significantly worse overall survival (OS) than those with low risk scores. Moreover, the multivariate Cox analysis results indicated that risk score, tumor grade, and tumor stage were significantly correlated with patient OS. A nomogram was constructed based on the three RBP genes and showed a good ability to predict outcomes in ccRCC patients. In conclusion, this study identified a three-RBP gene risk model for predicting the prognosis of patients, which is conducive to the identification of novel diagnostic and prognostic molecular markers.

Systematic analysis of enhancer regulatory circuit perturbation driven by copy number variations in malignant glioma

Background: Enhancers are emerging regulatory regions controlling gene expression in diverse cancer types. However, the functions of enhancer regulatory circuit perturbations driven by copy number variations (CNVs) in malignant glioma are unclear. Therefore, we aimed to investigate the comprehensive enhancer regulatory perturbation and identify potential biomarkers in glioma.

Results: We performed a meta-analysis of the enhancer centered regulatory circuit perturbations in 683 gliomas by integrating CNVs, gene expression, and transcription factors (TFs) binding. We found widespread CNVs of enhancers during glioma progression, and CNVs were associated with the perturbations of enhancer activities. In particular, the degree of perturbations for amplified enhancers was much greater accompanied by the glioma malignant progression. In addition, CNVs and enhancers cooperatively regulated the expressions of cancer-related genes. Genome-wide TF binding profiles revealed that enhancers were pervasively regulated by TFs. A network-based analysis of TF-enhancer-gene regulatory circuits revealed a core TF-gene module (58 interactions including seven genes and 14 TFs) that was associated survival of patients with glioma (p < 0.001). Finally, we validated this prognosis-associated TF-gene regulatory module in an independent cohort. In summary, our analyses provided new molecular insights for enhancer-centered transcriptional perturbation in glioma therapy.

Conclusion: Integrative analysis revealed enhancer regulatory perturbations in glioma and also identified a network module that was associated with patient survival, thereby providing novel insights into enhancer-centered cancer therapy.

Molecular Characterization and Clinical Relevance of RNA Binding Proteins in Colorectal Cancer

Abnormal expression of RNA binding proteins (RBPs) has been reported across various cancers. However, the potential role of RBPs in colorectal cancer (CRC) remains unclear. In this study, we performed a systematic bioinformatics analysis of RBPs in CRC. We downloaded CRC data from The Cancer Genome Atlas (TCGA) database. Our analysis identified 242 differentially expressed RBPs between tumor and normal tissues, including 200 upregulated and 42 downregulated RBPs. Next, we found eight RBPs (RRS1, PABPC1L, TERT, SMAD6, UPF3B, RP9, NOL3, and PTRH1) related to the prognoses of CRC patients. Among these eight prognosis-related RBPs, four RBPs (NOL3, PTRH1, UPF3B, and SMAD6) were selected to construct a prognostic risk score model. Furthermore, our results indicated that the prognostic risk score model accurately predicted the prognosis of CRC patients [area under the receiver operating characteristic curve (AUC)for 3- and 5-year overall survival (OS) and was 0.645 and 0.672, respectively]. Furthermore, we developed a nomogram based on a prognostic risk score model. The nomogram was able to demonstrate the wonderful performance in predicting 3- and 5-year OS. Additionally, we validated the clinical value of four risk genes in the prognostic risk score model and identified that these risk genes were associated with tumorigenesis, lymph node metastasis, distant metastasis, clinical stage, and prognosis. Finally, we used the TIMER and Human Protein Atlas (HPA)database to validate the expression of four risk genes at the transcriptional and translational levels, respectively, and used a clinical cohort to validate the roles of NOL3 and UPF3B in predicting the prognosis of CRC patients. In summary, our study demonstrated that RBPs have an effect on CRC tumor progression and might be potential prognostic biomarkers for CRC patients.

Complex impact of DNA methylation on transcriptional dysregulation across 22 human cancer types

Accumulating evidence has demonstrated that transcriptional regulation is affected by DNA methylation. Understanding the perturbation of DNA methylation-mediated regulation between transcriptional factors (TFs) and targets is crucial for human diseases. However, the global landscape of DNA methylation-mediated transcriptional dysregulation (DMTD) across cancers has not been portrayed. Here, we systematically identified DMTD by integrative analysis of transcriptome, methylome and regulatome across 22 human cancer types. Our results revealed that transcriptional regulation was affected by DNA methylation, involving hundreds of methylation-sensitive TFs (MethTFs). In addition, pan-cancer MethTFs, the regulatory activity of which is generally affected by DNA methylation across cancers, exhibit dominant functional characteristics and regulate several cancer hallmarks. Moreover, pan-cancer MethTFs were found to be affected by DNA methylation in a complex pattern. Finally, we investigated the cooperation among MethTFs and identified a network module that consisted of 43 MethTFs with prognostic potential. In summary, we systematically dissected the transcriptional dysregulation mediated by DNA methylation across cancer types, and our results provide a valuable resource for both epigenetic and transcriptional regulation communities.

Integrated Analysis of the Functions and Prognostic Values of RNA Binding Proteins in Lung Squamous Cell Carcinoma

Lung cancer is the leading cause of cancer-related deaths worldwide. Dysregulation of RNA binding proteins (RBPs) has been found in a variety of cancers and is related to oncogenesis and progression. However, the functions of RBPs in lung squamous cell carcinoma (LUSC) remain unclear. In this study, we obtained gene expression data and corresponding clinical information for LUSC from The Cancer Genome Atlas (TCGA) database, identified aberrantly expressed RBPs between tumors and normal tissue, and conducted a series of bioinformatics analyses to explore the expression and prognostic value of these RBPs. A total of 300 aberrantly expressed RBPs were obtained, comprising 59 downregulated and 241 upregulated RBPs. Functional enrichment analysis indicated that the differentially expressed RBPs were mainly associated with mRNA metabolic processes, RNA processing, RNA modification, regulation of translation, the TGF-beta signaling pathway, and the Toll-like receptor signaling pathway. Nine RBP genes (A1CF, EIF2B5, LSM1, LSM7, MBNL2, RSRC1, TRMU, TTF2, and ZCCHC5) were identified as prognosis-associated hub genes by univariate, least absolute shrinkage and selection operator (LASSO), Kaplan–Meier survival, and multivariate Cox regression analyses, and were used to construct the prognostic model. Further analysis demonstrated that high risk scores for patients were significantly related to poor overall survival according to the model. The area under the time-dependent receiver operator characteristic curve of the prognostic model was 0.712 at 3 years and 0.696 at 5 years. We also developed a nomogram based on nine RBP genes, with internal validation in the TCGA cohort, which showed a favorable predictive efficacy for prognosis in LUSC. Our results provide novel insights into the pathogenesis of LUSC. The nine-RBP gene signature showed predictive value for LUSC prognosis, with potential applications in clinical decision-making and individualized treatment.

RBP EIF2S2 Promotes Tumorigenesis and Progression by Regulating MYC-Mediated Inhibition via FHIT-Related Enhancers

RNA-binding proteins (RBPs) play fundamental roles in cancer; however, we still lack knowledge about to what extent RBPs are dysregulated, as well as about perturbed signaling pathways in cancer. In this study, we integrated analysis of multidimensional data across >10,000 cancer patients and >1,000 cell lines. We identified a top candidate RBP: eukaryotic translation initiation factor 2 subunit beta (EIF2S2). EIF2S2 is highly expressed in tumors and is associated with malignant features as well as patient prognosis. Functional assays performed in cancer cells revealed that EIF2S2 promotes cancer cell proliferation, migration, and invasion in vitro as well as tumor growth and metastasis in vivo. Mechanistic investigations further demonstrated that EIF2S2 promotes tumorigenesis and progression by directly binding to a long non-coding RNA, LINC01600, which physically interacts with the MYC protein and increases its stability. Interestingly, we revealed that the EIF2S2-LINC01600-MYC axis can activate the Wnt/β-catenin pathway by inhibiting the activity of FHIT-related enhancers and FHIT expression. Finally, EIF2S2 knockdown combined with oxaliplatin treatment could be a potential combination therapy in cancer. Our integrated analysis provided detailed knowledge of the function of the EIF2S2-LINC01600-MYC axis, which will facilitate the development of rational combination therapies for cancer.

Systematically Dissecting the Function of RNA-Binding Proteins During Glioma Progression

RNA-binding proteins (RBPs) play important roles in regulating gene expression and dysregulation of RBPs have been observed in various types of cancer. However, the role of RBPs during glioma progression, and particular in Chinese patients, is only starting to be unveiled. Here, we systematically analyzed the somatic mutation, gene expression patterns of 2949 RBPs during glioma progression. Our comprehensive study reveals several of highly mutated genes (such as ATRX, TTN and SETD2) and differentially expressed genes (such as KIF4A, TTK and CEP55). Integration of the expression of RBPs and genes, we constructed a regulatory network in glioma and revealed the functional links between RBPs and cancer-related genes. Moreover, we identified the prognosis spectrum of RBPs during glioma progression. The expression of a number of RBPs, such as SNRPN and IGF2BP3, are significantly associated with overall survival of patients in all grades. Taken together, our analyses provided a valuable RBP resource during glioma progression, and revealed several candidates that potentially contribute to development of therapeutic targets for glioma.

The role of the protein-RNA recognition code in neurodegeneration

MicroRNAs are small endogenous RNAs that pair and bind to sites on mRNAs to direct post-transcriptional repression. However, there is a possibility that microRNAs directly influence protein structure and activity, and this influence can be termed post-translational riboregulation. This conceptual review explores the literature on neurodegenerative disorders. Research on the association between neurodegeneration and RNA-repeat toxicity provides data that support a protein-RNA recognition code. For example, this code explains why hnRNP H and SFPQ proteins, which are involved in amyotrophic lateral sclerosis, are sequestered by the (GGGGCC)n repeat sequence. Similarly, it explains why MNBL proteins and (CTG)n repeats in RNA, which are involved in myotonic dystrophy, are sequestered into RNA foci. Using this code, proteins involved in diseases can be identified. A simple protein BLAST search of the human genome for amino acid repeats that correspond to the nucleotide repeats reveals new proteins among already known proteins that are involved in diseases. For example, the (CAG)n repeat sequence, when transcribed into possible peptide sequences, leads to the identification of PTCD3, Rem2, MESP2, SYPL2, WDR33, COL23A1, and others. After confirming this approach on RNA repeats, in the next step, the code was used in the opposite manner. Proteins that are involved in diseases were compared with microRNAs involved in those diseases. For example, a reasonable correspondence of microRNA 9 and 107 with amyloid-β-peptide (Aβ42) was identified. In the last step, a miRBase search for micro-nucleotides, obtained by transcription of a prion amino acid sequence, revealed new microRNAs and microRNAs that have previously been identified as involved in prion diseases. This concept provides a useful key for designing RNA or peptide probes.

Gain-of-Function Mutations: An Emerging Advantage for Cancer Biology

Advances in next-generation sequencing have identified thousands of genomic variants that perturb the normal functions of proteins, further contributing to diverse phenotypic consequences in cancer. Elucidating the functional pathways altered by loss-of-function (LOF) or gain-of-function (GOF) mutations will be crucial for prioritizing cancer-causing variants and their resultant therapeutic liabilities. In this review, we highlight the fundamental function of GOF mutations and discuss the potential mechanistic effects in the context of signaling networks. We also summarize advances in experimental and computational resources, which will dramatically help with studies on the functional and phenotypic consequences of mutations. Together, systematic investigations of the function of GOF mutations will provide an important missing piece for cancer biology and precision therapy.

Splicing alterations contributing to cancer hallmarks in the liver: central role of dedifferentiation and genome instability

Hepatocellular carcinoma (HCC) is a major cause of cancer-related death worldwide. HCCs are molecularly heterogeneous tumors, and this complexity is to a great extent responsible for their poor response to conventional and targeted therapies. In this review we summarize recent evidence indicating that imbalanced expression of mRNA splicing factors can be a relevant source for this heterogeneity. We also discuss how these alterations may play a driver role in hepatocarcinogenesis by impinging on the general hallmarks of cancer. Considering the natural history of HCC, we focused on two pathogenic features that are characteristic of liver tumors: chromosomal instability and phenotypic de-differentiation. We highlight mechanisms connecting splicing derangement with these two processes and the enabling capacities acquired by liver cells along their neoplastic transformation. A thorough understanding of the alterations in the splicing machinery may also help to identify new HCC biomarkers and to design novel therapeutic strategies.