U1-like snRNAs lacking complementarity to canonical 5' splice sites

Examining the capacity of human U1 snRNA variants to facilitate pre-mRNA splicing

The human U1 snRNA is encoded by a multigene family consisting of transcribed variants and defective pseudogenes. Many variant U1 (vU1) snRNAs have been demonstrated to not only be transcribed but also processed by the addition of a trimethylated guanosine cap, packaged into snRNPs, and assembled into spliceosomes; however, their capacity to facilitate pre-mRNA splicing has, so far, not been tested. A recent systematic analysis of the human snRNA genes identified 178 U1 snRNA genes that are present in the genome as either tandem arrays or single genes on multiple chromosomes. Of these, 15 were found to be expressed in human tissues and cell lines, although at significantly low levels from their endogenous loci, <0.001% of the canonical U1 snRNA. In this study, we found that placing the variants in the context of the regulatory elements of the RNU1-1 gene improves the expression of many variants to levels comparable to the canonical U1 snRNA. Application of a previously established HeLa cell-based minigene reporter assay to examine the capacity of the vU1 snRNAs to support pre-mRNA splicing revealed that even though the exogenously expressed variant snRNAs were enriched in the nucleus, only a few had a measurable effect on splicing.

Deletions of singular U1 snRNA gene significantly interfere with transcription and 3'-end mRNA formation

Small nuclear RNAs (snRNAs) are structural and functional cores of the spliceosome. In metazoan genomes, each snRNA has multiple copies/variants, up to hundreds in mammals. However, the expressions and functions of each copy/variant in one organism have not been systematically studied. Focus on U1 snRNA genes, we investigated all five copies in Drosophila melanogaster using two series of constructed strains. Analyses of transgenic flies that each have a U1 promoter-driven gfp revealed that U1:21D is the major and ubiquitously expressed copy, and the other four copies have specificities in developmental stages and tissues. Mutant strains that each have a precisely deleted copy of U1-gene exhibited various extents of defects in fly morphology or mobility, especially deletion of U1:82Eb. Interestingly, splicing was changed at limited levels in the deletion strains, while large amounts of differentially-expressed genes and alternative polyadenylation events were identified, showing preferences in the down-regulation of genes with 1-2 introns and selection of proximal sites for 3'-end polyadenylation. In vitro assays suggested that Drosophila U1 variants pulled down fewer SmD2 proteins compared to the canonical U1. This study demonstrates that all five U1-genes in Drosophila have physiological functions in development and play regulatory roles in transcription and 3'-end formation.

PanCancer analysis of somatic mutations in repetitive regions reveals recurrent mutations in snRNA U2

Current somatic mutation callers are biased against repetitive regions, preventing the identification of potential driver alterations in these loci. We developed a mutation caller for repetitive regions, and applied it to study repetitive non protein-coding genes in more than 2200 whole-genome cases. We identified a recurrent mutation at position c.28 in the gene encoding the snRNA U2. This mutation is present in B-cell derived tumors, as well as in prostate and pancreatic cancer, suggesting U2 c.28 constitutes a driver candidate associated with worse prognosis. We showed that the GRCh37 reference genome is incomplete, lacking the U2 cluster in chromosome 17, preventing the identification of mutations in this gene. Furthermore, the 5′-flanking region of WDR74, previously described as frequently mutated in cancer, constitutes a functional copy of U2. These data reinforce the relevance of non-coding mutations in cancer, and highlight current challenges of cancer genomic research in characterizing mutations affecting repetitive genes.

Human spliceosomal snRNA sequence variants generate variant spliceosomes

Human pre-mRNA splicing is primarily catalyzed by the major spliceosome, comprising five small nuclear ribonucleoprotein complexes, U1, U2, U4, U5, and U6 snRNPs, each of which contains the corresponding U-rich snRNA. These snRNAs are encoded by large gene families exhibiting significant sequence variation, but it remains unknown if most human snRNA genes are untranscribed pseudogenes or produce variant snRNAs with the potential to differentially influence splicing. Since gene duplication and variation are powerful mechanisms of evolutionary adaptation, we sought to address this knowledge gap by systematically profiling human U1, U2, U4, and U5 snRNA variant gene transcripts. We identified 55 transcripts that are detectably expressed in human cells, 38 of which incorporate into snRNPs and spliceosomes in 293T cells. All U1 snRNA variants are more than 1000-fold less abundant in spliceosomes than the canonical U1, whereas at least 1% of spliceosomes contain a variant of U2 or U4. In contrast, eight U5 snRNA sequence variants occupy spliceosomes at levels of 1% to 46%. Furthermore, snRNA variants display distinct expression patterns across five human cell lines and adult and fetal tissues. Different RNA degradation rates contribute to the diverse steady state levels of snRNA variants. Our findings suggest that variant spliceosomes containing noncanonical snRNAs may contribute to different tissue- and cell-type-specific alternative splicing patterns.

