Tissue-specific alternative splicing remodels protein-protein interaction networks

Widespread Expansion of Protein Interaction Capabilities by Alternative Splicing

While alternative splicing is known to diversify the functional characteristics of some genes, the extent to which protein isoforms globally contribute to functional complexity on a proteomic scale remains unknown. To address this systematically, we cloned full-length open reading frames of alternatively spliced transcripts for a large number of human genes and used protein-protein interaction profiling to functionally compare hundreds of protein isoform pairs. The majority of isoform pairs share less than 50% of their interactions. In the global context of interactome network maps, alternative isoforms tend to behave like distinct proteins rather than minor variants of each other. Interaction partners specific to alternative isoforms tend to be expressed in a highly tissue-specific manner and belong to distinct functional modules. Our strategy, applicable to other functional characteristics, reveals a widespread expansion of protein interaction capabilities through alternative splicing and suggests that many alternative "isoforms" are functionally divergent (i.e., "functional alloforms").

A new view of transcriptome complexity and regulation through the lens of local splicing variations

Alternative splicing (AS) can critically affect gene function and disease, yet mapping splicing variations remains a challenge. Here, we propose a new approach to define and quantify mRNA splicing in units of local splicing variations (LSVs). LSVs capture previously defined types of alternative splicing as well as more complex transcript variations. Building the first genome wide map of LSVs from twelve mouse tissues, we find complex LSVs constitute over 30% of tissue dependent transcript variations and affect specific protein families. We show the prevalence of complex LSVs is conserved in humans and identify hundreds of LSVs that are specific to brain subregions or altered in Alzheimer's patients. Amongst those are novel isoforms in the Camk2 family and a novel poison exon in Ptbp1, a key splice factor in neurogenesis. We anticipate the approach presented here will advance the ability to relate tissue-specific splice variation to genetic variation, phenotype, and disease.

A proteogenomic approach to understand splice isoform functions through sequence and expression-based computational modeling

The products of multi-exon genes are a mixture of alternatively spliced isoforms, from which the translated proteins can have similar, different or even opposing functions. It is therefore essential to differentiate and annotate functions for individual isoforms. Computational approaches provide an efficient complement to expensive and time-consuming experimental studies. The input data of these methods range from DNA sequence, to RNA selection pressure, to expressed sequence tags, to full-length complementary DNA, to exon array, to RNA-seq expression, to proteomic data. Notably, RNA-seq technology generates quantitative profiling of transcript expression at the genome scale, with an unprecedented amount of expression data available for developing isoform function prediction methods. Integrative analysis of these data at different molecular levels enables a proteogenomic approach to systematically interrogate isoform functions. Here, we briefly review the state-of-the-art methods according to their input data sources, discuss their advantages and limitations and point out potential ways to improve prediction accuracies.

Molecular recognition features (MoRFs) in three domains of life

MoRFs are widespread intrinsically disordered protein-binding regions that have similar abundance and amino acid composition across the three domains of life.

Roles of RNA-binding Proteins and Post-transcriptional Regulation in Driving Male Germ Cell Development in the Mouse

Tissue development and homeostasis are dependent on highly regulated gene expression programs in which cell-specific combinations of regulatory factors determine which genes are expressed and the post-transcriptional fate of the resulting RNA transcripts. Post-transcriptional regulation of gene expression by RNA-binding proteins has critical roles in tissue development-allowing individual genes to generate multiple RNA and protein products, and the timing, location, and abundance of protein synthesis to be finely controlled. Extensive post-transcriptional regulation occurs during mammalian gametogenesis, including high levels of alternative mRNA expression, stage-specific expression of mRNA variants, broad translational repression, and stage-specific activation of mRNA translation. In this chapter, an overview of the roles of RNA-binding proteins and the importance of post-transcriptional regulation in male germ cell development in the mouse is presented.

Aberrant Splicing of Estrogen Receptor, HER2, and CD44 Genes in Breast Cancer

Breast cancer (BC) is the most common cause of cancer-related death among women under the age of 50 years. Established biomarkers, such as hormone receptors (estrogen receptor [ER]/progesterone receptor) and human epidermal growth factor receptor 2 (HER2), play significant roles in the selection of patients for endocrine and trastuzumab therapies. However, the initial treatment response is often followed by tumor relapse with intrinsic resistance to the first-line therapy, so it has been expected to identify novel molecular markers to improve the survival and quality of life of patients. Alternative splicing of pre-messenger RNAs is a ubiquitous and flexible mechanism for the control of gene expression in mammalian cells. It provides cells with the opportunity to create protein isoforms with different, even opposing, functions from a single genomic locus. Aberrant alternative splicing is very common in cancer where emerging tumor cells take advantage of this flexibility to produce proteins that promote cell growth and survival. While a number of splicing alterations have been reported in human cancers, we focus on aberrant splicing of ER, HER2, and CD44 genes from the viewpoint of BC development. ERα36, a splice variant from the ER1 locus, governs nongenomic membrane signaling pathways triggered by estrogen and confers 4-hydroxytamoxifen resistance in BC therapy. The alternative spliced isoform of HER2 lacking exon 20 (Δ16HER2) has been reported in human BC; this isoform is associated with transforming ability than the wild-type HER2 and recapitulates the phenotypes of endocrine therapy-resistant BC. Although both CD44 splice isoforms (CD44s, CD44v) play essential roles in BC development, CD44v is more associated with those with favorable prognosis, such as luminal A subtype, while CD44s is linked to those with poor prognosis, such as HER2 or basal cell subtypes that are often metastatic. Hence, the detection of splice variants from these loci will provide keys to understand the pathogenesis, predict the prognosis, and choose specific therapies for BC.

