Computational comparison of common event-based differential splicing tools: practical considerations for laboratory researchers

rMATS-turbo: an efficient and flexible computational tool for alternative splicing analysis of large-scale RNA-seq data

Pre-mRNA alternative splicing is a prevalent mechanism for diversifying eukaryotic transcriptomes and proteomes. Regulated alternative splicing plays a role in many biological processes, and dysregulated alternative splicing is a feature of many human diseases. Short-read RNA sequencing (RNA-seq) is now the standard approach for transcriptome-wide analysis of alternative splicing. Since 2011, our laboratory has developed and maintained Replicate Multivariate Analysis of Transcript Splicing (rMATS), a computational tool for discovering and quantifying alternative splicing events from RNA-seq data. Here we provide a protocol for the contemporary version of rMATS, rMATS-turbo, a fast and scalable re-implementation that maintains the statistical framework and user interface of the original rMATS software, while incorporating a revamped computational workflow with a substantial improvement in speed and data storage efficiency. The rMATS-turbo software scales up to massive RNA-seq datasets with tens of thousands of samples. To illustrate the utility of rMATS-turbo, we describe two representative application scenarios. First, we describe a broadly applicable two-group comparison to identify differential alternative splicing events between two sample groups, including both annotated and novel alternative splicing events. Second, we describe a quantitative analysis of alternative splicing in a large-scale RNA-seq dataset (~1,000 samples), including the discovery of alternative splicing events associated with distinct cell states. We detail the workflow and features of rMATS-turbo that enable efficient parallel processing and analysis of large-scale RNA-seq datasets on a compute cluster. We anticipate that this protocol will help the broad user base of rMATS-turbo make the best use of this software for studying alternative splicing in diverse biological systems.

Loss-of-function mutation in PRMT9 causes abnormal synapse development by dysregulation of RNA alternative splicing

Protein arginine methyltransferase 9 (PRMT9) is a recently identified member of the PRMT family, yet its biological function remains largely unknown. Here, by characterizing an intellectual disability associated PRMT9 mutation (G189R) and establishing a Prmt9 conditional knockout (cKO) mouse model, we uncover an important function of PRMT9 in neuronal development. The G189R mutation abolishes PRMT9 methyltransferase activity and reduces its protein stability. Knockout of Prmt9 in hippocampal neurons causes alternative splicing of ~1900 genes, which likely accounts for the aberrant synapse development and impaired learning and memory in the Prmt9 cKO mice. Mechanistically, we discover a methylation-sensitive protein-RNA interaction between the arginine 508 (R508) of the splicing factor 3B subunit 2 (SF3B2), the site that is exclusively methylated by PRMT9, and the pre-mRNA anchoring site, a cis-regulatory element that is critical for RNA splicing. Additionally, using human and mouse cell lines, as well as an SF3B2 arginine methylation-deficient mouse model, we provide strong evidence that SF3B2 is the primary methylation substrate of PRMT9, thus highlighting the conserved function of the PRMT9/SF3B2 axis in regulating pre-mRNA splicing.

Widespread alternative splicing dysregulation occurs presymptomatically in CAG expansion spinocerebellar ataxias

The spinocerebellar ataxias (SCAs) are a group of dominantly inherited neurodegenerative diseases, several of which are caused by CAG expansion mutations (SCAs 1, 2, 3, 6, 7 and 12) and more broadly belong to the large family of over 40 microsatellite expansion diseases. While dysregulation of alternative splicing is a well defined driver of disease pathogenesis across several microsatellite diseases, the contribution of alternative splicing in CAG expansion SCAs is poorly understood. Furthermore, despite extensive studies on differential gene expression, there remains a gap in our understanding of presymptomatic transcriptomic drivers of disease. We sought to address these knowledge gaps through a comprehensive study of 29 publicly available RNA-sequencing datasets. We identified that dysregulation of alternative splicing is widespread across CAG expansion mouse models of SCAs 1, 3 and 7. These changes were detected presymptomatically, persisted throughout disease progression, were repeat length-dependent, and were present in brain regions implicated in SCA pathogenesis including the cerebellum, pons and medulla. Across disease progression, changes in alternative splicing occurred in genes that function in pathways and processes known to be impaired in SCAs, such as ion channels, synaptic signalling, transcriptional regulation and the cytoskeleton. We validated several key alternative splicing events with known functional consequences, including Trpc3 exon 9 and Kcnma1 exon 23b, in the Atxn1154Q/2Q mouse model. Finally, we demonstrated that alternative splicing dysregulation is responsive to therapeutic intervention in CAG expansion SCAs with Atxn1 targeting antisense oligonucleotide rescuing key splicing events. Taken together, these data demonstrate that widespread presymptomatic dysregulation of alternative splicing in CAG expansion SCAs may contribute to disease onset, early neuronal dysfunction and may represent novel biomarkers across this devastating group of neurodegenerative disorders.