Croes C, Petersen N, Bosma P, Riccardi F, Pagani F, Pinotti M, van de Graaf SFJ. A Compensatory U1snRNA Partially Rescues FAH Splicing and Protein Expression in a Splicing-Defective Mouse Model of Tyrosinemia Type I

The elucidation of aberrant splicing mechanisms, frequently associated with disease has led to the development of RNA therapeutics based on the U1snRNA, which is involved in 5' splice site (5'ss) recognition. Studies in cellular models have demonstrated that engineered U1snRNAs can rescue different splicing mutation types. However, the assessment of their correction potential in vivo is limited by the scarcity of animal models with the targetable splicing defects. Here, we challenged the U1snRNA in the FAH5961SB mouse model of hepatic fumarylacetoacetate hydrolase (FAH) deficiency (Hereditary Tyrosinemia type I, HT1) due to the FAH c.706G>A splicing mutation. Through minigene expression studies we selected a compensatory U1snRNA (U1F) that was able to rescue this mutation. Intriguingly, adeno-associated virus-mediated delivery of U1F (AAV8-U1F), but not of U1wt, partially rescued FAH splicing in mouse hepatocytes. Consistently, FAH protein was detectable only in the liver of AAV8-U1F treated mice, which displayed a slightly prolonged survival. Moreover, RNA sequencing revealed the negligible impact of the U1F on the splicing profile and overall gene expression, thus pointing toward gene specificity. These data provide early in vivo proof-of-principle of the correction potential of compensatory U1snRNAs in HTI and encourage further optimization on a therapeutic perspective, and translation to other splicing-defective forms of metabolic diseases.

Cutting a Long Intron Short: Recursive Splicing and Its Implications

Over time eukaryotic genomes have evolved to host genes carrying multiple exons separated by increasingly larger intronic, mostly non-protein-coding, sequences. Initially, little attention was paid to these intronic sequences, as they were considered not to contain regulatory information. However, advances in molecular biology, sequencing, and computational tools uncovered that numerous segments within these genomic elements do contribute to the regulation of gene expression. Introns are differentially removed in a cell type-specific manner to produce a range of alternatively-spliced transcripts, and many span tens to hundreds of kilobases. Recent work in human and fruitfly tissues revealed that long introns are extensively processed cotranscriptionally and in a stepwise manner, before their two flanking exons are spliced together. This process, called "recursive splicing," often involves non-canonical splicing elements positioned deep within introns, and different mechanisms for its deployment have been proposed. Still, the very existence and widespread nature of recursive splicing offers a new regulatory layer in the transcript maturation pathway, which may also have implications in human disease.

Splicing of many human genes involves sites embedded within introns

The conventional model for splicing involves excision of each intron in one piece; we demonstrate this inaccurately describes splicing in many human genes. First, after switching on transcription of SAMD4A, a gene with a 134 kb-long first intron, splicing joins the 3′ end of exon 1 to successive points within intron 1 well before the acceptor site at exon 2 is made. Second, genome-wide analysis shows that >60% of active genes yield products generated by such intermediate intron splicing. These products are present at ∼15% the levels of primary transcripts, are encoded by conserved sequences similar to those found at canonical acceptors, and marked by distinctive structural and epigenetic features. Finally, using targeted genome editing, we demonstrate that inhibiting the formation of these splicing intermediates affects efficient exon–exon splicing. These findings greatly expand the functional and regulatory complexity of the human transcriptome.

Functional repeat-derived RNAs often originate from retrotransposon-propagated ncRNAs

The human genome is scattered with repetitive sequences, and the ENCODE project revealed that 60–70% of the genomic DNA is transcribed into RNA. As a consequence, the human transcriptome contains a large portion of repeat-derived RNAs (repRNAs). Here, we present a hypothesis for the evolution of novel functional repeat-derived RNAs from non-coding RNAs (ncRNAs) by retrotransposition. Upon amplification, the ncRNAs can diversify in sequence and subsequently evolve new activities, which can result in novel functions. Non-coding transcripts derived from highly repetitive regions can therefore serve as a reservoir for the evolution of novel functional RNAs. We base our hypothetical model on observations reported for short interspersed nuclear elements derived from 7SL RNA and tRNAs, α satellites derived from snoRNAs and SL RNAs derived from U1 small nuclear RNA. Furthermore, we present novel putative human repeat-derived ncRNAs obtained by the comparison of the Dfam and Rfam databases, as well as several examples in other species. We hypothesize that novel functional ncRNAs can derive also from other repetitive regions and propose Genomic SELEX as a tool for their identification.

