PUF60-activated exons uncover altered 3' splice-site selection by germline missense mutations in a single RRM

MAPP unravels frequent co-regulation of splicing and polyadenylation by RNA-binding proteins and their dysregulation in cancer

Maturation of eukaryotic pre-mRNAs via splicing and polyadenylation is modulated across cell types and conditions by a variety of RNA-binding proteins (RBPs). Although there exist over 1,500 RBPs in human cells, their binding motifs and functions still remain to be elucidated, especially in the complex environment of tissues and in the context of diseases. To overcome the lack of methods for the systematic and automated detection of sequence motif-guided pre-mRNA processing regulation from RNA sequencing (RNA-Seq) data we have developed MAPP (Motif Activity on Pre-mRNA Processing). Applying MAPP to RBP knock-down experiments reveals that many RBPs regulate both splicing and polyadenylation of nascent transcripts by acting on similar sequence motifs. MAPP not only infers these sequence motifs, but also unravels the position-dependent impact of the RBPs on pre-mRNA processing. Interestingly, all investigated RBPs that act on both splicing and 3' end processing exhibit a consistently repressive or activating effect on both processes, providing a first glimpse on the underlying mechanism. Applying MAPP to normal and malignant brain tissue samples unveils that the motifs bound by the PTBP1 and RBFOX RBPs coordinately drive the oncogenic splicing program active in glioblastomas demonstrating that MAPP paves the way for characterizing pre-mRNA processing regulators under physiological and pathological conditions.

RNA binding proteins in cardiovascular development and disease

Congenital heart disease (CHD) is the most common birth defect affecting>1.35 million newborn babies worldwide. CHD can lead to prenatal, neonatal, postnatal lethality or life-long cardiac complications. RNA binding protein (RBP) mutations or variants are emerging as contributors to CHDs. RBPs are wizards of gene regulation and are major contributors to mRNA and protein landscape. However, not much is known about RBPs in the developing heart and their contributions to CHD. In this chapter, we will discuss our current knowledge about specific RBPs implicated in CHDs. We are in an exciting era to study RBPs using the currently available and highly successful RNA-based therapies and methodologies. Understanding how RBPs shape the developing heart will unveil their contributions to CHD. Identifying their target RNAs in the embryonic heart will ultimately lead to RNA-based treatments for congenital heart disease.

Identification of a novel de novo PUF60 variant causing Verheij syndrome in a fetus

Verheij syndrome (VRJS) is a craniofacial spliceosomopathy with a wide phenotypic spectrum. Haploinsufficiency of the poly-uridine binding splicing factor 60 gene (PUF60) and its loss-of-function (LOF) variants are involved in VRJS. We evaluated a human fetus with congenital heart defects and preaxial polydactyly. Clinical data were obtained from the medical record. Whole-exome sequencing (WES) was used to explore the potential genetic etiology, and the detected variant verified using Sanger sequencing. Functional studies were performed to validate the pathogenic effects of the variant. Using trio-WES, we identified a novel PUF60 variant (NM_078480.2; c.1678 T > A, p.*560Argext*204) in the pedigree. Bioinformatic analyses revealed that the variant is potentially pathogenic, and functional studies indicated that it leads to degradation of the elongated protein and subsequently PUF60 LOF, producing some VRJS phenotypes. These findings confirmed the pathogenicity of the variant. This study implicates PUF60 LOF in the etiopathogenesis of VRJS. It not only expands the PUF60 variant spectrum, and also provides a basis for genetic counseling and the diagnosis of VRJS. Although trio-WES is a well-established approach for identifying the genetic etiology of rare multisystemic conditions, functional studies could aid in verifying the pathogenicity of novel variants.

Exonic splicing code and coordination of divalent metals in proteins

Exonic sequences contain both protein-coding and RNA splicing information but the interplay of the protein and splicing code is complex and poorly understood. Here, we have studied traditional and auxiliary splicing codes of human exons that encode residues coordinating two essential divalent metals at the opposite ends of the Irving-Williams series, a universal order of relative stabilities of metal-organic complexes. We show that exons encoding Zn2+-coordinating amino acids are supported much less by the auxiliary splicing motifs than exons coordinating Ca2+. The handicap of the former is compensated by stronger splice sites and uridine-richer polypyrimidine tracts, except for position -3 relative to 3' splice junctions. However, both Ca2+ and Zn2+ exons exhibit close-to-constitutive splicing in multiple tissues, consistent with their critical importance for metalloprotein function and a relatively small fraction of expendable, alternatively spliced exons. These results indicate that constraints imposed by metal coordination spheres on RNA splicing have been efficiently overcome by the plasticity of exon-intron architecture to ensure adequate metalloprotein expression.

Novel Genetic and Phenotypic Expansion in Ameliorated PUF60-Related Disorders

