An integrative analysis of colon cancer identifies an essential function for PRPF6 in tumor growth

PAIP1 binds to pre-mRNA and regulates alternative splicing of cancer pathway genes including VEGFA

BACKGROUND

Poly (A) binding protein interacting protein 1 (PAIP1) has been shown to causally contribute to the development and progression of cancer. However, the mechanisms of the PAIP1 regulation in tumor cells remain poorly understood.

RESULTS

Here, we used a recently developed UV cross-linking and RNA immunoprecipitation method (iRIP-seq) to map the direct and indirect interaction sites between PAIP1 and RNA on a transcriptome-wide level in HeLa cells. We found that PAIP1 not only binds to 3'UTRs, but also to pre-mRNAs/mRNAs with a strong bias towards the coding region and intron. PAIP1 binding sites are enriched in splicing enhancer consensus GA-rich motifs. RNA-seq analysis revealed that PAIP1 selectively modulates the alternative splicing of genes in some cancer hallmarks including cell migration, the mTOR signaling pathway and the HIF-1 signaling pathway. PAIP1-regulated alternative splicing events were strongly associated with PAIP1 binding, demonstrating that the binding may promote selection of the nearby splice sites. Deletion of a PAIP1 binding site containing seven repeats of GA motif reduced the PAIP1-mediated suppression of the exon 6 inclusion in a VEGFA mRNA isoform. Proteomic analysis of the PAIP1-interacted proteins revealed the enrichment of the spliceosome components and splicing factors.

CONCLUSIONS

These findings suggest that PAIP1 is both a polyadenylation and alternative splicing regulator, that may play a large role in RNA processing via its role in alternative splicing regulation.

Spliceosomic dysregulation in pancreatic cancer uncovers splicing factors PRPF8 and RBMX as novel candidate actionable targets

Pancreatic ductal adenocarcinoma (PDAC) is a highly lethal cancer, characterized by late diagnosis and poor treatment response. Surgery is the only curative approach, only available to early-diagnosed patients. Current therapies have limited effects, cause severe toxicities, and minimally improve overall survival. Understanding of splicing machinery alterations in PDAC remains incomplete. Here, we comprehensively examined 59 splicing machinery components, uncovering dysregulation in pre-mRNA processing factor 8 (PRPF8) and RNA-binding motif protein X-linked (RBMX). Their downregulated expression was linked to poor prognosis and malignancy features, including tumor stage, invasion and metastasis, and associated with poorer survival and the mutation of key PDAC genes. Experimental modulation of these splicing factors in pancreatic cancer cell lines reverted their expression to non-tumor levels and resulted in decreased key tumor-related features. These results provide evidence that the splicing machinery is altered in PDAC, wherein PRPF8 and RBMX emerge as candidate actionable therapeutic targets.

Dysregulated ribosome quality control in human diseases

Precise regulation of mRNA translation is of fundamental importance for maintaining homeostasis. Conversely, dysregulated general or transcript-specific translation, as well as abnormal translation events, have been linked to a multitude of diseases. However, driven by the misconception that the transient nature of mRNAs renders their abnormalities inconsequential, the importance of mechanisms that monitor the quality and fidelity of the translation process has been largely overlooked. In recent years, there has been a dramatic shift in this paradigm, evidenced by several seminal discoveries on the role of a key mechanism in monitoring the quality of mRNA translation - namely, Ribosome Quality Control (RQC) - in the maintenance of homeostasis and the prevention of diseases. Here, we will review recent advances in the field and emphasize the biological significance of the RQC mechanism, particularly its implications in human diseases.

Identification of novel molecular subtypes to improve the classification framework of nasopharyngeal carcinoma

BACKGROUND

Nasopharyngeal carcinoma (NPC) treatment is largely based on a 'one-drug-fits-all' strategy in patients with similar pathological characteristics. However, given its biological heterogeneity, patients at the same clinical stage or similar therapies exhibit significant clinical differences. Thus, novel molecular subgroups based on these characteristics may better therapeutic outcomes.

METHODS

Herein, 192 treatment-naïve NPC samples with corresponding clinicopathological information were obtained from Fujian Cancer Hospital between January 2015 and January 2018. The gene expression profiles of the samples were obtained by RNA sequencing. Molecular subtypes were identified by consensus clustering. External NPC cohorts were used as the validation sets.

RESULTS

Patients with NPC were classified into immune, metabolic, and proliferative molecular subtypes with distinct clinical features. Additionally, this classification was repeatable and predictable as validated by the external NPC cohorts. Metabolomics has shown that arachidonic acid metabolites were associated with NPC malignancy. We also identified several key genes in each subtype using a weighted correlation network analysis. Furthermore, a prognostic risk model based on these key genes was developed and was significantly associated with disease-free survival (hazard ratio, 1.11; 95% CI, 1.07-1.16; P < 0.0001), which was further validated by an external NPC cohort (hazard ratio, 7.71; 95% CI, 1.39-42.73; P < 0.0001). Moreover, the 1-, 3-, and 5-year areas under the curve were 0.84 (95% CI, 0.74-0.94), 0.81 (95% CI, 0.73-0.89), and 0.82 (95% CI, 0.73-0.90), respectively, demonstrating a high predictive value.

CONCLUSIONS

Overall, we defined a novel classification of nasopharyngeal carcinoma (immune, metabolism, and proliferation subtypes). Among these subtypes, metabolism and proliferation subtypes were associated with advanced stage and poor prognosis of NPC patients, whereas the immune subtype was linked to early stage and favorable prognosis.

Xrp1 governs the stress response program to spliceosome dysfunction

Co-transcriptional processing of nascent pre-mRNAs by the spliceosome is vital to regulating gene expression and maintaining genome integrity. Here, we show that the deficiency of functional U5 small nuclear ribonucleoprotein particles (snRNPs) in Drosophila imaginal cells causes extensive transcriptome remodeling and accumulation of highly mutagenic R-loops, triggering a robust stress response and cell cycle arrest. Despite compromised proliferative capacity, the U5 snRNP-deficient cells increased protein translation and cell size, causing intra-organ growth disbalance before being gradually eliminated via apoptosis. We identify the Xrp1-Irbp18 heterodimer as the primary driver of transcriptional and cellular stress program downstream of U5 snRNP malfunction. Knockdown of Xrp1 or Irbp18 in U5 snRNP-deficient cells attenuated JNK and p53 activity, restored normal cell cycle progression and growth, and inhibited cell death. Reducing Xrp1-Irbp18, however, did not rescue the splicing defects, highlighting the requirement of accurate splicing for cellular and tissue homeostasis. Our work provides novel insights into the crosstalk between splicing and the DNA damage response and defines the Xrp1-Irbp18 heterodimer as a critical sensor of spliceosome malfunction and mediator of the stress-induced cellular senescence program.

Integrated bioinformatics approach to unwind key genes and pathways involved in colorectal cancer

BACKGROUND

Colorectal cancer (CRC) is the fifth leading cause of death in India. Until now, the exact pathogenesis concerning CRC signaling pathways is largely unknown; however, the diseased condition is believed to deteriorate with lifestyle, aging, and inherited genetic disorders. Hence, the identification of hub genes and therapeutic targets is of great importance for disease monitoring.

OBJECTIVE

Identification of hub genes and targets for identification of candidate hub genes for CRC diagnosis and monitoring.

MATERIALS AND METHODS

The present study applied gene expression analysis by integrating two profile datasets (GSE20916 and GSE33113) from NCBI-GEO database to elucidate the potential key candidate genes and pathways in CRC. Differentially expressed genes (DEGs) between CRC (195 CRC tissues) and healthy control (46 normal mucosal tissue) were sorted using GEO2R tool. Further, Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analysis were performed using Cluster Profiler in Rv. 3.6.1. Moreover, protein-protein interactions (PPI), module detection, and hub gene identification were accomplished and visualized through the Search Tool for the Retrieval of Interacting Genes, Molecular Complex Detection (MCODE) plug-in of Cytoscape v3.8.0. Further hub genes were imported into ToppGene webserver for pathway analysis and prognostic expression analysis was conducted using Gene Expression Profiling Interactive Analysis webserver.

RESULTS

