Genome-wide RNA-binding analysis of the trypanosome U1 snRNP proteins U1C and U1-70K reveals cis/trans-spliceosomal network

Computational analysis of sense-antisense chimeric transcripts reveals their potential regulatory features and the landscape of expression in human cells

Many human genes are transcribed from both strands and produce sense-antisense gene pairs. Sense-antisense (SAS) chimeric transcripts are produced upon the coalescing of exons/introns from both sense and antisense transcripts of the same gene. SAS chimera was first reported in prostate cancer cells. Subsequently, numerous SAS chimeras have been reported in the ChiTaRS-2.1 database. However, the landscape of their expression in human cells and functional aspects are still unknown. We found that longer palindromic sequences are a unique feature of SAS chimeras. Structural analysis indicates that a long hairpin-like structure formed by many consecutive Watson-Crick base pairs appears because of these long palindromic sequences, which possibly play a similar role as double-stranded RNA (dsRNA), interfering with gene expression. RNA-RNA interaction analysis suggested that SAS chimeras could significantly interact with their parental mRNAs, indicating their potential regulatory features. Here, 267 SAS chimeras were mapped in RNA-seq data from 16 healthy human tissues, revealing their expression in normal cells. Evolutionary analysis suggested the positive selection favoring sense-antisense fusions that significantly impacted the evolution of their function and structure. Overall, our study provides detailed insight into the expression landscape of SAS chimeras in human cells and identifies potential regulatory features.

The pan-cancer lncRNA PLANE regulates an alternative splicing program to promote cancer pathogenesis

Genomic amplification of the distal portion of chromosome 3q, which encodes a number of oncogenic proteins, is one of the most frequent chromosomal abnormalities in malignancy. Here we functionally characterise a non-protein product of the 3q region, the long noncoding RNA (lncRNA) PLANE, which is upregulated in diverse cancer types through copy number gain as well as E2F1-mediated transcriptional activation. PLANE forms an RNA-RNA duplex with the nuclear receptor co-repressor 2 (NCOR2) pre-mRNA at intron 45, binds to heterogeneous ribonucleoprotein M (hnRNPM) and facilitates the association of hnRNPM with the intron, thus leading to repression of the alternative splicing (AS) event generating NCOR2-202, a major protein-coding NCOR2 AS variant. This is, at least in part, responsible for PLANE-mediated promotion of cancer cell proliferation and tumorigenicity. These results uncover the function and regulation of PLANE and suggest that PLANE may constitute a therapeutic target in the pan-cancer context.

SNRPC promotes hepatocellular carcinoma cell motility by inducing epithelial-mesenchymal transition

The therapeutic outcome of hepatocellular carcinoma (HCC) remains unsatisfactory because of poor response and acquired drug resistance. To better elucidate the molecular mechanisms of HCC, here we used three Gene Expression Omnibus datasets to identify potential oncogenes, and thereby identified small nuclear ribonucleoprotein polypeptide C (SNRPC). We report that SNRPC is highly up‐regulated in HCC tissues as determined using immunohistochemistry assays of samples from a cohort of 224 patients with HCC, and overexpression of SNRPC was correlated with multiple tumors, advanced stage, and poor outcome. Kaplan–Meier analysis confirmed that patients with high SNRPC expression exhibited shorter survival in four independent HCC cohorts (all P < 0.05). Furthermore, SNRPC mutations are significantly more frequent in HCC tissues than in normal liver tissues and are an early event in the development of HCC. Functional network analysis suggested that SNRPC is linked to the regulation of ribosome, spliceosome, and proteasome signaling. Subsequently, gain‐ and loss‐of‐function assays showed that SNRPC promotes the motility and epithelial–mesenchymal transition of HCC cells in vitro. SNRPC expression was negatively correlated with the infiltration of CD4⁺ T cells, macrophage cells, and neutrophil cells (all P < 0.05), as determined by analyzing the TIMER (Tumor IMmune Estimation Resource) database. In conclusion, our findings suggest that SNRPC has a potential role in epithelial–mesenchymal transition and motility in HCC.

Nuclear mRNA maturation and mRNA export control: from trypanosomes to opisthokonts

The passage of mRNAs through the nuclear pores into the cytoplasm is essential in all eukaryotes. For regulation, mRNA export is tightly connected to the full machinery of nuclear mRNA processing, starting at transcription. Export competence of pre-mRNAs gradually increases by both transient and permanent interactions with multiple RNA processing and export factors. mRNA export is best understood in opisthokonts, with limited knowledge in plants and protozoa. Here, I review and compare nuclear mRNA processing and export between opisthokonts and Trypanosoma brucei. The parasite has many unusual features in nuclear mRNA processing, such as polycistronic transcription and trans-splicing. It lacks several nuclear complexes and nuclear-pore-associated proteins that in opisthokonts play major roles in mRNA export. As a consequence, trypanosome mRNA export control is not tight and export can even start co-transcriptionally. Whether trypanosomes regulate mRNA export at all, or whether leakage of immature mRNA to the cytoplasm is kept to a low level by a fast kinetics of mRNA processing remains to be investigated. mRNA export had to be present in the last common ancestor of eukaryotes. Trypanosomes are evolutionary very distant from opisthokonts and a comparison helps understanding the evolution of mRNA export.

iCLIP analysis of RNA substrates of the archaeal exosome

Background