Heterozygous variants in the Poly(U) Binding Splicing Factor 60kDa gene (PUF60) have been associated with Verheij syndrome, which has the key features of coloboma, short stature, skeletal abnormalities, developmental delay, palatal abnormalities, and congenital heart and kidney defects. Here, we report five novel patients from unrelated families with PUF60-related disorders exhibiting novel genetic and clinical findings with three truncating variants, one splice-site variant with likely reduced protein expression, and one missense variant. Protein modeling of the patient's missense variant in the PUF60 AlphaFold structure revealed a loss of polar bonds to the surrounding residues. Neurodevelopmental disorders were present in all patients, with variability in speech, motor, cognitive, social-emotional and behavioral features. Novel phenotypic expansions included movement disorders as well as immunological findings with recurrent respiratory, urinary and ear infections, atopic diseases, and skin abnormalities. We discuss the role of PUF60 in immunity with and without infection based on recent organismic and cellular studies. As our five patients showed less-severe phenotypes than classical Verheij syndrome, particularly with the absence of key features such as coloboma or palatal abnormalities, we propose a reclassification as PUF60-related neurodevelopmental disorders with multi-system involvement. These findings will aid in the genetic counseling of patients and families.

Arginine reprograms metabolism in liver cancer via RBM39

Metabolic reprogramming is a hallmark of cancer. However, mechanisms underlying metabolic reprogramming and how altered metabolism in turn enhances tumorigenicity are poorly understood. Here, we report that arginine levels are elevated in murine and patient hepatocellular carcinoma (HCC), despite reduced expression of arginine synthesis genes. Tumor cells accumulate high levels of arginine due to increased uptake and reduced arginine-to-polyamine conversion. Importantly, the high levels of arginine promote tumor formation via further metabolic reprogramming, including changes in glucose, amino acid, nucleotide, and fatty acid metabolism. Mechanistically, arginine binds RNA-binding motif protein 39 (RBM39) to control expression of metabolic genes. RBM39-mediated upregulation of asparagine synthesis leads to enhanced arginine uptake, creating a positive feedback loop to sustain high arginine levels and oncogenic metabolism. Thus, arginine is a second messenger-like molecule that reprograms metabolism to promote tumor growth.

PUF60 promotes cell cycle and lung cancer progression by regulating alternative splicing of CDC25C

Alternative splicing (AS) has been implicated in cell cycle regulation and cancer, but the underlying mechanisms are poorly understood. The poly(U)-binding splicing factor 60 (PUF60) is essential for embryonic development and is overexpressed in multiple types of cancer. Here, we report that PUF60 promotes mitotic cell cycle and lung cancer progression by controlling AS of the cell division cycle 25C (CDC25C). Systematic analysis of splicing factors deregulated in lung adenocarcinoma (LUAD) identifies that elevated copy number and expression of PUF60 correlate with poor prognosis. PUF60 depletion inhibits LUAD cell-cycle G2/M transition, cell proliferation, and tumor development. Mechanistically, PUF60 knockdown leads to exon skipping enriched in mitotic cell cycle genes, including CDC25C. Exon 3 skipping in the full-length CDC25C results in nonsense-mediated mRNA decay and a decrease of CDC25C protein, thereby inhibiting cell proliferation. This study establishes PUF60 as a cell cycle regulator and an oncogenic splicing factor in lung cancer.

Molecular basis of RNA-binding and autoregulation by the cancer-associated splicing factor RBM39

Pharmacologic depletion of RNA-binding motif 39 (RBM39) using aryl sulfonamides represents a promising anti-cancer therapy but requires high levels of the adaptor protein DCAF15. Consequently, novel approaches to deplete RBM39 in an DCAF15-independent manner are required. Here, we uncover that RBM39 autoregulates via the inclusion of a poison exon into its own pre-mRNA and identify the cis-acting elements that govern this regulation. We also determine the NMR solution structures of RBM39's tandem RNA recognition motifs (RRM1 and RRM2) bound to their respective RNA targets, revealing how RRM1 recognises RNA stem loops whereas RRM2 binds specifically to single-stranded N(G/U)NUUUG. Our results support a model where RRM2 selects the 3'-splice site of a poison exon and the RRM3 and RS domain stabilise the U2 snRNP at the branchpoint. Our work provides molecular insights into RBM39-dependent 3'-splice site selection and constitutes a solid basis to design alternative anti-cancer therapies.

Alu RNA fold links splicing with signal recognition particle proteins

Transcriptomic diversity in primates was considerably expanded by exonizations of intronic Alu elements. To better understand their cellular mechanisms we have used structure-based mutagenesis coupled with functional and proteomic assays to study the impact of successive primate mutations and their combinations on inclusion of a sense-oriented AluJ exon in the human F8 gene. We show that the splicing outcome was better predicted by consecutive RNA conformation changes than by computationally derived splicing regulatory motifs. We also demonstrate an involvement of SRP9/14 (signal recognition particle) heterodimer in splicing regulation of Alu-derived exons. Nucleotide substitutions that accumulated during primate evolution relaxed the conserved left-arm AluJ structure including helix H1 and reduced the capacity of SRP9/14 to stabilize the closed Alu conformation. RNA secondary structure-constrained mutations that promoted open Y-shaped conformations of the Alu made the Alu exon inclusion reliant on DHX9. Finally, we identified additional SRP9/14 sensitive Alu exons and predicted their functional roles in the cell. Together, these results provide unique insights into architectural elements required for sense Alu exonization, identify conserved pre-mRNA structures involved in exon selection and point to a possible chaperone activity of SRP9/14 outside the mammalian signal recognition particle.

