Developmental analysis of spliceosomal snRNA isoform expression

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.

sRNAfrag: A pipeline and suite of tools to analyze fragmentation in small RNA sequencing data

Fragments derived from small RNAs such as small nucleolar RNAs hold biological relevance. However, they remain poorly understood, calling for more comprehensive methods for analysis. We developed sRNAfrag, a standardized workflow and set of scripts to quantify and analyze sRNA fragmentation of any biotype. In a benchmark, it is able to detect loci of mature microRNAs fragmented from precursors and, utilizing multi-mapping events, the conserved 5' seed sequence of miRNAs which we believe may extraoplate to other small RNA fragments. The tool detected 1411 snoRNA fragment conservation events between 2/4 eukaryotic species, providing the opportunity to explore motifs and fragmentation patterns not only within species, but between. Availability: https://github.com/kenminsoo/sRNAfrag.

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.

Chemical crosslinking and ligation methods for in vivo analysis of RNA structures and interactions

RNA structures and interactions in living cells drive a variety of biological processes and play critical roles in physiology and disease states. However, studies of RNA structures and interactions have been challenging due to limitations in available technologies. Direct determination of structures in vitro has been only possible to a small number of RNAs with limited sizes and conformations. We recently introduced two chemical crosslink-ligation techniques that enabled studies of transcriptome-wide secondary and tertiary structures and their dynamics. In a dramatically improved version of the psoralen analysis of RNA interactions and structures (PARIS2) method, we detailed the synthesis and use of amotosalen, a highly soluble psoralen analogue, and enhanced enzymology for higher efficiency duplex capture. We also introduced spatial 2'-hydroxyl acylation reversible crosslinking (SHARC) with exonuclease (exo) trimming, a method which utilizes a novel crosslinker class that targets the 2'-OH to capture three-dimensional (3D) structures. Both are powerful orthogonal approaches for solving in vivo RNA structure and interactions, integrating crosslinking, exo trimming, proximity ligation, and high throughput sequencing. In this chapter, we present a detailed protocol for the methods and highlight steps that outperform existing crosslink-ligation approaches.

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.

SnoRNA copy regulation affects family size, genomic location and family abundance levels

Background

Small nucleolar RNAs (snoRNAs) are an abundant class of noncoding RNAs present in all eukaryotes and best known for their involvement in ribosome biogenesis. In mammalian genomes, many snoRNAs exist in multiple copies, resulting from recombination and retrotransposition from an ancestral snoRNA. To gain insight into snoRNA copy regulation, we used Rfam classification and normal human tissue expression datasets generated using low structure bias RNA-seq to characterize snoRNA families.

Results

We found that although box H/ACA families are on average larger than box C/D families, the number of expressed members is similar for both types. Family members can cover a wide range of average abundance values, but importantly, expression variability of individual members of a family is preferred over the total variability of the family, especially for box H/ACA snoRNAs, suggesting that while members are likely differentially regulated, mechanisms exist to ensure uniformity of the total family abundance across tissues. Box C/D snoRNA family members are mostly embedded in the same host gene while box H/ACA family members tend to be encoded in more than one different host, supporting a model in which box C/D snoRNA duplication occurred mostly by cis recombination while box H/ACA snoRNA families have gained copy members through retrotransposition. And unexpectedly, snoRNAs encoded in the same host gene can be regulated independently, as some snoRNAs within the same family vary in abundance in a divergent way between tissues.

Conclusions

SnoRNA copy regulation affects family sizes, genomic location of the members and controls simultaneously member and total family abundance to respond to the needs of individual tissues.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07757-1.

The Integrator Complex in Transcription and Development

The Integrator complex is conserved across metazoans and controls the fate of many nascent RNAs transcribed by RNA polymerase II (RNAPII). Among the 14 subunits of Integrator is an RNA endonuclease that is crucial for the biogenesis of small nuclear RNAs and enhancer RNAs. Integrator is further employed to trigger premature transcription termination at many protein-coding genes, thereby attenuating gene expression. Integrator thus helps to shape the transcriptome and ensure that genes can be robustly induced when needed. The molecular functions of Integrator subunits beyond the RNA endonuclease remain poorly understood, but some can act independently of the multisubunit complex. We highlight recent molecular insights into Integrator and propose how misregulation of this complex may lead to developmental defects and disease.

The cellular landscape of mid-size noncoding RNA

Noncoding RNA plays an important role in all aspects of the cellular life cycle, from the very basic process of protein synthesis to specialized roles in cell development and differentiation. However, many noncoding RNAs remain uncharacterized and the function of most of them remains unknown. Mid-size noncoding RNAs (mncRNAs), which range in length from 50 to 400 nucleotides, have diverse regulatory functions but share many fundamental characteristics. Most mncRNAs are produced from independent promoters although others are produced from the introns of other genes. Many are found in multiple copies in genomes. mncRNAs are highly structured and carry many posttranscriptional modifications. Both of these facets dictate their RNA-binding protein partners and ultimately their function. mncRNAs have already been implicated in translation, catalysis, as guides for RNA modification, as spliceosome components and regulatory RNA. However, recent studies are adding new mncRNA functions including regulation of gene expression and alternative splicing. In this review, we describe the different classes, characteristics and emerging functions of mncRNAs and their relative expression patterns. Finally, we provide a portrait of the challenges facing their detection and annotation in databases. This article is categorized under: Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs RNA Structure and Dynamics > RNA Structure, Dynamics, and Chemistry RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems RNA Evolution and Genomics > RNA and Ribonucleoprotein Evolution.

The Integrator complex regulates differential snRNA processing and fate of adult stem cells in the highly regenerative planarian Schmidtea mediterranea

