Crosslinking-immunoprecipitation (iCLIP) analysis reveals global regulatory roles of hnRNP L

Ribosome profiling reveals that post-transcriptional control of Nalf1 by heterogeneous nuclear ribonucleoprotein L is required for paclitaxel-induced neuropathic pain

ABSTRACT

Sensory neurons are integral to the genesis and maintenance of neuropathic pain. The molecular mechanisms that mediate long-lived changes in their excitability are unclear. Here, we leverage functional genomics approaches to survey changes in RNA abundance and translation in dorsal root ganglion neurons from a mouse model of paclitaxel-induced neuropathic pain. We focus specifically on females as paclitaxel is a first-line therapy for breast cancer. The sequencing data indicate that substantially more changes occur at the level of translation (n = 404) than transcription and decay (n = 109). We discovered that a core subunit of the sodium leak channel (NALCN) channel, auxiliary factor 1 (NALF1), is preferentially translated in response to paclitaxel. This effect is mediated by the RNA-binding protein heterogeneous nuclear ribonucleoprotein L (HNRNP L). Heterogeneous nuclear ribonucleoprotein L binds a 14 base CA-rich element (CARE) in the Nalf1 3' untranslated region (3'UTR). Genetic elimination of either HNRNP L, the Nalf1 CARE motif, or the pore-forming subunit of the nonselective NALCN diminishes pain amplification in vivo. Collectively, these results illustrate that an element situated in a 3'UTR is required for neuropathic pain in female mice.

The glycine-rich domain of GRP7 plays a crucial role in binding long RNAs and facilitating phase separation

Microscale thermophoresis (MST) is a well-established method to quantify protein-RNA interactions. In this study, we employed MST to analyze the RNA binding properties of glycine-rich RNA binding protein 7 (GRP7), which is known to have multiple biological functions related to its ability to bind different types of RNA. However, the exact mechanism of GRP7's RNA binding is not fully understood. While the RNA-recognition motif of GRP7 is known to be involved in RNA binding, the glycine-rich region (known as arginine-glycine-glycine-domain or RGG-domain) also influences this interaction. To investigate to which extend the RGG-domain of GRP7 is involved in RNA binding, mutation studies on putative RNA interacting or modulating sites were performed. In addition to MST experiments, we examined liquid-liquid phase separation of GRP7 and its mutants, both with and without RNA. Furthermore, we systemically investigated factors that might affect RNA binding selectivity of GRP7 by testing RNAs of different sizes, structures, and modifications. Consequently, our study revealed that GRP7 exhibits a high affinity for a variety of RNAs, indicating a lack of pronounced selectivity. Moreover, we established that the RGG-domain plays a crucial role in binding longer RNAs and promoting phase separation.

TDP-43 and other hnRNPs regulate cryptic exon inclusion of a key ALS/FTD risk gene, UNC13A

A major function of TAR DNA-binding protein-43 (TDP-43) is to repress the inclusion of cryptic exons during RNA splicing. One of these cryptic exons is in UNC13A, a genetic risk factor for amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). The accumulation of cryptic UNC13A in disease is heightened by the presence of a risk haplotype located within the cryptic exon itself. Here, we revealed that TDP-43 extreme N-terminus is important to repress UNC13A cryptic exon inclusion. Further, we found hnRNP L, hnRNP A1, and hnRNP A2B1 bind UNC13A RNA and repress cryptic exon inclusion, independently of TDP-43. Finally, higher levels of hnRNP L protein associate with lower burden of UNC13A cryptic RNA in ALS/FTD brains. Our findings suggest that while TDP-43 is the main repressor of UNC13A cryptic exon inclusion, other hnRNPs contribute to its regulation and may potentially function as disease modifiers.

CHD8 suppression impacts on histone H3 lysine 36 trimethylation and alters RNA alternative splicing

Disruptive mutations in the chromodomain helicase DNA-binding protein 8 gene (CHD8) have been recurrently associated with autism spectrum disorders (ASDs). Here we investigated how chromatin reacts to CHD8 suppression by analyzing a panel of histone modifications in induced pluripotent stem cell-derived neural progenitors. CHD8 suppression led to significant reduction (47.82%) in histone H3K36me3 peaks at gene bodies, particularly impacting on transcriptional elongation chromatin states. H3K36me3 reduction specifically affects highly expressed, CHD8-bound genes and correlates with altered alternative splicing patterns of 462 genes implicated in 'regulation of RNA splicing' and 'mRNA catabolic process'. Mass spectrometry analysis uncovered a novel interaction between CHD8 and the splicing regulator heterogeneous nuclear ribonucleoprotein L (hnRNPL), providing the first mechanistic insights to explain the CHD8 suppression-derived splicing phenotype, partly implicating SETD2, a H3K36me3 methyltransferase. In summary, our results point toward broad molecular consequences of CHD8 suppression, entailing altered histone deposition/maintenance and RNA processing regulation as important regulatory processes in ASD.

Positional motif analysis reveals the extent of specificity of protein-RNA interactions observed by CLIP

BACKGROUND

Crosslinking and immunoprecipitation (CLIP) is a method used to identify in vivo RNA-protein binding sites on a transcriptome-wide scale. With the increasing amounts of available data for RNA-binding proteins (RBPs), it is important to understand to what degree the enriched motifs specify the RNA-binding profiles of RBPs in cells.

RESULTS

We develop positionally enriched k-mer analysis (PEKA), a computational tool for efficient analysis of enriched motifs from individual CLIP datasets, which minimizes the impact of technical and regional genomic biases by internal data normalization. We cross-validate PEKA with mCross and show that the use of input control for background correction is not required to yield high specificity of enriched motifs. We identify motif classes with common enrichment patterns across eCLIP datasets and across RNA regions, while also observing variations in the specificity and the extent of motif enrichment across eCLIP datasets, between variant CLIP protocols, and between CLIP and in vitro binding data. Thereby, we gain insights into the contributions of technical and regional genomic biases to the enriched motifs, and find how motif enrichment features relate to the domain composition and low-complexity regions of the studied proteins.

CONCLUSIONS

