Splice-site pairing is an intrinsically high fidelity process

Preliminary insights into RNA in CSF of pediatric SMA patients after 6 months of nusinersen

BACKGROUND

Spinal muscular atrophy (SMA) is a rare autosomal-recessive neurodegenerative disorder caused by mutations in survival motor neuron 1 (SMN1) gene, and consequent loss of function of SMN protein, which results in progressive loss of lower motor neurons, and muscular wasting. Antisense oligonucleotide (ASO) nusinersen (Spinraza®) modulates the pre-mRNA splicing of the SMN2 gene, allowing rebalance of biologically active SMN. It is administered intrathecally via lumbar puncture after removing an equal amount of cerebrospinal fluid (CSF). Its effect was proven beneficial and approved since 2017 for SMA treatment. Given the direct effect of nusinersen on RNA metabolism, the aim of this project was to evaluate cell-free RNA (cfRNA) in CSF of SMA patients under ASOs treatment for biomarker discovery.

METHODS

By RNA-sequencing approach, RNA obtained from CSF of pediatric SMA type 2 and 3 patients was processed after 6 months of nusinersen treatment, at fifth intrathecal injection (T6), and compared to baseline (T0).

RESULTS

We observed the deregulation of cfRNAs in patients at T6 and we were able to classify these RNAs into disease specific, treatment specific and treatment dependent. Moreover, we subdivided patients into "homogeneous" and "heterogeneous" according to their gene expression pattern. The "heterogeneous" group showed peculiar activation of genes coding for ribosomal components, meaning that in these patients a different molecular effect of nusinersen is observable, even if this specific molecular response was not referable to a clinical pattern.

CONCLUSIONS

This study provides preliminary insights into modulation of gene expression dependent on nusinersen treatment and lays the foundation for biomarkers discovery.

Re-evaluating the impact of alternative RNA splicing on proteomic diversity

Alternative splicing (AS) constitutes a mechanism by which protein-coding genes and long non-coding RNA (lncRNA) genes produce more than a single mature transcript. From plants to humans, AS is a powerful process that increases transcriptome complexity. Importantly, splice variants produced from AS can potentially encode for distinct protein isoforms which can lose or gain specific domains and, hence, differ in their functional properties. Advances in proteomics have shown that the proteome is indeed diverse due to the presence of numerous protein isoforms. For the past decades, with the help of advanced high-throughput technologies, numerous alternatively spliced transcripts have been identified. However, the low detection rate of protein isoforms in proteomic studies raised debatable questions on whether AS contributes to proteomic diversity and on how many AS events are really functional. We propose here to assess and discuss the impact of AS on proteomic complexity in the light of the technological progress, updated genome annotation, and current scientific knowledge.

Deep sequencing of pre-translational mRNPs reveals hidden flux through evolutionarily conserved alternative splicing nonsense-mediated decay pathways

BACKGROUND

Alternative splicing, which generates multiple mRNA isoforms from single genes, is crucial for the regulation of eukaryotic gene expression. The flux through competing splicing pathways cannot be determined by traditional RNA-Seq, however, because different mRNA isoforms can have widely differing decay rates. Indeed, some mRNA isoforms with extremely short half-lives, such as those subject to translation-dependent nonsense-mediated decay (AS-NMD), may be completely overlooked in even the most extensive RNA-Seq analyses.

RESULTS

RNA immunoprecipitation in tandem (RIPiT) of exon junction complex components allows for purification of post-splicing mRNA-protein particles (mRNPs) not yet subject to translation (pre-translational mRNPs) and, therefore, translation-dependent mRNA decay. Here we compare exon junction complex RIPiT-Seq to whole cell RNA-Seq data from HEK293 cells. Consistent with expectation, the flux through known AS-NMD pathways is substantially higher than that captured by RNA-Seq. Our RIPiT-Seq also definitively demonstrates that the splicing machinery itself has no ability to detect reading frame. We identify thousands of previously unannotated splicing events; while many can be attributed to splicing noise, others are evolutionarily conserved events that produce new AS-NMD isoforms likely involved in maintenance of protein homeostasis. Several of these occur in genes whose overexpression has been linked to poor cancer prognosis.

CONCLUSIONS

Deep sequencing of RNAs in post-splicing, pre-translational mRNPs provides a means to identify and quantify splicing events without the confounding influence of differential mRNA decay. For many known AS-NMD targets, the nonsense-mediated decay-linked alternative splicing pathway predominates. Exon junction complex RIPiT-Seq also revealed numerous conserved but previously unannotated AS-NMD events.

Two sides of the same medal: Noncoding mutations reveal new pathological mechanisms and insights into the regulation of gene expression

Noncoding sequences constitute the major part of the human genome and also of pre-mRNAs. Single nucleotide variants in these regions are often overlooked, but may be responsible for much of the variation of phenotypes observed. Mutations in the noncoding part of pre-mRNAs often reveal new and meaningful insights into the regulation of cellular gene expression. Thus, the mechanistic analysis of the pathological mechanism of such mutations will both foster a deeper understanding of the disease and the underlying cellular pathways. Even synonymous mutations can cause diseases, since the primary mRNA sequence not only encodes amino acids, but also encrypts information on RNA-binding proteins and secondary structure. In fact, the RNA sequence directs assembly of a specific mRNP complex, which in turn dictates the fate of the mRNA or regulates its biogenesis. The accumulation of genomic sequence information is increasing at a rapid pace. However, much of the diversity uncovered may not explain the phenotype of a certain syndrome or disease. For this reason, we also emphasize the value of mechanistic studies on pathological mechanisms being complementary to genome-wide studies and bioinformatic approaches. This article is categorized under: RNA Processing > Splicing Regulation/Alternative Splicing RNA Processing > 3' End Processing RNA in Disease and Development > RNA in Disease.

Ultra-deep sequencing reveals pre-mRNA splicing as a sequence driven high-fidelity process