The archaeal exosome is an exoribonucleolytic multiprotein complex, which degrades single-stranded RNA in 3′ to 5′ direction phosphorolytically. In a reverse reaction, it can add A-rich tails to the 3′-end of RNA. The catalytic center of the exosome is in the aRrp41 subunit of its hexameric core. Its RNA-binding subunits aRrp4 and aDnaG confer poly(A) preference to the complex. The archaeal exosome was intensely characterized in vitro, but still little is known about its interaction with natural substrates in the cell, particularly because analysis of the transcriptome-wide interaction of an exoribonuclease with RNA is challenging.

Results

To determine binding sites of the exosome to RNA on a global scale, we performed individual-nucleotide resolution UV crosslinking and immunoprecipitation (iCLIP) analysis with antibodies directed against aRrp4 and aRrp41 of the chrenarchaeon Sulfolobus solfataricus. A relatively high proportion (17–19%) of the obtained cDNA reads could not be mapped to the genome. Instead, they corresponded to adenine-rich RNA tails, which are post-transcriptionally synthesized by the exosome, and to circular RNAs (circRNAs). We identified novel circRNAs corresponding to 5′ parts of two homologous, transposase-related mRNAs. To detect preferred substrates of the exosome, the iCLIP reads were compared to the transcript abundance using RNA-Seq data. Among the strongly enriched exosome substrates were RNAs antisense to tRNAs, overlapping 3′-UTRs and RNAs containing poly(A) stretches. The majority of the read counts and crosslink sites mapped in mRNAs. Furthermore, unexpected crosslink sites clustering at 5′-ends of RNAs was detected.

Conclusions

In this study, RNA targets of an exoribonuclease were analyzed by iCLIP. The data documents the role of the archaeal exosome as an exoribonuclease and RNA-tailing enzyme interacting with all RNA classes, and underlines its role in mRNA turnover, which is important for adaptation of prokaryotic cells to changing environmental conditions. The clustering of crosslink sites near 5′-ends of genes suggests simultaneous binding of both RNA ends by the S. solfataricus exosome. This may serve to prevent translation of mRNAs dedicated to degradation in 3′-5′ direction.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07200-x.

Regulation of gene expression in trypanosomatids: living with polycistronic transcription

In trypanosomes, RNA polymerase II transcription is polycistronic and individual mRNAs are excised by trans-splicing and polyadenylation. The lack of individual gene transcription control is compensated by control of mRNA processing, translation and degradation. Although the basic mechanisms of mRNA decay and translation are evolutionarily conserved, there are also unique aspects, such as the existence of six cap-binding translation initiation factor homologues, a novel decapping enzyme and an mRNA stabilizing complex that is recruited by RNA-binding proteins. High-throughput analyses have identified nearly a hundred regulatory mRNA-binding proteins, making trypanosomes valuable as a model system to investigate post-transcriptional regulation.

Ranking noncanonical 5' splice site usage by genome-wide RNA-seq analysis and splicing reporter assays

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.

The polyadenylation complex of Trypanosoma brucei: Characterization of the functional poly(A) polymerase

The generation of mature mRNA in the protozoan parasite Trypanosoma brucei requires coupled polyadenylation and trans splicing. In contrast to other eukaryotes, we still know very little on components, mechanisms, and dynamics of the 3' end-processing machinery in trypanosomes. To characterize the catalytic core of the polyadenylation complex in T. brucei, we first identified the poly(A) polymerase [Tb927.7.3780] as the major functional, nuclear-localized enzyme in trypanosomes. In contrast, another poly(A) polymerase, encoded by an intron-containing gene [Tb927.3.3160], localizes mainly in the cytoplasm and appears not to be functional in general 3' end processing of mRNAs. Based on tandem-affinity purification with tagged CPSF160 and mass spectrometry, we identified ten associated components of the trypanosome polyadenylation complex, including homologues to all four CPSF subunits, Fip1, CstF50/64, and Symplekin, as well as two hypothetical proteins. RNAi-mediated knockdown revealed that most of these factors are essential for growth and required for both in vivo polyadenylation and trans splicing, arguing for a general coupling of these two mRNA-processing reactions.

Co-evolution of SNF spliceosomal proteins with their RNA targets in trans-splicing nematodes

Although the mechanism of pre-mRNA splicing has been well characterized, the evolution of spliceosomal proteins is poorly understood. The U1A/U2B″/SNF family (hereafter referred to as the SNF family) of RNA binding spliceosomal proteins participates in both the U1 and U2 small interacting nuclear ribonucleoproteins (snRNPs). The highly constrained nature of this system has inhibited an analysis of co-evolutionary trends between the proteins and their RNA binding targets. Here we report accelerated sequence evolution in the SNF protein family in Phylum Nematoda, which has allowed an analysis of protein:RNA co-evolution. In a comparison of SNF genes from ecdysozoan species, we found a correlation between trans-splicing species (nematodes) and increased phylogenetic branch lengths of the SNF protein family, with respect to their sister clade Arthropoda. In particular, we found that nematodes (~70-80 % of pre-mRNAs are trans-spliced) have experienced higher rates of SNF sequence evolution than arthropods (predominantly cis-spliced) at both the nucleotide and amino acid levels. Interestingly, this increased evolutionary rate correlates with the reliance on trans-splicing by nematodes, which would alter the role of the SNF family of spliceosomal proteins. We mapped amino acid substitutions to functionally important regions of the SNF protein, specifically to sites that are predicted to disrupt protein:RNA and protein:protein interactions. Finally, we investigated SNF's RNA targets: the U1 and U2 snRNAs. Both are more divergent in nematodes than arthropods, suggesting the RNAs have co-evolved with SNF in order to maintain the necessarily high affinity interaction that has been characterized in other species.