Alternative splicing in bladder cancer: potential strategies for cancer diagnosis, prognosis, and treatment

Bladder cancer is the most common malignancy of the urinary tract worldwide. The therapeutic options to tackle this disease comprise surgery, intravesical or systemic chemotherapy, and immunotherapy. Unfortunately, a wide number of patients ultimately become resistant to these treatments and develop aggressive metastatic disease, presenting a poor prognosis. Therefore, the identification of novel therapeutic approaches to tackle this devastating pathology is urgently needed. However, a significant limitation is that the progression and drug response of bladder cancer is strongly associated with its intrinsic molecular heterogeneity. In this sense, RNA splicing is recently gaining importance as a critical hallmark of cancer since can have a significant clinical value. In fact, a profound dysregulation of the splicing process has been reported in bladder cancer, especially in the expression of certain key splicing variants and circular RNAs with a potential clinical value as diagnostic/prognostic biomarkers or therapeutic targets in this pathology. Indeed, some authors have already evidenced a profound antitumor effect by targeting some splicing factors (e.g., PTBP1), mRNA splicing variants (e.g., PKM2, HYAL4-v1), and circular RNAs (e.g., circITCH, circMYLK), which illustrates new possibilities to significantly improve the management of this pathology. This review represents the first detailed overview of the splicing process and its alterations in bladder cancer, and highlights opportunities for the development of novel diagnostic/prognostic biomarkers and their clinical potential for the treatment of this devastating cancer type. This article is categorized under: RNA Processing > Splicing Regulation/Alternative Splicing RNA in Disease and Development > RNA in Disease.

Insights on the biological functions and diverse regulation of RNA-binding protein 39 and their implication in human diseases

RNA-binding protein 39 (RBM39) involves in pre-mRNA splicing and transcriptional regulation. RBM39 is dysregulated in many cancers and its upregulation enhances cancer cell proliferation. Recently, it has been discovered that aryl sulfonamides act as molecular glues to recruit RBM39 to the CRL4^DCAF15 E3 ubiquitin ligase complex for its ubiquitination and proteasomal degradation. Therefore, various studies have focused on the degradation of RBM39 by aryl sulfonamides in the aim of finding new cancer therapeutics. These discoveries also attracted focus for thorough study on the biological functions of RBM39. RBM39 was found to regulate the splicing and transcription of genes mainly involved in pre-mRNA splicing, cell cycle regulation, DNA damage response, and metabolism, but the understanding of these regulations is still in its infancy. This article reviews the advances of the current literature and discusses the remaining key issues on the biological function and dynamic regulation of RBM39 at the post-translational level.

RNA splicing dysregulation and the hallmarks of cancer

Dysregulated RNA splicing is a molecular feature that characterizes almost all tumour types. Cancer-associated splicing alterations arise from both recurrent mutations and altered expression of trans-acting factors governing splicing catalysis and regulation. Cancer-associated splicing dysregulation can promote tumorigenesis via diverse mechanisms, contributing to increased cell proliferation, decreased apoptosis, enhanced migration and metastatic potential, resistance to chemotherapy and evasion of immune surveillance. Recent studies have identified specific cancer-associated isoforms that play critical roles in cancer cell transformation and growth and demonstrated the therapeutic benefits of correcting or otherwise antagonizing such cancer-associated mRNA isoforms. Clinical-grade small molecules that modulate or inhibit RNA splicing have similarly been developed as promising anticancer therapeutics. Here, we review splicing alterations characteristic of cancer cell transcriptomes, dysregulated splicing's contributions to tumour initiation and progression, and existing and emerging approaches for targeting splicing for cancer therapy. Finally, we discuss the outstanding questions and challenges that must be addressed to translate these findings into the clinic.

The transcription factor c-Jun inhibits RBM39 to reprogram pre-mRNA splicing during genotoxic stress

Genotoxic agents, that are used in cancer therapy, elicit the reprogramming of the transcriptome of cancer cells. These changes reflect the cellular response to stress and underlie some of the mechanisms leading to drug resistance. Here, we profiled genome-wide changes in pre-mRNA splicing induced by cisplatin in breast cancer cells. Among the set of cisplatin-induced alternative splicing events we focused on COASY, a gene encoding a mitochondrial enzyme involved in coenzyme A biosynthesis. Treatment with cisplatin induces the production of a short isoform of COASY lacking exons 4 and 5, whose depletion impedes mitochondrial function and decreases sensitivity to cisplatin. We identified RBM39 as a major effector of the cisplatin-induced effect on COASY splicing. RBM39 also controls a genome-wide set of alternative splicing events partially overlapping with the cisplatin-mediated ones. Unexpectedly, inactivation of RBM39 in response to cisplatin involves its interaction with the AP-1 family transcription factor c-Jun that prevents RBM39 binding to pre-mRNA. Our findings therefore uncover a novel cisplatin-induced interaction between a splicing regulator and a transcription factor that has a global impact on alternative splicing and contributes to drug resistance.

