AS-CMC: a pan-cancer database of alternative splicing for molecular classification of cancer

Alternative splicing (AS) is a post-transcriptional regulation that leads to the complexity of the transcriptome. Despite the growing importance of AS in cancer research, the role of AS has not been systematically studied, especially in understanding cancer molecular classification. Herein, we analyzed the molecular subtype-specific regulation of AS using The Cancer Genome Atlas data and constructed a web-based database, named Alternative Splicing for Cancer Molecular Classification (AS-CMC). Our system harbors three analysis modules for exploring subtype-specific AS events, evaluating their phenotype association, and performing pan-cancer comparison. The number of subtype-specific AS events was found to be diverse across cancer types, and some differentially regulated AS events were recurrently found in multiple cancer types. We analyzed a subtype-specific AS in exon 11 of mitogen-activated protein kinase kinase 7 (MAP3K7) as an example of a pan-cancer AS biomarker. This AS marker showed significant association with the survival of patients with stomach adenocarcinoma. Our analysis revealed AS as an important determinant for cancer molecular classification. AS-CMC is the first web-based resource that provides a comprehensive tool to explore the biological implications of AS events, facilitating the discovery of novel AS biomarkers.

Matrix factorization-based data fusion for the prediction of RNA-binding proteins and alternative splicing event associations during epithelial-mesenchymal transition

MOTIVATION

The epithelial-mesenchymal transition (EMT) is a cellular-developmental process activated during tumor metastasis. Transcriptional regulatory networks controlling EMT are well studied; however, alternative RNA splicing also plays a critical regulatory role during this process. Unfortunately, a comprehensive understanding of alternative splicing (AS) and the RNA-binding proteins (RBPs) that regulate it during EMT remains largely unknown. Therefore, a great need exists to develop effective computational methods for predicting associations of RBPs and AS events. Dramatically increasing data sources that have direct and indirect information associated with RBPs and AS events have provided an ideal platform for inferring these associations.

RESULTS

In this study, we propose a novel method for RBP-AS target prediction based on weighted data fusion with sparse matrix tri-factorization (WDFSMF in short) that simultaneously decomposes heterogeneous data source matrices into low-rank matrices to reveal hidden associations. WDFSMF can select and integrate data sources by assigning different weights to those sources, and these weights can be assigned automatically. In addition, WDFSMF can identify significant RBP complexes regulating AS events and eliminate noise and outliers from the data. Our proposed method achieves an area under the receiver operating characteristic curve (AUC) of $90.78\%$, which shows that WDFSMF can effectively predict RBP-AS event associations with higher accuracy compared with previous methods. Furthermore, this study identifies significant RBPs as complexes for AS events during EMT and provides solid ground for further investigation into RNA regulation during EMT and metastasis. WDFSMF is a general data fusion framework, and as such it can also be adapted to predict associations between other biological entities.

RBFOX2 alters splicing outcome in distinct binding modes with multiple protein partners

RBFOX2 controls the splicing of a large number of transcripts implicated in cell differentiation and development. Parsing RNA-binding protein datasets, we uncover that RBFOX2 can interact with hnRNPC, hnRNPM and SRSF1 to regulate splicing of a broad range of splicing events using different sequence motifs and binding modes. Using immunoprecipitation, specific RBP knockdown, RNA-seq and splice-sensitive PCR, we show that RBFOX2 can target splice sites using three binding configurations: single, multiple or secondary modes. In the single binding mode RBFOX2 is recruited to its target splice sites through a single canonical binding motif, while in the multiple binding mode RBFOX2 binding sites include the adjacent binding of at least one other RNA binding protein partner. Finally, in the secondary binding mode RBFOX2 likely does not bind the RNA directly but is recruited to splice sites lacking its canonical binding motif through the binding of one of its protein partners. These dynamic modes bind distinct sets of transcripts at different positions and distances relative to alternative splice sites explaining the heterogeneity of RBFOX2 targets and splicing outcomes.

A combinatorially regulated RNA splicing signature predicts breast cancer EMT states and patient survival