KDM3A and KDM3B Maintain Naïve Pluripotency Through the Regulation of Alternative Splicing

Histone modifying enzymes play a central role in maintaining cell identity by establishing a conducive chromatin environment for lineage specific transcription factor activity. Pluripotent embryonic stem cell (ESC) identity is characterized by a lower abundance of gene repression associated histone modifications that enables rapid response to differentiation cues. The KDM3 family of histone demethylases removes the repressive histone H3 lysine 9 dimethylation (H3K9me2). Here we uncover a surprising role for the KDM3 proteins in the maintenance of the pluripotent state through post-transcriptional regulation. We find that KDM3A and KDM3B interact with RNA processing factors such as EFTUD2 and PRMT5. Acute selective degradation of the endogenous KDM3A and KDM3B proteins resulted in altered splicing independent of H3K9me2 status or catalytic activity. These splicing changes partially resemble the splicing pattern of the more blastocyst-like ground state of pluripotency and occurred in important chromatin and transcription factors such as Dnmt3b, Tbx3 and Tcf12. Our findings reveal non-canonical roles of histone demethylating enzymes in splicing to regulate cell identity.

Identification of Splice Variants and Isoforms in Transcriptomics and Proteomics

Alternative splicing is pivotal to the regulation of gene expression and protein diversity in eukaryotic cells. The detection of alternative splicing events requires specific omics technologies. Although short-read RNA sequencing has successfully supported a plethora of investigations on alternative splicing, the emerging technologies of long-read RNA sequencing and top-down mass spectrometry open new opportunities to identify alternative splicing and protein isoforms with less ambiguity. Here, we summarize improvements in short-read RNA sequencing for alternative splicing analysis, including percent splicing index estimation and differential analysis. We also review the computational methods used in top-down proteomics analysis regarding proteoform identification, including the construction of databases of protein isoforms and statistical analyses of search results. While many improvements in sequencing and computational methods will result from emerging technologies, there should be future endeavors to increase the effectiveness, integration, and proteome coverage of alternative splicing events.

Methotrexate Provokes Disparate Folate Metabolism Gene Expression and Alternative Splicing in Ex Vivo Monocytes and GM-CSF- and M-CSF-Polarized Macrophages

Macrophages constitute important immune cell targets of the antifolate methotrexate (MTX) in autoimmune diseases, including rheumatoid arthritis. Regulation of folate/MTX metabolism remains poorly understood upon pro-inflammatory (M1-type/GM-CSF-polarized) and anti-inflammatory (M2-type/M-CSF-polarized) macrophages. MTX activity strictly relies on the folylpolyglutamate synthetase (FPGS) dependent intracellular conversion and hence retention to MTX-polyglutamate (MTX-PG) forms. Here, we determined FPGS pre-mRNA splicing, FPGS enzyme activity and MTX-polyglutamylation in human monocyte-derived M1- and M2-macrophages exposed to 50 nmol/L MTX ex vivo. Moreover, RNA-sequencing analysis was used to investigate global splicing profiles and differential gene expression in monocytic and MTX-exposed macrophages. Monocytes displayed six-eight-fold higher ratios of alternatively-spliced/wild type FPGS transcripts than M1- and M2-macrophages. These ratios were inversely associated with a six-ten-fold increase in FPGS activity in M1- and M2-macrophages versus monocytes. Total MTX-PG accumulation was four-fold higher in M1- versus M2-macrophages. Differential splicing after MTX-exposure was particularly apparent in M2-macrophages for histone methylation/modification genes. MTX predominantly induced differential gene expression in M1-macrophages, involving folate metabolic pathway genes, signaling pathways, chemokines/cytokines and energy metabolism. Collectively, macrophage polarization-related differences in folate/MTX metabolism and downstream pathways at the level of pre-mRNA splicing and gene expression may account for variable accumulation of MTX-PGs, hence possibly impacting MTX treatment efficacy.

