Competing upstream 5' splice sites enhance the rate of proximal splicing

Identified five variants in CFTR gene that alter RNA splicing by minigene assay

Background

Cystic fibrosis (CF) is a common monogenic multisystem disease caused primarily by variants in the CFTR gene. Emerging evidence suggests that some variants, which are described as missense, synonymous or nonsense variants in the literature or databases, may be deleterious by affecting the pre-mRNA splicing process.

Methods

We analyzed 27 exonic variants in the CFTR gene utilizing bioinformatics tools and identified candidate variants that could lead to splicing changes through minigene assays. Ultimately, we selected eight candidate variants to assess their effects on pre-mRNA splicing. The numbering of DNA variants is based on the complementary DNA (cDNA)sequence of CFTR (Ref Seq NM_000492.4).

Results

This study assessed the impact of CFTR variants on exon splicing by combining predictive bioinformatics tools with minigene assays. Among the eight candidate single nucleotide alterations, five variants (c.488A>T,c.1117G>T, c.1209G>T, c.3239A>G and c.3367G>C) were identified as causing exon skipping.

Conclusion

Our study employed a minigene system, which offers great flexibility for assessing aberrant splicing patterns when patient mRNA samples are not accessible, to investigate the effects of exonic variants on pre-mRNA splicing. Our experimental outcomes highlight the importance of analyzing exonic variations at the mRNA level.

Identification of seven variants in the col4a1 gene that alter RNA splicing by minigene assay

Type IV collagen is an integral component of basement membranes. Mutations in COL4A1, one of the key genes encoding Type IV collagen, can result in a variety of diseases. It is clear that a significant proportion of mutations that affect splicing can cause disease directly or contribute to the susceptibility or severity of disease. Here, we analyzed exonic mutations and intronic mutations described in the COL4A1 gene using bioinformatics programs and identified candidate mutations that may alter the normal splicing pattern through a minigene system. We identified seven variants that induce splicing alterations by disrupting normal splice sites, creating new ones, or altering splice regulatory elements. These mutations are predicted to impact protein function. Our results help in the correct molecular characterization of variants in COL4A1 and may help develop more personalized treatment options.

U1 snRNA interactions with deep intronic sequences regulate splicing of multiple exons of spinal muscular atrophy genes

Introduction

The U1 small nuclear RNA (snRNA) forms ribonucleoprotein particles (RNPs) such as U1 snRNP and U1-TAF15 snRNP. U1 snRNP is one of the most studied RNPs due to its critical role in pre-mRNA splicing in defining the 5' splice site (5'ss) of every exon through direct interactions with sequences at exon/intron junctions. Recent reports support the role of U1 snRNP in all steps of transcription, namely initiation, elongation, and termination. Functions of U1-TAF15 snRNP are less understood, though it associates with the transcription machinery and may modulate pre-mRNA splicing by interacting with the 5'ss and/or 5'ss-like sequences within the pre-mRNA. An anti-U1 antisense oligonucleotide (ASO) that sequesters the 5' end of U1 snRNA inhibits the functions of U1 snRNP, including transcription and splicing. However, it is not known if the inhibition of U1 snRNP influences post-transcriptional regulation of pre-mRNA splicing through deep intronic sequences.

Methods

We examined the effect of an anti-U1 ASO that sequesters the 5' end of U1 snRNA on transcription and splicing of all internal exons of the spinal muscular atrophy (SMA) genes, SMN1 and SMN2. Our study was enabled by the employment of a multi-exon-skipping detection assay (MESDA) that discriminates against prematurely terminated transcripts. We employed an SMN2 super minigene to determine if anti-U1 ASO differently affects splicing in the context of truncated introns.

Results

We observed substantial skipping of multiple internal exons of SMN1 and SMN2 triggered by anti-U1 treatment. Suggesting a role for U1 snRNP in interacting with deep intronic sequences, early exons of the SMN2 super minigene with truncated introns were resistant to anti-U1 induced skipping. Consistently, overexpression of engineered U1 snRNAs targeting the 5'ss of early SMN1 and SMN2 exons did not prevent exon skipping caused by anti-U1 treatment.

Discussion

Our results uncover a unique role of the U1 snRNA-associated RNPs in splicing regulation executed through deep intronic sequences. Findings are significant for developing novel therapies for SMA based on deep intronic targets.

SRSF1 acts as an IFN-I-regulated cellular dependency factor decisively affecting HIV-1 post-integration steps

Efficient HIV-1 replication depends on balanced levels of host cell components including cellular splicing factors as the family of serine/arginine-rich splicing factors (SRSF, 1-10). Type I interferons (IFN-I) play a crucial role in the innate immunity against HIV-1 by inducing the expression of IFN-stimulated genes (ISGs) including potent host restriction factors. The less well known IFN-repressed genes (IRepGs) might additionally affect viral replication by downregulating host dependency factors that are essential for the viral life cycle; however, so far, the knowledge about IRepGs involved in HIV-1 infection is very limited. In this work, we could demonstrate that HIV-1 infection and the associated ISG induction correlated with low SRSF1 levels in intestinal lamina propria mononuclear cells (LPMCs) and peripheral blood mononuclear cells (PBMCs) during acute and chronic HIV-1 infection. In HIV-1-susceptible cell lines as well as primary monocyte-derived macrophages (MDMs), expression levels of SRSF1 were transiently repressed upon treatment with specific IFNα subtypes in vitro. Mechanically, 4sU labeling of newly transcribed mRNAs revealed that IFN-mediated SRSF1 repression is regulated on early RNA level. SRSF1 knockdown led to an increase in total viral RNA levels, but the relative proportion of the HIV-1 viral infectivity factor (Vif) coding transcripts, which is essential to counteract APOBEC3G-mediated host restriction, was significantly reduced. In the presence of high APOBEC3G levels, however, increased LTR activity upon SRSF1 knockdown facilitated the overall replication, despite decreased vif mRNA levels. In contrast, SRSF1 overexpression significantly impaired HIV-1 post-integration steps including LTR transcription, alternative splice site usage, and virus particle production. Since balanced SRSF1 levels are crucial for efficient viral replication, our data highlight the so far undescribed role of SRSF1 acting as an IFN-modulated cellular dependency factor decisively regulating HIV-1 post-integration steps.

