Cryptic splice site usage in exon 7 of the human fibrinogen Bbeta-chain gene is regulated by a naturally silent SF2/ASF binding site within this exon

IL-17D-induced inhibition of DDX5 expression in keratinocytes amplifies IL-36R-mediated skin inflammation

Aberrant RNA splicing in keratinocytes drives inflammatory skin disorders. In the present study, we found that the RNA helicase DDX5 was downregulated in keratinocytes from the inflammatory skin lesions in patients with atopic dermatitis and psoriasis, and that mice with keratinocyte-specific deletion of Ddx5 (Ddx5^∆KC) were more susceptible to cutaneous inflammation. Inhibition of DDX5 expression in keratinocytes was induced by the cytokine interleukin (IL)-17D through activation of the CD93-p38 MAPK-AKT-SMAD2/3 signaling pathway and led to pre-messenger RNA splicing events that favored the production of membrane-bound, intact IL-36 receptor (IL-36R) at the expense of soluble IL-36R (sIL-36R) and to the selective amplification of IL-36R-mediated inflammatory responses and cutaneous inflammation. Restoration of sIL-36R in Ddx5^∆KC mice with experimental atopic dermatitis or psoriasis suppressed skin inflammation and alleviated the disease phenotypes. These findings indicate that IL-17D modulation of DDX5 expression controls inflammation in keratinocytes during inflammatory skin diseases.

Splicing Enhancers at Intron-Exon Borders Participate in Acceptor Splice Sites Recognition

Acceptor splice site recognition (3′ splice site: 3′ss) is a fundamental step in precursor messenger RNA (pre-mRNA) splicing. Generally, the U2 small nuclear ribonucleoprotein (snRNP) auxiliary factor (U2AF) heterodimer recognizes the 3′ss, of which U2AF35 has a dual function: (i) It binds to the intron–exon border of some 3′ss and (ii) mediates enhancer-binding splicing activators’ interactions with the spliceosome. Alternative mechanisms for 3′ss recognition have been suggested, yet they are still not thoroughly understood. Here, we analyzed 3′ss recognition where the intron–exon border is bound by a ubiquitous splicing regulator SRSF1. Using the minigene analysis of two model exons and their mutants, BRCA2 exon 12 and VARS2 exon 17, we showed that the exon inclusion correlated much better with the predicted SRSF1 affinity than 3′ss quality, which were assessed using the Catalog of Inferred Sequence Binding Preferences of RNA binding proteins (CISBP-RNA) database and maximum entropy algorithm (MaxEnt) predictor and the U2AF35 consensus matrix, respectively. RNA affinity purification proved SRSF1 binding to the model 3′ss. On the other hand, knockdown experiments revealed that U2AF35 also plays a role in these exons’ inclusion. Most probably, both factors stochastically bind the 3′ss, supporting exon recognition, more apparently in VARS2 exon 17. Identifying splicing activators as 3′ss recognition factors is crucial for both a basic understanding of splicing regulation and human genetic diagnostics when assessing variants’ effects on splicing.

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.

Ranking noncanonical 5' splice site usage by genome-wide RNA-seq analysis and splicing reporter assays

Most human pathogenic mutations in 5' splice sites affect the canonical GT in positions +1 and +2, leading to noncanonical dinucleotides. On the other hand, noncanonical dinucleotides are observed under physiological conditions in ∼1% of all human 5'ss. It is therefore a challenging task to understand the pathogenic mutation mechanisms underlying the conditions under which noncanonical 5'ss are used. In this work, we systematically examined noncanonical 5' splice site selection, both experimentally using splicing competition reporters and by analyzing a large RNA-seq data set of 54 fibroblast samples from 27 subjects containing a total of 2.4 billion gapped reads covering 269,375 exon junctions. From both approaches, we consistently derived a noncanonical 5'ss usage ranking GC > TT > AT > GA > GG > CT. In our competition splicing reporter assay, noncanonical splicing was strictly dependent on the presence of upstream or downstream splicing regulatory elements (SREs), and changes in SREs could be compensated by variation of U1 snRNA complementarity in the competing 5'ss. In particular, we could confirm splicing at different positions (i.e., -1, +1, +5) of a splice site for all noncanonical dinucleotides "weaker" than GC. In our comprehensive RNA-seq data set analysis, noncanonical 5'ss were preferentially detected in weakly used exon junctions of highly expressed genes. Among high-confidence splice sites, they were 10-fold overrepresented in clusters with a neighboring, more frequently used 5'ss. Conversely, these more frequently used neighbors contained only the dinucleotides GT, GC, and TT, in accordance with the above ranking.

