A General Framework for Interrogation of mRNA Stability Programs Identifies RNA-Binding Proteins that Govern Cancer Transcriptomes

Low RNA stability signifies strong expression regulatability of tumor suppressors

RNA expression of a gene is determined by not only transcriptional regulation, but also post-transcriptional regulation of RNA decay. The precise regulation of RNA stability in the cell plays an important role in normal development. Dysregulation of RNA stability can lead to diseases such as cancer. Here we found tumor suppressor RNAs tended to decay fast in normal cell types when compared with other RNAs. Consistent with a negative effect of m6A modification on RNA stability, we observed preferential deposition of m6A on tumor suppressor RNAs. Moreover, abundant m6A and fast decay of tumor suppressor RNAs both tended to be further enhanced in prostate cancer cells relative to normal prostate epithelial cells. Further, knockdown of m6A methyltransferase METTL3 and reader YTHDF2 in prostate cancer cells both posed stronger effect on tumor suppressor RNAs than on other RNAs. These results indicated a strong post transcriptional expression regulatability mediated by abundant m6A modification on tumor suppressor RNAs.

Modulation of insulin secretion by RBFOX2-mediated alternative splicing

Insulin secretion is a tightly regulated process that is vital for maintaining blood glucose homeostasis. Although the molecular components of insulin granule trafficking and secretion are well established, how they are regulated to rapidly fine-tune secretion in response to changing environmental conditions is not well characterized. Recent studies have determined that dysregulation of RNA-binding proteins (RBPs) and aberrant mRNA splicing occurs at the onset of diabetes. We demonstrate that the RBP, RBFOX2, is a critical regulator of insulin secretion through the alternative splicing of genes required for insulin granule docking and exocytosis. Conditional mutation of Rbfox2 in the mouse pancreas results in decreased insulin secretion and impaired blood glucose homeostasis. Consistent with defects in secretion, we observe reduced insulin granule docking and corresponding splicing defects in the SNARE complex components. These findings identify an additional mechanism for modulating insulin secretion in both healthy and dysfunctional pancreatic β cells.

Surmounting Cancer Drug Resistance: New Perspective on RNA-Binding Proteins

RNA-binding proteins (RBPs), being pivotal elements in both physiological and pathological processes, possess the ability to directly impact RNA, thereby exerting a profound influence on cellular life. Furthermore, the dysregulation of RBPs not only induces alterations in the expression levels of genes associated with cancer but also impairs the occurrence of post-transcriptional regulatory mechanisms. Consequently, these circumstances can give rise to aberrations in cellular processes, ultimately resulting in alterations within the proteome. An aberrant proteome can disrupt the equilibrium between oncogenes and tumor suppressor genes, promoting cancer progression. Given their significant role in modulating gene expression and post-transcriptional regulation, directing therapeutic interventions towards RBPs represents a viable strategy for combating drug resistance in cancer treatment. RBPs possess significant potential as diagnostic and prognostic markers for diverse cancer types. Gaining comprehensive insights into the structure and functionality of RBPs, along with delving deeper into the molecular mechanisms underlying RBPs in tumor drug resistance, can enhance cancer treatment strategies and augment the prognostic outcomes for individuals afflicted with cancer.

High glucose-upregulated PD-L1 expression through RAS signaling-driven downregulation of PTRH1 leads to suppression of T cell cytotoxic function in tumor environment

BACKGROUND

Nearly 80% of patients with pancreatic cancer suffer from glucose intolerance or diabetes. Pancreatic cancer complicated by diabetes has a more immunosuppressive tumor microenvironment (TME) and is associated with a worse prognosis. The relationship between glucose metabolism and programmed cell death-Ligand 1 (PD-L1) is close and complex. It is important to explore the regulation of high glucose on PD-L1 expression in pancreatic cancer and its effect on infiltrating immune effectors in the tumor microenvironment.

METHODS

Diabetic murine models (C57BL/6) were used to reveal different immune landscape in euglycemic and hyperglycemic pancreatic tumor microenvironment. Bioinformatics, WB, iRIP [Improved RNA Binding Protein (RBP) Immunoprecipitation]-seq were used to confirm the potential regulating role of peptidyl-tRNA hydrolase 1 homolog (PTRH1) on the stability of the PD-L1 mRNA. Postoperative specimens were used to identify the expression of PD-L1 and PTRH1 in pancreatic cancer. Co-culturing T cells with pancreatic cancer cells to explore the immunosuppressive effect of pancreatic tumor cells.

RESULTS

Our results revealed that a high dose of glucose enhanced the stability of the PD-L1 mRNA in pancreatic tumor cells by downregulating PTRH1 through RAS signaling pathway activation following epidermal growth factor receptor (EGFR) stimulation. PTRH1 overexpression significantly suppressed PD-L1 expression in pancreatic cells and improved the proportion and cytotoxic function of CD8⁺ T cells in the pancreatic TME of diabetic mice.

CONCLUSIONS

PTRH1, an RBP, plays a key role in the regulation of PD-L1 by high glucose and is closely related to anti-tumor immunity in the pancreatic TME.

Crosstalk of renal cell carcinoma cells and tumor-associated macrophages aggravates tumor progression by modulating muscleblind-like protein 2/B-cell lymphoma 2/beclin 1-mediated autophagy

BACKGROUND AIMS

M2-polarized tumor-associated macrophages contribute to the development of multiple human cancers, including renal cell carcinoma (RCC). However, the crosstalk mechanism between M2 macrophages and RCC remains unclear.

METHODS

The authors constructed a co-culture system of M2 macrophages differentiated from THP-1 and RCC cells. Microscopic examination and quantitative real‑time polymerase chain reaction (qRT-PCR) validated the morphology and types of macrophages. The proliferation, migration and invasion of RCC cells were assessed by Cell Counting Kit 8 (Dojindo Molecular Technologies, Inc, Santa Clara, CA, USA) and Transwell assay (Corning, Corning, NY, USA). Messenger RNA (mRNA) and protein expression of target molecules was detected by qRT‑PCR and western blotting. Expression of Ki-67, E-cadherin and N-cadherin was measured by immunofluorescence staining or immunohistochemistry. Molecular interaction was evaluated by RNA pull-down, RNA immunoprecipitation and co-immunoprecipitation. A xenograft model was established to determine tumor growth in vivo.

RESULTS

RCC cells triggered the activation of M2 macrophages. Functionally, M2-polarized macrophages facilitated the growth, migration, invasion and epithelial-mesenchymal transition of RCC cells by suppressing autophagy, whereas rapamycin, an activator of autophagy, significantly counteracted the tumor-promoting effects of M2 macrophages. Mechanistically, M2 macrophage-derived C-C motif chemokine 2 (CCL2) enhanced modulation of muscleblind-like protein 2 (MBNL2) expression. MBNL2 raised the stability of B-cell lymphoma 2 (Bcl-2) by directly binding to Bcl-2 mRNA, which endowed RCC cells with malignant properties via inhibition of beclin 1-dependent autophagy.

CONCLUSIONS