Intronic Polyadenylation in Acquired Cancer Drug Resistance Circumvented by Utilizing CRISPR/Cas9 with Homology-Directed Repair: The Tale of Human DNA Topoisomerase IIα

Simple Summary

DNA topoisomerase IIα (170 kDa, TOP2α/170) resolves nucleic acid topological entanglements by generating transient double-strand DNA breaks. TOP2α inhibitors/poisons stabilize TOP2α-DNA covalent complexes resulting in persistent DNA damage and are frequently utilized to treat a variety of cancers. Acquired resistance to these chemotherapeutic agents is often associated with decreased TOP2α/170 expression levels. Studies have demonstrated that a reduction in TOP2α/170 results from a type of alternative polyadenylation designated intronic polyadenylation (IPA). As a consequence of IPA, variant TOP2α mRNA transcripts have been characterized that have resulted in the translation of C-terminal truncated TOP2α isoforms with altered biological activities. In this paper, an example is discussed where circumvention of acquired TOP2α-mediated drug resistance was achieved by utilizing CRISPR/Cas9 specific gene editing of an exon/intron boundary through homology directed repair (HDR) to reduce TOP2α IPA. These results illustrate the therapeutic potential of CRISPR/Cas9/HDR to impact drug resistance associated with aberrant IPA.

Abstract

Intronic polyadenylation (IPA) plays a critical role in malignant transformation, development, progression, and cancer chemoresistance by contributing to transcriptome/proteome alterations. DNA topoisomerase IIα (170 kDa, TOP2α/170) is an established clinical target for anticancer agents whose efficacy is compromised by drug resistance often associated with a reduction of nuclear TOP2α/170 levels. In leukemia cell lines with acquired resistance to TOP2α-targeted drugs and reduced TOP2α/170 expression, variant TOP2α mRNA transcripts have been reported due to IPA that resulted in the translation of C-terminal truncated isoforms with altered nuclear-cytoplasmic distribution or heterodimerization with wild-type TOP2α/170. This review provides an overview of the various mechanisms regulating pre-mRNA processing and alternative polyadenylation, as well as the utilization of CRISPR/Cas9 specific gene editing through homology directed repair (HDR) to decrease IPA when splice sites are intrinsically weak or potentially mutated. The specific case of TOP2α exon 19/intron 19 splice site editing is discussed in etoposide-resistant human leukemia K562 cells as a tractable strategy to circumvent acquired TOP2α-mediated drug resistance. This example supports the importance of aberrant IPA in acquired drug resistance to TOP2α-targeted drugs. In addition, these results demonstrate the therapeutic potential of CRISPR/Cas9/HDR to impact drug resistance associated with aberrant splicing/polyadenylation.

Use of CRISPR/Cas9 with homology-directed repair to silence the human topoisomerase IIα intron-19 5’ splice site: Generation of etoposide resistance in human leukemia K562 cells

DNA Topoisomerase IIα (TOP2α/170) is an enzyme essential for proliferating cells. For rapidly multiplying malignancies, this has made TOP2α/170 an important target for etoposide and other clinically active anticancer drugs. Efficacy of these agents is often limited by chemoresistance related to alterations in TOP2α/170 expression levels. Our laboratory recently demonstrated reduced levels of TOP2α/170 and overexpression of a C-terminal truncated 90-kDa isoform, TOP2α/90, due to intronic polyadenylation (IPA; within intron 19) in an acquired etoposide-resistant K562 clonal cell line, K/VP.5. We previously reported that this isoform heterodimerized with TOP2α/170 and was a determinant of acquired resistance to etoposide. Optimization of the weak TOP2α exon 19/intron 19 5′ splice site in drug-resistant K/VP.5 cells by gene-editing restored TOP2α/170 levels, diminished TOP2α/90 expression, and circumvented drug resistance. Conversely, in the present study, silencing of the exon 19/intron 19 5′ splice site in parental K562 cells by CRISPR/Cas9 with homology-directed repair (HDR), and thereby forcing intron 19 retention, was used to induce resistance by disrupting normal RNA processing (i.e., gene knockout), and to further evaluate the role of TOP2α/170 and TOP2α/90 isoforms as resistance determinants. Gene-edited clones were identified by quantitative polymerase chain reaction (qPCR) and verified by Sanger sequencing. TOP2α/170 mRNA/protein expression levels were attenuated in the TOP2α gene-edited clones which resulted in resistance to etoposide as assessed by reduced etoposide-induced DNA damage (γH2AX, Comet assays) and growth inhibition. RNA-seq and qPCR studies suggested that intron 19 retention leads to decreased TOP2α/170 expression by degradation of the TOP2α edited mRNA transcripts. Forced expression of TOP2α/90 in the gene-edited K562 cells further decreased etoposide-induced DNA damage in support of a dominant negative role for this truncated isoform. Together results support the important role of both TOP2α/170 and TOP2α/90 as determinants of sensitivity/resistance to TOP2α-targeting agents.

Relative strength of 5’ splice-site strength defines functions of SRSF2 and SRSF6 in alternative splicing of Bcl-x pre-mRNA