Reads meet rotamers: structural biology in the age of deep sequencing

Structure has traditionally been interrelated with sequence, usually in the framework of comparing sequences across species sharing a common fold. However, the nature of information within the sequence and structure databases is evolving, changing the type of comparisons possible. In particular, we now have a vast amount of personal genome sequences from human populations and a greater fraction of new structures contain interacting proteins within large complexes. Consequently, we have to recast our conception of sequence conservation and its relation to structure-for example, focusing more on selection within the human population. Moreover, within structural biology there is less emphasis on the discovery of novel folds and more on relating structures to networks of protein interactions. We cover this changing mindset here.

The effect of a muscarinic receptor 1 gene variant on grey matter volume in schizophrenia

Previous research has demonstrated that individuals with schizophrenia who are homozygous at the c.267C>A single nucleotide polymorphism (rs2067477) within the cholinergic muscarinic M1 receptor (CHRM1) perform less well on the Wisconsin Card Sorting Test (WCST) than those who are heterozygous. This study sought to determine whether variation in the rs2067477 genotype was associated with differential changes in brain structure. Data from 227 patients with established schizophrenia or schizoaffective disorder were obtained from the Australian Schizophrenia Research Bank. Whole-brain voxel-based morphometry was performed to compare regional grey matter volume (GMV) between the 267C/C (N=191) and 267C/A (N=36) groups. Secondary analyses tested for an effect of genotype on cognition (the WCST was not available). Individuals who were homozygous (267C/C) demonstrated significantly reduced GMV in the right precentral gyrus compared to those who were heterozygous (267C/A). These preliminary results suggest that the rs2067477 genotype is associated with brain structure in the right precentral gyrus in individuals with schizophrenia/schizoaffective disorder. Future studies are required to replicate these results and directly link the volumetric reductions with specific cognitive processes.

Experimental detection of short regulatory motifs in eukaryotic proteins: tips for good practice as well as for bad

It has become clear in outline though not yet in detail how cellular regulatory and signalling systems are constructed. The essential machines are protein complexes that effect regulatory decisions by undergoing internal changes of state. Subcomponents of these cellular complexes are assembled into molecular switches. Many of these switches employ one or more short peptide motifs as toggles that can move between one or more sites within the switch system, the simplest being on-off switches. Paradoxically, these motif modules (termed short linear motifs or SLiMs) are both hugely abundant but difficult to research. So despite the many successes in identifying short regulatory protein motifs, it is thought that only the “tip of the iceberg” has been exposed. Experimental and bioinformatic motif discovery remain challenging and error prone. The advice presented in this article is aimed at helping researchers to uncover genuine protein motifs, whilst avoiding the pitfalls that lead to reports of false discovery.

PPIXpress: construction of condition-specific protein interaction networks based on transcript expression

Protein-protein interaction networks are an important component of modern systems biology. Yet, comparatively few efforts have been made to tailor their topology to the actual cellular condition being studied. Here, we present a network construction method that exploits expression data at the transcript-level and thus reveals alterations in protein connectivity not only caused by differential gene expression but also by alternative splicing. We achieved this by establishing a direct correspondence between individual protein interactions and underlying domain interactions in a complete but condition-unspecific protein interaction network. This knowledge was then used to infer the condition-specific presence of interactions from the dominant protein isoforms. When we compared contextualized interaction networks of matched normal and tumor samples in breast cancer, our transcript-based construction identified more significant alterations that affected proteins associated with cancerogenesis than a method that only uses gene expression data. The approach is provided as the user-friendly tool PPIXpress.

AVAILABILITY AND IMPLEMENTATION

PPIXpress is available at https://sourceforge.net/projects/ppixpress/.

RNA Binding Protein Ptbp2 Is Essential for Male Germ Cell Development

RNA binding proteins (RBPs) are increasingly recognized as essential factors in tissue development and homeostasis. The polypyrimidine tract binding (PTB) protein family of RBPs are important posttranscriptional regulators of gene expression. In the nervous system, the function and importance of PTB protein 2 (Ptbp2) as a key alternative splicing regulator is well established. Ptbp2 is also abundantly expressed during spermatogenesis, but its role in this developmental program has not been explored. Additionally, the importance of alternative splicing regulation in spermatogenesis is unclear. Here, we demonstrate that Ptbp2 is essential for spermatogenesis. We also describe an improved dual fluorescence flow cytometry strategy to discriminate, quantify, and collect germ cells in different stages of development. Using this approach, in combination with traditional histological methods, we show that Ptbp2 ablation results in germ cell loss due to increased apoptosis of meiotic spermatocytes and postmeiotic arrest of spermatid differentiation. Furthermore, we show that Ptbp2 is required for alternative splicing regulation in the testis, as in brain. Strikingly, not all of the alternatively spliced RNAs examined were sensitive to Ptbp2 loss in both tissues. Collectively, the data provide evidence for an important role for alternative splicing regulation in germ cell development and a central role for Ptbp2 in this process.

Human protein interaction networks across tissues and diseases

Protein interaction networks are an important framework for studying protein function, cellular processes, and genotype-to-phenotype relationships. While our view of the human interaction network is constantly expanding, less is known about networks that form in biologically important contexts such as within distinct tissues or in disease conditions. Here we review efforts to characterize these networks and to harness them to gain insights into the molecular mechanisms underlying human disease.

DULIP: A Dual Luminescence-Based Co-Immunoprecipitation Assay for Interactome Mapping in Mammalian Cells

