Search

Intron diversity facilitates regulated gene expression and alternative splicing. Spliceosomes excise introns after recognizing their splicing signals: the 5′-splice site (5′ss), branchpoint (BP) and 3′-splice site (3′ss). The latter two signals are recognized by U2 small nuclear ribonucleoprotein (snRNP) and its accessory factors (U2AFs), but longer spacings between them result in weaker splicing. Here, we show that excision of introns with a BP-distant 3′ss (e.g. rap1 intron 2) requires the ubiquitin-fold-activated splicing regulator Sde2 in Schizosaccharomyces pombe. By monitoring splicing-specific ura4 reporters in a collection of S. pombe mutants, Cay1 and Tls1 were identified as additional regulators of this process. The role of Sde2, Cay1 and Tls1 was further confirmed by increasing BP–3′ss spacings in a canonical tho5 intron. We also examined BP-distant exons spliced independently of these factors and observed that RNA secondary structures possibly bridged the gap between the two signals. These proteins may guide the 3′ss towards the spliceosome's catalytic centre by folding the RNA between the BP and 3′ss. Orthologues of Sde2, Cay1 and Tls1, although missing in the intron-poor Saccharomyces cerevisiae, are present in intron-rich eukaryotes, including humans. This type of intron-specific pre-mRNA splicing appears to have evolved for regulated gene expression and alternative splicing of key heterochromatin factors.

Correct pre-mRNA processing in higher eukaryotes vastly depends on splice site recognition. Beyond conserved 5′ss and 3′ss motifs, splicing regulatory elements (SREs) play a pivotal role in this recognition process. Here, we present in silico designed sequences with arbitrary a priori prescribed splicing regulatory HEXplorer properties that can be concatenated to arbitrary length without changing their regulatory properties. We experimentally validated in silico predictions in a massively parallel splicing reporter assay on more than 3000 sequences and exemplarily identified some SRE binding proteins. Aiming at a unified ‘functional splice site strength’ encompassing both U1 snRNA complementarity and impact from neighboring SREs, we developed a novel RNA-seq based 5′ss usage landscape, mapping the competition of pairs of high confidence 5′ss and neighboring exonic GT sites along HBond and HEXplorer score coordinate axes on human fibroblast and endothelium transcriptome datasets. These RNA-seq data served as basis for a logistic 5′ss usage prediction model, which greatly improved discrimination between strong but unused exonic GT sites and annotated highly used 5′ss. Our 5′ss usage landscape offers a unified view on 5′ss and SRE neighborhood impact on splice site recognition, and may contribute to improved mutation assessment in human genetics.

The spliceosome undergoes extensive rearrangements as it assembles onto precursor messenger RNAs. In the earliest assembly step, U1snRNA identifies the 5′ splice site. However, U1snRNA leaves the spliceosome relatively early in assembly, and 5′ splice site identity is subsequently maintained through interactions with U6snRNA, protein factor PRP8, and other components during the rearrangements that build the catalytic site. Using a forward genetic screen in Caenorhabditis elegans, we have identified suppressors of a locomotion defect caused by a 5′ss mutation. Here we report three new suppressor alleles from this screen, two in PRP8 and one in SNRNP200/BRR2. mRNASeq studies of these suppressor strains indicate that they also affect specific native alternative 5′ss, especially for suppressor PRP8 D1549N. A strong suppressor at the unstructured N-terminus of SNRNP200, N18K, indicates a novel role for this region. By examining distinct changes in the splicing of native genes, examining double mutants between suppressors, comparing these new suppressors to previously identified splicing suppressors from yeast, and mapping conserved suppressor residues onto cryoEM structural models of assembling human spliceosomes, we conclude that there are multiple interactions at multiple stages in spliceosome assembly responsible for maintaining the initial 5′ss identified by U1snRNA for entry into the catalytic core.

