Tissue-specific Alternative Splicing Research Articles

Abstract 46: Identification of Novel Alternate Splicing Events in Humans With RBFOX2 Mutations and Hypoplastic Left Heart Syndrome

Alternative splicing (AS) of protein isoforms is an integral mechanism for cardiac development. RNA Binding Protein, Fox-1 Homolog (C. Elegans) 2 (RBFOX2) is an RNA binding protein preferentially expressed in muscle and neuronal cells and regulates tissue-specific alternate exon splicing in ~2,100 target genes by binding the conserved RNA sequence motif (U)GCAUG. RBFOX2 was recently implicated in the pathogenesis of abnormal cardiac and cerebral development via loss-of-function studies in zebrafish and mouse. However, convincing evidence remains incomplete, as the full complement of RBFOX2 target genes and differential exon usage (DEU) in human cardiovascular cell lines are incompletely defined. We identified de novo mutations in RBFOX2 from four human cases of congenital heart disease (CHD) with hypoplastic left heart syndrome (HLHS) via whole exome sequencing. To test the hypothesis that RBFOX2 mutations alter DEU in known target genes, we performed RNA-seq on ductus arteriosus tissue from human CHD cases with and without RBFOX2 mutation. Analysis of RNA-seq for DEU was performed with DEXSeq. To limit the effect of differential gene expression, we restricted analysis to subjects with high global gene expression correlation (r2 > 0.9, case=1 vs. control=5). DEU in known RBFOX2 target genes were highly enriched compared to all known genes (115/2,100 vs. 589/26,310, p=5.78e-15). A high percentage of the DEU genes (60.0%, 69/115) have high heart expression (HHE) in the developing mouse (mean expression of four cardiac chambers at e14.5). DEU genes with HHE include VCL, TPM1, FN1, ACTN1, and CALD1. Functional annotation clustering reveals enrichment for several actin, cytoskeletal, and contractile Gene Ontology terms, suggesting a possible role in the epithelial-mesenchymal transition (EMT) developmental processes active during early cardiac formation. We also identified enrichment from genes implicated in CHD by allelic specific expression or de novo mutations (44/1,263 vs. 589/26,310, p=0.007). These results are the first in humans to identify differentially expressed exons associated with RBFOX2 mutations in CHD and suggests RBFOX2-mediated alternate splicing may influence EMT pathways implicated in the pathogenesis of HLHS.

Coordinated tissue-specific regulation of adjacent alternative 3' splice sites in C. elegans.

Adjacent alternative 3′ splice sites, those separated by ≤18 nucleotides, provide a unique problem in the study of alternative splicing regulation; there is overlap of the cis-elements that define the adjacent sites. Identification of the intron's 3′ end depends upon sequence elements that define the branchpoint, polypyrimidine tract, and terminal AG dinucleotide. Starting with RNA-seq data from germline-enriched and somatic cell-enriched Caenorhabditis elegans samples, we identify hundreds of introns with adjacent alternative 3′ splice sites. We identify 203 events that undergo tissue-specific alternative splicing. For these, the regulation is monodirectional, with somatic cells preferring to splice at the distal 3′ splice site (furthest from the 5′ end of the intron) and germline cells showing a distinct shift toward usage of the adjacent proximal 3′ splice site (closer to the 5′ end of the intron). Splicing patterns in somatic cells follow C. elegans consensus rules of 3′ splice site definition; a short stretch of pyrimidines preceding an AG dinucleotide. Splicing in germline cells occurs at proximal 3′ splice sites that lack a preceding polypyrimidine tract, and in three instances the germline-specific site lacks the AG dinucleotide. We provide evidence that use of germline-specific proximal 3′ splice sites is conserved across Caenorhabditis species. We propose that there are differences between germline and somatic cells in the way that the basal splicing machinery functions to determine the intron terminus.

Porphobilinogen deaminase gene mutations in Polish patients with non-erythroid acute intermittent porphyria.