Alternative splicing diversifies mRNA transcripts in human cells. While the spliceosome pairs exons with a high degree of accuracy, the rates of rare aberrant and non-canonical pre-mRNA splicing have not been evaluated at the nucleotide level to determine the quantity and identity of these events across splice junctions. Using ultra-deep sequencing the frequency of aberrant and non-canonical splicing events for three splice junctions flanking exon 7 of SMN1 were determined at single nucleotide resolution. After correction for background noise introduced by PCR amplification and sequencing steps, pre-mRNA splicing was shown to maintain a low overall rate of aberrant and non-canonically spliced events. Several previously unannotated splicing events across 3 exon|intron junctions in SMN1 were identified. Mutations within SMN exon 7 were shown to affect splicing fidelity by modulating RNA secondary structures, by altering the binding site of regulatory proteins and by changing the 5’ splice site strength. Mutations also create a truncated SMN1 exon 7 through the introduction of a de novo non-canonical 5’ splice site. The results from the ultra-deep sequencing approach highlight the impressive fidelity of pre-mRNA splicing and demonstrate that the immediate sequence context around splice sites is the main driving force behind non-canonical splice site pairing.

Deep Splicing Code: Classifying Alternative Splicing Events Using Deep Learning

Alternative splicing (AS) is the process of combining different parts of the pre-mRNA to produce diverse transcripts and eventually different protein products from a single gene. In computational biology field, researchers try to understand AS behavior and regulation using computational models known as “Splicing Codes”. The final goal of these algorithms is to make an in-silico prediction of AS outcome from genomic sequence. Here, we develop a deep learning approach, called Deep Splicing Code (DSC), for categorizing the well-studied classes of AS namely alternatively skipped exons, alternative 5’ss, alternative 3’ss, and constitutively spliced exons based only on the sequence of the exon junctions. The proposed approach significantly improves the prediction and the obtained results reveal that constitutive exons have distinguishable local characteristics from alternatively spliced exons. Using the motif visualization technique, we show that the trained models learned to search for competitive alternative splice sites as well as motifs of important splicing factors with high precision. Thus, the proposed approach greatly expands the opportunities to improve alternative splicing modeling. In addition, a web-server for AS events prediction has been developed based on the proposed method.

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.

Functional characterization and phenotypic spectrum of three recurrent disease-causing deep intronic variants of the CFTR gene

BACKGROUND

The CFTR genotype remains incomplete in 1% of Cystic Fibrosis (CF) cases, because only one or no disease-causing variants is detected after extended analysis. This fraction is probably higher in CFTR-Related Disorders (CFTR-RD). Deep-intronic CFTR variants are putative candidates to fill this gap. However, the recurrence, phenotypic spectrum and full molecular characterization of newly reported variants are unknown.

METHODS

Minigenes and analysis of CFTR transcripts in nasal epithelial cells were used to determine the impact on CFTR splicing of intronic variants that we previously identified by next generation sequencing of the whole CFTR locus. Phenotypic data were collected in 19 patients with CF and CFTR-RD, in whom one of the deep intronic variants has been detected.

RESULTS

Three deep-intronic variants promoted the inclusion of pseudo-exons (PE) in the CFTR transcript, hindering the synthesis of a functional protein. The c.2989-313A > T variant, detected in four patients with CF or CFTR-RD from three different families, led to the inclusion of a 118 bp PE. The c.3469-1304C > G variant promoted the inclusion of a 214 bp-PE and was identified in five patients with CF from four families. Haplotype analysis confirmed that this variant was associated with one CF chromosome of African origin. The most represented variant in our cohort was the c.3874-4522A > G, detected in 10 patients with various phenotypes, from male infertility to CF with pancreatic insufficiency.

CONCLUSION

These three deep intronic CFTR variants are associated with a large phenotypic spectrum, including typical CF. They should be included in CF diagnostic testing and carrier screening strategies.

Poor codon optimality as a signal to degrade transcripts with frameshifts

Frameshifting errors are common and mRNA quality control pathways, such as nonsense-mediated decay (NMD), exist to degrade these aberrant transcripts. Recent work has shown the existence of a genetic link between NMD and codon-usage mediated mRNA decay. Here we present computational evidence that these pathways are synergic for removing frameshifts.

Splicing heterogeneity: separating signal from noise

Single-cell analyses have revealed a tremendous variety among cells in the abundance and chemical composition of RNA. Much of this heterogeneity is due to alternative splicing by the spliceosome. Little is known about how many of the resulting isoforms are biologically functional or just provide noise with little to no impact. The dynamic nature of the spliceosome provides numerous opportunities for regulation but is also the source of stochastic fluctuations. We discuss possible origins of splicing stochasticity, the experimental approaches for studying heterogeneity in isoforms, and the potential biological significance of noisy splicing in development and disease.

Synergistic assembly of human pre-spliceosomes across introns and exons

Most human genes contain multiple introns, necessitating mechanisms to effectively define exons and ensure their proper connection by spliceosomes. Human spliceosome assembly involves both cross-intron and cross-exon interactions, but how these work together is unclear. We examined in human nuclear extracts dynamic interactions of single pre-mRNA molecules with individual fluorescently tagged spliceosomal subcomplexes to investigate how cross-intron and cross-exon processes jointly promote pre-spliceosome assembly. U1 subcomplex bound to the 5' splice site of an intron acts jointly with U1 bound to the 5' splice site of the next intron to dramatically increase the rate and efficiency by which U2 subcomplex is recruited to the branch site/3' splice site of the upstream intron. The flanking 5' splice sites have greater than additive effects implying distinct mechanisms facilitating U2 recruitment. This synergy of 5' splice sites across introns and exons is likely important in promoting correct and efficient splicing of multi-intron pre-mRNAs.

A gene is a segment of DNA that usually carries the information required to build a protein, the molecules responsible for most of life’s processes. This DNA segment is organized in modules, with coding sections separated by portions of non-coding DNA known as introns.

When a gene is ‘turned on’, it gets faithfully copied into a molecule of pre-messenger RNA (pre-mRNA), which contains the alternating coding and non-coding modules. Before it can serve as a template to create a protein, this pre-mRNA must be processed and all the introns removed by a structure called the spliceosome. If this delicate process goes wrong, inaccurate protein templates are produced that may be damaging for the cell.

Spliceosomes are precise molecular ‘scissors’ that can recognize where a coding module stops and an intron starts, and then make a snip in the pre-mRNA to remove the non-coding sequence. The spliceosome is a complex molecular machine formed of numerous parts – including one known as U1 snRNP – that must come together. When a pre-mRNA has several introns, a spliceosome assembles anew for each of them.