The diverse pleiotropic effects of spliceosomal protein PUF60: A case series of Verheij syndrome

Verheij syndrome (VRJS) is a rare craniofacial spliceosomopathy presenting with craniofacial dysmorphism, multiple congenital anomalies and variable neurodevelopmental delay. It is caused by single nucleotide variants (SNVs) in PUF60 or interstitial deletions of the 8q24.3 region. PUF60 encodes a splicing factor which forms part of the spliceosome. To date, 36 patients with a sole diagnosis of VRJS due to disease-causing PUF60 SNVs have been reported in peer-reviewed publications. Although the depth of their phenotyping has varied greatly, they exhibit marked phenotypic heterogeneity. We report 10 additional unrelated patients, including the first described patients of Khmer, Indian, and Vietnamese ethnicities, and the eldest patient to date, with 10 heterozygous PUF60 variants identified through exome sequencing, 8 previously unreported. All patients underwent deep phenotyping identifying variable dysmorphism, growth delay, neurodevelopmental delay, and multiple congenital anomalies, including several unique features. The eldest patient is the only reported individual with a germline variant and neither neurodevelopmental delay nor intellectual disability. In combining these detailed phenotypic data with that of previously reported patients (n = 46), we further refine the known frequencies of features associated with VRJS. These include neurodevelopmental delay/intellectual disability (98%), axial skeletal anomalies (74%), appendicular skeletal anomalies (73%), oral anomalies (68%), short stature (66%), cardiac anomalies (63%), brain malformations (48%), hearing loss (46%), microcephaly (41%), colobomata (38%), and other ocular anomalies (65%). This case series, incorporating three patients from previously unreported ethnic backgrounds, further delineates the broad pleiotropy and mutational spectrum of PUF60 pathogenic variants.

Decreased spliceosome fidelity and egl-8 intron retention inhibit mTORC1 signaling to promote longevity

Changes in splicing fidelity are associated with loss of homeostasis and aging, yet only a handful of splicing factors have been shown to be causally required to promote longevity, and the underlying mechanisms and downstream targets in these paradigms remain elusive. Surprisingly, we found a hypomorphic mutation within ribonucleoprotein RNP-6/poly(U)-binding factor 60 kDa (PUF60), a spliceosome component promoting weak 3′-splice site recognition, which causes aberrant splicing, elevates stress responses and enhances longevity in Caenorhabditis elegans. Through genetic suppressor screens, we identify a gain-of-function mutation within rbm-39, an RNP-6-interacting splicing factor, which increases nuclear speckle formation, alleviates splicing defects and curtails longevity caused by rnp-6 mutation. By leveraging the splicing changes induced by RNP-6/RBM-39 activities, we uncover intron retention in egl-8/phospholipase C β4 (PLCB4) as a key splicing target prolonging life. Genetic and biochemical evidence show that neuronal RNP-6/EGL-8 downregulates mammalian target of rapamycin complex 1 (mTORC1) signaling to control organismal lifespan. In mammalian cells, PUF60 downregulation also potently and specifically inhibits mTORC1 signaling. Altogether, our results reveal that splicing fidelity modulates lifespan through mTOR signaling.

Systematic exploration of dynamic splicing networks reveals conserved multistage regulators of neurogenesis

Alternative splicing (AS) is a critical regulatory layer; yet, factors controlling functionally coordinated splicing programs during developmental transitions are poorly understood. Here, we employ a screening strategy to identify factors controlling dynamic splicing events important for mammalian neurogenesis. Among previously unknown regulators, Rbm38 acts widely to negatively control neural AS, in part through interactions mediated by the established repressor of splicing, Ptbp1. Puf60, a ubiquitous factor, is surprisingly found to promote neural splicing patterns. This activity requires a conserved, neural-differential exon that remodels Puf60 co-factor interactions. Ablation of this exon rewires distinct AS networks in embryonic stem cells and at different stages of mouse neurogenesis. Single-cell transcriptome analyses further reveal distinct roles for Rbm38 and Puf60 isoforms in establishing neuronal identity. Our results describe important roles for previously unknown regulators of neurogenesis and establish how an alternative exon in a widely expressed splicing factor orchestrates temporal control over cell differentiation.

Exonic splicing code and protein binding sites for calcium

Auxilliary splicing sequences in exons, known as enhancers (ESEs) and silencers (ESSs), have been subject to strong selection pressures at the RNA and protein level. The protein component of this splicing code is substantial, recently estimated at ∼50% of the total information within ESEs, but remains poorly understood. The ESE/ESS profiles were previously associated with the Irving-Williams (I-W) stability series for divalent metals, suggesting that the ESE/ESS evolution was shaped by metal binding sites. Here, we have examined splicing activities of exonic sequences that encode protein binding sites for Ca2+, a weak binder in the I-W affinity order. We found that predicted exon inclusion levels for the EF-hand motifs and for Ca2+-binding residues in nonEF-hand proteins were higher than for average exons. For canonical EF-hands, the increase was centred on the EF-hand chelation loop and, in particular, on Ca2+-coordinating residues, with a 1>12>3∼5>9 hierarchy in the 12-codon loop consensus and usage bias at codons 1 and 12. The same hierarchy but a lower increase was observed for noncanonical EF-hands, except for S100 proteins. EF-hand loops preferentially accumulated exon splits in two clusters, one located in their N-terminal halves and the other around codon 12. Using splicing assays and published crosslinking and immunoprecipitation data, we identify candidate trans-acting factors that preferentially bind conserved GA-rich motifs encoding negatively charged amino acids in the loops. Together, these data provide evidence for the high capacity of codons for Ca2+-coordinating residues to be retained in mature transcripts, facilitating their exon-level expansion during eukaryotic evolution.