RCC-induced M2-polarized macrophages secrete CCL2 to promote the growth and metastasis of RCC cells via inhibition of MBNL2/Bcl-2/beclin 1-mediated autophagy, which provide a novel perspective for the development of a therapeutic strategy for -RCC.

LARP6 suppresses colorectal cancer progression through ZNF267/SGMS2-mediated imbalance of sphingomyelin synthesis

BACKGROUND

With increasing incidence and mortality, colorectal cancer (CRC) seriously endangers human health. LARP6, a member of La-related protein (LARP) family, is a RNA binding protein and probably associates with CRC progression, but its specific roles and mechanisms in CRC still remain unknown.

METHOD

Quantitative real-time PCR (qPCR), western blot, and immunohistochemistry were employed to examine LARP6 expression in CRC tissues. Using the stable LARP6 overexpression or interference CRC cell lines, the effect of LARP6 on CRC progression were evaluated. High-throughput RNA immunoprecipitation sequencing (RIP-seq) and a series of relevant experiments were conducted to explain how LARP6 functions. SPSS software was used for statistical analysis.

RESULT

In this study, we found that LARP6 expression is downregulated in CRC and correlates with patients' overall survival and relapse-free survival. Furthermore, altered LARP6 expression influences CRC cells invasion and metastasis. Mechanically, we discovered that LARP6 bind ZNF267 mRNA and regulated its stability and translation. LARP6 inhibited expression of SGMS2, a downstream target of ZNF267, resulting in ceramide and sphingomyelin imbalance in CRC cells. Interestingly, LARP6 also enhances autophagy activity of CRC cells, and the effect was at least partially determined by the inhibition of SGMS2-mediated sphingomyelin synthesis.

CONCLUSION

Our study showed how LARP6/ZNF267/SGMS2 axis influence CRC progression, which contributes to further understanding of the molecular mechanisms underlying CRC development.

Spliceosomal profiling identifies EIF4A3 as a novel oncogene in hepatocellular carcinoma acting through the modulation of FGFR4 splicing

INTRODUCTION

Altered splicing landscape is an emerging cancer hallmark; however, the dysregulation and implication of the cellular machinery controlling this process (spliceosome components and splicing factors) in hepatocellular carcinoma (HCC) is poorly known. This study aimed to comprehensively characterize the spliceosomal profile and explore its role in HCC.

METHODS

Expression levels of 70 selected spliceosome components and splicing factors and clinical implications were evaluated in two retrospective and six in silico HCC cohorts. Functional, molecular and mechanistic studies were implemented in three cell lines (HepG2, Hep3B and SNU-387) and preclinical Hep3B-induced xenograft tumours.

RESULTS

Spliceosomal dysregulations were consistently found in retrospective and in silico cohorts. EIF4A3, RBM3, ESRP2 and SRPK1 were the most dysregulated spliceosome elements in HCC. EIF4A3 expression was associated with decreased survival and greater recurrence. Plasma EIF4A3 levels were significantly elevated in HCC patients. In vitro EIF4A3-silencing (or pharmacological inhibition) resulted in reduced aggressiveness, and hindered xenograft-tumours growth in vivo, whereas EIF4A3 overexpression increased tumour aggressiveness. EIF4A3-silencing altered the expression and splicing of key HCC-related genes, specially FGFR4. EIF4A3-silencing blocked the cellular response to the natural ligand of FGFR4, FGF19. Functional consequences of EIF4A3-silencing were mediated by FGFR4 splicing as the restoration of non-spliced FGFR4 full-length version blunted these effects, and FGFR4 inhibition did not exert further effects in EIF4A3-silenced cells.

CONCLUSIONS

Splicing machinery is strongly dysregulated in HCC, providing a source of new diagnostic, prognostic and therapeutic options in HCC. EIF4A3 is consistently elevated in HCC patients and associated with tumour aggressiveness and mortality, through the modulation of FGFR4 splicing.

Increased expression of the RNA-binding protein Musashi-2 is associated with immune infiltration and predicts better outcomes in ccRCC patients

RNA-binding proteins (RBPs) mainly contribute to abnormalities in posttranscriptional gene regulation. The RBP Musashi-2, an evolutionarily conserved protein, has been characterized as an oncoprotein in various tumors. However, the prognostic value and potential roles of Musashi-2 in clear cell renal cell carcinoma (ccRCC) have not yet been elucidated. In this study, we found that Musashi-2 was mainly expressed in the normal distal tubular cells and collecting duct cells of the kidneys, while its expression was significantly decreased in ccRCC. And higher expression levels of Musashi-2 indicated better overall survival (OS) in ccRCC. Furthermore, immunohistochemistry demonstrated that PD-L1 expression was negatively correlated with Musashi-2 expression, and Musashi-2 was found to be remarkably correlated with multiple immune cells and immune inhibitors, including CD8+ T cells, CD4+ T cells, regulatory T (Treg) cells, PDCD1, CTLA4, Foxp3, and LAG3. Functional enrichment analysis revealed that Musashi-2 might be involved in ccRCC metabolic reprogramming and immune infiltration and further predicted the therapeutic sensitivity of ccRCC. Taken together, Musashi-2 is a prognostic biomarker for ccRCC patients that may provide novel insights into individualized treatment strategies and guide effective immunotherapy.

Pan-cancer analysis of mRNA stability for decoding tumour post-transcriptional programs

Measuring mRNA decay in tumours is a prohibitive challenge, limiting our ability to map the post-transcriptional programs of cancer. Here, using a statistical framework to decouple transcriptional and post-transcriptional effects in RNA-seq data, we uncover the mRNA stability changes that accompany tumour development and progression. Analysis of 7760 samples across 18 cancer types suggests that mRNA stability changes are ~30% as frequent as transcriptional events, highlighting their widespread role in shaping the tumour transcriptome. Dysregulation of programs associated with >80 RNA-binding proteins (RBPs) and microRNAs (miRNAs) drive these changes, including multi-cancer inactivation of RBFOX and miR-29 families. Phenotypic activation or inhibition of RBFOX1 highlights its role in calcium signaling dysregulation, while modulation of miR-29 shows its impact on extracellular matrix organization and stemness genes. Overall, our study underlines the integral role of mRNA stability in shaping the cancer transcriptome, and provides a resource for systematic interrogation of cancer-associated stability pathways.

The role of mRNA stability in shaping the cancer transcriptome is revealed using a statistical analysis of transcriptomic data.

Construction and validation of a prognostic model of RNA binding proteins in clear cell renal carcinoma

BACKGROUND

The dysfunction of RNA binding proteins (RBPs) is associated with various inflammation and cancer. The occurrence and progression of tumors are closely related to the abnormal expression of RBPs. There are few studies on RBPs in clear cell renal carcinoma (ccRCC), which allows us to explore the role of RBPs in ccRCC.

METHODS

