Functional and dynamic profiling of transcript isoforms reveals essential roles of alternative splicing in interferon response.

Type I interferon (IFN-I) plays an important role in the innate immune response through inducing IFN-I-stimulated genes (ISGs). However, how alternative splicing (AS) events, especially over time, affect their function remains poorly understood. We generated an annotation (113,843 transcripts) for IFN-I-stimulated human B cells called isoISG using high-accuracy long-read sequencing data from PacBio Sequel II/IIe. Transcript isoform profiling using isoISG revealed that isoform switching occurred in the early response to IFN-I so that ISGs would gain functional domains (e.g., C4B) or higher protein production (e.g., IRF3). Conversely, isoforms lacking functional domains increased during the late phase of IFN-I response, mainly due to intron retention events. This suggests that isoform switching both triggers and terminates IFN-I responses at the translation and protein levels. Furthermore, genetic variants influencing the isoform ratio of ISGs were associated with immunological and infectious diseases. AS has essential roles in regulating innate immune response and associated diseases.

Isoform switching leads to downregulation of cytokine producing genes in estrogen receptor positive breast cancer.

Objective: Estrogen receptor breast cancer (BC) is characterized by the expression of estrogen receptors. It is the most common cancer among women, with an incidence rate of 2.26 million cases worldwide. The aim of this study was to identify differentially expressed genes and isoform switching between estrogen receptor positive and triple negative BC samples. Methods: The data were collected from ArrayExpress, followed by preprocessing and subsequent mapping from HISAT2. Read quantification was performed by StringTie, and then R package ballgown was used to perform differential expression analysis. Functional enrichment analysis was conducted using Enrichr, and then immune genes were shortlisted based on the ScType marker database. Isoform switch analysis was also performed using the IsoformSwitchAnalyzeR package. Results: A total of 9,771 differentially expressed genes were identified, of which 86 were upregulated and 117 were downregulated. Six genes were identified as mainly associated with estrogen receptor positive BC, while a novel set of ten genes were found which have not previously been reported in estrogen receptor positive BC. Furthermore, alternative splicing and subsequent isoform usage in the immune system related genes were determined. Conclusion: This study identified the differential usage of isoforms in the immune system related genes in cancer cells that suggest immunosuppression due to the dysregulation of CXCR chemokine receptor binding, iron ion binding, and cytokine activity.

Quantitative Analysis of Isoform Switching in Cancer.

Over the past 8 years, multiple studies examined the phenomenon of isoform switching in human cancers and discovered that isoform switching is widespread, with hundreds to thousands of such events per cancer type. Although all of these studies used slightly different definitions of isoform switching, which in part led to a rather poor overlap of their results, they all leveraged transcript usage, a proportion of the transcript's expression in the total expression level of the parent gene, to detect isoform switching. However, how changes in transcript usage correlate with changes in transcript expression is not sufficiently explored. In this article, we adopt the most common definition of isoform switching and use a state-of-the-art tool for the analysis of differential transcript usage, SatuRn, to detect isoform switching events in 12 cancer types. We analyze the detected events in terms of changes in transcript usage and the relationship between transcript usage and transcript expression on a global scale. The results of our analysis suggest that the relationship between changes in transcript usage and changes in transcript expression is far from straightforward, and that such quantitative information can be effectively used for prioritizing isoform switching events for downstream analyses.

The landscape of differential splicing and transcript alternations in severe COVID-19 infection.