Transcriptomic Complexity of Culm Growth and Development in Different Types of Moso Bamboo

Moso bamboo is capable of both sexual and asexual reproduction during natural growth, resulting in four distinct types of culms: the bamboo shoot-culm, the seedling stem, the leptomorph rhizome, and a long-ignored culm-the outward-rhizome. Sometimes, when the outward rhizomes break through the soil, they continue to grow longitudinally and develop into a new individual. However, the roles of alternative transcription start sites (aTSS) or termination sites (aTTS) as well as alternative splicing (AS) have not been comprehensively studied for their development. To re-annotate the moso bamboo genome and identify genome-wide aTSS, aTTS, and AS in growing culms, we utilized single-molecule long-read sequencing technology. In total, 169,433 non-redundant isoforms and 14,840 new gene loci were identified. Among 1311 lncRNAs, most of which showed a positive correlation with their target mRNAs, one-third of these IncRNAs were preferentially expressed in winter bamboo shoots. In addition, the predominant AS type observed in moso bamboo was intron retention, while aTSS and aTTS events occurred more frequently than AS. Notably, most genes with AS events were also accompanied by aTSS and aTTS events. Outward rhizome growth in moso bamboo was associated with a significant increase in intron retention, possibly due to changes in the growth environment. As different types of moso bamboo culms grow and develop, a significant number of isoforms undergo changes in their conserved domains due to the regulation of aTSS, aTTS, and AS. As a result, these isoforms may play different roles than their original functions. These isoforms then performed different functions from their original roles, contributing to the transcriptomic complexity of moso bamboo. Overall, this study provided a comprehensive overview of the transcriptomic changes underlying different types of moso bamboo culm growth and development.

Unveiling the Role of Dynamic Alternative Splicing Modulation After Infestation with Sea Lice (Caligus rogercresseyi) in Atlantic Salmon

Sea lice are pathogenic marine ectoparasite copepods that represent a severe risk to the worldwide salmon industry. Several transcriptomic investigations have characterized the regulation of gene expression response of Atlantic salmon to sea lice infestation. These studies have focused on the levels of transcript, overlooking the potentially relevant role of alternative splicing (AS), which corresponds to an essential control mechanism of gene expression through RNA processing. In the present study, we performed a genome-wide bioinformatics characterization of differential AS event dynamics in control and infested C. rogercresseyi Atlantic salmon and in resistant and susceptible phenotypes. We identified a significant rise of alternative splicing events and AS genes after infestation and 176 differential alternative splicing events (DASE) from 133 genes. In addition, a higher number of DASE and AS genes were observed among resistant and susceptible phenotypes. Functional annotation of AS genes shows several terms and pathways associated with behavior, RNA splicing, immune response, and RNA binding. Furthermore, three protein-coding genes were identified undergoing differential transcript usage events, among resistant and susceptible phenotypes. Our findings support AS performing a relevant regulatory role in the response of salmonids to sea lice infestation.

Reversal of splicing infidelity is a pre-activation step in B cell differentiation

Introduction

B cell activation and differentiation is central to the adaptive immune response. Changes in exon usage can have major impacts on cellular signaling and differentiation but have not been systematically explored in differentiating B cells.

Methods

We analyzed exon usage and intron retention in RNA-Seq data from subsets of human B cells at various stages of differentiation, and in an in vitro laboratory model of B cell activation and differentiation (Epstein Barr virus infection).

Results

Blood naïve B cells were found to have an unusual splicing profile, with unannotated splicing events in over 30% of expressed genes. Splicing changed substantially upon naïve B cell entry into secondary lymphoid tissue and before activation, involving significant increases in exon commitment and reductions in intron retention. These changes preferentially involved short introns with weak splice sites and were likely mediated by an overall increase in splicing efficiency induced by the lymphoid environment. The majority of transcripts affected by splicing changes showed restoration of encoded conserved protein domains and/or reduced targeting to the nonsense-mediated decay pathway. Affected genes were enriched in functionally important immune cell activation pathways such as antigen-mediated signaling, cell cycle control and mRNA processing and splicing.

Discussion