We obtained the gene expression data and clinical data of ccRCC from the Cancer Genome Atlas (TCGA) database and extracted all the information of RBPs. We performed differential expression analysis of RBPs. Risk model were constructed based on the differentially expressed RBPs (DERBPs). The expression levels of model markers were examined by reverse transcription-quantitative PCR (RT-qPCR) and analyzed for model-clinical relevance. Finally, we mapped the model's nomograms to predict the 1, 3 and 5-year survival rates for ccRCC patients.

RESULTS

The results showed that the five-year survival rate for the high-risk group was 40.2% (95% CI = 0.313 ~ 0.518), while the five-year survival rate for the low-risk group was 84.3% (95% CI = 0.767 ~ 0.926). The ROC curves (AUC = 0.748) also showed that our model had stable predictive power. Further RT-qPCR results were in accordance with our analysis (p < 0.05). The results of the independent prognostic analysis showed that the model could be an independent prognostic factor for ccRCC. The results of the correlation analysis also demonstrated the good predictive ability of the model.

CONCLUSION

In summary, the 4-RBPs (EZH2, RPL22L1, RNASE2, U2AF1L4) risk model could be used as a prognostic indicator of ccRCC. Our study provides a possibility for predicting the survival of ccRCC.

Pumilio protects Xbp1 mRNA from regulated Ire1-dependent decay

The unfolded protein response (UPR) maintains homeostasis of the endoplasmic reticulum (ER). Residing in the ER membrane, the UPR mediator Ire1 deploys its cytoplasmic kinase-endoribonuclease domain to activate the key UPR transcription factor Xbp1 through non-conventional splicing of Xbp1 mRNA. Ire1 also degrades diverse ER-targeted mRNAs through regulated Ire1-dependent decay (RIDD), but how it spares Xbp1 mRNA from this decay is unknown. Here, we identify binding sites for the RNA-binding protein Pumilio in the 3′UTR Drosophila Xbp1. In the developing Drosophila eye, Pumilio binds both the Xbp1^unspliced and Xbp1^spliced mRNAs, but only Xbp1^spliced is stabilized by Pumilio. Furthermore, Pumilio displays Ire1 kinase-dependent phosphorylation during ER stress, which is required for its stabilization of Xbp1^spliced. hIRE1 can phosphorylate Pumilio directly, and phosphorylated Pumilio protects Xbp1^spliced mRNA against RIDD. Thus, Ire1-mediated phosphorylation enables Pumilio to shield Xbp1^spliced from RIDD. These results uncover an unexpected regulatory link between an RNA-binding protein and the UPR.

In Drosophila, ER-targeted mRNAs are degraded by Ire1-dependent pathway. Here the authors report that the fly mRNA binding protein Pumilio is phosphorylated by Ire1 and binds to Xbp1 mRNA, protecting it from the non-canonical endoribonuclease activity of Ire1.

Distribution of alternative untranslated regions within the mRNA of the CELF1 splicing factor affects its expression

CUG-binding protein, ELAV-like Family Member 1 (CELF1) plays an important role during the development of different tissues, such as striated muscle and brain tissue. CELF1 is an RNA-binding protein that regulates RNA metabolism processes, e.g., alternative splicing, and antagonizes other RNA-binding proteins, such as Muscleblind-like proteins (MBNLs). Abnormal activity of both classes of proteins plays a crucial role in the pathogenesis of myotonic dystrophy type 1 (DM1), the most common form of muscular dystrophy in adults. In this work, we show that alternative splicing of exons forming both the 5' and 3' untranslated regions (UTRs) of CELF1 mRNA is efficiently regulated during development and tissue differentiation and is disrupted in skeletal muscles in the context of DM1. Alternative splicing of the CELF1 5'UTR leads to translation of two potential protein isoforms that differ in the lengths of their N-terminal domains. We also show that the MBNL and CELF proteins regulate the distribution of mRNA splicing isoforms with different 5'UTRs and 3'UTRs and affect the CELF1 expression by changing its sensitivity to specific microRNAs or RNA-binding proteins. Together, our findings show the existence of different mechanisms of regulation of CELF1 expression through the distribution of various 5' and 3' UTR isoforms within CELF1 mRNA.

Codon optimality in cancer

A key characteristic of cancer cells is their increased proliferative capacity, which requires elevated levels of protein synthesis. The process of protein synthesis involves the translation of codons within the mRNA coding sequence into a string of amino acids to form a polypeptide chain. As most amino acids are encoded by multiple codons, the nucleotide sequence of a coding region can vary dramatically without altering the polypeptide sequence of the encoded protein. Although mutations that do not alter the final amino acid sequence are often thought of as silent/synonymous, these can still have dramatic effects on protein output. Because each codon has a distinct translation elongation rate and can differentially impact mRNA stability, each codon has a different degree of 'optimality' for protein synthesis. Recent data demonstrates that the codon preference of a transcriptome matches the abundance of tRNAs within the cell and that this supply and demand between tRNAs and mRNAs varies between different cell types. The largest observed distinction is between mRNAs encoding proteins associated with proliferation or differentiation. Nevertheless, precisely how codon optimality and tRNA expression levels regulate cell fate decisions and their role in malignancy is not fully understood. This review describes the current mechanistic understanding on codon optimality, its role in malignancy and discusses the potential to target codon optimality therapeutically in the context of cancer.

RNA-binding Protein MBNL2 regulates Cancer Cell Metastasis through MiR-182-MBNL2-AKT Pathway

The aberrant expression of RNA-binding proteins (RBPs) plays important roles in the occurrence and progression of cancer. MBNL2 is a member of the RNA binding protein MBNL family that is widely expressed in mammalian cells. We report here that MBNL2 is downregulated in breast, lung and liver cancer tissues, the promoter methylation levels of MBNL2 are higher in cancer tissues than normal tissues. The enrichment analysis of MBNL2 correlated genes indicates the potential function of MBNL2 on cancer progression. MBNL2 regulates cancer cell migration and invasion by modulating PI3K/AKT-mediated epithelial-mesenchymal transition. PI3K/AKT inhibitor overcomes the promotive effect of shMBNL2 on metastasis. The expression of MBNL2 is directly targeted by miR-182. miR-182 is upregulated in breast, lung and liver cancers and has good potential for cancer diagnosis. miR-182 promotes cancer cell migration and invasion by inhibiting the expression of MBNL2. Re-introduction of exogenous MBNL2 reverses the promotive effect of miR-182 on metastasis. Collectively, these findings suggest that MBNL2 plays a tumor suppressive function through miR-182-MBNL2-AKT-EMT signaling pathways.

Novel insights into the molecular mechanisms underlying risk of colorectal cancer from smoking and red/processed meat carcinogens by modeling exposure in normal colon organoids

