Search

Rbm3 (RNA-binding motif protein 3) is a stress responsive gene, which maintains cellular homeostasis and promotes survival upon various harmful cellular stimuli. Rbm3 protein shows conserved structural and molecular similarities to heterogeneous nuclear ribonucleoproteins (hnRNPs), which regulate all steps of the mRNA metabolism. Growing evidence is pointing towards a broader role of Rbm3 in various steps of gene expression. Here, we demonstrate that Rbm3 deficiency is linked to transcriptome-wide pre-mRNA splicing alterations, which can be reversed through Rbm3 co-expression from a cDNA. Using an MS2 tethering assay, we show that Rbm3 regulates splice site selection similar to other hnRNP proteins when recruited between two competing 5   ' splice sites. Furthermore, we show that the N-terminal part of Rbm3 encompassing the RNA recognition motif (RRM), is sufficient to elicit changes in splice site selection. On the basis of these findings, we propose a novel, undescribed function of Rbm3 in RNA splicing that contributes to the preservation of transcriptome integrity.

The underlying molecular mechanism and their general effect on the replication capacity of HIV 1 drug-resistance-associated mutations is often poorly understood. To elucidate the effect of two such mutations located in a region with a high density of spicing regulatory elements on the HIV-1-splicing outcome, bioinformatic predictions were combined with transfection and infection experiments. Results show that the previously described R263K drug-resistance-associated integrase mutation has additionally a severe effect on the ESE2b splicing regulatory element (SRE) in exon 2b, which causes loss of SD2b recognition. This was confirmed by an R263R silent mutation with a similar predicted effect on the exon 2b SRE. In contrast, a V260I mutation and its silent counterpart with a lower effect on ESS2b did not exhibit any differences in the splicing pattern. Since HIV-1 highly relies on a balanced splicing reaction, changes in the splicing outcome can contribute to changes in viral replication and might add to the effect of escape mutations toward antiviral drugs. Thus, a classification of mutations purely addressing proteins is insufficient.

Recently, by combining transcriptomics with functional splicing reporter assays we were able to identify GT > GC > TT as the three highest ranked dinucleotides of human 5' splice sites (5'ss). Here, we have extended our investigations to the proteomic characterization of nuclear proteins that bind to canonical and noncanonical 5'ss. Surprisingly, we found that U1 snRNP binding to functional 5'ss sequences prevented components of the DNA damage response (DDR) from binding to the RNA, suggesting a close link between spliceosome arrangement and genome stability. We demonstrate that all tested noncanonical 5'ss sequences are bona-fide targets of the U2-type spliceosome and are bound by U1 snRNP, including U1-C, in the presence of splicing enhancers. The quantity of precipitated U1-C protein was similar for all noncanonical 5'ss dinucleotides, so that the highly different 5'ss usage was likely due to a later step after early U1 snRNP binding. In addition, we show that an internal GT at positions +5/+6 can be advantageous for splicing at position +1 of noncanonical splice sites. Likewise, and in agreement with previous observations, splicing inactive U1 snRNP binding sites could serve as splicing enhancers, which may also explain the higher abundance of U1 snRNPs compared to other U snRNPs. Finally, we observe that an arginine-serine (RS)-rich domain recruitment to stem loop I of the U1 snRNA is functionally sufficient to promote exon-definition and upstream 3'ss activation.