A total of 2221 DEGs, including 1286 up-regulated and 935down-regulated genes mainly enriched in signaling pathways of NOD-like receptor, FoxO, AMPK signalling and leishmaniasis. Three key modules were detected from PPI network using MCODE. Besides, top 20 high prioritized hub genes were selected. Further, prognostic expression analysis revealed ten of the hub genes, namely IL1B, CD44, Glyceraldehyde-3-phosphate dehydrogenase (GAPDH, MMP9, CREB1, STAT1, vascular endothelial growth factor (VEGFA), CDC5 L, Ataxia-telangiectasia mutated (ATM + and CDH1 to be differently expressed in normal and cancer patients.

CONCLUSION

The present study proposed five novel therapeutic targets, i.e., ATM, GAPDH, CREB1, VEGFA, and CDH1 genes that might provide new insights into molecular oncogenesis of CRC.

Data mining identifies novel RNA-binding proteins involved in colon and rectal carcinomas

Colorectal adenocarcinoma (COREAD) is the second most deadly cancer and third most frequently encountered malignancy worldwide. Despite efforts in molecular subtyping and subsequent personalized COREAD treatments, multidisciplinary evidence suggests separating COREAD into colon cancer (COAD) and rectal cancer (READ). This new perspective could improve diagnosis and treatment of both carcinomas. RNA-binding proteins (RBPs), as critical regulators of every hallmark of cancer, could fulfill the need to identify sensitive biomarkers for COAD and READ separately. To detect new RBPs involved in COAD and READ progression, here we used a multidata integration strategy to prioritize tumorigenic RBPs. We analyzed and integrated 1) RBPs genomic and transcriptomic alterations from 488 COAD and 155 READ patients, 2) ∼ 10,000 raw associations between RBPs and cancer genes, 3) ∼ 15,000 immunostainings, and 4) loss-of-function screens performed in 102 COREAD cell lines. Thus, we unraveled new putative roles of NOP56, RBM12, NAT10, FKBP1A, EMG1, and CSE1L in COAD and READ progression. Interestingly, FKBP1A and EMG1 have never been related with any of these carcinomas but presented tumorigenic features in other cancer types. Subsequent survival analyses highlighted the clinical relevance of FKBP1A, NOP56, and NAT10 mRNA expression to predict poor prognosis in COREAD and COAD patients. Further research should be performed to validate their clinical potential and to elucidate their molecular mechanisms underlying these malignancies.

Common genetic variations in telomere length genes and lung cancer: a Mendelian randomisation study and its novel application in lung tumour transcriptome

Background

Genome-wide association studies (GWASs) have identified genetic susceptibility variants for both leukocyte telomere length (LTL) and lung cancer susceptibility. Our study aims to explore the shared genetic basis between these traits and investigate their impact on somatic environment of lung tumours.

Methods

We performed genetic correlation, Mendelian randomisation (MR), and colocalisation analyses using the largest available GWASs summary statistics of LTL (N=464,716) and lung cancer (N=29,239 cases and 56,450 controls). Principal components analysis based on RNA-sequencing data was used to summarise gene expression profile in lung adenocarcinoma cases from TCGA (N=343).

Results

Although there was no genome-wide genetic correlation between LTL and lung cancer risk, longer LTL conferred an increased risk of lung cancer regardless of smoking status in the MR analyses, particularly for lung adenocarcinoma. Of the 144 LTL genetic instruments, 12 colocalised with lung adenocarcinoma risk and revealed novel susceptibility loci, including MPHOSPH6, PRPF6, and POLI. The polygenic risk score for LTL was associated with a specific gene expression profile (PC2) in lung adenocarcinoma tumours. The aspect of PC2 associated with longer LTL was also associated with being female, never smokers, and earlier tumour stages. PC2 was strongly associated with cell proliferation score and genomic features related to genome stability, including copy number changes and telomerase activity.

Conclusions

This study identified an association between longer genetically predicted LTL and lung cancer and sheds light on the potential molecular mechanisms related to LTL in lung adenocarcinomas.

Funding

Institut National du Cancer (GeniLuc2017-1-TABAC-03-CIRC-1-TABAC17-022), INTEGRAL/NIH (5U19CA203654-03), CRUK (C18281/A29019), and Agence Nationale pour la Recherche (ANR-10-INBS-09).

Alternative mRNA Splicing and Promising Therapies in Cancer

Cancer is among the leading causes of mortality worldwide. While considerable attention has been given to genetic and epigenetic sources of cancer-specific cellular activities, the role of alternative mRNA splicing has only recently received attention as a major contributor to cancer initiation and progression. The distribution of alternate mRNA splicing variants in cancer cells is different from their non-cancer counterparts, and cancer cells are more sensitive than non-cancer cells to drugs that target components of the splicing regulatory network. While many of the alternatively spliced mRNAs in cancer cells may represent "noise" from splicing dysregulation, certain recurring splicing variants have been shown to contribute to tumor progression. Some pathogenic splicing disruption events result from mutations in cis-acting splicing regulatory sequences in disease-associated genes, while others may result from shifts in balance among naturally occurring alternate splicing variants among mRNAs that participate in cell cycle progression and the regulation of apoptosis. This review provides examples of cancer-related alternate splicing events resulting from each step of mRNA processing and the promising therapies that may be used to address them.

Microsatellite Instability and Aberrant Pre-mRNA Splicing: How Intimate Is It?

Cancers that belong to the microsatellite instability (MSI) class can account for up to 15% of all cancers of the digestive tract. These cancers are characterized by inactivation, through the mutation or epigenetic silencing of one or several genes from the DNA MisMatch Repair (MMR) machinery, including MLH1, MLH3, MSH2, MSH3, MSH6, PMS1, PMS2 and Exo1. The unrepaired DNA replication errors turn into mutations at several thousand sites that contain repetitive sequences, mainly mono- or dinucleotides, and some of them are related to Lynch syndrome, a predisposition condition linked to a germline mutation in one of these genes. In addition, some mutations shortening the microsatellite (MS) stretch could occur in the 3'-intronic regions, i.e., in the ATM (ATM serine/threonine kinase), MRE11 (MRE11 homolog) or the HSP110 (Heat shock protein family H) genes. In these three cases, aberrant pre-mRNA splicing was observed, and it was characterized by the occurrence of selective exon skipping in mature mRNAs. Because both the ATM and MRE11 genes, which as act as players in the MNR (MRE11/NBS1 (Nibrin)/RAD50 (RAD50 double strand break repair protein) DNA damage repair system, participate in double strand breaks (DSB) repair, their frequent splicing alterations in MSI cancers lead to impaired activity. This reveals the existence of a functional link between the MMR/DSB repair systems and the pre-mRNA splicing machinery, the diverted function of which is the consequence of mutations in the MS sequences.

Aberrant Retinal Pigment Epithelial Cells Derived from Induced Pluripotent Stem Cells of a Retinitis Pigmentosa Patient with the PRPF6 Mutation

Pre-mRNA processing factors (PRPFs) are vital components of the spliceosome and are involved in the physiological process necessary for pre-mRNA splicing to mature mRNA. As an important member, PRPF6 mutation resulting in autosomal dominant retinitis pigmentosa (adRP) is not common. Recently, we reported the establishment of an induced pluripotent stem cells (iPSCs; CSUASOi004-A) model by reprogramming the peripheral blood mononuclear cells of a PRPF6-related adRP patient, which could recapitulate a consistent disease-specific genotype. In this study, a disease model of retinal pigment epithelial (RPE) cells was generated from the iPSCs of this patient to further investigate the underlying molecular and pathological mechanisms. The results showed the irregular morphology, disorganized apical microvilli and reduced expressions of RPE-specific genes in the patient’s iPSC-derived RPE cells. In addition, RPE cells carrying the PRPF6 mutation displayed a decrease in the phagocytosis of fluorescein isothiocyanate-labeled photoreceptor outer segments and exhibited impaired cell polarity and barrier function. This study will benefit the understanding of PRPF6-related RPE cells and future cell therapy.

Hereditary retinoblastoma iPSC model reveals aberrant spliceosome function driving bone malignancies

Rare human hereditary disorders provide unequivocal evidence of the role of gene mutations in human disease pathogenesis and offer powerful insights into their influence on human disease development. Using a hereditary retinoblastoma (RB) patient–derived induced pluripotent stem cell (iPSC) platform, we elucidate the role of pRB/E2F3a in regulating spliceosomal gene expression. Pharmacological inhibition of the spliceosome in RB1-mutant cells preferentially increases splicing abnormalities of genes involved in cancer-promoting signaling and impairs cell proliferation and tumorigenesis. Expression of pRB/E2F3a–regulated spliceosomal proteins is negatively associated with pRB expression and correlates with poor clinical outcomes of osteosarcoma (OS) patients. Our findings strongly indicate that the spliceosome is an “Achilles’ heel” of RB1-mutant OS.