Our study provides insights into the overall contributions of regional binding preferences, protein domains, and low-complexity regions to the specificity of protein-RNA interactions, and shows the value of cross-motif and cross-RBP comparison for data interpretation. Our results are presented for exploratory analysis via an online platform in an RBP-centric and motif-centric manner ( https://imaps.goodwright.com/apps/peka/ ).

Circular RNAs: Characterization, cellular roles, and applications

Most circular RNAs are produced from the back-splicing of exons of precursor mRNAs. Recent technological advances have in part overcome problems with their circular conformation and sequence overlap with linear cognate mRNAs, allowing a better understanding of their cellular roles. Depending on their localization and specific interactions with DNA, RNA, and proteins, circular RNAs can modulate transcription and splicing, regulate stability and translation of cytoplasmic mRNAs, interfere with signaling pathways, and serve as templates for translation in different biological and pathophysiological contexts. Emerging applications of RNA circles to interfere with cellular processes, modulate immune responses, and direct translation into proteins shed new light on biomedical research. In this review, we discuss approaches used in circular RNA studies and the current understanding of their regulatory roles and potential applications.

The SWI/SNF subunit BRG1 affects alternative splicing by changing RNA binding factor interactions with nascent RNA

BRG1 and BRM are ATPase core subunits of the human SWI/SNF chromatin remodelling complexes mainly associated with transcriptional initiation. They also have a role in alternative splicing, which has been shown for BRM-containing SWI/SNF complexes at a few genes. Here, we have identified a subset of genes which harbour alternative exons that are affected by SWI/SNF ATPases by expressing the ATPases BRG1 and BRM in C33A cells, a BRG1- and BRM-deficient cell line, and analysed the effect on splicing by RNA sequencing. BRG1- and BRM-affected sub-sets of genes favouring both exon inclusion and exon skipping, with only a minor overlap between the ATPase. Some of the changes in alternative splicing induced by BRG1 and BRM expression did not require the ATPase activity. The BRG1-ATPase independent included exons displayed an exon signature of a high GC content. By investigating three genes with exons affected by the BRG-ATPase-deficient variant, we show that these exons accumulated phosphorylated RNA pol II CTD, both serine 2 and serine 5 phosphorylation, without an enrichment of the RNA polymerase II. The ATPases were recruited to the alternative exons, together with both core and signature subunits of SWI/SNF complexes, and promoted the binding of RNA binding factors to chromatin and RNA at the alternative exons. The interaction with the nascent RNP, however, did not reflect the association to chromatin. The hnRNPL, hnRNPU and SAM68 proteins associated with chromatin in cells expressing BRG1 and BRM wild type, but the binding of hnRNPU to the nascent RNP was excluded. This suggests that SWI/SNF can regulate alternative splicing by interacting with splicing-RNA binding factor and influence their binding to the nascent pre-mRNA particle.

Expanding uncapped translation and emerging function of circular RNA in carcinomas and noncarcinomas

Circular RNAs (circRNAs) are classified as noncoding RNAs because they are devoid of a 5' end cap and a 3' end poly (A) tail necessary for cap-dependent translation. However, increasing numbers of translated circRNAs identified through high-throughput RNA sequencing overlapping with polysome profiling indicate that this rule is being broken. CircRNAs can be translated in cap-independent mechanism, including IRES (internal ribosome entry site)-initiated pattern, MIRES (m6A internal ribosome entry site) -initiated patterns, and rolling translation mechanism (RCA). CircRNA-encoded proteins harbour diverse functions similar to or different from host proteins. In addition, they are linked to the modulation of human disease including carcinomas and noncarcinomas. CircRNA-related translatomics and proteomics have attracted increasing attention. This review discusses the progress and exclusive characteristics of circRNA translation and highlights the latest mechanisms and regulation of circRNA translatomics. Furthermore, we summarize the extensive functions and mechanisms of circRNA-derived proteins in human diseases, which contribute to a better understanding of intricate noncanonical circRNA translatomics and proteomics and their therapeutic potential in human diseases.

The oncogenic kinase NEK2 regulates an RBFOX2-dependent pro-mesenchymal splicing program in triple-negative breast cancer cells

BACKGROUND

Triple-negative breast cancer (TNBC) is the most heterogeneous and malignant subtype of breast cancer (BC). TNBC is defined by the absence of expression of estrogen, progesterone and HER2 receptors and lacks efficacious targeted therapies. NEK2 is an oncogenic kinase that is significantly upregulated in TNBC, thereby representing a promising therapeutic target. NEK2 localizes in the nucleus and promotes oncogenic splice variants in different cancer cells. Notably, alternative splicing (AS) dysregulation has recently emerged as a featuring trait of TNBC that contributes to its aggressive phenotype.

METHODS

To investigate whether NEK2 modulates TNBC transcriptome we performed RNA-sequencing analyses in a representative TNBC cell line (MDA-MB-231) and results were validated in multiple TNBC cell lines. Bioinformatics and functional analyses were carried out to elucidate the mechanism of splicing regulation by NEK2. Data from The Cancer Genome Atlas were mined to evaluate the potential of NEK2-sensitive exons as markers to identify the TNBC subtype and to assess their prognostic value.

RESULTS

Transcriptome analysis revealed a widespread impact of NEK2 on the transcriptome of TNBC cells, with 1830 AS events that are susceptible to its expression. NEK2 regulates the inclusion of cassette exons in splice variants that discriminate TNBC from other BC and that correlate with poor prognosis, suggesting that this kinase contributes to the TNBC-specific splicing program. NEK2 elicits its effects by modulating the expression of the splicing factor RBFOX2, a well-known regulator of epithelial to mesenchymal transition (EMT). Accordingly, NEK2 splicing-regulated genes are enriched in functional terms related to cell adhesion and contractile cytoskeleton and NEK2 depletion in mesenchymal TNBC cells induces phenotypic and molecular traits typical of epithelial cells. Remarkably, depletion of select NEK2-sensitive splice-variants that are prognostic in TNBC patients is sufficient to interfere with TNBC cell morphology and motility, suggesting that NEK2 orchestrates a pro-mesenchymal splicing program that modulates migratory and invasive properties of TNBC cells.

CONCLUSIONS

Our study uncovers an extensive splicing program modulated by NEK2 involving splice variants that confer an invasive phenotype to TNBCs and that might represent, together with NEK2 itself, valuable therapeutic targets for this disease.

Improved Predicting of The Sequence Specificities of RNA Binding Proteins by Deep Learning

RNA-binding proteins (RBPs) have a significant role in various regulatory tasks. However, the mechanism by which RBPs identify the subsequence target RNAs is still not clear. In recent years, several machine and deep learning-based computational models have been proposed for understanding the binding preferences of RBPs. These methods required integrating multiple features with raw RNA sequences such as secondary structure and their performances can be further improved. In this paper, we propose an efficient and simple convolution neural network, RBPCNN, that relies on the combination of the raw RNA sequence and evolutionary information. We show that conservation scores (evolutionary information) for the RNA sequences can significantly improve the overall performance of the proposed predictor. In addition, the automatic extraction of the binding sequence motifs can enhance our understanding of the binding specificities of RBPs. The experimental results show that RBPCNN outperforms significantly the current state-of-the-art methods. More specifically, the average area under the receiver operator curve was improved by 2.67 percent and the mean average precision was improved by 8.03 percent. The datasets and results can be downloaded from https://home.jbnu.ac.kr/NSCL/RBPCNN.htm.

HP1γ binding pre-mRNA intronic repeats modulates RNA splicing decisions

HP1 proteins are best known as markers of heterochromatin and gene silencing. Yet, they are also RNA-binding proteins and the HP1γ/CBX3 family member is present on transcribed genes together with RNA polymerase II, where it regulates co-transcriptional processes such as alternative splicing. To gain insight in the role of the RNA-binding activity of HP1γ in transcriptionally active chromatin, we have captured and analysed RNAs associated with this protein. We find that HP1γ is specifically targeted to hexameric RNA motifs and coincidentally transposable elements of the SINE family. As these elements are abundant in introns, while essentially absent from exons, the HP1γ RNA association tethers unspliced pre-mRNA to chromatin via the intronic regions and limits the usage of intronic cryptic splice sites. Thus, our data unveil novel determinants in the relationship between chromatin and co-transcriptional splicing.

HNRNPL Circularizes ARHGAP35 to Produce an Oncogenic Protein

Circular RNAs (circRNAs) are an intriguing class of widely prevalent endogenous RNAs, the vast majority of which have not been characterized functionally. Here, we identified a novel oncogenic circRNA originating from the back-splicing of Exon2 and Exon3 of a tumor suppressor gene, ARHGAP35 (also known as P190-A), termed as circARHGAP35. have observe that circARHGAP35 and linear ARHGAP35 have antithetical expression and functions. Interestingly, circARHGAP35 contains a 3867 nt long ORF with an m6A-modified start codon and encodes a truncated protein comprising four FF domains and lacking the Rho GAP domain. Mechanistically, circARHGAP35 protein promotes cancer cell progression by interacting with TFII-I protein in the nucleus. The RNA binding protein, HNRNPL, facilitates the formation of circARHGAP35. Clinically, circARHGAP35 is associated with poor survival in cancer patients. Our findings characterize an oncogenic circRNA and demonstrate a novel mechanism of oncogene activation in cancer by circRNA through the production of a truncated protein.

The MTNR1A mRNA is stabilized by the cytoplasmic hnRNPL in renal tubular cells

The downregulation of melatonin receptor 1A (MTNR1A) is associated with a range of pathological conditions, including membranous nephropathy. Knowledge of the mechanism underlying MTNR1A expression has been limited to the transcriptional regulation level. Here, RNA interference screening in human kidney cells revealed that heterogeneous nuclear ribonucleoprotein L (hnRNPL) upregulated MTNR1A RNA post-transcriptionally. hnRNPL knockdown or overexpression led to increased or decreased levels of cyclic adenosine monophosphate-responsive element-binding protein phosphorylation, respectively. Molecular studies showed that cytoplasmic hnRNPL exerts a stabilizing effect on the MTNR1A transcript through CA-repeat elements in its coding region. Further studies revealed that the interaction between hnRNPL and MTNR1A serves to protect MNTR1A RNA degradation by the exosome component 10 protein. MTNR1A, but not hnRNPL, displays a diurnal rhythm in mouse kidneys. Enhanced levels of MTNR1A recorded at midnight correlated with robust binding activity between cytoplasmic hnRNPL and the MTNR1A transcript. Both hnRNPL and MTNR1A were decreased in the cytoplasm of tubular epithelial cells from experimental membranous nephropathy kidneys, supporting their clinical relevance. Collectively, our data identified cytoplasmic hnRNPL as a novel player in the upregulation of MTNR1A expression in renal tubular epithelial cells, and as a potential therapeutic target.

Design and application of circular RNAs with protein-sponge function

Circular RNAs (circRNAs) are a class of noncoding RNAs, generated from pre-mRNAs by circular splicing of exons and functionally largely uncharacterized. Here we report on the design, expression, and characterization of artificial circRNAs that act as protein sponges, specifically binding and functionally inactivating hnRNP (heterogeneous nuclear ribonucleoprotein) L. HnRNP L regulates alternative splicing, depending on short CA-rich RNA elements. We demonstrate that designer hnRNP L-sponge circRNAs with CA-repeat or CA-rich sequence clusters can efficiently and specifically modulate splicing-regulatory networks in mammalian cells, including alternative splicing patterns and the cellular distribution of a splicing factor. This new strategy can in principle be applied to any RNA-binding protein, opening up new therapeutic strategies in molecular medicine.

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.

The role of hnRNPs in frontotemporal dementia and amyotrophic lateral sclerosis

Dysregulated RNA metabolism is emerging as a crucially important mechanism underpinning the pathogenesis of frontotemporal dementia (FTD) and the clinically, genetically and pathologically overlapping disorder of amyotrophic lateral sclerosis (ALS). Heterogeneous nuclear ribonucleoproteins (hnRNPs) comprise a family of RNA-binding proteins with diverse, multi-functional roles across all aspects of mRNA processing. The role of these proteins in neurodegeneration is far from understood. Here, we review some of the unifying mechanisms by which hnRNPs have been directly or indirectly linked with FTD/ALS pathogenesis, including their incorporation into pathological inclusions and their best-known roles in pre-mRNA splicing regulation. We also discuss the broader functionalities of hnRNPs including their roles in cryptic exon repression, stress granule assembly and in co-ordinating the DNA damage response, which are all emerging pathogenic themes in both diseases. We then present an integrated model that depicts how a broad-ranging network of pathogenic events can arise from declining levels of functional hnRNPs that are inadequately compensated for by autoregulatory means. Finally, we provide a comprehensive overview of the most functionally relevant cellular roles, in the context of FTD/ALS pathogenesis, for hnRNPs A1-U.

RNA-Protein Binding Sites Prediction via Multi Scale Convolutional Gated Recurrent Unit Networks

RNA-Protein binding plays important roles in the field of gene expression. With the development of high throughput sequencing, several conventional methods and deep learning-based methods have been proposed to predict the binding preference of RNA-protein binding. These methods can hardly meet the need of consideration of the dependencies between subsequence and the various motif lengths of different translation factors (TFs). To overcome such limitations, we propose a predictive model that utilizes a combination of multi-scale convolutional layers and bidirectional gated recurrent unit (GRU) layer. Multi-scale convolution layer has the ability to capture the motif features of different lengths, and bidirectional GRU layer is able to capture the dependencies among subsequence. Experimental results show that the proposed method performs better than four state-of-the-art methods in this field. In addition, we investigate the effect of model structure on model performance by performing our proposed method with a different convolution layer and a different number of kernel size. We also demonstrate the effectiveness of bidirectional GRU in improving model performance through comparative experiments.

Aminoacyl-tRNA synthetases in cell signaling

Aminoacyl-tRNA synthetases (ARSs) are a family of essential "housekeeping" enzymes ubiquitous in the three major domains of life. ARSs uniquely connect the essential minimal units of both major oligomer classes-the 3-nucleotide codons of oligonucleotides and the amino acids of proteins. They catalyze the esterification of amino acids to the 3'-end of cognate transfer RNAs (tRNAs) bearing the correct anticodon triplet to ensure accurate transfer of information from mRNA to protein according to the genetic code. As an essential translation factor responsible for the first biochemical reaction in protein biosynthesis, ARSs control protein production by catalyzing aminoacylation, and by editing of mischarged aminoacyl-tRNAs to maintain translational fidelity. In addition to their primary enzymatic activities, many ARSs have noncanonical functions unrelated to their catalytic activity in protein synthesis. Among the ARSs with "moonlighting" activities, several, including GluProRS (or EPRS), LeuRS, LysRS, SerRS, TyrRS, and TrpRS, exhibit cell signaling-related activities that sense environmental signals, regulate gene expression, and modulate cellular functions. ARS signaling functions generally depend on catalytically-inactive, appended domains not present in ancient enzyme forms, and are activated by stimulus-dependent post-translational modification. Activation often results in cellular re-localization and gain of new interacting partners. The newly formed ARS-bearing complexes conduct a host of signal transduction functions, including immune response, mTORC1 pathway signaling, and fibrogenic and angiogenic signaling, among others. Because noncanonical functions of ARSs in signal transduction are uncoupled from canonical aminoacylation functions, function-specific inhibitors can be developed, thus providing promising opportunities and therapeutic targets for treatment of human disease.

DSCAM-AS1-Driven Proliferation of Breast Cancer Cells Involves Regulation of Alternative Exon Splicing and 3'-End Usage

DSCAM-AS1 is a cancer-related long noncoding RNA with higher expression levels in Luminal A, B, and HER2-positive Breast Carcinoma (BC), where its expression is strongly dependent on Estrogen Receptor Alpha (ERα). DSCAM-AS1 expression is analyzed in 30 public datasets and, additionally, by qRT-PCR in tumors from 93 BC patients, to uncover correlations with clinical data. Moreover, the effect of DSCAM-AS1 knockdown on gene expression and alternative splicing is studied by RNA-Seq in MCF-7 cells. We confirm DSCAM-AS1 overexpression in high grade Luminal A, B, and HER2+ BCs and find a significant correlation with disease relapse. In total, 908 genes are regulated by DSCAM-AS1-silencing, primarily involved in the cell cycle and inflammatory response. Noteworthily, the analysis of alternative splicing and isoform regulation reveals 2085 splicing events regulated by DSCAM-AS1, enriched in alternative polyadenylation sites, 3′UTR (untranslated region) shortening and exon skipping events. Finally, the DSCAM-AS1-interacting splicing factor heterogeneous nuclear ribonucleoprotein L (hnRNPL) is predicted as the most enriched RBP for exon skipping and 3′UTR events. The relevance of DSCAM-AS1 overexpression in BC is confirmed by clinical data and further enhanced by its possible involvement in the regulation of RNA processing, which is emerging as one of the most important dysfunctions in cancer.

Upstream analysis of alternative splicing: a review of computational approaches to predict context-dependent splicing factors

Alternative splicing (AS) has shown to play a pivotal role in the development of diseases, including cancer. Specifically, all the hallmarks of cancer (angiogenesis, cell immortality, avoiding immune system response, etc.) are found to have a counterpart in aberrant splicing of key genes. Identifying the context-specific regulators of splicing provides valuable information to find new biomarkers, as well as to define alternative therapeutic strategies. The computational models to identify these regulators are not trivial and require three conceptual steps: the detection of AS events, the identification of splicing factors that potentially regulate these events and the contextualization of these pieces of information for a specific experiment. In this work, we review the different algorithmic methodologies developed for each of these tasks. Main weaknesses and strengths of the different steps of the pipeline are discussed. Finally, a case study is detailed to help the reader be aware of the potential and limitations of this computational approach.

SEQing: web-based visualization of iCLIP and RNA-seq data in an interactive python framework

Background

RNA-binding proteins interact with their target RNAs at specific sites. These binding sites can be determined genome-wide through individual nucleotide resolution crosslinking immunoprecipitation (iCLIP). Subsequently, the binding sites have to be visualized. So far, no visualization tool exists that is easily accessible but also supports restricted access so that data can be shared among collaborators.

Results

Here we present SEQing, a customizable interactive dashboard to visualize crosslink sites on target genes of RNA-binding proteins that have been obtained by iCLIP. Moreover, SEQing supports RNA-seq data that can be displayed in a different window tab. This allows, e.g. crossreferencing the iCLIP data with genes differentially expressed in mutants of the RBP and thus obtain some insights into a potential functional relevance of the binding sites. Additionally, detailed information on the target genes can be incorporated in another tab.

Conclusion

SEQing is written in Python3 and runs on Linux. The web-based access makes iCLIP data easily accessible, even with mobile devices. SEQing is customizable in many ways and has also the option to be secured by a password. The source code is available at https://github.com/malewins/SEQing.

MyoD induced enhancer RNA interacts with hnRNPL to activate target gene transcription during myogenic differentiation

Emerging evidence supports roles of enhancer RNAs (eRNAs) in regulating target gene. Here, we study eRNA regulation and function during skeletal myoblast differentiation. We provide a panoramic view of enhancer transcription and categorization of eRNAs. Master transcription factor MyoD is crucial in activating eRNA production. Super enhancer (se) generated seRNA-1 and -2 promote myogenic differentiation in vitro and in vivo. seRNA-1 regulates expression levels of two nearby genes, myoglobin (Mb) and apolipoprotein L6 (Apol6), by binding to heterogeneous nuclear ribonucleoprotein L (hnRNPL). A CAAA tract on seRNA-1 is essential in mediating seRNA-1/hnRNPL binding and function. Disruption of seRNA-1-hnRNPL interaction attenuates Pol II and H3K36me3 deposition at the Mb locus, in coincidence with the reduction of its transcription. Furthermore, analyses of hnRNPL binding transcriptome-wide reveal its association with eRNAs is a general phenomenon in multiple cells. Collectively, we propose that eRNA-hnRNPL interaction represents a mechanism contributing to target mRNA activation.

Mammalian RNA switches: Molecular rheostats in gene regulation, disease, and medicine

Alteration of RNA structure by environmental signals is a fundamental mechanism of gene regulation. For example, the riboswitch is a noncoding RNA regulatory element that binds a small molecule and causes a structural change in the RNA, thereby regulating transcription, splicing, or translation of an mRNA. The role of riboswitches in metabolite sensing and gene regulation in bacteria and other lower species was reported almost two decades ago, but riboswitches have not yet been discovered in mammals. An analog of the riboswitch, the protein-directed RNA switch (PDRS), has been identified as an important regulatory mechanism of gene expression in mammalian cells. RNA-binding proteins and microRNAs are two major executors of PDRS via their interaction with target transcripts in mammals. These protein-RNA interactions influence cellular functions by integrating environmental signals and intracellular pathways from disparate stimuli to modulate stability or translation of specific mRNAs. The discovery of a riboswitch in eukaryotes that is composed of a single class of thiamine pyrophosphate (TPP) suggests that additional ligand-sensing RNAs may be present to control eukaryotic or mammalian gene expression. In this review, we focus on protein-directed RNA switch mechanisms in mammals. We offer perspectives on the potential discovery of mammalian protein-directed and compound-dependent RNA switches that are related to human disease and medicine.

CLIP and RNA interactome studies to unravel genome-wide RNA-protein interactions in vivo in Arabidopsis thaliana

Post-transcriptional regulation makes an important contribution to adjusting the transcriptome to environmental changes in plants. RNA-binding proteins are key players that interact specifically with mRNAs to co-ordinate their fate. While the regulatory interactions between proteins and RNA are well understood in animals, until recently little information was available on the global binding landscape of RNA-binding proteins in higher plants. This is not least due to technical challenges in plants. In turn, while numerous RNA-binding proteins have been identified through mutant analysis and homology-based searches in plants, only recently a full compendium of proteins with RNA-binding activity has been experimentally determined for the reference plant Arabidopsis thaliana. State-of-the-art techniques to determine RNA-protein interactions genome-wide in animals are based on the covalent fixation of RNA and protein in vivo by UV light. This has only recently been successfully applied to plants. Here, we present practical considerations on the application of UV irradiation based methods to comprehensively determine in vivo RNA-protein interactions in Arabidopsis thaliana, focussing on individual nucleotide resolution crosslinking immunoprecipitation (iCLIP) and mRNA interactome capture.

Heterogeneous nuclear ribonucleoprotein L facilitates recruitment of 53BP1 and BRCA1 at the DNA break sites induced by oxaliplatin in colorectal cancer

Although oxaliplatin is an effective chemotherapeutic drug for treatment of colorectal cancer (CRC), tumor cells can develop mechanisms to evade oxaliplatin-induced cell death and show high tolerance and acquired resistance to this drug. Heterogeneous nuclear ribonucleoprotein L (hnRNP L) has been proved to play a critical role in DNA repair during IgH class switch recombination (CSR) in B lymphocytes, while, its role in CRC and chemotherapeutic resistance remain unknown. Our study aims to uncover an unidentified mechanism of regulating DNA double-strand breaks (DSBs) by hnRNP L in CRC cells treated by oxaliplatin. In present study, we observed that knockdown of hnRNP L enhanced the level of DNA breakage and sensitivity of CRC cells to oxaliplatin. The expression of key DNA repair factors (BRCA1, 53BP1, and ATM) was unaffected by hnRNP L knockdown, thereby excluding the likelihood of hnRNP L mediation via mRNA regulation. Moreover, we observed that phosphorylation level of ATM changed oppositely to 53BP1 and BRCA1 in the CRC cells (SW620 and HCT116) which exhibit synergistic effect by oxaliplatin plus hnRNP L impairment. And similar phenomenon was observed in the foci formation of these critical repair factors. We also found that hnRNP L binds directly with these DNA repair factors through its RNA-recognition motifs (RRMs). Analysis of cell death indicated that the RRMs of hnRNP L are required for cell survival under incubation with oxaliplatin. In conclusion, hnRNP L is critical for the recruitment of the DNA repair factors in oxaliplatin-induced DSBs. Targeting hnRNP L is a promising new clinical approach that could enhance the effectiveness of current chemotherapeutic treatment in patients with resistance to oxaliplatin.

Tandem repeats of the 5′ flanking region of human MUC5AC have a role as a novel enhancer in MUC5AC gene expression

MUC5AC is a well-known gastric differentiation marker, which has been frequently used for the classification of stomach cancer. However, the molecular mechanism of regulation of MUC5AC expression remains to be elucidated. In previous studies, we have shown that Gli regulated MUC5AC transcription through the Gli-binding motif in the 5′ region of MUC5AC. Gli played important roles, but independently was not sufficient for MUC5AC expression. In this study, we analyzed a 4010 bp fragment of the 5′-flanking promoter region of the human MUC5AC gene by luciferase assay, and found a novel distal enhancer region located between −1434 bp to −3000 bp upstream from the first ATG initiation codon. This region is composed of repetitive DNA sequences 5′-TCACTCAC-3′. The strength of enhancer activities depended on the length of the repetitive region. The tandem repeats are conserved among primates, but not in other mammals. The tandem repeat regions enhanced promoter activities not only of MUC5AC but also of other genes. The enhancer effect of the tandem repeat regions was maintained even when inverted. ChIP analysis revealed that H3K9me3 binds to the tandem repeat regions. Together, our results suggest that the tandem repeat region in the MUC5AC promoter has the potential to act as a strong enhancer, and H3K9me3 may contribute to histone modifications of this region.

•

A novel distal enhancer region is located in the MUC5AC promoter.

•

The enhancer region is composed of repetitive DNA sequences 5′-TCACTCAC-3′.

•

H3K9me3 bound to the tandem repeat region in the MUC5AC promoter.

•

Length variants were observed in tandem repeats of the MUC5AC promoter.

hnRNP L-mediated RNA switches function as a hypoxia-induced translational regulon

The GAIT (gamma-interferon-activated inhibitor of translation) complex or miR-297-RISC (RNA-induced silencing complex), together with hnRNP L or hnRNP L-bearing complex, operates an RNA switch in myeloid cells that regulates stress-dependent expression of vascular endothelial growth factor-A (VEGFA). Here, we have shown that hnRNP L directs multiple hypoxia-inducible RNA switches simultaneously and regulates expression of these oncogenic genes in addition to VEGFA. Bioinformatic and polysome profiling-microarray screens have identified DNM1L (Dynamin 1-like) and PHF21A (PHD finger protein 21A) mRNAs as regulated at the translational level by GAIT-dependent, hnRNP L-directed RNA switches. We have also uncovered CDK6 (Cyclin dependent kinase 6), MKLN1 (Muskelin 1) and EIF5 (Eukaryotic initiation factor 5) as novel miR-297-dependent, hnRNP L-directed RNA switch transcripts. Src Kinase is required for the phosphorylation of hnRNP L and activation of the RNA switch pathway. Knockdown of hnRNP L sensitizes the human U937 monocytic cells under hypoxia stress but not in normoxia via inducing cell apoptosis partially due to the reduced translation of hnRNP L target mRNAs. Collectively, our findings suggest that commonly controlled genes by the hnRNP L-directed RNA switches form a translational regulon that promotes hypoxia resistance and cell survival.

RNA structural analysis of the MYC mRNA reveals conserved motifs that affect gene expression

The MYC gene encodes a human transcription factor and proto-oncogene that is dysregulated in over half of all known cancers. To better understand potential post-transcriptional regulatory features affecting MYC expression, we analyzed secondary structures in the MYC mRNA using a program that is optimized for finding small locally-folded motifs with a high propensity for function. This was accomplished by calculating folding metrics across the MYC sequence using a sliding analysis window and generating unique consensus base pairing models weighted by their lower-than-random predicted folding energy. A series of 30 motifs were identified, primarily in the 5' and 3' untranslated regions, which show evidence of structural conservation and compensating mutations across vertebrate MYC homologs. This analysis was able to recapitulate known elements found within an internal ribosomal entry site, as well as discover a novel element in the 3' UTR that is unusually stable and conserved. This novel motif was shown to affect MYC expression, potentially via the modulation of miRNA target accessibility or other trans-regulatory factors. In addition to providing basic insights into mechanisms that regulate MYC expression, this study provides numerous, potentially druggable RNA targets for the MYC gene, which is considered “undruggable” at the protein level.

A mechanism underlying position-specific regulation of alternative splicing

Many RNA-binding proteins including a master regulator of splicing in developing brain and muscle, polypyrimidine tract-binding protein 1 (PTBP1), can either activate or repress alternative exons depending on the pre-mRNA recruitment position. When bound upstream or within regulated exons PTBP1 tends to promote their skipping, whereas binding to downstream sites often stimulates inclusion. How this switch is orchestrated at the molecular level is poorly understood. Using bioinformatics and biochemical approaches we show that interaction of PTBP1 with downstream intronic sequences can activate natural cassette exons by promoting productive docking of the spliceosomal U1 snRNP to a suboptimal 5' splice site. Strikingly, introducing upstream PTBP1 sites to this circuitry leads to a potent splicing repression accompanied by the assembly of an exonic ribonucleoprotein complex with a tightly bound U1 but not U2 snRNP. Our data suggest a molecular mechanism underlying the transition between a better-known repressive function of PTBP1 and its role as a bona fide splicing activator. More generally, we argue that the functional outcome of individual RNA contacts made by an RNA-binding protein is subject to extensive context-specific modulation.

ALYREF links 3'-end processing to nuclear export of non-polyadenylated mRNAs

The RNA-binding protein ALYREF plays key roles in nuclear export and also 3'-end processing of polyadenylated mRNAs, but whether such regulation also extends to non-polyadenylated RNAs is unknown. Replication-dependent (RD)-histone mRNAs are not polyadenylated, but instead end in a stem-loop (SL) structure. Here, we demonstrate that ALYREF prevalently binds a region next to the SL on RD-histone mRNAs. SL-binding protein (SLBP) directly interacts with ALYREF and promotes its recruitment. ALYREF promotes histone pre-mRNA 3'-end processing by facilitating U7-snRNP recruitment through physical interaction with the U7-snRNP-specific component Lsm11. Furthermore, ALYREF, together with other components of the TREX complex, enhances histone mRNA export. Moreover, we show that 3'-end processing promotes ALYREF recruitment and histone mRNA export. Together, our results point to an important role of ALYREF in coordinating 3'-end processing and nuclear export of non-polyadenylated mRNAs.

HNRNPL Restrains miR-155 Targeting of BUB1 to Stabilize Aberrant Karyotypes of Transformed Cells in Chronic Lymphocytic Leukemia

Aneuploidy and overexpression of hsa-miR-155-5p (miR-155) characterize most solid and hematological malignancies. We recently demonstrated that miR-155 sustains aneuploidy at early stages of in vitro cellular transformation. During in vitro transformation of normal human fibroblast, upregulation of miR-155 downregulates spindle checkpoint proteins as the mitotic checkpoint serine/threonine kinase budding uninhibited by benzimidazoles 1 (BUB1), the centromere protein F (CENPF) and the zw10 kinetochore protein (ZW10), compromising the chromosome alignment at the metaphase plate and leading to aneuploidy in daughter cells. Here we show that the heterogeneous nuclear ribonucleoprotein L (HNRNPL) binds to the polymorphic marker D2S1888 at the 3'UTR of BUB1 gene, impairs the miR-155 targeting, and restores BUB1 expression in chronic lymphocytic leukemia. This mechanism occurs at advanced passages of cell transformation and allows the expansion of more favorable clones. Our findings have revealed, at least in part, the molecular mechanisms behind the chromosomal stabilization of cell lines and the concept that, to survive, tumor cells cannot continuously change their genetic heritage but need to stabilize the most suitable karyotype.

Interplay between coding and exonic splicing regulatory sequences

The inclusion of exons during the splicing process depends on the binding of splicing factors to short low-complexity regulatory sequences. The relationship between exonic splicing regulatory sequences and coding sequences is still poorly understood. We demonstrate that exons that are coregulated by any given splicing factor share a similar nucleotide composition bias and preferentially code for amino acids with similar physicochemical properties because of the nonrandomness of the genetic code. Indeed, amino acids sharing similar physicochemical properties correspond to codons that have the same nucleotide composition bias. In particular, we uncover that the TRA2A and TRA2B splicing factors that bind to adenine-rich motifs promote the inclusion of adenine-rich exons coding preferentially for hydrophilic amino acids that correspond to adenine-rich codons. SRSF2 that binds guanine/cytosine-rich motifs promotes the inclusion of GC-rich exons coding preferentially for small amino acids, whereas SRSF3 that binds cytosine-rich motifs promotes the inclusion of exons coding preferentially for uncharged amino acids, like serine and threonine that can be phosphorylated. Finally, coregulated exons encoding amino acids with similar physicochemical properties correspond to specific protein features. In conclusion, the regulation of an exon by a splicing factor that relies on the affinity of this factor for specific nucleotide(s) is tightly interconnected with the exon-encoded physicochemical properties. We therefore uncover an unanticipated bidirectional interplay between the splicing regulatory process and its biological functional outcome.

Role of Alternative Splicing in Prostate Cancer Aggressiveness and Drug Resistance in African Americans

Alternative splicing, the process of removing introns and joining exons of pre-mRNA, is critical for growth, development, tissue homeostasis, and species diversity. Dysregulation of alternative splicing can initiate and drive disease. Aberrant alternative splicing has been shown to promote the "hallmarks of cancer" in both hematological and solid cancers. Of interest, recent work has focused on the role of alternative splicing in prostate cancer and prostate cancer health disparities. We will provide a review of prostate cancer health disparities involving the African American population, alternative RNA splicing, and alternative splicing in prostate cancer. Lastly, we will summarize our work on differential alternative splicing in prostate cancer disparities and its implications for disparate health outcomes and therapeutic targets.

Platforms for Investigating LncRNA Functions

Prior to the sequencing of the human genome, it was presumed that most of the DNA coded for proteins. However, with the advent of next-generation sequencing, it has now been recognized that most complex eukaryotic genomes are in fact transcribed into noncoding RNAs (ncRNAs), including a family of transcripts referred to as long noncoding RNAs (lncRNAs). LncRNAs have been implicated in many biological processes ranging from housekeeping functions such as transcription to more specialized functions such as dosage compensation or genomic imprinting, among others. Interestingly, lncRNAs are not limited to a defined set of functions but can regulate varied activities such as messenger RNA degradation, translation, and protein kinetics or function as RNA decoys or scaffolds. Although still in its infancy, research into the biology of lncRNAs has demonstrated the importance of lncRNAs in development and disease. However, the specific mechanisms through which these lncRNAs act remain poorly defined. Focused research into a small number of these lncRNAs has provided important clues into the heterogeneous nature of this family of ncRNAs. Due to the complex diversity of lncRNA function, in this review, we provide an update on the platforms available for investigators to aid in the identification of lncRNA function.

Staufen1 Interacts with Multiple Components of the Ebola Virus Ribonucleoprotein and Enhances Viral RNA Synthesis

Ebola virus (EBOV) is a negative-strand RNA virus with significant public health importance. Currently, no therapeutics are available for Ebola, which imposes an urgent need for a better understanding of EBOV biology. Here we dissected the virus-host interplay between EBOV and host RNA-binding proteins. We identified novel EBOV host factors, including Staufen1, which interacts with multiple viral factors and is required for efficient viral RNA synthesis.

Ebola virus (EBOV) genome and mRNAs contain long, structured regions that could hijack host RNA-binding proteins to facilitate infection. We performed RNA affinity chromatography coupled with mass spectrometry to identify host proteins that bind to EBOV RNAs and identified four high-confidence proviral host factors, including Staufen1 (STAU1), which specifically binds both 3′ and 5′ extracistronic regions of the EBOV genome. We confirmed that EBOV infection rate and production of infectious particles were significantly reduced in STAU1-depleted cells. STAU1 was recruited to sites of EBOV RNA synthesis upon infection and enhanced viral RNA synthesis. Furthermore, STAU1 interacts with EBOV nucleoprotein (NP), virion protein 30 (VP30), and VP35; the latter two bridge the viral polymerase to the NP-coated genome, forming the viral ribonucleoprotein (RNP) complex. Our data indicate that STAU1 plays a critical role in EBOV replication by coordinating interactions between the viral genome and RNA synthesis machinery.

Coordinate regulation of alternative pre-mRNA splicing events by the human RNA chaperone proteins hnRNPA1 and DDX5

Alternative premessenger RNA (pre-mRNA) splicing is a post-transcriptional mechanism for controlling gene expression. Splicing patterns are determined by both RNA-binding proteins and nuclear pre-mRNA structure. Here, we analyzed pre-mRNA splicing patterns, RNA-binding sites, and RNA structures near these binding sites coordinately controlled by two splicing factors: the heterogeneous nuclear ribonucleoprotein hnRNPA1 and the RNA helicase DDX5. We identified thousands of alternative pre-mRNA splicing events controlled by these factors by RNA sequencing (RNA-seq) following RNAi. Enhanced cross-linking and immunoprecipitation (eCLIP) on nuclear extracts was used to identify protein-RNA-binding sites for both proteins in the nuclear transcriptome. We found a significant overlap between hnRNPA1 and DDX5 splicing targets and that they share many closely linked binding sites as determined by eCLIP analysis. In vivo SHAPE (selective 2'-hydroxyl acylation analyzed by primer extension) chemical RNA structure probing data were used to model RNA structures near several exons controlled and bound by both proteins. Both sequence motifs and in vivo UV cross-linking sites for hnRNPA1 and DDX5 were used to map binding sites in their RNA targets, and often these sites flanked regions of higher chemical reactivity, suggesting an organized nature of nuclear pre-mRNPs. This work provides a first glimpse into the possible RNA structures surrounding pre-mRNA splicing factor-binding sites.

HnRNP L represses cryptic exons

The fidelity of RNA splicing is regulated by a network of splicing enhancers and repressors, although the rules that govern this process are not yet fully understood. One mechanism that contributes to splicing fidelity is the repression of nonconserved cryptic exons by splicing factors that recognize dinucleotide repeats. We previously identified that TDP-43 and PTBP1/PTBP2 are capable of repressing cryptic exons utilizing UG and CU repeats, respectively. Here we demonstrate that hnRNP L (HNRNPL) also represses cryptic exons by utilizing exonic CA repeats, particularly near the 5'SS. We hypothesize that hnRNP L regulates CA repeat repression for both cryptic exon repression and developmental processes such as T cell differentiation.

Multilevel Differential Control of Hormone Gene Expression Programs by hnRNP L and LL in Pituitary Cells

The pituitary-derived somatolactotrophe GH₃ cells secrete both growth hormone (GH) and prolactin (PRL). We have found that the hnRNP L and L-like (LL) paralogs differentially regulate alternative splicing of genes in these cells. Here, we show that hnRNP L is essential for PRL only, but LL is essential for both PRL and GH production. Transcriptome-wide RNA sequencing (RNA-Seq) analysis indicates that they differentially control groups of hormone or hormone-related genes involved in hormone production/regulation at total transcript and alternative exon levels. Interestingly, hnRNP L also specifically binds and prevents the aberrant usage of a nonconserved CA-rich intron piece of Prl pre-mRNA transcripts, and many others involved in endocrine functions, to prevent mostly cryptic last exons and mRNA truncation. Essential for the full hnRNP L effect on specific exons is a proline-rich region that emerged during evolution in vertebrate hnRNP L only but not LL. Together, our data demonstrate that the hnRNP L and its paralog, LL, differentially control hormone gene expression programs at multiple levels, and hnRNP L in particular is critical for protecting the transcriptome from aberrant usage of intronic sequences. The multilevel differential control by hnRNPs likely tailors the transcriptome to help refine and safeguard the different gene expression programs for different hormones.

Plant Ribonomics: Proteins in Search of RNA Partners

Research into the regulation of gene expression underwent a shift from focusing on DNA-binding proteins as key transcriptional regulators to RNA-binding proteins (RBPs) that come into play once transcription has been initiated. RBPs orchestrate all RNA-processing steps in the cell. To obtain a global view of in vivo targets, the RNA complement associated with particular RBPs is determined via immunoprecipitation of the RBP and subsequent identification of bound RNAs via RNA-seq. Here, we describe technical advances in identifying RBP in vivo targets and their binding motifs. We provide an up-to-date view of targets of nucleocytoplasmic RBPs collected in arabidopsis. We also discuss current experimental limitations and provide an outlook on how the approaches may advance our understanding of post-transcriptional networks.

Alternative Splicing in Neurogenesis and Brain Development

Alternative splicing of precursor mRNA is an important mechanism that increases transcriptomic and proteomic diversity and also post-transcriptionally regulates mRNA levels. Alternative splicing occurs at high frequency in brain tissues and contributes to every step of nervous system development, including cell-fate decisions, neuronal migration, axon guidance, and synaptogenesis. Genetic manipulation and RNA sequencing have provided insights into the molecular mechanisms underlying the effects of alternative splicing in stem cell self-renewal and neuronal fate specification. Timely expression and perhaps post-translational modification of neuron-specific splicing regulators play important roles in neuronal development. Alternative splicing of many key transcription regulators or epigenetic factors reprograms the transcriptome and hence contributes to stem cell fate determination. During neuronal differentiation, alternative splicing also modulates signaling activity, centriolar dynamics, and metabolic pathways. Moreover, alternative splicing impacts cortical lamination and neuronal development and function. In this review, we focus on recent progress toward understanding the contributions of alternative splicing to neurogenesis and brain development, which has shed light on how splicing defects may cause brain disorders and diseases.

Systems Approaches to Map In Vivo RNA–Protein Interactions in Arabidopsis thaliana

Proteins that specifically interact with mRNAs orchestrate mRNA processing steps all the way from transcription to decay. Thus, these RNA-binding proteins represent an important control mechanism to double check which proportion of nascent pre-mRNAs is ultimately available for translation into distinct proteins. Here, we discuss recent progress to obtain a systems-level understanding of in vivo RNA–protein interactions in the reference plant Arabidopsis thaliana using protein-centric and RNA-centric methods as well as combined protein binding site and structure probing.

hnRNP L controls HPV16 RNA polyadenylation and splicing in an Akt kinase-dependent manner

Inhibition of the Akt kinase activates HPV16 late gene expression by reducing HPV16 early polyadenylation and by activating HPV16 late L1 mRNA splicing. We identified ‘hot spots’ for RNA binding proteins at the early polyA signal and at splice sites on HPV16 late mRNAs. We observed that hnRNP L was associated with sequences at all HPV16 late splice sites and at the early polyA signal. Akt kinase inhibition resulted in hnRNP L dephosphorylation and reduced association of hnRNP L with HPV16 mRNAs. This was accompanied by an increased binding of U2AF65 and Sam68 to HPV16 mRNAs. Furthermore, siRNA knock-down of hnRNP L or Akt induced HPV16 gene expression. Treatment of HPV16 immortalized keratinocytes with Akt kinase inhibitor reduced hnRNP L binding to HPV16 mRNAs and induced HPV16 L1 mRNA production. Finally, deletion of the hnRNP L binding sites in HPV16 subgenomic expression plasmids resulted in activation of HPV16 late gene expression. In conclusion, the Akt kinase inhibits HPV16 late gene expression at the level of RNA processing by controlling the RNA-binding protein hnRNP L. We speculate that Akt kinase activity upholds an intracellular milieu that favours HPV16 early gene expression and suppresses HPV16 late gene expression.

ALYREF mainly binds to the 5' and the 3' regions of the mRNA in vivo

The TREX complex (TREX) plays key roles in nuclear export of mRNAs. However, little is known about its transcriptome-wide binding targets. We used individual cross-linking and immunoprecipitation (iCLIP) to identify the binding sites of ALYREF, an mRNA export adaptor in TREX, in human cells. Consistent with previous in vitro studies, ALYREF binds to a region near the 5′ end of the mRNA in a CBP80-dependent manner. Unexpectedly, we identified PABPN1-dependent ALYREF binding near the 3′ end of the mRNA. Furthermore, the 3′ processing factor CstF64 directly interacts with ALYREF and is required for the overall binding of ALYREF on the mRNA. In addition, we found that numerous middle exons harbor ALYREF binding sites and identified ALYREF-binding motifs that promote nuclear export of intronless mRNAs. Together, our study defines enrichment of ALYREF binding sites at the 5′ and the 3′ regions of the mRNA in vivo, identifies export-promoting ALYREF-binding motifs, and reveals CstF64- and PABPN1-mediated coupling of mRNA nuclear export to 3′ processing.

Interplay between miR-574-3p and hnRNP L regulates VEGFA mRNA translation and tumorigenesis

MicroRNAs (miRNAs) and heterogeneous nuclear ribonucleoproteins (hnRNPs) are families of sequence-specific, posttranscriptional modulators of gene expression. Despite extensive mechanistic and functional studies on both regulatory classes, the interactions and crosstalk between them are largely unexplored. We have reported that competition between miR-297 and hnRNP L to bind a 3΄UTR-localized CA-rich element (CARE) of VEGFA mRNA regulates its translation. Here, we show that translation of VEGFA mRNA in human myeloid cells is dictated by a bi-directional interaction between miR-574-3p, a CA-rich microRNA, and hnRNP L. In normoxia, miR-574-3p, acting as a decoy, binds cytoplasmic hnRNP L and prevents its binding to the CARE and stimulation of VEGFA mRNA translation, simultaneously permitting miR-297-mediated translational silencing. However, in hypoxia, cytoplasmic accumulation of Tyr³⁵⁹-phosphorylated hnRNP L sequesters miR-574-3p, overcoming its decoy activity and seed sequence-dependent gene silencing activity. Ectopically expressed miR-574-3p binds multiple RNA recognition motif (RRM) domains of hnRNP L, synergizes with miR-297, reduces VEGFA mRNA translation, and triggers apoptosis, thereby suppressing tumorigenesis. Our studies establish a novel condition-dependent interplay between a miRNA and an hnRNP that regulates their functions in a bidirectional manner.

The determinants of alternative RNA splicing in human cells

Alternative splicing represents an important level of the regulation of gene function in eukaryotic organisms. It plays a critical role in virtually every biological process within an organism, including regulation of cell division and cell death, differentiation of tissues in the embryo and the adult organism, as well as in cellular response to diverse environmental factors. In turn, studies of the last decade have shown that alternative splicing itself is controlled by different mechanisms. Unfortunately, there is no clear understanding of how these diverse mechanisms, or determinants, regulate and constrain the set of alternative RNA species produced from any particular gene in every cell of the human body. Here, we provide a consolidated overview of alternative splicing determinants including RNA-protein interactions, epigenetic regulation via chromatin remodeling, coupling of transcription-to-alternative splicing, effect of secondary structures in pre-RNA, and function of the RNA quality control systems. We also extensively and critically discuss some mechanistic insights on coordinated inclusion/exclusion of exons during the formation of mature RNA molecules. We conclude that the final structure of RNA is pre-determined by a complex interplay between cis- and trans-acting factors. Altogether, currently available empirical data significantly expand our understanding of the functioning of the alternative splicing machinery of cells in normal and pathological conditions. On the other hand, there are still many blind spots that require further deep investigations.

Genome-wide CRISPR screen identifies HNRNPL as a prostate cancer dependency regulating RNA splicing

Alternative RNA splicing plays an important role in cancer. To determine which factors involved in RNA processing are essential in prostate cancer, we performed a genome-wide CRISPR/Cas9 knockout screen to identify the genes that are required for prostate cancer growth. Functional annotation defined a set of essential spliceosome and RNA binding protein (RBP) genes, including most notably heterogeneous nuclear ribonucleoprotein L (HNRNPL). We defined the HNRNPL-bound RNA landscape by RNA immunoprecipitation coupled with next-generation sequencing and linked these RBP-RNA interactions to changes in RNA processing. HNRNPL directly regulates the alternative splicing of a set of RNAs, including those encoding the androgen receptor, the key lineage-specific prostate cancer oncogene. HNRNPL also regulates circular RNA formation via back splicing. Importantly, both HNRNPL and its RNA targets are aberrantly expressed in human prostate tumors, supporting their clinical relevance. Collectively, our data reveal HNRNPL and its RNA clients as players in prostate cancer growth and potential therapeutic targets.

Perspectives of long non-coding RNAs in cancer

A recent advance in transcriptomics has spawned the 'Decade of non-coding RNAs' by potentiating the growing numbers of long non-coding RNA in cancer. LncRNA involvement in cancer denotes its significance beyond our perception as they participate in tumor suppression and promoting oncogenesis, which raises them as a mighty class of effectors or regulators. Aberrantly expressed lncRNAs interact with major protein and coding partners, which ultimately deregulate normal cellular processes and drive the cell towards malignant state. Identification of theses interactions are utmost important as lncRNAs can be ideal targets for therapy. Dysregulation of lncRNAs by genomic alterations like single nucleotide variations and gene fusions are also potential modulators of their secondary structure. In this review, we discuss the various molecular interactions of lncRNAs with major bio-molecules and genetic variations in lncRNA genes and their importance in cancer. This systematic review outlines the vivid role of lncRNAs in cancer context and opens up future conceptual applications.

Alternative splicing: the pledge, the turn, and the prestige : The key role of alternative splicing in human biological systems

Alternative pre-mRNA splicing is a tightly controlled process conducted by the spliceosome, with the assistance of several regulators, resulting in the expression of different transcript isoforms from the same gene and increasing both transcriptome and proteome complexity. The differences between alternative isoforms may be subtle but enough to change the function or localization of the translated proteins. A fine control of the isoform balance is, therefore, needed throughout developmental stages and adult tissues or physiological conditions and it does not come as a surprise that several diseases are caused by its deregulation. In this review, we aim to bring the splicing machinery on stage and raise the curtain on its mechanisms and regulation throughout several systems and tissues of the human body, from neurodevelopment to the interactions with the human microbiome. We discuss, on one hand, the essential role of alternative splicing in assuring tissue function, diversity, and swiftness of response in these systems or tissues, and on the other hand, what goes wrong when its regulatory mechanisms fail. We also focus on the possibilities that splicing modulation therapies open for the future of personalized medicine, along with the leading techniques in this field. The final act of the spliceosome, however, is yet to be fully revealed, as more knowledge is needed regarding the complex regulatory network that coordinates alternative splicing and how its dysfunction leads to disease.