Jumonji-domain-containing protein 6 (JMJD6) is a Fe(II) and 2-oxogluterate (2OG) dependent oxygenase involved in gene regulation through post-translationally modifying nuclear proteins. It is highly expressed in many cancer types and linked to tumor progression and metastasis. Four alternatively-spliced jmjd6 transcripts were annotated. Here, we focus on the two most abundantly expressed ones, which we call jmjd6-2 and jmjd6-Ex5 . TCGA SpliceSeq data revealed a significant decrease of jmjd6-Ex5 transcripts in patients and postmortem tissue of several tumors. The two protein isoforms are distinguished by their C-terminal sequences, which include a serine-rich region (polyS-domain) in JMJD6-2 that is not present in JMJD6-Ex5. Immunoprecipitation followed by LC-MS/MS for JMJD6-Ex5 shows that different sets of proteins interact with JMJD6-2 and JMJD6-Ex5 with only a few overlaps. In particular, we found TFIIF-associating CTD phosphatase (FCP1), proteins of the survival of motor neurons (SMN) complex, heterogeneous nuclear ribonucleoproteins (hnRNPs) and upstream binding factor (UBF) to interact with JMJD6-Ex5. Like JMJD6-2, both UBF and FCP1 comprise a polyS-domain. The polyS domain of JMJD6-2 might block the interaction with polyS-domains of other proteins. In contrast, JMJD6-2 interacts with many SR-like proteins with arginine/serine-rich (RS)-domains, including several splicing factors. In an HIV-based splicing reporter assay, co-expression of JMJD6-2 inhibited exon inclusion, whereas JMJD6-Ex5 did not have any effect. Furthermore, the silencing of jmjd6 by siRNAs favored jmjd6-Ex5 transcripts, suggesting that JMJD6 controls splicing of its own pre-mRNA. The distinct molecular properties of JMJD6-2 and JMJD6-Ex5 open a lead into the functional implications of the variations of their relative abundance in tumors.

Transcription of the HIV-1 provirus generates a viral pre-mRNA, which is alternatively spliced into more than 50 HIV-1 mRNAs encoding all viral proteins. Regulation of viral alternative splice site usage includes the presence of splicing regulatory elements (SREs) which can dramatically impact RNA expression and HIV-1 replication when mutated. Recently, we were able to show that two viral SREs, G I3 -2 and ESE tat , are important players in the generation of viral vif , vpr and tat mRNAs. Furthermore, we demonstrated that masking these SREs by transfected locked nucleic acid (LNA) mixmers affect the viral splicing pattern and viral particle production. With regard to the development of future therapeutic LNA mixmer-based antiretroviral approaches, we delivered the G I3 -2 and the ESE tat LNA mixmers "nakedly", without the use of transfection reagents (gymnosis) into HIV-1 infected cells. Surprisingly, we observed that gymnotically-delivered LNA mixmers accumulated in the cytoplasm, and seemed to co-localize with GW bodies and induced degradation of mRNAs containing their LNA target sequence. The G I3 -2 and the ESE tat LNA-mediated RNA degradation resulted in abrogation of viral replication in HIV-1 infected Jurkat and PM1 cells as well as in PBMCs.

Most human pathogenic mutations in 5' splice sites affect the canonical GT in positions +1 and +2, leading to noncanonical dinucleotides. On the other hand, noncanonical dinucleotides are observed under physiological conditions in ∼1% of all human 5'ss. It is therefore a challenging task to understand the pathogenic mutation mechanisms underlying the conditions under which noncanonical 5'ss are used. In this work, we systematically examined noncanonical 5' splice site selection, both experimentally using splicing competition reporters and by analyzing a large RNA-seq data set of 54 fibroblast samples from 27 subjects containing a total of 2.4 billion gapped reads covering 269,375 exon junctions. From both approaches, we consistently derived a noncanonical 5'ss usage ranking GC > TT > AT > GA > GG > CT. In our competition splicing reporter assay, noncanonical splicing was strictly dependent on the presence of upstream or downstream splicing regulatory elements (SREs), and changes in SREs could be compensated by variation of U1 snRNA complementarity in the competing 5'ss. In particular, we could confirm splicing at different positions (i.e., -1, +1, +5) of a splice site for all noncanonical dinucleotides "weaker" than GC. In our comprehensive RNA-seq data set analysis, noncanonical 5'ss were preferentially detected in weakly used exon junctions of highly expressed genes. Among high-confidence splice sites, they were 10-fold overrepresented in clusters with a neighboring, more frequently used 5'ss. Conversely, these more frequently used neighbors contained only the dinucleotides GT, GC, and TT, in accordance with the above ranking., (© 2018 Erkelenz et al.; Published by Cold Spring Harbor Laboratory Press.)

