Shihao Shen

Nonallele-specific Silencing of Mutant and Wild-type Huntingtin Demonstrates Therapeutic Efficacy in Huntington's Disease Mice

Huntingtons disease (HD) is a fatal neurodegenerative disease caused by mutant huntingtin (htt) protein, and there are currently no effective treatments. Recently, we and others demonstrated that silencing mutant htt via RNA interference (RNAi) provides therapeutic benefit in HD mice. We have since found that silencing wild-type htt in adult mouse striatum is tolerated for at least 4 months. However, given the role of htt in various cellular processes, it remains unknown whether nonallele-specific silencing of both wild-type and mutant htt is a viable therapeutic strategy for HD. Here, we tested whether cosilencing wild-type and mutant htt provides therapeutic benefit and is tolerable in HD mice. After treatment, HD mice showed significant reductions in wild-type and mutant htt, and demonstrated improved motor coordination and survival. We performed transcriptional profiling to evaluate the effects of reducing wild-type htt in adult mouse striatum. We identified gene expression changes that are concordant with previously described roles for htt in various cellular processes. Also, several abnormally expressed transcripts associated with early-stage HD were differentially expressed in our studies, but intriguingly, those involved in neuronal function changed in opposing directions. Together, these encouraging and surprising findings support further testing of nonallele-specific RNAi therapeutics for HD.

An ESRP-regulated splicing programme is abrogated during the epithelial-mesenchymal transition

Alternative splicing achieves coordinated changes in post‐transcriptional gene expression programmes through the activities of diverse RNA‐binding proteins. Epithelial splicing regulatory proteins 1 and 2 (ESRP1 and ESRP2) are cell‐type‐specific regulators of transcripts that switch splicing during the epithelial–mesenchymal transition (EMT). To define a comprehensive programme of alternative splicing that is regulated during the EMT, we identified an extensive ESRP‐regulated splicing network of hundreds of alternative splicing events within numerous genes with functions in cell–cell adhesion, polarity, and migration. Loss of this global ESRP‐regulated epithelial splicing programme induces the phenotypic changes in cell morphology that are observed during the EMT. Components of this splicing signature provide novel molecular markers that can be used to characterize the EMT. Bioinformatics and experimental approaches revealed a high‐affinity ESRP‐binding motif and a predictive RNA map that governs their activity. This work establishes the ESRPs as coordinators of a complex alternative splicing network that adds an important post‐transcriptional layer to the changes in gene expression that underlie epithelial–mesenchymal transitions during development and disease.

rMATS: Robust and flexible detection of differential alternative splicing from replicate RNA-Seq data

Significance Alternative splicing (AS) is an important mechanism of eukaryotic gene regulation. Deep RNA sequencing (RNA-Seq) has become a powerful approach for quantitative profiling of AS. With the increasing capacity of high-throughput sequencers, it has become common for RNA-Seq studies of AS to examine multiple biological replicates. We developed rMATS, a new statistical method for robust and flexible detection of differential AS from replicate RNA-Seq data. Besides the analysis of unpaired replicates, rMATS includes a model specifically designed for paired replicates, such as case–control matched pairs in clinical RNA-Seq datasets. We expect rMATS will be useful for genome-wide studies of AS in diverse research projects. Our data also provide new insights about the experimental design for RNA-Seq studies of AS. Ultra-deep RNA sequencing (RNA-Seq) has become a powerful approach for genome-wide analysis of pre-mRNA alternative splicing. We previously developed multivariate analysis of transcript splicing (MATS), a statistical method for detecting differential alternative splicing between two RNA-Seq samples. Here we describe a new statistical model and computer program, replicate MATS (rMATS), designed for detection of differential alternative splicing from replicate RNA-Seq data. rMATS uses a hierarchical model to simultaneously account for sampling uncertainty in individual replicates and variability among replicates. In addition to the analysis of unpaired replicates, rMATS also includes a model specifically designed for paired replicates between sample groups. The hypothesis-testing framework of rMATS is flexible and can assess the statistical significance over any user-defined magnitude of splicing change. The performance of rMATS is evaluated by the analysis of simulated and real RNA-Seq data. rMATS outperformed two existing methods for replicate RNA-Seq data in all simulation settings, and RT-PCR yielded a high validation rate (94%) in an RNA-Seq dataset of prostate cancer cell lines. Our data also provide guiding principles for designing RNA-Seq studies of alternative splicing. We demonstrate that it is essential to incorporate biological replicates in the study design. Of note, pooling RNAs or merging RNA-Seq data from multiple replicates is not an effective approach to account for variability, and the result is particularly sensitive to outliers. The rMATS source code is freely available at rnaseq-mats.sourceforge.net/. As the popularity of RNA-Seq continues to grow, we expect rMATS will be useful for studies of alternative splicing in diverse RNA-Seq projects.