Mapping of protein-protein interactions (PPIs) is critical for understanding protein function and complex biological processes. Here, we present DULIP, a dual luminescence-based co-immunoprecipitation assay, for systematic PPI mapping in mammalian cells. DULIP is a second-generation luminescence-based PPI screening method for the systematic and quantitative analysis of co-immunoprecipitations using two different luciferase tags. Benchmarking studies with positive and negative PPI reference sets revealed that DULIP allows the detection of interactions with high sensitivity and specificity. Furthermore, the analysis of a PPI reference set with known binding affinities demonstrated that both low- and high-affinity interactions can be detected with DULIP assays. Finally, using the well-characterized interaction between Syntaxin-1 and Munc18, we found that DULIP is capable of detecting the effects of point mutations on interaction strength. Taken together, our studies demonstrate that DULIP is a sensitive and reliable method of great utility for systematic interactome research. It can be applied for interaction screening and validation of PPIs in mammalian cells. Moreover, DULIP permits the specific analysis of mutation-dependent binding patterns.

Systems biology and gene networks in neurodevelopmental and neurodegenerative disorders

Genetic and genomic approaches have implicated hundreds of genetic loci in neurodevelopmental disorders and neurodegeneration, but mechanistic understanding continues to lag behind the pace of gene discovery. Understanding the role of specific genetic variants in the brain involves dissecting a functional hierarchy that encompasses molecular pathways, diverse cell types, neural circuits and, ultimately, cognition and behaviour. With a focus on transcriptomics, this Review discusses how high-throughput molecular, integrative and network approaches inform disease biology by placing human genetics in a molecular systems and neurobiological context. We provide a framework for interpreting network biology studies and leveraging big genomics data sets in neurobiology.

Pushing the annotation of cellular activities to a higher resolution: Predicting functions at the isoform level

In past decades, the experimental determination of protein functions was expensive and time-consuming, so numerous computational methods were developed to speed up and guide the process. However, most of these methods predict protein functions at the gene level and do not consider the fact that protein isoforms (translated from alternatively spliced transcripts), not genes, are the actual function carriers. Now, high-throughput RNA-seq technology is providing unprecedented opportunities to unravel protein functions at the isoform level. In this article, we review recent progress in the high-resolution functional annotations of protein isoforms, focusing on two methods developed by the authors. Both methods can integrate multiple RNA-seq datasets for comprehensively characterizing functions of protein isoforms.

Differential and coherent processing patterns from small RNAs

Post-transcriptional processing events related to short RNAs are often reflected in their read profile patterns emerging from high-throughput sequencing data. MicroRNA arm switching across different tissues is a well-known example of what we define as differential processing. Here, short RNAs from the nine cell lines of the ENCODE project, irrespective of their annotation status, were analyzed for genomic loci representing differential or coherent processing. We observed differential processing predominantly in RNAs annotated as miRNA, snoRNA or tRNA. Four out of five known cases of differentially processed miRNAs that were in the input dataset were recovered and several novel cases were discovered. In contrast to differential processing, coherent processing is observed widespread in both annotated and unannotated regions. While the annotated loci predominantly consist of ~24 nt short RNAs, the unannotated loci comparatively consist of ~17 nt short RNAs. Furthermore, these ~17 nt short RNAs are significantly enriched for overlap to transcription start sites and DNase I hypersensitive sites (p-value < 0.01) that are characteristic features of transcription initiation RNAs. We discuss how the computational pipeline developed in this study has the potential to be applied to other forms of RNA-seq data for further transcriptome-wide studies of differential and coherent processing.

How do disordered regions achieve comparable functions to structured domains?

The traditional structure to function paradigm conceives of a protein's function as emerging from its structure. In recent years, it has been established that unstructured, intrinsically disordered regions (IDRs) in proteins are equally crucial elements for protein function, regulation and homeostasis. In this review, we provide a brief overview of how IDRs can perform similar functions to structured proteins, focusing especially on the formation of protein complexes and assemblies and the mediation of regulated conformational changes. In addition to highlighting instances of such functional equivalence, we explain how differences in the biological and physicochemical properties of IDRs allow them to expand the functional and regulatory repertoire of proteins. We also discuss studies that provide insights into how mutations within functional regions of IDRs can lead to human diseases.

My road to alternative splicing control: from simple paths to loops and interconnections

With the functional importance of alternative splicing being validated in nearly every mammalian biological system and implicated in many human diseases, it is now crucial to identify the molecular programs that control the production of splice variants. In this article, I will survey how our knowledge of the basic principles of alternative splicing control evolved over the last 25 years. I will also describe how investigation of the splicing control of an apoptotic regulator led us to identify novel effectors and revealed the existence of converging pathways linking splicing decisions to DNA damage. Finally, I will review how our efforts at developing tools designed to monitor and redirect splicing helped assess the impact of misregulated splicing in cancer.

Revealing Pathway Dynamics in Heart Diseases by Analyzing Multiple Differential Networks

Development of heart diseases is driven by dynamic changes in both the activity and connectivity of gene pathways. Understanding these dynamic events is critical for understanding pathogenic mechanisms and development of effective treatment. Currently, there is a lack of computational methods that enable analysis of multiple gene networks, each of which exhibits differential activity compared to the network of the baseline/healthy condition. We describe the iMDM algorithm to identify both unique and shared gene modules across multiple differential co-expression networks, termed M-DMs (multiple differential modules). We applied iMDM to a time-course RNA-Seq dataset generated using a murine heart failure model generated on two genotypes. We showed that iMDM achieves higher accuracy in inferring gene modules compared to using single or multiple co-expression networks. We found that condition-specific M-DMs exhibit differential activities, mediate different biological processes, and are enriched for genes with known cardiovascular phenotypes. By analyzing M-DMs that are present in multiple conditions, we revealed dynamic changes in pathway activity and connectivity across heart failure conditions. We further showed that module dynamics were correlated with the dynamics of disease phenotypes during the development of heart failure. Thus, pathway dynamics is a powerful measure for understanding pathogenesis. iMDM provides a principled way to dissect the dynamics of gene pathways and its relationship to the dynamics of disease phenotype. With the exponential growth of omics data, our method can aid in generating systems-level insights into disease progression.

