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

Protocol for siRNA-mediated U1 snRNA knockdown using fluorinated α-helical polypeptide in vitro and in vivo

Here, we present a protocol for small interfering RNA (siRNA)-mediated U1 small nuclear RNA (snRNA) knockdown using fluorinated α-helical polypeptide in macrophages and mouse lungs, providing a dependable approach to silence U1 snRNA in vitro and in vivo. We describe steps for preparing P7F7/siRNA polyplexes and silencing U1 snRNA with P7F7/siRNA polyplexes in macrophages and mouse lungs. Knockdown efficiency is validated through reverse-transcription quantitative real-time PCR analysis. This protocol is applicable for studying the physiological or pathophysiological function of U1 snRNA. For complete details on the use and execution of this protocol, please refer to Zhang et al.1.

U1 snRNA interactions with deep intronic sequences regulate splicing of multiple exons of spinal muscular atrophy genes

Introduction

The U1 small nuclear RNA (snRNA) forms ribonucleoprotein particles (RNPs) such as U1 snRNP and U1-TAF15 snRNP. U1 snRNP is one of the most studied RNPs due to its critical role in pre-mRNA splicing in defining the 5' splice site (5'ss) of every exon through direct interactions with sequences at exon/intron junctions. Recent reports support the role of U1 snRNP in all steps of transcription, namely initiation, elongation, and termination. Functions of U1-TAF15 snRNP are less understood, though it associates with the transcription machinery and may modulate pre-mRNA splicing by interacting with the 5'ss and/or 5'ss-like sequences within the pre-mRNA. An anti-U1 antisense oligonucleotide (ASO) that sequesters the 5' end of U1 snRNA inhibits the functions of U1 snRNP, including transcription and splicing. However, it is not known if the inhibition of U1 snRNP influences post-transcriptional regulation of pre-mRNA splicing through deep intronic sequences.

Methods

We examined the effect of an anti-U1 ASO that sequesters the 5' end of U1 snRNA on transcription and splicing of all internal exons of the spinal muscular atrophy (SMA) genes, SMN1 and SMN2. Our study was enabled by the employment of a multi-exon-skipping detection assay (MESDA) that discriminates against prematurely terminated transcripts. We employed an SMN2 super minigene to determine if anti-U1 ASO differently affects splicing in the context of truncated introns.

Results

We observed substantial skipping of multiple internal exons of SMN1 and SMN2 triggered by anti-U1 treatment. Suggesting a role for U1 snRNP in interacting with deep intronic sequences, early exons of the SMN2 super minigene with truncated introns were resistant to anti-U1 induced skipping. Consistently, overexpression of engineered U1 snRNAs targeting the 5'ss of early SMN1 and SMN2 exons did not prevent exon skipping caused by anti-U1 treatment.

Discussion

Our results uncover a unique role of the U1 snRNA-associated RNPs in splicing regulation executed through deep intronic sequences. Findings are significant for developing novel therapies for SMA based on deep intronic targets.

Splicing analysis of STAT3 tandem donor suggests non-canonical binding registers for U1 and U6 snRNAs