MATS: a Bayesian framework for flexible detection of differential alternative splicing from RNA-Seq data

Ultra-deep RNA sequencing has become a powerful approach for genome-wide analysis of pre-mRNA alternative splicing. We develop MATS (multivariate analysis of transcript splicing), a Bayesian statistical framework for flexible hypothesis testing of differential alternative splicing patterns on RNA-Seq data. MATS uses a multivariate uniform prior to model the between-sample correlation in exon splicing patterns, and a Markov chain Monte Carlo (MCMC) method coupled with a simulation-based adaptive sampling procedure to calculate the P-value and false discovery rate (FDR) of differential alternative splicing. Importantly, the MATS approach is applicable to almost any type of null hypotheses of interest, providing the flexibility to identify differential alternative splicing events that match a given user-defined pattern. We evaluated the performance of MATS using simulated and real RNA-Seq data sets. In the RNA-Seq analysis of alternative splicing events regulated by the epithelial-specific splicing factor ESRP1, we obtained a high RT–PCR validation rate of 86% for differential exon skipping events with a MATS FDR of <10%. Additionally, over the full list of RT–PCR tested exons, the MATS FDR estimates matched well with the experimental validation rate. Our results demonstrate that MATS is an effective and flexible approach for detecting differential alternative splicing from RNA-Seq data.

The epithelial splicing factors ESRP1 and ESRP2 positively and negatively regulate diverse types of alternative splicing events

Cell-type and tissue-specific alternative splicing events are regulated by combinatorial control involving both abundant RNA binding proteins as well as those with more discrete expression and specialized functions. Epithelial Splicing Regulatory Proteins 1 and 2 (ESRP1 and ESRP2) are recently discovered epithelial-specific RNA binding proteins that promote splicing of the epithelial variant of the FGFR2, ENAH, CD44, and CTNND1 transcripts. To cataloge a larger set of splicing events under the regulation of the ESRPs we profiled splicing changes induced by RNA interference-mediated knockdown of ESRP1 and ESRP2 expression in a human epithelial cell line using the splicing sensitive Affymetrix Exon ST1.0 Arrays. Analysis of the microarray data resulted in the identification of over a hundred candidate ESRP regulated splicing events. We were able to independently validate 37 of these targets by RT-PCR. The ESRP regulated events encompass all known types of alternative splicing events, most prominent being alternative cassette exons and splicing events leading to alternative 3’ terminal exons. Importantly, a number of these regulated splicing events occur in gene transcripts that encode proteins with well-described roles in the regulation of actin cytoskeleton organization, cell-cell adhesion, cell polarity, and cell migration. In sum, this work reveals a novel list of transcripts differentially spliced in epithelial and mesenchymal cells, implying that coordinated alternative splicing plays a critical role in determination of cell type identity. These results further establish ESRP1 and ESRP2 as global regulators of an epithelial splicing regulatory network.

Widespread establishment and regulatory impact of Alu exons in human genes

The Alu element has been a major source of new exons during primate evolution. Thousands of human genes contain spliced exons derived from Alu elements. However, identifying Alu exons that have acquired genuine biological functions remains a major challenge. We investigated the creation and establishment of Alu exons in human genes, using transcriptome profiles of human tissues generated by high-throughput RNA sequencing (RNA-Seq) combined with extensive RT-PCR analysis. More than 25% of Alu exons analyzed by RNA-Seq have estimated transcript inclusion levels of at least 50% in the human cerebellum, indicating widespread establishment of Alu exons in human genes. Genes encoding zinc finger transcription factors have significantly higher levels of Alu exonization. Importantly, Alu exons with high splicing activities are strongly enriched in the 5′-UTR, and two-thirds (10/15) of 5′-UTR Alu exons tested by luciferase reporter assays significantly alter mRNA translational efficiency. Mutational analysis reveals the specific molecular mechanisms by which newly created 5′-UTR Alu exons modulate translational efficiency, such as the creation or elongation of upstream ORFs that repress the translation of the primary ORFs. This study presents genomic evidence that a major functional consequence of Alu exonization is the lineage-specific evolution of translational regulation. Moreover, the preferential creation and establishment of Alu exons in zinc finger genes suggest that Alu exonization may have globally affected the evolution of primate and human transcriptomes by regulating the protein production of master transcriptional regulators in specific lineages.

Genome-Wide Determination of a Broad ESRP-Regulated Posttranscriptional Network by High-Throughput Sequencing