Acute intermittent porphyria (AIP) is an metabolic disorder characterized by a partial deficiency of the porphobilinogen deaminase, the enzyme of heme biosynthesis. The metabolic defect in AIP involves an approximately half-normal activity of porphobilinogen deaminase (PBGD, EC 4.3.1.8), the enzyme catalyzing condensation of four porphobilinogen molecules to hydroxymethylbilane. Due to tissue-specific alternative transcript splicing, the PBGD gene mutations within the range of exons 3-15 may lead to classical AIP involving erythrocytes and all the other tissues. Mutations within intron and exon 1 may result in the so-called non-erythroid AIP in which the PBGD activity is normal in erythrocytes and diminished in other tissues. The aim of the present study was to characterise molecular errors in the PBGD gene in Polish patients with non-erythroid AIP and to evaluate the efficacy of the DNA sequencing method in the early diagnosis of this disorder. Twenty five members of nine non-erythroid AIP families were assessed. In each of them DNA sequencing was performed using the Big Dye Terminator Cycle Sequencing Kit v.1.1 on the Hitachi 3730 Analyzer (Applied Biosystem, USA). Four mutations were detected in intron 1 of the PBGD gene, including one unreported novel mutation, 33+(4-12) del AGTGCTGAG, of an unknown biological mechanism, and three previously described mutations, i.e. 33+1 G > A, 33+2 T > C, 33+5 G > C, responsible for abnormal transcript splicing in the area of exon 1. Of 14 asymptomatic members of proband families in 6 subjects were diagnosed with AIP, and in 8 the AIP was excluded based on the DNA sequencing method. DNA sequencing based analysis is the only reliable method for correct diagnosis of asymptomatic non-erythroid AIP patients with normal urinary excretion of heme precursors. The mutations found in Polish patients with non-erythroid AIP represented those of splice defect and resulted in abnormal exon 1 splicing.

Rbfox proteins regulate tissue-specific alternative splicing of Mef2D required for muscle differentiation.

Among the Mef2 family of transcription factors, Mef2D is unique in that it undergoes tissue-specific splicing to generate an isoform that is essential for muscle differentiation. However, the mechanisms mediating this muscle-specific processing of Mef2D remain unknown. Using bioinformatics, we identified Rbfox proteins as putative modulators of Mef2D muscle-specific splicing. Accordingly, we found direct and specific Rbfox1 and Rbfox2 binding to Mef2D pre-mRNA in vivo. Gain- and loss-of-function experiments demonstrated that Rbfox1 and Rbfox2 cooperate in promoting Mef2D splicing and subsequent myogenesis. Thus, our findings reveal a new role for Rbfox proteins in regulating myogenesis through activation of essential muscle-specific splicing events.

Splicing factor SUP-12 and the molecular complexity of apparent cooperativity.

The splicing factor SUP-12 from C. elegans, in combination with either ASD-1 or FOX-1 from the Fox-1 (RBFOX) family, is required for generating a muscle-specific isoform of the fibroblast growth factor receptor EGL-15. Biophysical techniques have revealed the sequence preference for the RNA Recognition Motif (RRM) domain from SUP-12 as well as the structural details of the RNA-bound complex. Detailed genetics have identified a requisite need for the presence of both SUP-12 and ASD-1/FOX-1 to regulate the alternative splicing event, prompting speculation of a cooperative mechanism between these proteins on binding RNA. In contrast, the interplay between SUP-12 and ASD-1 suggests that although the RRM domains from each protein are in direct contact on the egl-15 pre-mRNA, there is no simple contribution of binding cooperativity. Evidence for an independent binding mechanism by SUP-12 and ASD-1 will be discussed, including a model in which both positive and negative contributions are balanced during complex assembly. The ability to monitor tissue-specific alternative splicing in live nematodes will continue to provide a powerful method to test in vivo mechanistic models derived from atomic-level investigation.

Promiscuous mRNA splicing under the control of AIRE in medullary thymic epithelial cells.