The specific recognition of splice signals at or near exon-intron junctions is not explained by their weak conservation and instead is postulated to require a multitude of features embedded in the pre-mRNA strand. We explored the possibility of 3D structural scaffold of AdML—a model pre-mRNA substrate—guiding early spliceosomal components to the splice signal sequences. We find that mutations in the non-cognate splice signal sequences impede recruitment of early spliceosomal components due to disruption of the global structure of the pre-mRNA. We further find that the pre-mRNA segments potentially interacting with the early spliceosomal component U1 snRNP are distributed across the intron, that there is a spatial proximity of 5′ and 3′ splice sites within the pre-mRNA scaffold, and that an interplay exists between the structural scaffold and splicing regulatory elements in recruiting early spliceosomal components. These results suggest that early spliceosomal components can recognize a 3D structural scaffold beyond the short splice signal sequences, and that in our model pre-mRNA, this scaffold is formed across the intron involving the major splice signals. This provides a conceptual basis to analyze the contribution of recognizable 3D structural scaffolds to the splicing code across the mammalian transcriptome.

Here we present an update to MutationTaster, our DNA variant effect prediction tool. The new version uses a different prediction model and attains higher accuracy than its predecessor, especially for rare benign variants. In addition, we have integrated many sources of data that only became available after the last release (such as gnomAD and ExAC pLI scores) and changed the splice site prediction model. To more easily assess the relevance of detected known disease mutations to the clinical phenotype of the patient, MutationTaster now provides information on the diseases they cause. Further changes represent a major overhaul of the interfaces to increase user-friendliness whilst many changes under the hood have been designed to accelerate the processing of uploaded VCF files. We also offer an API for the rapid automated query of smaller numbers of variants from within other software. MutationTaster2021 integrates our disease mutation search engine, MutationDistiller, to prioritise variants from VCF files using the patient's clinical phenotype. The novel version is available at https://www.genecascade.org/MutationTaster2021/. This website is free and open to all users and there is no login requirement., Graphical Abstract Graphical AbstractIdentification of disease-causing variants with MutationTaster2021.

Efficient splicing requires a balance between high-fidelity splice-site (SS) selection and speed. In Saccharomyces cerevisiae, Pre-mRNA processing factor 8 (Prp8) helps to balance precise SS selection and rapid, efficient intron excision and exon joining. argonaute1-52 (ago1-52) and incurvata13 (icu13) are hypomorphic alleles of the Arabidopsis thaliana genes ARGONAUTE1 (AGO1) and AUXIN RESISTANT6 (AXR6) that harbor point mutations creating a novel 3′SS and 5′SS, respectively. The spliceosome recognizes these novel SSs, as well as the intact genuine SSs, producing a mixture of wild-type and aberrant mature mRNAs. Here, we characterized five novel mutant alleles of PRP8 (one of the two Arabidopsis co-orthologs of yeast Prp8), naming these alleles morphology of ago1-52 suppressed5 (mas5). In the mas5-1 background, the spliceosome preferentially recognizes the intact genuine 3′SS of ago1-52 and 5′SS of icu13. Since point mutations that damage genuine SSs make the spliceosome prone to recognizing cryptic SSs, we also tested alleles of four genes carrying damaged genuine SSs, finding that mas5-1 did not suppress their missplicing. The mas5-1 and mas5-3 mutations represent a novel class of missplicing suppressors that increase splicing fidelity by hampering the use of novel SSs, but do not alter general pre-mRNA splicing.

Splicing is an essential step of RNA processing for multi-exon genes, in which introns are removed from a precursor RNA, thereby producing mature RNAs containing splice junctions. Here, we develope the RJunBase (www.RJunBase.org), a web-accessible database of three types of RNA splice junctions (linear, back-splice, and fusion junctions) that are derived from RNA-seq data of non-cancerous and cancerous tissues. The RJunBase aims to integrate and characterize all RNA splice junctions of both healthy or pathological human cells and tissues. This new database facilitates the visualization of the gene-level splicing pattern and the junction-level expression profile, as well as the demonstration of unannotated and tumor-specific junctions. The first release of RJunBase contains 682 017 linear junctions, 225 949 back-splice junctions and 34 733 fusion junctions across 18 084 non-cancerous and 11 540 cancerous samples. RJunBase can aid researchers in discovering new splicing-associated targets and provide insights into the identification and assessment of potential neoepitopes for cancer treatment.