Essential roles for the splicing regulator nSR100/SRRM4 during nervous system development

Quesnel-Vallières et al. show that loss of the vertebrate- and neural-specific Ser/Arg repeat-related protein of 100 kDa (nSR100/SRRM4) impairs development of the central and peripheral nervous systems. Accompanying these developmental defects are widespread changes in alternative splicing (AS) that primarily result in shifts to nonneural patterns for different classes of splicing events. The main component of the altered AS program comprises 3- to 27-nt neural microexons, and inclusion of a 6-nt nSR100-activated microexon in Unc13b transcripts is sufficient to rescue a neuritogenesis defect in nSR100 mutant primary neurons.

Alternative splicing (AS) generates vast transcriptomic complexity in the vertebrate nervous system. However, the extent to which trans-acting splicing regulators and their target AS regulatory networks contribute to nervous system development is not well understood. To address these questions, we generated mice lacking the vertebrate- and neural-specific Ser/Arg repeat-related protein of 100 kDa (nSR100/SRRM4). Loss of nSR100 impairs development of the central and peripheral nervous systems in part by disrupting neurite outgrowth, cortical layering in the forebrain, and axon guidance in the corpus callosum. Accompanying these developmental defects are widespread changes in AS that primarily result in shifts to nonneural patterns for different classes of splicing events. The main component of the altered AS program comprises 3- to 27-nucleotide (nt) neural microexons, an emerging class of highly conserved AS events associated with the regulation of protein interaction networks in developing neurons and neurological disorders. Remarkably, inclusion of a 6-nt, nSR100-activated microexon in Unc13b transcripts is sufficient to rescue a neuritogenesis defect in nSR100 mutant primary neurons. These results thus reveal critical in vivo neurodevelopmental functions of nSR100 and further link these functions to a conserved program of neuronal microexon splicing.

Role of the Wnt/β-catenin pathway in gastric cancer: An in-depth literature review

Gastric cancer remains one of the most common cancers worldwide and one of the leading cause for cancer-related deaths. Gastric adenocarcinoma is a multifactorial disease that is genetically, cytologically and architecturally more heterogeneous than other gastrointestinal carcinomas. The aberrant activation of the Wnt/β-catenin signaling pathway is involved in the development and progression of a significant proportion of gastric cancer cases. This review focuses on the participation of the Wnt/β-catenin pathway in gastric cancer by offering an analysis of the relevant literature published in this field. Indeed, it is discussed the role of key factors in Wnt/β-catenin signaling and their downstream effectors regulating processes involved in tumor initiation, tumor growth, metastasis and resistance to therapy. Available data indicate that constitutive Wnt signalling resulting from Helicobacter pylori infection and inactivation of Wnt inhibitors (mainly by inactivating mutations and promoter hypermethylation) play an important role in gastric cancer. Moreover, a number of recent studies confirmed CTNNB1 and APC as driver genes in gastric cancer. The identification of specific membrane, intracellular, and extracellular components of the Wnt pathway has revealed potential targets for gastric cancer therapy. High-throughput “omics” approaches will help in the search for Wnt pathway antagonist in the near future.

Proteomic Validation of Transcript Isoforms, Including Those Assembled from RNA-Seq Data

Human proteome analysis now requires an understanding of protein isoforms. We recently published the PG Nexus pipeline, which facilitates high confidence validation of exons and splice junctions by integrating genomics and proteomics data. Here we comprehensively explore how RNA-seq transcriptomics data, and proteomic analysis of the same sample, can identify protein isoforms. RNA-seq data from human mesenchymal (hMSC) stem cells were analyzed with our new TranscriptCoder tool to generate a database of protein isoform sequences. MS/MS data from matching hMSC samples were then matched against the TranscriptCoder-derived database, along with Ensembl and the neXtProt database. Querying the TranscriptCoder-derived or Ensembl database could unambiguously identify ∼450 protein isoforms, with isoform-specific proteotypic peptides, including candidate hMSC-specific isoforms for the genes DPYSL2 and FXR1. Where isoform-specific peptides did not exist, groups of nonisoform-specific proteotypic peptides could specifically identify many isoforms. In both the above cases, isoforms will be detectable with targeted MS/MS assays. Unfortunately, our analysis also revealed that some isoforms will be difficult to identify unambiguously as they do not have peptides that are sufficiently distinguishing. We covisualize mRNA isoforms and peptides in a genome browser to illustrate the above situations. Mass spectrometry data is available via ProteomeXchange (PXD001449).

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.

The variation game: Cracking complex genetic disorders with NGS and omics data

Tremendous advances in Next Generation Sequencing (NGS) and high-throughput omics methods have brought us one step closer towards mechanistic understanding of the complex disease at the molecular level. In this review, we discuss four basic regulatory mechanisms implicated in complex genetic diseases, such as cancer, neurological disorders, heart disease, diabetes, and many others. The mechanisms, including genetic variations, copy-number variations, posttranscriptional variations, and epigenetic variations, can be detected using a variety of NGS methods. We propose that malfunctions detected in these mechanisms are not necessarily independent, since these malfunctions are often found associated with the same disease and targeting the same gene, group of genes, or functional pathway. As an example, we discuss possible rewiring effects of the cancer-associated genetic, structural, and posttranscriptional variations on the protein-protein interaction (PPI) network centered around P53 protein. The review highlights multi-layered complexity of common genetic disorders and suggests that integration of NGS and omics data is a critical step in developing new computational methods capable of deciphering this complexity.