Braun et al. designed a new method that allows them to ‘tag’ spliceosomes extracted from a human cell and follow them as they come together. The experiments show that spliceosomes working on different introns in the same pre-mRNA actually help each other out. As one assembles, this helps the spliceosome that processes the neighboring intron to get built. In particular, the U1 snRNPs processing nearby introns collaborate to promote the assembly and activity of the spliceosomes. This teamwork is likely important to guarantee that multiple introns are cut out quickly and accurately.

Transcriptional mutagenesis reduces splicing fidelity in mammalian cells

Splicing fidelity is essential to the maintenance of cellular functions and viability, and mutations or natural variations in pre-mRNA sequences and consequent alteration of splicing have been implicated in the etiology and progression of numerous diseases. The extent to which transcriptional errors or lesion-induced transcriptional mutagenesis (TM) influences splicing fidelity is not currently known. To investigate this, we employed site-specific DNA lesions on the transcribed strand of a minigene splicing reporter in normal mammalian cells. These were the common mutagenic lesions O⁶-methylguanine (O⁶-meG) and 8-oxoguanine (8-oxoG). The minigene splicing reporters were derived from lamin A (LMNA) and proteolipid protein 1 (PLP1), both with known links to human diseases that result from deregulated splicing. In cells with active DNA repair, 1–4% misincorporation occurred opposite the lesions, which increased to 20–40% when repair was compromised. Furthermore, our results reveal that TM at a splice site significantly reduces in vivo splicing fidelity, thereby changing the relative expression of alternative splicing forms in mammalian cells. These findings suggest that splicing defects caused by transcriptional errors can potentially lead to phenotypic cellular changes and increased susceptibility to the development of disease.

The determinants of alternative RNA splicing in human cells

Alternative splicing represents an important level of the regulation of gene function in eukaryotic organisms. It plays a critical role in virtually every biological process within an organism, including regulation of cell division and cell death, differentiation of tissues in the embryo and the adult organism, as well as in cellular response to diverse environmental factors. In turn, studies of the last decade have shown that alternative splicing itself is controlled by different mechanisms. Unfortunately, there is no clear understanding of how these diverse mechanisms, or determinants, regulate and constrain the set of alternative RNA species produced from any particular gene in every cell of the human body. Here, we provide a consolidated overview of alternative splicing determinants including RNA-protein interactions, epigenetic regulation via chromatin remodeling, coupling of transcription-to-alternative splicing, effect of secondary structures in pre-RNA, and function of the RNA quality control systems. We also extensively and critically discuss some mechanistic insights on coordinated inclusion/exclusion of exons during the formation of mature RNA molecules. We conclude that the final structure of RNA is pre-determined by a complex interplay between cis- and trans-acting factors. Altogether, currently available empirical data significantly expand our understanding of the functioning of the alternative splicing machinery of cells in normal and pathological conditions. On the other hand, there are still many blind spots that require further deep investigations.

Linking amyotrophic lateral sclerosis and spinal muscular atrophy through RNA-transcriptome homeostasis: a genomics perspective

In this review, we present our most recent understanding of key biomolecular processes that underlie two motor neuron degenerative disorders, amyotrophic lateral sclerosis, and spinal muscular atrophy. We focus on the role of four multifunctional proteins involved in RNA metabolism (TDP-43, FUS, SMN, and Senataxin) that play a causal role in these diseases. Recent results have led to a novel scenario of intricate connections between these four proteins, bringing transcriptome homeostasis into the spotlight as a common theme in motor neuron degeneration. We review reported functional and physical interactions between these four proteins, highlighting their common association with nuclear bodies and small nuclear ribonucleoprotein particle biogenesis and function. We discuss how these interactions are turning out to be particularly relevant for the control of transcription and chromatin homeostasis, including the recent identification of an association between SMN and Senataxin required to ensure the resolution of DNA-RNA hybrid formation and proper termination by RNA polymerase II. These connections strongly support the existence of common pathways underlying the spinal muscular atrophy and amyotrophic lateral sclerosis phenotype. We also discuss the potential of genome-wide expression profiling, in particular RNA sequencing derived data, to contribute to unravelling the underlying mechanisms. We provide a review of publicly available datasets that have addressed both diseases using these approaches, and highlight the value of investing in cross-disease studies to promote our understanding of the pathways leading to neurodegeneration.

Genetics and biochemistry remain essential in the structural era of the spliceosome

The spliceosome is not a single macromolecular machine. Rather it is a collection of dynamic heterogeneous subcomplexes that rapidly interconvert throughout the course of a typical splicing cycle. Because of this, for many years the only high resolution structures of the spliceosome available were of smaller, isolated protein or RNA components. Consequently much of our current understanding of the spliceosome derives from biochemical and genetic techniques. Now with the publication of multiple, high resolution structures of the spliceosome, some question the relevance of traditional biochemical and genetic techniques to the splicing field. We argue such techniques are not only relevant, but vital for an in depth mechanistic understanding of pre-mRNA splicing.

Selection on Position of Nonsense Codons in Introns

Introns occasionally remain in mature messenger RNAs (mRNAs) due to splicing errors and the translated, aberrant proteins that result represent a metabolic cost and may have other deleterious consequences. The nonsense-mediated decay (NMD) pathway degrades aberrant mRNAs, which it recognizes by the presence of an in-frame premature termination codon (PTC). We investigated whether selection has shaped the location of PTCs in introns to reduce waste and facilitate NMD. We found across seven model organisms, that in both first and last introns, PTCs occur earlier in introns than expected by chance, suggesting that selection favors earlier position. This pattern is more pronounced in species with larger effective population sizes. The pattern does not hold for last introns in the two mammal species, however, perhaps because in these species NMD is not initiated from 3'-terminal introns. We conclude that there is compelling evidence that the location of PTCs is shaped by selection for reduced waste and efficient degradation of aberrant mRNAs.

Structural and mechanistic insights into human splicing factor SF3b complex derived using an integrated approach guided by the cryo-EM density maps