In multicellular organisms, cell type diversity and fate depend on specific sets of transcript isoforms generated by post-transcriptional RNA processing. Here, we used Schmidtea mediterranea, a flatworm with extraordinary regenerative abilities and a large pool of adult stem cells, as an in vivo model to study the role of Uridyl-rich small nuclear RNAs (UsnRNAs), which participate in multiple RNA processing reactions including splicing, in stem cell regulation. We characterized the planarian UsnRNA repertoire, identified stem cell-enriched variants and obtained strong evidence for an increased rate of UsnRNA 3'-processing in stem cells compared to their differentiated counterparts. Consistently, components of the Integrator complex showed stem cell-enriched expression and their depletion by RNAi disrupted UsnRNA processing resulting in global changes of splicing patterns and reduced processing of histone mRNAs. Interestingly, loss of Integrator complex function disrupted both stem cell maintenance and regeneration of tissues. Our data show that the function of the Integrator complex in UsnRNA 3'-processing is conserved in planarians and essential for maintaining their stem cell pool. We propose that cell type-specific modulation of UsnRNA composition and maturation contributes to in vivo cell fate choices, such as stem cell self-renewal in planarians.

PARIS: Psoralen Analysis of RNA Interactions and Structures with High Throughput and Resolution

RNA has the intrinsic propensity to form base pairs, leading to complex intramolecular and intermolecular helices. Direct measurement of base pairing interactions in living cells is critical to solving transcriptome structure and interactions, and investigating their functions (Lu and Chang, Curr Opin Struct Biol 36:142-148, 2016). Toward this goal, we developed an experimental method, PARIS (Psoralen Analysis of RNA Interactions and Structures), to directly determine transcriptome-wide base pairing interactions (Lu et al., Cell 165(5):1267-1279, 2016). PARIS combines four critical steps, in vivo cross-linking, 2D gel purification, proximity ligation, and high-throughput sequencing to achieve high-throughput and near-base pair resolution determination of the RNA structurome and interactome in living cells. In this chapter, we aim to provide a comprehensive discussion on the principles behind the experimental and computational strategies, and a step-by-step description of the experiment and analysis.

Functional Interplay between Small Non-Coding RNAs and RNA Modification in the Brain

Small non-coding RNAs are essential for transcription, translation and gene regulation in all cell types, but are particularly important in neurons, with known roles in neurodevelopment, neuroplasticity and neurological disease. Many small non-coding RNAs are directly involved in the post-transcriptional modification of other RNA species, while others are themselves substrates for modification, or are functionally modulated by modification of their target RNAs. In this review, we explore the known and potential functions of several distinct classes of small non-coding RNAs in the mammalian brain, focusing on the newly recognised interplay between the epitranscriptome and the activity of small RNAs. We discuss the potential for this relationship to influence the spatial and temporal dynamics of gene activation in the brain, and predict that further research in the field of epitranscriptomics will identify interactions between small RNAs and RNA modifications which are essential for higher order brain functions such as learning and memory.

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.

Integrative analysis unveils new functions for the Drosophila Cutoff protein in noncoding RNA biogenesis and gene regulation

Piwi-interacting RNAs (piRNAs) are central components of the piRNA pathway, which directs transposon silencing and guarantees genome integrity in the germ cells of several metazoans. In Drosophila, piRNAs are produced from discrete regions of the genome termed piRNA clusters, whose expression relies on the RDC complex comprised of the core proteins Rhino, Deadlock, and Cutoff. To date, the RDC complex has been exclusively implicated in the regulation of the piRNA loci. Here we further elucidate the function of Cutoff and the RDC complex by performing genome-wide ChIP-seq and RNA-seq assays in the Drosophila ovaries and analyzing these data together with other publicly available data sets. In agreement with previous studies, we confirm that Cutoff is involved in the transcriptional regulation of piRNA clusters and in the repression of transposable elements in germ cells. Surprisingly, however, we find that Cutoff is enriched at and affects the expression of other noncoding RNAs, including spliceosomal RNAs (snRNAs) and small nucleolar RNAs (snoRNAs). At least in some instances, Cutoff appears to act at a transcriptional level in concert with Rhino and perhaps Deadlock. Finally, we show that mutations in Cutoff result in the deregulation of hundreds of protein-coding genes in germ cells. Our study uncovers a broader function for the RDC complex in the Drosophila germline development.

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.

Transcriptomic comparison of Drosophila snRNP biogenesis mutants reveals mutant-specific changes in pre-mRNA processing: implications for spinal muscular atrophy

Survival motor neuron (SMN) functions in the assembly of spliceosomal small nuclear ribonucleoproteins (snRNPs) that catalyze pre-mRNA splicing. Here, we used disruptions in Smn and two additional snRNP biogenesis genes, Phax and Ars2, to classify RNA processing differences as snRNP-dependent or gene-specific in Drosophila Phax and Smn mutants exhibited comparable reductions in snRNAs, and comparison of their transcriptomes uncovered shared sets of RNA processing changes. In contrast, Ars2 mutants displayed only small decreases in snRNA levels, and RNA processing changes in these mutants were generally distinct from those identified in Phax and Smn animals. Instead, RNA processing changes in Ars2 mutants support the known interaction of Ars2 protein with the cap-binding complex, as splicing changes showed a clear bias toward the first intron. Bypassing disruptions in snRNP biogenesis, direct knockdown of spliceosomal proteins caused similar changes in the splicing of snRNP-dependent events. However, these snRNP-dependent events were largely unaltered in three Smn mutants expressing missense mutations that were originally identified in human spinal muscular atrophy (SMA) patients. Hence, findings here clarify the contributions of Phax, Smn, and Ars2 to snRNP biogenesis in Drosophila, and loss-of-function mutants for these proteins reveal differences that help disentangle cause and effect in SMA model flies.