Introns are mediators of cell response to starvation.

Introns are ubiquitous features of all eukaryotic cells. Introns need to be removed from nascent messenger RNA through the process of splicing to produce functional proteins. Here we show that the physical presence of introns in the genome promotes cell survival under starvation conditions. A systematic deletion set of all known introns in budding yeast genes indicates that, in most cases, cells with an intron deletion are impaired when nutrients are depleted. This effect of introns on growth is not linked to the expression of the host gene, and was reproduced even when translation of the host mRNA was blocked. Transcriptomic and genetic analyses indicate that introns promote resistance to starvation by enhancing the repression of ribosomal protein genes that are downstream of the nutrient-sensing TORC1 and PKA pathways. Our results reveal functions of introns that may help to explain their evolutionary preservation in genes, and uncover regulatory mechanisms of cell adaptations to starvation.

snoDB 2.0: an enhanced interactive database, specializing in human snoRNAs.

snoDB is an interactive database of human small nucleolar RNAs (snoRNAs) that includes up-to-date information on snoRNA features, genomic location, conservation, host gene, snoRNA-RNA targets and snoRNA abundance and provides links to other resources. In the second edition of this database (snoDB 2.0), we added an entirely new section on ribosomal RNA (rRNA) chemical modifications guided by snoRNAs with easy navigation between the different rRNA versions used in the literature and experimentally measured levels of modification. We also included new layers of information, including snoRNA motifs, secondary structure prediction, snoRNA-protein interactions, copy annotations and low structure bias expression data in a wide panel of tissues and cell lines to bolster functional probing of snoRNA biology. Version 2.0 features updated identifiers, more links to external resources and duplicate entry resolution. As a result, snoDB 2.0, which is freely available at https://bioinfo-scottgroup.med.usherbrooke.ca/snoDB/, represents a one-stop shop for snoRNA features, rRNA modification targets, functional impact and potential regulators.

The snoGloBe interaction predictor reveals a broad spectrum of C/D snoRNA RNA targets.

Box C/D small nucleolar RNAs (snoRNAs) are a conserved class of RNA known for their role in guiding ribosomal RNA 2'-O-ribose methylation. Recently, C/D snoRNAs were also implicated in regulating the expression of non-ribosomal genes through different modes of binding. Large scale RNA-RNA interaction datasets detect many snoRNAs binding messenger RNA, but are limited by specific experimental conditions. To enable a more comprehensive study of C/D snoRNA interactions, we created snoGloBe, a human C/D snoRNA interaction predictor based on a gradient boosting classifier. SnoGloBe considers the target type, position and sequence of the interactions, enabling it to outperform existing predictors. Interestingly, for specific snoRNAs, snoGloBe identifies strong enrichment of interactions near gene expression regulatory elements including splice sites. Abundance and splicing of predicted targets were altered upon the knockdown of their associated snoRNA. Strikingly, the predicted snoRNA interactions often overlap with the binding sites of functionally related RNA binding proteins, reinforcing their role in gene expression regulation. SnoGloBe is also an excellent tool for discovering viral RNA targets, as shown by its capacity to identify snoRNAs targeting the heavily methylated SARS-CoV-2 RNA. Overall, snoGloBe is capable of identifying experimentally validated binding sites and predicting novel sites with shared regulatory function.

RBFOX2 alters splicing outcome in distinct binding modes with multiple protein partners.

RBFOX2 controls the splicing of a large number of transcripts implicated in cell differentiation and development. Parsing RNA-binding protein datasets, we uncover that RBFOX2 can interact with hnRNPC, hnRNPM and SRSF1 to regulate splicing of a broad range of splicing events using different sequence motifs and binding modes. Using immunoprecipitation, specific RBP knockdown, RNA-seq and splice-sensitive PCR, we show that RBFOX2 can target splice sites using three binding configurations: single, multiple or secondary modes. In the single binding mode RBFOX2 is recruited to its target splice sites through a single canonical binding motif, while in the multiple binding mode RBFOX2 binding sites include the adjacent binding of at least one other RNA binding protein partner. Finally, in the secondary binding mode RBFOX2 likely does not bind the RNA directly but is recruited to splice sites lacking its canonical binding motif through the binding of one of its protein partners. These dynamic modes bind distinct sets of transcripts at different positions and distances relative to alternative splice sites explaining the heterogeneity of RBFOX2 targets and splicing outcomes.

SAPFIR: A webserver for the identification of alternative protein features.