U1 small nuclear RNA variants differentially form ribonucleoprotein particles in vitro

The U1 small nuclear (sn)RNA participates in splicing of pre-mRNAs by recognizing and binding to 5' splice sites at exon/intron boundaries. U1 snRNAs associate with 5' splice sites in the form of ribonucleoprotein particles (snRNPs) that are comprised of the U1 snRNA and 10 core components, including U1A, U1-70K, U1C and the 'Smith antigen', or Sm, heptamer. The U1 snRNA is highly conserved across a wide range of taxa; however, a number of reports have identified the presence of expressed U1-like snRNAs in multiple species, including humans. While numerous U1-like molecules have been shown to be expressed, it is unclear whether these variant snRNAs have the capacity to form snRNPs and participate in splicing. The purpose of the present study was to further characterize biochemically the ability of previously identified human U1-like variants to form snRNPs and bind to U1 snRNP proteins. A bioinformatics analysis provided support for the existence of multiple expressed variants. In vitro gel shift assays, competition assays, and immunoprecipitations (IPs) revealed that the variants formed high molecular weight assemblies to varying degrees and associated with core U1 snRNP proteins to a lesser extent than the canonical U1 snRNA. Together, these data suggest that the human U1 snRNA variants analyzed here are unable to efficiently bind U1 snRNP proteins. The current work provides additional biochemical insights into the ability of the variants to assemble into snRNPs.

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

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

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

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

Differentially expressed, variant U1 snRNAs regulate gene expression in human cells

Human U1 small nuclear (sn)RNA, required for splicing of pre-mRNA, is encoded by genes on chromosome 1 (1p36). Imperfect copies of these U1 snRNA genes, also located on chromosome 1 (1q12-21), were thought to be pseudogenes. However, many of these "variant" (v)U1 snRNA genes produce fully processed transcripts. Using antisense oligonucleotides to block the activity of a specific vU1 snRNA in HeLa cells, we have identified global transcriptome changes following interrogation of the Affymetrix Human Exon ST 1.0 array. Our results indicate that this vU1 snRNA regulates expression of a subset of target genes at the level of pre-mRNA processing. This is the first indication that variant U1 snRNAs have a biological function in vivo. Furthermore, some vU1 snRNAs are packaged into unique ribonucleoproteins (RNPs), and many vU1 snRNA genes are differentially expressed in human embryonic stem cells (hESCs) and HeLa cells, suggesting developmental control of RNA processing through expression of different sets of vU1 snRNPs.

Widespread recognition of 5' splice sites by noncanonical base-pairing to U1 snRNA involving bulged nucleotides

An established paradigm in pre-mRNA splicing is the recognition of the 5' splice site (5'ss) by canonical base-pairing to the 5' end of U1 small nuclear RNA (snRNA). We recently reported that a small subset of 5'ss base-pair to U1 in an alternate register that is shifted by 1 nucleotide. Using genetic suppression experiments in human cells, we now demonstrate that many other 5'ss are recognized via noncanonical base-pairing registers involving bulged nucleotides on either the 5'ss or U1 RNA strand, which we term "bulge registers." By combining experimental evidence with transcriptome-wide free-energy calculations of 5'ss/U1 base-pairing, we estimate that 10,248 5'ss (∼5% of human 5'ss) in 6577 genes use bulge registers. Several of these 5'ss occur in genes with mutations causing genetic diseases and are often associated with alternative splicing. These results call for a redefinition of an essential element for gene expression that incorporates these registers, with important implications for the molecular classification of splicing mutations and for alternative splicing.

Correct mRNA processing at a mutant TT splice donor in FANCC ameliorates the clinical phenotype in patients and is enhanced by delivery of suppressor U1 snRNAs

The U1 small nuclear RNA (U1 snRNA) as a component of the major U2-dependent spliceosome recognizes 5' splice sites (5'ss) containing GT as the canonical dinucleotide in the intronic positions +1 and +2. The c.165+1G>T germline mutation in the 5'ss of exon 2 of the Fanconi anemia C (FANCC) gene commonly predicted to prevent correct splicing was identified in nine FA patients from three pedigrees. RT-PCR analysis of the endogenous FANCC mRNA splicing pattern of patient-derived fibroblasts revealed aberrant mRNA processing, but surprisingly also correct splicing at the TT dinucleotide, albeit with lower efficiency. This consequently resulted in low levels of correctly spliced transcript and minute levels of normal posttranslationally processed FANCD2 protein, indicating that this naturally occurring TT splicing might contribute to the milder clinical manifestations of the disease in these patients. Functional analysis of this FANCC 5'ss within splicing reporters revealed that both the noncanonical TT dinucleotide and the genomic context of FANCC were required for the residual correct splicing at this mutant 5'ss. Finally, use of lentiviral vectors as a delivery system to introduce expression cassettes for TT-adapted U1 snRNAs into primary FANCC patient fibroblasts allowed the correction of the DNA-damage-induced G2 cell-cycle arrest in these cells, thus representing an alternative transcript-targeting approach for genetic therapy of inherited splice-site mutations.