Tobacco smoke and red/processed meats are well-known risk factors for colorectal cancer (CRC). Most research has focused on studies of normal colon biopsies in epidemiologic studies or treatment of CRC cell lines in vitro. These studies are often constrained by challenges with accuracy of self-report data or, in the case of CRC cell lines, small sample sizes and lack of relationship to normal tissue at risk. In an attempt to address some of these limitations, we performed a 24-hour treatment of a representative carcinogens cocktail in 37 independent organoid lines derived from normal colon biopsies. Machine learning algorithms were applied to bulk RNA-sequencing and revealed cellular composition changes in colon organoids. We identified 738 differentially expressed genes in response to carcinogens exposure. Network analysis identified significantly different modules of co-expression, that included genes related to MSI-H tumor biology, and genes previously implicated in CRC through genome-wide association studies. Our study helps to better define the molecular effects of representative carcinogens from smoking and red/processed meat in normal colon epithelial cells and in the etiology of the MSI-H subtype of CRC, and suggests an overlap between molecular mechanisms involved in inherited and environmental CRC risk.

lncRNA SNHG6 promotes hepatocellular carcinoma progression by interacting with HNRNPL/PTBP1 to facilitate SETD7/LZTFL1 mRNA destabilization

The lncRNA SNHG6 (small nucleolar RNA host gene 6) plays vital roles in tumorigenesis and the progression of hepatocellular carcinoma (HCC). However, the regulatory mechanisms of SNHG6 are largely unknown. In this study, we identified, via quantitative proteomics, specific cytoskeleton-associated proteins and enzyme modulators to be potential targets of SNHG6. SNHG6 reduced the mRNA levels of lysine methyltransferase, SET domain containing 7 (SETD7) and leucine zipper transcription factor-like 1 (LZTFL1) by posttranscriptional destabilization. Silencing of SETD7 or LZTFL1 reversed the suppressive effects of SNHG6 knockdown on HCC progression. Heterogeneous nuclear ribonucleoprotein L (HNRNPL) and polypyrimidine tract binding protein 1 (PTBP1) were identified as SNHG6-interacting proteins that bind to SETD7 or LZTFL1 mRNA. Forced expression of SNHG6 led to HNRNPL being competitively adsorbed by SNHG6, thereby removing its stabilizing effect on SETD7. Concurrently, the functional SNHG6-PTBP1 complex facilitated the degradation of LZTFL1 mRNA in hepatoma cells. These results indicated that SNHG6 promotes HCC progression by functioning as a "decoy plus guide" for HNRNPL and PTBP1 to facilitate mRNA decay of SETD7 and LZTFL1, thereby serving as a novel therapeutic target for HCC.

Systematic analysis of the function and prognostic value of RNA binding protein in head and neck squamous cell carcinoma

OBJECTIVE

Dysregulation of RNA binding proteins (RBPs) plays an important role in controlling processes in cancer development. However, the function of RBPs has not been thoroughly and systematically documented in head and neck cancer. We aim to explore the role of RPB in the pathogenesis of HNSC.

METHODS

We obtained HNSC gene expression data and corresponding clinical information from The Cancer Genome Atlas (TCGA) and the GEO databases, and identified aberrantly expressed RBPs between tumors and normal tissues. Meanwhile, we performed a series of bioinformatics to explore the function and prognostic value of these RBPs.

RESULTS

A total of 249 abnormally expressed RBPs were identified, including 101 downregulated RBPs and 148 upregulated RBPs. Using least absolute shrinkage and selection operator (LASSO) and univariate Cox regression analysis, the 15 RPBs were identified as hub genes. With the 15 RPBS, the prognostic prediction model was constructed. Further analysis showed that the high-risk score of the patients expressed a better survival outcome. The prediction model was validated in another HNSC dataset, and similar findings were observed.

CONCLUSIONS

Our results provide novel insights into the pathogenesis of HNSC. The fifteen RBP gene signature exhibited the predictive value of moderate HNSC prognosis, and have potential application value in clinical decision-making and individualized treatment.

Integrated analysis of the functions and prognostic values of RNA binding proteins in hepatocellular carcinoma

Background

Hepatocellular carcinoma (HCC), one of the most common malignant tumors worldwide, ranks as the fifth most common cancer and has been the second most frequent cause of cancer-related death. RNA binding proteins (RBPs) are proteins that interact with different classes of RNA and are commonly detected in cells.

Methods

We used RNA sequencing data from TCGA to display dysfunctional RBPs microenvironments and provide potential useful biomarkers for HCC diagnosis and prognosis.

Results

330 differently expressed RBPs (208 upregulated and 122 downregulated) were identified. KEGG were mainly enriched in RNA degradation, Influenza A, Hepatitis C, RIG-I-like receptor signaling pathway, Herpes simplex virus 1 infection and RNA transport. CBioPortal results demonstrated that these genes were altered in 50 samples out of 357 HCC patients (14%) and the amplification of BRCA1 was the largest frequent copy-number alteration.

Conclusion

Based on the online database, we identified novel RBPs markers for the prognosis of hepatocellular carcinoma.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12876-021-01843-0.

Transite: A Computational Motif-Based Analysis Platform That Identifies RNA-Binding Proteins Modulating Changes in Gene Expression

RNA-binding proteins (RBPs) play critical roles in regulating gene expression by modulating splicing, RNA stability, and protein translation. Stimulus-induced alterations in RBP function contribute to global changes in gene expression, but identifying which RBPs are responsible for the observed changes remains an unmet need. Here, we present Transite, a computational approach that systematically infers RBPs influencing gene expression through changes in RNA stability and degradation. As a proof of principle, we apply Transite to RNA expression data from human patients with non-small-cell lung cancer whose tumors were sampled at diagnosis or after recurrence following treatment with platinum-based chemotherapy. Transite implicates known RBP regulators of the DNA damage response and identifies hnRNPC as a new modulator of chemotherapeutic resistance, which we subsequently validated experimentally. Transite serves as a framework for the identification of RBPs that drive cell-state transitions and adds additional value to the vast collection of publicly available gene expression datasets.

Krismer et al. present a computational approach to identify RNA-binding proteins (RBPs) that modulate post-transcriptional control of gene expression using RNA expression data as inputs. By applying this approach to publicly available patient datasets, they identify and experimentally confirm that the RBP hnRNPC contributes to chemotherapy resistance in lung cancer.

Integrated Analysis of Nine Prognostic RNA-Binding Proteins in Soft Tissue Sarcoma

RNA-binding proteins (RBPs) have been shown to be dysregulated in cancer transcription and translation, but few studies have investigated their mechanism of action in soft tissue sarcoma (STS). Here, The Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) databases were used to identify differentially expressed RBPs in STS and normal tissues. Through a series of biological information analyses, 329 differentially expressed RBPs were identified. Functional enrichment analysis showed that differentially expressed RBPs were mainly involved in RNA transport, RNA splicing, mRNA monitoring pathways, ribosome biogenesis and translation regulation. Through Cox regression analyses, 9 RBPs (BYSL, IGF2BP3, DNMT3B, TERT, CD3EAP, SRSF12, TLR7, TRIM21 and MEX3A) were all up-regulated in STS as prognosis-related genes, and a prognostic model was established. The model calculated a risk score based on the expression of 9 hub RBPs. The risk score could be used for risk stratification of patients and had a high prognostic value based on the receiver operating characteristic (ROC) curve. We also established a nomogram containing risk scores and 9 key RBPs to predict the 1-year, 3-year, and 5-year survival rates of patients in STS. Afterwards, methylation analysis showed significant changes in the methylation degree of BYSL, CD3EAP and MEX2A. Furthermore, the expression of 9 hub RBPs was closely related to immune infiltration rather than tumor purity. Based on the above studies, these findings may provide new insights into the pathogenesis of STS and will provide candidate biomarkers for the prognosis of STS.

