Intronic splicing enhancers, cognate splicing factors and context-dependent regulation rules

Placental co-transcriptional activator Vestigial-like 1 (VGLL1) drives tumorigenesis via increasing transcription of proliferation and invasion genes

Introduction

Vestigial-like 1 (VGLL1) is a co-transcriptional activator that binds to TEA domain-containing transcription factors (TEADs). Its expression is upregulated in a variety of aggressive cancer types, including pancreatic and basal-like breast cancer, and increased transcription of VGLL1 is strongly correlated with poor prognosis and decreased overall patient survival. In normal tissues, VGLL1 is most highly expressed within placental trophoblast cells, which share the common attributes of rapid cellular proliferation and invasion with tumor cells. The impact of VGLL1 in cancer has not been fully elucidated and no VGLL1-targeted therapy currently exists.

Methods

The aim of this study was to evaluate the cellular function and downstream genomic targets of VGLL1 in placental, pancreatic, and breast cancer cells. Functional assays were employed to assess the role of VGLL1 in cellular invasion and proliferation, and ChIP-seq and RNAseq assays were performed to identify VGLL1 target genes and potential impact using pathway analysis.

Results

ChIP-seq analysis identified eight transcription factors with a VGLL1-binding motif that were common between all three cell types, including TEAD1-4, AP-1, and GATA6, and revealed ~3,000 shared genes with which VGLL1 interacts. Furthermore, increased VGLL1 expression led to an enhancement of cell invasion and proliferation, which was supported by RNAseq analysis showing transcriptional changes in several genes known to be involved in these processes.

Discussion

This work expands our mechanistic understanding of VGLL1 function in tumor cells and provides a strong rationale for developing VGLL1-targeted therapies for treating cancer patients.

The Many Roads from Alternative Splicing to Cancer: Molecular Mechanisms Involving Driver Genes

Cancer driver genes are either oncogenes or tumour suppressor genes that are classically activated or inactivated, respectively, by driver mutations. Alternative splicing-which produces various mature mRNAs and, eventually, protein variants from a single gene-may also result in driving neoplastic transformation because of the different and often opposed functions of the variants of driver genes. The present review analyses the different alternative splicing events that result in driving neoplastic transformation, with an emphasis on their molecular mechanisms. To do this, we collected a list of 568 gene drivers of cancer and revised the literature to select those involved in the alternative splicing of other genes as well as those in which its pre-mRNA is subject to alternative splicing, with the result, in both cases, of producing an oncogenic isoform. Thirty-one genes fall into the first category, which includes splicing factors and components of the spliceosome and splicing regulators. In the second category, namely that comprising driver genes in which alternative splicing produces the oncogenic isoform, 168 genes were found. Then, we grouped them according to the molecular mechanisms responsible for alternative splicing yielding oncogenic isoforms, namely, mutations in cis splicing-determining elements, other causes involving non-mutated cis elements, changes in splicing factors, and epigenetic and chromatin-related changes. The data given in the present review substantiate the idea that aberrant splicing may regulate the activation of proto-oncogenes or inactivation of tumour suppressor genes and details on the mechanisms involved are given for more than 40 driver genes.

An interpretable model of pre-mRNA splicing for animal and plant genes

Pre-mRNA splicing is a fundamental step in gene expression, conserved across eukaryotes, in which the spliceosome recognizes motifs at the 3' and 5' splice sites (SSs), excises introns, and ligates exons. SS recognition and pairing is often influenced by protein splicing factors (SFs) that bind to splicing regulatory elements (SREs). Here, we describe SMsplice, a fully interpretable model of pre-mRNA splicing that combines models of core SS motifs, SREs, and exonic and intronic length preferences. We learn models that predict SS locations with 83 to 86% accuracy in fish, insects, and plants and about 70% in mammals. Learned SRE motifs include both known SF binding motifs and unfamiliar motifs, and both motif classes are supported by genetic analyses. Our comparisons across species highlight similarities between non-mammals, increased reliance on intronic SREs in plant splicing, and a greater reliance on SREs in mammalian splicing.

Conditional splicing system for tight control of viral overlapping genes

Viral genomes frequently harbor overlapping genes, complicating the development of virus-vectored vaccines and gene therapies. This study introduces a novel conditional splicing system to precisely control the expression of such overlapping genes through recombinase-mediated conditional splicing. We refined site-specific recombinase (SSR) conditional splicing systems and explored their mechanisms. The systems demonstrated exceptional inducibility (116,700-fold increase) with negligible background expression, facilitating the conditional expression of overlapping genes in adenovirus-associated virus (AAV) and human immunodeficiency virus type 1. Notably, this approach enabled the establishment of stable AAV producer cell lines, encapsulating all necessary packaging genes. Our findings underscore the potential of the SSR-conditional splicing system to significantly advance vector engineering, enhancing the efficacy and scalability of viral-vector-based therapies and vaccines.

IMPORTANCE

Regulating overlapping genes is vital for gene therapy and vaccine development using viral vectors. The regulation of overlapping genes presents challenges, including cytotoxicity and impacts on vector capacity and genome stability, which restrict stable packaging cell line development and broad application. To address these challenges, we present a "loxp-splice-loxp"-based conditional splicing system, offering a novel solution for conditional expression of overlapping genes and stable cell line establishment. This system may also regulate other cytotoxic genes, representing a significant advancement in cell engineering and gene therapy as well as biomass production.

Role of epithelial splicing regulatory protein 1 in cancer progression

As aberrant alternative splicing by either dysregulation or mutations of splicing factors contributes to cancer initiation and progression, splicing factors are emerging as potential therapeutic targets for cancer therapy. Therefore, pharmacological modulators targeting splicing factors have been under development. Epithelial splicing regulatory protein 1 (ESRP1) is an epithelial cell-specific splicing factor, whose downregulation is associated with epithelial-mesenchymal transition (EMT) by regulating alternative splicing of multiple genes, such as CD44, CTNND1, ENAH, and FGFR2. Consistent with the downregulation of ESRP1 during EMT, it has been initially revealed that high ESRP1 expression is associated with favorable prognosis and ESRP1 plays a tumor-suppressive role in cancer progression. However, ESRP1 has been found to promote cancer progression in some cancers, such as breast and ovarian cancers, indicating that it plays a dual role in cancer progression depending on the type of cancer. Furthermore, recent studies have reported that ESRP1 affects tumor growth by regulating the metabolism of tumor cells or immune cell infiltration in the tumor microenvironment, suggesting the novel roles of ESRP1 in addition to EMT. ESRP1 expression was also associated with response to anticancer drugs. This review describes current understanding of the roles and mechanisms of ESRP1 in cancer progression, and further discusses the emerging novel roles of ESRP1 in cancer and recent attempts to target splicing factors for cancer therapy.

Adenine base editor-mediated splicing remodeling activates noncanonical splice sites