Tandem donor splice sites (5'ss) are unique regions with at least two GU dinucleotides serving as splicing cleavage sites. The Δ3 tandem 5'ss are a specific subclass of 5'ss separated by 3 nucleotides which can affect protein function by inserting/deleting a single amino acid. One 5'ss is typically preferred, yet factors governing particular 5'ss choice are not fully understood. A highly conserved exon 21 of the STAT3 gene was chosen as a model to study Δ3 tandem 5'ss splicing mechanisms. Based on multiple lines of experimental evidence, endogenous U1 snRNA most likely binds only to the upstream 5'ss. However, the downstream 5'ss is used preferentially, and the splice site choice is not dependent on the exact U1 snRNA binding position. Downstream 5'ss usage was sensitive to exact nucleotide composition and dependent on the presence of downstream regulatory region. The downstream 5'ss usage could be best explained by two novel interactions with endogenous U6 snRNA. U6 snRNA enables the downstream 5'ss usage in STAT3 exon 21 by two mechanisms: (i) binding in a novel non-canonical register and (ii) establishing extended Watson-Crick base pairing with the downstream regulatory region. This study suggests that U6:5'ss interaction is more flexible than previously thought.

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.

Sm complex assembly and 5' cap trimethylation promote selective processing of snRNAs by the 3' exonuclease TOE1

Competing exonucleases that promote 3' end maturation or degradation direct quality control of small non-coding RNAs, but how these enzymes distinguish normal from aberrant RNAs is poorly understood. The Pontocerebellar Hypoplasia 7 (PCH7)-associated 3' exonuclease TOE1 promotes maturation of canonical small nuclear RNAs (snRNAs). Here, we demonstrate that TOE1 achieves specificity toward canonical snRNAs through their Sm complex assembly and cap trimethylation, two features that distinguish snRNAs undergoing correct biogenesis from other small non-coding RNAs. Indeed, disruption of Sm complex assembly via snRNA mutations or protein depletions obstructs snRNA processing by TOE1, and in vitro snRNA processing by TOE1 is stimulated by a trimethylated cap. An unstable snRNA variant that normally fails to undergo maturation becomes fully processed by TOE1 when its degenerate Sm binding motif is converted into a canonical one. Our findings uncover the molecular basis for how TOE1 distinguishes snRNAs from other small non-coding RNAs and explain how TOE1 promotes maturation specifically of canonical snRNAs undergoing proper processing.

The forgotten variable? Does the euthanasia method and sample storage condition influence an organisms transcriptome - a gene expression analysis on multiple tissues in pigs

BACKGROUND

Transcriptomic studies often require collection of fresh tissues post euthanasia. The chosen euthanasia method might have the potential to induce variations in gene expressions that are unlinked with the experimental design. The present study compared the suitability of 'nitrogen gas in foam' (ANOXIA) in comparison to a non-barbiturate anaesthetic, T-61® (T61), for euthanizing piglets used in transcriptome research. Further, the effect of common tissue storage conditions, RNAlater™ (RL) and snap freezing in liquid nitrogen (LN2), on gene expression profiles were also analysed.

RESULTS

On comparison of the 3'mRNA-Seq data generated from pituitary, hypothalamus, liver and lung tissues, no significant differential expression in the protein coding genes were detected between the euthanasia methods. This implies that the nitrogen anoxia method could be a suitable alternative for euthanasia of piglets used in transcriptomic research. However, small nuclear RNAs (snRNAs) that constitute the eukaryotic spliceosomal machinery were found to be significantly higher (log2fold change ≥ 2.0, and adjusted p value ≤ 0.1) in pituitary samples collected using ANOXIA. Non-protein coding genes like snRNAs that play an important role in pre-mRNA splicing can subsequently modify gene expression. Storage in RL was found to be superior in preserving RNA compared to LN2 storage, as evidenced by the significantly higher RIN values in representative samples. However, storage in RL as opposed to LN2, also influenced differential gene expression in multiple tissues, perhaps as a result of its inability to inhibit biological activity during storage. Hence such external sources of variations should be carefully considered before arriving at research conclusions.

CONCLUSIONS

Source of biological variations like euthanasia method and storage condition can confound research findings. Even if we are unable to prevent the effect of these external factors, it will be useful to identify the impact of these variables on the parameter under observation and thereby prevent misinterpretation of our results.

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.

Non-small Cell Lung Cancer Epigenomes Exhibit Altered DNA Methylation in Smokers and Never-smokers

Epigenetic alterations are widespread in cancer and can complement genetic alterations to influence cancer progression and treatment outcome. To determine the potential contribution of DNAmethylation alterations to tumor phenotype in non-small cell lung cancer (NSCLC) in both smoker and never-smoker patients, we performed genome-wide profiling of DNA methylation in 17 primary NSCLC tumors and 10 matched normal lung samples using the complementary assays, methylated DNA immunoprecipitation sequencing (MeDIP-seq) and methylation sensitive restriction enzyme sequencing (MRE-seq). We reported recurrent methylation changes in the promoters of several genes, many previously implicated in cancer, including FAM83A and SEPT9 (hypomethylation), as well as PCDH7, NKX2-1, and SOX17 (hypermethylation). Although many methylation changes between tumors and their paired normal samples were shared across patients, several were specific to a particular smoking status. For example, never-smokers displayed a greater proportion of hypomethylated differentially methylated regions (hypoDMRs) and a greater number of recurrently hypomethylated promoters, including those of ASPSCR1, TOP2A, DPP9, and USP39, all previously linked to cancer. Changes outside of promoters were also widespread and often recurrent, particularly methylation loss over repetitive elements, highly enriched for ERV1 subfamilies. Recurrent hypoDMRs were enriched for several transcription factor binding motifs, often for genes involved in signaling and cell proliferation. For example, 71% of recurrent promoter hypoDMRs contained a motif for NKX2-1. Finally, the majority of DMRs were located within an active chromatin state in tissues profiled using the Roadmap Epigenomics data, suggesting that methylation changes may contribute to altered regulatory programs through the adaptation of cell type-specific expression programs.

Development of Engineered-U1 snRNA Therapies: Current Status

Splicing of pre-mRNA is a crucial regulatory stage in the pathway of gene expression. The majority of human genes that encode proteins undergo alternative pre-mRNA splicing and mutations that affect splicing are more prevalent than previously thought. Targeting aberrant RNA(s) may thus provide an opportunity to correct faulty splicing and potentially treat numerous genetic disorders. To that purpose, the use of engineered U1 snRNA (either modified U1 snRNAs or exon-specific U1s-ExSpeU1s) has been applied as a potentially therapeutic strategy to correct splicing mutations, particularly those affecting the 5' splice-site (5'ss). Here we review and summarize a vast panoply of studies that used either modified U1 snRNAs or ExSpeU1s to mediate gene therapeutic correction of splicing defects underlying a considerable number of genetic diseases. We also focus on the pre-clinical validation of these therapeutic approaches both in vitro and in vivo, and summarize the main obstacles that need to be overcome to allow for their successful translation to clinic practice in the future.

The 3' exonuclease TOE1 selectively processes snRNAs through recognition of Sm complex assembly and 5' cap trimethylation

Competing exonucleases that promote 3' end maturation or degradation direct quality control of small non-coding RNAs, but how these enzymes distinguish normal from aberrant RNAs is poorly understood. The Pontocerebellar Hypoplasia 7 (PCH7)-associated 3' exonuclease TOE1 promotes maturation of canonical small nuclear RNAs (snRNAs). Here, we demonstrate that TOE1 achieves specificity towards canonical snRNAs by recognizing Sm complex assembly and cap trimethylation, two features that distinguish snRNAs undergoing correct biogenesis from other small non-coding RNAs. Indeed, disruption of Sm complex assembly via snRNA mutations or protein depletions obstructs snRNA processing by TOE1, and in vitro snRNA processing by TOE1 is stimulated by a trimethylated cap. An unstable snRNA variant that normally fails to undergo maturation becomes fully processed by TOE1 when its degenerate Sm binding motif is converted into a canonical one. Our findings uncover the molecular basis for how TOE1 distinguishes snRNAs from other small non-coding RNAs and explain how TOE1 promotes maturation specifically of canonical snRNAs undergoing proper processing.

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.

Spliceosome assembly and regulation: insights from analysis of highly reduced spliceosomes