BACKGROUND: Alternative splicing can increase the diversity of gene functions by generating multiple isoforms with different sequences and functions. However, the extent to which splicing events have functional consequences remains unclear and predicting the impact of splicing events on protein activity is limited to gene-specific analysis. RESULTS: To accelerate the identification of functionally relevant alternative splicing events we created SAPFIR, a predictor of protein features associated with alternative splicing events. This webserver tool uses InterProScan to predict protein features such as functional domains, motifs and sites in the human and mouse genomes and link them to alternative splicing events. Alternative protein features are displayed as functions of the transcripts and splice sites. SAPFIR could be used to analyze proteins generated from a single gene or a group of genes and can directly identify alternative protein features in large sequence data sets. The accuracy and utility of SAPFIR was validated by its ability to rediscover previously validated alternative protein domains. In addition, our de novo analysis of public datasets using SAPFIR indicated that only a small portion of alternative protein domains was conserved between human and mouse, and that in human, genes involved in nervous system process, regulation of DNA-templated transcription and aging are more likely to produce isoforms missing functional domains due to alternative splicing. CONCLUSION: Overall SAPFIR represents a new tool for the rapid identification of functional alternative splicing events and enables the identification of cellular functions affected by a defined splicing program. SAPFIR is freely available at https://bioinfo-scottgroup.med.usherbrooke.ca/sapfir/ , a website implemented in Python, with all major browsers supported. The source code is available at https://github.com/DelongZHOU/SAPFIR .

Introns: Good Day Junk Is Bad Day Treasure.

Introns are ubiquitous in eukaryotic transcripts. They are often viewed as junk RNA but the huge energetic burden of transcribing, removing, and degrading them suggests a significant evolutionary advantage. Ostensibly, an intron functions within the host pre-mRNA to regulate its splicing, transport, and degradation. However, recent studies have revealed an entirely new class of trans-acting functions where the presence of intronic RNA in the cell impacts the expression of other genes in trans. Here, we review possible new mechanisms of intron functions, with a focus on the role of yeast introns in regulating the cell growth response to starvation.

Structure of a eukaryotic RNase III postcleavage complex reveals a double-ruler mechanism for substrate selection.

Ribonuclease III (RNase III) enzymes are a family of double-stranded RNA (dsRNA)-specific endoribonucleases required for RNA maturation and gene regulation. Prokaryotic RNase III enzymes have been well characterized, but how eukaryotic RNase IIIs work is less clear. Here, we describe the structure of the Saccharomyces cerevisiae RNase III (Rnt1p) postcleavage complex and explain why Rnt1p binds to RNA stems capped with an NGNN tetraloop. The structure shows specific interactions between a structural motif located at the end of the Rnt1p dsRNA-binding domain (dsRBD) and the guanine nucleotide in the second position of the loop. Strikingly, structural and biochemical analyses indicate that the dsRBD and N-terminal domains (NTDs) of Rnt1p function as two rulers that measure the distance between the tetraloop and the cleavage site. These findings provide a framework for understanding eukaryotic RNase IIIs.

Differential expression of duplicated ribosomal protein genes modifies ribosome composition in response to stress.

In Saccharomyces cerevisiae, most ribosomal proteins are synthesized from duplicated genes, increasing the potential for ribosome heterogeneity. However, the contribution of these duplicated genes to ribosome production and the mechanism determining their relative expression remain unclear. Here we demonstrate that in most cases, one of the two gene copies generate the bulk of the active ribosomes under normal growth conditions, while the other copy is favored only under stress. To understand the origin of these differences in paralog expression and their contribution to ribosome heterogeneity we used RNA polymerase II ChIP-Seq, RNA-seq, polyribosome association and peptide-based mass-spectrometry to compare their transcription potential, splicing, mRNA abundance, translation potential, protein abundance and incorporation into ribosomes. In normal conditions a post-transcriptional expression hierarchy of the duplicated ribosomal protein genes is the product of the efficient splicing, high stability and efficient translation of the major paralog mRNA. Exposure of the cell to stress modifies the expression ratio of the paralogs by repressing the expression of the major paralog and thus increasing the number of ribosomes carrying the minor paralog. Together the data indicate that duplicated ribosomal protein genes underlie a modular network permitting the modification of ribosome composition in response to changing growth conditions.

Reovirus µ2 Protein Impairs Translation to Reduce U5 snRNP Protein Levels.