The expression of tissue-restricted antigens (TRAs) in the thymus is required to ensure efficient negative selection of potentially auto-reactive T lymphocytes and avoid autoimmune disease. This promiscuous expression is under the control of the autoimmune regulator (AIRE), a transcription factor expressed in medullary thymic epithelial cells (mTECs). Tissue-specific alternative splicing may also produce TRAs but the extent to which splice isoforms that are restricted to specific tissues are expressed in mTECs is yet to be investigated. We reanalyzed microarray and RNA-Seq datasets from mouse mTECs and other epithelial and non-epithelial cell types and found that the diversity of splice isoforms in mTECs was greater than in any of the other cell types or tissues studied. We identified tissue-specific isoforms from a panel of mouse tissues and found several examples of such isoforms that are expressed in mTECs. The number of isoforms with restricted expression found in mTECs was significantly higher than for comparable cell types. Furthermore, we found evidence that AIRE influences the increased splicing diversity observed in mTECs as the genes for which tissue restricted isoforms are produced in mTECs were significantly more likely than other genes to be differentially spliced between AIRE knock-out and wild-type samples. Our results suggest that developing T lymphocytes are exposed to diverse tissue-restricted splice isoforms in the thymus and that AIRE has a direct or indirect role in this process, representing a novel aspect of its role in the maintenance of immune self-tolerance. Supplementary data are available at Bioinformatics online.

P01.08 * LINEAGE-SPECIFIC SPLICING OF AN ALTERNATIVE EXON OF ANXA7 PROMOTES EGFR SIGNALING ACTIVATION AND TUMOR PROGRESSION IN GLIOBLASTOMA

Tissue-specific alternative splicing is critical to the emergence of tissue identity during development, yet its role in malignant transformation is undefined. Tissue-specific splicing involves evolutionarily-conserved, alternative exons, which represent only a minority of total alternative exons. Many, however, have functional features that influence signaling pathways to profound biological effect. In the brain, Annexin A7 isoform 1 (ANXA7-I1) is exclusively expressed in mature neurons, while isoform 2 (ANXA7-I2) in which exon 6 is skipped, is expressed in glial and progenitor cells. We show that lineage-specific splicing of the cassette exon 6 in the membrane-binding tumor suppressor ANXA7diminishes endosomal targeting and consequent signal termination of the EGFR oncoprotein during brain tumor progression. Splicing of this exon is mediated by Polypyrimidine Tract-Binding Protein 1 (PTBP1), a ribonucleoprotein normally repressed during neuronal development but which we found to be highly expressed also in glioblastomas through loss of a brain-enriched microRNA, miR-124, and gene amplification. Here, we show that the PTBP1-ANXA7 splicing-EGFR signal activation axis promotes in vitro cell migration and invasion, and tumor angiogenesis in vivo. In glioblastoma, ANXA7 splicing is likely inherited from a potential tumor-initiating ancestor but this trait is further exploited through accumulation of mutations that enhance EGFR signaling. Our data illustrate how lineage-specific splicing of a tissue-regulated alternative exon in a constituent of an oncogenic pathway eliminates its tumor suppressor function and promotes glioblastoma progression. This paradigm may offer a general model as to how tissue-specific regulatory mechanisms can contribute to reprogramming normal development to oncogenesis.

RBFOX and SUP-12 sandwich a G base to cooperatively regulate tissue-specific splicing.

Tissue-specific alternative pre-mRNA splicing is often cooperatively regulated by multiple splicing factors, but the structural basis of cooperative RNA recognition is poorly understood. In Caenorhabditis elegans, ligand binding specificity of fibroblast growth factor receptors (FGFRs) is determined by mutually exclusive alternative splicing of the sole FGFR gene, egl-15. Here we determined the solution structure of a ternary complex of the RNA-recognition motif (RRM) domains from the RBFOX protein ASD-1, SUP-12 and their target RNA from egl-15. The two RRM domains cooperatively interact with the RNA by sandwiching a G base to form the stable complex. Multichromatic fluorescence splicing reporters confirmed the requirement of the G and the juxtaposition of the respective cis elements for effective splicing regulation in vivo. Moreover, we identified a new target for the heterologous complex through an element search, confirming the functional significance of the intermolecular coordination.

A comprehensive survey of non-canonical splice sites in the human transcriptome.

