Alternative Splicing of EZH2 pre-mRNA by SF3B3 Contributes to the Tumorigenic Potential of Renal Cancer

Landscape analysis of alternative splicing in kidney renal clear cell carcinoma and their clinical significance

A growing number of studies reveal that alternative splicing (AS) is associated with tumorigenesis, progression, and metastasis. Systematic analysis of alternative splicing signatures in renal cancer is lacking. In our study, we investigated the AS landscape of kidney renal clear cell carcinoma (KIRC) and identified AS predictive model to improve the prognostic prediction of KIRC. We obtained clinical data and gene expression profiles of KIRC patients from the TCGA database to evaluate AS events. The calculation results for seven types of AS events indicated that 46276 AS events from 10577 genes were identified. Next, we applied Cox regression analysis to identify 5864 prognostic-associated AS events. We used the Metascape database to verify the potential pathways of prognostic-associated AS. Moreover, we constructed KIRC prediction systems with prognostic-associated AS events by the LASSO Cox regression model. AUCs demonstrated that these prediction systems had excellent prognostic accuracy simultaneously. We identified 34 prognostic associated splicing factors (SFs) and constructed homologous regulatory networks. Furthermore, in vitro experiments were performed to validate the favorable effect of SFs FMR1 in KIRC. In conclusion, we overviewed AS events in KIRC and identified AS-based prognostic models to assist the survival prediction of KIRC patients. Our study may provide a novel predictive signature to improve the prognostic prediction of KIRC, which might facilitate KIRC patient counseling and individualized management.

The Many Roads from Alternative Splicing to Cancer: Molecular Mechanisms Involving Driver Genes

Cancer driver genes are either oncogenes or tumour suppressor genes that are classically activated or inactivated, respectively, by driver mutations. Alternative splicing-which produces various mature mRNAs and, eventually, protein variants from a single gene-may also result in driving neoplastic transformation because of the different and often opposed functions of the variants of driver genes. The present review analyses the different alternative splicing events that result in driving neoplastic transformation, with an emphasis on their molecular mechanisms. To do this, we collected a list of 568 gene drivers of cancer and revised the literature to select those involved in the alternative splicing of other genes as well as those in which its pre-mRNA is subject to alternative splicing, with the result, in both cases, of producing an oncogenic isoform. Thirty-one genes fall into the first category, which includes splicing factors and components of the spliceosome and splicing regulators. In the second category, namely that comprising driver genes in which alternative splicing produces the oncogenic isoform, 168 genes were found. Then, we grouped them according to the molecular mechanisms responsible for alternative splicing yielding oncogenic isoforms, namely, mutations in cis splicing-determining elements, other causes involving non-mutated cis elements, changes in splicing factors, and epigenetic and chromatin-related changes. The data given in the present review substantiate the idea that aberrant splicing may regulate the activation of proto-oncogenes or inactivation of tumour suppressor genes and details on the mechanisms involved are given for more than 40 driver genes.

Epigenetic regulation of FOXI2 promotes clear cell renal cell carcinoma progression

In recent decades, substantial advancements in epigenetics have unveiled a profound understanding of its mechanisms in tumorigenesis and have offered promising strategies for epigenetic therapy in cancer patients. In our study, through bioinformatics analysis, we discovered a significant downregulation and hypermethylation of FOXI2 in clear cell renal cell carcinoma (ccRCC), while the expression in chromophobe cell carcinoma (chRCC) exhibited the opposite trend. Moreover, we established a strong correlation between FOXI2 expression levels and the prognosis of ccRCC. Gene enrichment analysis and cell function experiments unequivocally demonstrate that FOXI2 possesses the capability to induce cell cycle arrest and inhibit cell proliferation. Our research findings demonstrate that the expression of FOXI2 in ccRCC is under the regulation of promoter hypermethylation. Furthermore, in vitro experiments have conclusively shown that the overexpression of FOXI2 induces cell cycle arrest and inhibits cell proliferation.

SF3B3-regulated mTOR alternative splicing promotes colorectal cancer progression and metastasis

BACKGROUND

Aberrant alternative splicing (AS) is a pervasive event during colorectal cancer (CRC) development. SF3B3 is a splicing factor component of U2 small nuclear ribonucleoproteins which are crucial for early stages of spliceosome assembly. The role of SF3B3 in CRC remains unknown.

METHODS

SF3B3 expression in human CRCs was analyzed using publicly available CRC datasets, immunohistochemistry, qRT-PCR, and western blot. RNA-seq, RNA immunoprecipitation, and lipidomics were performed in SF3B3 knockdown or overexpressing CRC cell lines. CRC cell xenografts, patient-derived xenografts, patient-derived organoids, and orthotopic metastasis mouse models were utilized to determine the in vivo role of SF3B3 in CRC progression and metastasis.

RESULTS

SF3B3 was upregulated in CRC samples and associated with poor survival. Inhibition of SF3B3 by RNA silencing suppressed the proliferation and metastasis of CRC cells in vitro and in vivo, characterized by mitochondria injury, increased reactive oxygen species (ROS), and apoptosis. Mechanistically, silencing of SF3B3 increased mTOR exon-skipped splicing, leading to the suppression of lipogenesis via mTOR-SREBF1-FASN signaling. The combination of SF3B3 shRNAs and mTOR inhibitors showed synergistic antitumor activity in patient-derived CRC organoids and xenografts. Importantly, we identified SF3B3 as a critical regulator of mTOR splicing and autophagy in multiple cancers.

CONCLUSIONS

Our findings revealed that SF3B3 promoted CRC progression and metastasis by regulating mTOR alternative splicing and SREBF1-FASN-mediated lipogenesis, providing strong evidence to support SF3B3 as a druggable target for CRC therapy.

Identification of splicing factors signature predicting prognosis risk and the mechanistic roles of novel oncogenes in HNSCC

Head and neck squamous cell carcinoma (HNSCC) is the most frequent subtype of head and neck cancer, generally with a poor prognosis and limited therapeutic options due to its highly heterogeneous malignancy. In this study, we screened functional splicing regulatory RNA binding proteins (RBPs) that were closely related with the prognosis of HNSCC patients and showed significant expression differences between HNSCC tumors and normal tissues. Based on this finding, we chose six candidate genes (HNRNPC, ZCRB1, RBM12B, SF3A2, SF3B3, and SRSF11) to generate a prognostic prediction model and validated the accuracy of the prognostic model for predicting patient survival outcomes. We found that the risk score predicted by our model can serve as an independent prognostic predictor. Notably, HNSCC tumors showing higher expression of SF3B3, HNRNPC, or ZCRB1 possessed higher risk scores in the discovered prediction model. The investigation of the underlying mechanism validated that knockdown of SF3B3, HNRNPC, and ZCRB1 separately induced a substantial impairment of HNSCC cell survival. Conversely, overexpression of each of the three genes promoted tumor cellular proliferation. High throughput RNA sequencing analysis revealed that changes in the expression of SF3B3 and HNRNPC remarkably affected alternative splicing of genes related to cell cycle regulation, whereas the depletion of ZCRB1 contributed to aberrant splicing events involving in DNA damage response. In addition, the prognostic prediction model's risk score was demonstrated to be related with the immune infiltration score. Particularly, SF3B3 has a negative correlation with CD8A expression. Therefore, our findings provide promising prognosis predictors and potential therapeutic targets for better treatment efficacy of HNSCC.

Alternative splicing decouples local from global PRC2 activity

The Polycomb repressive complex 2 (PRC2) mediates epigenetic maintenance of gene silencing in eukaryotes via methylation of histone H3 at lysine 27 (H3K27). Accessory factors define two distinct subtypes, PRC2.1 and PRC2.2, with different actions and chromatin-targeting mechanisms. The mechanisms orchestrating PRC2 assembly are not fully understood. Here, we report that alternative splicing (AS) of PRC2 core component SUZ12 generates an uncharacterized isoform SUZ12-S, which co-exists with the canonical SUZ12-L isoform in virtually all tissues and developmental stages. SUZ12-S drives PRC2.1 formation and favors PRC2 dimerization. While SUZ12-S is necessary and sufficient for the repression of target genes via promoter-proximal H3K27me3 deposition, SUZ12-L maintains global H3K27 methylation levels. Mouse embryonic stem cells (ESCs) lacking either isoform exit pluripotency more slowly and fail to acquire neuronal cell identity. Our findings reveal a physiological mechanism regulating PRC2 assembly and higher-order interactions in eutherians, with impacts on H3K27 methylation and gene repression.

An intricate rewiring of cancer metabolism via alternative splicing

