Juw Won Park

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.

BS69/ZMYND11 Reads and Connects Histone H3.3 Lysine 36 Trimethylation-Decorated Chromatin to Regulated Pre-mRNA Processing

BS69 (also called ZMYND11) contains tandemly arranged PHD, BROMO, and PWWP domains, which are chromatin recognition modalities. Here, we show that BS69 selectively recognizes histone variant H3.3 lysine 36 trimethylation (H3.3K36me3) via its chromatin-binding domains. We further identify BS69 association with RNA splicing regulators, including the U5 snRNP components of the spliceosome, such as EFTUD2. Remarkably, RNA sequencing shows that BS69 mainly regulates intron retention (IR), which is the least understood RNA alternative splicing event in mammalian cells. Biochemical and genetic experiments demonstrate that BS69 promotes IR by antagonizing EFTUD2 through physical interactions. We further show that regulation of IR by BS69 also depends on its binding to H3K36me3-decorated chromatin. Taken together, our study identifies an H3.3K36me3-specific reader and a regulator of IR and reveals that BS69 connects histone H3.3K36me3 to regulated RNA splicing, providing significant, important insights into chromatin regulation of pre-mRNA processing.

GLiMMPS: robust statistical model for regulatory variation of alternative splicing using RNA-seq data.

To characterize the genetic variation of alternative splicing, we develop GLiMMPS, a robust statistical method for detecting splicing quantitative trait loci (sQTLs) from RNA-seq data. GLiMMPS takes into account the individual variation in sequencing coverage and the noise prevalent in RNA-seq data. Analyses of simulated and real RNA-seq datasets demonstrate that GLiMMPS outperforms competing statistical models. Quantitative RT-PCR tests of 26 randomly selected GLiMMPS sQTLs yielded a validation rate of 100%. As population-scale RNA-seq studies become increasingly affordable and popular, GLiMMPS provides a useful tool for elucidating the genetic variation of alternative splicing in humans and model organisms.

The splicing regulators Esrp1 and Esrp2 direct an epithelial splicing program essential for mammalian development

Tissue- and cell-type-specific regulators of alternative splicing (AS) are essential components of posttranscriptional gene regulation, necessary for normal cellular function, patterning, and development. Mice with ablation of Epithelial splicing regulatory protein (Esrp1) develop cleft lip and palate. Loss of both Esrp1 and its paralog Esrp2 results in widespread developmental defects with broad implications to human disease. Deletion of the Esrps in the epidermis revealed their requirement for establishing a proper skin barrier, a primary function of epithelial cells comprising the epidermis. We profiled the global Esrp-mediated splicing regulatory program in epidermis, which revealed large-scale programs of epithelial cell-type-specific splicing required for epithelial cell functions. These mice represent a valuable model for evaluating the essential role for AS in development and function of epithelial cells, which play essential roles in tissue homeostasis in numerous organs, and provide a genetic tool to evaluate important functional properties of epithelial-specific splice variants in vivo. DOI: http://dx.doi.org/10.7554/eLife.08954.001

HDAC inhibition impedes epithelial–mesenchymal plasticity and suppresses metastatic, castration-resistant prostate cancer

PI3K (phosphoinositide 3-kinase)/AKT and RAS/MAPK (mitogen-activated protein kinase) pathway coactivation in the prostate epithelium promotes both epithelial–mesenchymal transition (EMT) and metastatic castration-resistant prostate cancer (mCRPC), which is currently incurable. To study the dynamic regulation of the EMT process, we developed novel genetically defined cellular and in vivo model systems from which epithelial, EMT and mesenchymal-like tumor cells with Pten deletion and Kras activation can be isolated. When cultured individually, each population has the capacity to regenerate all three tumor cell populations, indicative of epithelial–mesenchymal plasticity. Despite harboring the same genetic alterations, mesenchymal-like tumor cells are resistant to PI3K and MAPK pathway inhibitors, suggesting that epigenetic mechanisms may regulate the EMT process, as well as dictate the heterogeneous responses of cancer cells to therapy. Among differentially expressed epigenetic regulators, the chromatin remodeling protein HMGA2 is significantly upregulated in EMT and mesenchymal-like tumors cells, as well as in human mCRPC. Knockdown of HMGA2, or suppressing HMGA2 expression with the histone deacetylase inhibitor LBH589, inhibits epithelial–mesenchymal plasticity and stemness activities in vitro and markedly reduces tumor growth and metastasis in vivo through successful targeting of EMT and mesenchymal-like tumor cells. Importantly, LBH589 treatment in combination with castration prevents mCRPC development and significantly prolongs survival following castration by enhancing p53 and androgen receptor acetylation and in turn sensitizing castration-resistant mesenchymal-like tumor cells to androgen deprivation therapy. Taken together, these findings demonstrate that cellular plasticity is regulated epigenetically, and that mesenchymal-like tumor cell populations in mCRPC that are resistant to conventional and targeted therapies can be effectively treated with the epigenetic inhibitor LBH589.

Determination of a Comprehensive Alternative Splicing Regulatory Network and Combinatorial Regulation by Key Factors during the Epithelial-to-Mesenchymal Transition