Mammalian orthoreovirus (MRV) is a double-stranded RNA virus from the Reoviridae family that infects a large range of mammals, including humans. Recently, studies have shown that MRV alters cellular alternative splicing (AS) during viral infection. The structural protein µ2 appears to be the main determinant of these AS modifications by decreasing the levels of U5 core components EFTUD2, PRPF8, and SNRNP200 during infection. In the present study, we investigated the mechanism by which µ2 exerts this effect on the U5 components. Our results revealed that µ2 has no impact on steady-state mRNA levels, RNA export, and protein stability of these U5 snRNP proteins. However, polysome profiling and metabolic labeling of newly synthesized proteins revealed that µ2 exerts an inhibitory effect on global translation. Moreover, we showed that µ2 mutants unable to accumulate in the nucleus retain most of the ability to reduce PRPF8 protein levels, indicating that the effect of µ2 on U5 snRNP components mainly occurs in the cytoplasm. Finally, co-expression experiments demonstrated that µ2 suppresses the expression of U5 snRNP proteins in a dose-dependent manner, and that the expression of specific U5 snRNP core components have different sensitivities to µ2's presence. Altogether, these results suggest a novel mechanism by which the µ2 protein reduces the levels of U5 core components through translation inhibition, allowing this viral protein to alter cellular AS during infection.

Simultaneous sequencing of coding and noncoding RNA reveals a human transcriptome dominated by a small number of highly expressed noncoding genes.

Comparing the abundance of one RNA molecule to another is crucial for understanding cellular functions but most sequencing techniques can target only specific subsets of RNA. In this study, we used a new fragmented ribodepleted TGIRT sequencing method that uses a thermostable group II intron reverse transcriptase (TGIRT) to generate a portrait of the human transcriptome depicting the quantitative relationship of all classes of nonribosomal RNA longer than 60 nt. Comparison between different sequencing methods indicated that FRT is more accurate in ranking both mRNA and noncoding RNA than viral reverse transcriptase-based sequencing methods, even those that specifically target these species. Measurements of RNA abundance in different cell lines using this method correlate with biochemical estimates, confirming tRNA as the most abundant nonribosomal RNA biotype. However, the single most abundant transcript is 7SL RNA, a component of the signal recognition particle. Structured noncoding RNAs (sncRNAs) associated with the same biological process are expressed at similar levels, with the exception of RNAs with multiple functions like U1 snRNA. In general, sncRNAs forming RNPs are hundreds to thousands of times more abundant than their mRNA counterparts. Surprisingly, only 50 sncRNA genes produce half of the non-rRNA transcripts detected in two different cell lines. Together the results indicate that the human transcriptome is dominated by a small number of highly expressed sncRNAs specializing in functions related to translation and splicing.

CoCo: RNA-seq read assignment correction for nested genes and multimapped reads.

MOTIVATION: Next-generation sequencing techniques revolutionized the study of RNA expression by permitting whole transcriptome analysis. However, sequencing reads generated from nested and multi-copy genes are often either misassigned or discarded, which greatly reduces both quantification accuracy and gene coverage. RESULTS: Here we present count corrector (CoCo), a read assignment pipeline that takes into account the multitude of overlapping and repetitive genes in the transcriptome of higher eukaryotes. CoCo uses a modified annotation file that highlights nested genes and proportionally distributes multimapped reads between repeated sequences. CoCo salvages over 15% of discarded aligned RNA-seq reads and significantly changes the abundance estimates for both coding and non-coding RNA as validated by PCR and bedgraph comparisons. AVAILABILITY AND IMPLEMENTATION: The CoCo software is an open source package written in Python and available from http://gitlabscottgroup.med.usherbrooke.ca/scott-group/coco. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

The catalytic efficiency of yeast ribonuclease III depends on substrate specific product release rate.

Members of the ribonuclease III (RNase III) family regulate gene expression by triggering the degradation of double stranded RNA (dsRNA). Hundreds of RNase III cleavage targets have been identified and their impact on RNA maturation and stability is now established. However, the mechanism defining substrates' reactivity remains unclear. In this study, we developed a real-time FRET assay for the detection of dsRNA degradation by yeast RNase III (Rnt1p) and characterized the kinetic bottlenecks controlling the reactivity of different substrates. Surprisingly, the results indicate that Rnt1p cleavage reaction is not only limited by the rate of catalysis but can also depend on base-pairing of product termini. Cleavage products terminating with paired nucleotides, like the degradation signals found in coding mRNA sequence, were less reactive and more prone to inhibition than products having unpaired nucleotides found in non-coding RNA substrates. Mutational analysis of U5 snRNA and Mig2 mRNA confirms the pairing of the cleavage site as a major determinant for the difference between cleavage rates of coding and non-coding RNA. Together the data indicate that the base-pairing of Rnt1p substrates encodes reactivity determinants that permit both constitutive processing of non-coding RNA while limiting the rate of mRNA degradation.