CLOCK regulates Drp1 mRNA stability and mitochondrial homeostasis by interacting with PUF60

Circadian genes such as Clock, Bmal1, Cryptochrome1/2, and Period1/2/3 constitute the precise circadian system. Clock^Δ19 is a commonly used mouse model harboring a circadian clock gene mutation, which lacks the EXON-19-encoded 51 amino acids. Previous reports have shown that Clock^Δ19 mice have severe metabolic abnormalities. Here, we report that the mitochondria of Clock^Δ19 mice exhibit excessive fission and dysfunction. We also demonstrate that CLOCK binds to the RNA-binding protein PUF60 through its EXON 19. Further, we find that PUF60 directly maintains mitochondrial homeostasis through regulating Drp1 mRNA stability, while the association with CLOCK can competitively inhibit this function. In Clock^Δ19 mice, CLOCK^Δ19 releases PUF60, leading to enhanced Drp1 mRNA stability and persistent mitochondrial fission. Our results reveal a direct post-transcriptional role of CLOCK in regulating mitochondrial homeostasis via Drp1 mRNA stability and that the loss of EXON 19 of CLOCK in Clock^Δ19 mice leads to severe mitochondrial homeostasis disorders.

Integrative proteome analysis implicates aberrant RNA splicing in impaired developmental potential of aged mouse oocytes

Aging has many effects on the female reproductive system, among which decreased oocyte quality and impaired embryo developmental potential are the most important factors affecting female fertility. However, the mechanisms underlying oocyte aging are not yet fully understood. Here, we selected normal reproductively aging female mice and constructed a protein expression profile of metaphase II (MII) oocytes from three age groups. A total of 187 differentially expressed (DE) proteins were identified, and bioinformatics analyses showed that these DE proteins were highly enriched in RNA splicing. Next, RNA‐seq was performed on 2‐cell embryos from these three age groups, and splicing analysis showed that a large number of splicing events and genes were discovered at this stage. Differentially spliced genes (DSGs) in the two reproductively aging groups versus the younger group were enriched in biological processes related to DNA damage repair/response. Binding motif analysis suggested that PUF60 might be one of the core splicing factors causing a decline in DNA repair capacity in the subsequent development of oocytes from reproductively aging mice, and changing the splicing pattern of its potential downstream DSG Cdk9 could partially mimic phenotypes in the reproductively aging groups. Taken together, our study suggested that the abnormal expression of splicing regulation proteins in aged MII oocytes would affect the splicing of nascent RNA after zygotic genome activation in 2‐cell embryos, leading to the production of abnormally spliced transcripts of some key genes associated with DNA damage repair/response, thus affecting the developmental potential of aged oocytes.

RNA-binding protein 39: a promising therapeutic target for cancer

RNA-binding motif protein 39 (RBM39), as a key factor in tumor-targeted mRNA and protein expression, not only plays a vital role in tumorigenesis, but also has broad development prospects in clinical treatment and drug research. Moreover, since RBM39 was identified as a target of sulfonamides, it has played a key role in the emerging field of molecule drug development. Hence, it is of great significance to study the interaction between RBM39 and tumors and the clinical application of drug-targeted therapy. In this paper, we describe the possible multi-level regulation of RBM39, including gene transcription, protein translation, and alternative splicing. Importantly, the molecular function of RBM39 as an important splicing factor in most common tumors is systematically outlined. Furthermore, we briefly introduce RBM39’s tumor-targeted drug research and its clinical application, hoping to give reference significance for the molecular mechanism of RBM39 in tumors, and provide reliable ideas for in-depth research for future therapeutic strategies.

Inhibiting an RBM39/MLL1 epigenomic regulatory complex with dominant-negative peptides disrupts cancer cell transcription and proliferation

RBM39 is a known splicing factor and coactivator. Here, we report that RBM39 functions as a master transcriptional regulator that interacts with the MLL1 complex to facilitate chromatin binding and H3K4 trimethylation in breast cancer cells. We identify RBM39 functional domains required for DNA and complex binding and show that the loss of RBM39 has widespread effects on H3K4me3 and gene expression, including key oncogenic pathways. RBM39's RNA recognition motif 3 (RRM3) functions as a dominant-negative domain; namely, it disrupts the complex and H3K4me trimethylation and expression of RBM/MLL1 target genes. RRM3-derived cell-penetrating peptides phenocopy the effects of the loss of RBM39 to decrease growth and survival of all major subtypes of breast cancer and yet are nontoxic to normal cells. These findings establish RBM39/MLL1 as a major contributor to the abnormal epigenetic landscape in breast cancer and lay the foundation for peptide-mediated cancer-specific therapy based on disruption of RBM39 epigenomic functions.