We uncovered the diversity of non-canonical splice sites at the human transcriptome using deep transcriptome profiling. We mapped a total of 3.7 billion human RNA-seq reads and developed a set of stringent filters to avoid false non-canonical splice site detections. We identified 184 splice sites with non-canonical dinucleotides and U2/U12-like consensus sequences. We selected 10 of the herein identified U2/U12-like non-canonical splice site events and successfully validated 9 of them via reverse transcriptase-polymerase chain reaction and Sanger sequencing. Analyses of the 184 U2/U12-like non-canonical splice sites indicate that 51% of them are not annotated in GENCODE. In addition, 28% of them are conserved in mouse and 76% are involved in alternative splicing events, some of them with tissue-specific alternative splicing patterns. Interestingly, our analysis identified some U2/U12-like non-canonical splice sites that are converted into canonical splice sites by RNA A-to-I editing. Moreover, the U2/U12-like non-canonical splice sites have a differential distribution of splicing regulatory sequences, which may contribute to their recognition and regulation. Our analysis provides a high-confidence group of U2/U12-like non-canonical splice sites, which exhibit distinctive features among the total human splice sites.

RBFOX2 protein domains and cellular activities.

RBFOX2 (RNA-binding protein, Fox-1 homologue 2)/RBM9 (RNA-binding-motif protein 9)/RTA (repressor of tamoxifen action)/HNRBP2 (hexaribonucleotide-binding protein 2) encodes an RNA-binding protein involved in tissue specific alternative splicing regulation and steroid receptors transcriptional activity. Its ability to regulate specific splicing profiles depending on context has been related to different expression levels of the RBFOX2 protein itself and that of other splicing regulatory proteins involved in the shared modulation of specific genes splicing. However, this cannot be the sole explanation as to why RBFOX2 plays a widespread role in numerous cellular mechanisms from development to cell survival dependent on cell/tissue type. RBFOX2 isoforms with altered protein domains exist. In the present article, we describe the main RBFOX2 protein domains, their importance in the context of splicing and transcriptional regulation and we propose that RBFOX2 isoform distribution may play a fundamental role in RBFOX2-specific cellular effects.

ABERRANT SPLICING OF A BRAIN-ENRICHED ALTERNATIVE EXON ELIMINATES TUMOR SUPPRESSOR FUNCTION AND PROMOTES ONCOGENE FUNCTION DURING BRAIN TUMORIGENESIS

BACKGROUND: Tissue-specific alternative splicing is known to be critical to emergence of tissue identity during development, yet its role in malignant transformation is undefined. Tissue-specific splicing involves evolutionary-conserved, alternative exons, which represent only a minority of total alternative exons. Many, however, have functional features that influence activity in signaling pathways to profound biological effect. Given that tissue-specific splicing has a determinative role in brain development and the enrichment of genes containing tissue-specific exons for proteins with roles in signaling and development, it is thus plausible that changes in such exons could rewire normal neurogenesis towards malignant transformation. METHODS: We used integrated molecular genetic and cell biology analyses, computational biology, animal modeling, and clinical patient profiles to characterize the effect of aberrant splicing of a brain-enriched alternative exon in the membrane-binding tumor suppressor Annexin A7 (ANXA7) on oncogene regulation and brain tumorigenesis. RESULTS: We show that aberrant splicing of a tissue-specific cassette exon in ANXA7 diminishes endosomal targeting and consequent termination of the signal of the EGFR oncoprotein during brain tumorigenesis. Splicing of this exon is mediated by the ribonucleoprotein Polypyrimidine Tract-Binding Protein 1 (PTBP1), which is normally repressed during brain development but, we find, is excessively expressed in glioblastomas through either gene amplification or loss of a neuron-specific microRNA, miR-124. Silencing of PTBP1 attenuates both malignancy and angiogenesis in a stem cell-derived glioblastoma animal model characterized by a high native propensity to generate tumor endothelium or vascular pericytes to support tumor growth. We show that EGFR amplification and PTBP1 overexpression portend a similarly poor clinical outcome, further highlighting the importance of PTBP1-mediated activation of EGFR. CONCLUSIONS: Our data illustrate how anomalous splicing of a tissue-regulated exon in a constituent of an oncogenic signaling pathway eliminates its tumor suppressor function and promotes tumorigenesis. This paradigm of malignant glial transformation as a consequence of tissue-specific alternative exon splicing in a tumor suppressor, may have widespread applicability in explaining how changes in critical tissue-specific regulatory mechanisms reprogram normal development to oncogenesis. SECONDARY CATEGORY: n/a.