ABSTRACT The epithelial-to-mesenchymal transition (EMT) is an essential biological process during embryonic development that is also implicated in cancer metastasis. While the transcriptional regulation of EMT has been well studied, the role of alternative splicing (AS) regulation in EMT remains relatively uncharacterized. We previously showed that the epithelial cell-type-specific proteins epithelial splicing regulatory proteins 1 (ESRP1) and ESRP2 are important for the regulation of many AS events that are altered during EMT. However, the contributions of the ESRPs and other splicing regulators to the AS regulatory network in EMT require further investigation. Here, we used a robust in vitro EMT model to comprehensively characterize splicing switches during EMT in a temporal manner. These investigations revealed that the ESRPs are the major regulators of some but not all AS events during EMT. We determined that the splicing factor RBM47 is downregulated during EMT and also regulates numerous transcripts that switch splicing during EMT. We also determined that Quaking (QKI) broadly promotes mesenchymal splicing patterns. Our study highlights the broad role of posttranscriptional regulation during the EMT and the important role of combinatorial regulation by different splicing factors to fine tune gene expression programs during these physiological and developmental transitions.

Transcriptome sequencing reveals aberrant alternative splicing in Huntington's disease

Huntingtons disease (HD) is an autosomal dominant neurodegenerative disorder caused by a CAG expansion in the gene-encoding Huntingtin (HTT). Transcriptome dysregulation is a major feature of HD pathogenesis, as revealed by a large body of work on gene expression profiling of tissues from human HD patients and mouse models. These studies were primarily focused on transcriptional changes affecting steady-state overall gene expression levels using microarray based approaches. A major missing component, however, has been the study of transcriptome changes at the post-transcriptional level, such as alternative splicing. Alternative splicing is a critical mechanism for expanding regulatory and functional diversity from a limited number of genes, and is particularly complex in the mammalian brain. Here we carried out a deep RNA-seq analysis of the BA4 (Brodmann area 4) motor cortex from seven human HD brains and seven controls to systematically discover aberrant alternative splicing events and characterize potential associated splicing factors in HD. We identified 593 differential alternative splicing events between HD and control brains. Using two expanded panels with a total of 108 BA4 tissues from patients and controls, we identified four splicing factors exhibiting significantly altered expression levels in HD patient brains. Moreover, follow-up molecular analyses of one splicing factor PTBP1 revealed its impact on disease-associated splicing patterns in HD. Collectively, our data provide genomic evidence for widespread splicing dysregulation in HD brains, and suggest the role of aberrant alternative splicing in the pathogenesis of HD.

Transcriptome-wide landscape of pre-mRNA alternative splicing associated with metastatic colonization.

Metastatic colonization is an ominous feature of cancer progression. Recent studies have established the importance of pre-mRNA alternative splicing (AS) in cancer biology. However, little is known about the transcriptome-wide landscape of AS associated with metastatic colonization. Both in vitro and in vivo models of metastatic colonization were utilized to study AS regulation associated with cancer metastasis. Transcriptome profiling of prostate cancer cells and derivatives crossing in vitro or in vivo barriers of metastasis revealed splicing factors with significant gene expression changes associated with metastatic colonization. These include splicing factors known to be differentially regulated in epithelial–mesenchymal transition (ESRP1, ESRP2, and RBFOX2), a cellular process critical for cancer metastasis, as well as novel findings (NOVA1 and MBNL3). Finally, RNA-seq indicated a large network of AS events regulated by multiple splicing factors with altered gene expression or protein activity. These AS events are enriched for pathways important for cell motility and signaling, and affect key regulators of the invasive phenotype such as CD44 and GRHL1. Implications: Transcriptome-wide remodeling of AS is an integral regulatory process underlying metastatic colonization, and AS events affect the metastatic behavior of cancer cells. Mol Cancer Res; 13(2); 305–18. ©2014 AACR.

Identifying differential alternative splicing events from RNA sequencing data using RNASeq-MATS.

RNA sequencing (RNA-Seq) has emerged as a powerful and increasingly cost-effective technology for analysis of transcriptomes. RNA-Seq has several significant advantages over gene expression microarrays, including its high sensitivity and accuracy, broad dynamic range, nucleotide-level resolution, ability to detect novel mRNA transcripts, and ability to analyze pre-mRNA alternative splicing. A major application of RNA-Seq is to detect differential alternative splicing, i.e., differences in exon splicing patterns among different biological conditions. We recently developed a statistical method multivariate analysis of transcript splicing (MATS) for detecting differential alternative splicing events from RNA-Seq data. Here, we describe a computational pipeline RNASeq-MATS based on the MATS algorithm. This pipeline automatically detects and analyzes differential alternative splicing events corresponding to all major types of alternative splicing patterns from RNA-Seq data.

The contribution of Alu exons to the human proteome

BackgroundAlu elements are major contributors to lineage-specific new exons in primate and human genomes. Recent studies indicate that some Alu exons have high transcript inclusion levels or tissue-specific splicing profiles, and may play important regulatory roles in modulating mRNA degradation or translational efficiency. However, the contribution of Alu exons to the human proteome remains unclear and controversial. The prevailing view is that exons derived from young repetitive elements, such as Alu elements, are restricted to regulatory functions and have not had adequate evolutionary time to be incorporated into stable, functional proteins.ResultsWe adopt a proteotranscriptomics approach to systematically assess the contribution of Alu exons to the human proteome. Using RNA sequencing, ribosome profiling, and proteomics data from human tissues and cell lines, we provide evidence for the translational activities of Alu exons and the presence of Alu exon derived peptides in human proteins. These Alu exon peptides represent species-specific protein differences between primates and other mammals, and in certain instances between humans and closely related primates. In the case of the RNA editing enzyme ADARB1, which contains an Alu exon peptide in its catalytic domain, RNA sequencing analyses of A-to-I editing demonstrate that both the Alu exon skipping and inclusion isoforms encode active enzymes. The Alu exon derived peptide may fine tune the overall editing activity and, in limited cases, the site selectivity of ADARB1 protein products.ConclusionsOur data indicate that Alu elements have contributed to the acquisition of novel protein sequences during primate and human evolution.