ABSTRACT Tissue-specific alternative splicing is achieved through the coordinated assembly of RNA binding proteins at specific sites to enhance or silence splicing at nearby splice sites. We used high-throughput sequencing (RNA-Seq) to investigate the complete spectrum of alternative splicing events that are regulated by the epithelium-specific splicing regulatory proteins ESRP1 and ESRP2. We also combined this analysis with direct RNA sequencing (DRS) to reveal ESRP-mediated regulation of alternative polyadenylation. To define binding motifs that mediate direct regulation of splicing and polyadenylation by ESRP, SELEX-Seq analysis was performed, coupling traditional SELEX with high-throughput sequencing. Identification and scoring of high-affinity ESRP1 binding motifs within ESRP target genes allowed the generation of RNA maps that define the position-dependent activity of the ESRPs in regulating cassette exons and alternative 3′ ends. These extensive analyses provide a comprehensive picture of the functions of the ESRPs in an epithelial posttranscriptional gene expression program.

Diverse Splicing Patterns of Exonized Alu Elements in Human Tissues

Exonization of Alu elements is a major mechanism for birth of new exons in primate genomes. Prior analyses of expressed sequence tags show that almost all Alu-derived exons are alternatively spliced, and the vast majority of these exons have low transcript inclusion levels. In this work, we provide genomic and experimental evidence for diverse splicing patterns of exonized Alu elements in human tissues. Using Exon array data of 330 Alu-derived exons in 11 human tissues and detailed RT-PCR analyses of 38 exons, we show that some Alu-derived exons are constitutively spliced in a broad range of human tissues, and some display strong tissue-specific switch in their transcript inclusion levels. Most of such exons are derived from ancient Alu elements in the genome. In SEPN1, mutations of which are linked to a form of congenital muscular dystrophy, the muscle-specific inclusion of an Alu-derived exon may be important for regulating SEPN1 activity in muscle. Realtime qPCR analysis of this SEPN1 exon in macaque and chimpanzee tissues indicates human-specific increase in its transcript inclusion level and muscle specificity after the divergence of humans and chimpanzees. Our results imply that some Alu exonization events may have acquired adaptive benefits during the evolution of primate transcriptomes.

The nuclear experience of CPEB: Implications for RNA processing and translational control

CPEB is a sequence-specific RNA binding protein that promotes polyadenylation-induced translation in early development, during cell cycle progression and cellular senescence, and following neuronal synapse stimulation. It controls polyadenylation and translation through other interacting molecules, most notably the poly(A) polymerase Gld2, the deadenylating enzyme PARN, and the eIF4E-binding protein Maskin. Here, we report that CPEB shuttles between the nucleus and cytoplasm and that its export occurs via the CRM1-dependent pathway. In the nucleus of Xenopus oocytes, CPEB associates with lampbrush chromosomes and several proteins involved in nuclear RNA processing. CPEB also interacts with Maskin in the nucleus as well as with CPE-containing mRNAs. Although the CPE does not regulate mRNA export, it influences the degree to which mRNAs are translationally repressed in the cytoplasm. Moreover, CPEB directly or indirectly mediates the alternative splicing of at least one pre-mRNA in mouse embryo fibroblasts as well as certain mouse tissues. We propose that CPEB, together with Maskin, binds mRNA in the nucleus to ensure tight translational repression upon export to the cytoplasm. In addition, we propose that nuclear CPEB regulates specific pre-mRNA alternative splicing.

Evolution of alternative splicing in primate brain transcriptomes

Alternative splicing is a predominant form of gene regulation in higher eukaryotes. The evolution of alternative splicing provides an important mechanism for the acquisition of novel gene functions. In this work, we carried out a genome-wide phylogenetic survey of lineage-specific splicing patterns in the primate brain, via high-density exon junction array profiling of brain transcriptomes of humans, chimpanzees and rhesus macaques. We identified 509 genes showing splicing differences among these species. RT-PCR analysis of 40 exons confirmed the predicted splicing evolution of 33 exons. Of these 33 exons, outgroup analysis using rhesus macaques confirmed 13 exons with human-specific increase or decrease in transcript inclusion levels after humans diverged from chimpanzees. Some of the human-specific brain splicing patterns disrupt domains critical for protein-protein interactions, and some modulate translational efficiency of their host genes. Strikingly, for exons showing splicing differences across species, we observed a significant increase in the rate of silent substitutions within exons, coupled with accelerated sequence divergence in flanking introns. This indicates that evolution of cis-regulatory signals is a major contributor to the emergence of human-specific splicing patterns. In one gene (MAGOH), using minigene reporter assays, we demonstrated that the combination of two human-specific cis-sequence changes created its human-specific splicing pattern. Together, our data reveal widespread human-specific changes of alternative splicing in the brain and suggest an important role of splicing in the evolution of neuronal gene regulation and functions.