Viral infections can modulate the widespread alternations of cellular splicing, favouring viral replication within the host cells by overcoming host immune responses. However, how SARS-CoV-2 induces host cell differential splicing and affects the landscape of transcript alternation in severe COVID-19 infection remains elusive. Understanding the differential splicing and transcript alternations in severe COVID-19 infection may improve our molecular insights into the SARS-CoV-2 pathogenesis. In this study, we analysed the publicly available blood and lung transcriptome data of severe COVID-19 patients, blood transcriptome data of recovered COVID-19 patients at 12-, 16- and 24-week postinfection and healthy controls. We identified a significant transcript isoform switching in the individual blood and lung RNA-seq data of severe COVID-19-infected patients and 25 common genes that alter their transcript isoform in both blood and lung samples. Altered transcripts show significant loss of the open reading frame, functional domains and switch from coding to noncoding transcript, impacting normal cellular functions. Furthermore, we identified the expression of several novel recurrent chimeric transcripts in the blood samples from severe COVID-19 patients. Moreover, the analysis of the isoform switching into blood samples from recovered COVID-19 patients highlights that there is no significant isoform switching in 16- and 24-week postinfection, and the levels of expressed chimeric transcripts are reduced. This finding emphasizes that SARS-CoV-2 severe infection induces widespread splicing in the host cells, which could help the virus alter the host immune responses and facilitate the viral replication within the host and the efficient translation of viral proteins.

The coordinated regulation of early meiotic stages is dominated by non-coding RNAs and stage-specific transcription in wheat.

Reproductive success hinges on precisely coordinated meiosis, yet our understanding of how structural rearrangements of chromatin and phase transitions during meiosis are transcriptionally regulated is limited. In crop plants, detailed analysis of the meiotic transcriptome could identify regulatory genes and epigenetic regulators that can be targeted to increase recombination rates and broaden genetic variation, as well as provide a resource for comparison among eukaryotes of different taxa to answer outstanding questions about meiosis. We conducted a meiotic stage-specific analysis of messenger RNA (mRNA), small non-coding RNA (sncRNA), and long intervening/intergenic non-coding RNA (lincRNA) in wheat (Triticum aestivum L.) and revealed novel mechanisms of meiotic transcriptional regulation and meiosis-specific transcripts. Amidst general repression of mRNA expression, significant enrichment of ncRNAs was identified during prophase I relative to vegetative cells. The core meiotic transcriptome was comprised of 9309 meiosis-specific transcripts, 48 134 previously unannotated meiotic transcripts, and many known and novel ncRNAs differentially expressed at specific stages. The abundant meiotic sncRNAs controlled the reprogramming of central metabolic pathways by targeting genes involved in photosynthesis, glycolysis, hormone biosynthesis, and cellular homeostasis, and lincRNAs enhanced the expression of nearby genes. Alternative splicing was not evident in this polyploid species, but isoforms were switched at phase transitions. The novel, stage-specific regulatory controls uncovered here challenge the conventional understanding of this crucial biological process and provide a new resource of requisite knowledge for those aiming to directly modulate meiosis to improve crop plants. The wheat meiosis transcriptome dataset can be queried for genes of interest using an eFP browser located at https://bar.utoronto.ca/efp_wheat/cgi-bin/efpWeb.cgi?dataSource=Wheat_Meiosis.

Spatial and temporal requirement of Mlp60A isoforms during muscle development and function in Drosophila melanogaster.

Many myofibrillar proteins undergo isoform switching in a spatio-temporal manner during muscle development. The biological significance of the variants of several of these myofibrillar proteins remains elusive. One such myofibrillar protein, the Muscle LIM Protein (MLP), is a vital component of the Z-discs. In this paper, we show that one of the Drosophila MLP encoding genes, Mlp60A, gives rise to two isoforms: a short (279 bp, 10 kDa) and a long (1461 bp, 54 kDa) one. The short isoform is expressed throughout development, but the long isoform is adult-specific, being the dominant of the two isoforms in the indirect flight muscles (IFMs). A concomitant, muscle-specific knockdown of both isoforms leads to partial developmental lethality, with most of the surviving flies being flight defective. A global loss of both isoforms in a Mlp60A-null background also leads to developmental lethality, with muscle defects in the individuals that survive to the third instar larval stage. This lethality could be rescued partially by a muscle-specific overexpression of the short isoform. Genetic perturbation of only the long isoform, through a P-element insertion in the long isoform-specific coding sequence, leads to defective flight, in around 90% of the flies. This phenotype was completely rescued when the P-element insertion was precisely excised from the locus. Hence, our data show that the two Mlp60A isoforms are functionally specialized: the short isoform being essential for normal embryonic muscle development and the long isoform being necessary for normal adult flight muscle function.