Merging data from multiple samples is required to detect low expressed transcripts or splicing events that might be present only in a subset of samples. However, the exact number of required replicates enabling the detection of such rare events often remains a mystery but can be approached through probability theory. Here, we describe a probabilistic model, relating the number of observed events in a batch of samples with observation probabilities. Therein, samples appear as a heterogeneous collection of events, which are observed with some probability. The model is evaluated in a batch of 54 transcriptomes of human dermal fibroblast samples. The majority of putative splice-sites (alignment gap-sites) are detected in (almost) all samples or only sporadically, resulting in an U-shaped pattern for observation probabilities. The probabilistic model systematically underestimates event numbers due to a bias resulting from finite sampling. However, using an additional assumption, the probabilistic model can predict observed event numbers within a <10% deviation from the median. Single samples contain a considerable amount of uniquely observed putative splicing events (mean 7122 in alignments from TopHat alignments and 86,215 in alignments from STAR). We conclude that the probabilistic model provides an adequate description for observation of gap-sites in transcriptome data. Thus, the calculation of required sample sizes can be done by application of a simple binomial model to sporadically observed random events. Due to the large number of uniquely observed putative splice-sites and the known stochastic noise in the splicing machinery, it appears advisable to include observation of rare splicing events into analysis objectives. Therefore, it is beneficial to take scores for the validation of gap-sites into account., Competing Interests: The authors declare no conflict of interest.

Genomic alignments of sequenced cellular messenger RNA contain gapped alignments which are interpreted as consequence of intron removal. The resulting gap-sites, genomic locations of alignment gaps, are landmarks representing potential splice-sites. As alignment algorithms report gap-sites with a considerable false discovery rate, validations are required. We describe two quality scores, gap quality score ( gqs ) and weighted gap information score ( wgis ), developed for validation of putative splicing events: While gqs solely relies on alignment data wgis additionally considers information from the genomic sequence. FASTQ files obtained from 54 human dermal fibroblast samples were aligned against the human genome (GRCh38) using TopHat and STAR aligner. Statistical properties of gap-sites validated by gqs and wgis were evaluated by their sequence similarity to known exon-intron borders. Within the 54 samples, TopHat identifies 1,000,380 and STAR reports 6,487,577 gap-sites. Due to the lack of strand information, however, the percentage of identified GT-AG gap-sites is rather low. While gap-sites from TopHat contain ≈89% GT-AG, gap-sites from STAR only contain ≈42% GT-AG dinucleotide pairs in merged data from 54 fibroblast samples. Validation with gqs yields 156,251 gap-sites from TopHat alignments and 166,294 from STAR alignments. Validation with wgis yields 770,327 gap-sites from TopHat alignments and 1,065,596 from STAR alignments. Both alignment algorithms, TopHat and STAR, report gap-sites with considerable false discovery rate, which can drastically be reduced by validation with gqs and wgis ., Competing Interests: The authors declare no conflict of interest.

The open source environment R isf the most widely used software to statistically explore biological data sets including sequence alignments. BAM is the de facto standard file format for sequence alignment. With rbamtools, we provide now a full spectrum of accessibility to BAM for R users such as reading, writing, extraction of subsets and plotting of alignment depth where the script syntax closely follows the SAM/BAM format. Additionally, rbamtools enables fast accumulative tabulation of splicing events over multiple BAM files., (© The Author 2015. Published by Oxford University Press. All rights reserved. For Permissions, please e-mail: journals.permissions@oup.com.)

