Identification and analysis of U5 snRNA variants in Drosophila

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.

Introns: the "dark matter" of the eukaryotic genome

The emergence of introns was a significant evolutionary leap that is a major distinguishing feature between prokaryotic and eukaryotic genomes. While historically introns were regarded merely as the sequences that are removed to produce spliced transcripts encoding functional products, increasingly data suggests that introns play important roles in the regulation of gene expression. Here, we use an intron-centric lens to review the role of introns in eukaryotic gene expression. First, we focus on intron architecture and how it may influence mechanisms of splicing. Second, we focus on the implications of spliceosomal snRNAs and their variants on intron splicing. Finally, we discuss how the presence of introns and the need to splice them influences transcription regulation. Despite the abundance of introns in the eukaryotic genome and their emerging role regulating gene expression, a lot remains unexplored. Therefore, here we refer to introns as the "dark matter" of the eukaryotic genome and discuss some of the outstanding questions in the field.

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.

Cytogenomic analysis unveils mixed molecular evolution and recurrent chromosomal rearrangements shaping the multigene families on Schistocerca grasshopper genomes

Multigene families are essential components of eukaryotic genomes and play key roles either structurally and functionally. Their modes of evolution remain elusive even in the era of genomics, because multiple multigene family sequences coexist in genomes, particularly in large repetitive genomes. Here, we investigate how the multigene families 18S rDNA, U2 snDNA, and H3 histone evolved in 10 species of Schistocerca grasshoppers with very large and repeat-enriched genomes. Using sequenced genomes and fluorescence in situ hybridization mapping, we find substantial differences between species, including the number of chromosomal clusters, changes in sequence abundance and nucleotide composition, pseudogenization, and association with transposable elements (TEs). The intragenomic analysis of Schistocerca gregaria using long-read sequencing and genome assembly unveils conservation for H3 histone and recurrent pseudogenization for 18S rDNA and U2 snDNA, likely promoted by association with TEs and sequence truncation. Remarkably, TEs were frequently associated with truncated copies, were also among the most abundant in the genome, and revealed signatures of recent activity. Our findings suggest a combined effect of concerted and birth-and-death models driving the evolution of multigene families in Schistocerca over the last 8 million years, and the occurrence of intra- and interchromosomal rearrangements shaping their chromosomal distribution. Despite the conserved karyotype in Schistocerca, our analysis highlights the extensive reorganization of repetitive DNAs in Schistocerca, contributing to the advance of comparative genomics for this important grasshopper genus.

Competition between maturation and degradation drives human snRNA 3' end quality control

Polymerases and exonucleases act on 3' ends of nascent RNAs to promote their maturation or degradation but how the balance between these activities is controlled to dictate the fates of cellular RNAs remains poorly understood. Here, we identify a central role for the human DEDD deadenylase TOE1 in distinguishing the fates of small nuclear (sn)RNAs of the spliceosome from unstable genome-encoded snRNA variants. We found that TOE1 promotes maturation of all regular RNA polymerase II transcribed snRNAs of the major and minor spliceosomes by removing posttranscriptional oligo(A) tails, trimming 3' ends, and preventing nuclear exosome targeting. In contrast, TOE1 promotes little to no maturation of tested U1 variant snRNAs, which are instead targeted by the nuclear exosome. These observations suggest that TOE1 is positioned at the center of a 3' end quality control pathway that selectively promotes maturation and stability of regular snRNAs while leaving snRNA variants unprocessed and exposed to degradation in what could be a widespread mechanism of RNA quality control given the large number of noncoding RNAs processed by DEDD deadenylases.

Regulation of expression of human RNA polymerase II-transcribed snRNA genes

In addition to protein-coding genes, RNA polymerase II (pol II) transcribes numerous genes for non-coding RNAs, including the small-nuclear (sn)RNA genes. snRNAs are an important class of non-coding RNAs, several of which are involved in pre-mRNA splicing. The molecular mechanisms underlying expression of human pol II-transcribed snRNA genes are less well characterized than for protein-coding genes and there are important differences in expression of these two gene types. Here, we review the DNA features and proteins required for efficient transcription of snRNA genes and co-transcriptional 3′ end formation of the transcripts.

Chromatin-associated RNA sequencing (ChAR-seq) maps genome-wide RNA-to-DNA contacts

RNA is a critical component of chromatin in eukaryotes, both as a product of transcription, and as an essential constituent of ribonucleoprotein complexes that regulate both local and global chromatin states. Here, we present a proximity ligation and sequencing method called Chromatin-Associated RNA sequencing (ChAR-seq) that maps all RNA-to-DNA contacts across the genome. Using Drosophila cells, we show that ChAR-seq provides unbiased, de novo identification of targets of chromatin-bound RNAs including nascent transcripts, chromosome-specific dosage compensation ncRNAs, and genome-wide trans-associated RNAs involved in co-transcriptional RNA processing.

Variant snRNPs: New players within the spliceosome system

Much evidence is now accumulating that, in addition to their general role in splicing, the components of the core splicing machinery have extensive regulatory potential. In particular, recent evidence has demonstrated that de-regulation of these factors cause the highest extent of alternative splicing changes compared to de-regulation of the classical splicing regulators. This lack of a general inhibition of splicing resonates the differential splicing effects observed in different disease pathologies associated with specific mutations targeting core spliceosomal components. In this review we will summarize what is currently known regarding the involvement of core spliceosomal U-snRNP complexes in perturbed tissue development and human diseases and argue for the existence of a compensatory mechanism enabling cells to cope with drastic perturbations in core splicing components. This system maintains the correct balance of spliceosomal snRNPs through differential expression of variant (v)U-snRNPs.