Bcl-x, a member of the Bcl-2 family, plays a key role in apoptosis. Alternative splicing of Bcl-x pre-mRNA through alternative 5’ splice-site selection produces an anti-apoptotic mRNA isoform that includes exon 2b and a pro-apoptotic Bcl-x mRNA isoform that excludes exon 2b. Here we used Bcl-x minigene and identified SRSF2 and SRSF6 as two regulatory factors of 5’ splice-site selection of Bcl-x pre-mRNA. We selected binding clusters closer to 5’ splice-sites from multiple potential binding sites of SRSF2 and SRSF6 to perform loss of functions analysis through site-directed mutagenesis. Our results demonstrated that these mutations did not abolish regulatory functions of SRSF2 or SRSF6, indicating that a single binding motif or a cluster was not a functional target of these proteins in Bcl-x pre-mRNA splicing. Random deletion mutagenesis did not disrupt the role of SRSF2 and SRSF6. Importantly, mutagenesis of 5’ splice-site to a conserved or a weaker score demonstrated that the weaker strength of the target 5’ splice-site or higher strength of the other 5’ splice-site strength limited the role of SRSF2 and SRSF6 in 5’ splice-site activation.

CRISPR/Cas9 Genome Editing of the Human Topoisomerase IIα Intron 19 5' Splice Site Circumvents Etoposide Resistance in Human Leukemia K562 Cells

An essential function of DNA topoisomerase IIα (TOP2α; 170 kDa, TOP2α/170) is to resolve DNA topologic entanglements during chromosome disjunction by introducing transient DNA double-stranded breaks. TOP2α/170 is an important target for DNA damage-stabilizing anticancer drugs, whose clinical efficacy is compromised by drug resistance often associated with decreased TOP2α/170 expression. We recently demonstrated that an etoposide-resistant K562 clonal subline, K/VP.5, with reduced levels of TOP2α/170, expresses high levels of a novel C-terminal truncated TOP2α isoform (90 kDa, TOP2α/90). TOP2α/90, the translation product of a TOP2α mRNA that retains a processed intron 19 (I19), heterodimerizes with TOP2α/170 and is a resistance determinant through a dominant-negative effect on drug activity. We hypothesized that genome editing to enhance I19 removal would provide a tractable strategy to circumvent acquired TOP2α-mediated drug resistance. To enhance I19 removal in K/VP.5 cells, CRISPR/Cas9 was used to make changes (GAG//GTAA AC →GAG//GTAA GT ) in the TOP2α gene's suboptimal exon 19/intron 19 5' splice site (E19/I19 5' SS). Gene-edited clones were identified by quantitative polymerase chain reaction and verified by sequencing. Characterization of a clone with all TOP2α alleles edited revealed improved I19 removal, decreased TOP2α/90 mRNA/protein, and increased TOP2α/170 mRNA/protein. Sensitivity to etoposide-induced DNA damage (γH2AX, Comet assays) and growth inhibition was restored to levels comparable to those in parental K562 cells. Together, the results indicate that our gene-editing strategy for optimizing the TOP2α E19/I19 5' SS in K/VP.5 cells circumvents resistance to etoposide and other TOP2α-targeted drugs. SIGNIFICANCE STATEMENT: Results presented here indicate that CRISPR/Cas9 gene editing of a suboptimal exon 19/intron 19 5' splice site in the DNA topoisomerase IIα (TOP2α) gene results in circumvention of acquired drug resistance to etoposide and other TOP2α-targeted drugs in a clonal K562 cell line by enhancing removal of intron 19 and thereby decreasing formation of a truncated TOP2α 90 kDa isoform and increasing expression of full-length TOP2α 170 kDa in these resistant cells. Results demonstrate the importance of RNA processing in acquired drug resistance to TOP2α-targeted drugs.

Splicing at the phase-separated nuclear speckle interface: a model

Phase-separated membraneless bodies play important roles in nucleic acid biology. While current models for the roles of phase separation largely focus on the compartmentalization of constituent proteins, we reason that other properties of phase separation may play functional roles. Specifically, we propose that interfaces of phase-separated membraneless bodies could have functional roles in spatially organizing biochemical reactions. Here we propose such a model for the nuclear speckle, a membraneless body implicated in RNA splicing. In our model, sequence-dependent RNA positioning along the nuclear speckle interface coordinates RNA splicing. Our model asserts that exons are preferentially sequestered into nuclear speckles through binding by SR proteins, while introns are excluded through binding by nucleoplasmic hnRNP proteins. As a result, splice sites at exon-intron boundaries are preferentially positioned at nuclear speckle interfaces. This positioning exposes splice sites to interface-localized spliceosomes, enabling the subsequent splicing reaction. Our model provides a simple mechanism that seamlessly explains much of the complex logic of splicing. This logic includes experimental results such as the antagonistic duality between splicing factors, the position dependence of splicing sequence motifs, and the collective contribution of many motifs to splicing decisions. Similar functional roles for phase-separated interfaces may exist for other membraneless bodies.

Interplay of primary sequence, position and secondary RNA structure determines alternative splicing of LMNA in a pre-mature aging syndrome

Aberrant splicing in exon 11 of the LMNA gene causes the premature aging disorder Hutchinson-Gilford Progeria Syndrome. A de novo C1824T mutation activates an internal alternative 5' splice site, resulting in formation of the disease-causing progerin protein. The underlying mechanism for this 5' splice site selection is unknown. Here, we have applied a combination of targeted mutational analysis in a cell-based system and structural mapping by SHAPE-MaP to comprehensively probe the contributions of primary sequence, secondary RNA structure and linear splice site position in determining in vivo mechanisms of splice site choice in LMNA. While splice site choice is in part defined by sequence complementarity to U1 snRNA, we identify RNA secondary structural elements near the alternative 5' splice sites and show that splice site choice is significantly influenced by the structural context of the available splice sites. Furthermore, relative positioning of the competing sites within the primary sequence of the pre-mRNA is a predictor of 5' splice site usage, with the distal position favored over the proximal, regardless of sequence composition. Together, these results demonstrate that 5' splice site selection in LMNA is determined by an intricate interplay among RNA sequence, secondary structure and splice site position.

Dissecting splicing decisions and cell-to-cell variability with designed sequence libraries

Most human genes are alternatively spliced, allowing for a large expansion of the proteome. The multitude of regulatory inputs to splicing limits the potential to infer general principles from investigating native sequences. Here, we create a rationally designed library of >32,000 splicing events to dissect the complexity of splicing regulation through systematic sequence alterations. Measuring RNA and protein splice isoforms allows us to investigate both cause and effect of splicing decisions, quantify diverse regulatory inputs and accurately predict (R² = 0.73–0.85) isoform ratios from sequence and secondary structure. By profiling individual cells, we measure the cell-to-cell variability of splicing decisions and show that it can be encoded in the DNA and influenced by regulatory inputs, opening the door for a novel, single-cell perspective on splicing regulation.

Alternative splicing is regulated by multiple mechanisms. Here the authors employed designed splice site libraries and massively parallel reporter assays to dissect the regulatory complexity and cell-to-cell variability of splicing decisions and to build accurate predictive models.

Combinatorial regulation of alternative splicing

The generation of protein coding mRNAs from pre-mRNA is a fundamental biological process that is required for gene expression. Alternative pre-mRNA splicing is responsible for much of the transcriptomic and proteomic diversity observed in higher order eukaryotes. Aberrations that disrupt regular alternative splicing patterns are known to cause human diseases, including various cancers. Alternative splicing is a combinatorial process, meaning many factors affect which two splice sites are ligated together. The features that dictate exon inclusion are comprised of splice site strength, intron-exon architecture, RNA secondary structure, splicing regulatory elements, promoter use and transcription speed by RNA polymerase and the presence of post-transcriptional nucleotide modifications. A comprehensive view of all of the factors that influence alternative splicing decisions is necessary to predict splicing outcomes and to understand the molecular basis of disease. This article is part of a Special Issue entitled: RNA structure and splicing regulation edited by Francisco Baralle, Ravindra Singh and Stefan Stamm.

High-throughput analysis revealed mutations' diverging effects on SMN1 exon 7 splicing

Splicing-affecting mutations can disrupt gene function by altering the transcript assembly. To ascertain splicing dysregulation principles, we modified a minigene assay for the parallel high-throughput evaluation of different mutations by next-generation sequencing. In our model system, all exonic and six intronic positions of the SMN1 gene's exon 7 were mutated to all possible nucleotide variants, which amounted to 180 unique single-nucleotide mutants and 470 double mutants. The mutations resulted in a wide range of splicing aberrations. Exonic splicing-affecting mutations resulted either in substantial exon skipping, supposedly driven by predicted exonic splicing silencer or cryptic donor splice site (5'ss) and de novo 5'ss strengthening and use. On the other hand, a single disruption of exonic splicing enhancer was not sufficient to cause major exon skipping, suggesting these elements can be substituted during exon recognition. While disrupting the acceptor splice site led only to exon skipping, some 5'ss mutations potentiated the use of three different cryptic 5'ss. Generally, single mutations supporting cryptic 5'ss use displayed better pre-mRNA/U1 snRNA duplex stability and increased splicing regulatory element strength across the original 5'ss. Analyzing double mutants supported the predominating splicing regulatory elements' effect, but U1 snRNA binding could contribute to the global balance of splicing isoforms. Based on these findings, we suggest that creating a new splicing enhancer across the mutated 5'ss can be one of the main factors driving cryptic 5'ss use.