Background: The viral regulatory protein Tat is essential for establishing a productive transcription from the 5'-LTR promoter during the early phase of viral gene expression. Formation of the Tat-encoding mRNAs requires splicing at the viral 3'ss A3, which has previously been shown to be both negatively and positively regulated by the downstream splicing regulatory elements (SREs) ESS2p and ESE2/ESS2. However, using the novel RESCUE-type computational HEXplorer algorithm, we were recently able to identify another splicing enhancer (ESE(5807-5838), henceforth referred to as ESE tat ) located between ESS2p and ESE2/ESS2. Here we show that ESE tat has a great impact on viral tat-mRNA splicing and that it is fundamental for regulated 3'ss A3 usage., Results: Mutational inactivation or locked nucleic acid (LNA)-directed masking of the ESE tat sequence in the context of a replication-competent virus was associated with a failure (i) to activate viral 3'ss A3 and (ii) to accumulate Tat-encoding mRNA species. Consequently, due to insufficient amounts of Tat protein efficient viral replication was drastically impaired. RNA in vitro binding assays revealed SRSF2 and SRSF6 as candidate splicing factors acting through ESE tat and ESE2 for 3'ss A3 activation. This notion was supported by coexpression experiments, in which wild-type, but not ESE tat -negative provirus responded to higher levels of SRSF2 and SRSF6 proteins with higher levels of tat-mRNA splicing. Remarkably, we could also find that SRSF6 overexpression established an antiviral state within provirus-transfected cells, efficiently blocking virus particle production. For the anti-HIV-1 activity the arginine-serine (RS)-rich domain of the splicing factor was dispensable., Conclusions: Based on our results, we propose that splicing at 3'ss A3 is dependent on binding of the enhancing SR proteins SRSF2 and SRSF6 to the ESE tat and ESE2 sequence. Mutational inactivation or interference specifically with ESE tat activity by LNA-directed masking seem to account for an early stage defect in viral gene expression, probably by cutting off the supply line of Tat that HIV needs to efficiently transcribe its genome.

HIV-1 mediates pro-survival signals and prevents apoptosis via the phosphatidylinositol-3-kinase (PI3K) pathway. This pathway, however, also affects phosphorylation of serine-arginine (SR) proteins, a family of splicing regulatory factors balancing splice site selection. We now show that pharmacologic inhibition of PI3K signalling alters the HIV-1 splicing pattern of both minigene- and provirus-derived mRNAs. This indicates that HIV-1 might also promote PI3K signalling to balance processing of its transcripts by regulating phosphorylation of splicing regulatory proteins., (© 2014 The Authors.)

Within target T lymphocytes, human immunodeficiency virus type I (HIV-1) encounters the retroviral restriction factor APOBEC3G (apolipoprotein B mRNA-editing enzyme, catalytic polypeptide-like 3G; A3G), which is counteracted by the HIV-1 accessory protein Vif. Vif is encoded by intron-containing viral RNAs that are generated by splicing at 3' splice site (3'ss) A1 but lack splicing at 5'ss D2, which results in the retention of a large downstream intron. Hence, the extents of activation of 3'ss A1 and repression of D2, respectively, determine the levels of vif mRNA and thus the ability to evade A3G-mediated antiviral effects. The use of 3'ss A1 can be enhanced or repressed by splicing regulatory elements that control the recognition of downstream 5'ss D2. Here we show that an intronic G run (G(I2)-1) represses the use of a second 5'ss, termed D2b, that is embedded within intron 2 and, as determined by RNA deep-sequencing analysis, is normally inefficiently used. Mutations of G(I2)-1 and activation of D2b led to the generation of transcripts coding for Gp41 and Rev protein isoforms but primarily led to considerable upregulation of vif mRNA expression. We further demonstrate, however, that higher levels of Vif protein are actually detrimental to viral replication in A3G-expressing T cell lines but not in A3G-deficient cells. These observations suggest that an appropriate ratio of Vif-to-A3G protein levels is required for optimal virus replication and that part of Vif level regulation is effected by the novel G run identified here.

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

Usage of the HIV-1 major 5' splice site D1 is a prerequisite for generation of all spliced viral mRNAs encoding essential regulatory and structural proteins. We set out to determine whether flanking sequences ensure D1-activation. We found that an exonic splicing enhancer function is exerted by the region upstream of D1, which is crucially required for its activation. Additionally, we identified an intronic splicing regulatory element within the p17-instability element of the Gag-ORF enhancing D1-activation. Furthermore, our experimental data demonstrated that sequence motifs displaying high similarity to consensus binding sites for SR protein SC35 (SRSF2) overlapping with D1 fine-tune its activation. Our results reveal that D1-activation is safe-guarded by the interplay of upstream and downstream located splicing enhancer elements ensuring usage of D1 even if its strength is decreased upon mutation. The identification of sequence elements activating D1-usage sheds further light on the balanced expression of alternatively spliced HIV-1 mRNAs., (Copyright © 2012 Elsevier Inc. All rights reserved.)