Integrated analysis of the prognostic values of RNA-binding proteins in head and neck squamous cell carcinoma

Head and neck squamous cell carcinoma is a malignant tumor of the upper aerodigestive tract. These RNA-binding proteins (RBPs) influence post-transcriptional in cells and regulate cell physiology, participate in regulating RNA stability, alternative splicing, translation, modification, localization, and apoptosis. We used RNA sequencing data from The Cancer Genome Atlas to display dysfunctional RBPs microenvironments and provide potential useful biomarkers for head and neck squamous cell carcinoma (HNSCC) diagnosis and prognosis. Six RBPs (DNMT1, PCF11, EIF5A2, RNASE10, PSMA6, and IGF2BP2) were selected as independent prognosis factors of HNSCC patients. The Kyoto Encyclopedia of Genes and Genomes were mainly enriched in RNA transport, Spliceosome, RNA degradation, mRNA surveillance pathway, and Epstein-Barr virus infection. cBioPortal results demonstrated that these six genes were altered in 150 samples out of 504 HNSCC patients (30%) and the amplification of IGF2BP2 was the largest frequent copy-number alteration. Based on the online database, we identified novel RBPs markers for the prognosis of HNSCC.

RNA-Binding Proteins Hold Key Roles in Function, Dysfunction, and Disease

RNA-binding proteins (RBPs) are multi-faceted proteins in the regulation of RNA or its RNA splicing, localisation, stability, and translation. Amassing proof from many recent and dedicated studies reinforces the perception of RBPs exerting control through differing expression levels, cellular localization and post-transcriptional alterations. However, since the regulation of RBPs is reliant on the micro-environment and events like stress response and metabolism, their binding affinities and the resulting RNA-RBP networks may be affected. Therefore, any misregulation and disruption in the features of RNA and its related homeostasis can lead to a number of diseases that include diabetes, cardiovascular disease, and other disorders such as cancer and neurodegenerative diseases. As such, correct regulation of RNA and RBPs is crucial to good health as the effect RBPs exert through loss of function can cause pathogenesis. In this review, we will discuss the significance of RBPs and their typical function and how this can be disrupted in disease.

NFIB promotes the migration and progression of kidney renal clear cell carcinoma by regulating PINK1 transcription

Kidney renal clear cell carcinoma (KIRC) is the most common and aggressive type of renal cell carcinoma. Due to high mortality rate, high metastasis rate and chemical resistance, the prognosis of KIRC patients is poor. Therefore, it is necessary to study the mechanisms of KIRC development and to develop more effective prognostic molecular biomarkers to help clinical patients. In our study, we used The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases to investigate that the expression of nuclear factor I B (NFIB) is significantly higher in KIRC than in adjacent tissues. Moreover, NFIB expression levels are associated with multiple clinical pathological parameters of KIRC, and KIRC patients with high NFIB expression have poor prognosis, suggesting that NFIB may play vital roles in the malignant development of KIRC. Further studies demonstrated that NFIB could promote the progression and metastasis of KIRC and participate in the regulation of PTEN induced kinase 1 (PINK1). Furthermore, we used chromatin immunoprecipitation (ChIP) experiments to confirm that NFIB binds to the PINK1 promoter and regulates its expression at the transcriptional level. Further experiments also confirmed the important roles of PINK1 in promoting the development of tumors by NFIB. Hence, our data provide a new NFIB-mediated regulatory mechanism for the tumor progression of KIRC and suggest that NFIB can be applied as a new predictor and therapeutic target for KIRC.

A Prognostic Model Based on RNA Binding Protein Predicts Clinical Outcomes in Hepatocellular Carcinoma Patients

Dysregulation of RNA binding proteins (RBPs) is closely associated with tumor events. However, the function of RBPs in hepatocellular carcinoma (HCC) has not been fully elucidated. The RNA sequences and relevant clinical data of HCC were retrieved from the The Cancer Genome Atlas (TCGA) database to identify distinct RBPs. Subsequently, univariate and multivariate cox regression analysis was performed to evaluate the overall survival (OS)-associated RBPs. The expression levels of prognostic RBP genes and survival information were analyzed using a series of bioinformatics tool. A total of 365 samples with 1,542 RBPs were included in this study. One hundred and eighty-seven differently RBPs were screened, including 175 up-regulated and 12 down-regulated. The independent OS-associated RBPs of NHP2, UPF3B, and SMG5 were used to develop a prognostic model. Survival analysis showed that low-risk patients had a significantly longer OS and disease-free survival (DFS) when compared to high-risk patients (HR: 2.577, 95% CI: 1.793-3.704, P < 0.001 and HR: 1.599, 95% CI: 1.185-2.159, P = 0.001, respectively). The International Cancer Genome Consortium (ICGC) database was used to externally validate the model, and the OS of low-risk patients were found to be longer than that of high-risk patients (P < 0.001). The Nomograms of OS and DFS were plotted to help in clinical decision making. These results showed that the model was effective and may help in prognostic stratification of HCC patients. The prognostic prediction model based on RBPs provides new insights for HCC diagnosis and personalized treatment.

NFIB promotes the progression of gastric cancer by upregulating circMAP7D1 to stabilize HER2 mRNA

The present study evaluated the expression levels of nuclear factor I B (NFIB) in gastric cancer (GC) specimens and cells, and its regulatory roles were further elucidated. The expression levels of NFIB were examined in GC and paired normal specimens, and in human GC and normal gastric epithelial cells by reverse transcription‑quantitative PCR. A circular RNA (circRNA) microarray was performed to identify the novel downstream circRNA of NFIB. Cell proliferation was determined by Cell Counting Kit‑8 assay. Furthermore, cell cycle distribution and apoptosis were assessed using flow cytometry. Interactions between RNA were examined by RNA pulldown assay and the stability of target mRNA was evaluated using a mRNA stability assay. The results of the present study revealed that NFIB was upregulated in GC. Furthermore, silencing NFIB suppressed the proliferation of GC cells, whereas cell cycle arrest and apoptosis were enhanced. In addition, significant downregulation of circMAP7D1 (hsa_circ_0004093) was observed in GC cells infected with short hairpin RNA‑NFIB. These findings indicated that circMAP7D1 may be a promising downstream molecule of NFIB in GC, and further functional analyses indicated that circMAP7D1 was involved in NFIB‑modulated GC cell proliferation and apoptosis. Moreover, human epidermal growth factor receptor 2 (HER2) was identified as a novel target of circMAP7D1 in GC, and NFIB was able to increase the stability of HER2 mRNA through regulating circMAP7D1. In conclusion, the present findings indicated that NFIB expression was increased in GC. In addition, NFIB may promote the proliferation of GC cells and function through stabilizing HER2 mRNA by upregulating circMAP7D1. Notably, NFIB and its novel downstream signaling pathway may serve essential roles during the development of GC, and NFIB may be considered a promising candidate for the treatment of patients with GC.