Premessenger RNA splicing is catalyzed by the spliceosome, a multimegadalton RNA-protein complex that assembles in a highly regulated process on each intronic substrate. Most studies of splicing and spliceosomes have been carried out in human or S. cerevisiae model systems. There exists, however, a large diversity of spliceosomes, particularly in organisms with reduced genomes, that suggests a means of analyzing the essential elements of spliceosome assembly and regulation. In this review, we characterize changes in spliceosome composition across phyla, describing those that are most frequently observed and highlighting an analysis of the reduced spliceosome of the red alga Cyanidioschyzon merolae We used homology modeling to predict what effect splicing protein loss would have on the spliceosome, based on currently available cryo-EM structures. We observe strongly correlated loss of proteins that function in the same process, for example, in interacting with the U1 snRNP (which is absent in C. merolae), regulation of Brr2, or coupling transcription and splicing. Based on our observations, we predict splicing in C. merolae to be inefficient, inaccurate, and post-transcriptional, consistent with the apparent trend toward its elimination in this lineage. This work highlights the striking flexibility of the splicing pathway and the spliceosome when viewed in the context of eukaryotic diversity.

Adaptation of the Porcine Pituitary Transcriptome, Spliceosome and Editome during Early Pregnancy

The physiological mechanisms of the porcine reproduction are relatively well-known. However, transcriptomic changes and the mechanisms accompanying transcription and translation processes in various reproductive organs, as well as their dependence on hormonal status, are still poorly understood. The aim of this study was to gain a principal understanding of alterations within the transcriptome, spliceosome and editome occurring in the pituitary of the domestic pig (Sus scrofa domestica L.), which controls basic physiological processes in the reproductive system. In this investigation, we performed extensive analyses of data obtained by high-throughput sequencing of RNA from the gilts' pituitary anterior lobes during embryo implantation and the mid-luteal phase of the estrous cycle. During analyses, we obtained detailed information on expression changes of 147 genes and 43 long noncoding RNAs, observed 784 alternative splicing events and also found the occurrence of 8729 allele-specific expression sites and 122 RNA editing events. The expression profiles of the selected 16 phenomena were confirmed by PCR or qPCR techniques. As a final result of functional meta-analysis, we acquired knowledge regarding intracellular pathways that induce changes in the processes accompanying transcription and translation regulation, which may induce modifications in the secretory activity of the porcine adenohypophyseal cells.

A method for simultaneous detection of small and long RNA biotypes by ribodepleted RNA-Seq

RNA sequencing offers unprecedented access to the transcriptome. Key to this is the identification and quantification of many different species of RNA from the same sample at the same time. In this study we describe a novel protocol for simultaneous detection of coding and non-coding transcripts using modifications to the Ion Total RNA-Seq kit v2 protocol, with integration of QIASeq FastSelect rRNA removal kit. We report highly consistent sequencing libraries can be produced from both frozen high integrity mouse hippocampal tissue and the more challenging post-mortem human tissue. Removal of rRNA using FastSelect was extremely efficient, resulting in less than 1.5% rRNA content in the final library. We identified > 30,000 unique transcripts from all samples, including protein-coding genes and many species of non-coding RNA, in biologically-relevant proportions. Furthermore, the normalized sequencing read count for select genes significantly negatively correlated with Ct values from qRT-PCR analysis from the same samples. These results indicate that this protocol accurately and consistently identifies and quantifies a wide variety of transcripts simultaneously. The highly efficient rRNA depletion, coupled with minimized sample handling and without complicated and high-loss size selection protocols, makes this protocol useful to researchers wishing to investigate whole transcriptomes.

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.

An alternative D. melanogaster 7SK snRNP

BACKGROUND

The 7SK small nuclear RNA (snRNA) found in most metazoans is a key regulator of P-TEFb which in turn regulates RNA polymerase II elongation. Although its primary sequence varies in protostomes, its secondary structure and function are conserved across evolutionary distant taxa.

RESULTS

Here, we describe a novel ncRNA sharing many features characteristic of 7SK RNAs, in D. melanogaster. We examined the structure of the corresponding gene and determined the expression profiles of the encoded RNA, called snRNA:7SK:94F, during development. It is probably produced from the transcription of a lncRNA which is processed into a mature snRNA. We also addressed its biological function and we show that, like dm7SK, this alternative 7SK interacts in vivo with the different partners of the P-TEFb complex, i.e. HEXIM, LARP7 and Cyclin T. This novel RNA is widely expressed across tissues.

CONCLUSION

We propose that two distinct 7SK genes might contribute to the formation of the 7SK snRNP complex in D. melanogaster.

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 cap-specific m6A methyltransferase, PCIF1/CAPAM, is dynamically recruited to the gene promoter in a transcription-dependent manner

N 6-methyladenosine (m6A), the most abundant modification in eukaryotic mRNAs, plays an important role in mRNA metabolism and functions. When adenosine is transcribed as the first cap-adjacent nucleotide, it is methylated at the ribose 2'-O and N6 positions, thus generating N6, 2'-O-dimethyladenosine (m6Am). Phosphorylated C-terminal domain (CTD)-interacting factor 1 (PCIF1) is a novel cap-specific adenine N6-methyltransferase responsible for m6Am formation. As PCIF1 specifically interacts with the Ser5-phosphorylated CTD of RNA polymerase II (Pol II), which is a marker for the early phase of transcription, PCIF1 is speculated to be recruited to the early elongating Pol II. In this study, subcellular fractionation and immunofluorescence microscopy demonstrated that PCIF1 is mainly localized to the transcriptionally active chromatin regions in HeLa cells. Chromatin immunoprecipitation (ChIP) revealed that PCIF1 was predominantly localized to the promoter of a broad range of Pol II-transcribed genes, including several protein-coding genes and non-coding RNA genes. Moreover, PCIF1 accumulation on these promoters depended entirely on transcriptional activity and Ser5 phosphorylation of the CTD. These results suggest that PCIF1 dynamically localizes to the Pol II early in transcription and may efficiently catalyze N6-methylation of the first adenosine residue of nascent mRNAs cotranscriptionally.