Pre-mRNA splicing in eukaryotes is performed by the spliceosome, a highly complex macromolecular machine. SF3b is a multi-protein complex which recognizes the branch point adenosine of pre-mRNA as part of a larger U2 snRNP or U11/U12 di-snRNP in the dynamic spliceosome machinery. Although a cryo-EM map is available for human SF3b complex, the structure and relative spatial arrangement of all components in the complex are not yet known. We have recognized folds of domains in various proteins in the assembly and generated comparative models. Using an integrative approach involving structural and other experimental data, guided by the available cryo-EM density map, we deciphered a pseudo-atomic model of the closed form of SF3b which is found to be a "fuzzy complex" with highly flexible components and multiplicity of folds. Further, the model provides structural information for 5 proteins (SF3b10, SF3b155, SF3b145, SF3b130 and SF3b14b) and localization information for 4 proteins (SF3b10, SF3b145, SF3b130 and SF3b14b) in the assembly for the first time. Integration of this model with the available U11/U12 di-snRNP cryo-EM map enabled elucidation of an open form. This now provides new insights on the mechanistic features involved in the transition between closed and open forms pivoted by a hinge region in the SF3b155 protein that also harbors cancer causing mutations. Moreover, the open form guided model of the 5' end of U12 snRNA, which includes the branch point duplex, shows that the architecture of SF3b acts as a scaffold for U12 snRNA: pre-mRNA branch point duplex formation with potential implications for branch point adenosine recognition fidelity.

Well-characterized sequence features of eukaryote genomes and implications for ab initio gene prediction

In silico analysis of DNA sequences is an important area of computational biology in the post-genomic era. Over the past two decades, computational approaches for ab initio prediction of gene structure from genome sequence alone have largely facilitated our understanding on a variety of biological questions. Although the computational prediction of protein-coding genes has already been well-established, we are also facing challenges to robustly find the non-coding RNA genes, such as miRNA and lncRNA. Two main aspects of ab initio gene prediction include the computed values for describing sequence features and used algorithm for training the discriminant function, and by which different combinations are employed into various bioinformatic tools. Herein, we briefly review these well-characterized sequence features in eukaryote genomes and applications to ab initio gene prediction. The main purpose of this article is to provide an overview to beginners who aim to develop the related bioinformatic tools.

Selection for reduced translation costs at the intronic 5' end in fungi

It is generally believed that introns are not translated; therefore, the potential intronic features that may be related to the translation step (occurring after splicing) have yet to be thoroughly studied. Here, focusing on four fungi, we performed for the first time a comprehensive study aimed at characterizing how translation efficiency is encoded in introns and affects their evolution. By analysing their intronome we provide evidence of selection for STOP codons close to the intronic 5′ end, and show that the beginning of introns are selected for significantly high translation, presumably to reduce translation and metabolic costs in cases of non-spliced introns. Ribosomal profiling data analysis in Saccharomyces cerevisiae supports the conjecture that in this organism intron retention frequently occurs, introns are partially translated, and their translation efficiency affects organismal fitness. We show that the reported results are more significant in highly translated and highly spliced genes, but are not associated only with genes with a specific function. We also discuss the potential relation of the reported signals to efficient nonsense-mediated decay due to splicing errors. These new discoveries are supported by population-genetics considerations. In addition, they are contributory steps towards a broader understanding of intron evolution and the effect of silent mutations on gene expression and organismal fitness.

mRNA-Associated Processes and Their Influence on Exon-Intron Structure in Drosophila melanogaster

mRNA-associated processes and gene structure in eukaryotes are typically treated as separate research subjects. Here, we bridge this separation and leverage the extensive multidisciplinary work on Drosophila melanogaster to examine the roles that capping, splicing, cleavage/polyadenylation, and telescripting (i.e., the protection of nascent transcripts from premature cleavage/polyadenylation by the splicing factor U1) might play in shaping exon-intron architecture in protein-coding genes. Our findings suggest that the distance between subsequent internal 5′ splice sites (5′ss) in Drosophila genes is constrained such that telescripting effects are maximized, in theory, and thus nascent transcripts are less vulnerable to premature termination. Exceptionally weak 5′ss and constraints on intron-exon size at the gene 5′ end also indicate that capping might enhance the recruitment of U1 and, in turn, promote telescripting at this location. Finally, a positive correlation between last exon length and last 5′ss strength suggests that optimal donor splice sites in the proximity of the pre-mRNA tail may inhibit the processing of downstream polyadenylation signals more than weak donor splice sites do. These findings corroborate and build upon previous experimental and computational studies on Drosophila genes. They support the possibility, hitherto scantly explored, that mRNA-associated processes impose significant constraints on the evolution of eukaryotic gene structure.

Comparative analysis of housekeeping and tissue-selective genes in human based on network topologies and biological properties

Housekeeping genes are genes that are turned on most of the time in almost every tissue to maintain cellular functions. Tissue-selective genes are predominantly expressed in one or a few biologically relevant tissue types. Benefitting from the massive gene expression microarray data obtained over the past decades, the properties of housekeeping and tissue-selective genes can now be investigated on a large-scale manner. In this study, we analyzed the topological properties of housekeeping and tissue-selective genes in the protein-protein interaction (PPI) network. Furthermore, we compared the biological properties and amino acid usage between these two gene groups. The results indicated that there were significant differences in topological properties between housekeeping and tissue-selective genes in the PPI network, and housekeeping genes had higher centrality properties and may play important roles in the complex biological network environment. We also found that there were significant differences in multiple biological properties and many amino acid compositions. The functional genes enrichment and subcellular localizations analysis was also performed to investigate the characterization of housekeeping and tissue-selective genes. The results indicated that the two gene groups showed significant different enrichment in drug targets, disease genes and toxin targets, and located in different subcellular localizations. At last, the discriminations between the properties of two gene groups were measured by the F-score, and expression stage had the most discriminative index in all properties. These findings may elucidate the biological mechanisms for understanding housekeeping and tissue-selective genes and may contribute to better annotate housekeeping and tissue-selective genes in other organisms.

Estimation of the minimum mRNA splicing error rate in vertebrates