Functional observations from donor B cell subsets in progressive states of differentiation and from timecourse experiments using the in vitro model suggest that these widespread changes in mRNA splicing play a role in preparing naïve B cells for the decisive step of antigen-mediated activation and differentiation.

Transcriptome analysis provides critical answers to the "variants of uncertain significance" conundrum

While whole-genome and exome sequencing have transformed our collective understanding of genetics' role in disease pathogenesis, there are certain conditions and populations for whom DNA-level data fails to identify the underlying genetic etiology. Specifically, patients of non-White race and non-European ancestry are disproportionately affected by "variants of unknown/uncertain significance" (VUS), limiting the scope of precision medicine for minority patients and perpetuating health disparities. VUS often include deep intronic and splicing variants which are difficult to interpret from DNA data alone. RNA analysis can illuminate the consequences of VUS, thereby allowing for their reclassification as pathogenic versus benign. Here we review the critical role transcriptome analysis plays in clarifying VUS in both neoplastic and non-neoplastic diseases.

Epigenetic Alterations in Cryopreserved Human Spermatozoa: Suspected Potential Functional Defects

Background: Gene set enrichment analysis (GSEA) was conducted on raw data, and alternative splicing (AS) events were found after mRNA sequencing of human spermatozoa. In this study, we aimed to compare unknown micro-epigenetics alternations in fresh and cryopreserved spermatozoa to evaluate the effectivity of cryopreservation protocols. Methods: Spermatozoa were divided into three groups: fresh spermatozoa (group 1), cryoprotectant-free vitrified spermatozoa (group 2), and conventionally frozen spermatozoa (group 3). Nine RNA samples (three replicates in each group) were detected and were used for library preparation with an Illumina compatible kit and sequencing by the Illumina platform. Results: Three Gene Ontology (GO) terms were found to be enriched in vitrified spermatozoa compared with fresh spermatozoa: mitochondrial tRNA aminoacylation, ATP-dependent microtubule motor activity, and male meiotic nuclear division. In alternative splicing analysis, a number of unknown AS events were found, including functional gene exon skipping (SE), alternative 5′ splice sites (A5SS), alternative 3′ splice sites (A3SS), mutually exclusive exon (MXE), and retained intron (RI). Conclusions: Cryopreservation of spermatozoa from some patients can agitate epigenetic instability, including increased alternative splicing events and changes in crucial mitochondrial functional activities. For fertilization of oocytes, for such patients, it is recommended to use fresh spermatozoa whenever possible; cryopreservation of sperm is recommended to be used only in uncontested situations.

The role of alternative splicing in adaptation and evolution

Regulation of gene expression plays a central role in adaptive divergence and evolution. Although the role of gene regulation in microevolutionary processes is gaining wide acceptance, most studies have only investigated the evolution of transcript levels, ignoring the potentially significant role of transcript structures. We argue that variation in alternative splicing plays an important and widely unexplored role in adaptation (e.g., by increasing transcriptome and/or proteome diversity, or buffering potentially deleterious genetic variation). New studies increasingly highlight the potential for independent evolution in alternative splicing and transcript level, providing alternative paths for selection to act upon. We propose that alternative splicing and transcript levels can provide contrasting, nonredundant mechanisms of equal importance for adaptive diversification of gene function and regulation.

Probing Isoform Switching Events in Various Cancer Types: Lessons From Pan-Cancer Studies

Alternative splicing is an essential regulatory mechanism for gene expression in mammalian cells contributing to protein, cellular, and species diversity. In cancer, alternative splicing is frequently disturbed, leading to changes in the expression of alternatively spliced protein isoforms. Advances in sequencing technologies and analysis methods led to new insights into the extent and functional impact of disturbed alternative splicing events. In this review, we give a brief overview of the molecular mechanisms driving alternative splicing, highlight the function of alternative splicing in healthy tissues and describe how alternative splicing is disrupted in cancer. We summarize current available computational tools for analyzing differential transcript usage, isoform switching events, and the pathogenic impact of cancer-specific splicing events. Finally, the strategies of three recent pan-cancer studies on isoform switching events are compared. Their methodological similarities and discrepancies are highlighted and lessons learned from the comparison are listed. We hope that our assessment will lead to new and more robust methods for cancer-specific transcript detection and help to produce more accurate functional impact predictions of isoform switching events.