U1 snRNP telescripting: molecular mechanisms and beyond

U1 snRNP is one of the most abundant ribonucleoprotein (RNP) complexes in eukaryotic cells and is estimated to be approximately 1 million copies per cell. Apart from its canonical role in mRNA splicing, this complex has emerged as a key regulator of eukaryotic mRNA length via inhibition of mRNA 3'-end processing at numerous intronic polyadenylation sites, in a process that is also termed 'U1 snRNP telescripting'. Several reviews have extensively described the concept of U1 telescripting and subsequently highlighted its potential impacts in mRNA metabolism. Here, we review what is currently known regarding the underlying mechanisms of this important phenomenon and discuss open questions and future challenges.

Alcohol use disorder causes global changes in splicing in the human brain

Alcohol use disorder (AUD) is a widespread disease leading to the deterioration of cognitive and other functions. Mechanisms by which alcohol affects the brain are not fully elucidated. Splicing constitutes a nuclear process of RNA maturation, which results in the formation of the transcriptome. We tested the hypothesis as to whether AUD impairs splicing in the superior frontal cortex (SFC), nucleus accumbens (NA), basolateral amygdala (BLA), and central nucleus of the amygdala (CNA). To evaluate splicing, bam files from STAR alignments were indexed with samtools for use by rMATS software. Computational analysis of affected pathways was performed using Gene Ontology Consortium, Gene Set Enrichment Analysis, and LncRNA Ontology databases. Surprisingly, AUD was associated with limited changes in the transcriptome: expression of 23 genes was altered in SFC, 14 in NA, 102 in BLA, and 57 in CNA. However, strikingly, mis-splicing in AUD was profound: 1421 mis-splicing events were detected in SFC, 394 in NA, 1317 in BLA, and 469 in CNA. To determine the mechanism of mis-splicing, we analyzed the elements of the spliceosome: small nuclear RNAs (snRNAs) and splicing factors. While snRNAs were not affected by alcohol, expression of splicing factor heat shock protein family A (Hsp70) member 6 (HSPA6) was drastically increased in SFC, BLA, and CNA. Also, AUD was accompanied by aberrant expression of long noncoding RNAs (lncRNAs) related to splicing. In summary, alcohol is associated with genome-wide changes in splicing in multiple human brain regions, likely due to dysregulation of splicing factor(s) and/or altered expression of splicing-related lncRNAs.

Derivation of stable embryonic stem cell-like, but transcriptionally heterogenous, induced pluripotent stem cells from non-permissive mouse strains

Genetic background is known to play a role in the ability to derive pluripotent, embryonic stem cells (ESC), a trait referred to as permissiveness. Previously we demonstrated that induced pluripotent stem cells (iPSC) can be readily derived from non-permissive mouse strains by addition of serum-based media supplemented with GSK3B and MEK inhibitors, termed 2iS media, 3 days into reprogramming. Here, we describe the derivation of second type of iPSC colony from non-permissive mouse strains that can be stably maintained independently of 2iS media. The resulting cells display transcriptional heterogeneity similar to that observed in ESC from permissive genetic backgrounds derived in conventional serum containing media supplemented with leukemia inhibitor factor. However, unlike previous studies that report exclusive subpopulations, we observe both exclusive and simultaneous expression of naive and primed cell surface markers. Herein, we explore shifts in pluripotency in the presence of 2iS and characterize heterogenous subpopulations to determine their pluripotent state and role in heterogenous iPSCs derived from the non-permissive NOD/ShiLtJ strain. We conclude that heterogeneity is a naturally occurring, necessary quality of stem cells that allows for the maintenance of pluripotency. This study further demonstrates the efficacy of the 2iS reprogramming technique. It is also the first study to derive stable ESC-like stem cells from the non-permissive NOD/ShiLtJ and WSB/EiJ strains, enabling easier and broader research possibilities into pluripotency for these and similar non-permissive mouse strains and species.

ISG20 and nuclear exosome promote destabilization of nascent transcripts for spliceosomal U snRNAs and U1 variants

Primary RNA transcripts are processed in a plethora of ways to become mature functional forms. In one example, human spliceosomal U snRNAs are matured at their 3'-end by an exonuclease termed TOE1. This process is important because mutations in TOE1 gene can cause a human genetic disease, pontocerebellar hypoplasia (PCH). Nevertheless, TOE1 may not be the only maturation exonuclease for U snRNAs in the cell. Here, we biochemically identify two exonucleolytic factors, Interferon-stimulated gene 20-kDa protein (ISG20) and the nuclear exosome as such candidates, using a newly developed in vitro system that recapitulates 3'-end maturation of U1 snRNA. However, extensive 3'-end sequencing of endogenous U1 snRNA of the knockdown (KD) cells revealed that these factors are not the maturation factors per se. Instead, the nascent transcripts of the spliceosomal U snRNAs as well as of unstable U1 variants were found to increase in quantity upon KD of the factors. These results indicated that ISG20 and the nuclear exosome promote the degradation of nascent spliceosomal U snRNAs and U1 variants, and therefore implied their role in the quality control of newly synthesized U snRNAs.

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.

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.

5' and 3' splicing signals evolution in vertebrates: Analysis in a conserved gene family