The majority of protein coding genes in vertebrates contain several introns that are removed by the mRNA splicing machinery. Errors during splicing can generate aberrant transcripts and degrade the transmission of genetic information thus contributing to genomic instability and disease. However, estimating the error rate of constitutive splicing is complicated by the process of alternative splicing which can generate multiple alternative transcripts per locus and is particularly active in humans. In order to estimate the error frequency of constitutive mRNA splicing and avoid bias by alternative splicing we have characterized the frequency of splice variants at three loci, HPRT, POLB, and TRPV1 in multiple tissues of six vertebrate species. Our analysis revealed that the frequency of splice variants varied widely among loci, tissues, and species. However, the lowest observed frequency is quite constant among loci and approximately 0.1% aberrant transcripts per intron. Arguably this reflects the "irreducible" error rate of splicing, which consists primarily of the combination of replication errors by RNA polymerase II in splice consensus sequences and spliceosome errors in correctly pairing exons.

RNA polymerase errors cause splicing defects and can be regulated by differential expression of RNA polymerase subunits

Errors during transcription may play an important role in determining cellular phenotypes: the RNA polymerase error rate is >4 orders of magnitude higher than that of DNA polymerase and errors are amplified >1000-fold due to translation. However, current methods to measure RNA polymerase fidelity are low-throughout, technically challenging, and organism specific. Here I show that changes in RNA polymerase fidelity can be measured using standard RNA sequencing protocols. I find that RNA polymerase is error-prone, and these errors can result in splicing defects. Furthermore, I find that differential expression of RNA polymerase subunits causes changes in RNA polymerase fidelity, and that coding sequences may have evolved to minimize the effect of these errors. These results suggest that errors caused by RNA polymerase may be a major source of stochastic variability at the level of single cells.

The Silent Sway of Splicing by Synonymous Substitutions

Alternative splicing diversifies mRNA transcripts in human cells. This sequence-driven process can be influenced greatly by mutations, even those that do not change the protein coding potential of the transcript. Synonymous mutations have been shown to alter gene expression through modulation of splicing, mRNA stability, and translation. Using a synonymous position mutation library in SMN1 exon 7, we show that 23% of synonymous mutations across the exon decrease exon inclusion, suggesting that nucleotide identity across the entire exon has been evolutionarily optimized to support a particular exon inclusion level. Although phylogenetic conservation scores are insufficient to identify synonymous positions important for exon inclusion, an alignment of organisms filtered based on similar exon/intron architecture is highly successful. Although many of the splicing neutral mutations are observed to occur, none of the exon inclusion reducing mutants was found in the filtered alignment. Using the modified phylogenetic comparison as an approach to evaluate the impact on pre-mRNA splicing suggests that up to 45% of synonymous SNPs are likely to alter pre-mRNA splicing. These results demonstrate that coding and pre-mRNA splicing pressures co-evolve and that a modified phylogenetic comparison based on the exon/intron architecture is a useful tool in identifying splice altering SNPs.

Widespread alternative and aberrant splicing revealed by lariat sequencing

Alternative splicing is an important and ancient feature of eukaryotic gene structure, the existence of which has likely facilitated eukaryotic proteome expansions. Here, we have used intron lariat sequencing to generate a comprehensive profile of splicing events in Schizosaccharomyces pombe, amongst the simplest organisms that possess mammalian-like splice site degeneracy. We reveal an unprecedented level of alternative splicing, including alternative splice site selection for over half of all annotated introns, hundreds of novel exon-skipping events, and thousands of novel introns. Moreover, the frequency of these events is far higher than previous estimates, with alternative splice sites on average activated at ∼3% the rate of canonical sites. Although a subset of alternative sites are conserved in related species, implying functional potential, the majority are not detectably conserved. Interestingly, the rate of aberrant splicing is inversely related to expression level, with lowly expressed genes more prone to erroneous splicing. Although we validate many events with RNAseq, the proportion of alternative splicing discovered with lariat sequencing is far greater, a difference we attribute to preferential decay of aberrantly spliced transcripts. Together, these data suggest the spliceosome possesses far lower fidelity than previously appreciated, highlighting the potential contributions of alternative splicing in generating novel gene structures.

Splicing predictions reliably classify different types of alternative splicing

Alternative splicing is a key player in the creation of complex mammalian transcriptomes and its misregulation is associated with many human diseases. Multiple mRNA isoforms are generated from most human genes, a process mediated by the interplay of various RNA signature elements and trans-acting factors that guide spliceosomal assembly and intron removal. Here, we introduce a splicing predictor that evaluates hundreds of RNA features simultaneously to successfully differentiate between exons that are constitutively spliced, exons that undergo alternative 5' or 3' splice-site selection, and alternative cassette-type exons. Surprisingly, the splicing predictor did not feature strong discriminatory contributions from binding sites for known splicing regulators. Rather, the ability of an exon to be involved in one or multiple types of alternative splicing is dictated by its immediate sequence context, mainly driven by the identity of the exon's splice sites, the conservation around them, and its exon/intron architecture. Thus, the splicing behavior of human exons can be reliably predicted based on basic RNA sequence elements.

The low information content of Neurospora splicing signals: implications for RNA splicing and intron origin

When we expressed a small (0.9 kb) nonprotein-coding transcript derived from the mitochondrial VS plasmid in the nucleus of Neurospora we found that it was efficiently spliced at one or more of eight 5' splice sites and ten 3' splice sites, which are present apparently by chance in the sequence. Further experimental and bioinformatic analyses of other mitochondrial plasmids, random sequences, and natural nuclear genes in Neurospora and other fungi indicate that fungal spliceosomes recognize a wide range of 5' splice site and branchpoint sequences and predict introns to be present at high frequency in random sequence. In contrast, analysis of intronless fungal nuclear genes indicates that branchpoint, 5' splice site and 3' splice site consensus sequences are underrepresented compared with random sequences. This underrepresentation of splicing signals is sufficient to deplete the nuclear genome of splice sites at locations that do not comprise biologically relevant introns. Thus, the splicing machinery can recognize a wide range of splicing signal sequences, but splicing still occurs with great accuracy, not because the splicing machinery distinguishes correct from incorrect introns, but because incorrect introns are substantially depleted from the genome.

Stem-loop 4 of U1 snRNA is essential for splicing and interacts with the U2 snRNP-specific SF3A1 protein during spliceosome assembly

The pairing of 5′ and 3′ splice sites across an intron is a critical step in spliceosome formation and its regulation. Sharma et al. report a new interaction between stem–loop 4 (SL4) of the U1 snRNA, which recognizes the 5′ splice, and a component of the U2 snRNP complex, which assembles across the intron at the 3′ splice site. U1-SL4 interacts with the SF3A1 protein of the U2 snRNP, and this interaction occurs within prespliceosomal complexes assembled on the pre-mRNA.

