Regulation of constitutive and alternative mRNA splicing across the human transcriptome by PRPF8 is determined by 5' splice site strength

Pluripotent stem cell-derived models of retinal disease: Elucidating pathogenesis, evaluating novel treatments, and estimating toxicity

Blindness poses a growing global challenge, with approximately 26% of cases attributed to degenerative retinal diseases. While gene therapy, optogenetic tools, photosensitive switches, and retinal prostheses offer hope for vision restoration, these high-cost therapies will benefit few patients. Understanding retinal diseases is therefore key to advance effective treatments, requiring in vitro models replicating pathology and allowing quantitative assessments for drug discovery. Pluripotent stem cells (PSCs) provide a unique solution given their limitless supply and ability to differentiate into light-responsive retinal tissues encompassing all cell types. This review focuses on the history and current state of photoreceptor and retinal pigment epithelium (RPE) cell generation from PSCs. We explore the applications of this technology in disease modelling, experimental therapy testing, biomarker identification, and toxicity studies. We consider challenges in scalability, standardisation, and reproducibility, and stress the importance of incorporating vasculature and immune cells into retinal organoids. We advocate for high-throughput automation in data acquisition and analyses and underscore the value of advanced micro-physiological systems that fully capture the interactions between the neural retina, RPE, and choriocapillaris.

PRPF8-mediated dysregulation of hBrr2 helicase disrupts human spliceosome kinetics and 5´-splice-site selection causing tissue-specific defects

The carboxy-terminus of the spliceosomal protein PRPF8, which regulates the RNA helicase Brr2, is a hotspot for mutations causing retinitis pigmentosa-type 13, with unclear role in human splicing and tissue-specificity mechanism. We used patient induced pluripotent stem cells-derived cells, carrying the heterozygous PRPF8 c.6926 A > C (p.H2309P) mutation to demonstrate retinal-specific endophenotypes comprising photoreceptor loss, apical-basal polarity and ciliary defects. Comprehensive molecular, transcriptomic, and proteomic analyses revealed a role of the PRPF8/Brr2 regulation in 5'-splice site (5'SS) selection by spliceosomes, for which disruption impaired alternative splicing and weak/suboptimal 5'SS selection, and enhanced cryptic splicing, predominantly in ciliary and retinal-specific transcripts. Altered splicing efficiency, nuclear speckles organisation, and PRPF8 interaction with U6 snRNA, caused accumulation of active spliceosomes and poly(A)+ mRNAs in unique splicing clusters located at the nuclear periphery of photoreceptors. Collectively these elucidate the role of PRPF8/Brr2 regulatory mechanisms in splicing and the molecular basis of retinal disease, informing therapeutic approaches.

Critical Differential Expression Assessment for Individual Bulk RNA-Seq Projects

Finding the right balance of quality and quantity can be important, and it is essential that project quality does not drop below the level where important main conclusions are missed or misstated. We use knock-out and over-expression studies as a simplification to test recovery of a known causal gene in RNA-Seq cell line experiments. When single-end RNA-Seq reads are aligned with STAR and quantified with htseq-count, we found potential value in testing the use of the Generalized Linear Model (GLM) implementation of edgeR with robust dispersion estimation more frequently for either single-variate or multi-variate 2-group comparisons (with the possibility of defining criteria less stringent than |fold-change| > 1.5 and FDR < 0.05). When considering a limited number of patient sample comparisons with larger sample size, there might be some decreased variability between methods (except for DESeq1). However, at the same time, the ranking of the gene identified using immunohistochemistry (for ER/PR/HER2 in breast cancer samples from The Cancer Genome Atlas) showed as possible shift in performance compared to the cell line comparisons, potentially highlighting utility for standard statistical tests and/or limma-based analysis with larger sample sizes. If this continues to be true in additional studies and comparisons, then that could be consistent with the possibility that it may be important to allocate time for potential methods troubleshooting for genomics projects. Analysis of public data presented in this study does not consider all experimental designs, and presentation of downstream analysis is limited. So, any estimate from this simplification would be an underestimation of the true need for some methods testing for every project. Additionally, this set of independent cell line experiments has a limitation in being able to determine the frequency of missing a highly important gene if the problem is rare (such as 10% or lower). For example, if there was an assumption that only one method can be tested for "initial" analysis, then it is not completely clear to the extent that using edgeR-robust might perform better than DESeq2 in the cell line experiments. Importantly, we do not wish to cause undue concern, and we believe that it should often be possible to define a gene expression differential expression workflow that is suitable for some purposes for many samples. Nevertheless, at the same time, we provide a variety of measures that we believe emphasize the need to critically assess every individual project and maximize confidence in published results.

PRPF8 controls alternative splicing of PIRH2 to modulate the p53 pathway and survival of human ESCs

Human embryonic stem cells (hESCs) have great potential for developmental biology and regenerative medicine. However, extensive apoptosis often occurs when hESCs respond to various stresses or injuries. Understanding the molecular control and identifying new factors associated with hESC survival are fundamental to ensure the high quality of hESCs. In this study, we report that PRPF8, an RNA spliceosome component, is essential for hESC survival. PRPF8 knockdown (KD) induces p53 protein accumulation and activates the p53 pathway, leading to apoptosis in hESCs. Strikingly, silencing of p53 rescues PRPF8 KD-induced apoptosis, indicating that PRPF8 KD triggers hESC apoptosis through activating the p53 pathway. In search for the mechanism by which p53 pathway is activated by PRPF8 KD, we find that PRPF8 KD alters alternative splicing of many genes, including PIRH2 which encodes an E3 ubiquitin ligase of p53. PIRH2 has several isoforms such as PIRH2A, PIRH2B, and PIRH2C. Intriguingly, PRPF8 KD specifically increases the transcript level of the PIRH2B isoform, which lacks a RING domain and E3 ligase activity. Functionally, PIRH2B KD partially rescues the reduction in cell numbers and upregulation of P21 caused by PRPF8 KD in hESCs. The finding suggests that PRPF8 controls alternative splicing of PIRH2 to maintain the balance of p53 pathway activity and survival of hESCs. The PRPF8/PIRH2/p53 axis identified here provides new insights into how p53 pathway and hESC survival are precisely regulated at multiple layers, highlighting an important role of posttranscriptional machinery in supporting hESC survival.

Mouse nuclear RNAi-defective 2 promotes splicing of weak 5' splice sites

Removal of introns during pre-mRNA splicing, which is central to gene expression, initiates by base pairing of U1 snRNA with a 5' splice site (5'SS). In mammals, many introns contain weak 5'SSs that are not efficiently recognized by the canonical U1 snRNP, suggesting alternative mechanisms exist. Here, we develop a cross-linking immunoprecipitation coupled to a high-throughput sequencing method, BCLIP-seq, to identify NRDE2 (nuclear RNAi-defective 2), and CCDC174 (coiled-coil domain-containing 174) as novel RNA-binding proteins in mouse ES cells that associate with U1 snRNA and 5'SSs. Both proteins bind directly to U1 snRNA independently of canonical U1 snRNP-specific proteins, and they are required for the selection and effective processing of weak 5'SSs. Our results reveal that mammalian cells use noncanonical splicing factors bound directly to U1 snRNA to effectively select suboptimal 5'SS sequences in hundreds of genes, promoting proper splice site choice, and accurate pre-mRNA splicing.

Insplico: effective computational tool for studying splicing order of adjacent introns genome-wide with short and long RNA-seq reads

Although splicing occurs largely co-transcriptionally, the order by which introns are removed does not necessarily follow the order in which they are transcribed. Whereas several genomic features are known to influence whether or not an intron is spliced before its downstream neighbor, multiple questions related to adjacent introns' splicing order (AISO) remain unanswered. Here, we present Insplico, the first standalone software for quantifying AISO that works with both short and long read sequencing technologies. We first demonstrate its applicability and effectiveness using simulated reads and by recapitulating previously reported AISO patterns, which unveiled overlooked biases associated with long read sequencing. We next show that AISO around individual exons is remarkably constant across cell and tissue types and even upon major spliceosomal disruption, and it is evolutionarily conserved between human and mouse brains. We also establish a set of universal features associated with AISO patterns across various animal and plant species. Finally, we used Insplico to investigate AISO in the context of tissue-specific exons, particularly focusing on SRRM4-dependent microexons. We found that the majority of such microexons have non-canonical AISO, in which the downstream intron is spliced first, and we suggest two potential modes of SRRM4 regulation of microexons related to their AISO and various splicing-related features. Insplico is available on gitlab.com/aghr/insplico.

Genetic control of RNA editing in neurodegenerative disease