Intrinsically disordered proteins in cellular signalling and regulation

Intrinsically disordered proteins (IDPs) are important components of the cellular signalling machinery, allowing the same polypeptide to undertake different interactions with different consequences. IDPs are subject to combinatorial post-translational modifications and alternative splicing, adding complexity to regulatory networks and providing a mechanism for tissue-specific signalling. These proteins participate in the assembly of signalling complexes and in the dynamic self-assembly of membrane-less nuclear and cytoplasmic organelles. Experimental, computational and bioinformatic analyses combine to identify and characterize disordered regions of proteins, leading to a greater appreciation of their widespread roles in biological processes.

Oncogenic fusion protein EWS-FLI1 is a network hub that regulates alternative splicing

The synthesis and processing of mRNA, from transcription to translation initiation, often requires splicing of intragenic material. The final mRNA composition varies based on proteins that modulate splice site selection. EWS-FLI1 is an Ewing sarcoma (ES) oncoprotein with an interactome that we demonstrate to have multiple partners in spliceosomal complexes. We evaluate the effect of EWS-FLI1 on posttranscriptional gene regulation using both exon array and RNA-seq. Genes that potentially regulate oncogenesis, including CLK1, CASP3, PPFIBP1, and TERT, validate as alternatively spliced by EWS-FLI1. In a CLIP-seq experiment, we find that EWS-FLI1 RNA-binding motifs most frequently occur adjacent to intron-exon boundaries. EWS-FLI1 also alters splicing by directly binding to known splicing factors including DDX5, hnRNP K, and PRPF6. Reduction of EWS-FLI1 produces an isoform of γ-TERT that has increased telomerase activity compared with wild-type (WT) TERT. The small molecule YK-4-279 is an inhibitor of EWS-FLI1 oncogenic function that disrupts specific protein interactions, including helicases DDX5 and RNA helicase A (RHA) that alters RNA-splicing ratios. As such, YK-4-279 validates the splicing mechanism of EWS-FLI1, showing alternatively spliced gene patterns that significantly overlap with EWS-FLI1 reduction and WT human mesenchymal stem cells (hMSC). Exon array analysis of 75 ES patient samples shows similar isoform expression patterns to cell line models expressing EWS-FLI1, supporting the clinical relevance of our findings. These experiments establish systemic alternative splicing as an oncogenic process modulated by EWS-FLI1. EWS-FLI1 modulation of mRNA splicing may provide insight into the contribution of splicing toward oncogenesis, and, reciprocally, EWS-FLI1 interactions with splicing proteins may inform the splicing code.

Rethinking gene regulatory networks in light of alternative splicing, intrinsically disordered protein domains, and post-translational modifications

Models for genetic regulation and cell fate specification characteristically assume that gene regulatory networks (GRNs) are essentially deterministic and exhibit multiple stable states specifying alternative, but pre-figured cell fates. Mounting evidence shows, however, that most eukaryotic precursor RNAs undergo alternative splicing (AS) and that the majority of transcription factors contain intrinsically disordered protein (IDP) domains whose functionalities are context dependent as well as subject to post-translational modification (PTM). Consequently, many transcription factors do not have fixed cis-acting regulatory targets, and developmental determination by GRNs alone is untenable. Modeling these phenomena requires a multi-scale approach to explain how GRNs operationally interact with the intra- and intercellular environments. Evidence shows that AS, IDP, and PTM complicate gene expression and act synergistically to facilitate and promote time- and cell-specific protein modifications involved in cell signaling and cell fate specification and thereby disrupt a strict deterministic GRN-phenotype mapping. The combined effects of AS, IDP, and PTM give proteomes physiological plasticity, adaptive responsiveness, and developmental versatility without inefficiently expanding genome size. They also help us understand how protein functionalities can undergo major evolutionary changes by buffering mutational consequences.

A highly conserved program of neuronal microexons is misregulated in autistic brains

Alternative splicing (AS) generates vast transcriptomic and proteomic complexity. However, which of the myriad of detected AS events provide important biological functions is not well understood. Here, we define the largest program of functionally coordinated, neural-regulated AS described to date in mammals. Relative to all other types of AS within this program, 3-15 nucleotide "microexons" display the most striking evolutionary conservation and switch-like regulation. These microexons modulate the function of interaction domains of proteins involved in neurogenesis. Most neural microexons are regulated by the neuronal-specific splicing factor nSR100/SRRM4, through its binding to adjacent intronic enhancer motifs. Neural microexons are frequently misregulated in the brains of individuals with autism spectrum disorder, and this misregulation is associated with reduced levels of nSR100. The results thus reveal a highly conserved program of dynamic microexon regulation associated with the remodeling of protein-interaction networks during neurogenesis, the misregulation of which is linked to autism.

IIIDB: a database for isoform-isoform interactions and isoform network modules

BACKGROUND

Protein-protein interactions (PPIs) are key to understanding diverse cellular processes and disease mechanisms. However, current PPI databases only provide low-resolution knowledge of PPIs, in the sense that "proteins" of currently known PPIs generally refer to "genes." It is known that alternative splicing often impacts PPI by either directly affecting protein interacting domains, or by indirectly impacting other domains, which, in turn, impacts the PPI binding. Thus, proteins translated from different isoforms of the same gene can have different interaction partners.

RESULTS