Human U1 snRNA forms a new chromatin-associated snRNP with TAF15

The U1 small nuclear RNA (snRNA)--in the form of the U1 spliceosomal Sm small nuclear ribonucleoprotein particle (snRNP) that contains seven Sm and three U1-specific RNP proteins-has a crucial function in the recognition and removal of pre-messenger RNA introns. Here, we show that a fraction of human U1 snRNA specifically associates with the nuclear RNA-binding protein TBP-associated factor 15 (TAF15). We show that none of the known protein components of the spliceosomal U1-Sm snRNP interacts with the newly identified U1-TAF15 snRNP. In addition, the U1-TAF15 snRNP tightly associates with chromatin in an RNA-dependent manner and accumulates in nucleolar caps upon transcriptional inhibition. The Sm-binding motif of U1 snRNA is essential for the biogenesis of both U1-Sm and U1-TAF15 snRNPs, suggesting that the U1-TAF15 particle is produced by remodelling of the U1-Sm snRNP. A demonstration that human U1 snRNA forms at least two structurally distinct snRNPs supports the idea that the U1 snRNA has many nuclear functions.

Evolution of spliceosomal snRNA genes in metazoan animals

While studies of the evolutionary histories of protein families are commonplace, little is known on noncoding RNAs beyond microRNAs and some snoRNAs. Here we investigate in detail the evolutionary history of the nine spliceosomal snRNA families (U1, U2, U4, U5, U6, U11, U12, U4atac, and U6atac) across the completely or partially sequenced genomes of metazoan animals. Representatives of the five major spliceosomal snRNAs were found in all genomes. None of the minor splicesomal snRNAs were detected in nematodes or in the shotgun traces of Oikopleura dioica, while in all other animal genomes at most one of them is missing. Although snRNAs are present in multiple copies in most genomes, distinguishable paralogue groups are not stable over long evolutionary times, although they appear independently in several clades. In general, animal snRNA secondary structures are highly conserved, albeit, in particular, U11 and U12 in insects exhibit dramatic variations. An analysis of genomic context of snRNAs reveals that they behave like mobile elements, exhibiting very little syntenic conservation.

Low U1 snRNP dependence at the NF1 exon 29 donor splice site

Many disease-causing splicing mutations described in the literature produce changes in splice sites (SS) or in exon-regulatory sequences. The delineation of these splice aberrations can provide important insights into novel regulation mechanisms. In this study, we evaluated the effect of patient variations in neurofibromatosis type 1 (NF1) exon 29 and its 5'SS surrounding area on its splicing process. Only two of all nonsense, missense, synonymous and intronic variations analyzed in this study clearly altered exon 29 inclusion/exclusion levels. In particular, the intronic mutation +5g>a had the strongest effect, resulting in total exon exclusion. This finding prompted us to evaluate the exon 29 5'SS in relation to its ability to bind U1 snRNP. This was performed by direct analysis of the ability of U1 to bind to wild-type and mutant donor sites, by engineering an in vitro splicing system to directly evaluate the functional importance of U1 snRNA base pairing with the exon 29 donor site, and by coexpression of mutant U1 snRNP molecules to try to rescue exon 29 inclusion in vivo. The results revealed a low dependency on the presence of U1 snRNP, and suggest that exon 29 donor site definition may depend on alternative mechanisms of 5'SS recognition.

Recognition of atypical 5' splice sites by shifted base-pairing to U1 snRNA