During breast cancer metastasis, the developmental process epithelial-mesenchymal transition (EMT) is abnormally activated. Transcriptional regulatory networks controlling EMT are well-studied; however, alternative RNA splicing also plays a critical regulatory role during this process. A comprehensive understanding of alternative splicing (AS) and the RNA binding proteins (RBPs) that regulate it during EMT and their impact on breast cancer remains largely unknown. In this study, we annotated AS in the breast cancer TCGA data set and identified an AS signature that is capable of distinguishing epithelial and mesenchymal states of the tumors. This AS signature contains 25 AS events, among which nine showed increased exon inclusion and 16 showed exon skipping during EMT. This AS signature accurately assigns the EMT status of cells in the CCLE data set and robustly predicts patient survival. We further developed an effective computational method using bipartite networks to identify RBP-AS networks during EMT. This network analysis revealed the complexity of RBP regulation and nominated previously unknown RBPs that regulate EMT-associated AS events. This study highlights the importance of global AS regulation during EMT in cancer progression and paves the way for further investigation into RNA regulation in EMT and metastasis.

Regeneration of Damaged Tendon-Bone Junctions (Entheses)-TAK1 as a Potential Node Factor

Musculoskeletal dysfunctions are highly prevalent due to increasing life expectancy. Consequently, novel solutions to optimize treatment of patients are required. The current major research focus is to develop innovative concepts for single tissues. However, interest is also emerging to generate applications for tissue transitions where highly divergent properties need to work together, as in bone-cartilage or bone-tendon transitions. Finding medical solutions for dysfunctions of such tissue transitions presents an added challenge, both in research and in clinics. This review aims to provide an overview of the anatomical structure of healthy adult entheses and their development during embryogenesis. Subsequently, important scientific progress in restoration of damaged entheses is presented. With respect to enthesis dysfunction, the review further focuses on inflammation. Although molecular, cellular and tissue mechanisms during inflammation are well understood, tissue regeneration in context of inflammation still presents an unmet clinical need and goes along with unresolved biological questions. Furthermore, this review gives particular attention to the potential role of a signaling mediator protein, transforming growth factor beta-activated kinase-1 (TAK1), which is at the node of regenerative and inflammatory signaling and is one example for a less regarded aspect and potential important link between tissue regeneration and inflammation.

Rbfox-Splicing Factors Maintain Skeletal Muscle Mass by Regulating Calpain3 and Proteostasis

Maintenance of skeletal muscle mass requires a dynamic balance between protein synthesis and tightly controlled protein degradation by the calpain, autophagy-lysosome, and ubiquitin-proteasome systems (proteostasis). Several sensing and gene-regulatory mechanisms act together to maintain this balance in response to changing conditions. Here, we show that deletion of the highly conserved Rbfox1 and Rbfox2 alternative splicing regulators in adult mouse skeletal muscle causes rapid, severe loss of muscle mass. Rbfox deletion did not cause a reduction in global protein synthesis, but it led to altered splicing of hundreds of gene transcripts, including capn3, which produced an active form of calpain3 protease. Rbfox knockout also led to a reduction in autophagy flux, likely producing a compensatory increase in general protein degradation by the proteasome. Our results indicate that the Rbfox-splicing factors are essential for the maintenance of skeletal muscle mass and proteostasis.

Conserved functions of RNA-binding proteins in muscle

Animals require different types of muscle for survival, for example for circulation, motility, reproduction and digestion. Much emphasis in the muscle field has been placed on understanding how transcriptional regulation generates diverse types of muscle during development. Recent work indicates that alternative splicing and RNA regulation are as critical to muscle development, and altered function of RNA-binding proteins causes muscle disease. Although hundreds of genes predicted to bind RNA are expressed in muscles, many fewer have been functionally characterized. We present a cross-species view summarizing what is known about RNA-binding protein function in muscle, from worms and flies to zebrafish, mice and humans. In particular, we focus on alternative splicing regulated by the CELF, MBNL and RBFOX families of proteins. We discuss the systemic nature of diseases associated with loss of RNA-binding proteins in muscle, focusing on mis-regulation of CELF and MBNL in myotonic dystrophy. These examples illustrate the conservation of RNA-binding protein function and the marked utility of genetic model systems in understanding mechanisms of RNA regulation.