The mitochondrial solute carrier genes (SLC25) are highly conserved during vertebrate evolution. In most SLC25 genes of zebrafish, chicken, mouse, and human, the introns are located at exactly superimposable positions. In these topographically corresponding introns we studied the composition of the initial and terminal hexanucleotides (5'ss and 3'ss) which are instrumental in splicing signaling, focusing on the evolutionary conservation/mutation dynamics of these genetically related sequences. At each position, the per cent conservation of zebrafish individual nucleotides in chicken, mouse and human is proportional to their percent frequency in zebrafish; furthermore, nucleotide mutations are biased in favor of the more represented nucleotides, thus compensating for those highly represented zebrafish nucleotides which have not been conserved. As a result of these evolutionary dynamics, the general nucleotide composition at each position has remained relatively conserved throughout vertebrates. At 5'ss, following the canonical GT, A and G are largely prevailing at position +3, A at +4 and G at +5 (GT[A/G]AGx). At 3'ss, T and C are largely prevailing at positions -6, -5 and -3, preceding the canonical intron terminal AG ([C/T] [C/T]x[C/T]AG). However, the actual composition of the tetranucleotides at 5' and 3' often does not conform to the above scheme. At 5'ss the more canonical sequence is completely expressed in 63% of cases and partially (2 or 1 matches) in 37 % of cases. At 3'ss the more canonical sequence is completely expressed in 71 % of cases and partially (2 or 1 matches) in 29 % of cases. The nucleotide conservation loss (nucleotide mutation) is higher in the evolution from fish to the last common ancestor of birds and mammals (58 %), then diminishes in the successive evolution steps up to the mammalian common ancestor (10 %), and becomes still lower at the divergence of rodents and primates (5 %).

A guide to naming human non-coding RNA genes

Research on non-coding RNA (ncRNA) is a rapidly expanding field. Providing an official gene symbol and name to ncRNA genes brings order to otherwise potential chaos as it allows unambiguous communication about each gene. The HUGO Gene Nomenclature Committee (HGNC, www.genenames.org) is the only group with the authority to approve symbols for human genes. The HGNC works with specialist advisors for different classes of ncRNA to ensure that ncRNA nomenclature is accurate and informative, where possible. Here, we review each major class of ncRNA that is currently annotated in the human genome and describe how each class is assigned a standardised nomenclature.

Genome-wide evolution of wobble base-pairing nucleotides of branchpoint motifs with increasing organismal complexity

How have the branchpoint motifs evolved in organisms of different complexity? Here we identified and examined the consensus motifs (R₁C₂T₃R₄A₅Y₆, R: A or G, Y: C or T) of 898 fungal genomes. In Ascomycota unicellular yeasts, the G₄/A₄ ratio is mostly (98%) below 0.125 but increases sharply in multicellular species by about 40 times on average, and in the more complex Basidiomycota, it increases further by about 7 times. The global G₄ increase is consistent with A₄ to G₄ transitions in evolution. Of the G₄/A₄-interacting amino acids of the branchpoint binding protein MSL5 (SF1) and the HSH155 (SF3B1), as well as the 5' splice sites (SS) and U2 snRNA genes, the 5' SS G³/A³ co-vary with the G₄ to some extent. However, corresponding increase of the G₄-complementary GCAGTA-U2 gene is rare, suggesting wobble-base pairing between the G₄-containing branchpoint motif and GTAGTA-U2 in most of these species. Interestingly, the G₄/A₄ ratio correlates well with the abundance of alternative splicing in the two phyla, and G₄ enriched significantly at the alternative 3' SS of genes in RNA metabolism, kinases and membrane proteins. Similar wobble nucleotides also enriched at the 3' SS of multicellular fungi with only thousands of protein-coding genes. Thus, branchpoint motifs have evolved U2-complementarity in unicellular Ascomycota yeasts, but have gradually gained more wobble base-pairing nucleotides in fungi of higher complexity, likely to destabilize branchpoint motif-U2 interaction and/or branchpoint A protrusion for alternative splicing. This implies an important role of relaxing the branchpoint signals in the multicellularity and further complexity of fungi.

SCAF4 and SCAF8, mRNA Anti-Terminator Proteins

Accurate regulation of mRNA termination is required for correct gene expression. Here, we describe a role for SCAF4 and SCAF8 as anti-terminators, suppressing the use of early, alternative polyadenylation (polyA) sites. The SCAF4/8 proteins bind the hyper-phosphorylated RNAPII C-terminal repeat domain (CTD) phosphorylated on both Ser2 and Ser5 and are detected at early, alternative polyA sites. Concomitant knockout of human SCAF4 and SCAF8 results in altered polyA selection and subsequent early termination, leading to expression of truncated mRNAs and proteins lacking functional domains and is cell lethal. While SCAF4 and SCAF8 work redundantly to suppress early mRNA termination, they also have independent, non-essential functions. SCAF8 is an RNAPII elongation factor, whereas SCAF4 is required for correct termination at canonical, distal transcription termination sites in the presence of SCAF8. Together, SCAF4 and SCAF8 coordinate the transition between elongation and termination, ensuring correct polyA site selection and RNAPII transcriptional termination in human cells.

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.

U1 snRNP Telescripting: Suppression of Premature Transcription Termination in Introns as a New Layer of Gene Regulation

Recent observations showed that nascent RNA polymerase II transcripts, pre-mRNAs, and noncoding RNAs are highly susceptible to premature 3'-end cleavage and polyadenylation (PCPA) from numerous intronic cryptic polyadenylation signals (PASs). The importance of this in gene regulation was not previously appreciated as PASs, despite their prevalence, were thought to be active in terminal exons at gene ends. Unexpectedly, antisense oligonucleotide interference with U1 snRNA base-pairing to 5' splice sites, which is necessary for U1 snRNP's (U1) function in splicing, caused widespread PCPA in metazoans. This uncovered U1's PCPA suppression activity, termed telescripting, as crucial for full-length transcription in thousands of vertebrate genes, providing a general role in transcription elongation control. Progressive intron-size expansion in metazoan evolution greatly increased PCPA vulnerability and dependence on U1 telescripting. We describe how these observations unfolded and discuss U1 telescripting's role in shaping the transcriptome.

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.

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.

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.