Variant U1 snRNAs are implicated in human pluripotent stem cell maintenance and neuromuscular disease

The U1 small nuclear (sn)RNA (U1) is a multifunctional ncRNA, known for its pivotal role in pre-mRNA splicing and regulation of RNA 3' end processing events. We recently demonstrated that a new class of human U1-like snRNAs, the variant (v)U1 snRNAs (vU1s), also participate in pre-mRNA processing events. In this study, we show that several human vU1 genes are specifically upregulated in stem cells and participate in the regulation of cell fate decisions. Significantly, ectopic expression of vU1 genes in human skin fibroblasts leads to increases in levels of key pluripotent stem cell mRNA markers, including NANOG and SOX2. These results reveal an important role for vU1s in the control of key regulatory networks orchestrating the transitions between stem cell maintenance and differentiation. Moreover, vU1 expression varies inversely with U1 expression during differentiation and cell re-programming and this pattern of expression is specifically de-regulated in iPSC-derived motor neurons from Spinal Muscular Atrophy (SMA) type 1 patient's. Accordingly, we suggest that an imbalance in the vU1/U1 ratio, rather than an overall reduction in Uridyl-rich (U)-snRNAs, may contribute to the specific neuromuscular disease phenotype associated with SMA.

Developmental analysis of spliceosomal snRNA isoform expression

Pre-mRNA splicing is a critical step in eukaryotic gene expression that contributes to proteomic, cellular, and developmental complexity. Small nuclear (sn)RNAs are core spliceosomal components; however, the extent to which differential expression of snRNA isoforms regulates splicing is completely unknown. This is partly due to difficulties in the accurate analysis of the spatial and temporal expression patterns of snRNAs. Here, we use high-throughput RNA-sequencing (RNA-seq) data to profile expression of four major snRNAs throughout Drosophila development. This analysis shows that individual isoforms of each snRNA have distinct expression patterns in the embryo, larva, and pharate adult stages. Expression of these isoforms is more heterogeneous during embryogenesis; as development progresses, a single isoform from each snRNA subtype gradually dominates expression. Despite the lack of stable snRNA orthologous groups during evolution, this developmental switching of snRNA isoforms also occurs in distantly related vertebrate species, such as Xenopus, mouse, and human. Our results indicate that expression of snRNA isoforms is regulated and lays the foundation for functional studies of individual snRNA isoforms.

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.

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.

A Drosophila model of spinal muscular atrophy uncouples snRNP biogenesis functions of survival motor neuron from locomotion and viability defects

The spinal muscular atrophy (SMA) protein, survival motor neuron (SMN), functions in the biogenesis of small nuclear ribonucleoproteins (snRNPs). SMN has also been implicated in tissue-specific functions; however, it remains unclear which of these is important for the etiology of SMA. Smn null mutants display larval lethality and show significant locomotion defects as well as reductions in minor-class spliceosomal snRNAs. Despite these reductions, we found no appreciable defects in the splicing of mRNAs containing minor-class introns. Transgenic expression of low levels of either wild-type or an SMA patient-derived form of SMN rescued the larval lethality and locomotor defects; however, snRNA levels were not restored. Thus, the snRNP biogenesis function of SMN is not a major contributor to the phenotype of Smn null mutants. These findings have major implications for SMA etiology because they show that SMN's role in snRNP biogenesis can be uncoupled from the organismal viability and locomotor defects.

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.

Spliceosomal small nuclear RNA genes in 11 insect genomes

The removal of introns from the primary transcripts of protein-coding genes is accomplished by the spliceosome, a large macromolecular complex of which small nuclear RNAs (snRNAs) are crucial components. Following the recent sequencing of the honeybee (Apis mellifera) genome, we used various computational methods, ranging from sequence similarity search to RNA secondary structure prediction, to search for putative snRNA genes (including their promoters) and to examine their pattern of conservation among 11 available insect genomes (A. mellifera, Tribolium castaneum, Bombyx mori, Anopheles gambiae, Aedes aegypti, and six Drosophila species). We identified candidates for all nine spliceosomal snRNA genes in all the analyzed genomes. All the species contain a similar number of snRNA genes, with the exception of A. aegypti, whose genome contains more U1, U2, and U5 genes, and A. mellifera, whose genome contains fewer U2 and U5 genes. We found that snRNA genes are generally more closely related to homologs within the same genus than to those in other genera. Promoter regions for all spliceosomal snRNA genes within each insect species share similar sequence motifs that are likely to correspond to the PSEA (proximal sequence element A), the binding site for snRNA activating protein complex, but these promoter elements vary in sequence among the five insect families surveyed here. In contrast to the other insect species investigated, Dipteran genomes are characterized by a rapid evolution (or loss) of components of the U12 spliceosome and a striking loss of U12-type introns.

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

We have detected a surprising heterogeneity among human spliceosomal U1 small nuclear RNA (snRNA). Most interestingly, we have identified three U1 snRNA variants that lack complementarity to the canonical 5' splice site (5'SS) GU dinucleotide. Furthermore, we have observed heterogeneity among the identified variant U1 snRNA genes caused by single nucleotide polymorphism (SNP). The identified snRNAs were ubiquitously expressed in a variety of human tissues representing different stages of development and displayed features of functional spliceosomal snRNAs, i.e., trimethylated cap structures, association with Sm proteins and presence in nuclear RNA-protein complexes. The unanticipated heterogeneity among spliceosomal snRNAs could contribute to the complexity of vertebrates by expanding the coding capacity of their genomes.