All human genes undergo alternative splicing leading to the diversity of the proteins. However, in some cases, abnormal regulation of alternative splicing can result in diseases that trigger defects in metabolism, reduced apoptosis, increased proliferation, and progression in almost all tumor types. Metabolic dysregulations and immune dysfunctions are crucial factors in cancer. In this respect, alternative splicing in tumors could be a potential target for therapeutic cancer strategies. Dysregulation of alternative splicing during mRNA maturation promotes carcinogenesis and drug resistance in many cancer types. Alternative splicing (changing the target mRNA 3'UTR binding site) can result in a protein with altered drug affinity, ultimately leading to drug resistance.. Here, we will highlight the function of various alternative splicing factors, how it regulates the reprogramming of cancer cell metabolism, and their contribution to tumor initiation and proliferation. Also, we will discuss emerging therapeutics for treating tumors via abnormal alternative splicing. Finally, we will discuss the challenges associated with these therapeutic strategies for clinical applications.

Adverse developmental impacts in progeny of zebrafish exposed to the agricultural herbicide atrazine during embryogenesis

Atrazine (ATZ) is an herbicide commonly used on crops in the Midwestern US and other select global regions. The US Environmental Protection Agency ATZ regulatory limit is 3 parts per billion (ppb; µg/L), but this limit is often exceeded. ATZ has a long half-life, is a common contaminant of drinking water sources, and is indicated as an endocrine disrupting chemical in multiple species. The zebrafish was used to test the hypothesis that an embryonic parental ATZ exposure alters protein levels leading to modifications in morphology and behavior in developing progeny. Zebrafish embryos (F1) were collected from adults (F0) exposed to 0, 0.3, 3, or 30 ppb ATZ during embryogenesis. Differential proteomics, morphology, and behavior assays were completed with offspring aged 120 or 144 h with no additional chemical treatment. Proteomic analysis identified differential expression of proteins associated with neurological development and disease; and organ and organismal morphology, development, and injury, specifically the skeletomuscular system. Head length and ratio of head length to total length was significantly increased in the F1 of 0.3 and 30 ppb ATZ groups (p < 0.05). Based on molecular pathway alterations, further craniofacial morphology assessment found decreased distance for cartilaginous structures, decreased surface area and distance between saccular otoliths, and a more posteriorly positioned notochord (p < 0.05), indicating delayed ossification and skeletal growth. The visual motor response assay showed hyperactivity in progeny of the 30 ppb treatment group for distance moved and of the 0.3 and 30 ppb treatment groups for time spent moving (p < 0.05). Due to the changes in saccular otoliths, an acoustic startle assay was completed and showed decreased response in the 0.3 and 30 ppb treatments (p < 0.05). These findings suggest that a single embryonic parental exposure alters cellular pathways in their progeny that lead to perturbations in craniofacial development and behavior.

ROCK2-RNA interaction map reveals multiple biological mechanisms underlying tumor progression in renal cell carcinoma

Renal cell carcinoma (RCC) is the most common form of kidney cancer in adults. Despite new therapeutic modalities, the outcomes for RCC patients remain unsatisfactory. Rho-associated coiled-coil forming protein kinase 2 (ROCK2) has previously been shown to be upregulated in RCC, and its expression was negatively correlated with patient survival. However, the precise molecular function of ROCK2 has remained unclear. Herein, using RNA-seq analysis of ROCK2 knockdown and control cells, we identified 464 differentially expressed genes, and 1287 alternative splicing events in 786-O RCC cells. Furthermore, mapping of iRIP-seq reads in 786-O cells showed a biased distribution at 5' UTR, intronic and intergenic regions. By comparing ROCK2-regulated alternative splicing and iRIP-seq data, we found 292 overlapping genes that are enriched in multiple tumorigenic pathways. Taken together, our work defined a complex ROCK2-RNA interaction map on a genomic scale in a human RCC cell line, which deepens our understanding of the molecular function of ROCK2 in cancer development.

The interplay between lncRNA NR_030777 and SF3B3 in neuronal damage caused by paraquat

Paraquat (PQ) has been widely acknowledged as an environmental risk factor for Parkinson's disease (PD). However, the interaction between splicing factor and long non-coding RNA (lncRNA) in the process of PQ-induced PD has rarely been studied. Based on previous research, this study focused on splicing factor 3 subunit 3 (SF3B3) and lncRNA NR_030777. After changing the target gene expression level by lentiviral transfection technology, the related gene expression was detected by western blot and qRT-PCR. The expression of SF3B3 protein was reduced in Neuro-2a cells after PQ exposure, and the reactive oxygen species (ROS) scavenger N-acetylcysteine prevented this decline. Knockdown of SF3B3 reduced the PQ-triggered NR_030777 expression increase, and overexpression of NR_030777 reduced the transcriptional and translational level of Sf3b3. Then, knockdown of SF3B3 exacerbated the PQ-induced decrease in cell viability and aggravated the reduction of tyrosine hydroxylase (TH) protein expression. Overexpressing SF3B3 reversed the reduction of TH expression caused by PQ. Moreover, after intervention with the autophagy inhibitor Bafilomycin A1, LC3B-II protein expression was further increased in Neuro-2a cells with the knockdown of SF3B3, indicating that autophagy was enhanced. In conclusion, PQ modulated the interplay between NR_030777 and SF3B3 through ROS production, thereby impairing autophagic flux and causing neuronal damage.

Multi-omics analysis of the Indian ovarian cancer cohort revealed histotype-specific mutation and gene expression patterns

Introduction: In India, OVCa is women’s third most common and lethal cancer type, accounting for 6.7% of observed cancer incidences. The contribution of somatic mutations, aberrant expression of gene and splice forms in determining the cell fate, gene networks, tumour-specific variants, and the role of immune fraction infiltration have been proven essential in understanding tumorigenesis. However, their interplay in OVCa in a histotype-specific manner remains unclear in the Indian context. In the present study, we aimed to unravel the Indian population histotype-specific exome variants, differentially expressed gene modules, splice events and immune profiles of OVCa samples.

Methods: We analysed 10 tumour samples across 4 ovarian cancer histotypes along with 2 normal patient samples. This included BCFtool utilities and CNVkit for exome, WGCNA and DESeq2 for obtaining differential module hub genes and dysregulated miRNA targets, CIBERSORTx for individual immune profiles and rMATS for tumour specific splice variants.

Result: We identified population-specific novel mutations in Cancer Gene Census Tier1 and Tier2 genes. MUC16, MUC4, CIITA, and NCOR2 were among the most mutated genes, along with TP53. Transcriptome analysis showed significant overexpression of mutated genes MUC16, MUC4, and CIITA, whereas NCOR2 was downregulated. WGCNA revealed histotype-specific gene hubs and networks. Among the significant pathways, alteration in the immune system was one of the pathways, and immune profiling using CIBERSORTx revealed histotype-specific immune cell fraction. miRNA analysis revealed miR-200 family, miR-200a and miR-429 were upregulated in HGSOCs.Splice factor abrasion caused splicing perturbations, with the most abundant alternative splice event being exon skipping and the most spliced gene, SNHG17. Pathway analysis of spliced genes revealed translational elongation and Base excision repair as the pathways altered in OVCa.

Conclusion: Integrated exome, transcriptome, and splicing patterns revealed different population-specific molecular signatures of ovarian cancer histotypes in the Indian Cohort.

Novel long noncoding RNA LINC02820 augments TNF signaling pathway to remodel cytoskeleton and potentiate metastasis in esophageal squamous cell carcinoma

Esophageal squamous cell carcinoma (ESCC) is one of the most common malignant tumors in China. However, there are no targets to treat ESCC because the molecular mechanism behind the cancer is still unclear. Here, we found a novel long noncoding RNA LINC02820 was upregulated in ESCC and associated with the ESCC clinicopathological stage. Through a series of functional experiments, we observed that LINC02820 only promoted the migration and invasion capabilities of ESCC cell lines. Mechanically, we found that LINC02820 may affect the cytoskeletal remodeling, interact with splice factor 3B subunit 3 (SF3B3), and cooperate with TNFα to amplify the NF-κB signaling pathway, which can lead to ESCC metastasis. Overall, our findings revealed that LINC02820 is a potential biomarker and therapeutic target for the diagnosis and treatment of ESCC.

Splicing factor SF3B3, a NS5-binding protein, restricts ZIKV infection by targeting GCH1

Zika virus (ZIKV), a positive-sense single-stranded RNA virus, causes congenital ZIKV syndrome in children and Guillain-Barré Syndrome (GBS) in adults. ZIKV expresses nonstructural protein 5 (NS5), a large protein that is essential for viral replication. ZIKV NS5 confers the ability to evade interferon (IFN) signalling; however, the exact mechanism remains unclear. In this study, we employed affinity pull-down and liquid chromatography-tandem mass spectrometry (LC-MS/MS) analyses and found that splicing factor 3b subunit 3 (SF3B3) is associated with the NS5-Flag pull-down complex through interaction with NS5. Functional assays showed that SF3B3 overexpression inhibited ZIKV replication by promoting IFN-stimulated gene (ISG) expression whereas silencing of SF3B3 inhibited expression of ISGs to promote ZIKV replication. GTP cyclohydrolase I (GCH1) is the first and rate-limiting enzyme in tetrahydrobiopterin (BH4) biosynthesis. NS5 upregulates the expression of GCH1 during ZIKV infection. And GCH1 marginally promoted ZIKV replication via the IFN pathway. Additionally, GCH1 expression is related to the regulation of SF3B3. Overexpression of the SF3B3 protein effectively reduced GCH1 protein levels, whereas SF3B3 knockdown increased its levels. These findings indicated that ZIKV NS5 binding protein SF3B3 contributed to the host immune response against ZIKV replication by modulating the expression of GCH1.

The Role of Alternative Splicing in Cancer: Regulatory Mechanism, Therapeutic Strategy, and Bioinformatics Application

[Formula: see text] Alternative splicing (AS) can generate distinct transcripts and subsequent isoforms that play differential functions from the same pre-mRNA. Recently, increasing numbers of studies have emerged, unmasking the association between AS and cancer. In this review, we arranged AS events that are closely related to cancer progression and presented promising treatments based on AS for cancer therapy. Obtaining proliferative capacity, acquiring invasive properties, gaining angiogenic features, shifting metabolic ability, and getting immune escape inclination are all splicing events involved in biological processes. Spliceosome-targeted and antisense oligonucleotide technologies are two novel strategies that are hopeful in tumor therapy. In addition, bioinformatics applications based on AS were summarized for better prediction and elucidation of regulatory routines mingled in. Together, we aimed to provide a better understanding of complicated AS events associated with cancer biology and reveal AS a promising target of cancer treatment in the future.

The splicing factor SF3B4 drives proliferation and invasion in cervical cancer by regulating SPAG5

Regulation of alternative splicing (AS) by the splicing factor 3b (SF3B) family plays an essential role in cancer. However, the biological function of SF3B family members in cervical cancer (CC) needs to be further elucidated. In this study, we found that splicing factor 3b subunit 4 (SF3B4) was highly expressed in CC by bioinformatics analysis using cervical squamous cell carcinoma and endocervical adenocarcinoma (CESC) data from The Cancer Genome Atlas (TCGA). Then, we demonstrated that high expression of SF3B4 promoted proliferation and invasion abilities of CC cells in vitro and in vivo and that reduced expression of SF3B4 performed the opposite effect. Further RNA-seq and AS analysis showed that sperm-associated antigen 5 (SPAG5) was a downstream target gene of SF3B4. Interestingly, SPAG5 expression was decreased after SF3B4 knockdown because of retained introns (RIs) and reduced maturation of SPAG5 pre-mRNA. Importantly, SPAG5 deficiency impaired the oncogenic effects of SF3B4 overexpression on CC cells. In conclusion, SF3B4 promotes CC progression by regulating the effective splicing of SPAG5. SF3B4 could be a promising target for CC.

RNA binding protein RALY activates the cholesterol synthesis pathway through an MTA1 splicing switch in hepatocellular carcinoma

Alternative splicing is an important RNA processing event that contributes to RNA complexity and protein diversity in cancer. Accumulating evidence demonstrates the essential roles of some alternatively spliced genes in carcinogenesis. However, the potential roles of alternatively spliced genes in hepatocellular carcinoma (HCC) are still largely unknown. Here we showed that the HnRNP Associated with Lethal Yellow Protein Homolog (RALY) gene is upregulated and associated with poor outcomes in HCC patients. RALY acts as a tumor-promoting factor by cooperating with splicing factor 3b subunit 3 (SF3B3) and modulating the splicing switch of Metastasis Associated 1 (MTA1) from MTA-S to MTA1-L. Normally, MTA1-S inhibits cell proliferation by reducing the transcription of cholesterol synthesis genes. In HCC, RALY and SF3B3 cooperate to regulate the MTA1 splicing switch, leading to a reduction in the MTA1-S level, and alleviating the inhibitory effect of MTA1-S on cholesterol synthesis genes, thus promoting HCC cell proliferation. In conclusion, our results revealed that the RALY-SF3B3/MTA1/cholesterol synthesis pathway contributes essentially to hepatic carcinogenesis and could serve as a promising therapeutic target for HCC.

ASpediaFI: Functional Interaction Analysis of Alternative Splicing Events

Alternative splicing (AS) regulates biological processes governing phenotypes and diseases. Differential AS (DAS) gene test methods have been developed to investigate important exonic expression from high-throughput datasets. However, the DAS events extracted using statistical tests are insufficient to delineate relevant biological processes. In this study, we developed a novel application, Alternative Splicing Encyclopedia: Functional Interaction (ASpediaFI), to systemically identify DAS events and co-regulated genes and pathways. ASpediaFI establishes a heterogeneous interaction network of genes and their feature nodes (i.e., AS events and pathways) connected by co-expression or pathway gene set knowledge. Next, ASpediaFI explores the interaction network using the random walk with restart algorithm and interrogates the proximity from a query gene set. Finally, ASpediaFI extracts significant AS events, genes, and pathways. To evaluate the performance of our method, we simulated RNA sequencing (RNA-seq) datasets to consider various conditions of sequencing depth and sample size. The performance was compared with that of other methods. Additionally, we analyzed three public datasets of cancer patients or cell lines to evaluate how well ASpediaFI detects biologically relevant candidates. ASpediaFI exhibits strong performance in both simulated and public datasets. Our integrative approach reveals that DAS events that recognize a global co-expression network and relevant pathways determine the functional importance of spliced genes in the subnetwork. ASpediaFI is publicly available at https://bioconductor.org/packages/ASpediaFI.

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.

The exon 12-containing LHX6 isoforms promote cervical cancer cell proliferation by regulating the MAPK signaling pathway

LIM homeobox 6 (LHX6) has been reported to be downregulated and inhibits cell proliferation in various cancers. Alternative splicing of LHX6 leads to six annotated isoforms, which can be found in the NCBI database. However, the expression patterns and potential roles of these isoforms remain poorly characterized in cervical cancer. Here, we demonstrated that the LHX6 isoforms containing exon 12 (LHX6^EX(+12) group) and isoforms lacking exon 12 (LHX6^EX(-12) group) were differentially expressed in cervical tissue by qRT-PCR. The mRNA expression level of LHX6^EX(+12) group was higher than that of LHX6^EX(-12) group in cervical cancer tissue. Knockdown of LHX6^EX(+12) group and all LHX6 isoforms (LHX6^All group) inhibited cell growth, increased cell apoptosis, and induced cell cycle arrest from G0/G1 phase to S phase in vitro. Consistently, overexpression of the LHX6^EX(+12) group promoted cervical cancer cell proliferation in vitro. In contrast, no significant differences in cell proliferation were found between LHX6^EX(-12) isoform knockdown group and its control. RNA-sequencing suggested that the LHX6^EX(+12) isoform group might exert its cancer-promoting effects in cervical cancer via regulating MAPK signaling pathway. Downregulation of the LHX6^EX(+12) group significantly suppressed the phosphorylation of MRK, ERK, JNK, and P38 at the protein level. We also identified some unique biological processes and signaling pathways in which each isoform group might be involved. In summary, our results indicated that LHX6^EX(+12) isoform group was the dominant oncogenic type of LHX6 in cervical cancer, which may be a new biomarker and a potential precise therapeutic target for cervical cancer in the future.

Key pituitary miRNAs mediate the expression of pig GHRHR splice variants by regulating splice factors

The growth hormone releasing hormone receptor (GHRHR) is well documented in organism growth and its alternative splicing may generate multiple functional GHRHR splice variants (SVs). Our previous study has demonstrated the key pituitary miRNAs (let-7e and miR-328-5p) in pig regulated the expression of GHRHR SVs by directly targeting to them. And according to recent reports, the interplay between miRNA-based silencing of mRNAs and alternative splicing of pre-mRNAs is a crucial post-transcriptional mechanism. In this study, SF3B3 and CPSF4 were firstly excavated as the splice factors that involved in the formation of GHRHR SVs mediated by let-7e and miR-328-5p through the comparation of the expression relations of GHRHR SVs, let-7e/miR-328-5p and SF3B3/CPSF4 in pituitary tissues between Landrace pigs and BaMa pigs, as well as the prediction of the target relations of let-7e/miR-328-5p with SF3B3 and/or CPSF4. SF3B3 and CPSF4 targeted by let-7e and miR-328-5p were further verified by performing dual-luciferase reporter assays and detecting the expression of target transcripts. Then the RT-PCR, RT-qPCR and Western blot assays were used to confirm SF3B3 and CPSF4 were involved in the formation of the GHRHR SVs, and in this process, let-7e and miR-328-5p mediated GHRHR SVs by regulating SF3B3 and CPSF4. Finally, the target site of SF3B3 on pre-GHRHR was on the Exon 12 to Exon14, while CPSF4 acted on the other fragments of the pre-GHRHR, which were explored by dual-luciferase reporter system preliminarily. To the best of our knowledge, this paper is the first to report the miRNAs regulate GHRHR SVs indirectly by splice factors.

SF3B4 promotes ovarian cancer progression by regulating alternative splicing of RAD52

Many studies have proven that splicing factors are crucial for human malignant tumor development. However, as a classical splicing factor, the expression of SF3B4 is not clear, and its biological function needs to be further clarified in ovarian cancer (OC). We determined that SF3B4 was obviously upregulated and its high expression was associated with poor prognosis in OC patients. In vitro and in vivo assays suggested that SF3B4 overexpression promoted OC cell proliferation and mobility, and downregulation of SF3B4 had the opposite effect. Further studies found that miR-509–3p decreased SF3B4 mRNA expression by binding to the 3’ -UTR of SF3B4 directly. Importantly, we revealed that RAD52 was a potential target of SF3B4 through alternative splicing events analysis. Loss of SF3B4 led to decreased expression of RAD52, owing to intron 8 retention and generation of premature termination codons. Moreover, decreased expression of RAD52 partially counteracted the tumor-promoting effect of SF3B4 overexpression. In conclusion, our results suggested that SF3B4, negatively regulated by miR-509–3p, promoted OC progression through effective splicing of RAD52. Therefore, SF3B4 may be a promising biomarker and effective therapeutic target for OC.

Alternative Splicing, Epigenetic Modifications and Cancer: A Dangerous Triangle, or a Hopeful One?

Simple Summary

Epigenetics studies the alteration of gene expression without changing DNA sequence and very often, epigenetic dysregulation causes cancer. Alternative splicing is a mechanism that results in the production of several mRNA isoforms from a single gene and aberrant splicing is also a frequent cause of cancer. The present review is built on the interrelations of epigenetics and alternative splicing. In an intuitive way, we say that epigenetic modifications and alternative splicing are at two vertices of a triangle, the third vertex being occupied by cancer. Interconnection between alternative splicing and epigenetic modifications occurs backward and forward and the mechanisms involved are widely reviewed. These connections also provide novel diagnostic or prognostic tools, which are listed. Finally, as epigenetic alterations are reversible and aberrant alternative splicing may be corrected, the therapeutic possibilities to break the triangle are discussed.

Abstract

The alteration of epigenetic modifications often causes cancer onset and development. In a similar way, aberrant alternative splicing may result in oncogenic products. These issues have often been individually reviewed, but there is a growing body of evidence for the interconnection of both causes of cancer. Actually, aberrant splicing may result from abnormal epigenetic signalization and epigenetic factors may be altered by alternative splicing. In this way, the interrelation between epigenetic marks and alternative splicing form the base of a triangle, while cancer may be placed at the vertex. The present review centers on the interconnections at the triangle base, i.e., between alternative splicing and epigenetic modifications, which may result in neoplastic transformations. The effects of different epigenetic factors, including DNA and histone modifications, the binding of non-coding RNAs and the alterations of chromatin organization on alternative splicing resulting in cancer are first considered. Other less-frequently considered questions, such as the epigenetic regulation of the splicing machinery, the aberrant splicing of epigenetic writers, readers and erasers, etc., are next reviewed in their connection with cancer. The knowledge of the above-mentioned relationships has allowed increasing the collection of biomarkers potentially useful as cancer diagnostic and/or prognostic tools. Finally, taking into account on one hand that epigenetic changes are reversible, and some epigenetic drugs already exist and, on the other hand, that drugs intended for reversing aberrations in alternative splicing, therapeutic possibilities for breaking the mentioned cancer-related triangle are discussed.

Depleted HDAC3 attenuates hyperuricemia-induced renal interstitial fibrosis via miR-19b-3p/SF3B3 axis

Dysfunctional histone deacetylases (HDACs) elicit unrestrained fibrosis and damage to organs. With regard to the link between HDACs and fibrosis, this research is practiced to decipher the concrete mechanism of HDAC3 in hyperuricemia (HN)-induced renal interstitial fibrosis (RIF) from microRNA-19b-3p/splicing factor 3b subunit 3 (miR-19b-3p/SF3B3) axis.The HN model was established on rats to induce RIF by oral administration of adenine and potassium oxalate. HN rats were injected with miR-19b-3p- or HDAC3-related vectors to figure out their effects on RIF through detecting 24-h urine protein, uric acid (UA), blood urea nitrogen (BUN) and serum creatinine (Scr) contents and α-smooth muscle actin (α-SMA), transforming growth factor β1 (TGF-β1) and fibronectin (FN) contents in renal tissues and observing pathological damages and RIF index of renal tissues. HDAC3, miR-19b-3p and SF3B3 expression in renal tissues were tested, along with their interactions.Elevated HDAC3 and SF3B3 and reduced miR-19b-3p were displayed in renal tissues of HN rats. Suppressed HDAC3 or promoted miR-19b-3p relieved HN-induced RIF, as reflected by their inhibitory effects on 24 h urine protein, UA, BUN, Scr, α-SMA, TGF-β1, and FN contents and RIF index and their ameliorated effects on pathological damages of renal tissues. HDAC3 bound to the promoter of miR-19b-3p to regulate SF3B3. MiR-19b-3p depletion abrogated down-regulated HDAC3-induced effects on HN-induced RIF.It is delineated that depressed HDAC3 relives HN-induced RIF through restoring miR-19b-3p and knocking down SF3B3, replenishing the references for RIF curing.

Systematic analyses of the role of prognostic and immunological EIF3A, a reader protein, in clear cell renal cell carcinoma

Background

Eukaryotic initiation factor 3a (EIF3A), a “reader” protein for RNA methylation, has been found to be involved in promoting tumorigenesis in a variety of cancers. The impact of EIF3A in clear cell renal cell carcinoma (ccRCC) has yet to be reported. This study aimed to identify the prognostic value of EIF3A in ccRCC and investigate the relationship between EIF3A expression and immune infiltration.

Methods

We collected 29 m6A-related mRNA data and clinicopathological parameters from The Cancer Genome Atlas (TCGA) database. Logistic regression analyses were used to analyse the correlation between EIF3A expression and clinical characteristics. Immunohistochemistry (IHC) was applied to examine EIF3A levels in normal and ccRCC tissues. Univariate and multivariate analyses were conducted to recognize independent factors associated with overall survival (OS) and disease-free survival (DFS). The nomogram aimed to predict the 1-, 3- and 5-year survival probabilities. Gene set enrichment analysis (GSEA) was carried out to determine the potential functions and related signalling pathways of EIF3A expression. To investigate EIF3A of coexpressed genes, we used LinkedOmics, and the results were subjected to enrichment analysis. Simultaneously, LinkedOmics and STRING datasets were used to identify EIF3A coexpressed genes that were visualized via Cytoscape. Finally, we evaluated whether EIF3A expression correlated with the infiltration of immune cells and the expression of marker genes in ccRCC by Tumour Immune Estimation Resource (TIMER) and Gene Expression Profiling Interactive Analysis (GEPIA).

Result

EIF3A expression was significantly different between ccRCC tissues and normal tissues. EIF3A expression was correlated with poor prognostic clinicopathological factors, and K–M analyses revealed that low EIF3A expression was correlated with a poor prognosis. The results of univariate and multivariate analyses proved that EIF3A was a prognostic factor in ccRCC patients. GSEA results indicated that EIF3A high expression was enriched in the renal cell carcinoma pathway. EIF3A expression was significantly positively correlated with B cells, CD8 + T cells, CD4 + T cells, neutrophils, macrophages, and dendritic cells. Furthermore, EIF3A expression was associated with most marker genes of immune cells.

Conclusions

EIF3A could serve as a potential biomarker for prognostic and diagnostic stratification of ccRCC and is related to immune cell infiltrates.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12935-021-02364-2.

Impact of alternative splicing on mechanisms of resistance to anticancer drugs

A shared characteristic of many tumors is the lack of response to anticancer drugs. Multiple mechanisms of pharmacoresistance (MPRs) are involved in permitting cancer cells to overcome the effect of these agents. Pharmacoresistance can be primary (intrinsic) or secondary (acquired), i.e., triggered or enhanced in response to the treatment. Moreover, MPRs usually result in the lack of sensitivity to several agents, which accounts for diverse multidrug-resistant (MDR) phenotypes. MPRs are based on the dynamic expression of more than one hundred genes, constituting the so-called resistome. Alternative splicing (AS) during pre-mRNA maturation results in changes affecting proteins involved in the resistome. The resulting splicing variants (SVs) reduce the efficacy of anticancer drugs by lowering the intracellular levels of active agents, altering molecular targets, enhancing both DNA repair ability and defensive mechanism of tumors, inducing changes in the balance between pro-survival and pro-apoptosis signals, modifying interactions with the tumor microenvironment, and favoring malignant phenotypic transitions. Reasons accounting for cancer-associated aberrant splicing include mutations that create or disrupt splicing sites or splicing enhancers or silencers, abnormal expression of splicing factors, and impaired signaling pathways affecting the activity of the splicing machinery. Here we have reviewed the impact of AS on MPR in cancer cells.

RNA-Sequencing-Based Transcriptomic Score with Prognostic and Theranostic Values in Multiple Myeloma

Multiple myeloma (MM) is the second most frequent hematological cancer and is characterized by the clonal proliferation of malignant plasma cells. Genome-wide expression profiling (GEP) analysis with DNA microarrays has emerged as a powerful tool for biomedical research, generating a huge amount of data. Microarray analyses have improved our understanding of MM disease and have led to important clinical applications. In MM, GEP has been used to stratify patients, define risk, identify therapeutic targets, predict treatment response, and understand drug resistance. In this study, we built a gene risk score for 267 genes using RNA-seq data that demonstrated a prognostic value in two independent cohorts (n = 674 and n = 76) of newly diagnosed MM patients treated with high-dose Melphalan and autologous stem cell transplantation. High-risk patients were associated with the expression of genes involved in several major pathways implicated in MM pathophysiology, including interferon response, cell proliferation, hypoxia, IL-6 signaling pathway, stem cell genes, MYC, and epigenetic deregulation. The RNA-seq-based risk score was correlated with specific MM somatic mutation profiles and responses to targeted treatment including EZH2, MELK, TOPK/PBK, and Aurora kinase inhibitors, outlining potential utility for precision medicine strategies in MM.

A multi-omics study delineates new molecular features and therapeutic targets for esophageal squamous cell carcinoma

Esophageal squamous cell carcinoma (ESCC) is a major histological subtype of esophageal cancer with inferior prognosis. Here, we conducted comprehensive transcriptomic, proteomic, phosphoproteomic, and metabolomic characterization of human, treatment-naive ESCC and paired normal adjacent tissues (cohort 1, n = 24) in an effort to identify new molecular vulnerabilities for ESCC and potential therapeutic targets. Integrative analysis revealed a small group of genes that were related to the active posttranscriptional and posttranslational regulation of ESCC. By using proteomic, phosphoproteomic, and metabolomic data, networks of ESCC-related signaling and metabolic pathways that were closely linked to cancer etiology were unraveled. Notably, integrative analysis of proteomic and phosphoproteomic data pinpointed that certain pathways involved in RNA transcription, processing, and metabolism were stimulated in ESCC. Importantly, proteins with close linkage to ESCC prognosis were identified. By enrolling an ESCC patient cohort 2 (n = 41), three top-ranked prognostic proteins X-prolyl aminopeptidase 3 (XPNPEP3), bromodomain PHD finger transcription factor (BPTF), and fibrillarin (FBL) were verified to have increased expression in ESCC. Among these prognostic proteins, only FBL, a well-known nucleolar methyltransferase, was essential for ESCC cell growth in vitro and in vivo. Furthermore, a validation study using an ESCC patient cohort 3 (n = 100) demonstrated that high FBL expression predicted unfavorable patient survival. Finally, common cancer/testis antigens and established cancer drivers and kinases, all of which could direct therapeutic decisions, were characterized. Collectively, our multi-omics analyses delineated new molecular features associated with ESCC pathobiology involving epigenetic, posttranscriptional, posttranslational, and metabolic characteristics, and unveiled new molecular vulnerabilities with therapeutic potential for ESCC.

Comprehensive analysis of alternative splicing profiling reveals novel events associated with prognosis and the infiltration of immune cells in prostate cancer

Background

Alternative splicing (AS) is believed to play a vital role in tumor development. Therefore, comprehensive investigation of AS and its biological function in prostate cancer (PCa) is crucial.

Methods

The AS profiling of 489 patients with PCa was obtained from The Cancer Genome Atlas (TCGA) SpliceSeq database. Bioinformatics tools were used to describe splicing associations and build prognostic models. Unsupervised clustering of the determined prognostic AS events and the relationship with immune characteristics were also explored.

Results

In total, 20,723 AS events were detected and 2,805 were identified in PCa. In the regulatory networks, the data suggested a significant correlation between splicing factor (SF) expression and AS events. To stratify the progression risk of PCa patients, prognostic models were constructed using splicing patterns. Six AS events were screened out as independent prognostic factors for progression-free survival. Based on the gene features, we constructed the combined prognostic predictors model, and the receiver operating characteristic (ROC) curve for this model reached a high area under the ROC curve (AUC) of 0.729793, indicating a favorable ability to predict patient outcomes. Through unsupervised clustering analysis, the correlations between AS-based clusters and prognosis as well as immune characteristics were revealed. The correlation analysis on TIMER revealed the relationship between gene expression and immune cell infiltration.

Conclusions

This in-depth genome-wide analysis of the AS profiling in PCa revealed unique AS events associated with cancer progression and the infiltration of immune cells, with potential for predicting outcomes and therapeutic responses.

Scutellarein induces apoptosis and inhibits proliferation, migration, and invasion in ovarian cancer via inhibition of EZH2/FOXO1 signaling

Scutellarein, a flavone found in the perennial herb Scutellaria baicalensis, has antitumorigenic activity in multiple human cancers. However, whether scutellarein can attenuate ovarian cancer (OC) is unclear. This study investigated the effects of scutellarein in OC. In vitro cell viability was assessed using MTT assay whereas proliferation was assessed using 5-ethynyl-2'-deoxyuridine and colony formation assays. Cell apoptosis was detected by an Annexin V-fluorescein isothiocyanate/propidium iodide assay. Wound-healing and Transwell assays were used to determine cell migration and invasion. The differential expression of enhancer of zeste homolog 2 (EZH2) and forkhead box protein O1 (FOXO1) was measured by Quantitative real-time PCR and western blot analysis. We found that scutellarein inhibited viability, migration, invasion of A2780 and SKOV-3 cells, and reduced the expression of EZH2 in OC cells. In addition, FOXO1 was downregulated in OC tissues and cells and negatively regulated by EZH2. Also, scutellarein inhibited tumor growth and metastasis in vivo. In conclusion, scutellarein alleviates OC by the regulation of EZH2/FOXO1 signaling.

Alternative RNA Splicing-The Trojan Horse of Cancer Cells in Chemotherapy

Almost all transcribed human genes undergo alternative RNA splicing, which increases the diversity of the coding and non-coding cellular landscape. The resultant gene products might have distinctly different and, in some cases, even opposite functions. Therefore, the abnormal regulation of alternative splicing plays a crucial role in malignant transformation, development, and progression, a fact supported by the distinct splicing profiles identified in both healthy and tumor cells. Drug resistance, resulting in treatment failure, still remains a major challenge for current cancer therapy. Furthermore, tumor cells often take advantage of aberrant RNA splicing to overcome the toxicity of the administered chemotherapeutic agents. Thus, deciphering the alternative RNA splicing variants in tumor cells would provide opportunities for designing novel therapeutics combating cancer more efficiently. In the present review, we provide a comprehensive outline of the recent findings in alternative splicing in the most common neoplasms, including lung, breast, prostate, head and neck, glioma, colon, and blood malignancies. Molecular mechanisms developed by cancer cells to promote oncogenesis as well as to evade anticancer drug treatment and the subsequent chemotherapy failure are also discussed. Taken together, these findings offer novel opportunities for future studies and the development of targeted therapy for cancer-specific splicing variants.

LINC01348 suppresses hepatocellular carcinoma metastasis through inhibition of SF3B3-mediated EZH2 pre-mRNA splicing

Long non-coding RNAs (lncRNA) play crucial roles in hepatocellular carcinoma (HCC) progression. However, the specific functions of lncRNAs in alternative splicing (AS) and the metastatic cascade in liver cancer remain largely unclear. In this study, we identified a novel lncRNA, LINC01348, which was significantly downregulated in HCC and correlated with survival functions in HCC patients. Ectopic expression of LINC01348 induced marked inhibition of cell growth, and metastasis in vitro and in vivo. Conversely, these phenotypes were reversed upon knockdown of LINC01348. Mechanistically, LINC01348 complexed with splicing factor 3b subunit 3 (SF3B3) acted as a modulator of EZH2 pre-mRNA AS, and induced alterations in JNK/c-Jun activity and expression of Snail. Notably, C-terminal truncated HBx (Ct-HBx) negatively regulated LINC01348 through c-Jun signaling. Our data collectively highlight those novel regulatory associations involving LINC01348/SF3B3/EZH2/JNK/c-Jun/Snail are an important determinant of metastasis in HCC cells and support the potential utility of targeting LINC01348 as a therapeutic strategy for HCC.

Rewards of divergence in sequences, 3-D structures and dynamics of yeast and human spliceosome SF3b complexes

The evolution of homologous and functionally equivalent multiprotein assemblies is intriguing considering sequence divergence of constituent proteins. Here, we studied the implications of protein sequence divergence on the structure, dynamics and function of homologous yeast and human SF3b spliceosomal subcomplexes. Human and yeast SF3b comprise of 7 and 6 proteins respectively, with all yeast proteins homologous to their human counterparts at moderate sequence identity. SF3b6, an additional component in the human SF3b, interacts with the N-terminal extension of SF3b1 while the yeast homologue Hsh155 lacks the equivalent region. Through detailed homology studies, we show that SF3b6 is absent not only in yeast but in multiple lineages of eukaryotes implying that it is critical in specific organisms. We probed for the potential role of SF3b6 in the spliceosome assembled form through structural and flexibility analyses. By analysing normal modes derived from anisotropic network models of SF3b1, we demonstrate that when SF3b1 is bound to SF3b6, similarities in the magnitude of residue motions (0.86) and inter-residue correlated motions (0.94) with Hsh155 are significantly higher than when SF3b1 is considered in isolation (0.21 and 0.89 respectively). We observed that SF3b6 promotes functionally relevant 'open-to-close' transition in SF3b1 by enhancing concerted residue motions. Such motions are found to occur in the Hsh155 without SF3b6. The presence of SF3b6 influences motions of 16 residues that interact with U2 snRNA/branchpoint duplex and supports the participation of its interface residues in long-range communication in the SF3b1. These results advocate that SF3b6 potentially acts as an allosteric regulator of SF3b1 for BPS selection and might play a role in alternative splicing. Furthermore, we observe variability in the relative orientation of SF3b4 and in the local structure of three β-propeller domains of SF3b3 with reference to their yeast counterparts. Such differences influence the inter-protein interactions of SF3b between these two organisms. Together, our findings highlight features of SF3b evolution and suggests that the human SF3b may have evolved sophisticated mechanisms to fine tune its molecular function.

Age-Associated Proteomic Signatures and Potential Clinically Actionable Targets of Colorectal Cancer

The occurrence and prevalence of colorectal cancer (CRC) is closely associated with age. More than 90% of patients with CRC are diagnosed after 50 years of age. However, CRC incidence of young individuals has been increasing since 1990s, whereas the overall CRC frequency is declining. Distinct overall survival rates between young and aged patients with CRC have been established. Tremendous efforts have been made to clarify the underlying mechanisms of age-dependent clinical differences, but it still remains elusive. Here, we performed proteomic profiling of 50 patients with CRC and revealed proteomic signatures of CRC across age groups. Gene set enrichment analysis showed that distinct age-dependent clinical outcomes might mainly attribute to varied MYC targets V1/V2, E2F targets and G2M checkpoint gene sets, which were associated with cancer cell proliferation, cell apoptosis, tumor growth, and tumor metastasis. Multiple linear regression analysis revealed a large number of functional proteins, such as NOP2, CSE1L, NHP2, NOC2L and CDK1, with adjusted expression significantly correlated with age (p < 0.05). Among them, NHP2 is a core component of the telomerase complex associated with age. High NHP2 expression predicted poor overall survival, with a more significant correlation in aged patients with CRC. Knockdown of NHP2 significantly suppressed cancer cell proliferation. In addition, we revealed some age-related potential clinically actionable targets, such as PSEN1, TSPO, and CDK1, which might be more suitable for patients with late-onset CRC. Collectively, this study identifies age-associated proteomic signatures and potential therapeutic targets of CRC and may help make a precise decision on CRC treatment.

•

The proteomic signatures of early-onset CRC are disclosed.

•

Alterations of some proteins between cancerous and normal tissues are age-correlated.

•

NHP2, overexpressed in tumors especially in aged patients, predicts poor prognosis.

•

Potential age-dependent druggable targets and their inhibitors are summarized.

CD46 splice variant enhances translation of specific mRNAs linked to an aggressive tumor cell phenotype in bladder cancer

CD46 is well known to be involved in diverse biological processes. Although several splice variants of CD46 have been identified, little is known about the contribution of alternative splicing to its tumorigenic functions. In this study, we found that exclusion of CD46 exon 13 is significantly increased in bladder cancer (BCa) samples. In BCa cell lines, enforced expression of CD46-CYT2 (exon 13-skipping isoform) promoted, and CD46-CYT1 (exon 13-containing isoform) attenuated, cell growth, migration, and tumorigenicity in a xenograft model. We also applied interaction proteomics to identify exhaustively the complexes containing the CYT1 or CYT2 domain in EJ-1 cells. 320 proteins were identified that interact with the CYT1 and/or CYT2 domain, and most of them are new interactors. Using an internal ribosome entry site (IRES)-dependent reporter system, we established that CD46 could regulate mRNA translation through an interaction with the translation machinery. We also identified heterogeneous nuclear ribonucleoprotein (hnRNP)A1 as a novel CYT2 binding partner, and this interaction facilitates the interaction of hnRNPA1 with IRES RNA to promote IRES-dependent translation of HIF1a and c-Myc. Strikingly, the splicing factor SRSF1 is highly correlated with CD46 exon 13 exclusion in clinical BCa samples. Taken together, our findings contribute to understanding the role of CD46 in BCa development.

Noncanonical Functions of the Polycomb Group Protein EZH2 in Breast Cancer

Enhancer of Zeste Homologue 2 (EZH2) is the catalytic subunit of the polycomb repressive complex 2 (PRC2) that is critical for determining cell identity. An epigenetic writer, EZH2 has a well-defined role in transcriptional repression by depositing trimethyl marks on lysine 27 of histone H3. However, there is mounting evidence that histone methyltransferases like EZH2 exert histone methyltransferase-independent functions. The relevance of these functions to breast cancer progression and their regulatory mechanisms are only beginning to become understood. Here, we review the current understanding of EZH2 H3K27me3-independent, noncanonical, functions and their regulation in breast cancer.

Suppression of Heterogeneous Nuclear Ribonucleoprotein C Inhibit Hepatocellular Carcinoma Proliferation, Migration, and Invasion via Ras/MAPK Signaling Pathway

Hepatocellular carcinoma (HCC), the most common malignant tumor, has high fatality and recurrence rates. Accumulating evidence shows that heterogeneous nuclear ribonucleoprotein C (HNRNPC), which is mainly involved in RNA splicing, export, and translation, promotes progression and metastasis of multiple tumor types; however, the effects of HNRNPC in HCC are unknown. In the present study, high levels of HNRNPC were detected in tumor tissues compared with para-tumor tissues by immunohistochemical and western blot assays. Furthermore, Cox proportional hazards regression models, the Kaplan–Meier method, and clinicopathologic features analysis showed that HNRNPC was not only an independent prognostic factor for both overall and disease-free survival in HCC but also a predictor of large tumor size and advanced tumor stage. Functional experiments revealed that silencing of HNRNPC not only led to arrest of more HCC cells at G0/G1 phase to inhibit their proliferation, but also suppressed EMT process to block their invasion, and migration in vitro; this was related to the Ras/MAPK signaling pathway. In addition, blocking of HCC cell proliferation regulated by HNRNPC silencing was observed in vivo. Finally, rescue tests showed that after recovery of Ras/MAPK signaling pathway activity by treatment with Ras agonists, the proliferation, migration, and invasion suppression of Huh-7 and Hep 3B cell lines caused by HNRNPC knockdown was partially reversed. Taken together, these results indicate that HNRNPC knockdown inhibits HCC cell proliferation, migration and invasion, in part via the Ras/MAPK signaling pathway. Thus, HNRNPC may have an important role in the progression of HCC and represents a promising biomarker for evaluation of prognosis and a potential therapeutic target in HCC patients.

A rare variant in EZH2 is associated with prostate cancer risk

Prostate cancer (PrCa) is highly heritable, and although rare variants contribute significantly to PrCa risk, few have been identified to date. Herein, whole-genome sequencing was performed in a large PrCa family featuring multiple affected relatives spanning several generations. A rare, predicted splice site EZH2 variant, rs78589034 (G > A), was identified as segregating with disease in all but two individuals in the family, one of whom was affected with lymphoma and bowel cancer and a female relative. This variant was significantly associated with disease risk in combined familial and sporadic PrCa datasets (n = 1551; odds ratio [OR] = 3.55, P = 1.20 × 10^-5 ). Transcriptome analysis was performed on prostate tumour needle biopsies available for two rare variant carriers and two wild-type cases. Although no allele-dependent differences were detected in EZH2 transcripts, a distinct differential gene expression signature was observed when comparing prostate tissue from the rare variant carriers with the wild-type samples. The gene expression signature comprised known downstream targets of EZH2 and included the top-ranked genes, DUSP1, FOS, JUNB and EGR1, which were subsequently validated by qPCR. These data provide evidence that rs78589034 is associated with increased PrCa risk in Tasmanian men and further, that this variant may be associated with perturbed EZH2 function in prostate tissue. Disrupted EZH2 function is a driver of tumourigenesis in several cancers, including prostate, and is of significant interest as a therapeutic target.

Novel insights into clear cell renal cell carcinoma prognosis by comprehensive characterization of aberrant alternative splicing signature: a study based on large-scale sequencing data

Clear cell renal cell carcinoma (ccRCC) is the most common type with poor prognosis in kidney tumor. Growing evidence has indicated that aberrant alternative splicing (AS) events are efficacious signatures for tumor prognosis prediction and therapeutic targets. However, the detailed roles of AS events in ccRCC are largely unknown. In our study, level 3 RNA-seq data was acquired from The Cancer Genome Atlas dataset and corresponding AS profiles were detected with the assistance of SpliceSeq software. A total of 2100 aberrant survival-associated AS events were identified via differential expression and univariate cox regression analysis. The final prognostic panel formed by 17 specific events was developed by stepwise least absolute shrinkage and selection operator (LASSO) penalty, with the area under curve (AUC) values of receiver operator characteristic (ROC) curves keeping above 0.7 spanning 1 year to 5 years. And the results from functional enrichment analyses are unanimous that autophagy could be a potential mechanism of splicing regulation in ccRCC. Furthermore, splicing regulatory network was constructed via Spearman correlation between splicing factors and AS events. Finally, unsupervised clustering analysis revealed three clusters with distinct survival patterns, and associated with specific clinicopathological phenotypes. In overall, we developed a robust and individualized predictive model based on large-scale sequencing data. The identified AS events and splicing network may be valuable in deciphering the crucial posttranscriptional mechanisms on tumorigenesis of ccRCC.

Genome-Wide Analyses of Prognostic and Therapeutic Alternative Splicing Signatures in Bladder Urothelial Carcinoma

Background

Alternative splicing (AS) is an indispensable post-transcriptional modification applied during the maturation of mRNA, and AS defects have been associated with many cancers. This study was designed to thoroughly analyze AS events in bladder urothelial carcinoma (BLCA) at the genome-wide level.

Methods

We adopted a gap analysis to screen for significant differential AS events (DASEs) associated with BLCA. DASEs with prognostic value for OS and the disease-free interval (DFI) were identified by Cox analysis. In addition, a differential AS network and AS clusters were identified using unsupervised cluster analysis. We examined differences in the sensitivity to chemotherapy and immunotherapy between BLCA patients with high and low overall survival (OS) risk.

Results

An extensive number of DASEs (296) were found to be clinically relevant in BLCA. A prognosis model was established based prognostic value of OS and DFI. CUGBP elav-like family member 2 (CELF2) was identified as a hub splicing factor for AS networks. We also identified AS clusters associated with OS using unsupervised cluster analysis, and we predicted that the effects of cisplatin and gemcitabine chemotherapy would be different between high- and low-risk groups based on OS prognosis.

Conclusion

We completed a comprehensive analysis of AS events in BLCA at the genome-wide level. The present findings revealed that DASEs and splicing factors tended to impact BLCA patient survival and sensitivity to chemotherapy drugs, which may provide novel prospects for BLCA therapies.

MicroRNA-214 enriched exosomes from human cerebral endothelial cells (hCEC) sensitize hepatocellular carcinoma to anti-cancer drugs

Hepatocellular carcinoma (HCC) is the most common primary liver tumor worldwide. Current medical therapy for HCC has limited efficacy. The present study tests the hypothesis that human cerebral endothelial cell-derived exosomes carrying elevated miR-214 (hCEC-Exo-214) can amplify the efficacy of anti-cancer drugs on HCC cells. Treatment of HepG2 and Hep3B cells with hCEC-Exo-214 in combination with anti-cancer agents, oxaliplatin or sorafenib, significantly reduced cancer cell viability and invasion compared with monotherapy with either drug. Additionally, the therapeutic effect of the combination therapy was detected in primary tumor cells derived from patients with HCC. The ability of hCEC-Exo-214 in sensitizing HCC cells to anti-cancer drugs was specific, in that combination therapy did not affect the viability and invasion of human liver epithelial cells and non-cancer primary cells. Furthermore, compared to monotherapy with oxaliplatin and sorafenib, hCEC-Exo-214 in combination with either drug substantially reduced protein levels of P-glycoprotein (P-gp) and splicing factor 3B subunit 3 (SF3B3) in HCC cells. P-gp and SF3B3 are among miR-214 target genes and are known to mediate drug resistance and cancer cell proliferation, respectively. In conclusion, the present in vitro study provides evidence that hCEC-Exo-214 significantly enhances the anti-tumor efficacy of oxaliplatin and sorafenib on HCC cells.

Prp19 Is an Independent Prognostic Marker and Promotes Neuroblastoma Metastasis by Regulating the Hippo-YAP Signaling Pathway

Pre-mRNA processing factor 19 (Prp19) was previously reported to be involved in tumor progression. However, Prp19 expression and its functions remain elusive in neuroblastoma. Here, we aim to identify the functions and mechanisms of Prp19 in neuroblastoma. Neuroblastic tumor tissue microarrays and two independent validation data sets indicate that Prp19 is associated with high-risk markers and bone marrow metastasis and serves as a prognostic marker for worse clinical outcomes with neuroblastoma. Gain- and loss-of-expression assays reveal that Prp19 promotes invasion, migration, and epithelial-mesenchymal transition (EMT) of neuroblastoma cells in vitro. Bioinformatics analysis of RNA-seq data shows that the expressions of YAP and its downstream genes are significantly inhibited after downregulation of Prp19. Prp19 and YAP expression in metastatic lymph nodes is higher than in situ neuroblastoma tissue. Further experiments show that Prp19 regulates YAP expression and consequently affects cell invasion, migration, and EMT in neuroblastoma by pre-mRNA splicing of YAP. In conclusion, our findings provide the first evidence that Prp19 is a potential therapeutic target and prognostic biomarker for patients with neuroblastoma.

Identification and Expression Pattern of EZH2 in Pig Developing Fetuses

The proper methylation status of histones is essential for appropriate cell lineage and organogenesis. EZH2, a methyltransferase catalyzing H3K27me3, has been abundantly studied in human and mouse embryonic development. The pig is an increasing important animal model for molecular study and pharmaceutical research. However, the transcript variant and temporal expression pattern of EZH2 in the middle and late porcine fetus are still unknown. Here, we identified the coding sequence of the EZH2 gene and characterized its expression pattern in fetal tissues of Duroc pigs at 65- and 90-day postcoitus (dpc). Our results showed that the coding sequence of EZH2 was 2241 bp, encoding 746 amino acids. There were 9 amino acid insertions and an amino acid substitution in this transcript compared with the validated reference sequence in NCBI. EZH2 was ubiquitously expressed in the fetal tissues of two time points with different expression levels. These results validated a different transcript in pigs and characterized its expression profile in fetal tissues of different gestation stages, which indicated that EZH2 played important roles during porcine embryonic development.

Roles of N6 -methyladenosine (m6 A) RNA modifications in urological cancers

Epigenetics has long been a hot topic in the field of scientific research. The scope of epigenetics usually includes chromatin remodelling, DNA methylation, histone modifications, non-coding RNAs and RNA modifications. In recent years, RNA modifications have emerged as important regulators in a variety of physiological processes and in disease progression, especially in human cancers. Among the various RNA modifications, m⁶ A is the most common. The function of m⁶ A modifications is mainly regulated by 3 types of proteins: m⁶ A methyltransferases (writers), m⁶ A demethylases (erasers) and m⁶ A-binding proteins (readers). In this review, we focus on RNA m⁶ A modification and its relationship with urological cancers, particularly focusing on its roles and potential clinical applications.

Identification of PDE7B as a Potential Core Gene Involved in the Metastasis of Clear Cell Renal Cell Carcinoma

Background

Metastasis is the main cause of treatment failure in various cancer, including ccRCC. However, the key genes involved in ccRCC metastasis remain largely unknown.

Purpose

The identification of the aberrant gene expression patterns associated with metastatic traits is of great clinical significance. The aim of this study was to investigate the clinical significance and function of PDE7B in ccRCC.

Materials and Methods

Expression profiling data for patient-matched primary and metastatic ccRCC tumors were obtained from GEO Dataset. Limma package was used to identify differentially expressed genes (DEGs) between the metastatic and the primary groups. Gene Ontology, Kyoto Encyclopedia of Genes Genomes (KEGG), and PPI network analysis were used to study the interacting activities and the interconnection of the DEGs. CCK-8 assays and Transwell assays were performed to detect the proliferation and migration of renal cancer cells.

Results

We obtained 163 DEGs, including 132 that were upregulated and 31 that were downregulated in metastatic ccRCC tissues. Both Gene Ontology function and KEGG pathway analysis showed that DEGs were involved in extracellular matrix (ECM) organization and cell adhesion. After utilizing PPI network to explore the interconnection among the DEGs, 22 genes were selected as the hub genes. Subsequently, survival analysis revealed that seven hub genes (SFN, NKX2-1, HP, MAPT, EPHA4, KCNAB1, and PDE7B) were significantly associated with overall survival disease-specific survival, and progression-free interval in ccRCC. Moreover, the low expression of PDE7B was found in clinical ccRCC samples and correlated with TNM stage and histologic grade. We further showed that knockdown of PDE7B increased cell growth and migration of renal cancer cells.

Conclusion

Our results implicated that PDE7B may play a key role in the development of metastatic RCC.

HnRNP A1 - mediated alternative splicing of CCDC50 contributes to cancer progression of clear cell renal cell carcinoma via ZNF395

BACKGROUND

Aberrant alternative splicing events play critical roles in carcinogenesis and progression of many cancers, while sparse studies regarding to alternative splicing are available for clear cell renal cell carcinoma (ccRCC). We identified that alternative splicing of coiled-coil domain containing 50 (CCDC50) was dysregulated in ccRCC, whereas the clinical significance of this splicing event and its splicing regulation mechanisms were still elusive.

METHODS

Bioinformatic algorithm was utilized to identify significant exon skipping events in ccRCC via exon sequencing data from The Cancer Genome Atlas. Semi-quantitative real-time polymerase chain reaction and western blot were used to validate the aberrant expression of different transcripts in renal cancer tissues, cell lines and corresponding noncancerous controls. Short hairpin RNA targeting CCDC50 and overexpressing plasmids for each transcript were introduced into ccRCC cell lines, followed by a series of in vitro and in vivo functional experiments. Moreover, a panel of splicing factors were identified and their roles on splicing regulation of CCDC50 precursor mRNA (pre-mRNA) were studied. Furthermore, RNAseq data were analyzed to elucidate downstream molecules of CCDC50. Two-way analysis of variance and unpaired Student t test were used in statistical analysis.

RESULTS

Pre-mRNA of CCDC50 generated two transcripts, full-length transcript (CCDC50-FL) and truncated transcript (CCDC50-S) with exon 6 skipped. CCDC50-S was overexpressed in ccRCC tissues and cell lines compared to noncancerous counterparts, but CCDC50-FL was only detected in noncancerous tissues and normal renal epithelial cells. Higher percent spliced-in index was associated with better survival in ccRCC patients. In vitro and in vivo functional experiments indicated that CCDC50-S transcript promoted the proliferation, migration, invasion and tumorigenesis of ccRCC, while CCDC50-FL exerted opposite tumor suppressive functions. Besides, we identified that heterogeneous nuclear ribonucleoprotein A1 (HnRNP A1) could promote the skipping of exon 6, which resulted in higher portion of CCDC50-S and oncogenic transformation. Moreover, zinc finger protein 395 (ZNF395) was identified as a downstream protein of CCDC50-S, and the interaction initiated oncogenic pathways which were involved in ccRCC progression.

CONCLUSIONS

Aberrant alternative splicing of CCDC50 is regulated by HnRNP A1 in ccRCC. This splicing event contributes to cancer progression through the downstream pathway involving ZNF395.

Decreased miR-124 contributes to the epithelial-mesenchymal transition phenotype formation of lung adenocarcinoma cells via targeting enhancer of zeste homolog 2

INTRODUCTION

MiR-124, a tumor suppressor, is involved in regulating various cellular processes. The purpose of this study was to investigate the possible function of miR-124 in LA (lung adenocarcinoma) cells.

AIMS

MiR-124 expression levels in the 54 pairs of LA tissues (and corresponding non-tumor tissues) obtained at the Sixth People's Hospital of Yancheng City and in LA cells were assessed by qRT-PCR. Colony formation assay, wound healing assay, transwell assays, attachment/detachment, western blotting and immunofluorescence assays were performed to assess the function of miR-124 on proliferation, migration and epithelial-to-mesenchymal (EMT) phenotypes in LA cells in vitro. Enhancer of zeste homolog 2 (EZH2) is identified as a target of miR-124 by bioinformatics analysis and luciferase reporter assays. Rescue assays were applied to verify the relationship between miR-124 and EZH2.

RESULTS

MiR-124 was down-regulated in LA tissues (compared to adjacent non-tumor tissues), and was down-regulated in 3 out of 4 lung cancer cell lines compared to immortalized, non-tumorigenic bronchial epithelial cells. Forced expression of miR-124 significantly suppressed tumor cell proliferation, migration and inhibited the EMT process. On the contrary, deletion of miR-124 could obviously promote cell proliferation, migration and facilitate the formation of EMT phenotype. Bioinformatics analysis and luciferase reporter assays confirmed that EZH2 was a target gene of miR-124 and was negatively correlated with the level of miR-124 in cancer tissues.

CONCLUSION

Our current study suggested that miR-124 was a tumor suppressor in LA, and miR-124 was associated with LA cell EMT phenotype formation via targeting EZH2.

Identification of survival-associated alternative splicing events and signatures in adrenocortical carcinoma based on TCGA SpliceSeq data

Objective: To explore the correlations among alternative splicing (AS), splicing factors (SF) and survival outcome in adrenocortical carcinoma (ACC) patients.

Results: A total of 92 ACC patients were included. Univariate analysis identified 3919 AS events significantly associated with overall survival. Lasso method followed by multivariate analysis revealed that the prognostic capacity of these survival-related AS events is satisfactory. Interestingly, we found that the area under the curve (AUC) of AA, AD, AP and RI were more than 0.9, indicating that these four types of AS were of great significance. Independent prognostic analysis showed that only the risk score was the independent risk factor of ACC survival. Finally, we constructed an interesting interaction network between AS and SF.

Conclusions: This is the first and most comprehensive study to explore the aberrant AS variants in ACC, which might provide novel insights into molecular mechanism of ACC.

Methods: The transcriptome data, clinical information and Percent Spliced In (PSI) values of the ACC were obtained from TCGA database and TCGA SpliceSeq data portal. Lasso method and uni/multivariate Cox regression analysis were used to identify survival-related AS events and develop multi-AS-based signatures. The relationship between AS events and SFs was also investigated.

Germ Cell-Specific Gene 1-Like Protein Regulated by Splicing Factor CUGBP Elav-Like Family Member 5 and Primary Bile Acid Biosynthesis are Prognostic in Glioblastoma Multiforme

Background

Alternative splicing (AS) modifies 92-94% human genes, abnormal splicing events might relate to tumor development and invasion. Glioblastoma Multiforme (GBM) is a fatal, invasive, and malignant tumor in nervous system. The recurrence and development leads to poor prognosis. However, few studies have focused on AS in GBM.

Methods

RNA-seq and Alternative splicing events (ASEs) data of GBM samples were downloaded from The Cancer Genome Atlas (TCGA) and TCGASpliceSeq databases, respectively. Firstly, the Cox regression analysis was utilized to identify the overall survival splicing events (OS-SEs). Secondly, a multivariable model was applied to access the prognostic value of risk score. Then, we constructed a co-expressed network between splicing factors (SFs) and overall survival alternative splicing events (OS-SEs). Additionally, to explore the relationship between the potential prognostic signaling pathways and OS-SEs, we constructed a network between these pathways and OS-SEs. Ultimately, to better explain the results, validations from multi-dimension platforms were applied.

Results

In the first step, 1,062 OS-SEs were selected by Cox regression. Then, 11 OS-SEs were integrated in a multivariate model by Lasso regression. The area under the curve (AUC) of receiver operator characteristic (ROC) curve was 0.861. In addition, the risk score generated from the multivariate model was confirmed to be an independent prognostic factor (P < 0.001). What's more, in the network of SFs and ASEs, CELF5 significantly regulated GSG1L|35696|AP and GSG1L|35698|AP (P < 0.001, R = 0.511 and = -0.492). Additionally, GSG1L|35696|AP (P = 0.006) and GSG1L|35698|AP (P = 0.007) showed a significant relationship with cancer status. Eventually, KEGG pathways related to prognosis of GBM were selected by GSVA. The primary bile acid synthesis (P < 0.001, R = 0.420) was the significant pathway co-expressed with Germ Cell-Specific Gene 1-Like Protein (GSG1L).

Conclusions

Based on the comprehensive bioinformatics analysis, we proposed that aberrant splicing factor CUGBP Elav-like family member 5 (CELF5) significantly, positively and negatively, regulated ASE of GSG1L, and the primary bile acid synthesis pathway might play an important role in tumorigenesis and prognosis of GBM.