Due to the limitations of current experimental capacities, little data is available for PPIs at the resolution of isoforms, although such high-resolution data is crucial to map pathways and to understand protein functions. In fact, alternative splicing can often change the internal structure of a pathway by rearranging specific PPIs. To fill the gap, we systematically predicted genome-wide isoform-isoform interactions (IIIs) using RNA-seq datasets, domain-domain interaction and PPIs. Furthermore, we constructed an III database (IIIDB) that is a resource for studying PPIs at isoform resolution. To discover functional modules in the III network, we performed III network clustering, and then obtained 1025 isoform modules. To evaluate the module functionality, we performed the GO/pathway enrichment analysis for each isoform module.

CONCLUSIONS

The IIIDB provides predictions of human protein-protein interactions at the high resolution of transcript isoforms that can facilitate detailed understanding of protein functions and biological pathways. The web interface allows users to search for IIIs or III network modules. The IIIDB is freely available at http://syslab.nchu.edu.tw/IIIDB.

RBM24 is a major regulator of muscle-specific alternative splicing

Cell-type-specific splicing generates numerous alternatively spliced transcripts playing important roles for organ development and homeostasis, but only a few tissue-specific splicing factors have been identified. We found that RBM24 governs a large number of muscle-specific splicing events that are critically involved in cardiac and skeletal muscle development and disease. Targeted inactivation of RBM24 in mice disrupted cardiac development and impaired sarcomerogenesis in striated muscles. In vitro splicing assays revealed that recombinant RBM24 is sufficient to promote muscle-specific exon inclusion in nuclear extracts of nonmuscle cells. Furthermore, we demonstrate that binding of RBM24 to an intronic splicing enhancer (ISE) is essential and sufficient to overcome repression of exon inclusion by an exonic splicing silencer (ESS) containing PTB and hnRNP A1/A2 binding sites. Introduction of ESS and ISE converted a constitutive exon into an RMB24-dependent alternative exon. We reason that RBM24 is a major regulator of alternative splicing in striated muscles.

CLUMSY VEIN, the Arabidopsis DEAH-box Prp16 ortholog, is required for auxin-mediated development

Pre-messenger RNA (pre-mRNA) splicing is essential in eukaryotic cells. In animals and yeasts, the DEAH-box RNA-dependent ATPase Prp16 mediates conformational change of the spliceosome, thereby facilitating pre-mRNA splicing. In yeasts, Prp16 also plays an important role in splicing fidelity. Conversely, PRP16 orthologs in Chlamydomonas reinhardtii and nematode do not have an important role in general pre-mRNA splicing, but are required for gene silencing and sex determination, respectively. Functions of PRP16 orthologs in higher plants have not been described until now. Here we show that the CLUMSY VEIN (CUV) gene encoding the unique Prp16 ortholog in Arabidopsis thaliana facilitates auxin-mediated development including male-gametophyte transmission, apical-basal patterning of embryonic and gynoecium development, stamen development, phyllotactic flower positioning, and vascular development. cuv-1 mutation differentially affects splicing and expression of genes involved in auxin biosynthesis, polar auxin transport, auxin perception and auxin signaling. The cuv-1 mutation does not have an equal influence on pre-mRNA substrates. We propose that Arabidopsis PRP16/CUV differentially facilitates expression of genes, which include genes involved in auxin biosynthesis, transport, perception and signaling, thereby collectively influencing auxin-mediated development.

Alternative RNA structure-coupled gene regulations in tumorigenesis

Alternative RNA structures (ARSs), or alternative transcript isoforms, are critical for regulating cellular phenotypes in humans. In addition to generating functionally diverse protein isoforms from a single gene, ARS can alter the sequence contents of 5'/3' untranslated regions (UTRs) and intronic regions, thus also affecting the regulatory effects of these regions. ARS may introduce premature stop codon(s) into a transcript, and render the transcript susceptible to nonsense-mediated decay, which in turn can influence the overall gene expression level. Meanwhile, ARS can regulate the presence/absence of upstream open reading frames and microRNA targeting sites in 5'UTRs and 3'UTRs, respectively, thus affecting translational efficiencies and protein expression levels. Furthermore, since ARS may alter exon-intron structures, it can influence the biogenesis of intronic microRNAs and indirectly affect the expression of the target genes of these microRNAs. The connections between ARS and multiple regulatory mechanisms underline the importance of ARS in determining cell fate. Accumulating evidence indicates that ARS-coupled regulations play important roles in tumorigenesis. Here I will review our current knowledge in this field, and discuss potential future directions.

Revisiting the identification of canonical splice isoforms through integration of functional genomics and proteomics evidence

Canonical isoforms in different databases have been defined as the most prevalent, most conserved, most expressed, longest, or the one with the clearest description of domains or posttranslational modifications. In this article, we revisit these definitions of canonical isoforms based on functional genomics and proteomics evidence, focusing on mouse data. We report a novel functional relationship network-based approach for identifying the highest connected isoforms (HCIs). We show that 46% of these HCIs are not the longest transcripts. In addition, this approach revealed many genes that have more than one highly connected isoforms. Averaged across 175 RNA-seq datasets covering diverse tissues and conditions, 65% of the HCIs show higher expression levels than nonhighest connected isoforms at the transcript level. At the protein level, these HCIs highly overlap with the expressed splice variants, based on proteomic data from eight different normal tissues. These results suggest that a more confident definition of canonical isoforms can be made through integration of multiple lines of evidence, including HCIs defined by biological processes and pathways, expression prevalence at the transcript level, and relative or absolute abundance at the protein level. This integrative proteogenomics approach can successfully identify principal isoforms that are responsible for the canonical functions of genes.

BIN1 membrane curvature sensing and generation show autoinhibition regulated by downstream ligands and PI(4,5)P2