An essential heterodimer of the U2AF1 and U2AF2 pre-mRNA splicing factors nucleates spliceosome assembly at polypyrimidine (Py) signals preceding the major class of 3′ splice sites. U2AF1 frequently acquires an S34F-encoding mutation among patients with myelodysplastic syndromes (MDS). The influence of the U2AF1 subunit and its S34F mutation on the U2AF2 conformations remains unknown. Here, we employ single molecule Förster resonance energy transfer (FRET) to determine the influence of wild-type or S34F-substituted U2AF1 on the conformational dynamics of U2AF2 and its splice site RNA complexes. In the absence of RNA, the U2AF1 subunit stabilizes a high FRET value, which by structure-guided mutagenesis corresponds to a closed conformation of the tandem U2AF2 RNA recognition motifs (RRMs). When the U2AF heterodimer is bound to a strong, uridine-rich splice site, U2AF2 switches to a lower FRET value characteristic of an open, side-by-side arrangement of the RRMs. Remarkably, the U2AF heterodimer binds weak, uridine-poor Py tracts as a mixture of closed and open U2AF2 conformations, which are modulated by the S34F mutation. Shifts between open and closed U2AF2 may underlie U2AF1-dependent splicing of degenerate Py tracts and contribute to a subset of S34F-dysregulated splicing events in MDS patients.

Intronic splicing enhancers and silencers (ISEs and ISSs) are two groups of splicing-regulatory elements (SREs) that play critical roles in determining splice-site selection, particularly for alternatively spliced introns or exons. SREs are often short motifs; their mutation or dysregulation of their cognate proteins frequently causes aberrant splicing and results in disease. To date, however, knowledge about SRE sequences and how they regulate splicing remains limited. Here, using an SMN2 minigene, we generated a complete pentamer-sequence library that comprises all possible combinations of 5 nucleotides in intron 7, at a fixed site downstream of the 5′ splice site. We systematically analyzed the effects of all 1023 mutant pentamers on exon 7 splicing, in comparison to the wild-type minigene, in HEK293 cells. Our data show that the majority of pentamers significantly affect exon 7 splicing: 584 of them are stimulatory and 230 are inhibitory. To identify actual SREs, we utilized a motif set enrichment analysis (MSEA), from which we identified groups of stimulatory and inhibitory SRE motifs. We experimentally validated several strong SREs in SMN1/2 and other minigene settings. Our results provide a valuable resource for understanding how short RNA sequences regulate splicing. Many novel SREs can be explored further to elucidate their mechanism of action.

RBFOX2 controls the splicing of a large number of transcripts implicated in cell differentiation and development. Parsing RNA-binding protein datasets, we uncover that RBFOX2 can interact with hnRNPC, hnRNPM and SRSF1 to regulate splicing of a broad range of splicing events using different sequence motifs and binding modes. Using immunoprecipitation, specific RBP knockdown, RNA-seq and splice-sensitive PCR, we show that RBFOX2 can target splice sites using three binding configurations: single, multiple or secondary modes. In the single binding mode RBFOX2 is recruited to its target splice sites through a single canonical binding motif, while in the multiple binding mode RBFOX2 binding sites include the adjacent binding of at least one other RNA binding protein partner. Finally, in the secondary binding mode RBFOX2 likely does not bind the RNA directly but is recruited to splice sites lacking its canonical binding motif through the binding of one of its protein partners. These dynamic modes bind distinct sets of transcripts at different positions and distances relative to alternative splice sites explaining the heterogeneity of RBFOX2 targets and splicing outcomes.

Group II introns can self-splice from RNA transcripts through branching, hydrolysis and circularization, being released as lariats, linear introns and circles, respectively. In contrast to branching, the circularization pathway is mostly based on assumptions and has been largely overlooked. Here, we address the molecular details of both transesterification reactions of the group II intron circularization pathway in vivo. We show that free E1 is recruited by the intron through base pairing interactions and that released intron circles can generate free E1 by the spliced exon reopening reaction. The first transesterification reaction was found to be induced inaccurately by the 3′OH of the terminal residue of free E1 at the 3′ splice site, producing circularization intermediates with heterogeneous 3′ ends. Nevertheless, specific terminal 3′OH, selected by a molecular ruler, was shown to precisely attack the 5′ splice site and release intron circles with 3′–5′ rather than 2′–5′ bonds at their circularization junction. Our work supports a circularization model where the recruitment of free E1 and/or displacement of cis-E1 induce a conformational change of the intron active site from the pre-5′ to the pre-3′ splice site processing conformation, suggesting how circularization might initiate at the 3′ instead of the 5′ splice site.