Restriction of an intron size en route to endothermy

Ca2+-insensitive and -sensitive E1 subunits of the 2-oxoglutarate dehydrogenase complex (OGDHC) regulate tissue-specific NADH and ATP supply by mutually exclusive OGDH exons 4a and 4b. Here we show that their splicing is enforced by distant lariat branch points (dBPs) located near the 5' splice site of the intervening intron. dBPs restrict the intron length and prevent transposon insertions, which can introduce or eliminate dBP competitors. The size restriction was imposed by a single dominant dBP in anamniotes that expanded into a conserved constellation of four dBP adenines in amniotes. The amniote clusters exhibit taxon-specific usage of individual dBPs, reflecting accessibility of their extended motifs within a stable RNA hairpin rather than U2 snRNA:dBP base-pairing. The dBP expansion took place in early terrestrial species and was followed by a uridine enrichment of large downstream polypyrimidine tracts in mammals. The dBP-protected megatracts permit reciprocal regulation of exon 4a and 4b by uridine-binding proteins, including TIA-1/TIAR and PUF60, which promote U1 and U2 snRNP recruitment to the 5' splice site and BP, respectively, but do not significantly alter the relative dBP usage. We further show that codons for residues critically contributing to protein binding sites for Ca2+ and other divalent metals confer the exon inclusion order that mirrors the Irving-Williams affinity series, linking the evolution of auxiliary splicing motifs in exons to metallome constraints. Finally, we hypothesize that the dBP-driven selection for Ca2+-dependent ATP provision by E1 facilitated evolution of endothermy by optimizing the aerobic scope in target tissues.

RNA-binding motif protein 39 (RBM39): An emerging cancer target

RNA-binding motif protein 39 (RBM39) is an RNA-binding protein involved in transcriptional co-regulation and alternative RNA splicing. Recent studies have revealed that RBM39 is the unexpected target of aryl sulphonamides, which act as molecular glues between RBM39 and the DCAF15-associated E3 ubiquitin ligase complex leading to selective degradation of the target. Loss of RBM39 leads to aberrant splicing events and differential gene expression, thereby inhibiting cell cycle progression and causing tumour regression in a number of preclinical models. Many clinical studies have shown that aryl sulphonamides were well tolerated, but their clinical performance was limited due to an insufficient understanding of the target, RBM39 biology and a lack of predictive biomarkers. This review summarises the current knowledge of RBM39 function and discusses the therapeutic potential of this spliceosome target in cancer therapy.

PUF60/AURKA Axis Contributes to Tumor Progression and Malignant Phenotypes in Bladder Cancer

Abnormal expression or mutation of RNA splicing proteins are widely observed in human cancers. Here, we identified poly(U) binding splicing factor 60 (PUF60) as one of the most differentially expressed genes out of 97 RNA splicing proteins between normal and bladder cancer tissues by bioinformatics analysis of TCGA bladder cancer expression data. The expression of PUF60 was significantly higher in tumor tissues, while high PUF60 expression was associated with malignant phenotypes of bladder cancer and shorter survival time. Moreover, we identified aurora kinase A (AURKA) as a new downstream target of PUF60 in bladder cancer cells. PUF60 knockdown significantly inhibited cell viability and colony formation capacity in bladder cancer cells, whereas AURKA overexpression reversed this inhibition effect. Overexpression of PUF60 significantly promoted cell viability and colony formation in bladder cancer cells, while treatment with AURKA specific inhibitor reversed this promotive effect. Mechanistically, PUF60 specifically bound to the AURKA promoter, thereby activating its transcription and expression. Furthermore, we showed that there was a significant positive correlation between PUF60 and AURKA expression in bladder cancer tissues, and PUF60 and AURKA expression contributed to tumor progression and malignant phenotypes in the patients with bladder cancer. Collectively, these results indicate that the PUF60/AURKA axis plays a key role in regulating tumorigenesis and progression of bladder cancer, and may be a potential prognostic biomarker and therapeutic target for bladder cancer patients.

Identification of spliceosome components pivotal to breast cancer survival