In striated muscles, invaginations from the plasma membrane, termed transverse tubules (T-tubule), function in the excitation–contraction coupling machinery. BIN1 (isoform8) plays a critical role in the biogenesis of T-tubules. BIN1 contains an N-terminal BAR domain to sense and induce membrane curvature, an isoform8-specific polybasic motif (exon10) as the phosphoinositide binding module and a C-terminal Src homology 3 (SH3) domain for the recruitment of downstream proteins such as dynamin 2. Previous studies of N-BAR domains focused on elucidating mechanisms of membrane curvature sensing and generation (MC-S&G). Less is known about how MC-S&G is regulated. We found that the SH3 domain binds to the exon10 motif more strongly compared to the proline-rich domain (PRD) of dynamin 2. Furthermore, we found that the MC-S&G ability of full-length BIN1 is inhibited on membranes lacking PI(4,5)P₂. Addition of PI(4,5)P₂ in the membrane activates BIN1 to sense and induce membrane curvature. Co-presence of the SH3 domain and exon10 motif leads to the strongest phosphoinositide-mediated control of BIN1 function. Addition of SH3 domain ligand (such as PRD peptides), as well as addition of the water-soluble PI(4,5)P₂ analogue, can both enhance the MC-S&G ability of BIN1 on membranes without PI(4,5)P₂, indicating that the key to activate BIN1 is to disrupt the exon10–SH3 interaction. The nonsense mutation K436X, found in centronuclear myopathy (CNM) patients, abolishes SH3 domain binding with either exon10 or the PRD motif, resulting in increased membrane deformation capacity. Our results suggest an autoinhibition model for BIN1 that involves a synergistic regulation by membrane composition and protein–protein interactions.

A million peptide motifs for the molecular biologist

A molecular description of functional modules in the cell is the focus of many high-throughput studies in the postgenomic era. A large portion of biomolecular interactions in virtually all cellular processes is mediated by compact interaction modules, referred to as peptide motifs. Such motifs are typically less than ten residues in length, occur within intrinsically disordered regions, and are recognized and/or posttranslationally modified by structured domains of the interacting partner. In this review, we suggest that there might be over a million instances of peptide motifs in the human proteome. While this staggering number suggests that peptide motifs are numerous and the most understudied functional module in the cell, it also holds great opportunities for new discoveries.

The intrinsically disordered regions of the Drosophila melanogaster Hox protein ultrabithorax select interacting proteins based on partner topology

Interactions between structured proteins require a complementary topology and surface chemistry to form sufficient contacts for stable binding. However, approximately one third of protein interactions are estimated to involve intrinsically disordered regions of proteins. The dynamic nature of disordered regions before and, in some cases, after binding calls into question the role of partner topology in forming protein interactions. To understand how intrinsically disordered proteins identify the correct interacting partner proteins, we evaluated interactions formed by the Drosophila melanogaster Hox transcription factor Ultrabithorax (Ubx), which contains both structured and disordered regions. Ubx binding proteins are enriched in specific folds: 23 of its 39 partners include one of 7 folds, out of the 1195 folds recognized by SCOP. For the proteins harboring the two most populated folds, DNA-RNA binding 3-helical bundles and α-α superhelices, the regions of the partner proteins that exhibit these preferred folds are sufficient for Ubx binding. Three disorder-containing regions in Ubx are required to bind these partners. These regions are either alternatively spliced or multiply phosphorylated, providing a mechanism for cellular processes to regulate Ubx-partner interactions. Indeed, partner topology correlates with the ability of individual partner proteins to bind Ubx spliceoforms. Partners bind different disordered regions within Ubx to varying extents, creating the potential for competition between partners and cooperative binding by partners. The ability of partners to bind regions of Ubx that activate transcription and regulate DNA binding provides a mechanism for partners to modulate transcription regulation by Ubx, and suggests that one role of disorder in Ubx is to coordinate multiple molecular functions in response to tissue-specific cues.

Context-dependent control of alternative splicing by RNA-binding proteins

Sequence-specific RNA-binding proteins (RBPs) bind to pre-mRNA to control alternative splicing, but it is not yet possible to read the 'splicing code' that dictates splicing regulation on the basis of genome sequence. Each alternative splicing event is controlled by multiple RBPs, the combined action of which creates a distribution of alternatively spliced products in a given cell type. As each cell type expresses a distinct array of RBPs, the interpretation of regulatory information on a given RNA target is exceedingly dependent on the cell type. RBPs also control each other's functions at many levels, including by mutual modulation of their binding activities on specific regulatory RNA elements. In this Review, we describe some of the emerging rules that govern the highly context-dependent and combinatorial nature of alternative splicing regulation.

A global regulatory mechanism for activating an exon network required for neurogenesis

The vertebrate and neural-specific Ser/Arg (SR)-related protein nSR100/SRRM4 regulates an extensive program of alternative splicing with critical roles in nervous system development. However, the mechanism by which nSR100 controls its target exons is poorly understood. We demonstrate that nSR100-dependent neural exons are associated with a unique configuration of intronic cis-elements that promote rapid switch-like regulation during neurogenesis. A key feature of this configuration is the insertion of specialized intronic enhancers between polypyrimidine tracts and acceptor sites that bind nSR100 to potently activate exon inclusion in neural cells while weakening 3' splice site recognition and contributing to exon skipping in nonneural cells. nSR100 further operates by forming multiple interactions with early spliceosome components bound proximal to 3' splice sites. These multifaceted interactions achieve dominance over neural exon silencing mediated by the splicing regulator PTBP1. The results thus illuminate a widespread mechanism by which a critical neural exon network is activated during neurogenesis.

Nodes occupying central positions in human tissue specific PPI networks are enriched with many splice variants