Tumor microenvironment-associated modifications of alternative splicing.

Pre-mRNA alternative splicing is modified in cancer, but the origin and specificity of these changes remain unclear. Here, we probed ovarian tumors to identify cancer-associated splicing isoforms and define the mechanism by which splicing is modified in cancer cells. Using high-throughput quantitative PCR, we monitored the expression of splice variants in laser-dissected tissues from ovarian tumors. Surprisingly, changes in alternative splicing were not limited to the tumor tissues but were also found in the tumor microenvironment. Changes in the tumor-associated splicing events were found to be regulated by splicing factors that are differentially expressed in cancer tissues. Overall, ∼20% of the alternative splicing events affected by the down-regulation of the splicing factors QKI and RBFOX2 were altered in the microenvironment of ovarian tumors. Together, our results indicate that the tumor microenvironment undergoes specific changes in alternative splicing orchestrated by a limited number of splicing factors.

Identification of discrete classes of small nucleolar RNA featuring different ends and RNA binding protein dependency.

Small nucleolar RNAs (snoRNAs) are among the first discovered and most extensively studied group of small non-coding RNA. However, most studies focused on a small subset of snoRNAs that guide the modification of ribosomal RNA. In this study, we annotated the expression pattern of all box C/D snoRNAs in normal and cancer cell lines independent of their functions. The results indicate that C/D snoRNAs are expressed as two distinct forms differing in their ends with respect to boxes C and D and in their terminal stem length. Both forms are overexpressed in cancer cell lines but display a conserved end distribution. Surprisingly, the long forms are more dependent than the short forms on the expression of the core snoRNP protein NOP58, thought to be essential for C/D snoRNA production. In contrast, a subset of short forms are dependent on the splicing factor RBFOX2. Analysis of the potential secondary structure of both forms indicates that the k-turn motif required for binding of NOP58 is less stable in short forms which are thus less likely to mature into a canonical snoRNP. Taken together the data suggest that C/D snoRNAs are divided into at least two groups with distinct maturation and functional preferences.

Redirecting splicing with bifunctional oligonucleotides.

Ectopic modulators of alternative splicing are important tools to study the function of splice variants and for correcting mis-splicing events that cause human diseases. Such modulators can be bifunctional oligonucleotides made of an antisense portion that determines target specificity, and a non-hybridizing tail that recruits proteins or RNA/protein complexes that affect splice site selection (TOSS and TOES, respectively, for targeted oligonucleotide silencer of splicing and targeted oligonucleotide enhancer of splicing). The use of TOSS and TOES has been restricted to a handful of targets. To generalize the applicability and demonstrate the robustness of TOSS, we have tested this approach on more than 50 alternative splicing events. Moreover, we have developed an algorithm that can design active TOSS with a success rate of 80%. To produce bifunctional oligonucleotides capable of stimulating splicing, we built on the observation that binding sites for TDP-43 can stimulate splicing and improve U1 snRNP binding when inserted downstream from 5' splice sites. A TOES designed to recruit TDP-43 improved exon 7 inclusion in SMN2. Overall, our study shows that bifunctional oligonucleotides can redirect splicing on a variety of genes, justifying their inclusion in the molecular arsenal that aims to alter the production of splice variants.

Cell cycle-dependent transcription factors control the expression of yeast telomerase RNA.

Telomerase is a specialized ribonucleoprotein that adds repeated DNA sequences to the ends of eukaryotic chromosomes to preserve genome integrity. Some secondary structure features of the telomerase RNA are very well conserved, and it serves as a central scaffold for the binding of associated proteins. The Saccharomyces cerevisiae telomerase RNA, TLC1, is found in very low copy number in the cell and is the limiting component of the known telomerase holoenzyme constituents. The reasons for this low abundance are unclear, but given that the RNA is very stable, transcriptional control mechanisms must be extremely important. Here we define the sequences forming the TLC1 promoter and identify the elements required for its low expression level, including enhancer and repressor elements. Within an enhancer element, we found consensus sites for Mbp1/Swi4 association, and chromatin immunoprecipitation (ChIP) assays confirmed the binding of Mbp1 and Swi4 to these sites of the TLC1 promoter. Furthermore, the enhancer element conferred cell cycle-dependent regulation to a reporter gene, and mutations in the Mbp1/Swi4 binding sites affected the levels of telomerase RNA and telomere length. Finally, ChIP experiments using a TLC1 RNA-binding protein as target showed cell cycle-dependent transcription of the TLC1 gene. These results indicate that the budding yeast TLC1 RNA is transcribed in a cell cycle-dependent fashion late in G1 and may be part of the S phase-regulated group of genes involved in DNA replication.