Characterizing isoform switching events in esophageal adenocarcinoma.

Isoform switching events with predicted functional consequences are common in many cancers, but characterization of switching events in esophageal adenocarcinoma (EAC) is lacking. Next-generation sequencing was used to detect levels of RNA transcripts and identify specific isoforms in treatment-naïve esophageal tissues ranging from premalignant Barrett's esophagus (BE), BE with low- or high-grade dysplasia (BE.LGD, BE.HGD), and EAC. Samples were stratified by histopathology and TP53 mutation status, identifying significant isoform switching events with predicted functional consequences. Comparing BE.LGD with BE.HGD, a histopathology linked to cancer progression, isoform switching events were identified in 75 genes including KRAS, RNF128, and WRAP53. Stratification based on TP53 status increased the number of significant isoform switches to 135, suggesting switching events affect cellular functions based on TP53 mutation and tissue histopathology. Analysis of isoforms agnostic, exclusive, and shared with mutant TP53 revealed unique signatures including demethylation, lipid and retinoic acid metabolism, and glucuronidation, respectively. Nearly half of isoform switching events were identified without significant gene-level expression changes. Importantly, two TP53-interacting isoforms, RNF128 and WRAP53, were significantly linked to patient survival. Thus, analysis of isoform switching events may provide new insight for the identification of prognostic markers and inform new potential therapeutic targets for EAC.

circRAB3IP modulates cell proliferation by reorganizing gene expression and mRNA processing in a paracrine manner.

Circular RNAs are an endogenous long-lived and abundant noncoding species. Despite their prevalence, only a few circRNAs have been dissected mechanistically to date. Here, we cataloged nascent RNA-enriched circRNAs from primary human cells and functionally assigned a role to circRAB3IP in sustaining cellular homeostasis. We combined "omics" and functional experiments to show how circRAB3IP depletion deregulates hundreds of genes, suppresses cell cycle progression, and induces senescence-associated gene expression changes. Conversely, excess circRAB3IP delivered to endothelial cells via extracellular vesicles suffices for accelerating their division. We attribute these effects to an interplay between circRAB3IP and the general splicing factor SF3B1, which can affect transcript variant expression levels of cell cycle-related genes. Together, our findings link the maintenance of cell homeostasis to the presence of a single circRNA.

Target isoforms are an overlooked challenge and opportunity in chimeric antigen receptor cell therapy.

The development of novel chimeric antigen receptor (CAR) cell therapies is rapidly growing, with 299 new agents being reported and 109 new clinical trials initiated so far this year. One critical lesson from approved CD19-specific CAR therapies is that target isoform switching has been shown to cause tumour relapse, but little is known about the isoforms of CAR targets in solid cancers. Here we assess the protein isoform landscape and identify both the challenges and opportunities protein isoform switching present as CAR therapy is applied to solid cancers.

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.

Neural Precursor Cells Expanded Inside the 3D Micro-Scaffold Nichoid Present Different Non-Coding RNAs Profiles and Transcript Isoforms Expression: Possible Epigenetic Modulation by 3D Growth.