The RB1 gene is frequently mutated in human cancers but its role in tumorigenesis remains incompletely defined. Using an induced pluripotent stem cell (iPSC) model of hereditary retinoblastoma (RB), we report that the spliceosome is an up-regulated target responding to oncogenic stress in RB1-mutant cells. By investigating transcriptomes and genome occupancies in RB iPSC–derived osteoblasts (OBs), we discover that both E2F3a, which mediates spliceosomal gene expression, and pRB, which antagonizes E2F3a, coregulate more than one-third of spliceosomal genes by cobinding to their promoters or enhancers. Pharmacological inhibition of the spliceosome in RB1-mutant cells leads to global intron retention, decreased cell proliferation, and impaired tumorigenesis. Tumor specimen studies and genome-wide TCGA (The Cancer Genome Atlas) expression profile analyses support the clinical relevance of pRB and E2F3a in modulating spliceosomal gene expression in multiple cancer types including osteosarcoma (OS). High levels of pRB/E2F3a–regulated spliceosomal genes are associated with poor OS patient survival. Collectively, these findings reveal an undiscovered connection between pRB, E2F3a, the spliceosome, and tumorigenesis, pointing to the spliceosomal machinery as a potentially widespread therapeutic vulnerability of pRB-deficient cancers.

Probing Isoform Switching Events in Various Cancer Types: Lessons From Pan-Cancer Studies

Alternative splicing is an essential regulatory mechanism for gene expression in mammalian cells contributing to protein, cellular, and species diversity. In cancer, alternative splicing is frequently disturbed, leading to changes in the expression of alternatively spliced protein isoforms. Advances in sequencing technologies and analysis methods led to new insights into the extent and functional impact of disturbed alternative splicing events. In this review, we give a brief overview of the molecular mechanisms driving alternative splicing, highlight the function of alternative splicing in healthy tissues and describe how alternative splicing is disrupted in cancer. We summarize current available computational tools for analyzing differential transcript usage, isoform switching events, and the pathogenic impact of cancer-specific splicing events. Finally, the strategies of three recent pan-cancer studies on isoform switching events are compared. Their methodological similarities and discrepancies are highlighted and lessons learned from the comparison are listed. We hope that our assessment will lead to new and more robust methods for cancer-specific transcript detection and help to produce more accurate functional impact predictions of isoform switching events.

Alternative splicing of mRNA in colorectal cancer: new strategies for tumor diagnosis and treatment

Alternative splicing (AS) is an important event that contributes to posttranscriptional gene regulation. This process leads to several mature transcript variants with diverse physiological functions. Indeed, disruption of various aspects of this multistep process, such as cis- or trans- factor alteration, promotes the progression of colorectal cancer. Therefore, targeting some specific processes of AS may be an effective therapeutic strategy for treating cancer. Here, we provide an overview of the AS events related to colorectal cancer based on research done in the past 5 years. We focus on the mechanisms and functions of variant products of AS that are relevant to malignant hallmarks, with an emphasis on variants with clinical significance. In addition, novel strategies for exploiting the therapeutic value of AS events are discussed.

Elucidation of the BMI1 interactome identifies novel regulatory roles in glioblastoma

Glioblastoma (GBM) is the most common and aggressive intrinsic brain tumour in adults. Epigenetic mechanisms controlling normal brain development are often dysregulated in GBM. Among these, BMI1, a structural component of the Polycomb Repressive Complex 1 (PRC1), which promotes the H2AK119ub catalytic activity of Ring1B, is upregulated in GBM and its tumorigenic role has been shown in vitro and in vivo. Here, we have used protein and chromatin immunoprecipitation followed by mass spectrometry (MS) analysis to elucidate the protein composition of PRC1 in GBM and transcriptional silencing of defining interactors in primary patient-derived GIC lines to assess their functional impact on GBM biology. We identify novel regulatory functions in mRNA splicing and cholesterol transport which could represent novel targetable mechanisms in GBM.

Emerging roles of spliceosome in cancer and immunity

Precursor messenger RNA (pre-mRNA) splicing is catalyzed by an intricate ribonucleoprotein complex called the spliceosome. Although the spliceosome is considered to be general cell “housekeeping” machinery, mutations in core components of the spliceosome frequently correlate with cell- or tissue-specific phenotypes and diseases. In this review, we expound the links between spliceosome mutations, aberrant splicing, and human cancers. Remarkably, spliceosome-targeted therapies (STTs) have become efficient anti-cancer strategies for cancer patients with splicing defects. We also highlight the links between spliceosome and immune signaling. Recent studies have shown that some spliceosome gene mutations can result in immune dysregulation and notable phenotypes due to mis-splicing of immune-related genes. Furthermore, several core spliceosome components harbor splicing-independent immune functions within the cell, expanding the functional repertoire of these diverse proteins.

CD44 alternative splicing senses intragenic DNA methylation in tumors via direct and indirect mechanisms

DNA methylation (meDNA) is a modulator of alternative splicing, and splicing perturbations are involved in tumorigenesis nearly as frequently as DNA mutations. However, the impact of meDNA on tumorigenesis via splicing-mediated mechanisms has not been thoroughly explored. Here, we found that HCT116 colon carcinoma cells inactivated for the DNA methylases DNMT1/3b undergo a partial epithelial to mesenchymal transition associated with increased CD44 variant exon skipping. These skipping events are directly mediated by the loss of intragenic meDNA and the chromatin factors MBD1/2/3 and HP1γ and are also linked to phosphorylation changes in elongating RNA polymerase II. The role of meDNA in alternative splicing was confirmed by using the dCas9/DNMT3b tool. We further tested whether the meDNA level could have predictive value in the MCF10A model for breast cancer progression and in patients with acute lymphoblastic leukemia (B ALL). We found that a small number of differentially spliced genes, mostly involved in splicing and signal transduction, are correlated with the local modulation of meDNA. Our observations suggest that, although DNA methylation has multiple avenues to affect alternative splicing, its indirect effect may also be mediated through alternative splicing isoforms of these meDNA sensors.

Epigenetic Mechanisms Are Involved in the Oncogenic Properties of ZNF518B in Colorectal Cancer

Simple Summary

The ZNF518B gene, which is up-regulated in colorectal cancer, plays a role in metastasis, but neither the mechanisms involved in this process nor the role of the different isoforms of the gene are known. Here we show that the ratio of these isoforms is related to the relapsing of the disease, and that the protein ZNF518B interacts with enzymes able to introduce epigenetic changes, which may affect the activity of many genes. We also report a list of genes affected in common by ZNF518B and by two of those related enzymes, namely, G9A and EZH2. An in-depth analysis of five of those genes revealed that ZNF518B is involved in the recruitment of the enzymes and in the deposition of the corresponding epigenetic marks. The results highlight the relevance of epigenetic changes in cancer development, and open the possibility of developing therapeutic approaches, as the introduction of epigenetic modifications is reversible.

Abstract

The ZNF518B gene, which is up-regulated in colorectal cancer, plays a role in cell dissemination and metastasis. It encodes a zinc-finger protein, which interacts with histone methyltransferases G9A and EZH2. The expression of the two major mRNA isoforms 1 (coding for the full protein) and 2 was quantified by RT-qPCR in a cohort of 66 patients. The effects of silencing ZNF518B on the transcriptome of DLD1 and HCT116 cells were analysed by Clariom-S assays and validated by RT-qPCR. The recruitment of methyltransferases and the presence of H3K27me3 were studied by chromatin immunoprecipitation (ChIP). The ratio (isoform 2)/(isoform 1) negatively correlated with the relapsing of disease. The study of the transcriptome of DLD1 and HCT116 cells revealed that many genes affected by silencing ZNF518B are related to cancer. After crossing these results with the list of genes affected by silencing the histone methyltransferases (retrieved in silico), five genes were selected. ChIP analysis revealed that the recruitment of EZH2 is ZNF518B-dependent in KAT2B, RGS4 and EFNA5; the level of H3K27me3 changes in accordance. G9A also binds RGS4 and PADI3 in a ZNF518B-dependent manner. The results highlight the importance of epigenetics in cancer and open a novel therapeutic possibility, as inhibition of histone methyltransferases may reverse the disease-linked histone marks.