The essential splicing factor Cwc24 contains a zinc-finger (ZF) domain required for its function in splicing. Cwc24 binds over the 5′ splice site after the spliceosome is activated, and its binding prior to Prp2-mediated spliceosome remodeling is important for proper interactions of U5 and U6 with the 5′ splice site sequence and selection of the 5′ splice site. Here, we show that Cwc24 transiently interacts with the 5′ splice site in formation of the functional RNA catalytic core during spliceosome remodeling, and the ZF-motif is required for specific interaction of Cwc24 with the 5′ splice site. Deletion of the ZF domain or mutation of the conserved ZF residues greatly weakened the association of Cwc24 with the spliceosome, and lowered the affinity and specificity of its interaction with the 5′ splice site, resulting in atypical interactions of U5, U6 and Prp8 with the 5′ splice site, and aberrant cleavage at the 5′ splice site. Our results reveal a crucial role of the Cwc24 ZF-motif for defining 5′ splice site selection in the first splicing step.

Chromatin organization and epigenetic markers influence splicing, though the magnitudes of these effects and the mechanisms are largely unknown. Here, we demonstrate that linker histone H1.5 influences mRNA splicing. We observed that linker histone H1.5 binds DNA over splice sites of short exons in human lung fibroblasts (IMR90 cells). We found that association of H1.5 with these splice sites correlated with the level of inclusion of alternatively spliced exons. Exons marked by H1.5 had more RNA polymerase II (RNAP II) stalling near the 3′ splice site than did exons not associated with H1.5. In cells depleted of H1.5, we showed that the inclusion of five exons evaluated decreased and that RNAP II levels over these exons were also reduced. Our findings indicate that H1.5 is involved in regulation of splice site selection and alternative splicing, a function not previously demonstrated for linker histones.

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.

Precise identification of correct exon–intron boundaries is a prerequisite to analyze the location and structure of genes. The existing framework for genomic signals, delineating exon and introns in a genomic segment, seems insufficient, predominantly due to poor sequence consensus as well as limitations of training on available experimental data sets. We present here a novel concept for characterizing exon–intron boundaries in genomic segments on the basis of structural and energetic properties. We analyzed boundary junctions on both sides of all the exons (3 28 368) of protein coding genes from human genome (GENCODE database) using 28 structural and three energy parameters. Study of sequence conservation at these sites shows very poor consensus. It is observed that DNA adopts a unique structural and energy state at the boundary junctions. Also, signals are somewhat different for housekeeping and tissue specific genes. Clustering of 31 parameters into four derived vectors gives some additional insights into the physical mechanisms involved in this biological process. Sites of structural and energy signals correlate well to the positions playing important roles in pre-mRNA splicing.

The rMAPS2 (RNA Map Analysis and Plotting Server 2) web server, freely available at http://rmaps.cecsresearch.org/, has provided the high-throughput sequencing data research community with curated tools for the identification of RNA binding protein sites. rMAPS2 analyzes differential alternative splicing or CLIP peak data obtained from high-throughput sequencing data analysis tools like MISO, rMATS, Piranha, PIPE-CLIP and PARalyzer, and then, graphically displays enriched RNA-binding protein target sites. The initial release of rMAPS focused only on the most common alternative splicing event, skipped exon or exon skipping. However, there was a high demand for the analysis of other major types of alternative splicing events, especially for retained intron events since this is the most common type of alternative splicing in plants, such as Arabidopsis thaliana. Here, we expanded the implementation of rMAPS2 to facilitate analyses for all five major types of alternative splicing events: skipped exon, mutually exclusive exons, alternative 5′ splice site, alternative 3′ splice site and retained intron. In addition, by employing multi-threading, rMAPS2 has vastly improved the user experience with significant reductions in running time, ∼3.5 min for the analysis of all five major alternative splicing types at once.