The pairing of 5′ and 3′ splice sites across an intron is a critical step in spliceosome formation and its regulation. Interactions that bring the two splice sites together during spliceosome assembly must occur with a high degree of specificity and fidelity to allow expression of functional mRNAs and make particular alternative splicing choices. Here, we report a new interaction between stem–loop 4 (SL4) of the U1 snRNA, which recognizes the 5′ splice site, and a component of the U2 small nuclear ribonucleoprotein particle (snRNP) complex, which assembles across the intron at the 3′ splice site. Using a U1 snRNP complementation assay, we found that SL4 is essential for splicing in vivo. The addition of free U1-SL4 to a splicing reaction in vitro inhibits splicing and blocks complex assembly prior to formation of the prespliceosomal A complex, indicating a requirement for a SL4 contact in spliceosome assembly. To characterize the interactions of this RNA structure, we used a combination of stable isotope labeling by amino acids in cell culture (SILAC), biotin/Neutravidin affinity pull-down, and mass spectrometry. We show that U1-SL4 interacts with the SF3A1 protein of the U2 snRNP. We found that this interaction between the U1 snRNA and SF3A1 occurs within prespliceosomal complexes assembled on the pre-mRNA. Thus, SL4 of the U1 snRNA is important for splicing, and its interaction with SF3A1 mediates contact between the 5′ and 3′ splice site complexes within the assembling spliceosome.

smFRET studies of the 'encounter' complexes and subsequent intermediate states that regulate the selectivity of ligand binding

The selectivity with which a biomolecule can bind its cognate ligand when confronted by the vast array of structurally similar, competing ligands that are present in the cell underlies the fidelity of some of the most fundamental processes in biology. Because they collectively comprise one of only a few methods that can sensitively detect the 'encounter' complexes and subsequent intermediate states that regulate the selectivity of ligand binding, single-molecule fluorescence, and particularly single-molecule fluorescence resonance energy transfer (smFRET), approaches have revolutionized studies of ligand-binding reactions. Here, we describe a widely used smFRET strategy that enables investigations of a large variety of ligand-binding reactions, and discuss two such reactions, aminoacyl-tRNA selection during translation elongation and splice site selection during spliceosome assembly, that highlight both the successes and challenges of smFRET studies of ligand-binding reactions. We conclude by reviewing a number of emerging experimental and computational approaches that are expanding the capabilities of smFRET approaches for studies of ligand-binding reactions and that promise to reveal the mechanisms that control the selectivity of ligand binding with unprecedented resolution.

Human Tra2 proteins jointly control a CHEK1 splicing switch among alternative and constitutive target exons

Alternative splicing--the production of multiple messenger RNA isoforms from a single gene--is regulated in part by RNA binding proteins. While the RBPs transformer2 alpha (Tra2α) and Tra2β have both been implicated in the regulation of alternative splicing, their relative contributions to this process are not well understood. Here we find simultaneous--but not individual--depletion of Tra2α and Tra2β induces substantial shifts in splicing of endogenous Tra2β target exons, and that both constitutive and alternative target exons are under dual Tra2α-Tra2β control. Target exons are enriched in genes associated with chromosome biology including CHEK1, which encodes a key DNA damage response protein. Dual Tra2 protein depletion reduces expression of full-length CHK1 protein, results in the accumulation of the DNA damage marker γH2AX and decreased cell viability. We conclude Tra2 proteins jointly control constitutive and alternative splicing patterns via paralog compensation to control pathways essential to the maintenance of cell viability.

Illuminating the Transcriptome through the Genome

Sequencing the human genome was a huge milestone in genetic research that revealed almost the total DNA sequence required to create a human being. However, in order to function, the DNA genome needs to be expressed as an RNA transcriptome. This article reviews how knowledge of genome sequence information has led to fundamental discoveries in how the transcriptome is processed, with a focus on new system-wide insights into how pre-mRNAs that are encoded by split genes in the genome are rearranged by splicing into functional mRNAs. These advances have been made possible by the development of new post-genome technologies to probe splicing patterns. Transcriptome-wide approaches have characterised a "splicing code" that is embedded within and has a significant role in deciphering the genome, and is deciphered by RNA binding proteins. These analyses have also found that most human genes encode multiple mRNA isoforms, and in some cases proteins, leading in turn to a re-assessment of what exactly a gene is. Analysis of the transcriptome has given insights into how the genome is packaged and transcribed, and is helping to explain important aspects of genome evolution.

Global dissection of alternative splicing in paleopolyploid soybean

Alternative splicing (AS) is common in higher eukaryotes and plays an important role in gene posttranscriptional regulation. It has been suggested that AS varies dramatically among species, tissues, and duplicated gene families of different sizes. However, the genomic forces that govern AS variation remain poorly understood. Here, through genome-wide identification of AS events in the soybean (Glycine max) genome using high-throughput RNA sequencing of 28 samples from different developmental stages, we found that more than 63% of multiexonic genes underwent AS. More AS events occurred in the younger developmental stages than in the older developmental stages for the same type of tissue, and the four main AS types, exon skipping, intron retention, alternative donor sites, and alternative acceptor sites, exhibited different characteristics. Global computational analysis demonstrated that the variations of AS frequency and AS types were significantly correlated with the changes of gene features and gene transcriptional level. Further investigation suggested that the decrease of AS within the genome-wide duplicated genes were due to the diminution of intron length, exon number, and transcriptional level. Altogether, our study revealed that a large number of genes were alternatively spliced in the soybean genome and that variations in gene structure and transcriptional level may play important roles in regulating AS.

Exon identity crisis: disease-causing mutations that disrupt the splicing code

Cis-acting RNA elements control the accurate expression of human multi-exon protein coding genes. Single nucleotide variants altering the fidelity of this regulatory code and, consequently, pre-mRNA splicing are expected to contribute to the etiology of numerous human diseases.

Fungal alternative splicing is associated with multicellular complexity and virulence: a genome-wide multi-species study