TGF-β-induced alternative splicing of TAK1 promotes EMT and drug resistance

Transforming growth factor-β (TGF-β) is major inducer of epithelial to mesenchymal transition (EMT), which associates with cancer cell metastasis and resistance to chemotherapy and targeted drugs, through both transcriptional and non-transcriptional mechanisms. We previously reported that in cancer cells, heightened mitogenic signaling allows TGF-β-activated Smad3 to interact with poly(RC) binding protein 1 (PCBP1) and together they regulate many alternative splicing events that favors expression of protein isoforms essential for EMT, cytoskeletal rearrangement, and adherens junction signaling. Here, we show that the exclusion of TGF-β-activated kinase 1 (TAK1) variable exon 12 requires another RNA-binding protein, Fox-1 homolog 2 (Rbfox2), which binds intronic sequences in front of exon 12 independently of the Smad3-PCBP1 complex. Functionally, exon 12-excluded TAK1∆E12 and full length TAK1FL are distinct. The short isoform TAK1∆E12 is constitutively active and supports TGF-β-induced EMT and nuclear factor kappa B (NF-κB) signaling, whereas the full-length isoform TAK1FL promotes TGF-β-induced apoptosis. These observations offer a harmonious explanation for how a single TAK1 kinase can mediate the opposing responses of cell survival and apoptosis in response to TGF-β. They also reveal a propensity of the alternatively spliced TAK1 isoform TAK1∆E12 to cause drug resistance due to its activity in supporting EMT and NF-κB survival signaling.

Regulation of Tak1 alternative splicing by splice-switching oligonucleotides

Alternative splicing (AS) generates multiple isoforms from a single precursor mRNA, and these isoforms usually exhibit different tissue distributions and functions. Aberrant protein isoforms can lead to abnormalities in protein function and may even result in genetic disorders or cancer. In recent years, splice-switching oligonucleotides (SSOs) have emerged as a promising therapeutic strategy for several neurological diseases, but the efficacy of this strategy in other organs is less reported. In this study, we designed and synthesized SSOs targeting the splicing regulators of exon 12 of the Tak1 gene, inducing variant switching between Tak1-A and Tak1-B. We also designed SSOs capable of knockdown both Tak1 variants by inducing the aberrant splicing of exon 4. The Vivo-morpholino SSOs showed significant splice-switching of Tak1 in mouse liver, with a persistence of at least 10 days after initial SSOs delivery. Bioinformatics analysis indicated a lipid metabolism-related function for Tak1-B but not Tak1-A. The conversion of Tak1-B to Tak1-A consistently led to significant accumulation of lipids in cultured AML12 cells, as well as the dysregulation of several lipid metabolism-related genes in mouse liver. Different functional properties of the two isoforms may explain the conflicting functions previously reported for Tak1. In conclusion, our research clarified the different functions of Tak1 isoforms, and provided an efficient strategy for the functional research of the AS isoforms.

New inter-correlated genes targeted by diatom-derived polyunsaturated aldehydes in the sea urchin Paracentrotus lividus

The marine environment is continually subjected to the action of stressors (including natural toxins), which represent a constant danger for benthic communities. In the present work using network analysis we identified ten genes on the basis of associated functions (FOXA, FoxG, GFI-1, nodal, JNK, OneCut/Hnf6, TAK1, tcf4, TCF7, VEGF) in the sea urchin Paracentrotus lividus, having key roles in different processes, such as embryonic development and asymmetry, cell fate specification, cell differentiation and morphogenesis, and skeletogenesis. These genes are correlated with three HUB genes, Foxo, Jun and HIF1A. Real Time qPCR revealed that during sea urchin embryonic development the expression levels of these genes were modulated by three diatom-derived polyunsaturated aldehydes (PUAs), decadienal, heptadienal and octadienal. Our findings show how changes in gene expression levels may be used as an early indicator of stressful conditions in the marine environment. The identification of key genes and the molecular pathways in which they are involved represents a fundamental tool in understanding how marine organisms try to afford protection against toxicants, to avoid deleterious consequences and irreversible damages. The genes identified in this work as targets for PUAs can be considered as possible biomarkers to detect exposure to different environmental pollutants.