Cancer-associated mutations of the core splicing factor 3 B1 (SF3B1) result in selection of novel 3′ splice sites (3′SS), but precise molecular mechanisms of oncogenesis remain unclear. SF3B1 stabilizes the interaction between U2 snRNP and branch point (BP) on the pre-mRNA. It has hence been speculated that a change in BP selection is the basis for novel 3′SS selection. Direct quantitative determination of BP utilization is however technically challenging. To define BP utilization by SF3B1-mutant spliceosomes, we used an overexpression approach in human cells as well as a complementary strategy using isogenic murine embryonic stem cells with monoallelic K700E mutations constructed via CRISPR/Cas9-based genome editing and a dual vector homology-directed repair methodology. A synthetic minigene library with degenerate regions in 3′ intronic regions (3.4 million individual minigenes) was used to compare BP usage of SF3B1K700E and SF3B1WT. Using this model, we show that SF3B1K700E spliceosomes utilize non-canonical sequence variants (at position −1 relative to BP adenosine) more frequently than wild-type spliceosomes. These predictions were confirmed using minigene splicing assays. Our results suggest a model of BP utilization by mutant SF3B1 wherein it is able to utilize non-consensus alternative BP sequences by stabilizing weaker U2-BP interactions.

The spliceosome is assembled via sequential interactions of pre-mRNA with five small nuclear RNAs and many proteins. Recent determination of cryo-EM structures for several spliceosomal complexes has provided deep insights into interactions between spliceosomal components and structural changes of the spliceosome between steps, but information on how the proteins interact with pre-mRNA to mediate the reaction is scarce. By systematic analysis of proteins interacting with the splice sites (SSs), we have identified many previously unknown interactions of spliceosomal components with the pre-mRNA. Prp8 directly binds over the 5′SS and the branch site (BS) for the first catalytic step, and the 5′SS and 3′SS for the second step. Switching the Prp8 interaction from the BS to the 3′SS requires Slu7, which interacts dynamically with pre-mRNA first, and then interacts stably with the 3′-exon after Prp16-mediated spliceosome remodeling. Our results suggest that Prp8 plays a key role in positioning the 5′SS and 3′SS, facilitated by Slu7 through interactions with Prp8 and substrate RNA to advance exon ligation. We also provide evidence that Prp16 first docks on the intron 3′ tail, then translocates in the 3′ to 5′ direction on remodeling the spliceosome.

RNA splicing is a highly regulated process dependent on sequences near splice sites. Insertions of Alu retrotransposons can disrupt splice sites or bind splicing regulators. We hypothesized that some common inherited polymorphic Alu insertions are responsible for splicing QTLs (sQTL). We focused on intronic Alu variants mapping within 100 bp of an alternatively used exon and screened for those that alter splicing. We identify five loci, 21.7% of those assayed, where the polymorphic Alu alters splicing. While in most cases the Alu promotes exon skipping, at one locus the Alu increases exon inclusion. Of particular interest is an Alu polymorphism in the CD58 gene. Reduced CD58 expression is associated with risk for developing multiple sclerosis. We show that the Alu insertion promotes skipping of CD58 exon 3 and results in a frameshifted transcript, indicating that the Alu may be the causative variant for increased MS risk at this locus. Using RT-PCR analysis at the endogenous locus, we confirm that the Alu variant is a sQTL for CD58. In summary, altered splicing efficiency is a common functional consequence of Alu polymorphisms including at least one instance where the variant is implicated in disease risk. This work broadens our understanding of splicing regulatory sequences around exons.

Splice-switching antisense oligonucleotides (ASOs), which bind specific RNA-target sequences and modulate pre-mRNA splicing by sterically blocking the binding of splicing factors to the pre-mRNA, are a promising therapeutic modality to treat a range of genetic diseases. ASOs are typically 15–25 nt long and considered to be highly specific towards their intended target sequence, typically elements that control exon definition and/or splice-site recognition. However, whether or not splice-modulating ASOs also induce hybridization-dependent mis-splicing of unintended targets has not been systematically studied. Here, we tested the in vitro effects of splice-modulating ASOs on 108 potential off-targets predicted on the basis of sequence complementarity, and identified 17 mis-splicing events for one of the ASOs tested. Based on analysis of data from two overlapping ASO sequences, we conclude that off-target effects are difficult to predict, and the choice of ASO chemistry influences the extent of off-target activity. The off-target events caused by the uniformly modified ASOs tested in this study were significantly reduced with mixed-chemistry ASOs of the same sequence. Furthermore, using shorter ASOs, combining two ASOs, and delivering ASOs by free uptake also reduced off-target activity. Finally, ASOs with strategically placed mismatches can be used to reduce unwanted off-target splicing events.