Context matters: Regulation of splice donor usage

Elaborate research on splicing, starting in the late seventies, evolved from the discovery that 5' splice sites are recognized by their complementarity to U1 snRNA towards the realization that RNA duplex formation cannot be the sole basis for 5'ss selection. Rather, their recognition is highly influenced by a number of context factors including transcript architecture as well as splicing regulatory elements (SREs) in the splice site neighborhood. In particular, proximal binding of splicing regulatory proteins highly influences splicing outcome. The importance of SRE integrity especially becomes evident in the light of human pathogenic mutations where single nucleotide changes in SREs can severely affect the resulting transcripts. Bioinformatics tools nowadays greatly assist in the computational evaluation of 5'ss, their neighborhood and the impact of pathogenic mutations. Although predictions are already quite robust, computational evaluation of the splicing regulatory landscape still faces challenges to increase future reliability. This article is part of a Special Issue entitled: RNA structure and splicing regulation edited by Francisco Baralle, Ravindra Singh and Stefan Stamm.

Global donor and acceptor splicing site kinetics in human cells

RNA splicing is an essential part of eukaryotic gene expression. Although the mechanism of splicing has been extensively studied in vitro, in vivo kinetics for the two-step splicing reaction remain poorly understood. Here, we combine transient transcriptome sequencing (TT-seq) and mathematical modeling to quantify RNA metabolic rates at donor and acceptor splice sites across the human genome. Splicing occurs in the range of minutes and is limited by the speed of RNA polymerase elongation. Splicing kinetics strongly depends on the position and nature of nucleotides flanking splice sites, and on structural interactions between unspliced RNA and small nuclear RNAs in spliceosomal intermediates. Finally, we introduce the 'yield' of splicing as the efficiency of converting unspliced to spliced RNA and show that it is highest for mRNAs and independent of splicing kinetics. These results lead to quantitative models describing how splicing rates and yield are encoded in the human genome.