Lineage-specific splicing of a brain-enriched alternative exon promotes glioblastoma progression.

Tissue-specific alternative splicing is critical for the emergence of tissue identity during development, yet the role of this process in malignant transformation is undefined. Tissue-specific splicing involves evolutionarily conserved, alternative exons that represent only a minority of the total alternative exons identified. Many of these conserved exons have functional features that influence signaling pathways to profound biological effect. Here, we determined that lineage-specific splicing of a brain-enriched cassette exon in the membrane-binding tumor suppressor annexin A7 (ANXA7) diminishes endosomal targeting of the EGFR oncoprotein, consequently enhancing EGFR signaling during brain tumor progression. ANXA7 exon splicing was mediated by the ribonucleoprotein PTBP1, which is normally repressed during neuronal development. PTBP1 was highly expressed in glioblastomas due to loss of a brain-enriched microRNA (miR-124) and to PTBP1 amplification. The alternative ANXA7 splicing trait was present in precursor cells, suggesting that glioblastoma cells inherit the trait from a potential tumor-initiating ancestor and that these cells exploit this trait through accumulation of mutations that enhance EGFR signaling. Our data illustrate that lineage-specific splicing of a tissue-regulated alternative exon in a constituent of an oncogenic pathway eliminates tumor suppressor functions and promotes glioblastoma progression. This paradigm may offer a general model as to how tissue-specific regulatory mechanisms can reprogram normal developmental processes into oncogenic ones.

Expression of the rodent-specific alternative splice variant of tryptophanyl-tRNA synthetase in murine tissues and cells

Tryptophanyl-tRNA synthetase (TrpRS) catalyzes the aminoacylation of tRNATrp. mRNA of a rodent-specific alternative splice variant of TrpRS (SV-TrpRS), which results in the inclusion of an additional hexapeptide at the C-terminus of full-length TrpRS (FL-TrpRS), has been identified in murine embryonic stem (ES) cells. In the present study, we evaluated the expression of mouse TrpRS mRNA by real-time reverse transcription PCR. We show that SV-TrpRS and FL-TrpRS mRNAs are highly expressed in murine ES cells, embryo, spleen, lung, liver and uterus, and that the relative expression of SV-TrpRS compared to FL-TrpRS is significantly less in the brain. Moreover, we found that interferon-γ increases the expression of TrpRS in a mouse cell line. These results provide the first evidence for tissue-specific expression and alternative splicing of mouse TrpRS.

RBM5 Is a Male Germ Cell Splicing Factor and Is Required for Spermatid Differentiation and Male Fertility

Alternative splicing of precursor messenger RNA (pre-mRNA) is common in mammalian cells and enables the production of multiple gene products from a single gene, thus increasing transcriptome and proteome diversity. Disturbance of splicing regulation is associated with many human diseases; however, key splicing factors that control tissue-specific alternative splicing remain largely undefined. In an unbiased genetic screen for essential male fertility genes in the mouse, we identified the RNA binding protein RBM5 (RNA binding motif 5) as an essential regulator of haploid male germ cell pre-mRNA splicing and fertility. Mice carrying a missense mutation (R263P) in the second RNA recognition motif (RRM) of RBM5 exhibited spermatid differentiation arrest, germ cell sloughing and apoptosis, which ultimately led to azoospermia (no sperm in the ejaculate) and male sterility. Molecular modelling suggested that the R263P mutation resulted in compromised mRNA binding. Within the adult mouse testis, RBM5 localises to somatic and germ cells including spermatogonia, spermatocytes and round spermatids. Through the use of RNA pull down coupled with microarrays, we identified 11 round spermatid-expressed mRNAs as putative RBM5 targets. Importantly, the R263P mutation affected pre-mRNA splicing and resulted in a shift in the isoform ratios, or the production of novel spliced transcripts, of most targets. Microarray analysis of isolated round spermatids suggests that altered splicing of RBM5 target pre-mRNAs affected expression of genes in several pathways, including those implicated in germ cell adhesion, spermatid head shaping, and acrosome and tail formation. In summary, our findings reveal a critical role for RBM5 as a pre-mRNA splicing regulator in round spermatids and male fertility. Our findings also suggest that the second RRM of RBM5 is pivotal for appropriate pre-mRNA splicing.