The functional repertoire of genes in the eukaryotic organisms is enhanced by the phenomenon of alternative splicing. Hence, a node in a tissue specific protein-protein interaction (TS PPIN) network can be thought of as an ensemble of various spliced protein products of the corresponding gene expressed in that tissue. Here we demonstrate that the nodes that occupy topologically central positions characterized by high degree, betweenness, closeness, and eigenvector centrality values in TS PPINs of Homo sapiens are associated with high number of splice variants. We also show that the high "centrality" of these genes/nodes could in part be explained by the presence of a large number of promiscuous domains.

Structured and disordered facets of the GPCR fold

The seven-transmembrane (7TM) helix fold of G-protein coupled receptors (GPCRs) has been adapted for a wide variety of physiologically important signaling functions. Here, we discuss the diversity in the structured and disordered regions of GPCRs based on the recently published crystal structures and sequence analysis of all human GPCRs. A comparison of the structures of rhodopsin-like receptors (class A), secretin-like receptors (class B), metabotropic receptors (class C) and frizzled receptors (class F) shows that the relative arrangement of the transmembrane helices is conserved across all four GPCR classes although individual receptors can be activated by ligand binding at varying positions within and around the transmembrane helical bundle. A systematic analysis of GPCR sequences reveals the presence of disordered segments in the cytoplasmic side, abundant post-translational modification sites, evidence for alternative splicing and several putative linear peptide motifs that have the potential to mediate interactions with cytosolic proteins. While the structured regions permit the receptor to bind diverse ligands, the disordered regions appear to have an underappreciated role in modulating downstream signaling in response to the cellular state. An integrated paradigm combining the knowledge of structured and disordered regions is imperative for gaining a holistic understanding of the GPCR (un)structure-function relationship.

Intracellular clusterin interacts with brain isoforms of the bridging integrator 1 and with the microtubule-associated protein Tau in Alzheimer's disease

Sporadic or late-onset Alzheimer's disease (AD) is expected to affect 50% of individuals reaching 85 years of age. The most significant genetic risk factor for late-onset AD is the e4 allele of APOE gene encoding apolipoprotein E, a lipid carrier shown to modulate brain amyloid burden. Recent genome-wide association studies have uncovered additional single nucleotide polymorphisms (SNPs) linked to AD susceptibility, including those in the CLU and BIN1 genes encoding for clusterin (CLU) and the bridging integrator 1 (BIN1) proteins, respectively. Because CLU has been implicated in brain amyloid-β (Aβ) clearance in mouse models of amyloid deposition, we sought to investigate whether an AD-linked SNP in the CLU gene altered Aβ42 biomarker levels in the cerebrospinal fluid (CSF). Instead, we found that the CLU rs11136000 SNP modified CSF levels of the microtubule-associated protein Tau in AD patients. We also found that an intracellular form of CLU (iCLU) was upregulated in the brain of Tau overexpressing Tg4510 mice, but not in Tg2576 amyloid mouse model. By overexpressing iCLU and Tau in cell culture systems we discovered that iCLU was a Tau-interacting protein and that iCLU associated with brain-specific isoforms of BIN1, also recently identified as a Tau-binding protein. Through expression analysis of CLU and BIN1 variants, we found that CLU and BIN1 interacted via their coiled-coil motifs. In co-immunoprecipitation studies using human brain tissue, we showed that iCLU and the major BIN1 isoform expressed in neurons were associated with modified Tau species found in AD. Finally, we showed that expression of certain coding CLU variants linked to AD risk led to increased levels of iCLU. Together, our findings suggest that iCLU and BIN1 interaction might impact Tau function in neurons and uncover potential new mechanisms underlying the etiology of Tau pathology in AD.

Alternative splicing of MBD2 supports self-renewal in human pluripotent stem cells

Alternative RNA splicing (AS) regulates proteome diversity, including isoform-specific expression of several pluripotency genes. Here, we integrated global gene expression and proteomic analyses and identified a molecular signature suggesting a central role for AS in maintaining human pluripotent stem cell (hPSC) self-renewal. We demonstrate that the splicing factor SFRS2 is an OCT4 target gene required for pluripotency. SFRS2 regulates AS of the methyl-CpG binding protein MBD2, whose isoforms play opposing roles in maintenance of and reprogramming to pluripotency. Although both MDB2a and MBD2c are enriched at the OCT4 and NANOG promoters, MBD2a preferentially interacts with repressive NuRD chromatin remodeling factors and promotes hPSC differentiation, whereas overexpression of MBD2c enhances reprogramming of fibroblasts to pluripotency. The miR-301 and miR-302 families provide additional regulation by targeting SFRS2 and MDB2a. These data suggest that OCT4, SFRS2, and MBD2 participate in a positive feedback loop, regulating proteome diversity in support of hPSC self-renewal and reprogramming.

The emerging era of genomic data integration for analyzing splice isoform function

The vast majority of multi-exon genes in humans undergo alternative splicing, which greatly increases the functional diversity of protein species. Predicting functions at the isoform level is essential to further our understanding of developmental abnormalities and cancers, which frequently exhibit aberrant splicing and dysregulation of isoform expression. However, determination of isoform function is very difficult, and efforts to predict isoform function have been limited in the functional genomics field. Deep sequencing of RNA now provides an unprecedented amount of expression data at the transcript level. We describe here emerging computational approaches that integrate such large-scale whole-transcriptome sequencing (RNA-seq) data for predicting the functions of alternatively spliced isoforms, and we discuss their applications in developmental and cancer biology. We outline future directions for isoform function prediction, emphasizing the need for heterogeneous genomic data integration and tissue-specific, dynamic isoform-level network modeling, which will allow the field to realize its full potential.