A-to-I RNA editing diversifies human transcriptome to confer its functional effects on the downstream genes or regulations, potentially involving in neurodegenerative pathogenesis. Its variabilities are attributed to multiple regulators, including the key factor of genetic variants. To comprehensively investigate the potentials of neurodegenerative disease-susceptibility variants from the view of A-to-I RNA editing, we analyzed matched genetic and transcriptomic data of 1596 samples across nine brain tissues and whole blood from two large consortiums, Accelerating Medicines Partnership-Alzheimer's Disease and Parkinson's Progression Markers Initiative. The large-scale and genome-wide identification of 95 198 RNA editing quantitative trait loci revealed the preferred genetic effects on adjacent editing events. Furthermore, to explore the underlying mechanisms of the genetic controls of A-to-I RNA editing, several top RNA-binding proteins were pointed out, such as EIF4A3, U2AF2, NOP58, FBL, NOP56 and DHX9, since their regulations on multiple RNA-editing events were probably interfered by these genetic variants. Moreover, these variants may also contribute to the variability of other molecular phenotypes associated with RNA editing, including the functions of 3 proteins, expressions of 277 genes and splicing of 449 events. All the analyses results shown in NeuroEdQTL (https://relab.xidian.edu.cn/NeuroEdQTL/) constituted a unique resource for the understanding of neurodegenerative pathogenesis from genotypes to phenotypes related to A-to-I RNA editing.

Predictive model of transcriptional elongation control identifies trans regulatory factors from chromatin signatures

Promoter-proximal Polymerase II (Pol II) pausing is a key rate-limiting step for gene expression. DNA and RNA-binding trans-acting factors regulating the extent of pausing have been identified. However, we lack a quantitative model of how interactions of these factors determine pausing, therefore the relative importance of implicated factors is unknown. Moreover, previously unknown regulators might exist. Here we address this gap with a machine learning model that accurately predicts the extent of promoter-proximal Pol II pausing from large-scale genome and transcriptome binding maps and gene annotation and sequence composition features. We demonstrate high accuracy and generalizability of the model by validation on an independent cell line which reveals the model's cell line agnostic character. Model interpretation in light of prior knowledge about molecular functions of regulatory factors confirms the interconnection of pausing with other RNA processing steps. Harnessing underlying feature contributions, we assess the relative importance of each factor, quantify their predictive effects and systematically identify previously unknown regulators of pausing. We additionally identify 16 previously unknown 7SK ncRNA interacting RNA-binding proteins predictive of pausing. Our work provides a framework to further our understanding of the regulation of the critical early steps in transcriptional elongation.

PRPF8 increases the aggressiveness of hepatocellular carcinoma by regulating FAK/AKT pathway via fibronectin 1 splicing

Hepatocellular carcinoma (HCC) pathogenesis is associated with alterations in splicing machinery components (spliceosome and splicing factors) and aberrant expression of oncogenic splice variants. We aimed to analyze the expression and potential role of the spliceosome component PRPF8 (pre-mRNA processing factor 8) in HCC. PRPF8 expression (mRNA/protein) was analyzed in a retrospective cohort of HCC patients (n = 172 HCC and nontumor tissues) and validated in two in silico cohorts (TCGA and CPTAC). PRPF8 expression was silenced in liver cancer cell lines and in xenograft tumors to understand the functional and mechanistic consequences. In silico RNAseq and CLIPseq data were also analyzed. Our results indicate that PRPF8 is overexpressed in HCC and associated with increased tumor aggressiveness (patient survival, etc.), expression of HCC-related splice variants, and modulation of critical genes implicated in cancer-related pathways. PRPF8 silencing ameliorated aggressiveness in vitro and decreased tumor growth in vivo. Analysis of in silico CLIPseq data in HepG2 cells demonstrated that PRPF8 binds preferentially to exons of protein-coding genes, and RNAseq analysis showed that PRPF8 silencing alters splicing events in multiple genes. Integrated and in vitro analyses revealed that PRPF8 silencing modulates fibronectin (FN1) splicing, promoting the exclusion of exon 40.2, which is paramount for binding to integrins. Consistent with this finding, PRPF8 silencing reduced FAK/AKT phosphorylation and blunted stress fiber formation. Indeed, HepG2 and Hep3B cells exhibited a lower invasive capacity in membranes treated with conditioned medium from PRPF8-silenced cells compared to medium from scramble-treated cells. This study demonstrates that PRPF8 is overexpressed and associated with aggressiveness in HCC and plays important roles in hepatocarcinogenesis by altering FN1 splicing, FAK/AKT activation and stress fiber formation.

Microtubule modification defects underlie cilium degeneration in cell models of retinitis pigmentosa associated with pre-mRNA splicing factor mutations

Retinitis pigmentosa (RP) is the most common cause of hereditary blindness, and may occur in isolation as a non-syndromic condition or alongside other features in a syndromic presentation. Biallelic or monoallelic mutations in one of eight genes encoding pre-mRNA splicing factors are associated with non-syndromic RP. The molecular mechanism of disease remains incompletely understood, limiting opportunities for targeted treatment. Here we use CRISPR and base edited PRPF6 and PRPF31 mutant cell lines, and publicly-available data from human PRPF31^+/− patient derived retinal organoids and PRPF31 siRNA-treated organotypic retinal cultures to confirm an enrichment of differential splicing of microtubule, centrosomal, cilium and DNA damage response pathway genes in these cells. We show that genes with microtubule/centrosome/centriole/cilium gene ontology terms are enriched for weak 3′ and 5′ splice sites, and that subtle defects in spliceosome activity predominantly affect efficiency of splicing of these exons. We suggest that the primary defect in PRPF6 or PRPF31 mutant cells is microtubule and centrosomal defects, leading to defects in cilium and mitotic spindle stability, with the latter leading to DNA damage, triggering differential splicing of DNA damage response genes to activate this pathway. Finally, we expand understanding of “splicing factor RP” by investigating the function of TTLL3, one of the most statistically differentially expressed genes in PRPF6 and PRPF31 mutant cells. We identify that TTLL3 is the only tubulin glycylase expressed in the human retina, essential for monoglycylation of microtubules of the cilium, including the retinal photoreceptor cilium, to prevent cilium degeneration and retinal degeneration. Our preliminary data suggest that rescue of tubulin glycylation through overexpression of TTLL3 is sufficient to rescue cilium number in PRPF6 and PRPF31 mutant cells, suggesting that this defect underlies the cellular defect and may represent a potential target for therapeutic intervention in this group of disorders.

The Role of Alternative Splicing in Cancer: Regulatory Mechanism, Therapeutic Strategy, and Bioinformatics Application

[Formula: see text] Alternative splicing (AS) can generate distinct transcripts and subsequent isoforms that play differential functions from the same pre-mRNA. Recently, increasing numbers of studies have emerged, unmasking the association between AS and cancer. In this review, we arranged AS events that are closely related to cancer progression and presented promising treatments based on AS for cancer therapy. Obtaining proliferative capacity, acquiring invasive properties, gaining angiogenic features, shifting metabolic ability, and getting immune escape inclination are all splicing events involved in biological processes. Spliceosome-targeted and antisense oligonucleotide technologies are two novel strategies that are hopeful in tumor therapy. In addition, bioinformatics applications based on AS were summarized for better prediction and elucidation of regulatory routines mingled in. Together, we aimed to provide a better understanding of complicated AS events associated with cancer biology and reveal AS a promising target of cancer treatment in the future.

Whole cell response to receptor stimulation involves many deep and distributed subcellular biochemical processes

Neurite outgrowth is an integrated whole cell response triggered by the cannabinoid-1 receptor. We sought to identify the many different biochemical pathways that contribute to this whole cell response. To understand underlying mechanisms, we identified subcellular processes (SCPs) composed of one or more biochemical pathways and their interactions required for this response. Differentially expressed genes and proteins were obtained from bulk transcriptomics and proteomic analysis of extracts from cells stimulated with a cannabinoid-1 receptor agonist. We used these differentially expressed genes and proteins to build networks of interacting SCPs by combining the expression data with prior pathway knowledge. From these SCP networks, we identified additional genes that when ablated, experimentally validated the SCP involvement in neurite outgrowth. Our experiments and informatics modeling allowed us to identify diverse SCPs such as those involved in pyrimidine metabolism, lipid biosynthesis, and mRNA splicing and stability, along with more predictable SCPs such as membrane vesicle transport and microtubule dynamics. We find that SCPs required for neurite outgrowth are widely distributed among many biochemical pathways required for constitutive cellular functions, several of which are termed ‘deep’, since they are distal to signaling pathways and the key SCPs directly involved in extension of the neurite. In contrast, ‘proximal’ SCPs are involved in microtubule growth and membrane vesicle transport dynamics required for neurite outgrowth. From these bioinformatics and dynamical models based on experimental data, we conclude that receptor-mediated regulation of subcellular functions for neurite outgrowth is both distributed, that is, involves many different biochemical pathways, and deep.

Intronic Polyadenylation in Acquired Cancer Drug Resistance Circumvented by Utilizing CRISPR/Cas9 with Homology-Directed Repair: The Tale of Human DNA Topoisomerase IIα

Simple Summary

DNA topoisomerase IIα (170 kDa, TOP2α/170) resolves nucleic acid topological entanglements by generating transient double-strand DNA breaks. TOP2α inhibitors/poisons stabilize TOP2α-DNA covalent complexes resulting in persistent DNA damage and are frequently utilized to treat a variety of cancers. Acquired resistance to these chemotherapeutic agents is often associated with decreased TOP2α/170 expression levels. Studies have demonstrated that a reduction in TOP2α/170 results from a type of alternative polyadenylation designated intronic polyadenylation (IPA). As a consequence of IPA, variant TOP2α mRNA transcripts have been characterized that have resulted in the translation of C-terminal truncated TOP2α isoforms with altered biological activities. In this paper, an example is discussed where circumvention of acquired TOP2α-mediated drug resistance was achieved by utilizing CRISPR/Cas9 specific gene editing of an exon/intron boundary through homology directed repair (HDR) to reduce TOP2α IPA. These results illustrate the therapeutic potential of CRISPR/Cas9/HDR to impact drug resistance associated with aberrant IPA.

Abstract

Intronic polyadenylation (IPA) plays a critical role in malignant transformation, development, progression, and cancer chemoresistance by contributing to transcriptome/proteome alterations. DNA topoisomerase IIα (170 kDa, TOP2α/170) is an established clinical target for anticancer agents whose efficacy is compromised by drug resistance often associated with a reduction of nuclear TOP2α/170 levels. In leukemia cell lines with acquired resistance to TOP2α-targeted drugs and reduced TOP2α/170 expression, variant TOP2α mRNA transcripts have been reported due to IPA that resulted in the translation of C-terminal truncated isoforms with altered nuclear-cytoplasmic distribution or heterodimerization with wild-type TOP2α/170. This review provides an overview of the various mechanisms regulating pre-mRNA processing and alternative polyadenylation, as well as the utilization of CRISPR/Cas9 specific gene editing through homology directed repair (HDR) to decrease IPA when splice sites are intrinsically weak or potentially mutated. The specific case of TOP2α exon 19/intron 19 splice site editing is discussed in etoposide-resistant human leukemia K562 cells as a tractable strategy to circumvent acquired TOP2α-mediated drug resistance. This example supports the importance of aberrant IPA in acquired drug resistance to TOP2α-targeted drugs. In addition, these results demonstrate the therapeutic potential of CRISPR/Cas9/HDR to impact drug resistance associated with aberrant splicing/polyadenylation.

Snrpb is required in murine neural crest cells for proper splicing and craniofacial morphogenesis

Heterozygous mutations in SNRPB, an essential core component of the five small ribonucleoprotein particles of the spliceosome, are responsible for cerebrocostomandibular syndrome (CCMS). We show that Snrpb heterozygous mouse embryos arrest shortly after implantation. Additionally, heterozygous deletion of Snrpb in the developing brain and neural crest cells models craniofacial malformations found in CCMS, and results in death shortly after birth. RNAseq analysis of mutant heads prior to morphological defects revealed increased exon skipping and intron retention in association with increased 5′ splice site strength. We found increased exon skipping in negative regulators of the P53 pathway, along with increased levels of nuclear P53 and P53 target genes. However, removing Trp53 in Snrpb heterozygous mutant neural crest cells did not completely rescue craniofacial development. We also found a small but significant increase in exon skipping of several transcripts required for head and midface development, including Smad2 and Rere. Furthermore, mutant embryos exhibited ectopic or missing expression of Fgf8 and Shh, which are required to coordinate face and brain development. Thus, we propose that mis-splicing of transcripts that regulate P53 activity and craniofacial-specific genes contributes to craniofacial malformations.

This article has an associated First Person interview with the first author of the paper.

Summary: We report the first mouse model for cerebrocostomandibular syndrome, showing that mis-splicing of transcripts that regulate P53 activity and craniofacial-specific genes contributes to craniofacial malformations.

Use of CRISPR/Cas9 with homology-directed repair to silence the human topoisomerase IIα intron-19 5’ splice site: Generation of etoposide resistance in human leukemia K562 cells

DNA Topoisomerase IIα (TOP2α/170) is an enzyme essential for proliferating cells. For rapidly multiplying malignancies, this has made TOP2α/170 an important target for etoposide and other clinically active anticancer drugs. Efficacy of these agents is often limited by chemoresistance related to alterations in TOP2α/170 expression levels. Our laboratory recently demonstrated reduced levels of TOP2α/170 and overexpression of a C-terminal truncated 90-kDa isoform, TOP2α/90, due to intronic polyadenylation (IPA; within intron 19) in an acquired etoposide-resistant K562 clonal cell line, K/VP.5. We previously reported that this isoform heterodimerized with TOP2α/170 and was a determinant of acquired resistance to etoposide. Optimization of the weak TOP2α exon 19/intron 19 5′ splice site in drug-resistant K/VP.5 cells by gene-editing restored TOP2α/170 levels, diminished TOP2α/90 expression, and circumvented drug resistance. Conversely, in the present study, silencing of the exon 19/intron 19 5′ splice site in parental K562 cells by CRISPR/Cas9 with homology-directed repair (HDR), and thereby forcing intron 19 retention, was used to induce resistance by disrupting normal RNA processing (i.e., gene knockout), and to further evaluate the role of TOP2α/170 and TOP2α/90 isoforms as resistance determinants. Gene-edited clones were identified by quantitative polymerase chain reaction (qPCR) and verified by Sanger sequencing. TOP2α/170 mRNA/protein expression levels were attenuated in the TOP2α gene-edited clones which resulted in resistance to etoposide as assessed by reduced etoposide-induced DNA damage (γH2AX, Comet assays) and growth inhibition. RNA-seq and qPCR studies suggested that intron 19 retention leads to decreased TOP2α/170 expression by degradation of the TOP2α edited mRNA transcripts. Forced expression of TOP2α/90 in the gene-edited K562 cells further decreased etoposide-induced DNA damage in support of a dominant negative role for this truncated isoform. Together results support the important role of both TOP2α/170 and TOP2α/90 as determinants of sensitivity/resistance to TOP2α-targeting agents.

C9ORF78 partially localizes to centromeres and plays a role in chromosome segregation

C9ORF78 is a poorly characterized protein found in diverse eukaryotes. Previous work indicated overexpression of C9ORF78 in malignant tissues indicating a possible involvement in growth regulatory pathways. Additional studies in fission yeast and humans uncover a potential function in regulating the spliceosome. In studies of GFP-tagged C9ORF78 we observed a dramatic reduction in protein abundance in cells grown to confluence and/or deprived of serum growth factors. Serum stimulation induced synchronous re-expression of the protein in HeLa cells. This effect was also observed with the endogenous protein. Overexpressing either E2F1 or N-Myc resulted in elevated C9ORF78 expression potentially explaining the serum-dependent upregulation of the protein. Immunofluorescence analysis indicates that C9ORF78 localizes to nuclei in interphase but does not appear to concentrate in speckles as would be expected for a splicing protein. Surprisingly, a subpopulation of C9ORF78 co-localizes with ACA, Mad1 and Ndc80 in mitotic cells suggesting that this protein associates with kinetochores or centromeres. Levels of C9ORF78 at the centromere/kinetochore also increased upon activation of the mitotic checkpoint. Furthermore, knocking-down C9ORF78 caused mitotic defects. These studies uncover novel mitotic function and subcellular localization of C9ORF78.

Cwc27, associated with retinal degeneration, functions as a splicing factor in vivo

Previous in vitro studies indicate that CWC27 functions as a splicing factor in the Bact spliceosome complex, interacting with CWC22 to form a landing platform for eIF4A3, a core component of the exon junction complex. However, the function of CWC27 as a splicing factor has not been validated in any in vivo systems. CWC27 variants have been shown to cause autosomal recessive retinal degeneration, in both syndromic and non-syndromic forms. The Cwc27K338fs/K338fs mouse model was shown to have significant retinal dysfunction and degeneration by 6 months of age. In this report, we have taken advantage of the Cwc27K338fs/K338fs mouse model to show that Cwc27 is involved in splicing in vivo in the context of the retina. Bulk RNA and single cell RNA-sequencing of the mouse retina showed that there were gene expression and splicing pattern changes, including alternative splice site usage and intron retention. Positive staining for CHOP suggests that ER stress may be activated in response to the splicing pattern changes and is a likely contributor to the disease mechanism. Our results provide the first evidence that CWC27 functions as a splicing factor in an in vivo context. The splicing defects and gene expression changes observed in the Cwc27K338fs/K338fs mouse retina provide insight to the potential disease mechanisms, paving the way for targeted therapeutic development.

The SON RNA splicing factor is required for intracellular trafficking structures that promote centriole assembly and ciliogenesis

Control of centrosome assembly is critical for cell division, intracellular trafficking, and cilia. Regulation of centrosome number occurs through the precise duplication of centrioles that reside in centrosomes. Here we explored transcriptional control of centriole assembly and find that the RNA splicing factor SON is specifically required for completing procentriole assembly. Whole genome mRNA sequencing identified genes whose splicing and expression are affected by the reduction of SON, with an enrichment in genes involved in the microtubule (MT) cytoskeleton, centrosome, and centriolar satellites. SON is required for the proper splicing and expression of CEP131, which encodes a major centriolar satellite protein and is required to organize the trafficking and MT network around the centrosomes. This study highlights the importance of the distinct MT trafficking network that is intimately associated with nascent centrioles and is responsible for procentriole development and efficient ciliogenesis.

What's Wrong in a Jump? Prediction and Validation of Splice Site Variants

Alternative splicing (AS) is a crucial process to enhance gene expression driving organism development. Interestingly, more than 95% of human genes undergo AS, producing multiple protein isoforms from the same transcript. Any alteration (e.g., nucleotide substitutions, insertions, and deletions) involving consensus splicing regulatory sequences in a specific gene may result in the production of aberrant and not properly working proteins. In this review, we introduce the key steps of splicing mechanism and describe all different types of genomic variants affecting this process (splicing variants in acceptor/donor sites or branch point or polypyrimidine tract, exonic, and deep intronic changes). Then, we provide an updated approach to improve splice variants detection. First, we review the main computational tools, including the recent Machine Learning-based algorithms, for the prediction of splice site variants, in order to characterize how a genomic variant interferes with splicing process. Next, we report the experimental methods to validate the predictive analyses are defined, distinguishing between methods testing RNA (transcriptomics analysis) or proteins (proteomics experiments). For both prediction and validation steps, benefits and weaknesses of each tool/procedure are accurately reported, as well as suggestions on which approaches are more suitable in diagnostic rather than in clinical research.

BRCA2 deficiency reveals that oxidative stress impairs RNaseH1 function to cripple mitochondrial DNA maintenance

Oxidative stress is a ubiquitous cellular challenge implicated in aging, neurodegeneration, and cancer. By studying pathogenic mutations in the tumor suppressor BRCA2, we identify a general mechanism by which oxidative stress restricts mitochondrial (mt)DNA replication. BRCA2 inactivation induces R-loop accumulation in the mtDNA regulatory region and diminishes mtDNA replication initiation. In BRCA2-deficient cells, intracellular reactive oxygen species (ROS) are elevated, and ROS scavengers suppress the mtDNA defects. Conversely, wild-type cells exposed to oxidative stress by pharmacologic or genetic manipulation phenocopy these defects. Mechanistically, we find that 8-oxoguanine accumulation in mtDNA caused by oxidative stress suffices to impair recruitment of the mitochondrial enzyme RNaseH1 to sites of R-loop accrual, restricting mtDNA replication initiation. Thus, oxidative stress impairs RNaseH1 function to cripple mtDNA maintenance. Our findings highlight a molecular mechanism that links oxidative stress to mitochondrial dysfunction and is elicited by the inactivation of genes implicated in neurodegeneration and cancer.

•

BRCA2-deficient cells accumulate mtDNA R-loops due to oxidative stress

•

This stress creates 8-oxoguanine lesions impairing RNaseH1 recruitment to mtDNA

•

RNaseH1 impairment triggers R-loop formation and restricts mtDNA replication

•

Other sources of oxidative stress also cripple mtDNA maintenance via this mechanism

Pre-mRNA Processing Factors and Retinitis Pigmentosa: RNA Splicing and Beyond

Retinitis pigmentosa (RP) is the most common inherited retinal disease characterized by progressive degeneration of photoreceptors and/or retinal pigment epithelium that eventually results in blindness. Mutations in pre-mRNA processing factors (PRPF3, 4, 6, 8, 31, SNRNP200, and RP9) have been linked to 15–20% of autosomal dominant RP (adRP) cases. Current evidence indicates that PRPF mutations cause retinal specific global spliceosome dysregulation, leading to mis-splicing of numerous genes that are involved in a variety of retina-specific functions and/or general biological processes, including phototransduction, retinol metabolism, photoreceptor disk morphogenesis, retinal cell polarity, ciliogenesis, cytoskeleton and tight junction organization, waste disposal, inflammation, and apoptosis. Importantly, additional PRPF functions beyond RNA splicing have been documented recently, suggesting a more complex mechanism underlying PRPF-RPs driven disease pathogenesis. The current review focuses on the key RP-PRPF genes, depicting the current understanding of their roles in RNA splicing, impact of their mutations on retinal cell’s transcriptome and phenome, discussed in the context of model species including yeast, zebrafish, and mice. Importantly, information on PRPF functions beyond RNA splicing are discussed, aiming at a holistic investigation of PRPF-RP pathogenesis. Finally, work performed in human patient-specific lab models and developing gene and cell-based replacement therapies for the treatment of PRPF-RPs are thoroughly discussed to allow the reader to get a deeper understanding of the disease mechanisms, which we believe will facilitate the establishment of novel and better therapeutic strategies for PRPF-RP patients.

Mutant PRPF8 Causes Widespread Splicing Changes in Spliceosome Components in Retinitis Pigmentosa Patient iPSC-Derived RPE Cells

Retinitis pigmentosa (RP) is a rare, progressive disease that affects photoreceptors and retinal pigment epithelial (RPE) cells with blindness as a final outcome. Despite high medical and social impact, there is currently no therapeutic options to slow down the progression of or cure the disease. The development of effective therapies was largely hindered by high genetic heterogeneity, inaccessible disease tissue, and unfaithful model organisms. The fact that components of ubiquitously expressed splicing factors lead to the retina-specific disease is an additional intriguing question. Herein, we sought to correlate the retinal cell-type-specific disease phenotype with the splicing profile shown by a patient with autosomal recessive RP, caused by a mutation in pre-mRNA splicing factor 8 (PRPF8). In order to get insight into the role of PRPF8 in homeostasis and disease, we capitalize on the ability to generate patient-specific RPE cells and reveal differentially expressed genes unique to RPE cells. We found that spliceosomal complex and ribosomal functions are crucial in determining cell-type specificity through differential expression and alternative splicing (AS) and that PRPF8 mutation causes global changes in splice site selection and exon inclusion that particularly affect genes involved in these cellular functions. This finding corroborates the hypothesis that retinal tissue identity is conferred by a specific splicing program and identifies retinal AS events as a framework toward the design of novel therapeutic opportunities.

Uncovering transcriptional dark matter via gene annotation independent single-cell RNA sequencing analysis

Conventional scRNA-seq expression analyses rely on the availability of a high quality genome annotation. Yet, as we show here with scRNA-seq experiments and analyses spanning human, mouse, chicken, mole rat, lemur and sea urchin, genome annotations are often incomplete, in particular for organisms that are not routinely studied. To overcome this hurdle, we created a scRNA-seq analysis routine that recovers biologically relevant transcriptional activity beyond the scope of the best available genome annotation by performing scRNA-seq analysis on any region in the genome for which transcriptional products are detected. Our tool generates a single-cell expression matrix for all transcriptionally active regions (TARs), performs single-cell TAR expression analysis to identify biologically significant TARs, and then annotates TARs using gene homology analysis. This procedure uses single-cell expression analyses as a filter to direct annotation efforts to biologically significant transcripts and thereby uncovers biology to which scRNA-seq would otherwise be in the dark.

Conventional single-cell RNA sequencing analysis rely on genome annotations that may be incomplete or inaccurate especially for understudied organisms. Here the authors present a bioinformatic tool that leverages single-cell data to uncover biologically relevant transcripts beyond the best available genome annotation.

Prpf31 is essential for the survival and differentiation of retinal progenitor cells by modulating alternative splicing

Dysfunction of splicing factors often result in abnormal cell differentiation and apoptosis, especially in neural tissues. Mutations in pre-mRNAs processing factor 31 (PRPF31) cause autosomal dominant retinitis pigmentosa, a progressive retinal degeneration disease. The transcriptome-wide splicing events specifically regulated by PRPF31 and their biological roles in the development and maintenance of retina are still unclear. Here, we showed that the differentiation and viability of retinal progenitor cells (RPCs) are severely perturbed in prpf31 knockout zebrafish when compared with other tissues at an early embryonic stage. At the cellular level, significant mitotic arrest and DNA damage were observed. These defects could be rescued by the wild-type human PRPF31 rather than the disease-associated mutants. Further bioinformatic analysis and experimental verification uncovered that Prpf31 deletion predominantly causes the skipping of exons with a weak 5′ splicing site. Moreover, genes necessary for DNA repair and mitotic progression are most enriched among the differentially spliced events, which may explain the cellular and tissular defects in prpf31 mutant retinas. This is the first time that Prpf31 is demonstrated to be essential for the survival and differentiation of RPCs during retinal neurogenesis by specifically modulating the alternative splicing of genes involved in DNA repair and mitosis.

Mutation in Eftud2 causes craniofacial defects in mice via mis-splicing of Mdm2 and increased P53

EFTUD2 is mutated in patients with mandibulofacial dysostosis with microcephaly (MFDM). We generated a mutant mouse line with conditional mutation in Eftud2 and used Wnt1-Cre2 to delete it in neural crest cells. Homozygous deletion of Eftud2 causes brain and craniofacial malformations, affecting the same precursors as in MFDM patients. RNAseq analysis of embryonic heads revealed a significant increase in exon skipping and increased levels of an alternatively spliced Mdm2 transcript lacking exon 3. Exon skipping in Mdm2 was also increased in O9-1 mouse neural crest cells after siRNA knock-down of Eftud2 and in MFDM patient cells. Moreover, we found increased nuclear P53, higher expression of P53-target genes and increased cell death. Finally, overactivation of the P53 pathway in Eftud2 knockdown cells was attenuated by overexpression of non-spliced Mdm2, and craniofacial development was improved when Eftud2-mutant embryos were treated with Pifithrin-α, an inhibitor of P53. Thus, our work indicates that the P53-pathway can be targeted to prevent craniofacial abnormalities and shows a previously unknown role for alternative splicing of Mdm2 in the etiology of MFDM.

Isoform-resolved correlation analysis between mRNA abundance regulation and protein level degradation

Profiling of biological relationships between different molecular layers dissects regulatory mechanisms that ultimately determine cellular function. To thoroughly assess the role of protein post‐translational turnover, we devised a strategy combining pulse stable isotope‐labeled amino acids in cells (pSILAC), data‐independent acquisition mass spectrometry (DIA‐MS), and a novel data analysis framework that resolves protein degradation rate on the level of mRNA alternative splicing isoforms and isoform groups. We demonstrated our approach by the genome‐wide correlation analysis between mRNA amounts and protein degradation across different strains of HeLa cells that harbor a high grade of gene dosage variation. The dataset revealed that specific biological processes, cellular organelles, spatial compartments of organelles, and individual protein isoforms of the same genes could have distinctive degradation rate. The protein degradation diversity thus dissects the corresponding buffering or concerting protein turnover control across cancer cell lines. The data further indicate that specific mRNA splicing events such as intron retention significantly impact the protein abundance levels. Our findings support the tight association between transcriptome variability and proteostasis and provide a methodological foundation for studying functional protein degradation.

Identification and analysis of splicing quantitative trait loci across multiple tissues in the human genome

Alternative splicing (AS) is a fundamental step in eukaryotic mRNA biogenesis. Here, we develop an efficient and reproducible pipeline for the discovery of genetic variants that affect AS (splicing QTLs, sQTLs). We use it to analyze the GTEx dataset, generating a comprehensive catalog of sQTLs in the human genome. Downstream analysis of this catalog provides insight into the mechanisms underlying splicing regulation. We report that a core set of sQTLs is shared across multiple tissues. sQTLs often target the global splicing pattern of genes, rather than individual splicing events. Many also affect the expression of the same or other genes, uncovering regulatory loci that act through different mechanisms. sQTLs tend to be located in post-transcriptionally spliced introns, which would function as hotspots for splicing regulation. While many variants affect splicing patterns by altering the sequence of splice sites, many more modify the binding sites of RNA-binding proteins. Genetic variants affecting splicing can have a stronger phenotypic impact than those affecting gene expression.

The Role of the U5 snRNP in Genetic Disorders and Cancer

Pre-mRNA splicing is performed by the spliceosome, a dynamic macromolecular complex consisting of five small uridine-rich ribonucleoprotein complexes (the U1, U2, U4, U5, and U6 snRNPs) and numerous auxiliary splicing factors. A plethora of human disorders are caused by genetic variants affecting the function and/or expression of splicing factors, including the core snRNP proteins. Variants in the genes encoding proteins of the U5 snRNP cause two distinct and tissue-specific human disease phenotypes – variants in PRPF6, PRPF8, and SNRP200 are associated with retinitis pigmentosa (RP), while variants in EFTUD2 and TXNL4A cause the craniofacial disorders mandibulofacial dysostosis Guion-Almeida type (MFDGA) and Burn-McKeown syndrome (BMKS), respectively. Furthermore, recurrent somatic mutations or changes in the expression levels of a number of U5 snRNP proteins (PRPF6, PRPF8, EFTUD2, DDX23, and SNRNP40) have been associated with human cancers. How and why variants in ubiquitously expressed spliceosome proteins required for pre-mRNA splicing in all human cells result in tissue-restricted disease phenotypes is not clear. Additionally, why variants in different, yet interacting, proteins making up the same core spliceosome snRNP result in completely distinct disease outcomes – RP, craniofacial defects or cancer – is unclear. In this review, we define the roles of different U5 snRNP proteins in RP, craniofacial disorders and cancer, including how disease-associated genetic variants affect pre-mRNA splicing and the proposed disease mechanisms. We then propose potential hypotheses for how U5 snRNP variants cause tissue specificity resulting in the restricted and distinct human disorders.

CRISPR/Cas9 Genome Editing of the Human Topoisomerase IIα Intron 19 5' Splice Site Circumvents Etoposide Resistance in Human Leukemia K562 Cells

An essential function of DNA topoisomerase IIα (TOP2α; 170 kDa, TOP2α/170) is to resolve DNA topologic entanglements during chromosome disjunction by introducing transient DNA double-stranded breaks. TOP2α/170 is an important target for DNA damage-stabilizing anticancer drugs, whose clinical efficacy is compromised by drug resistance often associated with decreased TOP2α/170 expression. We recently demonstrated that an etoposide-resistant K562 clonal subline, K/VP.5, with reduced levels of TOP2α/170, expresses high levels of a novel C-terminal truncated TOP2α isoform (90 kDa, TOP2α/90). TOP2α/90, the translation product of a TOP2α mRNA that retains a processed intron 19 (I19), heterodimerizes with TOP2α/170 and is a resistance determinant through a dominant-negative effect on drug activity. We hypothesized that genome editing to enhance I19 removal would provide a tractable strategy to circumvent acquired TOP2α-mediated drug resistance. To enhance I19 removal in K/VP.5 cells, CRISPR/Cas9 was used to make changes (GAG//GTAA AC →GAG//GTAA GT ) in the TOP2α gene's suboptimal exon 19/intron 19 5' splice site (E19/I19 5' SS). Gene-edited clones were identified by quantitative polymerase chain reaction and verified by sequencing. Characterization of a clone with all TOP2α alleles edited revealed improved I19 removal, decreased TOP2α/90 mRNA/protein, and increased TOP2α/170 mRNA/protein. Sensitivity to etoposide-induced DNA damage (γH2AX, Comet assays) and growth inhibition was restored to levels comparable to those in parental K562 cells. Together, the results indicate that our gene-editing strategy for optimizing the TOP2α E19/I19 5' SS in K/VP.5 cells circumvents resistance to etoposide and other TOP2α-targeted drugs. SIGNIFICANCE STATEMENT: Results presented here indicate that CRISPR/Cas9 gene editing of a suboptimal exon 19/intron 19 5' splice site in the DNA topoisomerase IIα (TOP2α) gene results in circumvention of acquired drug resistance to etoposide and other TOP2α-targeted drugs in a clonal K562 cell line by enhancing removal of intron 19 and thereby decreasing formation of a truncated TOP2α 90 kDa isoform and increasing expression of full-length TOP2α 170 kDa in these resistant cells. Results demonstrate the importance of RNA processing in acquired drug resistance to TOP2α-targeted drugs.

RNA in cancer

While the processing of mRNA is essential for gene expression, recent findings have highlighted that RNA processing is systematically altered in cancer. Mutations in RNA splicing factor genes and the shortening of 3' untranslated regions are widely observed. Moreover, evidence is accumulating that other types of RNAs, including circular RNAs, can contribute to tumorigenesis. In this Review, we highlight how altered processing or activity of coding and non-coding RNAs contributes to cancer. We introduce the regulation of gene expression by coding and non-coding RNA and discuss both established roles (microRNAs and long non-coding RNAs) and emerging roles (selective mRNA processing and circular RNAs) for RNAs, highlighting the potential mechanisms by which these RNA subtypes contribute to cancer. The widespread alteration of coding and non-coding RNA demonstrates that altered RNA biogenesis contributes to multiple hallmarks of cancer.

Alternative Splicing of Putative Stroke/Vascular Risk Factor Genes Expressed in Blood Following Ischemic Stroke Is Sexually Dimorphic and Cause-Specific

Genome-wide association studies have identified putative ischemic stroke risk genes, yet, their expression after stroke is unexplored in spite of growing interest in elucidating their specific role and identifying candidate genes for stroke treatment. Thus, we took an exploratory approach to investigate sexual dimorphism, alternative splicing, and etiology in putative risk gene expression in blood following cardioembolic, atherosclerotic large vessel disease and small vessel disease/lacunar causes of ischemic stroke in each sex compared to controls. Whole transcriptome arrays assessed 71 putative stroke/vascular risk factor genes for blood RNA expression at gene-, exon-, and alternative splicing-levels. Male (n = 122) and female (n = 123) stroke and control volunteers from three university medical centers were matched for race, age, vascular risk factors, and blood draw time since stroke onset. Exclusion criteria included: previous stroke, drug abuse, subarachnoid or intracerebral hemorrhage, hemorrhagic transformation, infection, dialysis, cancer, hematological abnormalities, thrombolytics, anticoagulants or immunosuppressants. Significant differential gene expression (fold change > |1.2|, p < 0.05, partial correlation > |0.4|) and alternative splicing (false discovery rate p < 0.3) were assessed. At gene level, few were differentially expressed: ALDH2, ALOX5AP, F13A1, and IMPA2 (males, all stroke); ITGB3 (females, cardioembolic); ADD1 (males, atherosclerotic); F13A1, IMPA2 (males, lacunar); and WNK1 (females, lacunar). GP1BA and ITGA2B were alternatively spliced in both sexes (all patients vs. controls). Six genes in males, five in females, were alternatively spliced in all stroke compared to controls. Alternative splicing and exon-level analyses associated many genes with specific etiology in either sex. Of 71 genes, 70 had differential exon-level expression in stroke patients compared to control subjects. Among stroke patients, 24 genes represented by differentially expressed exons were male-specific, six were common between sexes, and two were female-specific. In lacunar stroke, expression of 19 differentially expressed exons representing six genes (ADD1, NINJ2, PCSK9, PEMT, SMARCA4, WNK1) decreased in males and increased in females. Results demonstrate alternative splicing and sexually dimorphic expression of most putative risk genes in stroke patients' blood. Since expression was also often cause-specific, sex, and etiology are factors to consider in stroke treatment trials and genetic association studies as society trends toward more personalized medicine.

Spliceosomal factor mutations and mis-splicing in MDS

Somatic, heterozygous missense and nonsense mutations in at least seven proteins that function in the spliceosome are found at high frequency in MDS patients. These proteins act at various steps in the process of splicing by the spliceosome and lead to characteristic alterations in the alternative splicing of a subset of genes. Several studies have investigated the effects of these mutations and have attempted to identify a commonly affected gene or pathway. Here, we summarize what is known about the normal function of these proteins and how the mutations alter the splicing landscape of the genome. We also summarize the commonly mis-spliced gene targets and discuss the state of mechanistic unification that has been achieved. Finally, we discuss alternative mechanisms by which these mutations may lead to disease.

WBP11 is required for splicing the TUBGCP6 pre-mRNA to promote centriole duplication

Centriole duplication occurs once in each cell cycle to maintain centrosome number. A previous genome-wide screen revealed that depletion of 14 RNA splicing factors leads to a specific defect in centriole duplication, but the cause of this deficit remains unknown. Here, we identified an additional pre-mRNA splicing factor, WBP11, as a novel protein required for centriole duplication. Loss of WBP11 results in the retention of ∼200 introns, including multiple introns in TUBGCP6, a central component of the γ-TuRC. WBP11 depletion causes centriole duplication defects, in part by causing a rapid decline in the level of TUBGCP6. Several additional splicing factors that are required for centriole duplication interact with WBP11 and are required for TUBGCP6 expression. These findings provide insight into how the loss of a subset of splicing factors leads to a failure of centriole duplication. This may have clinical implications because mutations in some spliceosome proteins cause microcephaly and/or growth retardation, phenotypes that are strongly linked to centriole defects.

Epigenome-based splicing prediction using a recurrent neural network

Alternative RNA splicing provides an important means to expand metazoan transcriptome diversity. Contrary to what was accepted previously, splicing is now thought to predominantly take place during transcription. Motivated by emerging data showing the physical proximity of the spliceosome to Pol II, we surveyed the effect of epigenetic context on co-transcriptional splicing. In particular, we observed that splicing factors were not necessarily enriched at exon junctions and that most epigenetic signatures had a distinctly asymmetric profile around known splice sites. Given this, we tried to build an interpretable model that mimics the physical layout of splicing regulation where the chromatin context progressively changes as the Pol II moves along the guide DNA. We used a recurrent-neural-network architecture to predict the inclusion of a spliced exon based on adjacent epigenetic signals, and we showed that distinct spatio-temporal features of these signals were key determinants of model outcome, in addition to the actual nucleotide sequence of the guide DNA strand. After the model had been trained and tested (with >80% precision-recall curve metric), we explored the derived weights of the latent factors, finding they highlight the importance of the asymmetric time-direction of chromatin context during transcription.

Spliceosomopathies and neurocristopathies: Two sides of the same coin?

Mutations in core components of the spliceosome are responsible for a group of syndromes collectively known as spliceosomopathies. Patients exhibit microcephaly, micrognathia, malar hypoplasia, external ear anomalies, eye anomalies, psychomotor delay, intellectual disability, limb, and heart defects. Craniofacial malformations in these patients are predominantly found in neural crest cells-derived structures of the face and head. Mutations in eight genes SNRPB, RNU4ATAC, SF3B4, PUF60, EFTUD2, TXNL4, EIF4A3, and CWC27 are associated with craniofacial spliceosomopathies. In this review, we provide a brief description of the normal development of the head and the face and an overview of mutations identified in genes associated with craniofacial spliceosomopathies. We also describe a model to explain how and when these mutations are most likely to impact neural crest cells. We speculate that mutations in a subset of core splicing factors lead to disrupted splicing in neural crest cells because these cells have increased sensitivity to inefficient splicing. Hence, disruption in splicing likely activates a cellular stress response that includes increased skipping of regulatory exons in genes such as MDM2 and MDM4, key regulators of P53. This would result in P53-associated death of neural crest cells and consequently craniofacial malformations associated with spliceosomopathies.

Food-Derived Compounds Apigenin and Luteolin Modulate mRNA Splicing of Introns with Weak Splice Sites

Cancer cells often exhibit extreme sensitivity to splicing inhibitors. We identified food-derived flavonoids, apigenin and luteolin, as compounds that modulate mRNA splicing at the genome-wide level, followed by proliferation inhibition. They bind to mRNA splicing-related proteins to induce a widespread change of splicing patterns in treated cells. Their inhibitory activity on splicing is relatively moderate, and introns with weak splice sites tend to be sensitive to them. Such introns remain unspliced, and the resulting intron-containing mRNAs are retained in the nucleus, resulting in the nuclear accumulation of poly(A)⁺ RNAs in these flavonoid-treated cells. Tumorigenic cells are more susceptible to these flavonoids than nontumorigenic cells, both for the nuclear poly(A)⁺ RNA-accumulating phenotype and cell viability. This study illustrates the possible mechanism of these flavonoids to suppress tumor progression in vivo that were demonstrated by previous studies and provides the potential of daily intake of moderate splicing inhibitors to prevent cancer development.

•

Food-derived compounds, apigenin and luteolin, modulate mRNA splicing

•

The treatment of these flavonoids causes numerous alternative splicing events

•

Splicing of introns with weak splice sites tend to be inhibited by these flavonoids

•

Tumorigenic cells are more sensitive to these flavonoids than non-tumorigenic cells

The role of deubiquitinating enzymes in gastric cancer

The epigenetic regulation of gene expression (via DNA methylation, histone modification and microRNA interference) contributes to a variety of diseases, particularly cancer. Protein deubiquitination serves a key role in the mechanism underlying histone modification, and consequently influences tumor development and progression. Improved characterization of the role of ubiquitinating enzymes has led to the identification of numerous deubiquitinating enzymes (DUBs) with various functions. Gastric cancer (GC) is a highly prevalent cancer type that exhibits a high mortality rate. Latest analysis about cancer patient revealed that GC is sixth deadliest cancer type, which frequently occur in male (7.2%) than female (4.1%). Complex associations between DUBs and GC progression have been revealed in multiple studies; however, the molecular mechanism underpinning the metastasis and recurrence of GC is yet to be elucidated. Generally, DUBs were upregulated in gastric cancer. The relation of DUBs and tumor size, classification and staging was observed in GC. Besides, 5-yar survival rate of patients with GC is effeccted by expression level of DUBs. Among the highly expressed DUBs, specifically six DUBs namely UCHs, USPs, OTUs, MJDs, JAMMs and MCPIPs effect on this survival rate. Consequently, the association between GC and DUBs has received increasing attention in recent years. Therefore, in the present review, literature investigating the association between DUBs and GC pathophysiology was analyzed and critically appraised.

A missense mutation in SNRPE linked to non-syndromal microcephaly interferes with U snRNP assembly and pre-mRNA splicing

Malfunction of pre-mRNA processing factors are linked to several human diseases including cancer and neurodegeneration. Here we report the identification of a de novo heterozygous missense mutation in the SNRPE gene (c.65T>C (p.Phe22Ser)) in a patient with non-syndromal primary (congenital) microcephaly and intellectual disability. SNRPE encodes SmE, a basal component of pre-mRNA processing U snRNPs. We show that the microcephaly-linked SmE variant is unable to interact with the SMN complex and as a consequence fails to assemble into U snRNPs. This results in widespread mRNA splicing alterations in fibroblast cells derived from this patient. Similar alterations were observed in HEK293 cells upon SmE depletion that could be rescued by the expression of wild type but not mutant SmE. Importantly, the depletion of SmE in zebrafish causes aberrant mRNA splicing alterations and reduced brain size, reminiscent of the patient microcephaly phenotype. We identify the EMX2 mRNA, which encodes a protein required for proper brain development, as a major mis-spliced down stream target. Together, our study links defects in the SNRPE gene to microcephaly and suggests that alterations of cellular splicing of specific mRNAs such as EMX2 results in the neurological phenotype of the disease.

In higher eukaryotes, the protein coding genes are first transcribed as precursor mRNAs (pre-mRNAs) and further processed by the spliceosome to form the mature mRNA for translation. Malfunction of pre-mRNA processing factors are linked to several human diseases including cancer and neurodegeneration. Here we report the identification of a de novo heterozygous missense mutation in the SNRPE/SmE gene in a patient with non-syndromal primary (congenital) microcephaly and intellectual disability. The effect of identified de novo mutation on SNRPE/SmE was characterized in vitro. The zebrafish was used as in vivo model to further dissect the physiological consequence and pathomechanism. Finally, the EMX2 gene was identified as one of the major down stream target genes responsible for the phenotype. Our study links defects in the SNRPE/SmE gene to microcephaly and provides the new pathogenic mechanism for microcephaly.

PRPF31 reduction causes mis-splicing of the phototransduction genes in human organotypic retinal culture

PRPF31 is ubiquitously expressed splicing factor and has an essential role in the pre-mRNA splicing in all tissues. However, it is not clear how reduced expression of this general splicing factor leads to retinal restricted disease, retinitis pigmentosa (RP). In this study, we used RNA interference and RNA-sequencing to mimic the PRPF31 haploinsufficiency in human organotypic retinal cultures (HORCs). We examined the effects of PRPF31 deficiency on splicing by analyzing the differential exon usages (DEUs) and intron retentions of the retinal transcriptome. Our results revealed that the PRPF31 deficiency causes mis-splicing of genes involved in RNA processing (PRPF3, PRPF8, PRPF4, and PRPF19) and phototransduction (RHO, ROM1, FSCN2, GNAT2, and GNAT1) in the retina in the PRPF31 reduced samples. Mis-splicing of genes implicated in phototransduction was associated with photoreceptor degeneration observed in RP patients. Our data revealed that PRPF31 deficiency leads to the mis-splicing of a distinct subset of pre-mRNAs with a widespread effect on phototransduction and RNA processing.

Splicing events in the control of genome integrity: role of SLU7 and truncated SRSF3 proteins

Genome instability is related to disease development and carcinogenesis. DNA lesions are caused by genotoxic compounds but also by the dysregulation of fundamental processes like transcription, DNA replication and mitosis. Recent evidence indicates that impaired expression of RNA-binding proteins results in mitotic aberrations and the formation of transcription-associated RNA-DNA hybrids (R-loops), events strongly associated with DNA injury. We identify the splicing regulator SLU7 as a key mediator of genome stability. SLU7 knockdown results in R-loops formation, DNA damage, cell-cycle arrest and severe mitotic derangements with loss of sister chromatid cohesion (SCC). We define a molecular pathway through which SLU7 keeps in check the generation of truncated forms of the splicing factor SRSF3 (SRp20) (SRSF3-TR). Behaving as dominant negative, or by gain-of-function, SRSF3-TR impair the correct splicing and expression of the splicing regulator SRSF1 (ASF/SF2) and the crucial SCC protein sororin. This unique function of SLU7 was found in cancer cells of different tissue origin and also in the normal mouse liver, demonstrating a conserved and fundamental role of SLU7 in the preservation of genome integrity. Therefore, the dowregulation of SLU7 and the alterations of this pathway that we observe in the cirrhotic liver could be involved in the process of hepatocarcinogenesis.

A systems view of spliceosomal assembly and branchpoints with iCLIP

Studies of spliceosomal interactions are challenging due to their dynamic nature. Here we used spliceosome iCLIP, which immunoprecipitates SmB along with small nuclear ribonucleoprotein particles and auxiliary RNA binding proteins, to map spliceosome engagement with pre-messenger RNAs in human cell lines. This revealed seven peaks of spliceosomal crosslinking around branchpoints (BPs) and splice sites. We identified RNA binding proteins that crosslink to each peak, including known and candidate splicing factors. Moreover, we detected the use of over 40,000 BPs with strong sequence consensus and structural accessibility, which align well to nearby crosslinking peaks. We show how the position and strength of BPs affect the crosslinking patterns of spliceosomal factors, which bind more efficiently upstream of strong or proximally located BPs and downstream of weak or distally located BPs. These insights exemplify spliceosome iCLIP as a broadly applicable method for transcriptomic studies of splicing mechanisms.

RNA components of the spliceosome regulate tissue- and cancer-specific alternative splicing

Alternative splicing of pre-mRNAs plays a pivotal role during the establishment and maintenance of human cell types. Characterizing the trans-acting regulatory proteins that control alternative splicing has therefore been the focus of much research. Recent work has established that even core protein components of the spliceosome, which are required for splicing to proceed, can nonetheless contribute to splicing regulation by modulating splice site choice. We here show that the RNA components of the spliceosome likewise influence alternative splicing decisions. Although these small nuclear RNAs (snRNAs), termed U1, U2, U4, U5, and U6 snRNA, are present in equal stoichiometry within the spliceosome, we found that their relative levels vary by an order of magnitude during development, across tissues, and across cancer samples. Physiologically relevant perturbation of individual snRNAs drove widespread gene-specific differences in alternative splicing but not transcriptome-wide splicing failure. Genes that were particularly sensitive to variations in snRNA abundance in a breast cancer cell line model were likewise preferentially misspliced within a clinically diverse cohort of invasive breast ductal carcinomas. As aberrant mRNA splicing is prevalent in many cancers, we propose that a full understanding of such dysregulated pre-mRNA processing requires study of snRNAs, as well as protein splicing factors. Together, our data show that the RNA components of the spliceosome are not merely basal factors, as has long been assumed. Instead, these noncoding RNAs constitute a previously uncharacterized layer of regulation of alternative splicing, and contribute to the establishment of global splicing programs in both healthy and malignant cells.

Integration of CLIP experiments of RNA-binding proteins: a novel approach to predict context-dependent splicing factors from transcriptomic data

Background

Splicing is a genetic process that has important implications in several diseases including cancer. Deciphering the complex rules of splicing regulation is crucial to understand and treat splicing-related diseases. Splicing factors and other RNA-binding proteins (RBPs) play a key role in the regulation of splicing. The specific binding sites of an RBP can be measured using CLIP experiments. However, to unveil which RBPs regulate a condition, it is necessary to have a priori hypotheses, as a single CLIP experiment targets a single protein.

Results

In this work, we present a novel methodology to predict context-specific splicing factors from transcriptomic data. For this, we systematically collect, integrate and analyze more than 900 CLIP experiments stored in four CLIP databases: POSTAR2, CLIPdb, DoRiNA and StarBase. The analysis of these experiments shows the strong coherence between the binding sites of RBPs of similar families. Augmenting this information with expression changes, we are able to correctly predict the splicing factors that regulate splicing in two gold-standard experiments in which specific splicing factors are knocked-down.

Conclusions

The methodology presented in this study allows the prediction of active splicing factors in either cancer or any other condition by only using the information of transcript expression. This approach opens a wide range of possible studies to understand the splicing regulation of different conditions. A tutorial with the source code and databases is available at https://gitlab.com/fcarazo.m/sfprediction.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5900-1) contains supplementary material, which is available to authorized users.

Early splicing functions of fission yeast Prp16 and its unexpected requirement for gene Silencing is governed by intronic features

Prp16 is a DEAH box pre-mRNA splicing factor that triggers a key spliceosome conformational switch to facilitate second step splicing in Saccharomyces cerevisiae. However, Prp16 functions are largely unexplored in Schizosaccharomyces pombe, an attractive model with exon-intron architecture more relevant to several other eukaryotes. Here, we generated mis-sense alleles in SpPrp16 whose consequences on genome-wide splicing uncover its nearly global splicing role with only a small subset of unaffected introns. Prp16 dependent and independent intron categories displayed a striking difference in the strength of intronic 5' splice site (5'SS)-U6 snRNA and branch site (BS)-U2 snRNA interactions. Selective weakening of these interactions could convert a Prp16 dependent intron into an independent one. These results point to the role of SpPrp16 in destabilizing 5'SS-U6snRNA and BS-U2snRNA interactions which plausibly trigger structural alterations in the spliceosome to facilitate first step catalysis. Our data suggest that SpPrp16 interactions with early acting factors, its enzymatic activities and association with intronic elements collectively account for efficient and accurate first step catalysis. In addition to splicing derangements in the spprp16F528S mutant, we show that SpPrp16 influences cell cycle progression and centromeric heterochromatinization. We propose that strong 5'SS-U6 snRNA and BS-U2 snRNA complementarity of intron-like elements in non-coding RNAs which lead to complete splicing arrest and impaired Seb1 functions at the pericentromeric loci may cumulatively account for the heterochromatin defects in spprp16F528S cells. These findings suggest that the diverse Prp16 functions within a genome are likely governed by its intronic features that influence splice site-snRNA interaction strength.

Deregulated Splicing Is a Major Mechanism of RNA-Induced Toxicity in Huntington's Disease

Huntington's disease (HD) is caused by an expanded CAG repeat in the huntingtin (HTT) gene, translating into an elongated polyglutamine stretch. In addition to the neurotoxic mutant HTT protein, the mutant CAG repeat RNA can exert toxic functions by trapping RNA-binding proteins. While few examples of proteins that aberrantly bind to mutant HTT RNA and execute abnormal function in conjunction with the CAG repeat RNA have been described, an unbiased approach to identify the interactome of mutant HTT RNA is missing. Here, we describe the analysis of proteins that preferentially bind mutant HTT RNA using a mass spectrometry approach. We show that (I) the majority of proteins captured by mutant HTT RNA belong to the spliceosome pathway, (II) expression of mutant CAG repeat RNA induces mis-splicing in a HD cell model, (III) overexpression of one of the splice factors trapped by mutant HTT ameliorates the HD phenotype in a fly model and (VI) deregulated splicing occurs in human HD brain. Our data suggest that deregulated splicing is a prominent mechanism of RNA-induced toxicity in HD.

Alternative Splicing and Protein Diversity: Plants Versus Animals

Plants, unlike animals, exhibit a very high degree of plasticity in their growth and development and employ diverse strategies to cope with the variations during diurnal cycles and stressful conditions. Plants and animals, despite their remarkable morphological and physiological differences, share many basic cellular processes and regulatory mechanisms. Alternative splicing (AS) is one such gene regulatory mechanism that modulates gene expression in multiple ways. It is now well established that AS is prevalent in all multicellular eukaryotes including plants and humans. Emerging evidence indicates that in plants, as in animals, transcription and splicing are coupled. Here, we reviewed recent evidence in support of co-transcriptional splicing in plants and highlighted similarities and differences between plants and humans. An unsettled question in the field of AS is the extent to which splice isoforms contribute to protein diversity. To take a critical look at this question, we presented a comprehensive summary of the current status of research in this area in both plants and humans, discussed limitations with the currently used approaches and suggested improvements to current methods and alternative approaches. We end with a discussion on the potential role of epigenetic modifications and chromatin state in splicing memory in plants primed with stresses.

Spliceosome factors target timeless (tim) mRNA to control clock protein accumulation and circadian behavior in Drosophila

Transcription-translation feedback loops that comprise eukaryotic circadian clocks rely upon temporal delays that separate the phase of active transcription of clock genes, such as Drosophila period (per) and timeless (tim), from negative feedback by the two proteins. However, our understanding of the mechanisms involved is incomplete. Through an RNA interference screen, we found that pre-mRNA processing 4 (PRP4) kinase, a component of the U4/U5.U6 triple small nuclear ribonucleoprotein (tri-snRNP) spliceosome, and other tri-snRNP components regulate cycling of the molecular clock as well as rest:activity rhythms. Unbiased RNA-Sequencing uncovered an alternatively spliced intron in tim whose increased retention upon prp4 downregulation leads to decreased TIM levels. We demonstrate that the splicing of tim is rhythmic with a phase that parallels delayed accumulation of the protein in a 24 hr cycle. We propose that alternative splicing constitutes an important clock mechanism for delaying the daily accumulation of clock proteins, and thereby negative feedback by them.

Editorial note: This article has been through an editorial process in which the authors decide how to respond to the issues raised during peer review. The Reviewing Editor's assessment is that all the issues have been addressed (see decision letter).

Targeting the spliceosome for cutaneous squamous cell carcinoma therapy: a role for c-MYC and wild-type p53 in determining the degree of tumour selectivity

We show that suppression of the spliceosome has potential for the treatment of cutaneous squamous cell carcinoma (cSCC). The small-molecule inhibitors of the spliceosome at the most advanced stage of development target the splicing factor SF3B1/SF3b155. The majority of cSCC cell lines are more sensitive than normal skin cells to death induced by the SF3B1 inhibitor pladienolide B. Knockdown of SF3B1 and a range of other splicing factors with diverse roles in the spliceosome can also selectively kill cSCC cells. We demonstrate that endogenous c-MYC participates in conferring sensitivity to spliceosome inhibition. c-MYC expression is elevated in cSCC lines and its knockdown reduces alterations in mRNA splicing and attenuates cell death caused by interference with the spliceosome. In addition, this study provides further support for a key role of the p53 pathway in the response to spliceosome disruption. SF3B1 inhibition causes wild-type p53 upregulation associated with altered mRNA splicing and reduced protein expression of both principal p53 negative regulators MDMX/MDM4 and MDM2. We observed that wild-type p53 can promote pladienolide B-induced death in tumour cells. However, p53 is commonly inactivated by mutation in cSCCs and p53 participates in killing normal skin cells at high concentrations of pladienolide B. This may limit the therapeutic window of SF3B1 inhibitors for cSCC. We provide evidence that, while suppression of SF3B1 has promise for treating cSCCs with mutant p53, inhibitors which target the spliceosome through SF3B1-independent mechanisms could have greater cSCC selectivity as a consequence of reduced p53 upregulation in normal cells.

Neurodevelopmental Genetic Diseases Associated With Microdeletions and Microduplications of Chromosome 17p13.3

Chromosome 17p13.3 is a region of genomic instability that is linked to different rare neurodevelopmental genetic diseases, depending on whether a deletion or duplication of the region has occurred. Chromosome microdeletions within 17p13.3 can result in either isolated lissencephaly sequence (ILS) or Miller-Dieker syndrome (MDS). Both conditions are associated with a smooth cerebral cortex, or lissencephaly, which leads to developmental delay, intellectual disability, and seizures. However, patients with MDS have larger deletions than patients with ILS, resulting in additional symptoms such as poor muscle tone, congenital anomalies, abnormal spasticity, and craniofacial dysmorphisms. In contrast to microdeletions in 17p13.3, recent studies have attracted considerable attention to a condition known as a 17p13.3 microduplication syndrome. Depending on the genes involved in their microduplication, patients with 17p13.3 microduplication syndrome may be categorized into either class I or class II. Individuals in class I have microduplications of the YWHAE gene encoding 14-3-3ε, as well as other genes in the region. However, the PAFAH1B1 gene encoding LIS1 is never duplicated in these patients. Class I microduplications generally result in learning disabilities, autism, and developmental delays, among other disorders. Individuals in class II always have microduplications of the PAFAH1B1 gene, which may include YWHAE and other genetic microduplications. Class II microduplications generally result in smaller body size, developmental delays, microcephaly, and other brain malformations. Here, we review the phenotypes associated with copy number variations (CNVs) of chromosome 17p13.3 and detail their developmental connection to particular microdeletions or microduplications. We also focus on existing single and double knockout mouse models that have been used to study human phenotypes, since the highly limited number of patients makes a study of these conditions difficult in humans. These models are also crucial for the study of brain development at a mechanistic level since this cannot be accomplished in humans. Finally, we emphasize the usefulness of the CRISPR/Cas9 system and next generation sequencing in the study of neurodevelopmental diseases.