Adenine base editors (ABEs) are genome-editing tools that have been harnessed to introduce precise A•T to G•C conversion. The discovery of split genes revealed that all introns contain two highly conserved dinucleotides, canonical "AG" (acceptor) and "GT" (donor) splice sites. ABE can directly edit splice acceptor sites of the adenine (A) base, leading to aberrant gene splicing, which may be further adopted to remodel splicing. However, spliced isoforms triggered with ABE have not been well explored. To address it, we initially generated a cell line harboring C-terminal enhanced GFP (eGFP)-tagged β-actin (ACTB), in which the eGFP signal can track endogenous β-actin expression. Expectedly, after the editing of splice acceptor sites, we observed a dramatical decrease in the percentage of eGFP-positive cells and generation of splicing products with the noncanonical splice site. Furthermore, we manipulated Peroxidasin in mouse embryos with ABE, in which a noncanonical acceptor was activated to remodel splicing, successfully generating a mouse disease model of anophthalmia and severely malformed microphthalmia. Collectively, we demonstrate that ABE-mediated splicing remodeling can activate a noncanonical acceptor to manipulate human and mouse genomes, which will facilitate the investigation of basic and translational medicine studies.

The human genome contains over a million autonomous exons

Mammalian mRNA and lncRNA exons are often small compared to introns. The exon definition model predicts that exons splice autonomously, dependent on proximal exon sequence features, explaining their delineation within large introns. This model has not been examined on a genome-wide scale, however, leaving open the question of how often mRNA and lncRNA exons are autonomous. It is also unknown how frequently such exons can arise by chance. Here, we directly assayed large fragments (500-1000 bp) of the human genome by exon trapping, which detects exons spliced into a heterologous transgene, here designed with a large intron context. We define the trapped exons as "autonomous." We obtained ∼1.25 million trapped exons, including most known mRNA and well-annotated lncRNA internal exons, demonstrating that human exons are predominantly autonomous. mRNA exons are trapped with the highest efficiency. Nearly a million of the trapped exons are unannotated, most located in intergenic regions and antisense to mRNA, with depletion from the forward strand of introns. These exons are not conserved, suggesting they are nonfunctional and arose from random mutations. They are nonetheless highly enriched with known splicing promoting sequence features that delineate known exons. Novel autonomous exons are more numerous than annotated lncRNA exons, and computational models also indicate they will occur with similar frequency in any randomly generated sequence. These results show that most human coding exons splice autonomously, and provide an explanation for the existence of many unconserved lncRNAs, as well as a new annotation and inclusion levels of spliceable loci in the human genome.

The role of alternative pre-mRNA splicing in cancer progression

Alternative pre-mRNA splicing is a critical mechanism that generates multiple mRNA from a single gene, thereby increasing the diversity of the proteome. Recent research has highlighted the significance of specific splicing isoforms in cellular processes, particularly in regulating cell numbers. In this review, we examine the current understanding of the role of alternative splicing in controlling cancer cell growth and discuss specific splicing factors and isoforms and their molecular mechanisms in cancer progression. These isoforms have been found to intricately control signaling pathways crucial for cell cycle progression, proliferation, and apoptosis. Furthermore, studies have elucidated the characteristics and functional importance of splicing factors that influence cell numbers. Abnormal expression of oncogenic splicing isoforms and splicing factors, as well as disruptions in splicing caused by genetic mutations, have been implicated in the development and progression of tumors. Collectively, these findings provide valuable insights into the complex interplay between alternative splicing and cell proliferation, thereby suggesting the potential of alternative splicing as a therapeutic target for cancer.

The nexus between RNA-binding proteins and their effectors

RNA-binding proteins (RBPs) regulate essentially every event in the lifetime of an RNA molecule, from its production to its destruction. Whereas much has been learned about RNA sequence specificity and general functions of individual RBPs, the ways in which numerous RBPs instruct a much smaller number of effector molecules, that is, the core engines of RNA processing, as to where, when and how to act remain largely speculative. Here, we survey the known modes of communication between RBPs and their effectors with a particular focus on converging RBP-effector interactions and their roles in reducing the complexity of RNA networks. We discern the emerging unifying principles and discuss their utility in our understanding of RBP function, regulation of biological processes and contribution to human disease.

An iterative gene-editing strategy broadens eIF4E1 genetic diversity in Solanum lycopersicum and generates resistance to multiple potyvirus isolates

Resistance to potyviruses in plants has been largely provided by the selection of natural variant alleles of eukaryotic translation initiation factors (eIF) 4E in many crops. However, the sources of such variability for breeding can be limited for certain crop species, while new virus isolates continue to emerge. Different methods of mutagenesis have been applied to inactivate the eIF4E genes to generate virus resistance, but with limited success due to the physiological importance of translation factors and their redundancy. Here, we employed genome editing approaches at the base level to induce non-synonymous mutations in the eIF4E1 gene and create genetic diversity in cherry tomato (Solanum lycopersicum var. cerasiforme). We sequentially edited the genomic sequences coding for two regions of eIF4E1 protein, located around the cap-binding pocket and known to be important for susceptibility to potyviruses. We show that the editing of only one of the two regions, by gene knock-in and base editing, respectively, is not sufficient to provide resistance. However, combining amino acid mutations in both regions resulted in resistance to multiple potyviruses without affecting the functionality in translation initiation. Meanwhile, we report that extensive base editing in exonic region can alter RNA splicing pattern, resulting in gene knockout. Altogether our work demonstrates that precision editing allows to design plant factors based on the knowledge on evolutionarily selected alleles and enlarge the gene pool to potentially provide advantageous phenotypes such as pathogen resistance.

Mechanism and modeling of human disease-associated near-exon intronic variants that perturb RNA splicing

It is estimated that 10%-30% of disease-associated genetic variants affect splicing. Splicing variants may generate deleteriously altered gene product and are potential therapeutic targets. However, systematic diagnosis or prediction of splicing variants is yet to be established, especially for the near-exon intronic splice region. The major challenge lies in the redundant and ill-defined branch sites and other splicing motifs therein. Here, we carried out unbiased massively parallel splicing assays on 5,307 disease-associated variants that overlapped with branch sites and collected 5,884 variants across the 5' splice region. We found that strong splice sites and exonic features preserve splicing from intronic sequence variation. Whereas the splice-altering mechanism of the 3' intronic variants is complex, that of the 5' is mainly splice-site destruction. Statistical learning combined with these molecular features allows precise prediction of altered splicing from an intronic variant. This statistical model provides the identity and ranking of biological features that determine splicing, which serves as transferable knowledge and out-performs the benchmarking predictive tool. Moreover, we demonstrated that intronic splicing variants may associate with disease risks in the human population. Our study elucidates the mechanism of splicing response of intronic variants, which classify disease-associated splicing variants for the promise of precision medicine.

hnRNP A1 in RNA metabolism regulation and as a potential therapeutic target