The Role of the U5 snRNP in Genetic Disorders and Cancer

Pre-mRNA splicing is performed by the spliceosome, a dynamic macromolecular complex consisting of five small uridine-rich ribonucleoprotein complexes (the U1, U2, U4, U5, and U6 snRNPs) and numerous auxiliary splicing factors. A plethora of human disorders are caused by genetic variants affecting the function and/or expression of splicing factors, including the core snRNP proteins. Variants in the genes encoding proteins of the U5 snRNP cause two distinct and tissue-specific human disease phenotypes – variants in PRPF6, PRPF8, and SNRP200 are associated with retinitis pigmentosa (RP), while variants in EFTUD2 and TXNL4A cause the craniofacial disorders mandibulofacial dysostosis Guion-Almeida type (MFDGA) and Burn-McKeown syndrome (BMKS), respectively. Furthermore, recurrent somatic mutations or changes in the expression levels of a number of U5 snRNP proteins (PRPF6, PRPF8, EFTUD2, DDX23, and SNRNP40) have been associated with human cancers. How and why variants in ubiquitously expressed spliceosome proteins required for pre-mRNA splicing in all human cells result in tissue-restricted disease phenotypes is not clear. Additionally, why variants in different, yet interacting, proteins making up the same core spliceosome snRNP result in completely distinct disease outcomes – RP, craniofacial defects or cancer – is unclear. In this review, we define the roles of different U5 snRNP proteins in RP, craniofacial disorders and cancer, including how disease-associated genetic variants affect pre-mRNA splicing and the proposed disease mechanisms. We then propose potential hypotheses for how U5 snRNP variants cause tissue specificity resulting in the restricted and distinct human disorders.

Construction of an RNA-Binding Protein-Related Prognostic Model for Pancreatic Adenocarcinoma Based on TCGA and GTEx Databases

Background: Recently, RNA-binding proteins (RBPs) were reported to interact with target mRNA to regulate gene posttranscriptional expression, and RBP-mediated RNA modification can regulate the expression and function of proto-oncogenes and tumor suppressor genes. We systematically analyzed the expression of RBPs in pancreatic adenocarcinoma (PAAD) and constructed an RBP-associated prognostic risk model. Methods: Gene expression data of normal pancreatic samples as well as PAAD samples were downloaded from TCGA-PAAD and GTEx databases. Wilcoxon test and univariate Cox analysis were, respectively, applied to screen differential expression RBPs (DE-RBPs) and prognostic-associated RBPs (pRBPs). Functional enrichment was analyzed by GO, KEGG, and GSEA. Protein-protein interaction (PPI) network was constructed by STRING online database. Modeling RBPs were selected by multivariate Cox analysis. Kaplan-Meier survival and Cox analysis were applied to evaluate the effects of risk score on the overall survival of PAAD patients. ROC curves and validation cohort were applied to verify the accuracy of the model. Nomogram was applied for predicting 1-, 3-, and 5-year overall survival (OS) of PAAD patients. At last, modeling RBPs were further analyzed to explore their differential expression, prognostic value, as well as enrichment pathways in PAAD. Results: RBPs (453) were differentially expressed in normal and tumor samples, besides, 28 of which were prognostic associated. DE-RBPs (453) are functionally associated with ribosome, ribonuclease, spliceosome, etc. Eight RBPs (PABPC1, PRPF6, OAS1, RBM5, LSM12, IPO7, FXR1, and RBM6) were identified to construct a prognostic risk model. Higher risk score not only predicted poor prognosis but also was an independent poor prognostic indicator, which was verified by ROC curves and validation cohort. Eight modeling RBPs were confirmed to be significantly differentially expressed between normal and tumor samples from RNA and protein level. Besides, all of eight RBPs were related with overall survival of PAAD patients. Conclusions: We successfully constructed an RBP-associated prognostic risk model in PAAD, which has a potential clinical application prospect.

RNA in cancer

While the processing of mRNA is essential for gene expression, recent findings have highlighted that RNA processing is systematically altered in cancer. Mutations in RNA splicing factor genes and the shortening of 3' untranslated regions are widely observed. Moreover, evidence is accumulating that other types of RNAs, including circular RNAs, can contribute to tumorigenesis. In this Review, we highlight how altered processing or activity of coding and non-coding RNAs contributes to cancer. We introduce the regulation of gene expression by coding and non-coding RNA and discuss both established roles (microRNAs and long non-coding RNAs) and emerging roles (selective mRNA processing and circular RNAs) for RNAs, highlighting the potential mechanisms by which these RNA subtypes contribute to cancer. The widespread alteration of coding and non-coding RNA demonstrates that altered RNA biogenesis contributes to multiple hallmarks of cancer.

PRPF6 promotes androgen receptor/androgen receptor-variant 7 actions in castration-resistant prostate cancer cells

Androgen receptor (AR) and its variants play vital roles in development and progression of prostate cancer. To clarify the mechanisms involved in the enhancement of their actions would be crucial for understanding the process in prostate cancer and castration-resistant prostate cancer transformation. Here, we provided the evidence to show that pre-mRNA processing factor 6 (PRPF6) acts as a key regulator for action of both AR full length (AR-FL) and AR variant 7 (AR-V7), thereby participating in the enhancement of AR-FL and AR-V7-induced transactivation in prostate cancer. In addition, PRPF6 is recruited to cis-regulatory elements in AR target genes and associates with JMJD1A to enhance AR-induced transactivation. PRPF6 also promotes expression of AR-FL and AR-V7. Moreover, PRPF6 depletion reduces tumor growth in prostate cancer-derived cell lines and results in significant suppression of xenograft tumors even under castration condition in mouse model. Furthermore, PRPF6 is obviously highly expressed in human prostate cancer samples. Collectively, our results suggest PRPF6 is involved in enhancement of oncogenic AR signaling, which support a previously unknown role of PRPF6 during progression of prostate cancer and castration-resistant prostate cancers.

Splicing factor PRPF6 upregulates oncogenic androgen receptor signaling pathway in hepatocellular carcinoma

Androgen receptor (AR) signaling is considered to be crucial for the pathogenesis of hepatocellular carcinoma (HCC) with obvious sexual dimorphism. Pre‐mRNA processing factor 6 (PRPF6) was identified as a coactivator of AR. However, the molecular mechanism underlying the modulation function of PRPF6 on AR‐mediated transcriptional activity in HCC needs to be further clarified. In this study, we analyzed data from The Cancer Genome Atlas to show that PRPF6 is highly expressed in HCC. . Our data indicated that PRPF6 interacts with AR/AR splice variants (AR‐Vs) and upregulates AR/AR splice variant 7‐mediated transcriptional activity even without dihydrotestosterone treatment. We observed that AR is obviously induced by androgen treatment and is mainly expressed in the nucleus in HCC‐derived cell lines. Moreover, overexpression of PRPF6 enhances AR expression accompanied with the increase of AR‐Vs expression. We provided evidence that PRPF6 participates in upregulating AR self‐transcription. PRPF6 facilitates the recruitment of AR to the androgen responsive element region of the AR gene. Finally, PRPF6 depletion inhibits cell proliferation in HCC cells and mouse xenografts. Taken together, our results suggest that PRPF6 as a splicing factor enhances AR self‐transcription, thereby coactivating oncogenic AR/AR‐Vs actions in HCC.

PCDH7 interacts with GluN1 and regulates dendritic spine morphology and synaptic function

The N-terminal domain (NTD) of the GluN1 subunit (GluN1-NTD) is important for NMDA receptor structure and function, but the interacting proteins of the GluN1-NTD are not well understood. Starting with an unbiased screen of ~ 1,500 transmembrane proteins using the purified GluN1-NTD protein as a bait, we identify Protocadherin 7 (PCDH7) as a potential interacting protein. PCDH7 is highly expressed in the brain and has been linked to CNS disorders, including epilepsy. Using primary neurons and brain slice cultures, we find that overexpression and knockdown of PCDH7 induce opposing morphological changes of dendritic structures. We also find that PCDH7 overexpression reduces synaptic NMDA receptor currents. These data show that PCDH7 can regulate dendritic spine morphology and synaptic function, possibly via interaction with the GluN1 subunit.