Splicing regulatory elements (SREs) are sequences bound by proteins that influence splicing of nearby splice sites. Constitutively spliced introns have evolved to utilize many different splicing factors. The evolutionary processes that influenced which splicing factors are used for splicing of individual introns are generally unclear. We demonstrate that in the lineage that gave rise to mammals, many introns lost U-rich sequences and gained G-rich sequences, both of which resemble known SREs. The apparent conversion of U-rich to G-rich SREs suggests that the associated splicing factors are functionally equivalent. In support of this we demonstrated that U-rich and G-rich SREs are both capable of promoting splicing of an SRE-dependent splicing reporter. Furthermore, we demonstrate, using the heterologous MS2 tethering system (bacterial MS2 coat fusion-protein and its RNA stem-loop binding site), that both the U-rich SRE-binding protein (TIA1) and the G-rich SRE-binding protein (HNRNPF) can promote splicing of the same intron. We also observed that gain of G-rich SREs is significantly associated with G/C-rich genomic isochores, suggesting that gain or loss of SREs was driven by the same processes that ultimately resulted in the formation of mammalian genomic isochores. We propose the following model for the gain and loss of mammalian SREs. Ancestral U-rich SREs located in genomic regions that were experiencing high rates of A/T to G/C conversion would have suffered frequent deleterious mutations. However, this same process resulted in increased formation of functionally equivalent G-rich SREs, and acquisition of new G-rich SREs decreased purifying selection on the U-rich SREs, which were then free to decay.

Background: The HIV-1 Rev protein is a key component in the early to late switch in HIV-1 splicing from early intronless (e.g. tat, rev) to late intron-containing Rev-dependent (e.g. gag, vif, env) transcripts. Previous results suggested that cis-acting sequences and inefficient 5' and 3' splice sites are a prerequisite for Rev function. However, we and other groups have shown that two of the HIV-1 5' splice sites, D1 and D4, are efficiently used in vitro and in vivo. Here, we focus on the efficiency of the HIV-1 3' splice sites taking into consideration to what extent their intrinsic efficiencies are modulated by their downstream cis-acting exonic sequences. Furthermore, we delineate their role in RNA stabilization and Rev function., Results: In the presence of an efficient upstream 5' splice site the integrity of the 3' splice site is not essential for Rev function whereas an efficient 3' splice site impairs Rev function. The detrimental effect of a strong 3' splice site on the amount of Rev-dependent intron-containing HIV-1 glycoprotein coding (env) mRNA is not compensatable by weakening the strength of the upstream 5' splice site. Swapping the HIV-1 3' splice sites in an RRE-containing minigene, we found a 3' splice site usage which was variably dependent on the presence of the usual downstream exonic sequence. The most evident activation of 3' splice site usage by its usual downstream exonic sequence was observed for 3' splice site A1 which was turned from an intrinsic very weak 3' splice site into the most active 3' splice site, even abolishing Rev activity. Performing pull-down experiments with nuclear extracts of HeLa cells we identified a novel ASF/SF2-dependent exonic splicing enhancer (ESE) within HIV-1 exon 2 consisting of a heptameric sequence motif occurring twice (M1 and M2) within this short non-coding leader exon. Single point mutation of M1 within an infectious molecular clone is detrimental for HIV-1 exon 2 recognition without affecting Rev-dependent vif expression., Conclusion: Under the conditions of our assay, the rate limiting step of retroviral splicing, competing with Rev function, seems to be exclusively determined by the functional strength of the 3' splice site. The bipartite ASF/SF2-dependent ESE within HIV-1 exon 2 supports cross-talk between splice site pairs across exon 2 (exon definition) which is incompatible with processing of the intron-containing vif mRNA. We propose that Rev mediates a switch from exon to intron definition necessary for the expression of all intron-containing mRNAs.