Non-coding RNAs show relevant implications in various biological and pathological processes. Thus, understanding the biological implications of these molecules in stem cell biology still represents a major challenge. The aim of this work is to study the transcriptional dysregulation of 357 non-coding genes, found through RNA-Seq approach, in murine neural precursor cells expanded inside the 3D micro-scaffold Nichoid versus standard culture conditions. Through weighted co-expression network analysis and functional enrichment, we highlight the role of non-coding RNAs in altering the expression of coding genes involved in mechanotransduction, stemness, and neural differentiation. Moreover, as non-coding RNAs are poorly conserved between species, we focus on those with human homologue sequences, performing further computational characterization. Lastly, we looked for isoform switching as possible mechanism in altering coding and non-coding gene expression. Our results provide a comprehensive dissection of the 3D scaffold Nichoid's influence on the biological and genetic response of neural precursor cells. These findings shed light on the possible role of non-coding RNAs in 3D cell growth, indicating that also non-coding RNAs are implicated in cellular response to mechanical stimuli.

A dynamic intron retention program regulates the expression of several hundred genes during pollen meiosis.

KEY MESSAGE: Intron retention is a stage-specific mechanism of functional attenuation of a subset of co-regulated, functionally related genes during early stages of pollen development. To improve our understanding of the gene regulatory mechanisms that drive developmental processes, we performed a genome-wide study of alternative splicing and isoform switching during five key stages of pollen development in field mustard, Brassica rapa. Surprisingly, for several hundred genes (12.3% of the genes analysed), isoform switching results in stage-specific expression of intron-retaining transcripts at the meiotic stage of pollen development. In such cases, we report temporally regulated switching between expression of a canonical, translatable isoform and an intron-retaining transcript that is predicted to produce a truncated and presumably inactive protein. The results suggest a new pervasive mechanism underlying modulation of protein levels in a plant developmental program. The effect is not based on gene expression induction but on the type of transcript produced. We conclude that intron retention is a stage-specific mechanism of functional attenuation of a subset of co-regulated, functionally related genes during meiosis, especially genes related to ribosome biogenesis, mRNA transport and nuclear envelope architecture. We also propose that stage-specific expression of a non-functional isoform of Brassica rapa BrSDG8, a non-redundant member of histone methyltransferase gene family, linked to alternative splicing regulation, may contribute to the intron retention observed.

Alternative Splicing Regulator RBM20 and Cardiomyopathy.

RBM20 is a vertebrate-specific RNA-binding protein with two zinc finger (ZnF) domains, one RNA-recognition motif (RRM)-type RNA-binding domain and an arginine/serine (RS)-rich region. RBM20 has initially been identified as one of dilated cardiomyopathy (DCM)-linked genes. RBM20 is a regulator of heart-specific alternative splicing and Rbm20 ΔRRM mice lacking the RRM domain are defective in the splicing regulation. The Rbm20 ΔRRM mice, however, do not exhibit a characteristic DCM-like phenotype such as dilatation of left ventricles or systolic dysfunction. Considering that most of the RBM20 mutations identified in familial DCM cases were heterozygous missense mutations in an arginine-serine-arginine-serine-proline (RSRSP) stretch whose phosphorylation is crucial for nuclear localization of RBM20, characterization of a knock-in animal model is awaited. One of the major targets for RBM20 is the TTN gene, which is comprised of the largest number of exons in mammals. Alternative splicing of the TTN gene is exceptionally complicated and RBM20 represses >160 of its consecutive exons, yet detailed mechanisms for such extraordinary regulation are to be elucidated. The TTN gene encodes the largest known protein titin, a multi-functional sarcomeric structural protein specific to striated muscles. As titin is the most important factor for passive tension of cardiomyocytes, extensive heart-specific and developmentally regulated alternative splicing of the TTN pre-mRNA by RBM20 plays a critical role in passive stiffness and diastolic function of the heart. In disease models with diastolic dysfunctions, the phenotypes were rescued by increasing titin compliance through manipulation of the Ttn pre-mRNA splicing, raising RBM20 as a potential therapeutic target.

Characterization of kinase gene expression and splicing profile in prostate cancer with RNA-Seq data.