Comparative proteogenomics profiling of non-small and small lung carcinoma cell lines using mass spectrometry

Background

Evidences indicated that non-small-cell lung cancer (NSCLC) and small-cell lung cancer (SCLC) might originate from the same cell type, which however ended up to be two different subtypes of lung carcinoma, requiring different therapeutic regimens. We aimed to identify the differences between these two subtypes of lung cancer by using integrated proteome and genome approaches.

Methods and Materials

Two representative cell lines for each lung cancer subtype were comparatively analysed by quantitative proteomics, and their corresponding transcriptomics data were obtained from the Gene Expression Omnibus database. The integrated analyses of proteogenomic data were performed to determine key differentially expressed proteins that were positively correlated between proteomic and transcriptomic data.

Result

The proteomics analysis revealed 147 differentially expressed proteins between SCLC and NSCLC from a total of 3,970 identified proteins. Combined with available transcriptomics data, we further confirmed 14 differentially expressed proteins including six known and eight new lung cancer related proteins that were positively correlated with their transcriptomics data. These proteins are mainly involved in cell migration, proliferation, and invasion.

Conclusion

The proteogenomic data on both NSCLC and SCLC cell lines presented in this manuscript is complementary to existing genomic and proteomic data related to lung cancers and will be crucial for a systems biology-level understanding of the molecular mechanism of lung cancers. The raw mass spectrometry data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository with the dataset identifier PXD015270.

Exon Definition Facilitates Reliable Control of Alternative Splicing in the RON Proto-Oncogene

Alternative splicing is a key step in eukaryotic gene expression that allows for the production of multiple transcript and protein isoforms from the same gene. Even though splicing is perturbed in many diseases, we currently lack insights into regulatory mechanisms promoting its precision and efficiency. We analyze high-throughput mutagenesis data obtained for an alternatively spliced exon in the proto-oncogene RON and determine the functional units that control this splicing event. Using mathematical modeling of distinct splicing mechanisms, we show that alternative splicing is based in RON on a so-called "exon definition" mechanism. Here, the recognition of the adjacent exons by the spliceosome is required for removal of an intron. We use our model to analyze the differences between the exon and intron definition scenarios and find that exon definition prevents the accumulation of deleterious, partially spliced retention products during alternative splicing regulation. Furthermore, it modularizes splicing control, as multiple regulatory inputs are integrated into a common net input, irrespective of the location and nature of the corresponding cis-regulatory elements in the pre-messenger RNA. Our analysis suggests that exon definition promotes robust and reliable splicing outcomes in RON splicing.

Is Hydrogen Sulfide a Concern During Treatment of Lung Adenocarcinoma With Ammonium Tetrathiomolybdate?

Ammonium tetrathiomolybdate (ATTM) has been used in breast cancer therapy for copper chelation, as elevated copper promotes tumor growth. ATTM is also an identified H₂S donor and endogenous H₂S facilitates VitB₁₂-induced S-adenosylmethionine (SAM) generation, which have been confirmed in m⁶A methylation and lung cancer development. The m⁶A modification was recently shown to participate in lung adenocarcinoma (LUAD) progression. These conflicting analyses of ATTM's anticancer vs. H₂S's carcinogenesis suggest that H₂S should not be ignored during LUAD's treatment with ATTM. This study was aimed to explore ATTM's effects on LUAD cells and mechanisms associated with H₂S and m⁶A. It was found that treatment with ATTM inhibited cell growth at high concentrations, while enhanced cell growth at low concentrations in three LUAD cell lines (A549, HCC827, and PC9). However, another copper chelator triethylenetetramine, without H₂S releasing activity, was not found to induce cell growth. Low ATTM concentrations also elevated m⁶A content in A549 cells. Analysis of differentially expressed genes in TCGA cohort indicated that m⁶A writer METTL3 and reader YTHDF1 were upregulated while eraser FTO was downregulated in LUAD tissues, consistent with the findings of protein expression in patient tissues. ATTM treatment of A549 cells significantly increased METTL3/14 and YTHDF1 while decreased FTO expression. Furthermore, inhibition of m⁶A with shMETTL3 RNA significantly attenuated eukaryotic translation initiation factor (eIF) expressions in A549 cells. Correlation analysis indicated that small nuclear ribonucleic protein PRPF6 was positively expressed with YTHDF1 in LUAD tissues. Knockdown of YTHDF1 partially blocked both basal and ATTM-induced PRPF6 expression, as well as A549 cell growth. Lastly, ATTM treatment not only raised intracellular H₂S content but also upregulated H₂S-producing enzymes. Exogenous H₂S application mimicked ATTM's aforementioned effects, but the effects could be weakened by zinc-induced H₂S scavenging. Collectively, H₂S impedes ATTM-induced anticancer effects through YTHDF1-dependent PRPF6 m⁶A methylation in lung adenocarcinoma cells.

A functional network of gastric-cancer-associated splicing events controlled by dysregulated splicing factors

Comprehensive understanding of aberrant splicing in gastric cancer is lacking. We RNA-sequenced 19 gastric tumor–normal pairs and identified 118 high-confidence tumor-associated (TA) alternative splicing events (ASEs) based on high-coverage sequencing and stringent filtering, and also identified 8 differentially expressed splicing factors (SFs). The TA ASEs occurred in genes primarily involved in cytoskeletal organization. We constructed a correlative network between TA ASE splicing ratios and SF expression, replicated it in independent gastric cancer data from The Cancer Genome Atlas and experimentally validated it by knockdown of the nodal SFs (PTBP1, ESRP2 and MBNL1). Each SF knockdown drove splicing alterations in several corresponding TA ASEs and led to alterations in cellular migration consistent with the role of TA ASEs in cytoskeletal organization. We have therefore established a robust network of dysregulated splicing associated with tumor invasion in gastric cancer. Our work is a resource for identifying oncogenic splice forms, SFs and splicing-generated tumor antigens as biomarkers and therapeutic targets.

Therapeutic potential of targeting mixed lineage kinases in cancer and inflammation

Dysregulation of intracellular signaling pathways is a key attribute of diseases associated with chronic inflammation, including cancer. Mitogen activated protein kinases have emerged as critical conduits of intracellular signal transmission, yet due to their ubiquitous roles in cellular processes, their direct inhibition may lead to undesired effects, thus limiting their usefulness as therapeutic targets. Mixed lineage kinases (MLKs) are mitogen-activated protein kinase kinase kinases (MAP3Ks) that interact with scaffolding proteins and function upstream of p38, JNK, ERK, and NF-kappaB to mediate diverse cellular signals. Studies involving gene silencing, genetically engineered mouse models, and small molecule inhibitors suggest that MLKs are critical in tumor progression as well as in inflammatory processes. Recent advances indicate that they may be useful targets in some types of cancer and in diseases driven by chronic inflammation including neurodegenerative diseases and metabolic diseases such as nonalcoholic steatohepatitis. This review describes existing MLK inhibitors, the roles of MLKs in various aspects of tumor progression and in the control of inflammatory processes, and the potential for therapeutic targeting of MLKs.

An Intricate Connection between Alternative Splicing and Phenotypic Plasticity in Development and Cancer

During tumor progression, hypoxia, nutrient deprivation or changes in the extracellular environment (i.e., induced by anti-cancer drugs) elicit adaptive responses in cancer cells. Cellular plasticity increases the chance that tumor cells may survive in a challenging microenvironment, acquire new mechanisms of resistance to conventional drugs, and spread to distant sites. Re-activation of stem pathways appears as a significant cause of cellular plasticity because it promotes the acquisition of stem-like properties through a profound phenotypic reprogramming of cancer cells. In addition, it is a major contributor to tumor heterogeneity, depending on the coexistence of phenotypically distinct subpopulations in the same tumor bulk. Several cellular mechanisms may drive this fundamental change, in particular, high-throughput sequencing technologies revealed a key role for alternative splicing (AS). Effectively, AS is one of the most important pre-mRNA processes that increases the diversity of transcriptome and proteome in a tissue- and development-dependent manner. Moreover, defective AS has been associated with several human diseases. However, its role in cancer cell plasticity and tumor heterogeneity remains unclear. Therefore, unravelling the intricate relationship between AS and the maintenance of a stem-like phenotype may explain molecular mechanisms underlying cancer cell plasticity and improve cancer diagnosis and treatment.