Alternative splicing (AS) is a regulatory mechanism of gene expression that greatly expands the coding capacities of genomes by allowing the generation of multiple mRNAs from a single gene. In Drosophila, the mod(mdg4) locus is an extreme example of AS that produces more than 30 different mRNAs via trans-splicing that joins together the common exons and the 3′ variable exons generated from alternative promoters. To map the regions required for trans-splicing, we have developed an assay for measuring trans-splicing events and identified a 73-bp region in the last common intron that is critical for trans-splicing of three pre-mRNAs synthesized from different DNA strands. We have also found that conserved sequences in the distal part of the last common intron induce polyadenylation-independent transcription termination and are enriched by paused RNA polymerase II (RNAP II). These results suggest that all mod(mdg4) mRNAs are formed by joining in trans the 5′ splice site in the last common exon with the 3′ splice site in one of the alternative exons.

Variant interpretation is the key issue in molecular diagnosis. Spliceogenic variants exemplify this issue as each nucleotide variant can be deleterious via disruption or creation of splice site consensus sequences. Consequently, reliable in silico prediction of variant spliceogenicity would be a major improvement. Thanks to an international effort, a set of 395 variants studied at the mRNA level and occurring in 5′ and 3′ consensus regions (defined as the 11 and 14 bases surrounding the exon/intron junction, respectively) was collected for 11 different genes, including BRCA1, BRCA2, CFTR and RHD, and used to train and validate a new prediction protocol named Splicing Prediction in Consensus Elements (SPiCE). SPiCE combines in silico predictions from SpliceSiteFinder-like and MaxEntScan and uses logistic regression to define optimal decision thresholds. It revealed an unprecedented sensitivity and specificity of 99.5 and 95.2%, respectively, and the impact on splicing was correctly predicted for 98.8% of variants. We therefore propose SPiCE as the new tool for predicting variant spliceogenicity. It could be easily implemented in any diagnostic laboratory as a routine decision making tool to help geneticists to face the deluge of variants in the next-generation sequencing era. SPiCE is accessible at (https://sourceforge.net/projects/spicev2-1/).

Familial dysautonomia (FD) is a rare inherited neurodegenerative disorder caused by a point mutation in the IKBKAP gene that results in defective splicing of its pre-mRNA. The mutation weakens the 5′ splice site of exon 20, causing this exon to be skipped, thereby introducing a premature termination codon. Though detailed FD pathogenesis mechanisms are not yet clear, correcting the splicing defect in the relevant tissue(s), thus restoring normal expression levels of the full-length IKAP protein, could be therapeutic. Splice-switching antisense oligonucleotides (ASOs) can be effective targeted therapeutics for neurodegenerative diseases, such as nusinersen (Spinraza), an approved drug for spinal muscular atrophy. Using a two-step screen with ASOs targeting IKBKAP exon 20 or the adjoining intronic regions, we identified a lead ASO that fully restored exon 20 splicing in FD patient fibroblasts. We also characterized the corresponding cis-acting regulatory sequences that control exon 20 splicing. When administered into a transgenic FD mouse model, the lead ASO promoted expression of full-length human IKBKAP mRNA and IKAP protein levels in several tissues tested, including the central nervous system. These findings provide insights into the mechanisms of IKBKAP exon 20 recognition, and pre-clinical proof of concept for an ASO-based targeted therapy for FD.

Spinal muscular atrophy (SMA) is caused by deletions or mutations of the Survival Motor Neuron 1 (SMN1) gene coupled with predominant skipping of SMN2 exon 7. The only approved SMA treatment is an antisense oligonucleotide that targets the intronic splicing silencer N1 (ISS-N1), located downstream of the 5′ splice site (5′ss) of exon 7. Here, we describe a novel approach to exon 7 splicing modulation through activation of a cryptic 5′ss (Cr1). We discovered the activation of Cr1 in transcripts derived from SMN1 that carries a pathogenic G-to-C mutation at the first position (G1C) of intron 7. We show that Cr1-activating engineered U1 snRNAs (eU1s) have the unique ability to reprogram pre-mRNA splicing and restore exon 7 inclusion in SMN1 carrying a broad spectrum of pathogenic mutations at both the 3′ss and 5′ss of the exon 7. Employing a splicing-coupled translation reporter, we demonstrate that mRNAs generated by an eU1-induced activation of Cr1 produce full-length SMN. Our findings underscore a wider role for U1 snRNP in splicing regulation and reveal a novel approach for the restoration of SMN exon 7 inclusion for a potential therapy of SMA.

Inhibition of the Akt kinase activates HPV16 late gene expression by reducing HPV16 early polyadenylation and by activating HPV16 late L1 mRNA splicing. We identified ‘hot spots’ for RNA binding proteins at the early polyA signal and at splice sites on HPV16 late mRNAs. We observed that hnRNP L was associated with sequences at all HPV16 late splice sites and at the early polyA signal. Akt kinase inhibition resulted in hnRNP L dephosphorylation and reduced association of hnRNP L with HPV16 mRNAs. This was accompanied by an increased binding of U2AF65 and Sam68 to HPV16 mRNAs. Furthermore, siRNA knock-down of hnRNP L or Akt induced HPV16 gene expression. Treatment of HPV16 immortalized keratinocytes with Akt kinase inhibitor reduced hnRNP L binding to HPV16 mRNAs and induced HPV16 L1 mRNA production. Finally, deletion of the hnRNP L binding sites in HPV16 subgenomic expression plasmids resulted in activation of HPV16 late gene expression. In conclusion, the Akt kinase inhibits HPV16 late gene expression at the level of RNA processing by controlling the RNA-binding protein hnRNP L. We speculate that Akt kinase activity upholds an intracellular milieu that favours HPV16 early gene expression and suppresses HPV16 late gene expression.

Although the involvement of Ser/Arg-rich (SR) proteins in RNA metabolism is well documented, their role in vertebrate development remains elusive. We, therefore, elected to take advantage of the zebrafish model organism to study the SR genes' functions using the splicing morpholino (sMO) microinjection and the programmable site-specific nucleases. Consistent with previous research, we revealed discrepancies between the mutant and morphant phenotypes and we show that these inconsistencies may result from a large number of unsuspected inadvertent morpholino RNA targets. While microinjection of MOs directed against srsf5a (sMOsrsf5a) led to developmental defects, the corresponding homozygous mutants did not display any phenotypic traits. Furthermore, microinjection of sMOsrsf5a into srsf5a−/− led to the previously observed morphant phenotype. Similar findings were observed for other SR genes. sMOsrsf5a alternative target genes were identified using deep mRNA sequencing. We uncovered that only 11 consecutive bases complementary to sMOsrsf5a are sufficient for binding and subsequent blocking of splice sites. In addition, we observed that sMOsrsf5a secondary targets can be reduced by increasing embryos growth temperature after microinjection. Our data contribute to the debate about MO specificity, efficacy and the number of unknown targeted sequences.

U2 auxiliary factor 1 (U2AF1) functions in 3′-splice site selection during pre-mRNA processing. Alternative usage of duplicated tandem exons in U2AF1 produces two isoforms, U2AF1a and U2AF1b, but their functional differences are unappreciated due to their homology. Through integrative approaches of genome editing, customized-transcriptome profiling and crosslinking-mediated interactome analyses, we discovered that the expression of U2AF1 isoforms is controlled by mTOR and they exhibit a distinctive molecular profile for the splice site and protein interactomes. Mechanistic dissection of mutually exclusive alternative splicing events revealed that U2AF1 isoforms’ inherent differential preferences of nucleotide sequences and their stoichiometry determine the 3′-splice site. Importantly, U2AF1a-driven transcriptomes feature alternative splicing events in the 5′-untranslated region (5′-UTR) that are favorable for translation. These findings unveil distinct roles of duplicated tandem exon-derived U2AF1 isoforms in the regulation of the transcriptome and suggest U2AF1a-driven 5′-UTR alternative splicing as a molecular mechanism of mTOR-regulated translational control.