PHYSIOLOGY AND ENDOCRINOLOGY SYMPOSIUM:The effects of poor maternal nutrition during gestation on offspring postnatal growth and metabolism

Poor maternal nutrition during gestation has been linked to poor growth and development, metabolic dysfunction, impaired health, and reduced productivity of offspring in many species. Poor maternal nutrition can be defined as an excess or restriction of overall nutrients or specific macro- or micronutrients in the diet of the mother during gestation. Interestingly, there are several reports that both restricted- and over-feeding during gestation negatively affect offspring postnatal growth with reduced muscle and bone deposition, increased adipose accumulation, and metabolic dysregulation through reduced leptin and insulin sensitivity. Our laboratory and others have used experimental models of restricted- and over-feeding during gestation to evaluate effects on early postnatal growth of offspring. Restricted- and over-feeding during gestation alters body size, circulating growth factors, and metabolic hormones in offspring postnatally. Both restricted- and over-feeding alter muscle growth, increase lipid content in the muscle, and cause changes in expression of myogenic factors. Although the negative effects of poor maternal nutrition on offspring growth have been well characterized in recent years, the mechanisms contributing to these changes are not well established. Our laboratory has focused on elucidating these mechanisms by evaluating changes in gene and protein expression, and stem cell function. Through RNA-Seq analysis, we observed changes in expression of genes involved in protein synthesis, metabolism, cell function, and signal transduction in muscle tissue. We recently reported that satellite cells, muscle stem cells, have altered expression of myogenic factors in offspring from restricted-fed mothers. Bone marrow derived mesenchymal stem cells, multipotent cells that contribute to development and maintenance of several tissues including bone, muscle, and adipose, have a 50% reduction in cell proliferation and altered metabolism in offspring from both restricted- and over-fed mothers. These findings indicate that poor maternal nutrition may alter offspring postnatal growth by programming stem cell populations. In conclusion, poor maternal nutrition during gestation negatively affects offspring postnatal growth, potentially through impaired stem and satellite cell function. Therefore, determining the mechanisms that contribute to fetal programming is critical to identifying effective management interventions for these offspring and improving efficiency of production.

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.

Small Non-coding RNA Abundance in Adrenocortical Carcinoma: A Footprint of a Rare Cancer

BACKGROUND: Adrenocortical carcinoma (ACC) is a relatively rare, but aggressive type of cancer, which affects both children and adults. OBJECTIVE: Small non-coding RNAs (sncRNAs) play important roles and may serve as biomarkers for disease diagnosis, prognosis and treatment. METHODS: In our study, we sought to identify sncRNAs associated with malignant adrenal tumors. We obtained publicly available, small RNA sequencing data derived from 45 ACC and 30 benign tumors arising from the cortex of the adrenal gland, adrenocortical adenomas (ACA), and compared their sncRNA expression profiles. RESULTS: First, we remapped small RNA-seq to miRBase version 21 to check expression of miRNAs and found 147 miRNAs were aberrantly expressed (p<0.05) in ACC samples compared to ACA samples. Pathway analysis of differentially expressed miRNAs revealed p53 signaling pathways to be profoundly affected in ACC samples. Further examination for other types of small RNAs revealed 16 piRNAs, 48 lncRNAs and 19 sn/snoRNAs identified in ACC samples. Conclusions: Our data analysis suggests that publically available resources can be mined for biomarker development and improvements in-patient care; however, further research must be performed to correlate tumor grade with gene expression.

Th17 Polarization under Hypoxia Results in Increased IL-10 Production in a Pathogen-Independent Manner

The IL-17-producing CD4⁺ T helper cell (Th17) differentiation is affected by stimulation of the aryl hydrocarbon receptor (AhR) pathway and by hypoxia-inducible factor 1 alpha (HIF-1α). In some cases, Th17 become non-pathogenic and produce IL-10. However, the initiating events triggering this phenotype are yet to be fully understood. Here, we show that such cells may be differentiated at low oxygen and regardless of AhR ligand treatment such as cigarette smoke extract. Hypoxia led to marked alterations of the transcriptome of IL-10-producing Th17 cells affecting genes involved in metabolic, anti-apoptotic, cell cycle, and T cell functional pathways. Moreover, we show that oxygen regulates the expression of CD52, which is a cell surface protein that has been shown to suppress the activation of other T cells upon release. Taken together, these findings suggest a novel ability for Th17 cells to regulate immune responses in vivo in an oxygen-dependent fashion.

Integration of quantitated expression estimates from polyA-selected and rRNA-depleted RNA-seq libraries

BACKGROUND

The availability of fast alignment-free algorithms has greatly reduced the computational burden of RNA-seq processing, especially for relatively poorly assembled genomes. Using these approaches, previous RNA-seq datasets could potentially be processed and integrated with newly sequenced libraries. Confounding factors in such integration include sequencing depth and methods of RNA extraction and selection. Different selection methods (typically, either polyA-selection or rRNA-depletion) omit different RNAs, resulting in different fractions of the transcriptome being sequenced. In particular, rRNA-depleted libraries sample a broader fraction of the transcriptome than polyA-selected libraries. This study aimed to develop a systematic means of accounting for library type that allows data from these two methods to be compared.

RESULTS

The method was developed by comparing two RNA-seq datasets from ovine macrophages, identical except for RNA selection method. Gene-level expression estimates were obtained using a two-part process centred on the high-speed transcript quantification tool Kallisto. Firstly, a set of reference transcripts was defined that constitute a standardised RNA space, with expression from both datasets quantified against it. Secondly, a simple ratio-based correction was applied to the rRNA-depleted estimates. The outcome is an almost perfect correlation between gene expression estimates, independent of library type and across the full range of levels of expression.

CONCLUSION