Multi-omic Dissection of Oncogenically Active Epiproteomes Identifies Drivers of Proliferative and Invasive Breast Tumors

Proliferative and invasive breast tumors evolve heterogeneously in individual patients, posing significant challenges in identifying new druggable targets for precision, effective therapy. Here we present a functional multi-omics method, interaction-Correlated Multi-omic Aberration Patterning (iC-MAP), which dissects intra-tumor heterogeneity and identifies in situ the oncogenic consequences of multi-omics aberrations that drive proliferative and invasive tumors. First, we perform chromatin activity-based chemoproteomics (ChaC) experiments on breast cancer (BC) patient tissues to identify genetic/transcriptomic alterations that manifest as oncogenically active proteins. ChaC employs a biotinylated small molecule probe that specifically binds to the oncogenically active histone methyltransferase G9a, enabling sorting/enrichment of a G9a-interacting protein complex that represents the predominant BC subtype in a tissue. Second, using patient transcriptomic/genomic data, we retrospectively identified some G9a interactor-encoding genes that showed individualized iC-MAP. Our iC-MAP findings represent both new diagnostic/prognostic markers to identify patient subsets with incurable metastatic disease and targets to create individualized therapeutic strategies.

•

ChaC dissects tumor heterogeneity for identifying oncogenic-active proteins

•

An oncogenic-active G9a-interactome represents the invasive tumor in a tissue

•

iC-MAP identifies multi-omics aberrations that drive invasive tumors

•

Patient-specific iC-MAP of select interactor genes are of prognostic value

Small RNA expression from viruses, bacteria and human miRNAs in colon cancer tissue and its association with microsatellite instability and tumor location

BACKGROUND

MicroRNAs (miRNA) and other small RNAs are frequently dysregulated in cancer and are promising biomarkers for colon cancer. Here we profile human, virus and bacteria small RNAs in normal and tumor tissue from early stage colon cancer and correlate the expression with clinical parameters.

METHODS

Small RNAs from colon cancer tissue and adjacent normal mucosa of 48 patients were sequenced using Illumina high-throughput sequencing. Clinical parameters were correlated with the small RNA expression data using linear models. We performed a meta-analysis by comparing publicly available small RNA sequencing datasets with our original sequencing data to confirm the main findings.

RESULTS

We identified 331 differentially expressed miRNAs between tumor and normal samples. We found that the major changes in miRNA expression between left and right colon are due to miRNAs located within the Hox-developmental genes, including miR-10b, miR-196b and miR-615. Further, we identified new miRNAs associated with microsatellite instability (MSI), including miR-335, miR-26 and miR-625. We performed a meta-analysis on all publicly available miRNA-seq datasets and identified 117 common miRNAs that were differentially expressed between tumor and normal tissue. The miRNAs miR-135b and miR-31 were the most significant upregulated miRNA in tumor across all datasets. The miRNA miR-133a was the most strongly downregulated miRNA in our dataset and also showed consistent downregulation in the other datasets. The miRNAs associated with MSI and tumor location in our data showed similar changes in the other datasets. Finally, we show that small RNAs from Epstein-Barr virus and Fusobacterium nucleatum are differentially expressed between tumor and normal adjacent tissue.

CONCLUSIONS

Small RNA profiling in colon cancer tissue revealed novel RNAs associated with MSI and tumor location. We show that Fusobacterium nucleatum are detectable at the RNA-level in colon tissue, and that both Fusobacterium nucleatum and Epstein-Barr virus separate tumor and normal tissue.

Splicing and cancer: Challenges and opportunities

Cancer arises from alterations in several metabolic processes affecting proliferation, growth, replication and death of cells. A fundamental challenge in the study of cancer biology is to uncover molecular mechanisms that lead to malignant cellular transformation. Recent genomic analyses revealed that many molecular alterations observed in cancers come from modifications in the splicing process, including mutations in pre-mRNA regulatory sequences, mutations in spliceosome components, and altered ratio of specific splicing regulators. While alterations in splice site preferences might generate alternative isoforms enabling different biological functions, these might also be responsible for nonfunctional isoforms that can eventually cause dysregulation in cellular processes. Molecular characteristics of regulatory sequences and proteins might also be important prognostic tools revealing a cancer-specific splicing pattern and linking splicing control to cancer development. The connection between cancer biology and splicing regulation is of primary importance to understand the mechanisms leading to disease and also to improve development of therapeutic approaches. Splicing modulation is being explored in new anti-cancer therapies and further investigation of targeted splicing factors is critical for the success of these strategies. This article is categorized under: RNA Processing > Splicing Mechanisms RNA-Based Catalysis > RNA Catalysis in Splicing and Translation RNA Processing > Splicing Regulation/Alternative Splicing RNA in Disease and Development > RNA in Disease.

Prognostic value and potential function of splicing events in prostate adenocarcinoma

Prostate adenocarcinoma (PRAD) is one of the most common types of malignancy in males and at present, effective prognostic indicators are limited. The development of PRAD has been associated with abnormalities in alternative splicing (AS), a requisite biological process of gene expression in eukaryotic cells; however, the prognostic value of AS products and splicing events remains to be elucidated. In the present study, the data of splicing events and the clinical information of PRAD patients were obtained from The Cancer Genome Atlas (TCGA)SpliceSeq and TCGA databases, respectively. A prognostic index (PI) was generated from disease-free survival-associated splicing events (DFS-SEs), which were identified by univariate/multivariate Cox regression analysis. A total of 6,909 DFS-SEs were identified in PRAD. The corresponding genes for the DFS-SEs were significantly enriched in mitochondria and their associated pathways according to Gene Ontology annotation and in the pathways of fatty acid metabolism, oxidative phosphorylation and Huntington's disease according to a Kyoto Encyclopedia of Genes and Genomes pathway analysis. The PI for mutually exclusive exons had the greatest ability to predict the probability of five-year disease-free survival of patients with PRAD, with an area under the time-dependent receiver-operating characteristic curve of 0.7606. Patients with PRAD, when divided into a 'low' and a 'high' group based on their median PI for exon skip values, exhibited a marked difference in disease-free survival (low vs. high, 3,588.45±250.51 vs. 1,531.08±136.50 days; P=7.43×10⁻⁹). A correlation network between DFS-SEs of splicing factors and non-splicing factors was constructed to determine the potential mechanisms in PRAD, which included the potential regulatory interaction between the splicing event of splicing factor RNA binding motif protein 5-alternate terminator (AT)-64957 and the splicing event of non-splicing factor heterochromatin protein 1 binding protein 3-AT-939. In conclusion, the PIs derived from DFS-SEs are valuable prognostic factors for patients with PRAD, and the function of splicing events in PRAD deserves further exploration.

The Cancer Spliceome: Reprograming of Alternative Splicing in Cancer

Alternative splicing allows for the expression of multiple RNA and protein isoforms from one gene, making it a major contributor to transcriptome and proteome diversification in eukaryotes. Advances in next generation sequencing technologies and genome-wide analyses have recently underscored the fact that the vast majority of multi-exon genes under normal physiology engage in alternative splicing in tissue-specific and developmental-specific manner. On the other hand, cancer cells exhibit remarkable transcriptome alterations partly by adopting cancer-specific splicing isoforms. These isoforms and their encoded proteins are not insignificant byproducts of the abnormal physiology of cancer cells, but either drivers of cancer progression or small but significant contributors to specific cancer hallmarks. Thus, it is paramount that the pathways that regulate alternative splicing in cancer, including the splicing factors that bind to pre-mRNAs and modulate spliceosome recruitment. In this review, we present a few distinct cases of alternative splicing in cancer, with an emphasis on their regulation as well as their contribution to cancer cell phenotype. Several categories of splicing aberrations are highlighted, including alterations in cancer-related genes that directly affect their pre-mRNA splicing, mutations in genes encoding splicing factors or core spliceosomal subunits, and the seemingly mutation-free disruptions in the balance of the expression of RNA-binding proteins, including components of both the major (U2-dependent) and minor (U12-dependent) spliceosomes. Given that the latter two classes cause global alterations in splicing that affect a wide range of genes, it remains a challenge to identify the ones that contribute to cancer progression. These challenges necessitate a systematic approach to decipher these aberrations and their impact on cancer. Ultimately, a sufficient understanding of splicing deregulation in cancer is predicted to pave the way for novel and innovative RNA-based therapies.