Abnormal RNA metabolism, regulated by various RNA binding proteins, can have functional consequences for multiple diseases. Heterogeneous nuclear ribonucleoprotein A1 (hnRNP A1) is an important RNA binding protein, that regulates various RNA metabolic processes, including transcription, alternative splicing of pre-mRNA, translation, miRNA processing and mRNA stability. As a potent splicing factor, hnRNP A1 can regulate multiple splicing events, including itself, collaborating with other cooperative or antagonistical splicing factors by binding to splicing sites and regulatory elements in exons or introns. hnRNP A1 can modulate gene transcription by directly interacting with promoters or indirectly impacting Pol II activities. Moreover, by interacting with the internal ribosome entry site (IRES) or 3′-UTR of mRNAs, hnRNP A1 can affect mRNA translation. hnRNP A1 can alter the stability of mRNAs by binding to specific locations of 3′-UTR, miRNAs biogenesis and Nonsense-mediated mRNA decay (NMD) pathway. In this review, we conclude the selective sites where hnRNP A1 binds to RNA and DNA, and the co-regulatory factors that interact with hnRNP A1. Given the dysregulation of hnRNP A1 in diverse diseases, especially in cancers and neurodegeneration diseases, targeting hnRNP A1 for therapeutic treatment is extremely promising. Therefore, this review also provides the small-molecule drugs, biomedicines and novel strategies targeting hnRNP A1 for therapeutic purposes.

A widespread length-dependent splicing dysregulation in cancer

Dysregulation of alternative splicing is a key molecular hallmark of cancer. However, the common features and underlying mechanisms remain unclear. Here, we report an intriguing length-dependent splicing regulation in cancers. By systematically analyzing the transcriptome of thousands of cancer patients, we found that short exons are more likely to be mis-spliced and preferentially excluded in cancers. Compared to other exons, cancer-associated short exons (CASEs) are more conserved and likely to encode in-frame low-complexity peptides, with functional enrichment in GTPase regulators and cell adhesion. We developed a CASE-based panel as reliable cancer stratification markers and strong predictors for survival, which is clinically useful because the detection of short exon splicing is practical. Mechanistically, mis-splicing of CASEs is regulated by elevated transcription and alteration of certain RNA binding proteins in cancers. Our findings uncover a common feature of cancer-specific splicing dysregulation with important clinical implications in cancer diagnosis and therapies.

Pervasive translation of circular RNAs driven by short IRES-like elements

Some circular RNAs (circRNAs) were found to be translated through IRES-driven mechanism, however the scope and functions of circRNA translation are unclear because endogenous IRESs are rare. To determine the prevalence and mechanism of circRNA translation, we develop a cell-based system to screen random sequences and identify 97 overrepresented hexamers that drive cap-independent circRNA translation. These IRES-like short elements are significantly enriched in endogenous circRNAs and sufficient to drive circRNA translation. We further identify multiple trans-acting factors that bind these IRES-like elements to initiate translation. Using mass-spectrometry data, hundreds of circRNA-coded peptides are identified, most of which have low abundance due to rapid degradation. As judged by mass-spectrometry, 50% of translatable endogenous circRNAs undergo rolling circle translation, several of which are experimentally validated. Consistently, mutations of the IRES-like element in one circRNA reduce its translation. Collectively, our findings suggest a pervasive translation of circRNAs, providing profound implications in translation control.

Unbiased screen of random sequences identified many short IRES-like elements to drive circular RNA translation and hundreds of rolling circle translation events, suggesting a pervasive cap-independent translation in human transcriptome.

HNRNPH1 destabilizes the G-quadruplex structures formed by G-rich RNA sequences that regulate the alternative splicing of an oncogenic fusion transcript

In the presence of physiological monovalent cations, thousands of RNA G-rich sequences can form parallel G-quadruplexes (G4s) unless RNA-binding proteins inhibit, destabilize, or resolve the formation of such secondary RNA structures. Here, we have used a disease-relevant model system to investigate the biophysical properties of the RNA-binding protein HNRNPH1's interaction with G-rich sequences. We demonstrate the importance of two EWSR1-exon 8 G-rich regions in mediating the exclusion of this exon from the oncogenic EWS-FLI1 transcripts expressed in a subset of Ewing sarcomas, using complementary analysis of tumor data, long-read sequencing, and minigene studies. We determined that HNRNPH1 binds the EWSR1-exon 8 G-rich sequences with low nM affinities irrespective of whether in a non-G4 or G4 state but exhibits different kinetics depending on RNA structure. Specifically, HNRNPH1 associates and dissociates from G4-folded RNA faster than the identical sequences in a non-G4 state. Importantly, we demonstrate using gel shift and spectroscopic assays that HNRNPH1, particularly the qRRM1-qRRM2 domains, destabilizes the G4s formed by the EWSR1-exon 8 G-rich sequences in a non-catalytic fashion. Our results indicate that HNRNPH1's binding of G-rich sequences favors the accumulation of RNA in a non-G4 state and that this contributes to its regulation of RNA processing.

Quantitative prediction of variant effects on alternative splicing in MAPT using endogenous pre-messenger RNA structure probing

Splicing is highly regulated and is modulated by numerous factors. Quantitative predictions for how a mutation will affect precursor mRNA (pre-mRNA) structure and downstream function are particularly challenging. Here, we use a novel chemical probing strategy to visualize endogenous precursor and mature MAPT mRNA structures in cells. We used these data to estimate Boltzmann suboptimal structural ensembles, which were then analyzed to predict consequences of mutations on pre-mRNA structure. Further analysis of recent cryo-EM structures of the spliceosome at different stages of the splicing cycle revealed that the footprint of the B^act complex with pre-mRNA best predicted alternative splicing outcomes for exon 10 inclusion of the alternatively spliced MAPT gene, achieving 74% accuracy. We further developed a β-regression weighting framework that incorporates splice site strength, RNA structure, and exonic/intronic splicing regulatory elements capable of predicting, with 90% accuracy, the effects of 47 known and 6 newly discovered mutations on inclusion of exon 10 of MAPT. This combined experimental and computational framework represents a path forward for accurate prediction of splicing-related disease-causing variants.

Polymorphic differences within HLA-C alleles contribute to alternatively spliced transcripts lacking exon 5

The HLA genes are amongst the most polymorphic in the human genome. Alternative splicing could add an extra layer of complexity, but has not been studied extensively. Here, we applied an RNA based approach to study the influence of allele polymorphism on alternative splicing of HLA‐C in peripheral blood. RNA was isolated from these peripheral cells, converted into cDNA and amplified specifically for 12 common HLA‐C allele groups. Through subsequent sequencing of HLA‐C, we observed alternative splicing variants of HLA‐C*04 and *16 that resulted in exon 5 skipping and were co‐expressed with the mature transcript. Investigation of intron 4 sequences of HLA‐C*04 and *16 compared with other HLA‐C alleles demonstrated no effect on predicted splice sites and branch point. To further investigate if the unique polymorphic positions in exon 5 of HLA‐C*04 or *16 may facilitate alternative splicing by acting on splicing regulatory elements (SRE), in‐silico splicing analysis was performed. While the HLA‐C*04 specific SNP in exon 5 had no effect on predicted exonic SRE, the HLA‐C*16 specific exon 5 SNP did alter exonic SRE. Our findings provide experimental and theoretical support for the concept that polymorphisms within the HLA‐C alleles influence the alternative splicing of HLA‐C.

p73α1, a p73 C-terminal isoform, regulates tumor suppression and the inflammatory response via Notch1

Significance p73 is expressed as multiple C-terminal isoforms, but their expression and activity are largely unknown. Here, we identified p73α1 as a p73 C-terminal isoform that results from exon 12 (E12) exclusion. We showed that E12 deficiency in mice leads to systemic inflammation but not spontaneous tumors. We also showed that Notch1 is regulated by p73α1 and plays a critical role in p73-dependent tumor suppression and systemic inflammation.