RBM24 promotes U1 snRNP recognition of the mutated 5' splice site in the IKBKAP gene of familial dysautonomia

The 5' splice site mutation (IVS20+6T>C) of the inhibitor of κ light polypeptide gene enhancer in B cells, kinase complex-associated protein (IKBKAP) gene in familial dysautonomia (FD) is at the sixth intronic nucleotide of the 5' splice site. It is known to weaken U1 snRNP recognition and result in an aberrantly spliced mRNA product in neuronal tissue, but normally spliced mRNA in other tissues. Aberrantly spliced IKBKAP mRNA abrogates IKK complex-associated protein (IKAP)/elongator protein 1 (ELP1) expression and results in a defect of neuronal cell development in FD. To elucidate the tissue-dependent regulatory mechanism, we screened an expression library of major RNA-binding proteins (RBPs) with our mammalian dual-color splicing reporter system and identified RBM24 as a regulator. RBM24 functioned as a cryptic intronic splicing enhancer binding to an element (IVS20+13-29) downstream from the intronic 5' splice site mutation in the IKBKAP gene and promoted U1 snRNP recognition only to the mutated 5' splice site (and not the wild-type 5' splice site). Our results show that tissue-specific expression of RBM24 can explain the neuron-specific aberrant splicing of IKBKAP exon 20 in familial dysautonomia, and that ectopic expression of RBM24 in neuronal tissue could be a novel therapeutic target of the disease.

Succession of splicing regulatory elements determines cryptic 5΄ss functionality

A critical step in exon definition is the recognition of a proper splice donor (5΄ss) by the 5’ end of U1 snRNA. In the selection of appropriate 5΄ss, cis-acting splicing regulatory elements (SREs) are indispensable. As a model for 5΄ss recognition, we investigated cryptic 5΄ss selection within the human fibrinogen Bβ-chain gene (FGB) exon 7, where we identified several exonic SREs that simultaneously acted on up- and downstream cryptic 5΄ss. In the FGB exon 7 model system, 5΄ss selection iteratively proceeded along an alternating sequence of U1 snRNA binding sites and interleaved SREs which in principle supported different 3’ exon ends. Like in a relay race, SREs either suppressed a potential 5΄ss and passed the splicing baton on or splicing actually occurred. From RNA-Seq data, we systematically selected 19 genes containing exons with silent U1 snRNA binding sites competing with nearby highly used 5΄ss. Extensive SRE analysis by different algorithms found authentic 5΄ss significantly more supported by SREs than silent U1 snRNA binding sites, indicating that our concept may permit generalization to a model for 5΄ss selection and 3’ exon end definition.

Analysis of Competing HIV-1 Splice Donor Sites Uncovers a Tight Cluster of Splicing Regulatory Elements within Exon 2/2b

The HIV-1 accessory protein Vif is essential for viral replication by counteracting the host restriction factor APOBEC3G (A3G), and balanced levels of both proteins are required for efficient viral replication. Noncoding exons 2/2b contain the Vif start codon between their alternatively used splice donors 2 and 2b (D2 and D2b). For vif mRNA, intron 1 must be removed while intron 2 must be retained. Thus, splice acceptor 1 (A1) must be activated by U1 snRNP binding to either D2 or D2b, while splicing at D2 or D2b must be prevented. Here, we unravel the complex interactions between previously known and novel components of the splicing regulatory network regulating HIV-1 exon 2/2b inclusion in viral mRNAs. In particular, using RNA pulldown experiments and mass spectrometry analysis, we found members of the heterogeneous nuclear ribonucleoparticle (hnRNP) A/B family binding to a novel splicing regulatory element (SRE), the exonic splicing silencer ESS2b, and the splicing regulatory proteins Tra2/SRSF10 binding to the nearby exonic splicing enhancer ESE2b. Using a minigene reporter, we performed bioinformatics HEXplorer-guided mutational analysis to narrow down SRE motifs affecting splice site selection between D2 and D2b. Eventually, the impacts of these SREs on the viral splicing pattern and protein expression were exhaustively analyzed in viral particle production and replication experiments. Masking of these protein binding sites by use of locked nucleic acids (LNAs) impaired Vif expression and viral replication.IMPORTANCE Based on our results, we propose a model in which a dense network of SREs regulates vif mRNA and protein expression, crucial to maintain viral replication within host cells with varying A3G levels and at different stages of infection. This regulation is maintained by several serine/arginine-rich splicing factors (SRSF) and hnRNPs binding to those elements. Targeting this cluster of SREs with LNAs may lead to the development of novel effective therapeutic strategies.

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.

An RNA splicing enhancer that does not act by looping

Out of the loop: Do the proteins bound to an enhancer site on pre-mRNA interact directly with the splice site by diffusion (looping), as is generally accepted, or does the intervening RNA play a role? By inserting a PEG linker between an enhancer sequence and alternative splice sites, the interaction of these two elements can be studied. Intervening RNA was essential for the enhancer activity, which rules out the looping model.

In vitro transcription of compound heterozygous hypofibrinogenemia Matsumoto IX; first identification of FGB IVS6 deletion of 4 nucleotides and FGG IVS3-2A>G causing abnormal RNA splicing

BACKGROUND

We reported a case of hypofibrinogenemia Matsumoto IX (M IX) caused by a novel compound heterozygous mutation involving an FGB IVS6 deletion of 4 nucleotides (Delta4b) (three T, one G; between FGB IVS6-10 and -16) and FGG IVS3-2A/G, which are both identified for the first time. To examine the transcription of mRNA from the M IX gene, we cloned the wild-type and mutant genes into expression vectors.

METHODS

The vectors were transfected into CHO cells and transiently produced wild-type, Bbeta- or gamma-mRNA in the cells. The mRNAs amplified with RT-PCR were analyzed by agarose gel electrophoresis and nucleotide sequencing.

RESULTS

The RT-PCR product from FGB IVS6Delta4b showed aberrant mRNA that included both introns 6 and 7, and that from FGG IVS3-2G showed two aberrant mRNAs, a major one including intron 3 and a minor in which intron 3 was spliced by a cryptic splice site in exon 4. We speculated that the aberrant mRNAs are degraded before translation into proteins, and/or translated variant chains are subjected to quality control and degraded in the cytoplasm.

CONCLUSION

The reduced plasma fibrinogen level of the M IX patient was caused by abnormal RNA splicing of one or both of the FGB and FGG genes.

Polypyrimidine tract binding protein regulates alternative splicing of an aberrant pseudoexon in NF1

In disease-associated genes, understanding the functional significance of deep intronic nucleotide variants represents a difficult challenge. We previously reported that an NF1 intron 30 exonization event is triggered from a single correct nomenclature is 'c.293-279 A>G' mutation [Raponi M, Upadhyaya M & Baralle D (2006) Hum Mutat 27, 294-295]. In this paper, we investigate which characteristics play a role in regulating inclusion of the aberrant pseudoexon. Our investigation shows that pseudoexon inclusion levels are strongly downregulated by polypyrimidine tract binding protein and its homologue neuronal polypyrimidine tract binding protein. In particular, we provide evidence that the functional effect of polypyrimidine tract binding protein is proportional to its concentration, and map the cis-acting elements that are principally responsible for this negative regulation. These results highlight the importance of evaluating local sequence context for diagnostic purposes, and the utility of developing therapies to turn off activated pseudoexons.

Multifactorial interplay controls the splicing profile of Alu-derived exons

Exonization of Alu elements creates primate-specific genomic diversity. Here we combine bioinformatic and experimental methodologies to reconstruct the molecular changes leading to exon selection. Our analyses revealed an intricate network involved in Alu exonization. A typical Alu element contains multiple sites with the potential to serve as 5' splice sites (5'ss). First, we demonstrated the role of 5'ss strength in controlling exonization events. Second, we found that a cryptic 5'ss enhances the selection of a more upstream site and demonstrate that this is mediated by binding of U1 snRNA to the cryptic splice site, challenging the traditional role attributed to U1 snRNA of binding the 5'ss only. Third, we used a simple algorithm to identify specific sequences that determine splice site selection within specific Alu exons. Finally, by inserting identical exons within different sequences, we demonstrated the importance of flanking genomic sequences in determining whether an Alu exon will undergo exonization. Overall, our results demonstrate the complex interplay between at least four interacting layers that affect Alu exonization. These results shed light on the mechanism through which Alu elements enrich the primate transcriptome and allow a better understanding of the exonization process in general.