Alternative splicing events implicated in carcinogenesis and prognosis of colorectal cancer

Background: Emerging evidence suggested that aberrant alternative splicing (AS) is pervasive event in development and progression of cancer. However, the information of aberrant splicing events involved in colorectal carcinogenesis and progression is still elusive. Materials and Methods: In this study, splicing data of 499 colon adenocarcinoma cases (COAD) and 176 rectum adenocarcinoma (READ) with clinicopathological information were obtained from The Cancer Genome Atlas (TCGA) to explore the changes of alternative splicing events in relation to the carcinogenesis and prognosis of colorectal cancer (CRC). Gene interaction network construction, functional and pathway enrichment analysis were performed by multiple bioinformatics tools. Results: Overall, most AS patterns were more active in CRC tissues than adjacent normal ones. We detected altogether 35391 AS events of 9084 genes in COAD and 34900 AS events of 9032 genes in READ, some of which were differentially spliced between cancer tissues and normal tissues including genes of SULT1A2, CALD1, DTNA, COL12A1 and TTLL12. Differentially spliced genes were enriched in biological process including muscle organ development, cytoskeleton organization, actin cytoskeleton organization, biological adhesion, and cell adhesion. The integrated predictor model of COAD showed an AUC of 0.805 (sensitivity: 0.734; specificity: 0.756) while READ predictor had an AUC of 0.738 (sensitivity: 0.614; specificity: 0.900). In addition, a number of prognosis-associated AS events were discovered, including genes of PSMD2, NOL8, ALDH4A1, SLC10A7 and PPAT. Conclusion: We draw comprehensive profiles of alternative splicing events in the carcinogenesis and prognosis of CRC. The interaction network and functional connections were constructed to elucidate the underlying mechanisms of alternative splicing in CRC.

The Circular RNA Profiles of Colorectal Tumor Metastatic Cells

Circular RNAs (circRNAs) have been reported that can be used as biomarkers for colorectal cancers (CRC) and other types of tumors. However, a limited number of studies have been performed investigating the potential role of circRNAs in tumor metastasis. Here, we examined the circRNAs in two CRC cell lines (a primary tumor cell SW480 and its metastasis cell SW620), and found a large set of circRNA (2,919 ncDECs) with significantly differential expression patterns relative to normal cells (NCM460). In addition, we uncovered a set of 623 pmDECs that differ between the primary CRC cells and its metastasis cells. Both differentially expressed circRNA (DEC) sets contain many previously unknown putative CRC-related circRNAs, thereby providing many new circRNAs as candidate biomarkers for CRC development and metastasis. These studies are the first large-scale identification of metastasis-related circRNAs for CRC and provide valuable candidate biomarkers for diagnostic and a starting point for additional investigations of CRC metastasis.

Mixed lineage kinase ZAK promotes epithelial-mesenchymal transition in cancer progression

ZAK, a mixed lineage kinase, is often described as a positive or negative regulator of cell growth. We identified it as one of the top hits in our kinome cDNA screen for potent regulators of epithelial mesenchymal transition (EMT). Ectopic expression of ZAK promoted EMT phenotypes and apoptosis resistance in multiple epithelial cell lines, while having different impacts on cell growth in different cell lines. Conversely, depletion of ZAK in aggressive mesenchymal cancer cells reversed EMT phenotypes, increased sensitivity to conventional cytotoxic drugs, and attenuated bone metastasis potential, with little impact on primary tumor growth. Mechanistically, ZAK-mediated EMT is associated with activation of ZEB1 and suppression of epithelial splicing regulatory proteins (ESRPs), which results in a switch in CD44 expression from the epithelial CD44v8-9 isoform to the mesenchymal CD44s isoform. Of note, transcriptomic analysis showed that ZAK overexpression is significantly associated with poor survival in a number of human cancer types. Tissue microarray analysis on breast invasive carcinoma further supported that ZAK overexpression is an independent poor prognostic factor for overall survival in breast cancer. Through combination with ZAK, prognostic accuracy of other common clinicopathological markers in breast cancer is improved by up to 21%. Taken together, these results suggest that promoting EMT is the primary role for ZAK in cancer progression. They also highlight its potential as a biomarker to identify high-risk patients, and suggest its promise as a therapeutic target for inhibiting metastasis and overcoming drug resistance.

Regulatory Potential of the RNA Processing Machinery: Implications for Human Disease

Splicing and nuclear export of mRNA are critical steps in the gene expression pathway. While RNA processing factors can perform general, essential functions for intron removal and bulk export of mRNA, emerging evidence highlights that the core RNA splicing and export machineries also display regulatory potential. Here, we discuss recent insights into how this regulatory potential can selectively alter gene expression and regulate important biological processes. We also highlight the participation of RNA processing pathways in the cellular response to DNA damage at multiple levels. These findings have important implications for the contribution of selective mRNA processing and export to the development of human cancers and neurodegenerative disorders.

snRNP proteins in health and disease

Split gene architecture of most human genes requires removal of intervening sequences by mRNA splicing that occurs on large multiprotein complexes called spliceosomes. Mutations compromising several spliceosomal components have been recorded in degenerative syndromes and haematological neoplasia, thereby highlighting the importance of accurate splicing execution in homeostasis of assorted adult tissues. Moreover, insufficient splicing underlies defective development of craniofacial skeleton and upper extremities. This review summarizes recent advances in the understanding of splicing factor function deduced from cryo-EM structures. We combine these data with the characterization of splicing factors implicated in hereditary or somatic disorders, with a focus on potential functional consequences the mutations may elicit in spliceosome assembly and/or performance. Given aberrant splicing or perturbations in splicing efficiency substantially underpin disease pathogenesis, profound understanding of the mis-splicing principles may open new therapeutic vistas. In three major sections dedicated to retinal dystrophies, hereditary acrofacial syndromes, and haematological malignancies, we delineate the noticeable variety of conditions associated with dysfunctional splicing and accentuate recurrent patterns in splicing defects.

RNA processing as an alternative route to attack glioblastoma

Genomic analyses have become an important tool to identify new avenues for therapy. This is especially true for cancer types with extremely poor outcomes, since our lack of effective therapies offers no tangible clinical starting point to build upon. The highly malignant brain tumor glioblastoma (GBM) exemplifies such a refractory cancer, with only 15 month average patient survival. Analyses of several hundred GBM samples compiled by the TCGA (The Cancer Genome Atlas) have produced an extensive transcriptomic map, identified prevalent chromosomal alterations, and defined important driver mutations. Unfortunately, clinical trials based on these results have not yet delivered an improvement on outcome. It is, therefore, necessary to characterize other regulatory routes known for playing a role in tumor relapse and response to treatment. Alternative splicing affects more than 90% of the human coding genes and it is an important source for transcript variation and gene regulation. Mutations and alterations in splicing factors are highly prevalent in multiple cancers, demonstrating the potential for splicing to act as a tumor driver. As a result, numerous genes are expressed as cancer-specific splicing isoforms that are functionally distinct from the canonical isoforms found in normal tissue. These include genes that regulate cancer-critical pathways such as apoptosis, DNA repair, cell proliferation, and migration. Splicing defects can even induce genomic instability, a common characteristic of cancer, and a driver of tumor evolution. Importantly, components of the splicing machinery are targetable; multiple drugs can inhibit splicing factors or promote changes in splicing which could be exploited to begin improving clinical outcomes. Here, we review the current literature and present a case for exploring RNA processing as therapeutic route for the treatment of GBM.

RNA editing-dependent epitranscriptome diversity in cancer stem cells

Cancer stem cells (CSCs) can regenerate all facets of a tumour as a result of their stem cell-like capacity to self-renew, survive and become dormant in protective microenvironments. CSCs evolve during tumour progression in a manner that conforms to Charles Darwin's principle of natural selection. Although somatic DNA mutations and epigenetic alterations promote evolution, post-transcriptional RNA modifications together with RNA binding protein activity (the 'epitranscriptome') might also contribute to clonal evolution through dynamic determination of RNA function and gene expression diversity in response to environmental stimuli. Deregulation of these epitranscriptomic events contributes to CSC generation and maintenance, which governs cancer progression and drug resistance. In this Review, we discuss the role of malignant RNA processing in CSC generation and maintenance, including mechanisms of RNA methylation, RNA editing and RNA splicing, and the functional consequences of their aberrant regulation in human malignancies. Finally, we highlight the potential of these events as novel CSC biomarkers as well as therapeutic targets.