Cancer cells employ alternative splicing (AS) to acquire splicing isoforms favouring their survival. However, the causes of aberrant AS in breast cancer are poorly understood. In this study, the METABRIC (Molecular Taxonomy of Breast Cancer International Consortium) data were analysed with univariate feature selection. Of 122 analysed spliceosome components, U2SURP, PUF60, DDX41, HNRNPAB, EIF4A3, and PPIL3 were significantly associated with breast cancer survival. The top 4 four genes, U2SURP, PUF60, DDX41, and HNRNPAB, were chosen for further analyses. Their expression was significantly associated with cancer molecular subtype, tumour stage, tumour grade, overall survival (OS), and cancer-specific survival in the METABRIC data. These results were verifiable using other cohorts. The Cancer Genome Atlas data unveiled the elevated expression of PUF60, DDX41, and HNRNPAB in tumours compared with the normal tissue and confirmed the differential expression of the four genes among cancer molecular subtypes, as well as the associations of U2SURP, PUF60, and DDX41 expression with tumour stage. A meta-analysis data verified the associations of U2SURP, PUF60, and HNRNPAB expression with tumour grade, the associations of PUF60, DDX41, and HNRNPAB expression with OS and distant metastasis-free survival, and the associations of U2SURP and HNRNPAB expression with relapse-free survival. Experimentally, we demonstrated that inhibiting the expression of the four genes separately suppressed cell colony formation and slowed down cell growth considerably in breast cancer cells, but not in immortal breast epithelial cells. In conclusion, we have identified U2SURP, PUF60, DDX41, and HNRNPAB are spliceosome-related genes pivotal for breast cancer survival.

Structural disruption of exonic stem-loops immediately upstream of the intron regulates mammalian splicing

Recognition of highly degenerate mammalian splice sites by the core spliceosomal machinery is regulated by several protein factors that predominantly bind exonic splicing motifs. These are postulated to be single-stranded in order to be functional, yet knowledge of secondary structural features that regulate the exposure of exonic splicing motifs across the transcriptome is not currently available. Using transcriptome-wide RNA structural information we show that retained introns in mouse are commonly flanked by a short (≲70 nucleotide), highly base-paired segment upstream and a predominantly single-stranded exonic segment downstream. Splicing assays with select pre-mRNA substrates demonstrate that loops immediately upstream of the introns contain pre-mRNA-specific splicing enhancers, the substitution or hybridization of which impedes splicing. Additionally, the exonic segments flanking the retained introns appeared to be more enriched in a previously identified set of hexameric exonic splicing enhancer (ESE) sequences compared to their spliced counterparts, suggesting that base-pairing in the exonic segments upstream of retained introns could be a means for occlusion of ESEs. The upstream exonic loops of the test substrate promoted recruitment of splicing factors and consequent pre-mRNA structural remodeling, leading up to assembly of the early spliceosome. These results suggest that disruption of exonic stem-loop structures immediately upstream (but not downstream) of the introns regulate alternative splicing events, likely through modulating accessibility of splicing factors.

Protein elongation variant of PUF60: Milder phenotypic end of the Verheij syndrome

The PUF60 gene encodes a ubiquitously expressed essential splicing factor that is recruited to the U2snRNA complex. The complex binds to the 3' splice site of exons in specific target genes and regulates the inclusion or exclusion of such exons. Recently, pathogenic variants of PUF60 have been shown to cause a relatively specific and potentially recognizable pattern of malformation referred to as Verheij syndrome. Here, we report a 12-year-old female patient with a de novo mutation in PUF60 whose phenotype was representative of the milder end of the phenotypic spectrum of Verheij syndrome; the de novo mutation was a frameshift mutation p.(Ser558Cysfs*21) that resulted in the addition of 21 extra amino acids at the carboxy end of the protein. Among the frequent features of Verheij syndrome, the patient exhibited coloboma, cervical spinal segmentation defects, and borderline intellectual functioning, but lacked cardiac abnormalities, deafness, and urogenital abnormalities. The results of RNA analysis using peripheral blood showed the escape of the mutant allele from nonsense-mediated mRNA decay, possibly accounting for the mild phenotype in the presently reported patient. Based on our clinical observations, we inferred that two embryologic processes, closure of the ocular plate and cervical spinal segmentation, are particularly susceptible to deficient PUF60-mediated splicing regulation, compared with other embryogenetic processes leading to the central nervous system, heart, ear, and kidney.

Intrinsic Regulatory Role of RNA Structural Arrangement in Alternative Splicing Control

Alternative splicing is a highly sophisticated process, playing a significant role in posttranscriptional gene expression and underlying the diversity and complexity of organisms. Its regulation is multilayered, including an intrinsic role of RNA structural arrangement which undergoes time- and tissue-specific alterations. In this review, we describe the principles of RNA structural arrangement and briefly decipher its cis- and trans-acting cellular modulators which serve as crucial determinants of biological functionality of the RNA structure. Subsequently, we engage in a discussion about the RNA structure-mediated mechanisms of alternative splicing regulation. On one hand, the impairment of formation of optimal RNA structures may have critical consequences for the splicing outcome and further contribute to understanding the pathomechanism of severe disorders. On the other hand, the structural aspects of RNA became significant features taken into consideration in the endeavor of finding potential therapeutic treatments. Both aspects have been addressed by us emphasizing the importance of ongoing studies in both fields.

Cancer-Associated Substitutions in RNA Recognition Motifs of PUF60 and U2AF65 Reveal Residues Required for Correct Folding and 3' Splice-Site Selection