Mechanisms of non-coding RNA-modulated alternative splicing in cancer

Alternative splicing (AS) is a common and pivotal process for eukaryotic gene expression regulation, which enables a precursor RNA to produce multiple transcript variants with diverse cellular functions. Aberrant AS represents a hallmark of cancer, engaged in all stages of tumorigenesis from initiation to metastasis. Accumulating pieces of evidence have revealed the involvement of non-coding RNAs (ncRNAs) in regulating AS in human cancers. In this review, we overview the underlying mechanisms of non-coding RNAs, including microRNAs (miRNAs), long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs) modulated AS at diverse levels in human cancers, and summarize their regulatory functions in tumorigenesis.

Clinical implementation of RNA sequencing for Mendelian disease diagnostics

BACKGROUND

Lack of functional evidence hampers variant interpretation, leaving a large proportion of individuals with a suspected Mendelian disorder without genetic diagnosis after whole genome or whole exome sequencing (WES). Research studies advocate to further sequence transcriptomes to directly and systematically probe gene expression defects. However, collection of additional biopsies and establishment of lab workflows, analytical pipelines, and defined concepts in clinical interpretation of aberrant gene expression are still needed for adopting RNA sequencing (RNA-seq) in routine diagnostics.

METHODS

We implemented an automated RNA-seq protocol and a computational workflow with which we analyzed skin fibroblasts of 303 individuals with a suspected mitochondrial disease that previously underwent WES. We also assessed through simulations how aberrant expression and mono-allelic expression tests depend on RNA-seq coverage.

RESULTS

We detected on average 12,500 genes per sample including around 60% of all disease genes-a coverage substantially higher than with whole blood, supporting the use of skin biopsies. We prioritized genes demonstrating aberrant expression, aberrant splicing, or mono-allelic expression. The pipeline required less than 1 week from sample preparation to result reporting and provided a median of eight disease-associated genes per patient for inspection. A genetic diagnosis was established for 16% of the 205 WES-inconclusive cases. Detection of aberrant expression was a major contributor to diagnosis including instances of 50% reduction, which, together with mono-allelic expression, allowed for the diagnosis of dominant disorders caused by haploinsufficiency. Moreover, calling aberrant splicing and variants from RNA-seq data enabled detecting and validating splice-disrupting variants, of which the majority fell outside WES-covered regions.

CONCLUSION

Together, these results show that streamlined experimental and computational processes can accelerate the implementation of RNA-seq in routine diagnostics.

Pervasive misannotation of microexons that are evolutionarily conserved and crucial for gene function in plants

It is challenging to identify the smallest microexons (≤15-nt) due to their small size. Consequently, these microexons are often misannotated or missed entirely during genome annotation. Here, we develop a pipeline to accurately identify 2,398 small microexons in 10 diverse plant species using 990 RNA-seq datasets, and most of them have not been annotated in the reference genomes. Analysis reveals that microexons tend to have increased detained flanking introns that require post-transcriptional splicing after polyadenylation. Examination of 45 conserved microexon clusters demonstrates that microexons and associated gene structures can be traced back to the origin of land plants. Based on these clusters, we develop an algorithm to genome-wide model coding microexons in 132 plants and find that microexons provide a strong phylogenetic signal for plant organismal relationships. Microexon modeling reveals diverse evolutionary trajectories, involving microexon gain and loss and alternative splicing. Our work provides a comprehensive view of microexons in plants.

The small size (≤15-nt) of micorexons poses difficulties for genome annotation and identification using standard RNA sequence mapping approaches. Here, the authors develop computational pipelines to discover and predict microexons in plants and reveal diverse evolutionary trajectories via genomewide microexon modeling.

Systematic characterization of short intronic splicing-regulatory elements in SMN2 pre-mRNA

Intronic splicing enhancers and silencers (ISEs and ISSs) are two groups of splicing-regulatory elements (SREs) that play critical roles in determining splice-site selection, particularly for alternatively spliced introns or exons. SREs are often short motifs; their mutation or dysregulation of their cognate proteins frequently causes aberrant splicing and results in disease. To date, however, knowledge about SRE sequences and how they regulate splicing remains limited. Here, using an SMN2 minigene, we generated a complete pentamer-sequence library that comprises all possible combinations of 5 nucleotides in intron 7, at a fixed site downstream of the 5′ splice site. We systematically analyzed the effects of all 1023 mutant pentamers on exon 7 splicing, in comparison to the wild-type minigene, in HEK293 cells. Our data show that the majority of pentamers significantly affect exon 7 splicing: 584 of them are stimulatory and 230 are inhibitory. To identify actual SREs, we utilized a motif set enrichment analysis (MSEA), from which we identified groups of stimulatory and inhibitory SRE motifs. We experimentally validated several strong SREs in SMN1/2 and other minigene settings. Our results provide a valuable resource for understanding how short RNA sequences regulate splicing. Many novel SREs can be explored further to elucidate their mechanism of action.

ERG transcription factors have a splicing regulatory function involving RBFOX2 that is altered in the EWS-FLI1 oncogenic fusion

ERG family proteins (ERG, FLI1 and FEV) are a subfamily of ETS transcription factors with key roles in physiology and development. In Ewing sarcoma, the oncogenic fusion protein EWS-FLI1 regulates both transcription and alternative splicing of pre-messenger RNAs. However, whether wild-type ERG family proteins might regulate splicing is unknown. Here, we show that wild-type ERG proteins associate with spliceosomal components, are found on nascent RNAs, and induce alternative splicing when recruited onto a reporter minigene. Transcriptomic analysis revealed that ERG and FLI1 regulate large numbers of alternative spliced exons (ASEs) enriched with RBFOX2 motifs and co-regulated by this splicing factor. ERG and FLI1 are associated with RBFOX2 via their conserved carboxy-terminal domain, which is present in EWS-FLI1. Accordingly, EWS-FLI1 is also associated with RBFOX2 and regulates ASEs enriched in RBFOX2 motifs. However, in contrast to wild-type ERG and FLI1, EWS-FLI1 often antagonizes RBFOX2 effects on exon inclusion. In particular, EWS-FLI1 reduces RBFOX2 binding to the ADD3 pre-mRNA, thus increasing its long isoform, which represses the mesenchymal phenotype of Ewing sarcoma cells. Our findings reveal a RBFOX2-mediated splicing regulatory function of wild-type ERG family proteins, that is altered in EWS-FLI1 and contributes to the Ewing sarcoma cell phenotype.

Therapeutic Modulation of RNA Splicing in Malignant and Non-Malignant Disease

RNA splicing is the enzymatic process by which non-protein coding sequences are removed from RNA to produce mature protein-coding mRNA. Splicing is thereby a major mediator of proteome diversity as well as a dynamic regulator of gene expression. Genetic alterations disrupting splicing of individual genes or altering the function of splicing factors contribute to a wide range of human genetic diseases as well as cancer. These observations have resulted in the development of therapies based on oligonucleotides that bind to RNA sequences and modulate splicing for therapeutic benefit. In parallel, small molecules that bind to splicing factors to alter their function or modify RNA processing of individual transcripts are being pursued for monogenic disorders as well as for cancer.