PRPF overexpression induces drug resistance through actin cytoskeleton rearrangement and epithelial-mesenchymal transition

Pre-mRNA processing factor (PRPF) 4B kinase belongs to the CDK-like kinase family, and is involved in pre-mRNA splicing, and in signal transduction. In this study, we observed that PRPF overexpression decreased the intracellular levels of reactive oxygen species, and inhibited resveratrol-induced apoptosis by activating the cell survival signaling proteins NFκB, ERK, and c-MYC in HCT116 human colon cancer cells. PRPF overexpression altered cellular morphology, and rearranged the actin cytoskeleton, by regulating the activity of Rho family proteins. Moreover, it decreased the activity of RhoA, but increased the expression of Rac1. In addition, PRPF triggered the epithelial-mesenchymal transition (EMT), and decreased the invasiveness of HCT116, PC3 human prostate, and B16-F10 melanoma cells. The loss of E-cadherin, a hallmark of EMT, was observed in HCT116 cells overexpressing PRPF. Taken together, these results indicate that PRPF blocks the apoptotic effects of resveratrol by activating cell survival signaling pathways, rearranging the actin cytoskeleton, and inducing EMT. The elucidation of the mechanisms that underlie anticancer drug resistance and the anti-apoptosis effect of PRPF may provide a therapeutic basis for inhibiting tumor growth and preventing metastasis in various cancers.

Associations between genetic variants in mRNA splicing-related genes and risk of lung cancer: a pathway-based analysis from published GWASs

mRNA splicing is an important mechanism to regulate mRNA expression. Abnormal regulation of this process may lead to lung cancer. Here, we investigated the associations of 11,966 single-nucleotide polymorphisms (SNPs) in 206 mRNA splicing-related genes with lung cancer risk by using the summary data from six published genome-wide association studies (GWASs) of Transdisciplinary Research in Cancer of the Lung (TRICL) (12,160 cases and 16,838 controls) and another two lung cancer GWASs of Harvard University (984 cases and 970 controls) and deCODE (1,319 cases and 26,380 controls). We found that a total of 12 significant SNPs with false discovery rate (FDR) ≤0.05 were mapped to one novel gene PRPF6 and two previously reported genes (DHX16 and LSM2) that were also confirmed in this study. The six novel SNPs in PRPF6 were in high linkage disequilibrium and associated with PRPF6 mRNA expression in lymphoblastoid cells from 373 Europeans in the 1000 Genomes Project. Taken together, our studies shed new light on the role of mRNA splicing genes in the development of lung cancer.

Robust stratification of breast cancer subtypes using differential patterns of transcript isoform expression

Breast cancer, the second leading cause of cancer death of women worldwide, is a heterogenous disease with multiple different subtypes. These subtypes carry important implications for prognosis and therapy. Interestingly, it is known that these different subtypes not only have different biological behaviors, but also have distinct gene expression profiles. However, it has not been rigorously explored whether particular transcriptional isoforms are also differentially expressed among breast cancer subtypes, or whether transcript isoforms from the same sets of genes can be used to differentiate subtypes. To address these questions, we analyzed the patterns of transcript isoform expression using a small set of RNA-sequencing data for eleven Estrogen Receptor positive (ER+) subtype and fourteen triple negative (TN) subtype tumors. We identified specific sets of isoforms that distinguish these tumor subtypes with higher fidelity than standard mRNA expression profiles. We found that alternate promoter usage, alternative splicing, and alternate 3’UTR usage are differentially regulated in breast cancer subtypes. Profiling of isoform expression in a second, independent cohort of 68 tumors confirmed that expression of splice isoforms differentiates breast cancer subtypes. Furthermore, analysis of RNAseq data from 594 cases from the TCGA cohort confirmed the ability of isoform usage to distinguish breast cancer subtypes. Also using our expression data, we identified several RNA processing factors that were differentially expressed between tumor subtypes and/or regulated by estrogen receptor, including YBX1, YBX2, MAGOH, MAGOHB, and PCBP2. RNAi knock-down of these RNA processing factors in MCF7 cells altered isoform expression. These results indicate that global dysregulation of splicing in breast cancer occurs in a subtype-specific and reproducible manner and is driven by specific differentially expressed RNA processing factors.

Breast cancer, the second leading cause of cancer death of women worldwide, is a heterogenous disease. Different subtypes of breast cancer display very different expression programs, and these expression programs are associated with different patient outcomes and with different treatment protocols. However, little is known about what drives these subtype differences. By sequencing RNA in a discovery cohort of breast cancer patients, we demonstrate that different subtypes of breast cancer can be distinguished by simply using differential transcript isoform expression. We confirmed our findings using two additional patient cohorts. We also demonstrate that differential expression of RNA processing factors between subtypes can affect differences in isoform usage. Using RNAi we knock down differentially expressed RNA processing factors including YBX1, YBX2, MAGOH, MAGOHB, and PCBP2, and show that this knock-down results in differential isoform expression of the genes identified in our disease subtype panel. Taken together, our results indicate that global dysregulation of splicing occurs in a subtype-specific and reproducible manner in breast cancer, and is driven by specific differentially expressed RNA processing factors.

Potential effect of spliceosome inhibition in small cell lung cancer irrespective of the MYC status

Small cell lung cancer (SCLC) is a highly aggressive malignancy with few therapeutic advances in the treatment in recent decades. Based on a recent study that identified the spliceosome as a therapeutic vulnerability in MYC-driven breast cancers, we evaluated the efficacy of a spliceosome inhibitor in SCLC cell lines and analyzed the correlation with MYC status. Among 23 SCLC cell lines examined, eight showed high MYC protein expression (> 80% positive cells) by immunohistochemistry (IHC), while 10 cell lines demonstrated no staining for MYC. The remaining five cell lines showed weak staining (< 40% positive cells). All four cell lines that were previously demonstrated to have MYC gene amplification were positive for MYC by IHC. Four cell lines with high MYC expression and four with low MYC expression were used in further analysis. A spliceosome inhibitor, pladienolide B, showed high efficacy (IC₅₀ < 12nM) in all eight cell lines tested, irrespective of the MYC IHC or MYC gene amplification status. We observed that the four cell lines with higher sensitivity to the spliceosome inhibitor were established from patients with prior chemotherapy. Therefore we chronically treated H1048 cells, that were established from a treatment-naïve patient, with cisplatin for 4 weeks, and found that H1048-cisplatin treated cells became more sensitive to pladienolide B. In conclusion, our in vitro results indicate that spliceosome inhibitors would be promising molecular target drugs in SCLC irrespective of the MYC status, especially in the second-line settings after an effective front-line chemotherapy.

Transcriptomic evaluation of bovine blastocysts obtained from peri-pubertal oocyte donors

Assisted reproduction technologies (ART) and high selection pressure in the dairy industry are leading towards the use of younger females for reproduction, thereby reducing the interval between generations. This situation may have a negative impact on embryo quality, thus reducing the success rate of the procedures. This study aimed to document the effects of oocyte donor age on embryo quality, at the transcriptomic level, in order to characterize the effects of using young females for reproduction purpose. Young Holstein heifers (n = 10) were used at three different ages for ovarian stimulation protocols and oocyte collections (at 8, 11 and 14 months). All of the oocytes were fertilized in vitro with the semen of one adult bull, generating three lots of embryos per animal. Each animal was its own control for the evaluation of the effects of age. The EmbryoGENE platform was used for the assessment of gene expression patterns at the blastocyst stage. Embryos from animals at 8 vs 14 months and at 11 vs 14 months were used for microarray hybridization. Validation was done by performing RT-qPCR on seven candidate genes. Age-related contrast analysis (8 vs 14 mo and 11 vs 14 mo) identified 242 differentially expressed genes (DEGs) for the first contrast, and 296 for the second. The analysis of the molecular and biological functions of the DEGs suggests a metabolic cause to explain the differences that are observed between embryos from immature and adult subjects. The mTOR and PPAR signaling pathways, as well as the NRF2-mediated oxidative stress response pathways were among the gene expression pathways affected by donor age. In conclusion, the main differences between embryos produced at peri-pubertal ages are related to metabolic conditions resulting in a higher impact of in vitro conditions on blastocyts from younger heifers.