Alternative splicing (AS) is a cellular process that increases a cell's coding capacity from a limited set of genes. Although AS is common in higher plants and animals, its prevalence in other eukaryotes is mostly unknown. In fungi the involvement of AS in gene expression and its effect on multi-cellularity and virulence is of great medical and economic interest. We present a genome-wide comparative study of AS in 23 informative fungi of different taxa, based on alignments of public transcript sequences. Random sampling of expressed sequence tags allows for robust and comparable estimations of AS rates. We find that a greater fraction of fungal genes than previously expected is associated with AS. We estimate that on average, 6.4% of the annotated genes are affected by AS, with Cryptococcus neoformans showing an extraordinary rate of 18%. The investigated Basidiomycota show higher average AS rates (8.6%) than the Ascomycota (6.0%), although not significant. We find that multi-cellular complexity and younger evolutionary age associate with higher AS rates. Furthermore, AS affects genes involved in pathogenic lifestyle, particularly in functions of stress response and dimorphic switching. Together, our analysis strongly supports the view that AS is a rather common phenomenon in fungi and associates with higher multi-cellular complexity.

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.

HEXEvent: a database of Human EXon splicing Events

HEXEvent (http://hexevent.mmg.uci.edu) is a new database that permits the user to compile genome-wide exon data sets of human internal exons showing selected splicing events. User queries can be customized based on the type and the frequency of alternative splicing events. For each splicing version of an exon, an ESTs count is given, specifying the frequency of the event. A user-specific definition of constitutive exons can be entered to designate an exon exclusion level still acceptable for an exon to be considered as constitutive. Similarly, the user has the option to define a maximum inclusion level for an exon to be called an alternatively spliced exon. Unlike other existing splicing databases, HEXEvent permits the user to easily extract alternative splicing information for individual, multiple or genome-wide human internal exons. Importantly, the generated data sets are downloadable for further analysis.

Connections between chromatin signatures and splicing

Splicing and alternative splicing are involved in the expression of most human genes, playing key roles in differentiation, cell cycle progression, and development. Misregulation of splicing is frequently associated to disease, which imposes a better understanding of the mechanisms underlying splicing regulation. Accumulated evidence suggests that multiple trans-acting factors and cis-regulatory elements act together to determine tissue-specific splicing patterns. Besides, as splicing is often cotranscriptional, a complex picture emerges in which splicing regulation not only depends on the balance of splicing factor binding to their pre-mRNA target sites but also on transcription-associated features such as protein recruitment to the transcribing machinery and elongation kinetics. Adding more complexity to the splicing regulation network, recent evidence shows that chromatin structure is another layer of regulation that may act through various mechanisms. These span from regulation of RNA polymerase II elongation, which ultimately determines splicing decisions, to splicing factor recruitment by specific histone marks. Chromatin may not only be involved in alternative splicing regulation but in constitutive exon recognition as well. Moreover, splicing was found to be necessary for the proper 'writing' of particular chromatin signatures, giving further mechanistic support to functional interconnections between splicing, transcription and chromatin structure. These links between chromatin configuration and splicing raise the intriguing possibility of the existence of a memory for splicing patterns to be inherited through epigenetic modifications.

Alternative transcription and alternative splicing in cancer

In recent years, the notion of "one gene makes one protein that functions in one signaling pathway" in mammalian cells has been shown to be overly simplistic. Recent genome-wide studies suggest that at least half of the human genes, including many therapeutic target genes, produce multiple protein isoforms through alternative splicing and alternative usage of transcription initiation and/or termination. For example, alternative splicing of the vascular endothelial growth factor gene (VEGFA) produces multiple protein isoforms, which display either pro-angiogenic or anti-angiogenic activities. Similarly, for the majority of human genes, the inclusion or exclusion of exonic sequences enhances the generation of transcript variants and/or protein isoforms that can vary in structure and functional properties. Many of the isoforms produced in this manner are tightly regulated during normal development but are misregulated in cancer cells. Altered expression of transcript variants and protein isoforms for numerous genes is linked with disease and its prognosis, and cancer cells manipulate regulatory mechanisms to express specific isoforms that confer drug resistance and survival advantages. Emerging insights indicate that modulating the expression of transcript and protein isoforms of a gene may hold the key to impeding tumor growth and act as a model for efficient targeting of disease-associated genes at the isoform level. This review highlights the role and regulation of alternative transcription and splicing mechanisms in generating the transcriptome, and the misuse and diagnostic/prognostic potential of alternative transcription and splicing in cancer.

A novel splice variant of folate receptor 4 predominantly expressed in regulatory T cells

BACKGROUND

Regulatory T cells (Tregs) are required for proper maintenance of immunological self-tolerance and immune homeostasis. Folate receptor 4 (FR4) is expressed at high levels in transforming growth factor-beta (TGF-β)-induced Tregs and natural Tregs. Moreover, antibody-mediated targeting of FR4 is sufficient to mediate Treg depletion.

RESULTS

In this study, we describe a novel FR4 transcript variant, FR4D3, in which exon 3 is deleted. The mRNA of FR4D3 encodes a FR4 variant truncated by 189 bp. FR4D3 was found to be predominantly expressed in CD4(+)CD25(+) Treg cells. Overexpression of FR4D3 in CD4(+)CD25(+) Treg cells in vitro stimulated proliferation, which may modulate the ability of these cells to bind and incorporate folic acid.

CONCLUSIONS

Our results suggested that high levels of FR4D3 may be critical to support the substantial proliferative capacity of Treg cells.

Staying on message: ensuring fidelity in pre-mRNA splicing

The faithful expression of genes requires that cellular machinery select substrates with high specificity at each step in gene expression. High specificity is particularly important at the stage of nuclear pre-mRNA splicing, during which the spliceosome selects splice sites and excises intervening introns. With low specificity, the usage of alternative sites would yield insertions, deletions and frame shifts in mRNA. Recently, biochemical, genetic and genome-wide approaches have significantly advanced our understanding of splicing fidelity. In particular, we have learned that DExD/H-box ATPases play a general role in rejecting and discarding suboptimal substrates and that these factors serve as a paradigm for proofreading NTPases in other systems. Recent advances have also defined fundamental questions for future investigations.

Changes in intranuclear mobility of mature snRNPs provide a mechanism for splicing defects in spinal muscular atrophy

It is becoming increasingly clear that defects in RNA metabolism can lead to disease. Spinal muscular atrophy (SMA), a leading genetic cause of infant mortality, results from insufficient amounts of survival motor neuron (SMN) protein. SMN is required for the biogenesis of small nuclear ribonucleoproteins (snRNPs): essential components of the spliceosome. Splicing abnormalities have been detected in models of SMA but it is unclear how lowered SMN affects the fidelity of pre-mRNA splicing. We have examined the dynamics of mature snRNPs in cells depleted of SMN and demonstrated that SMN depletion increases the mobility of mature snRNPs within the nucleus. To dissect the molecular mechanism by which SMN deficiency affects intranuclear snRNP mobility, we employed a panel of inhibitors of different stages of pre-mRNA processing. This in vivo modelling demonstrates that snRNP mobility is altered directly as a result of impaired snRNP maturation. Current models of nuclear dynamics predict that subnuclear structures, including the spliceosome, form by self-organization mediated by stochastic interactions between their molecular components. Thus, alteration of the intranuclear mobility of snRNPs provides a molecular mechanism for splicing defects in SMA.

Error prevention and mitigation as forces in the evolution of genes and genomes

Why are short introns rarely a multiple of three nucleotides long? Why do essential genes cluster? Why are genes in operons often lined up in the order in which they are needed in the encoded pathway? In this Opinion article, we argue that these and many other - ostensibly disparate - observations are all pieces of an emerging picture in which multiple aspects of gene anatomy and genome architecture have evolved in response to error-prone gene expression.

The contribution of mouse models to understanding the pathogenesis of spinal muscular atrophy

Spinal muscular atrophy (SMA), which is caused by inactivating mutations in the survival motor neuron 1 (SMN1) gene, is characterized by loss of lower motor neurons in the spinal cord. The gene encoding SMN is very highly conserved in evolution, allowing the disease to be modeled in a range of species. The similarities in anatomy and physiology to the human neuromuscular system, coupled with the ease of genetic manipulation, make the mouse the most suitable model for exploring the basic pathogenesis of motor neuron loss and for testing potential treatments. Therapies that increase SMN levels, either through direct viral delivery or by enhancing full-length SMN protein expression from the SMN1 paralog, SMN2, are approaching the translational stage of development. It is therefore timely to consider the role of mouse models in addressing aspects of disease pathogenesis that are most relevant to SMA therapy. Here, we review evidence suggesting that the apparent selective vulnerability of motor neurons to SMN deficiency is relative rather than absolute, signifying that therapies will need to be delivered systemically. We also consider evidence from mouse models suggesting that SMN has its predominant action on the neuromuscular system in early postnatal life, during a discrete phase of development. Data from these experiments suggest that the timing of therapy to increase SMN levels might be crucial. The extent to which SMN is required for the maintenance of motor neurons in later life and whether augmenting its levels could treat degenerative motor neuron diseases, such as amyotrophic lateral sclerosis (ALS), requires further exploration.

Cyclic stretch increases splicing noise rate in cultured human fibroblasts

Background

Mechanical forces are known to alter the expression of genes, but it has so far not been reported whether they may influence the fidelity of nucleus-based processes. One experimental approach permitting to address this question is the application of cyclic stretch to cultured human fibroblasts. As a marker for the precision of nucleus-based processes, the number of errors that occur during co-transcriptional splicing can then be measured. This so-called splicing noise is found at low frequency in pre-mRNA splicing.

Findings

The amount of splicing noise was measured by RT-qPCR of seven exon skips from the test genes AATF, MAP3K11, NF1, PCGF2, POLR2A and RABAC1. In cells treated by altered uniaxial cyclic stretching for 18 h, a uniform and significant increase of splicing noise was found for all detectable exon skips.

Conclusion

Our data demonstrate that application of cyclic stretch to cultured fibroblasts correlates with a reduced transcriptional fidelity caused by increasing splicing noise.

Evolution of SR protein and hnRNP splicing regulatory factors

The splicing of pre-mRNAs is an essential step of gene expression in eukaryotes. Introns are removed from split genes through the activities of the spliceosome, a large ribonuclear machine that is conserved throughout the eukaryotic lineage. While unicellular eukaryotes are characterized by less complex splicing, pre-mRNA splicing of multicellular organisms is often associated with extensive alternative splicing that significantly enriches their proteome. The alternative selection of splice sites and exons permits multicellular organisms to modulate gene expression patterns in a cell type-specific fashion, thus contributing to their functional diversification. Alternative splicing is a regulated process that is mainly influenced by the activities of splicing regulators, such as SR proteins or hnRNPs. These modular factors have evolved from a common ancestor through gene duplication events to a diverse group of splicing regulators that mediate exon recognition through their sequence-specific binding to pre-mRNAs. Given the strong correlations between intron expansion, the complexity of pre-mRNA splicing, and the emergence of splicing regulators, it is argued that the increased presence of SR and hnRNP proteins promoted the evolution of alternative splicing through relaxation of the sequence requirements of splice junctions.

Aberrant Single Exon Skipping is not Altered by Age in Exons of NF1, RABAC1, AATF or PCGF2 in Human Blood Cells and Fibroblasts

In human pre-mRNA splicing, infrequent errors occur resulting in erroneous splice products as shown in a genome-wide approach. One characteristic subgroup consists of products lacking one cassette exon. The noise in the splicing process, represented by those misspliced products, can be increased by cold shock treatment or by inhibiting the nonsense mediated decay. Here, we investigated whether the splicing noise frequency increases with age in vivo in peripheral bloods cells or in vitro in cultured and aged fibroblasts from healthy donors. Splicing noise frequency was measured for four erroneously skipped NF1 exons and one exon of RABAC1, AATF and PCGF2 by RT-qPCR. Measurements were validated in cultured fibroblasts treated with cold shock or puromycin. Intragenic but not interpersonal differences were detected in splicing noise frequencies in vivo in peripheral blood cells of 11 healthy donors (15 y–85 y) and in in vitro senescent fibroblasts from three further donors. No correlation to the age of the donors was found in the splicing noise frequencies. Our data demonstrates that splicing error frequencies are not altered by age in peripheral blood cells or in vitro aged fibroblasts in the tested exons of the four investigated genes, indicating a high importance of correct splicing in these proliferating aged cells.

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.