A combination of reference transcriptome filtering and a ratio-based correction can create equivalent expression profiles from both polyA-selected and rRNA-depleted libraries. This approach will allow meta-analysis and integration of existing RNA-seq data into transcriptional atlas projects.

A Bioinformatics-Based Alternative mRNA Splicing Code that May Explain Some Disease Mutations Is Conserved in Animals

Deep sequencing of cDNAs made from spliced mRNAs indicates that most coding genes in many animals and plants have pre-mRNA transcripts that are alternatively spliced. In pre-mRNAs, in addition to invariant exons that are present in almost all mature mRNA products, there are at least 6 additional types of exons, such as exons from alternative promoters or with alternative polyA sites, mutually exclusive exons, skipped exons, or exons with alternative 5′ or 3′ splice sites. Our bioinformatics-based hypothesis is that, in analogy to the genetic code, there is an “alternative-splicing code” in introns and flanking exon sequences, analogous to the genetic code, that directs alternative splicing of many of the 36 types of introns. In humans, we identified 42 different consensus sequences that are each present in at least 100 human introns. 37 of the 42 top consensus sequences are significantly enriched or depleted in at least one of the 36 types of introns. We further supported our hypothesis by showing that 96 out of 96 analyzed human disease mutations that affect RNA splicing, and change alternative splicing from one class to another, can be partially explained by a mutation altering a consensus sequence from one type of intron to that of another type of intron. Some of the alternative splicing consensus sequences, and presumably their small-RNA or protein targets, are evolutionarily conserved from 50 plant to animal species. We also noticed the set of introns within a gene usually share the same splicing codes, thus arguing that one sub-type of splicesosome might process all (or most) of the introns in a given gene. Our work sheds new light on a possible mechanism for generating the tremendous diversity in protein structure by alternative splicing of pre-mRNAs.

A Comprehensive NGS Data Analysis of Differentially Regulated miRNAs, piRNAs, lncRNAs and sn/snoRNAs in Triple Negative Breast Cancer

Cancer is the second leading cause of death in the United States and is a major public health concern worldwide. Basic, clinical and epidemiological research is leading to improved cancer detection, prevention, and outcomes. Recent technological advances have allowed unbiased and comprehensive screening of genome-wide gene expression. Small non-coding RNAs (sncRNAs) have been shown to play an important role in biological processes and could serve as a diagnostic, prognostic and therapeutic biomarker for specific diseases. Recent findings have begun to reveal and enhance our understanding of the complex architecture of sncRNA expression including miRNAs, piRNAs, lncRNAs, sn/snoRNAs and their relationships with biological systems. We used publicly available small RNA sequencing data that was derived from 24 triple negative breast cancers (TNBC) and 14 adjacent normal tissue samples to remap various types of sncRNAs. We found a total of 55 miRNAs were aberrantly expressed (p<0.005) in TNBC samples (8 miRNAs upregulated; 47 downregulated) compared to adjacent normal tissues whereas the original study reported only 25 novel miRs. In this study, we used pathway analysis of differentially expressed miRNAs which revealed TGF-beta signaling pathways to be profoundly affected in the TNBC samples. Furthermore, our comprehensive re-mapping strategy allowed us to discover a number of other differentially expressed sncRNAs including piRNAs, lncRNAs, sn/snoRNAs, rRNAs, miscRNAs and nonsense-mediated decay RNAs. We believe that our sncRNA analysis workflow is extremely comprehensive and suitable for discovery of novel sncRNAs changes, which may lead to the development of innovative diagnostic and therapeutic tools for TNBC.

Biallelic mutations in the 3' exonuclease TOE1 cause pontocerebellar hypoplasia and uncover a role in snRNA processing

Deadenylases are best known for degrading the poly(A) tail during mRNA decay. The deadenylase family has expanded throughout evolution and, in mammals, consists of 12 Mg²⁺-dependent 3'-end RNases with substrate specificity that is mostly unknown. Pontocerebellar hypoplasia type 7 (PCH7) is a unique recessive syndrome characterized by neurodegeneration and ambiguous genitalia. We studied 12 human families with PCH7, uncovering biallelic, loss-of-function mutations in TOE1, which encodes an unconventional deadenylase. toe1-morphant zebrafish displayed midbrain and hindbrain degeneration, modeling PCH-like structural defects in vivo. Surprisingly, we found that TOE1 associated with small nuclear RNAs (snRNAs) incompletely processed spliceosomal. These pre-snRNAs contained 3' genome-encoded tails often followed by post-transcriptionally added adenosines. Human cells with reduced levels of TOE1 accumulated 3'-end-extended pre-snRNAs, and the immunoisolated TOE1 complex was sufficient for 3'-end maturation of snRNAs. Our findings identify the cause of a neurodegenerative syndrome linked to snRNA maturation and uncover a key factor involved in the processing of snRNA 3' ends.

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.

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.

Large Scale Gene Expression Meta-Analysis Reveals Tissue-Specific, Sex-Biased Gene Expression in Humans

The severity and prevalence of many diseases are known to differ between the sexes. Organ specific sex-biased gene expression may underpin these and other sexually dimorphic traits. To further our understanding of sex differences in transcriptional regulation, we performed meta-analyses of sex biased gene expression in multiple human tissues. We analyzed 22 publicly available human gene expression microarray data sets including over 2500 samples from 15 different tissues and 9 different organs. Briefly, by using an inverse-variance method we determined the effect size difference of gene expression between males and females. We found the greatest sex differences in gene expression in the brain, specifically in the anterior cingulate cortex, (1818 genes), followed by the heart (375 genes), kidney (224 genes), colon (218 genes), and thyroid (163 genes). More interestingly, we found different parts of the brain with varying numbers and identity of sex-biased genes, indicating that specific cortical regions may influence sexually dimorphic traits. The majority of sex-biased genes in other tissues such as the bladder, liver, lungs, and pancreas were on the sex chromosomes or involved in sex hormone production. On average in each tissue, 32% of autosomal genes that were expressed in a sex-biased fashion contained androgen or estrogen hormone response elements. Interestingly, across all tissues, we found approximately two-thirds of autosomal genes that were sex-biased were not under direct influence of sex hormones. To our knowledge this is the largest analysis of sex-biased gene expression in human tissues to date. We identified many sex-biased genes that were not under the direct influence of sex chromosome genes or sex hormones. These may provide targets for future development of sex-specific treatments for diseases.