HuR/ELAVL1 drives malignant peripheral nerve sheath tumor growth and metastasis

Cancer cells can develop a strong addiction to discrete molecular regulators, which control the aberrant gene expression programs that drive and maintain the cancer phenotype. Here, we report the identification of the RNA-binding protein HuR/ELAVL1 as a central oncogenic driver for malignant peripheral nerve sheath tumors (MPNSTs), which are highly aggressive sarcomas that originate from cells of the Schwann cell lineage. HuR was found to be highly elevated and bound to a multitude of cancer-associated transcripts in human MPNST samples. Accordingly, genetic and pharmacological inhibition of HuR had potent cytostatic and cytotoxic effects on tumor growth, and strongly suppressed metastatic capacity in vivo. Importantly, we linked the profound tumorigenic function of HuR to its ability to simultaneously regulate multiple essential oncogenic pathways in MPNST cells, including the Wnt/β-catenin, YAP/TAZ, RB/E2F, and BET pathways, which converge on key transcriptional networks. Given the exceptional dependency of MPNST cells on HuR for survival, proliferation, and dissemination, we propose that HuR represents a promising therapeutic target for MPNST treatment.

An APEX2 proximity ligation method for mapping interactions with the nuclear lamina

The nuclear lamina (NL) is a meshwork found beneath the inner nuclear membrane. The study of the NL is hindered by the insolubility of the meshwork and has driven the development of proximity ligation methods to identify the NL-associated/proximal proteins, RNA, and DNA. To simplify and improve temporal labeling, we fused APEX2 to the NL protein lamin-B1 to map proteins, RNA, and DNA. The identified NL-interacting/proximal RNAs show a long 3' UTR bias, a finding consistent with an observed bias toward longer 3' UTRs in genes deregulated in lamin-null cells. A C-rich motif was identified in these 3' UTR. Our APEX2-based proteomics identifies a C-rich motif binding regulatory protein that exhibits altered localization in lamin-null cells. Finally, we use APEX2 to map lamina-associated domains (LADs) during the cell cycle and uncover short, H3K27me3-rich variable LADs. Thus, the APEX2-based tools presented here permit identification of proteomes, transcriptomes, and genome elements associated with or proximal to the NL.

SFPQ promotes an oncogenic transcriptomic state in melanoma

The multifunctional protein, splicing factor, proline- and glutamine-rich (SFPQ) has been implicated in numerous cancers often due to interaction with coding and non-coding RNAs, however, its role in melanoma remains unclear. We report that knockdown of SFPQ expression in melanoma cells decelerates several cancer-associated cell phenotypes, including cell growth, migration, epithelial to mesenchymal transition, apoptosis, and glycolysis. RIP-seq analysis revealed that the SFPQ-RNA interactome is reprogrammed in melanoma cells and specifically enriched with key melanoma-associated coding and long non-coding transcripts, including SOX10, AMIGO2 and LINC00511 and in most cases SFPQ is required for the efficient expression of these genes. Functional analysis of two SFPQ-enriched lncRNA, LINC00511 and LINC01234, demonstrated that these genes independently contribute to the melanoma phenotype and a more detailed analysis of LINC00511 indicated that this occurs in part via modulation of the miR-625-5p/PKM2 axis. Importantly, analysis of a large clinical cohort revealed that elevated expression of SFPQ in primary melanoma tumours may have utility as a prognostic biomarker. Together, these data suggest that SFPQ is an important driver of melanoma, likely due to SFPQ-RNA interactions promoting the expression of numerous oncogenic transcripts.

RNA-binding protein SORBS2 suppresses clear cell renal cell carcinoma metastasis by enhancing MTUS1 mRNA stability

RNA-binding proteins (RBPs) predominantly contribute to abnormal posttranscriptional gene modulation and disease progression in cancer. Sorbin and SH3 domain-containing 2 (SORBS2), an RBP, has been reported to be a potent tumor suppressor in several cancer types. Through integrative analysis of clinical specimens, we disclosed that the expression level of SORBS2 was saliently decreased in metastatic tissues and positively correlated with overall survival. We observed that overexpression of SORBS2 brought about decreased metastatic capacity in ccRCC cell lines. Transcriptome-wide analysis revealed that SORBS2 notably increased microtubule-associated tumor-suppressor 1 gene (MTUS1) expression. In-depth mechanistic exploring discovered that the Cys2-His2 zinc finger (C2H2-ZnF) domain of SORBS2 directly bound to the 3′ untranslated region (3′UTR) of MTUS1 mRNA, which increased MTUS1 mRNA stability. In addition, we identified that MTUS1 regulated microtubule dynamics via promoting KIF2C^S192 phosphorylation by Aurora B. Together, our research identified SORBS2 as a suppressor of ccRCC metastasis by enhancing MTUS1 mRNA stability, providing a novel understanding of RBPs during ccRCC progression.

Integrated analysis of the functions of RNA binding proteins in clear cell renal cell carcinoma

RNA binding proteins (RBPs) dysregulation is involved in the processes of various tumors. However, the roles of RBPs in clear cell renal cell carcinoma (ccRCC) remain poorly understand. In present study, we first performed consensus clustering and identified two clusters, of which cluster 2 was closely correlated with the malignancy of ccRCC. Differentially expressed RBPs between normal and tumor tissues were obtained, comprising 71 up-regulated and 44 down-regulated ones. Then, ten hub genes were selected and validated using The Human Protein Atlas database and receiver operating characteristic curves, showing good diagnostic value for cancers. Besides, we identified ten RBPs with the most useful prognostic values, and were used to construct a risk score model. The model could be used to stratify patients with different prognosis and phenotype distributions. The model showed good performance and can be used as a complementation for clinical factors to guide clinical practice in the future.

Effects of RNA Binding Proteins on the Prognosis and Malignant Progression in Prostate Cancer

Prostate cancer (PCa) is a common lethal malignancy in men. RNA binding proteins (RBPs) have been proven to regulate the biological processes of various tumors, but their roles in PCa remain less defined. In the present study, we used bioinformatics analysis to identify RBP genes with prognostic and diagnostic values. A total of 59 differentially expressed RBPs in PCa were obtained, comprising 28 upregulated and 31 downregulated RBP genes, which may play important roles in PCa. Functional enrichment analyses showed that these RBPs were mainly involved in mRNA processing, RNA splicing, and regulation of RNA splicing. Additionally, we identified nine RBP genes (EXO1, PABPC1L, REXO2, MBNL2, MSI1, CTU1, MAEL, YBX2, and ESRP2) and their prognostic values by a protein-protein interaction network and Cox regression analyses. The expression of these nine RBPs was validated using immunohistochemical staining between the tumor and normal samples. Further, the associations between the expression of these nine RBPs and pathological T staging, Gleason score, and lymph node metastasis were evaluated. Moreover, these nine RBP genes showed good diagnostic values and could categorize the PCa patients into two clusters with different malignant phenotypes. Finally, we constructed a prognostic model based on these nine RBP genes and validated them using three external datasets. The model showed good efficiency in predicting patient survival and was independent of other clinical factors. Therefore, our model could be used as a supplement for clinical factors to predict patient prognosis and thereby improve patient survival.

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.