U2AF65 (U2AF2) and PUF60 (PUF60) are splicing factors important for recruitment of the U2 small nuclear ribonucleoprotein to lariat branch points and selection of 3′ splice sites (3′ss). Both proteins preferentially bind uridine-rich sequences upstream of 3′ss via their RNA recognition motifs (RRMs). Here, we examined 36 RRM substitutions reported in cancer patients to identify variants that alter 3′ss selection, RNA binding and protein properties. Employing PUF60- and U2AF65-dependent 3′ss previously identified by RNA-seq of depleted cells, we found that 43% (10/23) and 15% (2/13) of independent RRM mutations in U2AF65 and PUF60, respectively, conferred splicing defects. At least three RRM mutations increased skipping of internal U2AF2 (~9%, 2/23) or PUF60 (~8%, 1/13) exons, indicating that cancer-associated RRM mutations can have both cis- and trans-acting effects on splicing. We also report residues required for correct folding/stability of each protein and map functional RRM substitutions on to existing high-resolution structures of U2AF65 and PUF60. These results identify new RRM residues critical for 3′ss selection and provide relatively simple tools to detect clonal RRM mutations that enhance the mRNA isoform diversity.

Anti-parasite drug ivermectin can suppress ovarian cancer by regulating lncRNA-EIF4A3-mRNA axes

RELEVANCE

Ivermectin, as an old anti-parasite drug, can suppress almost completely the growth of various human cancers, including ovarian cancer (OC). However, its anticancer mechanism remained to be further studied at the molecular levels. Ivermectin-related molecule-panel changes will serve a useful tool for its personalized drug therapy and prognostic assessment in OCs.

PURPOSE

To explore the functional significance of ivermectin-mediated lncRNA-EIF4A3-mRNA axes in OCs and ivermectin-related molecule-panel for its personalized drug therapy monitoring.

METHODS

Based on our previous study, a total of 16 lncRNA expression patterns were analyzed using qRT-PCR before and after ivermectin-treated different OC cell lines (TOV-21G and A2780). Stable isotope labeling with amino acids in cell culture (SILAC)-based quantitative proteomics was used to analyze the protein expressions of EIF4A3 and EIF4A3-binding mRNAs in ovarian cancer cells treated with and without ivermectin. A total of 411 OC patients from the Cancer Genome Atlas (TCGA) database with the selected lncRNA expressions and the corresponding clinical data were included. Lasso regression was constructed to examine the relationship between lncRNA signature and OC survival risk. The overall survival analysis between high-risk and low-risk groups used the Kaplan-Meier method. Heatmap showed the correlation between risk groups and clinical characteristics. The univariate and multivariate models were established with Cox regression.

RESULTS

SILAC-based quantitative proteomics found the protein expression levels of EIF4A3 and 116 EIF4A3-binding mRNAs were inhibited by ivermectin in OC cells. Among the analyzed 16 lncRNAs (HCG15, KIF9-AS1, PDCD4-AS1, ZNF674-AS1, ZNRF3-AS1, SOS1-IT1, LINC00565, SNHG3, PLCH1-AS1, WWTR1-AS1, LINC00517, AL109767.1, STARD13-IT1, LBX2-AS1, LEMD1-AS1, and HOXC-AS3), only 7 lncRNAs (HCG15, KIF9-AS1, PDCD4-AS1, ZNF674-AS1, ZNRF3-AS1, SOS1-IT1, and LINC00565) were obtained for further lasso regression when combined with the results of drug testing and overall survival analysis. Lasso regression identified the prognostic model of ivermectin-related three-lncRNA signature (ZNRF3-AS1, SOS1-IT1, and LINC00565). The high-risk and low-risk groups based on the prognostic model were significantly related to overall survival and clinicopathologic characteristics (survival status, lymphatic invasion, cancer status, and clinical stage) in OC patients and remained independent risk factors according to multivariate COX analysis (p < 0.05).

CONCLUSION

Those findings provided the potential targeted lncRNA-EIF4A3-mRNA pathways of ivermectin in OC, and constructed the effective prognostic model, which benefits discovery of novel mechanism of ivermectin to suppress ovarian cancer cells, and the ivermectin-related molecule-panel changes benefit for its personalized drug therapy and prognostic assessment towards its predictive, preventive, and personalized medicine (PPPM) in OCs.

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.

U2AF65 assemblies drive sequence-specific splice site recognition

The essential splicing factor U2AF65 is known to help anchoring U2 snRNP at the branch site. Its C-terminal UHM domain interacts with ULM motifs of SF3b155, an U2 snRNP protein. Here, we report a cooperative binding of U2AF65 and the related protein CAPERα to the multi-ULM domain of SF3b155. In addition, we show that the RS domain of U2AF65 drives a liquid-liquid phase separation that is amplified by intronic RNA with repeated pyrimidine tracts. In cells, knockdown of either U2AF65 or CAPERα improves the inclusion of cassette exons that are preceded by such repeated pyrimidine-rich motifs. These results support a model in which liquid-like assemblies of U2AF65 and CAPERα on repetitive pyrimidine-rich RNA sequences are driven by their RS domains, and facilitate the recruitment of the multi-ULM domain of SF3b155. We anticipate that posttranslational modifications and proteins recruited in dynamical U2AF65 and CAPERα condensates may further contribute to the complex mechanisms leading to specific splice site choice that occurs in cells.