Identified of a novel cis-element regulating the alternative splicing of LcDREB2

Alternative splicing (AS) is an important gene regulation mechanism in plants. Despite the widespread use of AS in plant gene expression regulation, the identification of the cis-elements involved in the AS mechanism is rarely reported in plants. To explore the regulation mechanism of the AS of LcDREB2, a DREB2 ortholog from Sheepgrass (Leymus chinensis), the genomic sequences of LcDREB2 and its homologs in Poaceae were aligned, and six mutations were introduced in the conserved sequence of LcDREB2. By analyzing the distinct transcript patterns of the LcDREB2 mutants in transgenic Oryza sativa, a novel cis-element that affected the AS of LcDREB2 was identified as Exonic Splicing Enhancer 1 (ESE1). In addition, five serine-arginine rich (SR) proteins were confirmed to interact with ESE1 by electrophoretic mobility shift assay (EMSA). To further explore the expression regulation mechanism of the DREB subfamily, phylogenetic analysis of DREB2 paralogous genes was performed. The results strongly supported the hypothesis that AS is conserved in Poaceae plants and that it is an evolutionary strategy for the regulation of the functional expression of genes. The findings and methods of our study will promote a substantial step forward in understanding of the plant AS regulation mechanism.

Neural differentiation modulates the vertebrate brain specific splicing program

Alternative splicing patterns are known to vary between tissues but these patterns have been found to be predominantly peculiar to one species or another, implying only a limited function in fundamental neural biology. Here we used high-throughput RT-PCR to monitor the expression pattern of all the annotated simple alternative splicing events (ASEs) in the Reference Sequence Database, in different mouse tissues and identified 93 brain-specific events that shift from one isoform to another (switch-like) between brain and other tissues. Consistent with an important function, regulation of a core set of 9 conserved switch-like ASEs is highly conserved, as they have the same pattern of tissue-specific splicing in all vertebrates tested: human, mouse and zebrafish. Several of these ASEs are embedded within genes that encode proteins associated with the neuronal microtubule network, and show a dramatic and concerted shift within a short time window of human neural stem cell differentiation. Similarly these exons are dynamically regulated in zebrafish development. These data demonstrate that although alternative splicing patterns often vary between species, there is nonetheless a core set of vertebrate brain-specific ASEs that are conserved between species and associated with neural differentiation.

Analysis of pollen-specific alternative splicing in Arabidopsis thaliana via semi-quantitative PCR

Alternative splicing enables a single gene to produce multiple mRNA isoforms by varying splice site selection. In animals, alternative splicing of mRNA isoforms between cell types is widespread and supports cellular differentiation. In plants, at least 20% of multi-exon genes are alternatively spliced, but the extent and significance of tissue-specific splicing is less well understood, partly because it is difficult to isolate cells of a single type. Pollen is a useful model system to study tissue-specific splicing in higher plants because pollen grains contain only two cell types and can be collected in large amounts without damaging cells. Previously, we identified pollen-specific splicing patterns by comparing RNA-Seq data from Arabidopsis pollen and leaves. Here, we used semi-quantitative PCR to validate pollen-specific splicing patterns among genes where RNA-Seq data analysis indicated splicing was most different between pollen and leaves. PCR testing confirmed eight of nine alternative splicing patterns, and results from the ninth were inconclusive. In four genes, alternative transcriptional start sites coincided with alternative splicing. This study highlights the value of the low-cost PCR assay as a method of validating RNA-Seq results.

MBNL1 and RBFOX2 cooperate to establish a splicing programme involved in pluripotent stem cell differentiation