Accurate pre-mRNA splicing is critical for gene expression. The 5' splice site (5' ss) — the highly diverse element at the 5' end of introns — is initially recognized via base-pairing to the 5' end of U1 small nuclear RNA (snRNA). However, many natural 5' ss have a very poor match to the consensus sequence, and are predicted to be very weak. Using genetic suppression experiments in human cells, we demonstrate that some atypical 5' ss are actually efficiently recognized by U1, in an alternative base-pairing register that is shifted by one nucleotide. These atypical 5' ss are phylogenetically widespread, and many of them are conserved. Moreover, shifted base-pairing provides an explanation for the effect of a 5' ss mutation associated with pontocerebellar hypoplasia. The unexpected flexibility in 5' ss/U1 base-pairing challenges an established paradigm, and has broad implications for splice-site prediction algorithms and gene-annotation efforts in genome projects.

De novo search for non-coding RNA genes in the AT-rich genome of Dictyostelium discoideum: performance of Markov-dependent genome feature scoring

Genome data are increasingly important in the computational identification of novel regulatory non-coding RNAs (ncRNAs). However, most ncRNA gene-finders are either specialized to well-characterized ncRNA gene families or require comparisons of closely related genomes. We developed a method for de novo screening for ncRNA genes with a nucleotide composition that stands out against the background genome based on a partial sum process. We compared the performance when assuming independent and first-order Markov-dependent nucleotides, respectively, and used Karlin-Altschul and Karlin-Dembo statistics to evaluate the significance of hits. We hypothesized that a first-order Markov-dependent process might have better power to detect ncRNA genes since nearest-neighbor models have been shown to be successful in predicting RNA structures. A model based on a first-order partial sum process (analyzing overlapping dinucleotides) had better sensitivity and specificity than a zeroth-order model when applied to the AT-rich genome of the amoeba Dictyostelium discoideum. In this genome, we detected 94% of previously known ncRNA genes (at this sensitivity, the false positive rate was estimated to be 25% in a simulated background). The predictions were further refined by clustering candidate genes according to sequence similarity and/or searching for an ncRNA-associated upstream element. We experimentally verified six out of 10 tested ncRNA gene predictions. We conclude that higher-order models, in combination with other information, are useful for identification of novel ncRNA gene families in single-genome analysis of D. discoideum. Our generalizable approach extends the range of genomic data that can be searched for novel ncRNA genes using well-grounded statistical methods.

Requirements for gene silencing mediated by U1 snRNA binding to a target sequence

U1 interference (U1i) is a novel method to block gene expression. U1i requires expression of a 5'-end-mutated U1 snRNA designed to base pair to the 3'-terminal exon of the target gene's pre-mRNA that leads to inhibition of polyadenylation. Here, we show U1i is robust (> or =95%) and a 10-nt target length is sufficient for good silencing. Surprisingly, longer U1 snRNAs, which could increase annealing to the target, fail to improve silencing. Extensive mutagenesis of the 10-bp U1 snRNA:target duplex shows that any single mismatch different from GU at positions 3-8, destroys silencing. However, mismatches within the other positions give partial silencing, suggesting that off-target inhibition could occur. The specificity of U1i may be enhanced, however, by the fact that silencing is impaired by RNA secondary structure or by splicing factors binding nearby, the latter mediated by Arginine-Serine (RS) domains. U1i inhibition can be reconstituted in vivo by tethering of RS domains of U1-70K and U2AF65. These results help to: (i) define good target sites for U1i; (ii) identify and understand natural cellular examples of U1i; (iii) clarify the contribution of hydrogen bonding to U1i and to U1 snRNP binding to 5' splice sites and (iv) understand the mechanism of U1i.

Annotating noncoding RNA genes

Noncoding RNA genes produce a functional RNA product rather than a translated protein. More than 1500 homologs of known "classical" RNA genes can be annotated in the human genome sequence, and automatic homology-based methods predict up to 5000 related sequences. Methods to predict novel RNA genes on a whole-genome scale are immature at present, but their use hints at tens of thousands of such genes in the human genome. Messenger RNA-like transcripts with no protein-coding potential are routinely discovered by high-throughput transcriptome analyses. Meanwhile, various experimental studies have suggested that the vast majority of the human genome is transcribed, although the proportion of the detected RNAs that is functional remains unknown.

Spliceosome assembly and composition

Cells control alternative splicing by modulating assembly of the pre-mRNA splicing machinery at competing splice sites. Therefore, a working knowledge of spliceosome assembly is essential for understanding how alternative splice site choices are achieved. In this chapter, we review spliceosome assembly with particular emphasis on the known steps and factors subject to regulation during alternative splice site selection in mammalian cells. We also review recent advances regarding similarities and differences between the in vivo and in vitro assembly pathways, as well as proofreading mechanisms contributing to the fidelity of splice site selection.