HMGA1a Induces Alternative Splicing of the Estrogen Receptor-αlpha Gene by Trapping U1 snRNP to an Upstream Pseudo-5' Splice Site

Objectives: The high-mobility group A protein 1a (HMGA1a) protein is known as a transcription factor that binds to DNA, but recent studies have shown it exerts novel functions through RNA-binding. We were prompted to decipher the mechanism of HMGA1a-induced alternative splicing of the estrogen receptor alpha (ERα) that we recently reported would alter tamoxifen sensitivity in MCF-7 TAMR1 cells.

Methods: Endogenous expression of full length ERα66 and its isoform ERα46 were evaluated in MCF-7 breast cancer cells by transient expression of HMGA1a and an RNA decoy (2′-O-methylated RNA of the HMGA1a RNA-binding site) that binds to HMGA1a. RNA-binding of HMGA1a was checked by RNA-EMSA. In vitro splicing assay was performed to check the direct involvement of HMGA1a in splicing regulation. RNA-EMSA assay in the presence of purified U1 snRNP was performed with psoralen UV crosslinking to check complex formation of HMGA1a-U1 snRNP at the upstream pseudo-5′ splice site of exon 1.

Results: HMGA1a induced exon skipping of a shortened exon 1 of ERα in in vitro splicing assays that was blocked by the HMGA1a RNA decoy and sequence-specific RNA-binding was confirmed by RNA-EMSA. RNA-EMSA combined with psoralen UV crosslinking showed that HMGA1a trapped purified U1 snRNP at the upstream pseudo-5′ splice site.

Conclusions: Regulation of ERα alternative splicing by an HMGA1a-trapped U1 snRNP complex at the upstream 5′ splice site of exon 1 offers novel insight on 5′ splice site regulation by U1 snRNP as well as a promising target in breast cancer therapy where alternative splicing of ERα is involved.

The kinetics of pre-mRNA splicing in the Drosophila genome and the influence of gene architecture

Production of most eukaryotic mRNAs requires splicing of introns from pre-mRNA. The splicing reaction requires definition of splice sites, which are initially recognized in either intron-spanning (‘intron definition’) or exon-spanning (‘exon definition’) pairs. To understand how exon and intron length and splice site recognition mode impact splicing, we measured splicing rates genome-wide in Drosophila, using metabolic labeling/RNA sequencing and new mathematical models to estimate rates. We found that the modal intron length range of 60–70 nt represents a local maximum of splicing rates, but that much longer exon-defined introns are spliced even faster and more accurately. We observed unexpectedly low variation in splicing rates across introns in the same gene, suggesting the presence of gene-level influences, and we identified multiple gene level variables associated with splicing rate. Together our data suggest that developmental and stress response genes may have preferentially evolved exon definition in order to enhance the rate or accuracy of splicing.

Targeting Splicing in the Treatment of Myelodysplastic Syndromes and Other Myeloid Neoplasms

Genome sequencing of primary cells from patients with myelodysplastic syndromes (MDS) led to the identification of recurrent heterozygous mutations in gene encoding components of the spliceosome, the cellular machinery which processes pre-messenger RNA (mRNA) to mature mRNA during gene transcription. Splicing mutations are mutually exclusive with one another and collectively represent the most common mutation class in MDS, occurring in approximately 60 % of patients overall and more than 80 % of those with ring sideroblasts. Evidence from animal models suggests that homozygous splicing mutations are lethal, and that in heterozygously mutated models, any further disruption of splicing triggers apoptosis and cell death. MDS cells with spliceosome mutations are thus uniquely vulnerable to therapies targeting splicing, which may be tolerated by healthy cells. The spliceosome is emerging as a novel therapeutic target in MDS and related myeloid neoplasms, with the first clinical trial of a splicing modulator opening in 2016.

Activation of the intronic cryptic 5' splice site depends on its distance to the upstream cassette exon

Splice site selection is a key step that determines the mRNA isoforms generated from a single transcript. The large diversity in splice site sequences emphasizes the plasticity of splice site recognition and selection. In this report, a cell-based reporter system using a SMN1/2 cassette exon was applied to study the roles governing the activation of a cryptic 5'SS from the intron 4 of the CT/CGRP gene. We found that the cryptic site was activated when placed within 124nt downstream the cassette exon, and the level of activation was negatively correlated with its distance from the exon. In addition, activation was not affected by PTB but was eliminated by an insertion extending the exon length. Activated cryptic 5'SSs in intron or exon could override the original alternative 5'SS, obeying the U1 base-pairing rule. These results suggest that the exon length itself could represent a factor in determining the splice site selection.

Suppression of 5' splice-sites through multiple exonic motifs by hnRNP L

Selection of 5' splice-sites (5'SS) in alternative splicing plays an important role in gene regulation. Although regulatory mechanisms of heterogeneous nuclear ribonucleoprotein L (hnRNP L), a well-known splicing regulatory protein, have been studied in a substantial level, its role in 5'SS selection is not thoroughly defined. By using a KLF6 pre-mRNA alternative splicing model, we demonstrate in this report that hnRNP L inhibits proximal 5'SS but promotes two consecutive distal 5'SS splicing, antagonizing SRSF1 roles in KLF6 pre-mRNA splicing. In addition, three consecutive CA-rich sequences in a CA cassette immediately upstream of the proximal 5'SS are all required for hnRNP L functions. Importantly, the CA-cassette locations on the proximal exon do not affect hnRNP L roles. We further show that the proximal 5'SS but not the two distal 5'SSs are essential for hnRNP L activities. Notably, in a Bcl-x pre-mRNA model that contains two alternative 5'SS but includes CA-rich elements at distal exon, we demonstrate that hnRNP L also suppresses nearby 5'SS activation. Taken together, we conclude that hnRNP L suppresses 5'SS selection through multiple exonic motifs.