Analysis of the human immunodeficiency virus-1 RNA packageome

All retroviruses package cellular RNAs into virions. Studies of murine leukemia virus (MLV) revealed that the major host cell RNAs encapsidated by this simple retrovirus were LTR retrotransposons and noncoding RNAs (ncRNAs). Several classes of ncRNAs appeared to be packaged by MLV shortly after synthesis, as precursors to tRNAs, small nuclear RNAs, and small nucleolar RNAs were all enriched in virions. To determine the extent to which the human immunodeficiency virus (HIV-1) packages similar RNAs, we used high-throughput sequencing to characterize the RNAs within infectious HIV-1 virions produced in CEM-SS T lymphoblastoid cells. We report that the most abundant cellular RNAs in HIV-1 virions are 7SL RNA and transcripts from numerous divergent and truncated members of the long interspersed element (LINE) and short interspersed element (SINE) families of retrotransposons. We also detected precursors to several tRNAs and small nuclear RNAs as well as transcripts derived from the ribosomal DNA (rDNA) intergenic spacers. We show that packaging of a pre-tRNA requires the nuclear export receptor Exportin 5, indicating that HIV-1 recruits at least some newly made ncRNAs in the cytoplasm. Together, our work identifies the set of RNAs packaged by HIV-1 and reveals that early steps in HIV-1 assembly intersect with host cell ncRNA biogenesis pathways.

The alternative splicing program of differentiated smooth muscle cells involves concerted non-productive splicing of post-transcriptional regulators

Alternative splicing (AS) is a key component of gene expression programs that drive cellular differentiation. Smooth muscle cells (SMCs) are important in the function of a number of physiological systems; however, investigation of SMC AS has been restricted to a handful of events. We profiled transcriptome changes in mouse de-differentiating SMCs and observed changes in hundreds of AS events. Exons included in differentiated cells were characterized by particularly weak splice sites and by upstream binding sites for Polypyrimidine Tract Binding protein (PTBP1). Consistent with this, knockdown experiments showed that that PTBP1 represses many smooth muscle specific exons. We also observed coordinated splicing changes predicted to downregulate the expression of core components of U1 and U2 snRNPs, splicing regulators and other post-transcriptional factors in differentiated cells. The levels of cognate proteins were lower or similar in differentiated compared to undifferentiated cells. However, levels of snRNAs did not follow the expression of splicing proteins, and in the case of U1 snRNP we saw reciprocal changes in the levels of U1 snRNA and U1 snRNP proteins. Our results suggest that the AS program in differentiated SMCs is orchestrated by the combined influence of auxiliary RNA binding proteins, such as PTBP1, along with altered activity and stoichiometry of the core splicing machinery.

Cajal bodies are linked to genome conformation

The mechanisms underlying nuclear body (NB) formation and their contribution to genome function are unknown. Here we examined the non-random positioning of Cajal bodies (CBs), major NBs involved in spliceosomal snRNP assembly and their role in genome organization. CBs are predominantly located at the periphery of chromosome territories at a multi-chromosome interface. Genome-wide chromosome conformation capture analysis (4C-seq) using CB-interacting loci revealed that CB-associated regions are enriched with highly expressed histone genes and U small nuclear or nucleolar RNA (sn/snoRNA) loci that form intra- and inter-chromosomal clusters. In particular, we observed a number of CB-dependent gene-positioning events on chromosome 1. RNAi-mediated disassembly of CBs disrupts the CB-targeting gene clusters and suppresses the expression of U sn/snoRNA and histone genes. This loss of spliceosomal snRNP production results in increased splicing noise, even in CB-distal regions. Therefore, we conclude that CBs contribute to genome organization with global effects on gene expression and RNA splicing fidelity.

Nuclear bodies can nucleate at sites of active transcription and are beneficial for efficient gene expression. Here, the authors show that Cajal bodies, a prominent type of nuclear body, contribute to genome organization with global effects on gene expression and RNA splicing fidelity.

A U1 snRNP-specific assembly pathway reveals the SMN complex as a versatile hub for RNP exchange

Despite equal snRNP stoichiometry in spliceosomes, U1 snRNP (U1) is typically the most abundant vertebrate snRNP. Mechanisms regulating U1 overabundance and snRNP repertoire are unknown. In Sm-core assembly, a key snRNP-biogenesis step mediated by the SMN complex, the snRNA-specific RNA-binding protein (RBP) Gemin5 delivers pre-snRNAs, which join SMN-Gemin2-recruited Sm proteins. We show that the human U1-specific RBP U1-70K can bridge pre-U1 to SMN-Gemin2-Sm, in a Gemin5-independent manner, thus establishing an additional and U1-exclusive Sm core-assembly pathway. U1-70K hijacks SMN-Gemin2-Sm, enhancing Sm-core assembly on U1s and inhibiting that on other snRNAs, thereby promoting U1 overabundance and regulating snRNP repertoire. SMN-Gemin2's ability to facilitate transactions between different RBPs and RNAs explains its multi-RBP valency and the myriad transcriptome perturbations associated with SMN deficiency in neurodegenerative spinal muscular atrophy. We propose that SMN-Gemin2 is a versatile hub for RNP exchange that functions broadly in RNA metabolism.