BACKGROUND: Alternative splicing is a ubiquitous post-transcriptional regulation mechanism in most eukaryotic genes. Aberrant splicing isoforms and abnormal isoform ratios can contribute to cancer development. Kinase genes are key regulators of multiple cellular processes. Many kinases are found to be oncogenic and have been intensively investigated in the study of cancer and drugs. RNA-Seq provides a powerful technology for genome-wide study of alternative splicing in cancer besides the conventional gene expression profiling. But this potential has not been fully demonstrated yet. METHODS: We characterized the transcriptome profile of prostate cancer using RNA-Seq data from viewpoints of both differential expression and differential splicing, with an emphasis on kinase genes and their splicing variations. We built a pipeline to conduct differential expression and differential splicing analysis, followed by functional enrichment analysis. We performed kinase domain analysis to identify the functionally important candidate kinase gene in prostate cancer, and calculated the expression levels of isoforms to explore the function of isoform switching of kinase genes in prostate cancer. RESULTS: We identified distinct gene groups from differential expression and splicing analyses, which suggested that alternative splicing adds another level to gene expression regulation. Enriched GO terms of differentially expressed and spliced kinase genes were found to play different roles in regulation of cellular metabolism. Function analysis on differentially spliced kinase genes showed that differentially spliced exons of these genes are significantly enriched in protein kinase domains. Among them, we found that gene CDK5 has isoform switching between prostate cancer and benign tissues, which may affect cancer development by changing androgen receptor (AR) phosphorylation. The observation was validated in another RNA-Seq dataset of prostate cancer cell lines. CONCLUSIONS: Our work characterized the expression and splicing profiles of kinase genes in prostate cancer and proposed a hypothetical model on isoform switching of CDK5 and AR phosphorylation in prostate cancer. These findings bring new understanding to the role of alternatively spliced kinases in prostate cancer and also demonstrate the use of RNA-Seq data in studying alternative splicing in cancer.

Invention and Early History of Exon Skipping and Splice Modulation.

Since its discovery in 1977, much has been known about RNA splicing and how it plays a central role in human development, function, and, notably, disease. Defects in RNA splicing account for at least 10% of all genetic disorders, with the number expected to increase as more information is uncovered on the contribution of noncoding genomic regions to disease. Splice modulation through the use of antisense oligonucleotides (AOs) has emerged as a promising avenue for the treatment of these disorders. In fact, two splice-switching AOs have recently obtained approval from the US Food and Drug Administration: eteplirsen (Exondys 51) for Duchenne muscular dystrophy, and nusinersen (Spinraza) for spinal muscular atrophy. These work by exon skipping and exon inclusion, respectively. In this chapter, we discuss the early development of AO-based splice modulation therapy-its invention, first applications, and its evolution into the approach we are now familiar with. We give a more extensive history of exon skipping in particular, as it is the splice modulation approach given the most focus in this book.

Abiotic Stresses Modulate Landscape of Poplar Transcriptome via Alternative Splicing, Differential Intron Retention, and Isoform Ratio Switching.

Abiotic stresses affect plant physiology, development, growth, and alter pre-mRNA splicing. Western poplar is a model woody tree and a potential bioenergy feedstock. To investigate the extent of stress-regulated alternative splicing (AS), we conducted an in-depth survey of leaf, root, and stem xylem transcriptomes under drought, salt, or temperature stress. Analysis of approximately one billion of genome-aligned RNA-Seq reads from tissue- or stress-specific libraries revealed over fifteen millions of novel splice junctions. Transcript models supported by both RNA-Seq and single molecule isoform sequencing (Iso-Seq) data revealed a broad array of novel stress- and/or tissue-specific isoforms. Analysis of Iso-Seq data also resulted in the discovery of 15,087 novel transcribed regions of which 164 show AS. Our findings demonstrate that abiotic stresses profoundly perturb transcript isoform profiles and trigger widespread intron retention (IR) events. Stress treatments often increased or decreased retention of specific introns - a phenomenon described here as differential intron retention (DIR). Many differentially retained introns were regulated in a stress- and/or tissue-specific manner. A subset of transcripts harboring super stress-responsive DIR events showed persisting fluctuations in the degree of IR across all treatments and tissue types. To investigate coordinated dynamics of intron-containing transcripts in the study we quantified absolute copy number of isoforms of two conserved transcription factors (TFs) using Droplet Digital PCR. This case study suggests that stress treatments can be associated with coordinated switches in relative ratios between fully spliced and intron-retaining isoforms and may play a role in adjusting transcriptome to abiotic stresses.