Kinetic analysis of in vitro pre-mRNA splicing in HeLa nuclear extract

Kinetic analysis of in vitro splicing is a valuable technique for understanding splicing regulation. It allows the determination of specific contributions from functional elements for the efficient removal of introns. This chapter will describe the rationale and approach employed to use kinetic analysis to evaluate an in vitro splicing reaction using radiolabeled pre-mRNA incubated in splicing-competent HeLa nuclear extract (NE).

In silico to in vivo splicing analysis using splicing code models

With the growing appreciation of RNA splicing's role in gene regulation, development, and disease, researchers from diverse fields find themselves investigating exons of interest. Commonly, researchers are interested in knowing if an exon is alternatively spliced, if it is differentially included in specific tissues or in developmental stages, and what regulatory elements control its inclusion. An important step towards the ability to perform such analysis in silico was made with the development of computational splicing code models. Aimed as a practical how-to guide, we demonstrate how researchers can now use these code models to analyze a gene of interest, focusing on Bin1 as a case study. Bridging integrator 1 (BIN1) is a nucleocytoplasmic adaptor protein known to be functionally regulated through alternative splicing in a tissue-specific manner. Specific Bin1 isoforms have been associated with muscular diseases and cancers, making the study of its splicing regulation of wide interest. Using AVISPA, a recently released web tool based on splicing code models, we show that many Bin1 tissue-dependent isoforms are correctly predicted, along with many of its known regulators. We review the best practices and constraints of using the tool, demonstrate how AVISPA is used to generate high confidence novel regulatory hypotheses, and experimentally validate predicted regulators of Bin1 alternative splicing.

Misregulation of pre-mRNA alternative splicing in cancer

Alternative splicing of mRNA precursors enables one gene to produce multiple protein isoforms with differing functions. Under normal conditions, this mechanism is tightly regulated in order for the human genome to generate proteomic diversity sufficient for the functional requirements of complex tissues. When deregulated, however, cancer cells take advantage of this mechanism to produce aberrant proteins with added, deleted, or altered functional domains that contribute to tumorigenesis. Here, we discuss aspects of alternative splicing misregulation in cancer, focusing on splicing events affected by deregulation of regulatory splicing factors and also recent studies identifying mutated components of the splicing machinery.

SIGNIFICANCE

An increasing body of evidence indicates that aberrant splicing of mRNA precursors leads to production of aberrant proteins that contribute to tumorigenesis. Recent studies show that alterations in cellular concentrations of regulatory splicing factors and mutations in components of the core splicing machinery provide major mechanisms of misregulation of mRNA splicing in cancer. A better understanding of this misregulation will potentially reveal a group of novel drug targets for therapeutic intervention.

Splicing kinetics and transcript release from the chromatin compartment limit the rate of Lipid A-induced gene expression

The expression of eukaryotic mRNAs is achieved though an intricate series of molecular processes that provide many steps for regulating the production of a final gene product. However, the relationships between individual steps in mRNA biosynthesis and the rates at which they occur are poorly understood. By applying RNA-seq to chromatin-associated and soluble nucleoplasmic fractions of RNA from Lipid A-stimulated macrophages, we examined the timing of exon ligation and transcript release from chromatin relative to the induction of transcription. We find that for a subset of genes in the Lipid A response, the ligation of certain exon pairs is delayed relative to the synthesis of the complete transcript. In contrast, 3' end cleavage and polyadenylation occur rapidly once transcription extends through the cleavage site. Our data indicate that these transcripts with delayed splicing are not released from the chromatin fraction until all the introns have been excised. These unusual kinetics result in a chromatin-associated pool of completely transcribed and 3'-processed transcripts that are not yet fully spliced. We also find that long introns containing repressed exons that will be excluded from the final mRNA are excised particularly slowly relative to other introns in a transcript. These results indicate that the kinetics of splicing and transcript release contribute to the timing of expression for multiple genes of the inflammatory response.

Pick one, but be quick: 5' splice sites and the problems of too many choices

Splice site selection is fundamental to pre-mRNA splicing and the expansion of genomic coding potential. 5' Splice sites (5'ss) are the critical elements at the 5' end of introns and are extremely diverse, as thousands of different sequences act as bona fide 5'ss in the human transcriptome. Most 5'ss are recognized by base-pairing with the 5' end of the U1 small nuclear RNA (snRNA). Here we review the history of research on 5'ss selection, highlighting the difficulties of establishing how base-pairing strength determines splicing outcomes. We also discuss recent work demonstrating that U1 snRNA:5'ss helices can accommodate noncanonical registers such as bulged duplexes. In addition, we describe the mechanisms by which other snRNAs, regulatory proteins, splicing enhancers, and the relative positions of alternative 5'ss contribute to selection. Moreover, we discuss mechanisms by which the recognition of numerous candidate 5'ss might lead to selection of a single 5'ss and propose that protein complexes propagate along the exon, thereby changing its physical behavior so as to affect 5'ss selection.

Antisense oligonucleotide mediated therapy of spinal muscular atrophy