Integrative analysis of tissue-specific methylation and alternative splicing identifies conserved transcription factor binding motifs

The exact role of intragenic DNA methylation in regulating tissue-specific gene regulation is unclear. Recently, the DNA-binding protein CTCF has been shown to participate in the regulation of alternative splicing in a DNA methylation-dependent manner. To globally evaluate the relationship between DNA methylation and tissue-specific alternative splicing, we performed genome-wide DNA methylation profiling of mouse retina and brain. In protein-coding genes, tissue-specific differentially methylated regions (T-DMRs) were preferentially located in exons and introns. Gene ontology and evolutionary conservation analysis suggest that these T-DMRs are likely to be biologically relevant. More than 14% of alternatively spliced genes were associated with a T-DMR. T-DMR-associated genes were enriched for developmental genes, suggesting that a specific set of alternatively spliced genes may be regulated through DNA methylation. Novel DNA sequences motifs overrepresented in T-DMRs were identified as being associated with positive and/or negative regulation of alternative splicing in a position-dependent context. The majority of these evolutionarily conserved motifs contain a CpG dinucleotide. Some transcription factors, which recognize these motifs, are known to be involved in splicing. Our results suggest that DNA methylation-dependent alternative splicing is widespread and lay the foundation for further mechanistic studies of the role of DNA methylation in tissue-specific splicing regulation.

Alternative Splicing Regulates Biogenesis of miRNAs Located across Exon-Intron Junctions

The initial step in microRNA (miRNA) biogenesis requires processing of the precursor miRNA (pre-miRNA) from a longer primary transcript. Many pre-miRNAs originate from introns, and both a mature miRNA and a spliced RNA can be generated from the same transcription unit. We have identified a mechanism in which RNA splicing negatively regulates the processing of pre-miRNAs that overlap exon-intron junctions. Computational analysis identified dozens of such pre-miRNAs, and experimental validation demonstrated competitive interaction between the Microprocessor complex and the splicing machinery. Tissue-specific alternative splicing regulates maturation of one such miRNA, miR-412, resulting in effects on its targets that code a protein network involved in neuronal cell death processes. This mode of regulation specifically controls maturation of splice-site-overlapping pre-miRNAs but not pre-miRNAs located completely within introns or exons of the same transcript. Our data present a biological role of alternative splicing in regulation of miRNA biogenesis.

Preparation of Nuclear Extracts from HeLa Cells

HeLa cells are the archetypal tissue culture cell line for the preparation of mammalian cell-free systems. These cells, derived from a cervical carcinoma, have been maintained in culture since the late 1940s. They grow with a doubling time of ∼24 h and can be cultured in medium containing fetal bovine serum (optimal serum for growth) or medium containing horse serum (in which they grow more slowly, but this serum is considerably less expensive). Nuclear extracts prepared from these cells have been used to determine the mechanisms of splicing and polyadenylation, and such extracts have been characterized extensively. HeLa cells are usually the cell type of choice for initiating cell-free analysis of nearly any aspect of mammalian gene expression. In some instances (e.g., analysis of tissue-specific alternative splicing), it is necessary to use nuclei from a different cell type. We have found that the protocol described here can be used successfully to prepare active nuclear extracts from a wide variety of tissue culture cells, including Drosophila S2 cells.

Tissue-specific splicing of a ubiquitously expressed transcription factor is essential for muscle differentiation