RBM20, a potential target for treatment of cardiomyopathy via titin isoform switching.

Cardiomyopathy, also known as heart muscle disease, is an unfavorable condition leading to alterations in myocardial contraction and/or impaired ability of ventricular filling. The onset and development of cardiomyopathy have not currently been well defined. Titin is a giant multifunctional sarcomeric filament protein that provides passive stiffness to cardiomyocytes and has been implicated to play an important role in the origin and development of cardiomyopathy and heart failure. Titin-based passive stiffness can be mainly adjusted by isoform switching and post-translational modifications in the spring regions. Recently, genetic mutations of TTN have been identified that can also contribute to variable passive stiffness, though the detailed mechanisms remain unclear. In this review, we will discuss titin isoform switching as it relates to alternative splicing during development stages and differences between species and muscle types. We provide an update on the regulatory mechanisms of TTN splicing controlled by RBM20 and cover the roles of TTN splicing in adjusting the diastolic stiffness and systolic compliance of the healthy and the failing heart. Finally, this review attempts to provide future directions for RBM20 as a potential target for pharmacological intervention in cardiomyopathy and heart failure.

Mammalian Actins: Isoform-Specific Functions and Diseases.

Actin is the central building block of the actin cytoskeleton, a highly regulated filamentous network enabling dynamic processes of cells and simultaneously providing structure. Mammals have six actin isoforms that are very conserved and thus share common functions. Tissue-specific expression in part underlies their differential roles, but actin isoforms also coexist in various cell types and tissues, suggesting specific functions and preferential interaction partners. Gene deletion models, antibody-based staining patterns, gene silencing effects, and the occurrence of isoform-specific mutations in certain diseases have provided clues for specificity on the subcellular level and its consequences on the organism level. Yet, the differential actin isoform functions are still far from understood in detail. Biochemical studies on the different isoforms in pure form are just emerging, and investigations in cells have to deal with a complex and regulated system, including compensatory actin isoform expression.

CD44 single nucleotide polymorphism and isoform switching may predict gastric cancer recurrence.

BACKGROUND AND OBJECTIVES: The clinical implications of single nucleotide polymorphisms (SNPs) in CD44 remain unclear. This study examined the relationships of CD44 SNPs with clinicopathological parameters and prognosis in Japanese gastric cancer patients. METHODS: The CD44 SNPs were analyzed in 11 gastric cancer cell lines and 517 clinical specimens. The expression of CD44 standard (CD44s) and CD44 variant 9 isoform (CD44v9) transcripts were measured by quantitative real-time polymerase chain reaction. RESULTS: The CD44 rs187116 A/A, A/G, and G/G genotypes were present in 10.3%, 45.1%, and 44.7% of patients, respectively. The presence of CD44 rs187116 A/G or G/G genotypes was significantly associated with positive peritoneal washing cytology (P = 0.012). Disease-free survival of patients with these genotypes was significantly worse than in those with the A/A genotype (P = 0.039). Multivariate analysis showed that the CD44 rs187116 was independently prognostic of disease-free survival (P = 0.047). The CD44s/CD44v9 ratio was significantly lower in patients with the CD44 rs187116 A/A genotype than in those with the A/G (P = 0.046) and G/G (P = 0.047) genotypes. CONCLUSIONS: The CD44 rs187116 genotype could predict disease recurrence in Japanese gastric cancer patients, and the SNP was associated with CD44 isoform switching.