Regulatory roles and mechanisms of alternative RNA splicing in adipogenesis and human metabolic health

Alternative splicing (AS) regulates gene expression patterns at the post-transcriptional level and generates a striking expansion of coding capacities of genomes and cellular protein diversity. RNA splicing could undergo modulation and close interaction with genetic and epigenetic machinery. Notably, during the adipogenesis processes of white, brown and beige adipocytes, AS tightly interplays with the differentiation gene program networks. Here, we integrate the available findings on specific splicing events and distinct functions of different splicing regulators as examples to highlight the directive biological contribution of AS mechanism in adipogenesis and adipocyte biology. Furthermore, accumulating evidence has suggested that mutations and/or altered expression in splicing regulators and aberrant splicing alterations in the obesity-associated genes are often linked to humans’ diet-induced obesity and metabolic dysregulation phenotypes. Therefore, significant attempts have been finally made to overview novel detailed discussion on the prospects of splicing machinery with obesity and metabolic disorders to supply featured potential management mechanisms in clinical applicability for obesity treatment strategies.

Protocol for proteogenomic dissection of intronic splicing enhancer interactome for prediction of individualized cancer prognosis

Inter- or intra-patient tumor heterogeneity hinders the discovery of biomarkers for predicting individualized prognosis. Here, we present a protocol for an alternative splicing activity-based proteogenomic approach for identification of candidate prognostic markers in cancer cell lines and human breast cancer specimens. The pull-down of protein complexes with intronic splicing enhancer (ISE) probes is followed by tandem mass spectrometry (MS/MS) peptide sequencing. The proteogenomic analysis of data from these ISE-MS/MS assays identifies new prognostic markers that can be utilized to stratify patients with poor prognosis.

For complete details on the use and execution of this protocol, please refer to Wang et al. (2018).

•

Protocol for LC-MS/MS analysis of ISE interactomes from cancer cell lines or tissues

•

Protocol for proteogenomic identification of ISE interactors of prognostic significance

•

Protocol for protein extraction, ISE pull-downs, MS/MS, and proteogenomic analysis

•

Protocols applicable to cancer cell lines and clinical specimens

Alternative splicing of lncRNAs in human diseases

Alternative splicing (AS), a vital post-transcription process for eukaryote gene expression regulating, can efficiently improve gene utilization and increase the variety of RNA transcripts and proteins. However, AS of non-coding RNAs (ncRNAs) has not been paid enough attention to compared with that of protein-coding RNAs (mRNAs) for a long time. In fact, AS of ncRNAs, especially long noncoding RNAs (lncRNAs), also plays a significant regulatory role in the human disease. Recently, some bifunctional genes transcribed into both mRNA and lncRNA transcripts by AS have been observed. Here, we focus on the AS of lncRNAs and bifunctional genes producing lncRNA transcripts and propose a strategy for the future research of lncRNA AS.

Comprehensive characterisation of intronic mis-splicing mutations in human cancers

Previous studies studying mis-splicing mutations were based on exome data and thus our current knowledge is largely limited to exons and the canonical splice sites. To comprehensively characterise intronic mis-splicing mutations, we analysed 1134 pan-cancer whole genomes and transcriptomes together with 3022 normal control samples. The ratio-based splicing analysis resulted in 678 somatic intronic mutations, with 46% residing in deep introns. Among the 309 deep intronic single nucleotide variants, 245 altered core splicing codes, with 38% activating cryptic splice sites, 12% activating cryptic polypyrimidine tracts, and 36% and 12% disrupting authentic polypyrimidine tracts and branchpoints, respectively. All the intronic cryptic splice sites were created at pre-existing GT/AG dinucleotides or by GC-to-GT conversion. Notably, 85 deep intronic mutations indicated gain of splicing enhancers or loss of splicing silencers. We found that 64 tumour suppressors were affected by intronic mutations and blood cancers showed higher proportion of deep intronic mutations. In particular, a telomere maintenance gene, POT1, was recurrently mis-spliced by deep intronic mutations in blood cancers. We validated a pseudoexon activation involving a splicing silencer in POT1 by CRISPR/Cas9. Our results shed light on previously unappreciated mechanisms by which noncoding mutations acting on splicing codes in deep introns contribute to tumourigenesis.

Multiple Alternative Promoters and Alternative Splicing Enable Universal Transcription-Based Logic Computation in Mammalian Cells

Multi-input logic gene circuits can enable sophisticated control of cell function, yet large-scale synthetic circuitry in mammalian cells has relied on post-transcriptional regulation or recombinase-triggered state transitions. Large-scale transcriptional logic, on the other hand, has been challenging to implement. Inspired by a naturally found regulatory strategy of using multiple alternative promoters, followed by alternative splicing, we developed a scalable and compact platform for transcriptional OR logic using inputs to those promoters. The platform is extended to implement disjunctive normal form (DNF) computations capable of implementing arbitrary logic rules. Specifically, AND logic is implemented at individual promoters using synergistic transcriptional inputs, and NOT logic via microRNA inputs targeting unique exon sequences driven by those promoters. Together, these regulatory programs result in DNF-like logic control of output gene expression. The approach offers flexibility for building complex logic programs in mammalian cells.

Evaluation of both exonic and intronic variants for effects on RNA splicing allows for accurate assessment of the effectiveness of precision therapies

Elucidating the functional consequence of molecular defects underlying genetic diseases enables appropriate design of therapeutic options. Treatment of cystic fibrosis (CF) is an exemplar of this paradigm as the development of CFTR modulator therapies has allowed for targeted and effective treatment of individuals harboring specific genetic variants. However, the mechanism of these drugs limits effectiveness to particular classes of variants that allow production of CFTR protein. Thus, assessment of the molecular mechanism of individual variants is imperative for proper assignment of these precision therapies. This is particularly important when considering variants that affect pre-mRNA splicing, thus limiting success of the existing protein-targeted therapies. Variants affecting splicing can occur throughout exons and introns and the complexity of the process of splicing lends itself to a variety of outcomes, both at the RNA and protein levels, further complicating assessment of disease liability and modulator response. To investigate the scope of this challenge, we evaluated splicing and downstream effects of 52 naturally occurring CFTR variants (exonic = 15, intronic = 37). Expression of constructs containing select CFTR intronic sequences and complete CFTR exonic sequences in cell line models allowed for assessment of RNA and protein-level effects on an allele by allele basis. Characterization of primary nasal epithelial cells obtained from individuals harboring splice variants corroborated in vitro data. Notably, we identified exonic variants that result in complete missplicing and thus a lack of modulator response (e.g. c.2908G>A, c.523A>G), as well as intronic variants that respond to modulators due to the presence of residual normally spliced transcript (e.g. c.4242+2T>C, c.3717+40A>G). Overall, our data reveals diverse molecular outcomes amongst both exonic and intronic variants emphasizing the need to delineate RNA, protein, and functional effects of each variant in order to accurately assign precision therapies.

Genetic variants that impact pre-mRNA splicing are a common cause of genetic disease and have varying downstream molecular consequences. As a result, precision therapies that function at the protein level are not always effective for these variants and thus careful assessment is necessary. Here we evaluate RNA-level effects of 52 variants in the cystic fibrosis transmembrane conductance regulator (CFTR) gene and show that study of splicing and its consequences allows for more accurate assignment of precision therapies.