Alternate splicing contributes extensively to cellular complexity by generating protein isoforms with divergent functions. However, the role of alternate isoforms in development remains poorly understood. Mef2 transcription factors are essential transducers of cell signaling that modulate differentiation of many cell types. Among Mef2 family members, Mef2D is unique, as it undergoes tissue-specific splicing to generate a muscle-specific isoform. Since the ubiquitously expressed (Mef2Dα1) and muscle-specific (Mef2Dα2) isoforms of Mef2D are both expressed in muscle, we examined the relative contribution of each Mef2D isoform to differentiation. Using both in vitro and in vivo models, we demonstrate that Mef2D isoforms act antagonistically to modulate differentiation. While chromatin immunoprecipitation (ChIP) sequencing analysis shows that the Mef2D isoforms bind an overlapping set of genes, only Mef2Dα2 activates late muscle transcription. Mechanistically, the differential ability of Mef2D isoforms to activate transcription depends on their susceptibility to phosphorylation by protein kinase A (PKA). Phosphorylation of Mef2Dα1 by PKA provokes its association with corepressors. Conversely, exon switching allows Mef2Dα2 to escape this inhibitory phosphorylation, permitting recruitment of Ash2L for transactivation of muscle genes. Thus, our results reveal a novel mechanism in which a tissue-specific alternate splicing event has evolved that permits a ubiquitously expressed transcription factor to escape inhibitory signaling for temporal regulation of gene expression.

Prediction of clustered RNA-binding protein motif sites in the mammalian genome.

Sequence-specific interactions of RNA-binding proteins (RBPs) with their target transcripts are essential for post-transcriptional gene expression regulation in mammals. However, accurate prediction of RBP motif sites has been difficult because many RBPs recognize short and degenerate sequences. Here we describe a hidden Markov model (HMM)-based algorithm mCarts to predict clustered functional RBP-binding sites by effectively integrating the number and spacing of individual motif sites, their accessibility in local RNA secondary structures and cross-species conservation. This algorithm learns and quantifies rules of these features, taking advantage of a large number of in vivo RBP-binding sites obtained from cross-linking and immunoprecipitation data. We applied this algorithm to study two representative RBP families, Nova and Mbnl, which regulate tissue-specific alternative splicing through interacting with clustered YCAY and YGCY elements, respectively, and predicted their binding sites in the mouse transcriptome. Despite the low information content in individual motif elements, our algorithm made specific predictions for successful experimental validation. Analysis of predicted sites also revealed cases of extensive and distal RBP-binding sites important for splicing regulation. This algorithm can be readily applied to other RBPs to infer their RNA-regulatory networks. The software is freely available at http://zhanglab.c2b2.columbia.edu/index.php/MCarts.

Distinct Types of Disorder in the Human Proteome: Functional Implications for Alternative Splicing

Intrinsically disordered regions have been associated with various cellular processes and are implicated in several human diseases, but their exact roles remain unclear. We previously defined two classes of conserved disordered regions in budding yeast, referred to as “flexible” and “constrained” conserved disorder. In flexible disorder, the property of disorder has been positionally conserved during evolution, whereas in constrained disorder, both the amino acid sequence and the property of disorder have been conserved. Here, we show that flexible and constrained disorder are widespread in the human proteome, and are particularly common in proteins with regulatory functions. Both classes of disordered sequences are highly enriched in regions of proteins that undergo tissue-specific (TS) alternative splicing (AS), but not in regions of proteins that undergo general (i.e., not tissue-regulated) AS. Flexible disorder is more highly enriched in TS alternative exons, whereas constrained disorder is more highly enriched in exons that flank TS alternative exons. These latter regions are also significantly more enriched in potential phosphosites and other short linear motifs associated with cell signaling. We further show that cancer driver mutations are significantly enriched in regions of proteins associated with TS and general AS. Collectively, our results point to distinct roles for TS alternative exons and flanking exons in the dynamic regulation of protein interaction networks in response to signaling activity, and they further suggest that alternatively spliced regions of proteins are often functionally altered by mutations responsible for cancer.