Cryptic splicing sites are differentially utilized in vivo

It has long been considered that cryptic splice sites are ignored by the splicing machinery in the context of intact genuine splice sites. In the present study, it is shown that cryptic splice sites are utilized in all circumstances, when the authentic site is intact, partially functional or completely abolished. Their use would therefore contribute to a background lack of fidelity in the context of the wild-type sequence. We also found that a mutation at the 5' splice site of beta-globin intron 1 accommodates multiple cryptic splicing pathways, including three previously reported pathways. Focusing on the two major cryptic 5' splice sites within beta-globin exon 1, we show that cryptic splice site selection ex vivo varies depending upon: (a) the cell stage of development during terminal erythroid differentiation; (b) the nature of the mutation at the authentic 5' splice site; and (c) the nature of the promoter. Finally, we found that the two major cryptic 5' splice sites are utilized with differential efficiencies in two siblings sharing the same beta-globin chromosome haplotype in the homozygous state. Collectively, these data suggest that intrinsic, sequence specific factors and cell genetic background factors both contribute to promote a subtle differential use of cryptic splice sites in vivo.

Aberrant 5' splice sites in human disease genes: mutation pattern, nucleotide structure and comparison of computational tools that predict their utilization

Despite a growing number of splicing mutations found in hereditary diseases, utilization of aberrant splice sites and their effects on gene expression remain challenging to predict. We compiled sequences of 346 aberrant 5′splice sites (5′ss) that were activated by mutations in 166 human disease genes. Mutations within the 5′ss consensus accounted for 254 cryptic 5′ss and mutations elsewhere activated 92 de novo 5′ss. Point mutations leading to cryptic 5′ss activation were most common in the first intron nucleotide, followed by the fifth nucleotide. Substitutions at position +5 were exclusively G>A transitions, which was largely attributable to high mutability rates of C/G>T/A. However, the frequency of point mutations at position +5 was significantly higher than that observed in the Human Gene Mutation Database, suggesting that alterations of this position are particularly prone to aberrant splicing, possibly due to a requirement for sequential interactions with U1 and U6 snRNAs. Cryptic 5′ss were best predicted by computational algorithms that accommodate nucleotide dependencies and not by weight-matrix models. Discrimination of intronic 5′ss from their authentic counterparts was less effective than for exonic sites, as the former were intrinsically stronger than the latter. Computational prediction of exonic de novo 5′ss was poor, suggesting that their activation critically depends on exonic splicing enhancers or silencers. The authentic counterparts of aberrant 5′ss were significantly weaker than the average human 5′ss. The development of an online database of aberrant 5′ss will be useful for studying basic mechanisms of splice-site selection, identifying splicing mutations and optimizing splice-site prediction algorithms.

Defective splicing, disease and therapy: searching for master checkpoints in exon definition

The number of aberrant splicing processes causing human disease is growing exponentially and many recent studies have uncovered some aspects of the unexpectedly complex network of interactions involved in these dysfunctions. As a consequence, our knowledge of the various cis- and trans-acting factors playing a role on both normal and aberrant splicing pathways has been enhanced greatly. However, the resulting information explosion has also uncovered the fact that many splicing systems are not easy to model. In fact we are still unable, with certainty, to predict the outcome of a given genomic variation. Nonetheless, in the midst of all this complexity some hard won lessons have been learned and in this survey we will focus on the importance of the wide sequence context when trying to understand why apparently similar mutations can give rise to different effects. The examples discussed in this summary will highlight the fine 'balance of power' that is often present between all the various regulatory elements that define exon boundaries. In the final part, we shall then discuss possible therapeutic targets and strategies to rescue genetic defects of complex splicing systems.