The posttranscriptional regulatory element (PRE) is considered to enhance hepatitis B virus (HBV) gene expression by facilitating the nuclear export of intronless viral subgenomic RNAs. Its role in the RNA metabolism of the viral pregenomic RNA (pgRNA) is currently unknown. We identified a positively cis-acting splicing regulatory element (SRE-1) and present two lines of evidence for its functionality. Firstly, in a heterologous context SRE-1 functionally substitutes for a retroviral bidirectional exonic splicing enhancer (ESE). As expected, SRE-1 is a splicing enhancer also in its natural viral sequence context, since deletion of SRE-1 reduces splicing of pgRNA in cell culture experiments. Secondly, we show that stimulation of HBV RNA splicing by the splicing factor PSF was repressed by the PRE. Analysis of a variety of PSF mutants indicated that RNA-binding and protein-protein interaction were required to enhance splicing. In addition, we show that the PRE contributed to pgRNA stability, but has little influence on its nuclear export. Herein, we report for the first time that the PRE harbors splicing stimulating and inhibiting regulatory elements controlling processing of the viral pregenome. We discuss a model in which the regulation of pgRNA splicing depends on cellular factors interacting with the PRE.

Spliceosome formation is initiated by the recognition of the 5' splice site through formation of an RNA duplex between the 5' splice site and U1 snRNA. We have previously shown that RNA duplex formation between U1 snRNA and the 5' splice site can protect pre-mRNAs from degradation prior to splicing. This initial RNA duplex must be disrupted to expose the 5' splice site sequence for base pairing with U6 snRNA and to form the active spliceosome. Here, we investigated whether hyperstabilization of the U1 snRNA/5' splice site duplex interferes with splicing efficiency in human cell lines or nuclear extracts. Unlike observations in Saccharomyces cerevisiae, we demonstrate that an extended U1 snRNA/5' splice site interaction does not decrease splicing efficiency, but rather increases 5' splice site recognition and exon inclusion. However, low complementarity of the 5' splice site to U1 snRNA significantly increases exon skipping and RNA degradation. Although the splicing mechanisms are conserved between human and S.cerevisiae, these results demonstrate that distinct differences exist in the activation of the spliceosome.

The integrated human immunodeficiency virus type 1 (HIV-1) genome is transcribed in a single pre-mRNA that is alternatively spliced into more than 40 mRNAs. We characterized a novel bidirectional exonic splicing enhancer (ESE) that regulates the expression of the HIV-1 env, vpu, rev, and nef mRNAs. The ESE is localized downstream of the vpu-, env-, and nef-specific 3' splice site no. 5. SF2/ASF and SRp40 activate the ESE and are required for efficient 3' splice site usage and binding of the U1 snRNP to the downstream 5' splice site no. 4. U1 snRNP binding to the 5' splice site no. 4 is required for splicing of the rev and nef mRNAs and to increase expression of the partially spliced env mRNA. Finally, our results indicate that this ESE is necessary for the recruitment of the U1 snRNP to the 5' splice site no. 4, even when the 5' splice site and the U1 snRNA have been mutated to obtain a perfect complementary match. The ESE characterized here is highly conserved in most viral subtypes.

HIV-1 env expression from certain subgenomic vectors requires the viral regulatory protein Rev, its target sequence RRE, and a 5' splice site upstream of the env open reading frame. To determine the role of this splice site in the 5'-splice-site-dependent Rev-mediated env gene expression, we have subjected the HIV-1 5' splice site, SD4, to a mutational analysis and have analyzed the effect of those mutations on env expression. The results demonstrate that the overall strength of hydrogen bonding between the 5' splice site, SD4, and the free 5' end of the U1 snRNA correlates with env expression efficiency, as long as env expression is suboptimal, and that a continuous stretch of 14 hydrogen bonds can lead to full env expression, as a result of stabilizing the pre-mRNA. The U1 snRNA-mediated stabilization is independent of functional splicing, as a mismatch in position +1 of the 5' splice site that led to loss of detectable amounts of spliced transcripts did not preclude stabilization and expression of the unspliced env mRNA, provided that Rev enables its nuclear export. The nucleotides capable of participating in U1 snRNA:pre-mRNA interaction include positions -3 to +8 of the 5' splice site and all 11 nt constituting the single-stranded 5' end of U1 snRNA. Moreover, env gene expression is significantly decreased upon the introduction of point mutations in several upstream GAR nucleotide motifs, which are mediating SF2/ASF responsiveness in an in vitro splicing assay. This suggests that the GAR sequences may play a role in stabilizing the pre-mRNA by sequestering U1 snRNP to SD4.