Enigmatic rhodopsin mutation creates an exceptionally strong splice acceptor site

The c.620 T > G mutation in rhodopsin found in the first mapped autosomal dominant retinitis pigmentosa (adRP) locus is associated with severe, early-onset RP. Intriguingly, another mutation affecting the same nucleotide (c.620 T > A) is related to a mild, late-onset RP. Assuming that both mutations are missense mutations (Met207Arg and Met207Lys) hampering the ligand-binding pocket, previous work addressed how they might differentially impair rhodopsin function. Here, we investigated the impact of both mutations at the mRNA and protein level in HEK293 cells and in the mouse retina. We show that, in contrast to c.620 T > A, c.620 T > G is a splicing mutation, which generates an exceptionally strong splice acceptor site (SAS) resulting in a 90 bp in-frame deletion and protein mislocalization in vitro and in vivo. Moreover, we identified the core element underlying the c.620 T > G SAS strength. Finally, we demonstrate that the c.620 T > G SAS is very flexible in branch point choice, which might explain its remarkable performance. Based on these results, we suggest that (i) point mutations should be routinely tested for mRNA splicing to avoid dispensable analysis of mutations on protein level, which do not naturally exist. (ii) Puzzling disease courses of mutations in other genes might also correlate with their effects on mRNA splicing. (iii) Flexibility in branch point choice might be another factor influencing the SAS strength. (iv) The core splice element identified in this study could be useful for biotechnological applications requiring effective SAS.

Splicing Characteristics of Dystrophin Pseudoexons and Identification of a Novel Pathogenic Intronic Variant in the DMD Gene

Pseudoexon (PE) inclusion has been implicated in various dystrophinopathies; however, its splicing characteristics have not been fully investigated. This study aims to analyze the splicing characteristics of dystrophin PEs and compare them with those of dystrophin canonical exons (CEs). Forty-two reported dystrophin PEs were divided into a splice site (ss) group and a splicing regulatory element (SRE) group. Five dystrophin PEs with characteristics of poison exons were identified and categorized as the possible poison exon group. The comparative analysis of each essential splicing signal among different groups of dystrophin PEs and dystrophin CEs revealed that the possible poison exon group had a stronger 3′ ss compared to any other group. As for auxiliary SREs, different groups of dystrophin PEs were found to have a smaller density of diverse types of exonic splicing enhancers and a higher density of several types of exonic splicing silencers compared to dystrophin CEs. In addition, the possible poison exon group had a smaller density of 3′ ss intronic splicing silencers compared to dystrophin CEs. To our knowledge, our findings indicate for the first time that poison exons might exist in DMD (the dystrophin gene) and present with different splicing characteristics than other dystrophin PEs and CEs.

A polyclonal allelic expression assay for detecting regulatory effects of transcript variants

We present an assay to experimentally test the regulatory effects of genetic variants within transcripts using CRISPR/Cas9 followed by targeted sequencing. We applied the assay to 32 premature stop-gained variants across the genome and in two Mendelian disease genes, 33 putative causal variants of eQTLs, and 62 control variants in HEK293T cells, replicating a subset of variants in HeLa cells. We detected significant effects in the expected direction (in 60% of variants), demonstrating the ability of the assay to capture regulatory effects of eQTL variants and nonsense-mediated decay triggered by premature stop-gained variants. The results suggest a utility for validating transcript-level effects of genetic variants.

Abnormal PTBP1 Expression Sustains the Disease Progression of Multiple Myeloma

Multiple myeloma (MM) is a hematopoietic malignancy characterized by heterogeneity, which corresponds to alternative splicing (AS) profiles and disadjust gene expression. Bioinformatics analysis of AS factors possibly related to MM progression identified the polypyrimidine tract binding protein (PTBP1) as candidate. The purpose of this study was to confirm the incidence and prognostic value of PTBP1 in MM patients. Several cohorts of 2971 patients presenting newly diagnosed and relapsed MM were enrolled. Correlations between PTBP1 expression and clinicopathological characteristics, proliferative activity, and response to therapy of myeloma cells were analyzed. Moreover, the effect of PTBP1 on the AS pattern of specific aerobic glycolysis-related genes was explored in MM patients. Clinically, PTBP1 expression was present at all stages; it increased with disease progression and poor prognosis, which was even stronger elevated in patients with high tumor burden and drug resistance. Mechanistically, PTBP1 modulated AS of PKM2 and aerobic glycolysis-related genes in MM patients, which play synergistic or additive effects in clinical outcome. PTBP1 may be a novel marker for prognostic prediction and a promising therapeutic target for the development of anti-MM treatments.

Regulations on Messenger RNA: Wires and Nodes

Somatic cells of an organism virtually share the same DNA but it is the timely expression of specific genes that determine their phenotype and cellular identity. A series of complex molecular machinery allows for the regulated process of RNA transcription, splicing, and translation. In addition, microRNAs and specialized RNA binding proteins can trigger the degradation of mRNAs. Long non-coding RNAs can also regulate mRNA fate in multiple ways. In this chapter, we reviewed the RNA processing mechanisms directly regulating immune checkpoint genes. We also cover RNA-based therapeutic strategies aiming at restoring immunity by targeting immune checkpoint genes.

Starvation-induced suppression of DAZAP1 by miR-10b integrates splicing control into TSC2-regulated oncogenic autophagy in esophageal squamous cell carcinoma

Esophageal squamous cell carcinoma (ESCC) accounts for about 90% of all incident esophageal cancers, with a 5-year survival rate of < 20%. Autophagy is of particular importance in cancers; however, the detailed regulatory mechanisms of oncogenic autophagy in ESCC have not been fully elucidated. In the present study, we address how splicing control of TSC2 is involved in mTOR-regulated oncogenic autophagy. Methods: Alternative splicing events controlled by DAZAP1 in ESCC cells were identified via RNAseq. Differential phosphorylation of short or long TSC2 splicing variants by AKT and their impacts on mTOR signaling were also examined. Results: We found that starvation-induced miR-10b could enhance autophagy via silencing DAZAP1, a key regulator of pre-mRNA alternative splicing. Intriguingly, we observed a large number of significantly changed alternative splicing events, especially exon skipping, upon RNAi of DAZAP1. TSC2 was verified as one of the crucial target genes of DAZAP1. Silencing of DAZAP1 led to the exclusion of TSC2 exon 26 (from Leu947 to Arg988), producing a short TSC2 isoform. The short TSC2 isoform cannot be phosphorylated at Ser981 by AKT, which resulted in continuous activation of TSC2 in ESCC. The active TSC2 inhibited mTOR via RHEB, leading to continually stimulated oncogenic autophagy of ESCC cells. Conclusions: Our data revealed an important physiological function of tumor suppressor DAZAP1 in autophagy regulation and highlighted the potential of controlling mRNA alternative splicing as an effective therapeutic application for cancers.

A Depletion of Stop Codons in lincRNA is Owing to Transfer of Selective Constraint from Coding Sequences