Regulation of conditional gene expression by coupled transcription repression and RNA degradation.

Gene expression is determined by a combination of transcriptional and post-transcriptional regulatory events that were thought to occur independently. This report demonstrates that the genes associated with the Snf3p-Rgt2p glucose-sensing pathway are regulated by interconnected transcription repression and RNA degradation. Deletion of the dsRNA-specific ribonuclease III Rnt1p increased the expression of Snf3p-Rgt2p-associated transcription factors in vivo and the recombinant enzyme degraded their messenger RNA in vitro. Surprisingly, Rnt1ps effect on gene expression in vivo was both RNA and promoter dependent, thus linking RNA degradation to transcription. Strikingly, deletion of RNT1-induced promoter-specific transcription of the glucose sensing genes even in the absence of RNA cleavage signals. Together, the results presented here support a model in which co-transcriptional RNA degradation increases the efficiency of gene repression, thereby allowing an effective cellular response to the continuous changes in nutrient concentrations.

RNA-dependent regulation of the cell wall stress response.

Stress response requires the precise modulation of gene expression in response to changes in growth conditions. This report demonstrates that selective nuclear mRNA degradation is required for both the cell wall stress response and the regulation of the cell wall integrity checkpoint. More specifically, the deletion of the yeast nuclear dsRNA-specific ribonuclease III (Rnt1p) increased the expression of the mRNAs associated with both the morphogenesis checkpoint and the cell wall integrity pathway, leading to an attenuation of the stress response. The over-expression of selected Rnt1p substrates, including the stress associated morphogenesis protein kinase Hsl1p, in wild-type cells mimicked the effect of RNT1 deletion on cell wall integrity, and their mRNAs were directly cleaved by the recombinant enzyme in vitro. The data supports a model for gene regulation in which nuclear mRNA degradation optimizes the cell response to stress and links it to the cell cycle.

Budding yeast telomerase RNA transcription termination is dictated by the Nrd1/Nab3 non-coding RNA termination pathway.

The RNA component of budding yeast telomerase (Tlc1) occurs in two forms, a non-polyadenylated form found in functional telomerase and a rare polyadenylated version with unknown function. Previous work suggested that the functional Tlc1 polyA- RNA is processed from the polyA+ form, but the mechanisms regulating its transcription termination and 3'-end formation remained unclear. Here we examined transcription termination of Tlc1 RNA in the sequences 3' of the TLC1 gene and relate it to telomere maintenance. Strikingly, disruption of all probable or cryptic polyadenylation signals near the 3'-end blocked the accumulation of the previously reported polyA+ RNA without affecting the level, function or specific 3' nucleotide of the mature polyA- form. A genetic approach analysing TLC1 3'-end sequences revealed that transcription terminates upstream of the polyadenylation sites. Furthermore, the results also demonstrate that the function of this Tlc1 terminator depends on the Nrd1/Nab3 transcription termination pathway. The data thus show that transcription termination of the budding yeast telomerase RNA occurs as that of snRNAs and Tlc1 functions in telomere maintenance are not strictly dependent on a polyadenylated precursor, even if the polyA+ form can serve as intermediate in a redundant termination/maturation pathway.

Differential effects of hnRNP D/AUF1 isoforms on HIV-1 gene expression.

Control of RNA processing plays a major role in HIV-1 gene expression. To explore the role of several hnRNP proteins in this process, we carried out a siRNA screen to examine the effect of depletion of hnRNPs A1, A2, D, H, I and K on HIV-1 gene expression. While loss of hnRNPs H, I or K had little effect, depletion of A1 and A2 increased expression of viral structural proteins. In contrast, reduced hnRNP D expression decreased synthesis of HIV-1 Gag and Env. Loss of hnRNP D induced no changes in viral RNA abundance but reduced the accumulation of HIV-1 unspliced and singly spliced RNAs in the cytoplasm. Subsequent analyses determined that hnRNP D underwent relocalization to the cytoplasm upon HIV-1 infection and was associated with Gag protein. Screening of the four isoforms of hnRNP D determined that, upon overexpression, they had differential effects on HIV-1 Gag expression, p45 and p42 isoforms increased viral Gag synthesis while p40 and p37 suppressed it. The differential effect of hnRNP D isoforms on HIV-1 expression suggests that their relative abundance could contribute to the permissiveness of cell types to replicate the virus, a hypothesis subsequently confirmed by selective depletion of p45 and p42.