RNA-Binding Proteins in Cancer: Functional and Therapeutic Perspectives

Simple Summary

RNA-binding proteins (RBPs) play central roles in regulating posttranscriptional expression of genes. Many of them are known to be deregulated in a wide variety of cancers. Dysregulated RBPs influence the expression levels of target RNAs related to cancer phenotypes, such as proliferation, apoptosis, angiogenesis, senescence, and EMT/invasion/metastasis. Thus, understanding the molecular functions of RBPs and their roles in cancer-related phenotypes can lead to improved therapeutic strategies.

Abstract

RNA-binding proteins (RBPs) crucially regulate gene expression through post-transcriptional regulation, such as by modulating microRNA (miRNA) processing and the alternative splicing, alternative polyadenylation, subcellular localization, stability, and translation of RNAs. More than 1500 RBPs have been identified to date, and many of them are known to be deregulated in cancer. Alterations in the expression and localization of RBPs can influence the expression levels of oncogenes, tumor-suppressor genes, and genome stability-related genes. RBP-mediated gene regulation can lead to diverse cancer-related cellular phenotypes, such as proliferation, apoptosis, angiogenesis, senescence, and epithelial-mesenchymal transition (EMT)/invasion/metastasis. This regulation can also be associated with cancer prognosis. Thus, RBPs can be potential targets for the development of therapeutics for the cancer treatment. In this review, we describe the molecular functions of RBPs, their roles in cancer-related cellular phenotypes, and various approaches that may be used to target RBPs for cancer treatment.

Integrated analysis of RNA-binding proteins in human colorectal cancer

Background

Although RNA-binding proteins play an essential role in a variety of different tumours, there are still limited efforts made to systematically analyse the role of RNA-binding proteins (RBPs) in the survival of colorectal cancer (CRC) patients.

Methods

Analysis of CRC transcriptome data collected from the TCGA database was conducted, and RBPs were extracted from CRC. R software was applied to analyse the differentially expressed genes (DEGs) of RBPs. To identify related pathways and perform functional annotation of RBP DEGs, Gene Ontology (GO) function and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were carried out using the database for annotation, visualization and integrated discovery. Protein-protein interactions (PPIs) of these DEGs were analysed based on the Search Tool for the Retrieval of Interacting Genes (STRING) database and visualized by Cytoscape software. Based on the Cox regression analysis of the prognostic value of RBPs (from the PPI network) with survival time, the RBPs related to survival were identified, and a prognostic model was constructed. To verify the model, the data stored in the TCGA database were designated as the training set, while the chip data obtained from the GEO database were treated as the test set. Then, both survival analysis and ROC curve verification were conducted. Finally, the risk curves and nomograms of the two groups were generated to predict the survival period.

Results

Among RBP DEGs, 314 genes were upregulated while 155 were downregulated, of which twelve RBPs (NOP14, MRPS23, MAK16, TDRD6, POP1, TDRD5, TDRD7, PPARGC1A, LIN28B, CELF4, LRRFIP2, MSI2) with prognostic value were obtained.

Conclusions

The twelve identified genes may be promising predictors of CRC and play an essential role in the pathogenesis of CRC. However, further investigation of the underlying mechanism is needed.

Comprehensive Analysis of the Roles and Prognostic Value of RNA-Binding Proteins in Head and Neck Squamous Cell Carcinoma

Head and neck squamous cell carcinoma (HNSCC) is a malignancy with relatively high incidence and poor prognosis. RNA-binding proteins (RBPs) were reported to be dysregulated in multiple cancers and were closely associated with tumor initiation and progression. However, an integrated analysis of the roles of RBPs in HNSCC has not been conducted. In the present study, we obtained transcriptome data and corresponding clinical information of HNSCC patients from The Cancer Genome Atlas database and screened out differentially expressed RBPs between tumor and normal tissues. Subsequently, we utilized a series of bioinformatics analyses to elucidate the potential functions and prognostic value of these RBPs in HNSCC. As a result, a total of 88 aberrantly expressed RBPs were identified, including 63 downregulated and 25 upregulated RBPs. Functional enrichment analysis suggested that the differentially expressed RBPs mainly participated in mRNA metabolic processes, RNA processing, RNA transport, regulation of RNA stability, RNA degradation, and mRNA surveillance pathway. Three RBP genes (NOVA1, EZH2, and RBM24) were determined as prognosis-related hub genes from which EZH2 and NOVA1 were selected to construct a prognostic signature based on LASSO Cox regression algorithm. Further analysis demonstrated that the high-risk patient group stratified by the risk signature has advanced tumor grade and poorer overall survival when compared with low-risk group. Moreover, univariate analysis showed that the risk score, tumor stage, T stage, and N stage were significantly associated with patient overall survival and the multivariate analysis results indicated that the risk score and age were greatly correlated with patient prognosis. Overall, this study provided a comprehensive landscape of RBPs in HNSCC and identified an effective gene signature for predicting the clinical outcomes of HNSCC patient, which may contribute to clinical decision making and individualized cancer treatment.

Intrinsic Regulatory Role of RNA Structural Arrangement in Alternative Splicing Control

Alternative splicing is a highly sophisticated process, playing a significant role in posttranscriptional gene expression and underlying the diversity and complexity of organisms. Its regulation is multilayered, including an intrinsic role of RNA structural arrangement which undergoes time- and tissue-specific alterations. In this review, we describe the principles of RNA structural arrangement and briefly decipher its cis- and trans-acting cellular modulators which serve as crucial determinants of biological functionality of the RNA structure. Subsequently, we engage in a discussion about the RNA structure-mediated mechanisms of alternative splicing regulation. On one hand, the impairment of formation of optimal RNA structures may have critical consequences for the splicing outcome and further contribute to understanding the pathomechanism of severe disorders. On the other hand, the structural aspects of RNA became significant features taken into consideration in the endeavor of finding potential therapeutic treatments. Both aspects have been addressed by us emphasizing the importance of ongoing studies in both fields.

Identification of prognostic alternative splicing signatures in hepatitis B or/and C viruses related hepatocellular carcinoma

BACKGROUND

Alternative splicing (AS) takes a crucial part in tumor process. We aim to analyze AS in Hepatitis B virus (HBV) or/and hepatitis C virus (HCV) related hepatocellular carcinoma (HCC).

METHODS

Cox regression analysis was conducted to screen survival-associated AS events. The receiver operating characteristic curve used to evaluate the predictive accuracy. Splicing network was built to investigate the relationship between splicing factors and AS events.

RESULTS

Ninety-six survival-associated AS events were obtained by univariate Cox regression. Final prognostic model could significantly distinguish the prognosis. We identified RBFOX2 as the hub gene in splicing network based on differentially expressed splicing factors, and obtained MAP3K13_AT as the key AS event in survival-related splicing network.

