Identify Alternative Splicing Events Based on Position-Specific Evolutionary Conservation

The evolution of eukaryotes is accompanied by the increased complexity of alternative splicing which greatly expands genome information. One of the greatest challenges in the post-genome era is a complete revelation of human transcriptome with consideration of alternative splicing. Here, we introduce a comparative genomics approach to systemically identify alternative splicing events based on the differential evolutionary conservation between exons and introns and the high-quality annotation of the ENCODE regions. Specifically, we focus on exons that are included in some transcripts but are completely spliced out for others and we call them conditional exons. First, we characterize distinguishing features among conditional exons, constitutive exons and introns. One of the most important features is the position-specific conservation score. There are dramatic differences in conservation scores between conditional exons and constitutive exons. More importantly, the differences are position-specific. For flanking intronic regions, the differences between conditional exons and constitutive exons are also position-specific. Using the Random Forests algorithm, we can classify conditional exons with high specificities (97% for the identification of conditional exons from intron regions and 95% for the classification of known exons) and fair sensitivities (64% and 32% respectively). We applied the method to the human genome and identified 39,640 introns that actually contain conditional exons and classified 8,813 conditional exons from the current RefSeq exon list. Among those, 31,673 introns containing conditional exons and 5,294 conditional exons classified from known exons cannot be inferred from RefSeq, UCSC or Ensembl annotations. Some of these de novo predictions were experimentally verified

Studying alternative splicing regulatory networks through partial correlation analysis

The identification of links between exons and their regulators or targets and between co-spliced exons in human, mouse and rat provides novel insights into the alternative splicing regulatory network

A protein assembly mediates Xist localization and gene silencing

Nuclear compartments have diverse roles in regulating gene expression, yet the molecular forces and components that drive compartment formation remain largely unclear. The long non-coding RNA Xist establishes an intra-chromosomal compartment by localizing at a high concentration in a territory spatially close to its transcription locus and binding diverse proteins to achieve X-chromosome inactivation (XCI). The XCI process therefore serves as a paradigm for understanding how RNA-mediated recruitment of various proteins induces a functional compartment. The properties of the inactive X (Xi)-compartment are known to change over time, because after initial Xist spreading and transcriptional shutoff a state is reached in which gene silencing remains stable even if Xist is turned off. Here we show that the Xist RNA-binding proteins PTBP1, MATR3, TDP-43 and CELF1 assemble on the multivalent E-repeat element of Xist and, via self-aggregation and heterotypic protein–protein interactions, form a condensate in the Xi. This condensate is required for gene silencing and for the anchoring of Xist to the Xi territory, and can be sustained in the absence of Xist. Notably, these E-repeat-binding proteins become essential coincident with transition to the Xist-independent XCI phase, indicating that the condensate seeded by the E-repeat underlies the developmental switch from Xist-dependence to Xist-independence. Taken together, our data show that Xist forms the Xi compartment by seeding a heteromeric condensate that consists of ubiquitous RNA-binding proteins, revealing an unanticipated mechanism for heritable gene silencing

A hierarchical Bayesian model for comparing transcriptomes at the individual transcript isoform level

The complexity of mammalian transcriptomes is compounded by alternative splicing which allows one gene to produce multiple transcript isoforms. However, transcriptome comparison has been limited to differential analysis at the gene level instead of the individual transcript isoform level. High-throughput sequencing technologies and high-resolution tiling arrays provide an unprecedented opportunity to compare transcriptomes at the level of individual splice variants. However, sequence read coverage or probe intensity at each position may represent a family of splice variants instead of one single isoform. Here we propose a hierarchical Bayesian model, BASIS (Bayesian Analysis of Splicing IsoformS), to infer the differential expression level of each transcript isoform in response to two conditions. A latent variable was introduced to perform direct statistical selection of differentially expressed isoforms. Model parameters were inferred based on an ergodic Markov chain generated by our Gibbs sampler. BASIS has the ability to borrow information across different probes (or positions) from the same genes and different genes. BASIS can handle the heteroskedasticity of probe intensity or sequence read coverage. We applied BASIS to a human tiling-array data set and a mouse RNA-seq data set. Some of the predictions were validated by quantitative real-time RT–PCR experiments

A Cell-Based High-Throughput Method for Identifying Modulators of Alternative Splicing

Alternative splicing, a key regulatory process of gene expression, is controlled by trans-acting factors that recognize cis-elements in premature RNA transcripts to affect spliceosome assembly and splice site choices. Extracellular stimuli and signaling cascades can converge on RNA binding splicing regulators to affect alternative splicing. Defects in splicing regulation have been associated with various human diseases, and modification of disease-causing splicing events presents great therapeutic promise. Determining splicing regulators and/or upstream modulators has been difficult with low throughput, low sensitivity, and low specificity. IRAS (identifying regulators of alternative splicing) is a novel cell-based high-throughput screening strategy designed specifically to address these challenges and has achieved high throughput, high sensitivity, and high specificity. Here, we describe the IRAS method in detail with a pair of dual-fluorescence minigene reporters that produces GFP and RFP fluorescent signals to assay the two spliced isoforms exclusively. These two complementary mini-gene reporters alter GFP/RFP output ratios in the opposite direction in response to only a true splicing change. False positives from a signal screen do not stimulate opposite changes in GFP/RFP ratios. The reporter pair in conjunction with robotic liquid handlers and arrayed libraries allows IRAS to screen for both positive and negative splicing regulators with high sensitivity and specificity in a high-throughput manner

Putative Coiled-Coil Domain-Dependent Autoinhibition and Alternative Splicing Determine SHTN1’s Actin-Binding Activity

The actin cytoskeleton plays a pivotal role in cell development, morphogenesis, and other cellular functions. Precise control of actin dynamics requires actin-binding proteins. Here, we characterize multifarious regulation of SHTN1 (shootin1) and show that, unlike known actin-binding proteins, SHTN1's actin binding activity is intrinsically inhibited by a putative coiled-coil domain (CCD) and the autoinhibition is overcome by alternative splicing regulation. We found SHTN1 contains a noncanonical WH2 domain and an upstream proline-rich region (PRR) that by themselves are sufficient for actin interaction. Alternative splicing of Shtn1 at the C terminus and downstream of the WH2-PRR domain produces a long (SHTN1L or shootin1b) and a short (SHTN1S or shootin1a) isoform, which both contain the described PRR and WH2 domains. However, SHTN1S does not interact with actin due to inhibition mediated by an N-terminal CCD. A SHTN1L-specific C-terminal motif counters the intramolecular inhibition and allows SHNT1L to bind actin. A nuclear localization signal is embedded between PRR and WH2 and is subject to similar autoinhibition. SHTN1 would be the first WH2-containing molecule that adopts CCD-dependent autoinhibition and alternative splicing-dependent actin interaction