Reprogramming somatic cells into induced pluripotent stem cells (iPSCs) has provided huge insight into the pathways, mechanisms and transcription factors that control differentiation. Here we use high-throughput RT-PCR technology to take a snapshot of splicing changes in the full spectrum of high- and low-expressed genes during induction of fibroblasts, from several donors, into iPSCs and their subsequent redifferentiation. We uncover a programme of concerted alternative splicing changes involved in late mesoderm differentiation and controlled by key splicing regulators MBNL1 and RBFOX2. These critical splicing adjustments arise early in vertebrate evolution and remain fixed in at least 10 genes (including PLOD2, CLSTN1, ATP2A1, PALM, ITGA6, KIF13A, FMNL3, PPIP5K1, MARK2 and FNIP1), implying that vertebrates require alternative splicing to fully implement the instructions of transcriptional control networks.

Illuminating the Transcriptome through the Genome

Sequencing the human genome was a huge milestone in genetic research that revealed almost the total DNA sequence required to create a human being. However, in order to function, the DNA genome needs to be expressed as an RNA transcriptome. This article reviews how knowledge of genome sequence information has led to fundamental discoveries in how the transcriptome is processed, with a focus on new system-wide insights into how pre-mRNAs that are encoded by split genes in the genome are rearranged by splicing into functional mRNAs. These advances have been made possible by the development of new post-genome technologies to probe splicing patterns. Transcriptome-wide approaches have characterised a "splicing code" that is embedded within and has a significant role in deciphering the genome, and is deciphered by RNA binding proteins. These analyses have also found that most human genes encode multiple mRNA isoforms, and in some cases proteins, leading in turn to a re-assessment of what exactly a gene is. Analysis of the transcriptome has given insights into how the genome is packaged and transcribed, and is helping to explain important aspects of genome evolution.

RNA binding proteins in the regulation of heart development

In vivo, RNA molecules are constantly accompanied by RNA binding proteins (RBPs), which are intimately involved in every step of RNA biology, including transcription, editing, splicing, transport and localization, stability, and translation. RBPs therefore have opportunities to shape gene expression at multiple levels. This capacity is particularly important during development, when dynamic chemical and physical changes give rise to complex organs and tissues. This review discusses RBPs in the context of heart development. Since the targets and functions of most RBPs--in the heart and at large--are not fully understood, this review focuses on the expression and roles of RBPs that have been implicated in specific stages of heart development or developmental pathology. RBPs are involved in nearly every stage of cardiogenesis, including the formation, morphogenesis, and maturation of the heart. A fuller understanding of the roles and substrates of these proteins could ultimately provide attractive targets for the design of therapies for congenital heart defects, cardiovascular disease, or cardiac tissue repair.

RNA-seq of Arabidopsis pollen uncovers novel transcription and alternative splicing

Pollen grains of Arabidopsis (Arabidopsis thaliana) contain two haploid sperm cells enclosed in a haploid vegetative cell. Upon germination, the vegetative cell extrudes a pollen tube that carries the sperm to an ovule for fertilization. Knowing the identity, relative abundance, and splicing patterns of pollen transcripts will improve our understanding of pollen and allow investigation of tissue-specific splicing in plants. Most Arabidopsis pollen transcriptome studies have used the ATH1 microarray, which does not assay splice variants and lacks specific probe sets for many genes. To investigate the pollen transcriptome, we performed high-throughput sequencing (RNA-Seq) of Arabidopsis pollen and seedlings for comparison. Gene expression was more diverse in seedling, and genes involved in cell wall biogenesis were highly expressed in pollen. RNA-Seq detected at least 4,172 protein-coding genes expressed in pollen, including 289 assayed only by nonspecific probe sets. Additional exons and previously unannotated 5' and 3' untranslated regions for pollen-expressed genes were revealed. We detected regions in the genome not previously annotated as expressed; 14 were tested and 12 were confirmed by polymerase chain reaction. Gapped read alignments revealed 1,908 high-confidence new splicing events supported by 10 or more spliced read alignments. Alternative splicing patterns in pollen and seedling were highly correlated. For most alternatively spliced genes, the ratio of variants in pollen and seedling was similar, except for some encoding proteins involved in RNA splicing. This study highlights the robustness of splicing patterns in plants and the importance of ongoing annotation and visualization of RNA-Seq data using interactive tools such as Integrated Genome Browser.