Although the constraints on a gene’s sequence are often assumed to reflect the functioning of that gene, here we propose transfer selection, a constraint operating on one class of genes transferred to another, mediated by shared binding factors. We show that such transfer can explain an otherwise paradoxical depletion of stop codons in long intergenic noncoding RNAs (lincRNAs). Serine/arginine-rich proteins direct the splicing machinery by binding exonic splice enhancers (ESEs) in immature mRNA. As coding exons cannot contain stop codons in one reading frame, stop codons should be rare within ESEs. We confirm that the stop codon density (SCD) in ESE motifs is low, even accounting for nucleotide biases. Given that serine/arginine-rich proteins binding ESEs also facilitate lincRNA splicing, a low SCD could transfer to lincRNAs. As predicted, multiexon lincRNA exons are depleted in stop codons, a result not explained by open reading frame (ORF) contamination. Consistent with transfer selection, stop codon depletion in lincRNAs is most acute in exonic regions with the highest ESE density, disappears when ESEs are masked, is consistent with stop codon usage skews in ESEs, and is diminished in both single-exon lincRNAs and introns. Owing to low SCD, the maximum lengths of pseudo-ORFs frequently exceed null expectations. This has implications for ORF annotation and the evolution of de novo protein-coding genes from lincRNAs. We conclude that not all constraints operating on genes need be explained by the functioning of the gene but may instead be transferred owing to shared binding factors.

Programmable RNA manipulation in living cells

RNAs are responsible for mediating genetic information flow within the cell. RNA splicing, modification, trafficking, translation, and stability are all controlled at the transcript level. However, biological tools to study and manipulate them in a programmable fashion are currently limited. In this review, we summarize recent advances regarding available RNA-targeting systems discovered so far, including CRISPR-based technologies-Cas9 and Cas13, and programmable RNA-binding proteins-PUF and PPR. These tools allow transcript-specific manipulation in gene expression.

Splicing in the pathogenesis, diagnosis and treatment of ciliopathies

Primary cilia are essential signalling organelles found on the apical surface of epithelial cells, where they coordinate chemosensation, mechanosensation and light sensation. Motile cilia play a central role in establishing fluid flow in the respiratory tract, reproductive tract, brain ventricles and ear. Genetic defects affecting the structure or function of cilia can lead to a broad range of developmental and degenerative diseases known as ciliopathies. Splicing contributes to the pathogenesis, diagnosis and treatment of ciliopathies. Tissue-specific alternative splicing contributes to the tissue-specific manifestation of ciliopathy phenotypes, for example the retinal-specific effects of some genetic defects, due to specific transcript expression in the highly specialised ciliated cells of the retina, the photoreceptor cells. Ciliopathies can arise both as a result of genetic variants in spliceosomal proteins, or as a result of variants affecting splicing of specific cilia genes. Here we discuss the opportunities and challenges in diagnosing ciliopathies using RNA sequence analysis and the potential for treating ciliopathies in a relatively mutation-neutral way by targeting splicing. This article is part of a Special Issue entitled: RNA structure and splicing regulation edited by Francisco Baralle, Ravindra Singh and Stefan Stamm.

Context matters: Regulation of splice donor usage

Elaborate research on splicing, starting in the late seventies, evolved from the discovery that 5' splice sites are recognized by their complementarity to U1 snRNA towards the realization that RNA duplex formation cannot be the sole basis for 5'ss selection. Rather, their recognition is highly influenced by a number of context factors including transcript architecture as well as splicing regulatory elements (SREs) in the splice site neighborhood. In particular, proximal binding of splicing regulatory proteins highly influences splicing outcome. The importance of SRE integrity especially becomes evident in the light of human pathogenic mutations where single nucleotide changes in SREs can severely affect the resulting transcripts. Bioinformatics tools nowadays greatly assist in the computational evaluation of 5'ss, their neighborhood and the impact of pathogenic mutations. Although predictions are already quite robust, computational evaluation of the splicing regulatory landscape still faces challenges to increase future reliability. This article is part of a Special Issue entitled: RNA structure and splicing regulation edited by Francisco Baralle, Ravindra Singh and Stefan Stamm.

RNA-Binding Motif 4 (RBM4) Suppresses Tumor Growth and Metastasis in Human Gastric Cancer

BACKGROUND Dysregulation of the splicing activator, RNA-binding motif 4 (RBM4), has recently been reported to be involved in the progression of several cancers. However, the mechanisms that underpin the activity of RBM4 in gastric cancer (GC) remain unknown. The purpose of our study was to explore how RBM4 affects the biological behavior of GC through in vivo and in vitro experiments. MATERIAL AND METHODS Western blot and flow cytometry analyses were used to investigate the RBM4 protein levels in normal gastric epithelial cells and 5 types of GC cells. Cell Counting Kit-8 assay, flow cytometry analysis, wound-healing, and migration and invasion assays were evaluated in vitro in BGC823 and MGC803 GC cells. A xenograft tumor model was used to assess whether RBM4 inhibits GC growth in vivo. Mitogen-activated protein kinase (MAPK) protein levels were determined using western blot analyses. RESULTS Our study revealed that RBM4 protein was downregulated in GC cells. Re-expression of RBM4 inhibited the proliferation, migration, and invasion of GC cells, while promoting apoptosis. Thus, the overexpression of RBM4 can inhibit tumor growth in GC mouse models. We also report that RBM4 was involved in the activation of MAPK-dependent signaling pathways in human GC. CONCLUSIONS It is hoped that these findings will improve our understanding of GC pathogenesis while also helping us to explore the feasibility of RBM4-targeted therapy for GC treatment.

Alternative splicing in lung cancer

Lung cancer has the highest mortality rate of all cancers worldwide. Lung cancer is a very heterogeneous disease that is often diagnosed at later stages which have a poor prognosis. Aberrant alternative splicing patterns found in lung cancer contribute to important cell functions. These include changes in splicing for the BCL2L1, MDM2, MDM4, NUMB and MET genes during lung tumourigenesis, to affect pathways involved in apoptosis, cell proliferation and cellular cohesion. Global analyses of RNASeq datasets suggest there may be many more potentially influential aberrant splicing events that need to be investigated in lung cancer. Changes in expression of the splicing factors that regulate alternative splicing events have also been identified in lung cancer. Of these, changes in expression of QKI, RBM4, RBM5, RBM6, RBM10 and SRSF1 proteins regulate many of the most frequently referenced aberrant splicing events in lung cancer. The expanding list of genes known to be aberrantly spliced in lung cancer along with the altered expression of splicing factors that regulate them are providing new clues as to how lung cancer develops, and how these events can be exploited for better treatment. This article is part of a Special Issue entitled: RNA structure and splicing regulation edited by Francisco Baralle, Ravindra Singh and Stefan Stamm.

Splicing Factors Have an Essential Role in Prostate Cancer Progression and Androgen Receptor Signaling