CONCLUSION

Our results highlight the AS signatures in HCC patients with HBV or/and HCV infection. Meanwhile, AS events and splicing factors in different virus-infected HCC subgroups can provide novel perspectives as biomarkers and individualized therapeutic targets.

Muscleblind-like 2 controls the hypoxia response of cancer cells

Hypoxia is a hallmark of solid cancers, supporting proliferation, angiogenesis, and escape from apoptosis. There is still limited understanding of how cancer cells adapt to hypoxic conditions and survive. We analyzed transcriptome changes of human lung and breast cancer cells under chronic hypoxia. Hypoxia induced highly concordant changes in transcript abundance, but divergent splicing responses, underlining the cell type-specificity of alternative splicing programs. While RNA-binding proteins were predominantly reduced, hypoxia specifically induced muscleblind-like protein 2 (MBNL2). Strikingly, MBNL2 induction was critical for hypoxia adaptation by controlling the transcript abundance of hypoxia response genes, such as vascular endothelial growth factor A (VEGFA) MBNL2 depletion reduced the proliferation and migration of cancer cells, demonstrating an important role of MBNL2 as cancer driver. Hypoxia control is specific for MBNL2 and not shared by its paralog MBNL1. Thus, our study revealed MBNL2 as central mediator of cancer cell responses to hypoxia, regulating the expression and alternative splicing of hypoxia-induced genes.

Integrated Analysis of RNA-Binding Proteins in Glioma

RNA-binding proteins (RBPs) play important roles in many cancer types. However, RBPs have not been thoroughly and systematically studied in gliomas. Global analysis of the functional impact of RBPs will provide a better understanding of gliomagenesis and new insights into glioma therapy. In this study, we integrated a list of the human RBPs from six sources-Gerstberger, SONAR, Gene Ontology project, Poly(A) binding protein, CARIC, and XRNAX-which covered 4127 proteins with RNA-binding activity. The RNA sequencing data were downloaded from The Cancer Genome Atlas (TCGA) (n = 699) and Chinese Glioma Genome Atlas (CGGA) (n = 325 + 693). We examined the differentially expressed genes (DEGs) using the R package DESeq2, and constructed a weighted gene co-expression network analysis (WGCNA) of RBPs. Furthermore, survival analysis was also performed based on the univariate and multivariate Cox proportional hazards regression models. In the WGCNA analysis, we identified a key module involved in the overall survival (OS) of glioblastomas. Survival analysis revealed eight RBPs (PTRF, FNDC3B, SLC25A43, ZC3H12A, LRRFIP1, HSP90B1, HSPA5, and BNC2) are significantly associated with the survival of glioblastoma patients. Another 693 patients within the CGGA database were used to validate the findings. Additionally, 3564 RBPs were classified into canonical and non-canonical RBPs depending on the domains that they contain, and non-canonical RBPs account for the majority (72.95%). The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed that some non-canonical RBPs may have functions in glioma. Finally, we found that the knockdown of non-canonical RBPs, PTRF, or FNDC3B can alone significantly inhibit the proliferation of LN229 and U251 cells. Simultaneously, RNA Immunoprecipitation (RIP) analysis indicated that PTRF may regulate cell growth and death- related pathways to maintain tumor cell growth. In conclusion, our findings presented an integrated view to assess the potential death risks of glioblastoma at a molecular level, based on the expression of RBPs. More importantly, we identified non-canonical RNA-binding proteins PTRF and FNDC3B, showing them to be potential prognostic biomarkers for glioblastoma.

Characterizing the interplay between gene nucleotide composition bias and splicing

BACKGROUND

Nucleotide composition bias plays an important role in the 1D and 3D organization of the human genome. Here, we investigate the potential interplay between nucleotide composition bias and the regulation of exon recognition during splicing.

RESULTS

By analyzing dozens of RNA-seq datasets, we identify two groups of splicing factors that activate either about 3200 GC-rich exons or about 4000 AT-rich exons. We show that splicing factor-dependent GC-rich exons have predicted RNA secondary structures at 5' ss and are dependent on U1 snRNP-associated proteins. In contrast, splicing factor-dependent AT-rich exons have a large number of decoy branch points, SF1- or U2AF2-binding sites and are dependent on U2 snRNP-associated proteins. Nucleotide composition bias also influences local chromatin organization, with consequences for exon recognition during splicing. Interestingly, the GC content of exons correlates with that of their hosting genes, isochores, and topologically associated domains.

CONCLUSIONS

We propose that regional nucleotide composition bias over several dozens of kilobase pairs leaves a local footprint at the exon level and induces constraints during splicing that can be alleviated by local chromatin organization at the DNA level and recruitment of specific splicing factors at the RNA level. Therefore, nucleotide composition bias establishes a direct link between genome organization and local regulatory processes, like alternative splicing.

CCT3 acts upstream of YAP and TFCP2 as a potential target and tumour biomarker in liver cancer

Although Yes-associated protein (YAP) is very important to liver cancer, its nuclear localisation prevents consideration as a promising therapeutic target and a diagnostic biomarker. Recently, we reported that the protumourigenic roles of YAP in liver cancer are indispensable for transcription factor CP2 (TFCP2) in a Hippo-independent manner; however, proteins that act upstream to simultaneously control YAP and TFCP2 remain unclear. The aim of this study was to uncover such proteins and evaluate whether they are potential YAP-associated therapeutic targets and diagnostic biomarkers. Mass spectrometry revealed that chaperonin containing TCP1 subunit 3 (CCT3) co-interact with YAP and TFCP2, and notably, CCT3 is a non-nuclear protein. CCT3 was elevated in liver cancer, and its higher expression was associated with poorer overall survival. Inhibiting CCT3 resulted in a suppressed transformative phenotype in liver cancer cells, suggesting that CCT3 might be a potential therapeutic target. CCT3 prolonged half-life of YAP and TFCP2 by blocking their ubiquitination caused by poly(rC) binding protein 2 (PCBP2) in a beta-transducin repeat containing E3 ubiquitin protein ligase (βTrCP)-independent manner. Interestingly, PCBP2 directly interacted with YAP via a WB motif-WW domain interaction, whereas indirectly interacted with TFCP2 via the aid of YAP. Furthermore, CCT3 was capable of separating PCBP2-YAP interactions, thereby preventing YAP and TFCP2 from PCBP2-induced ubiquitination. Moreover, YAP and TFCP2 were downstream of CCT3 to positively control tumourigenesis, yet such effects were inhibited by PCBP2. Clinically, CCT3 was positively correlated with YAP and TFCP2, and elevated levels of the CCT3-YAP-TFCP2 axis might be critical for liver malignancy. In addition, seral-CCT3 was proven to be a potential biomarker, and its diagnostic capacity was better than that of alpha fetoprotein (AFP) to a certain extent. Together, CCT3 acts as a trigger of YAP and TFCP2 to affect tumourigenesis and serves as a potential therapeutic target and biomarker in liver cancer.