Spinal muscular atrophy (SMA) is the leading genetic cause of infant mortality. SMA results from deletions or mutations of survival motor neuron 1 (SMN1), an essential gene. SMN2, a nearly identical copy, can compensate for SMN1 loss if SMN2 exon 7 skipping is prevented. Among the many cis-elements involved in the splicing regulation of SMN exon 7, intronic splicing silencer N1 (ISS-N1) has emerged as the most effective target for an antisense oligonucleotide (ASO)-mediated splicing correction of SMN2 exon 7. Blocking of ISS-N1 by an ASO has been shown to fully restore SMN2 exon 7 inclusion in SMA patient cells as well as in vivo. Here we review how ISS-N1 targeting ASOs that use different chemistries respond differently in the various SMA mouse models. We also compare other ASO-based strategies for therapeutic splicing correction in SMA. Given that substantial progress on ASO-based strategies to promote SMN2 exon 7 inclusion in SMA has been made, and that similar approaches in a growing number of genetic diseases are possible, this report has wide implications.

Position-dependent splicing activation and repression by SR and hnRNP proteins rely on common mechanisms

Alternative splicing is regulated by splicing factors that modulate splice site selection. In some cases, however, splicing factors show antagonistic activities by either activating or repressing splicing. Here, we show that these opposing outcomes are based on their binding location relative to regulated 5' splice sites. SR proteins enhance splicing only when they are recruited to the exon. However, they interfere with splicing by simply relocating them to the opposite intronic side of the splice site. hnRNP splicing factors display analogous opposing activities, but in a reversed position dependence. Activation by SR or hnRNP proteins increases splice site recognition at the earliest steps of exon definition, whereas splicing repression promotes the assembly of nonproductive complexes that arrest spliceosome assembly prior to splice site pairing. Thus, SR and hnRNP splicing factors exploit similar mechanisms to positively or negatively influence splice site selection.

Global profiling of alternative splicing events and gene expression regulated by hnRNPH/F

In this study, we have investigated the global impact of heterogeneous nuclear Ribonuclear Protein (hnRNP) H/F-mediated regulation of splicing events and gene expression in oligodendrocytes. We have performed a genome-wide transcriptomic analysis at the gene and exon levels in Oli-neu cells treated with siRNA that targets hnRNPH/F compared to untreated cells using Affymetrix Exon Array. Gene expression levels and regulated exons were identified with the GenoSplice EASANA algorithm. Bioinformatics analyses were performed to determine the structural properties of G tracts that correlate with the function of hnRNPH/F as enhancers vs. repressors of exon inclusion. Different types of alternatively spliced events are regulated by hnRNPH/F. Intronic G tracts density, length and proximity to the 5′ splice site correlate with the hnRNPH/F enhancer function. Additionally, 6% of genes are differently expressed upon knock down of hnRNPH/F. Genes that regulate the transition of oligodendrocyte progenitor cells to oligodendrocytes are differentially expressed in hnRNPH/F depleted Oli-neu cells, resulting in a decrease of negative regulators and an increase of differentiation-inducing regulators. The changes were confirmed in developing oligodendrocytes in vivo. This is the first genome wide analysis of splicing events and gene expression regulated by hnRNPH/F in oligodendrocytes and the first report that hnRNPH/F regulate genes that are involved in the transition from oligodendrocyte progenitor cells to oligodendrocytes.

Intron definition and a branch site adenosine at nt 385 control RNA splicing of HPV16 E6*I and E7 expression

HPV16 E6 and E7, two viral oncogenes, are expressed from a single bicistronic pre-mRNA. In this report, we provide the evidence that the bicistronic pre-mRNA intron 1 contains three 5' splice sites (5' ss) and three 3' splice sites (3' ss) normally used in HPV16(+) cervical cancer and its derived cell lines. The choice of two novel alternative 5' ss (nt 221 5' ss and nt 191 5' ss) produces two novel isoforms of E6E7 mRNAs (E6*V and E6*VI). The nt 226 5' ss and nt 409 3' ss is preferentially selected over the other splice sites crossing over the intron to excise a minimal length of the intron in RNA splicing. We identified AACAAAC as the preferred branch point sequence (BPS) and an adenosine at nt 385 (underlined) in the BPS as a branch site to dictate the selection of the nt 409 3' ss for E6*I splicing and E7 expression. Introduction of point mutations into the mapped BPS led to reduced U2 binding to the BPS and thereby inhibition of the second step of E6E7 splicing at the nt 409 3' ss. Importantly, the E6E7 bicistronic RNA with a mutant BPS and inefficient splicing makes little or no E7 and the resulted E6 with mutations of (91)QYNK(94) to (91)PSFW(94) displays attenuate activity on p53 degradation. Together, our data provide structural basis of the E6E7 intron 1 for better understanding of how viral E6 and E7 expression is regulated by alternative RNA splicing. This study elucidates for the first time a mapped branch point in HPV16 genome involved in viral oncogene expression.

The transition in spliceosome assembly from complex E to complex A purges surplus U1 snRNPs from alternative splice sites

Spliceosomes are assembled in stages. The first stage forms complex E, which is characterized by the presence of U1 snRNPs base-paired to the 5′ splice site, components recognizing the 3′ splice site and proteins thought to connect them. The splice sites are held in close proximity and the pre-mRNA is committed to splicing. Despite this, the sites for splicing appear not to be fixed until the next complex (A) forms. We have investigated the reasons why 5′ splice sites are not fixed in complex E, using single molecule methods to determine the stoichiometry of U1 snRNPs bound to pre-mRNA with one or two strong 5′ splice sites. In complex E most transcripts with two alternative 5′ splice sites were bound by two U1 snRNPs. However, the surplus U1 snRNPs were displaced during complex A formation in an ATP-dependent process requiring an intact 3′ splice site. This process leaves only one U1 snRNP per complex A, regardless of the number of potential sites. We propose a mechanism for selection of the 5′ splice site. Our results show that constitutive splicing components need not be present in a fixed stoichiometry in a splicing complex.

DNA diagnostics and exon skipping