Although inhibition of the androgen–androgen receptor (AR) axis effectively represses the growth of prostate cancer, most of all cases eventually become castration-resistant prostate cancers (CRPCs). Enhancement of the expression of AR and its variants along with the downstream signals is important for disease progression. AR-V7, a constitutive active form of AR, is generated as a result of RNA splicing. RNA splicing creates multiple transcript variants from one pre-messenger RNA (mRNA) by removing introns/exons to allow mRNA translation. The molecular mechanisms leading to marked increases of AR and generation of AR-V7 have been unclear. However, recent papers highlighted the roles of RNA splicing factors which promote AR expression and production of variants. Notably, a broad range of splicing components were aberrantly regulated in CRPC tissues. Interestingly, expression of various spliceosome genes is enhanced by RNA-binding protein splicing factor proline- and glutamine-rich (PSF/SFPQ), leading to changes in the expression of AR transcript variants. Moreover, inhibition of several splicing factors repressed tumor growth in vivo. Altered expression of splicing factors is correlated to biochemical recurrence in prostate cancer patients. Thus, these findings suggest that splicing factors would be a potential therapeutic target. This review focuses on the emerging roles of splicing factors in prostate cancer progression and AR signaling.

RNA-Seq approach for accurate characterization of splicing efficiency of yeast introns

Introns in different genes, or even different introns within the same gene, often have different splice sites and differ in splicing efficiency (SE). One expects mass-transcribed genes to have introns with higher SE than weakly transcribed genes. However, such a simple expectation cannot be tested directly because variable SE for these genes is often not measured. Mechanistically, SE should depend on signal strength at key splice sites (SS) such as 5'SS, 3'SS and branchpoint site (BPS), i.e., SE = F(5'SS, 3'SS, BPS). However, without SE, we again cannot model how these splice sites contribute to SE. Here I present an RNA-Seq approach to quantify SE for each of the 304 introns in yeast (Saccharomyces cerevisiae) genes, including 24 in the 5'UTR, by measuring 1) number of reads mapped to exon-exon junctions (N_EE) as a proxy for the abundance of spliced form, and 2) number of reads mapped to exon-intron junction (N_EI5 and N_EI3 at 5' and 3' ends of intron) as a proxy for the abundance of unspliced form. The total mRNA is N_Total = N_EE + p * N_EI5 + (1-p) * N_EI3, with the simplest p = 0.5 but statistical methods were presented to estimate p from data. An estimated p is needed because N_EI5 is expected to be smaller than N_EI3 due to 1) step 1 splicing occurs before step 2 so EI5 is broken before EI3, 2) enrichment of poly(A) mRNA by oligo-dT, and 3) 5' degradation. SE is defined as the proportion (N_EE/N_Total). Application of the method shows that ribosomal protein messages are efficiently and mostly cotranscriptionally spliced. Yeast genes with long introns are also spliced efficiently. HAC1/YFL031W is poorly spliced partly because its splicing involves a nonspliceosome mechanism and partly because Ire1p, which participate in splicing HAC1, is hardly expressed. Many putative yeast genes have low SE, and some splice sites are incorrectly annotated.

Splicing and cancer: Challenges and opportunities

Cancer arises from alterations in several metabolic processes affecting proliferation, growth, replication and death of cells. A fundamental challenge in the study of cancer biology is to uncover molecular mechanisms that lead to malignant cellular transformation. Recent genomic analyses revealed that many molecular alterations observed in cancers come from modifications in the splicing process, including mutations in pre-mRNA regulatory sequences, mutations in spliceosome components, and altered ratio of specific splicing regulators. While alterations in splice site preferences might generate alternative isoforms enabling different biological functions, these might also be responsible for nonfunctional isoforms that can eventually cause dysregulation in cellular processes. Molecular characteristics of regulatory sequences and proteins might also be important prognostic tools revealing a cancer-specific splicing pattern and linking splicing control to cancer development. The connection between cancer biology and splicing regulation is of primary importance to understand the mechanisms leading to disease and also to improve development of therapeutic approaches. Splicing modulation is being explored in new anti-cancer therapies and further investigation of targeted splicing factors is critical for the success of these strategies. This article is categorized under: RNA Processing > Splicing Mechanisms RNA-Based Catalysis > RNA Catalysis in Splicing and Translation RNA Processing > Splicing Regulation/Alternative Splicing RNA in Disease and Development > RNA in Disease.

Splicing of long non-coding RNAs primarily depends on polypyrimidine tract and 5' splice-site sequences due to weak interactions with SR proteins

Many nascent long non-coding RNAs (lncRNAs) undergo the same maturation steps as pre-mRNAs of protein-coding genes (PCGs), but they are often poorly spliced. To identify the underlying mechanisms for this phenomenon, we searched for putative splicing inhibitory sequences using the ncRNA-a2 as a model. Genome-wide analyses of intergenic lncRNAs (lincRNAs) revealed that lincRNA splicing efficiency positively correlates with 5'ss strength while no such correlation was identified for PCGs. In addition, efficiently spliced lincRNAs have higher thymidine content in the polypyrimidine tract (PPT) compared to efficiently spliced PCGs. Using model lincRNAs, we provide experimental evidence that strengthening the 5'ss and increasing the T content in PPT significantly enhances lincRNA splicing. We further showed that lincRNA exons contain less putative binding sites for SR proteins. To map binding of SR proteins to lincRNAs, we performed iCLIP with SRSF2, SRSF5 and SRSF6 and analyzed eCLIP data for SRSF1, SRSF7 and SRSF9. All examined SR proteins bind lincRNA exons to a much lower extent than expression-matched PCGs. We propose that lincRNAs lack the cooperative interaction network that enhances splicing, which renders their splicing outcome more dependent on the optimality of splice sites.

Modeling and Predicting the Activities of Trans-Acting Splicing Factors with Machine Learning

Alternative splicing (AS) is generally regulated by trans-splicing factors that specifically bind to cis-elements in pre-mRNAs. The human genome encodes ∼1,500 RNA binding proteins (RBPs) that potentially regulate AS, yet their functions remain largely unknown. To explore their potential activities, we fused the putative functional domains of RBPs to a sequence-specific RNA-binding domain and systemically analyzed how these engineered factors affect splicing. We discovered that ∼80% of low-complexity domains in endogenous RBPs displayed distinct context-dependent activities in regulating splicing, indicating that AS is under more extensive regulation than previously expected. We developed a machine learning approach to classify and predict the activities of RBPs based on their sequence compositions and further validated this model using endogenous RBPs and synthetic polypeptides. These results represent a systematic inspection, modeling, prediction, and validation of how RBP sequences affect their activities in controlling splicing, paving the way for de novo engineering of artificial splicing factors.

Sequence, Structure, and Context Preferences of Human RNA Binding Proteins

RNA binding proteins (RBPs) orchestrate the production, processing, and function of mRNAs. Here, we present the affinity landscapes of 78 human RBPs using an unbiased assay that determines the sequence, structure, and context preferences of these proteins in vitro by deep sequencing of bound RNAs. These data enable construction of "RNA maps" of RBP activity without requiring crosslinking-based assays. We found an unexpectedly low diversity of RNA motifs, implying frequent convergence of binding specificity toward a relatively small set of RNA motifs, many with low compositional complexity. Offsetting this trend, however, we observed extensive preferences for contextual features distinct from short linear RNA motifs, including spaced "bipartite" motifs, biased flanking nucleotide composition, and bias away from or toward RNA structure. Our results emphasize the importance of contextual features in RNA recognition, which likely enable targeting of distinct subsets of transcripts by different RBPs that recognize the same linear motif.