The predominate form of DNA diagnostics remains nucleic acid sequencing in the research and clinical setting. While DNA sequencing allows a mutation to be correctly identified, only RNA sequencing can confirm the effect of that mutation on the resulting mRNA transcript. In the absence of RNA sequencing, predictions are reliant on either experimental studies or bioinformatic modelling. While each of these approaches provides insights into cellular splicing choices, of which exon skipping is but one, both possess inherent weaknesses. A method which is able to integrate and appropriately weigh the various factors influencing cellular splicing choices into an accurate, comprehensive modelling tool still remains elusive.In this overview chapter, the current methods utilised for DNA diagnostics and the impact of the emerging next-generation sequencing techniques are considered. We explore why RNA remains a problematic medium with which to work. To understand how exon skipping can be predicted from a DNA sequence, the key cis-acting elements influencing splicing are reviewed. Finally, the current methods used to predict exon skipping including RNA-based studies, experimental studies, and bioinformatic modelling approaches are outlined.

Identification and experimental validation of splicing regulatory elements in Drosophila melanogaster reveals functionally conserved splicing enhancers in metazoans

RNA sequence elements involved in the regulation of pre-mRNA splicing have previously been identified in vertebrate genomes by computational methods. Here, we apply such approaches to predict splicing regulatory elements in Drosophila melanogaster and compare them with elements previously found in the human, mouse, and pufferfish genomes. We identified 99 putative exonic splicing enhancers (ESEs) and 231 putative intronic splicing enhancers (ISEs) enriched near weak 5' and 3' splice sites of constitutively spliced introns, distinguishing between those found near short and long introns. We found that a significant proportion (58%) of fly enhancer sequences were previously reported in at least one of the vertebrates. Furthermore, 20% of putative fly ESEs were previously identified as ESEs in human, mouse, and pufferfish; while only two fly ISEs, CTCTCT and TTATAA, were identified as ISEs in all three vertebrate species. Several putative enhancer sequences are similar to characterized binding-site motifs for Drosophila and mammalian splicing regulators. To provide additional evidence for the function of putative ISEs, we separately identified 298 intronic hexamers significantly enriched within sequences phylogenetically conserved among 15 insect species. We found that 73 putative ISEs were among those enriched in conserved regions of the D. melanogaster genome. The functions of nine enhancer sequences were verified in a heterologous splicing reporter, demonstrating that these sequences are sufficient to enhance splicing in vivo. Taken together, these data identify a set of predicted positive-acting splicing regulatory motifs in the Drosophila genome and reveal regulatory sequences that are present in distant metazoan genomes.

The genomic signature of splicing-coupled selection differs between long and short introns

Understanding the function of noncoding regions in the genome, such as introns, is of central importance to evolutionary biology. One approach is to assay for the targets of natural selection. On one hand, the sequence of introns, especially short introns, appears to evolve in an almost neutral manner. Whereas on the other hand, a large proportion of intronic sequence is under selective constraint. This discrepancy is largely dependent on intron length and differences in the methods used to infer selection. We have used a method based on DNA strand asymmetery that does not require comparison with any putatively neutrally evolving sequence, nor sequence conservation between species, to detect selection within introns. The strongest signal we identify is associated with short introns. This signal comes from a family of motifs that could act as cryptic 5′ splice sites during mRNA processing, suggesting a mechanistic justification underlying this signal of selection. Together with an analysis of intron length and splice site strength, we observe that the genomic signature of splicing-coupled selection differs between long and short introns.

Efficient internal exon recognition depends on near equal contributions from the 3' and 5' splice sites

Pre-mRNA splicing is carried out by the spliceosome, which identifies exons and removes intervening introns. In vertebrates, most splice sites are initially recognized by the spliceosome across the exon, because most exons are small and surrounded by large introns. This gene architecture predicts that efficient exon recognition depends largely on the strength of the flanking 3' and 5' splice sites. However, it is unknown if the 3' or the 5' splice site dominates the exon recognition process. Here, we test the 3' and 5' splice site contributions towards efficient exon recognition by systematically replacing the splice sites of an internal exon with sequences of different splice site strengths. We show that the presence of an optimal splice site does not guarantee exon inclusion and that the best predictor for exon recognition is the sum of both splice site scores. Using a genome-wide approach, we demonstrate that the combined 3' and 5' splice site strengths of internal exons provide a much more significant separator between constitutive and alternative exons than either the 3' or the 5' splice site strength alone.

Compensatory signals associated with the activation of human GC 5' splice sites

GC 5′ splice sites (5′ss) are present in ∼1% of human introns, but factors promoting their efficient selection are poorly understood. Here, we describe a case of X-linked agammaglobulinemia resulting from a GC 5′ss activated by a mutation in BTK intron 3. This GC 5′ss was intrinsically weak, yet it was selected in >90% primary transcripts in the presence of a strong and intact natural GT counterpart. We show that efficient selection of this GC 5′ss required a high density of GAA/CAA-containing splicing enhancers in the exonized segment and was promoted by SR proteins 9G8, Tra2β and SC35. The GC 5′ss was efficiently inhibited by splice-switching oligonucleotides targeting either the GC 5′ss itself or the enhancer. Comprehensive analysis of natural GC-AG introns and previously reported pathogenic GC 5′ss showed that their efficient activation was facilitated by higher densities of splicing enhancers and lower densities of silencers than their GT 5′ss equivalents. Removal of the GC-AG introns was promoted to a minor extent by the splice-site strength of adjacent exons and inhibited by flanking Alu repeats, with the first downstream Alus located on average at a longer distance from the GC 5′ss than other transposable elements. These results provide new insights into the splicing code that governs selection of noncanonical splice sites.