Molecular impact of mutations in RNA splicing factors in cancer.

Somatic mutations in genes encoding components of the RNA splicing machinery occur frequently in multiple forms of cancer. The most frequently mutated RNA splicing factors in cancer impact intronic branch site and 3' splice site recognition. These include mutations in the core RNA splicing factor SF3B1 as well as mutations in the U2AF1/2 heterodimeric complex, which recruits the SF3b complex to the 3' splice site. Additionally, mutations in splicing regulatory proteins SRSF2 and RBM10 are frequent in cancer, and there has been a recent suggestion that variant forms of small nuclear RNAs (snRNAs) may contribute to splicing dysregulation in cancer. Here, we describe molecular mechanisms by which mutations in these factors alter splice site recognition and how studies of this process have yielded new insights into cancer pathogenesis and the molecular regulation of splicing. We also discuss data linking mutant RNA splicing factors to RNA metabolism beyond splicing.

The splicing factor CCAR1 regulates the Fanconi anemia/BRCA pathway.

The twenty-three Fanconi anemia (FA) proteins cooperate in the FA/BRCA pathway to repair DNA interstrand cross-links (ICLs). The cell division cycle and apoptosis regulator 1 (CCAR1) protein is also a regulator of ICL repair, though its possible function in the FA/BRCA pathway remains unknown. Here, we demonstrate that CCAR1 plays a unique upstream role in the FA/BRCA pathway and is required for FANCA protein expression in human cells. Interestingly, CCAR1 co-immunoprecipitates with FANCA pre-mRNA and is required for FANCA mRNA processing. Loss of CCAR1 results in retention of a poison exon in the FANCA transcript, thereby leading to reduced FANCA protein expression. A unique domain of CCAR1, the EF hand domain, is required for interaction with the U2AF heterodimer of the spliceosome and for excision of the poison exon. Taken together, CCAR1 is a splicing modulator required for normal splicing of the FANCA mRNA and other mRNAs involved in various cellular pathways.

Precision analysis of mutant U2AF1 activity reveals deployment of stress granules in myeloid malignancies.

Splicing factor mutations are common among cancers, recently emerging as drivers of myeloid malignancies. U2AF1 carries hotspot mutations in its RNA-binding motifs; however, how they affect splicing and promote cancer remain unclear. The U2AF1/U2AF2 heterodimer is critical for 3' splice site (3'SS) definition. To specifically unmask changes in U2AF1 function in vivo, we developed a crosslinking and immunoprecipitation procedure that detects contacts between U2AF1 and the 3'SS AG at single-nucleotide resolution. Our data reveal that the U2AF1 S34F and Q157R mutants establish new 3'SS contacts at -3 and +1 nucleotides, respectively. These effects compromise U2AF2-RNA interactions, resulting predominantly in intron retention and exon exclusion. Integrating RNA binding, splicing, and turnover data, we predicted that U2AF1 mutations directly affect stress granule components, which was corroborated by single-cell RNA-seq. Remarkably, U2AF1-mutant cell lines and patient-derived MDS/AML blasts displayed a heightened stress granule response, pointing to a novel role for biomolecular condensates in adaptive oncogenic strategies.

Activation of ERß hijacks the splicing machinery to trigger R-loop formation in triple-negative breast cancer.

Triple-negative breast cancer (TNBC) is a subtype of breast cancer with aggressive behavior and poor prognosis. Current therapeutic options available for TNBC patients are primarily chemotherapy. With our evolving understanding of this disease, novel targeted therapies, including poly ADP-ribose polymerase (PARP) inhibitors, antibody-drug conjugates, and immune-checkpoint inhibitors, have been developed for clinical use. Previous reports have demonstrated the essential role of estrogen receptor ß (ERß) in TNBC, but the detailed molecular mechanisms downstream ERß activation in TNBC are still far from elucidated. In this study, we demonstrated that a specific ERß agonist, LY500307, potently induces R-loop formation and DNA damage in TNBC cells. Subsequent interactome experiments indicated that the residues 151 to 165 of U2 small nuclear RNA auxiliary factor 1 (U2AF1) and the Trp439 and Lys443 of ERß were critical for the binding between U2AF1 and ERß. Combined RNA sequencing and ribosome sequencing analysis demonstrated that U2AF1-regulated downstream RNA splicing of 5-oxoprolinase (OPLAH) could affect its enzymatic activity and is essential for ERß-induced R-loop formation and DNA damage. In clinical samples including 115 patients from The Cancer Genome Atlas (TCGA) and 32 patients from an in-house cohort, we found a close correlation in the expression of ESR2 and U2AF1 in TNBC patients. Collectively, our study has unraveled the molecular mechanisms that explain the therapeutic effects of ERß activation in TNBC, which provides rationale for ERß activation-based single or combined therapy for patients with TNBC.

A systematic screen identifies Saf5 as a link between splicing and transcription in fission yeast.

Splicing is an important step of gene expression regulation in eukaryotes, as there are many mRNA precursors that can be alternatively spliced in different tissues, at different cell cycle phases or under different external stimuli. We have developed several integrated fluorescence-based in vivo splicing reporter constructs that allow the quantification of fission yeast splicing in vivo on intact cells, and we have compared their splicing efficiency in a wild type strain and in a prp2-1 (U2AF65) genetic background, showing a clear dependency between Prp2 and a consensus signal at 5' splicing site (5'SS). To isolate novel genes involved in regulated splicing, we have crossed the reporter showing more intron retention with the Schizosaccharomyces pombe knock out collection. Among the candidate genes involved in the regulation of splicing, we have detected strong splicing defects in two of the mutants -Δcwf12, a member of the NineTeen Complex (NTC) and Δsaf5, a methylosome subunit that acts together with the survival motor neuron (SMN) complex in small nuclear ribonucleoproteins (snRNP) biogenesis. We have identified that strains with mutations in cwf12 have inefficient splicing, mainly when the 5'SS differs from the consensus. However, although Δsaf5 cells also have some dependency on 5'SS sequence, we noticed that when one intron of a given pre-mRNA was affected, the rest of the introns of the same pre-mRNA had high probabilities of being also affected. This observation points Saf5 as a link between transcription rate and splicing.

The splicing factor U2AF1 contributes to cancer progression through a noncanonical role in translation regulation.

Somatic mutations in the genes encoding components of the spliceosome occur frequently in human neoplasms, including myeloid dysplasias and leukemias, and less often in solid tumors. One of the affected factors, U2AF1, is involved in splice site selection, and the most common change, S34F, alters a conserved nucleic acid-binding domain, recognition of the 3' splice site, and alternative splicing of many mRNAs. However, the role that this mutation plays in oncogenesis is still unknown. Here, we uncovered a noncanonical function of U2AF1, showing that it directly binds mature mRNA in the cytoplasm and negatively regulates mRNA translation. This splicing-independent role of U2AF1 is altered by the S34F mutation, and polysome profiling indicates that the mutation affects translation of hundreds of mRNA. One functional consequence is increased synthesis of the secreted chemokine interleukin 8, which contributes to metastasis, inflammation, and cancer progression in mice and humans.

The Augmented R-Loop Is a Unifying Mechanism for Myelodysplastic Syndromes Induced by High-Risk Splicing Factor Mutations.

Mutations in several general pre-mRNA splicing factors have been linked to myelodysplastic syndromes (MDSs) and solid tumors. These mutations have generally been assumed to cause disease by the resultant splicing defects, but different mutations appear to induce distinct splicing defects, raising the possibility that an alternative common mechanism is involved. Here we report a chain of events triggered by multiple splicing factor mutations, especially high-risk alleles in SRSF2 and U2AF1, including elevated R-loops, replication stress, and activation of the ataxia telangiectasia and Rad3-related protein (ATR)-Chk1 pathway. We further demonstrate that enhanced R-loops, opposite to the expectation from gained RNA binding with mutant SRSF2, result from impaired transcription pause release because the mutant protein loses its ability to extract the RNA polymerase II (Pol II) C-terminal domain (CTD) kinase-the positive transcription elongation factor complex (P-TEFb)-from the 7SK complex. Enhanced R-loops are linked to compromised proliferation of bone-marrow-derived blood progenitors, which can be partially rescued by RNase H overexpression, suggesting a direct contribution of augmented R-loops to the MDS phenotype.

Global analysis of binding sites of U2AF1 and ZRSR2 reveals RNA elements required for mutually exclusive splicing by the U2- and U12-type spliceosome.

Recurring mutations in genes encoding 3' splice-site recognition proteins, U2AF1 and ZRSR2 are associated with human cancers. Here, we determined binding sites of the proteins to reveal that U2-type and U12-type splice sites are recognized by U2AF1 and ZRSR2, respectively. However, some sites are spliced by both the U2-type and U12-type spliceosomes, indicating that well-conserved consensus motifs in some U12-type introns could be recognized by the U2-type spliceosome. Nucleotides flanking splice sites of U12-type introns are different from those flanking U2-type introns. Remarkably, the AG dinucleotide at the positions -1 and -2 of 5' splice sites of U12-type introns with GT-AG termini is not present. AG next to 5' splice site introduced by a single nucleotide substitution at the -2 position could convert a U12-type splice site to a U2-type site. The class switch of introns by a single mutation and the bias against G at the -1 position of U12-type 5' splice site support the notion that the identities of nucleotides in exonic regions adjacent to splice sites are fine-tuned to avoid recognition by the U2-type spliceosome. These findings may shed light on the mechanism of selectivity in U12-type intron splicing and the mutations that affect splicing.

Venetoclax abrogates the prognostic impact of splicing factor gene mutations in newly diagnosed acute myeloid leukemia.

Mutations in splicing factor (SF) genes SRSF2, U2AF1, SF3B1, and ZRSR2 are now considered adverse risk in the European LeukemiaNet 2022 acute myeloid leukemia (AML) risk stratification. The prognostic impact of SF mutations in AML has been predominantly derived from younger patients treated with intensive (INT) therapy. We evaluated 994 patients with newly diagnosed AML, including 266 (27%) with a SFmut. Median age was 67 years overall, with patients with SFmut being older at 72 years. SRSF2 (n = 140, 53%) was the most common SFmut. In patients treated with INT, median relapse-free survival (RFS) (9.6 vs 21.4 months, P = .04) and overall survival (OS) (15.9 vs 26.7 months, P = .06) were shorter for patients with SFmut than without SFwt, however this significance abrogated when evaluating patients who received venetoclax with INT therapy (RFS 15.4 vs 20.3 months, P = .36; OS 19.6 vs 30.7 months, P = .98). In patients treated with LI, median RFS (9.3 vs 7.7 months, P = .35) and OS (12.3 vs 8.5 months, P = .14) were similar for patients with and without SFmut , and outcomes improved in all groups with venetoclax. On multivariate analysis, SFmut did not affect hazards of relapse and death for INT arm but reduced both these hazards in LI arm. In a large AML data set with >60% of patients receiving venetoclax with LI/INT therapy, SFmut had no independent negative prognostic impact. Newer prognostic models that consider LI therapy and use of venetoclax among other factors are warranted.

Body Temperature Cycles Control Rhythmic Alternative Splicing in Mammals.

The core body temperature of all mammals oscillates with the time of the day. However, direct molecular consequences of small, physiological changes in body temperature remain largely elusive. Here we show that body temperature cycles drive rhythmic SR protein phosphorylation to control an alternative splicing (AS) program. A temperature change of 1°C is sufficient to induce a concerted splicing switch in a large group of functionally related genes, rendering this splicing-based thermometer much more sensitive than previously described temperature-sensing mechanisms. AS of two exons in the 5' UTR of the TATA-box binding protein (Tbp) highlights the general impact of this mechanism, as it results in rhythmic TBP protein levels with implications for global gene expression in vivo. Together our data establish body temperature-driven AS as a core clock-independent oscillator in mammalian peripheral clocks.

Recurrent Neurodevelopmentally Associated Variants of the Pre-mRNA Splicing Factor U2AF2 Alter RNA Binding Affinities and Interactions.

De novo mutations affecting the pre-mRNA splicing factor U2AF2 are associated with developmental delays and intellectual disabilities, yet the molecular basis is unknown. Here, we demonstrated by fluorescence anisotropy RNA binding assays that recurrent missense mutants (Arg149Trp, Arg150His, or Arg150Cys) decreased the binding affinity of U2AF2 for a consensus splice site RNA. Crystal structures at 1.4 Å resolutions showed that Arg149Trp or Arg150His disrupted hydrogen bonds between U2AF2 and the terminal nucleotides of the RNA site. Reanalysis of publicly available RNaseq data confirmed that U2AF2 depletion altered splicing of transcripts encoding RNA binding proteins (RBPs). These results confirmed that the impaired RNA interactions of Arg149Trp and Arg150His U2AF2 variants could contribute to dysregulating an RBP-governed neurodevelopmental program of alternative splicing.

U2AF2-SNORA68 promotes triple-negative breast cancer stemness through the translocation of RPL23 from nucleoplasm to nucleolus and c-Myc expression.

BACKGROUND: Small nucleolar RNAs (snoRNAs) play key roles in ribosome biosynthesis. However, the mechanism by which snoRNAs regulate cancer stemness remains to be fully elucidated. METHODS: SNORA68 expression was evaluated in breast cancer tissues by in situ hybridization and qRT‒PCR. Proliferation, migration, apoptosis and stemness analyses were used to determine the role of SNORA68 in carcinogenesis and stemness maintenance. Mechanistically, RNA pull-down, RNA immunoprecipitation (RIP), cell fractionation and coimmunoprecipitation assays were conducted. RESULTS: SNORA68 exhibited high expression in triple-negative breast cancer (TNBC) and was significantly correlated with tumor size (P = 0.048), ki-67 level (P = 0.037), and TNM stage (P = 0.015). The plasma SNORA68 concentration was significantly lower in patients who achieved clinical benefit. The SNORA68-high patients had significantly shorter disease-free survival (DFS) (P = 0.036). Functionally, SNORA68 was found to promote the cell stemness and carcinogenesis of TNBC in vitro and in vivo. Furthermore, elevated SNORA68 expression led to increased nucleolar RPL23 expression and retained RPL23 in the nucleolus by binding U2AF2. RPL23 in the nucleolus subsequently upregulated c-Myc expression. This pathway was validated using a xenograft model. CONCLUSION: U2AF2-SNORA68 promotes TNBC stemness by retaining RPL23 in the nucleolus and increasing c-Myc expression, which provides new insight into the regulatory mechanism of stemness.

CCAR-1 works together with the U2AF large subunit UAF-1 to regulate alternative splicing.

The Cell Division Cycle and Apoptosis Regulator (CCAR) protein family members have recently emerged as regulators of alternative splicing and transcription, as well as having other key physiological functions. For example, mammalian CCAR2/DBC1 forms a complex with the zinc factor protein ZNF326 to integrate alternative splicing with RNA polymerase II transcriptional elongation in AT-rich regions of the DNA. Additionally, Caenorhabditis elegans CCAR-1, a homolog to mammalian CCAR2, facilitates the alternative splicing of the perlecan unc-52 gene. However, much about the CCAR family's role in alternative splicing is unknown. Here, we have examined the role of CCAR-1 in genome-wide alternative splicing in Caenorhabditis elegans and have identified new alternative splicing targets of CCAR-1 using RNA sequencing. Also, we found that CCAR-1 interacts with the spliceosome factors UAF-1 and UAF-2 using mass spectrometry, and that knockdown of ccar-1 affects alternative splicing patterns, motility, and proteostasis of UAF-1 mutant worms. Collectively, we demonstrate the role of CCAR-1 in regulating global alternative splicing in C. elegans and in conjunction with UAF-1.

U2AF35(S34F) Promotes Transformation by Directing Aberrant ATG7 Pre-mRNA 3' End Formation.

Recurrent mutations in the splicing factor U2AF35 are found in several cancers and myelodysplastic syndrome (MDS). How oncogenic U2AF35 mutants promote transformation remains to be determined. Here we derive cell lines transformed by the oncogenic U2AF35(S34F) mutant and identify aberrantly processed pre-mRNAs by deep sequencing. We find that in U2AF35(S34F)-transformed cells the autophagy-related factor 7 (Atg7) pre-mRNA is abnormally processed, which unexpectedly is not due to altered splicing but rather selection of a distal cleavage and polyadenylation (CP) site. This longer Atg7 mRNA is translated inefficiently, leading to decreased ATG7 levels and an autophagy defect that predisposes cells to secondary mutations, resulting in transformation. MDS and acute myeloid leukemia patient samples harboring U2AF35(S34F) have a similar increased use of the ATG7 distal CP site, and previous studies have shown that mice with hematopoietic cells lacking Atg7 develop an MDS-like syndrome. Collectively, our results reveal a basis for U2AF35(S34F) oncogenic activity.

Molecular Architecture of SF3b and Structural Consequences of Its Cancer-Related Mutations.

SF3b is a heptameric protein complex of the U2 small nuclear ribonucleoprotein (snRNP) that is essential for pre-mRNA splicing. Mutations in the largest SF3b subunit, SF3B1/SF3b155, are linked to cancer and lead to alternative branch site (BS) selection. Here we report the crystal structure of a human SF3b core complex, revealing how the distinctive conformation of SF3b155's HEAT domain is maintained by multiple contacts with SF3b130, SF3b10, and SF3b14b. Protein-protein crosslinking enabled the localization of the BS-binding proteins p14 and U2AF65 within SF3b155's HEAT-repeat superhelix, which together with SF3b14b forms a composite RNA-binding platform. SF3b155 residues, the mutation of which leads to cancer, contribute to the tertiary structure of the HEAT superhelix and its surface properties in the proximity of p14 and U2AF65. The molecular architecture of SF3b reveals the spatial organization of cancer-related SF3b155 mutations and advances our understanding of their effects on SF3b structure and function.

Cooperative engagement and subsequent selective displacement of SR proteins define the pre-mRNA 3D structural scaffold for early spliceosome assembly.

We recently reported that serine-arginine-rich (SR) protein-mediated pre-mRNA structural remodeling generates a pre-mRNA 3D structural scaffold that is stably recognized by the early spliceosomal components. However, the intermediate steps between the free pre-mRNA and the assembled early spliceosome are not yet characterized. By probing the early spliceosomal complexes in vitro and RNA-protein interactions in vivo, we show that the SR proteins bind the pre-mRNAs cooperatively generating a substrate that recruits U1 snRNP and U2AF65 in a splice signal-independent manner. Excess U1 snRNP selectively displaces some of the SR protein molecules from the pre-mRNA generating the substrate for splice signal-specific, sequential recognition by U1 snRNP, U2AF65 and U2AF35. Our work thus identifies a novel function of U1 snRNP in mammalian splicing substrate definition, explains the need for excess U1 snRNP compared to other U snRNPs in vivo, demonstrates how excess SR proteins could inhibit splicing, and provides a conceptual basis to examine if this mechanism of splicing substrate definition is employed by other splicing regulatory proteins.

[Primary myelofibrosis with double mutation in U2AF1].

A 47-year-old woman presented with subcutaneous hemorrhage. Blood tests revealed leukoerythroblastosis, anemia, and thrombocytopenia. Bone marrow biopsy led to a diagnosis of primary myelofibrosis (aaDIPSS, DIPSS-plus: intermediate-II risk). JAK2, CALR, and MPL mutations were not detected in peripheral blood, but targeted sequencing of bone marrow specimens revealed a double mutation (Q157R, S34F) in U2AF1. Allo-PBSCT was performed using an HLA-matched related donor, and post-transplantation bone marrow examination showed complete donor chimerism on day 55. Two years after allogeneic transplantation, the patient remains relapse-free. Although U2AF1 gene abnormality is known as a poor prognostic factor in primary myelofibrosis, this patient had a favorable long-term prognosis due to prompt transplantation therapy. This case highlights the importance of detailed gene mutation analysis in patients with triple-negative MF.

Promoter R-Loops Recruit U2AF1 to Modulate Its Phase Separation and RNA Splicing.

R-loops and guanine quadruplexes (G4s) are secondary structures of nucleic acids that are ubiquitously present in cells and are enriched in promoter regions of genes. By employing a bioinformatic approach based on overlap analysis of transcription factor chromatin immunoprecipitation sequencing (ChIP-seq) data sets, we found that many splicing factors, including U2AF1 whose recognition of the 3' splicing site is crucial for pre-mRNA splicing, exhibit pronounced enrichment at endogenous R-loop- and DNA G4-structure loci in promoter regions of human genes. We also revealed that U2AF1 binds directly to R-loops and DNA G4 structures at a low-nM binding affinity. Additionally, we showed the ability of U2AF1 to undergo phase separation, which could be stimulated by binding with R-loops, but not duplex DNA, RNA/DNA hybrid, DNA G4, or single-stranded RNA. We also demonstrated that U2AF1 binds to promoter R-loops in human cells, and this binding competes with U2AF1's interaction with 3' splicing site and leads to augmented distribution of RNA polymerase II (RNAPII) to promoters over gene bodies, thereby modulating cotranscriptional pre-mRNA splicing. Together, we uncovered a group of candidate proteins that can bind to both R-loops and DNA G4s, revealed the direct and strong interactions of U2AF1 with these nucleic acid structures, and established a biochemical rationale for U2AF1's occupancy in gene promoters. We also unveiled that interaction with R-loops promotes U2AF1's phase separation, and our work suggests that U2AF1 modulates pre-mRNA splicing by regulating RNAPII's partition in transcription initiation versus elongation.

A presumed missense variant in the U2AF2 gene causes exon skipping in neurodevelopmental diseases.

U2 small nuclear RNA auxiliary factor 2 (U2AF2) is an indispensable pre-mRNA splicing factor in the early process of splicing. Recently, U2AF2 was reported as a novel candidate gene associated with neurodevelopmental disorders. Herein, we report a patient with a novel presumed heterozygous missense variant in the U2AF2 gene (c.603G>T), who has a similar clinical phenotype as the patient reported before, including epilepsy, intellectual disability, language delay, microcephaly, and hypoplastic corpus callosum. We reviewed the phenotypic and genetic spectrum of patients with U2AF2-related neurological diseases, both newly diagnosed and previously reported. To investigate the possible pathogenesis, EBV-immortalized lymphoblastoid cells were derived from the peripheral blood obtained from the patient and control groups. Furthermore, according to the results of WB, RT-PCR, Q-PCR, and cDNA sequencing of RT-PCR products, the presumed missense variant c.603G>T caused exon 6 skipping in the U2AF2 mRNA transcript and led to a truncated protein (p.E163_E201del). Cell Counting Kit-8 (CCK-8) and cell cycle detection demonstrated that the variant c.603G>T inhibited the proliferation of patient lymphocyte cells compared with the control group. This study is aimed at expanding the phenotypic and genetic spectrum of U2AF2-related neurodevelopmental diseases and investigating the potential effects. This is the first report of the possible pathogenesis of a U2AF2 gene pathogenic variant in a patient with neurodevelopmental diseases and shows that a novel presumed missense variant in the U2AF2 gene causes exon skipping.

Prediction of complete remission and survival in acute myeloid leukemia using supervised machine learning.

Achievement of complete remission signifies a crucial milestone in the therapy of acute myeloid leukemia (AML) while refractory disease is associated with dismal outcomes. Hence, accurately identifying patients at risk is essential to tailor treatment concepts individually to disease biology. We used nine machine learning (ML) models to predict complete remission and 2-year overall survival in a large multicenter cohort of 1,383 AML patients who received intensive induction therapy. Clinical, laboratory, cytogenetic and molecular genetic data were incorporated and our results were validated on an external multicenter cohort. Our ML models autonomously selected predictive features including established markers of favorable or adverse risk as well as identifying markers of so-far controversial relevance. De novo AML, extramedullary AML, double-mutated CEBPA, mutations of CEBPA-bZIP, NPM1, FLT3-ITD, ASXL1, RUNX1, SF3B1, IKZF1, TP53, and U2AF1, t(8;21), inv(16)/t(16;16), del(5)/del(5q), del(17)/del(17p), normal or complex karyotypes, age and hemoglobin concentration at initial diagnosis were statistically significant markers predictive of complete remission, while t(8;21), del(5)/del(5q), inv(16)/t(16;16), del(17)/del(17p), double-mutated CEBPA, CEBPA-bZIP, NPM1, FLT3-ITD, DNMT3A, SF3B1, U2AF1, and TP53 mutations, age, white blood cell count, peripheral blast count, serum lactate dehydrogenase level and hemoglobin concentration at initial diagnosis as well as extramedullary manifestations were predictive for 2-year overall survival. For prediction of complete remission and 2-year overall survival areas under the receiver operating characteristic curves ranged between 0.77-0.86 and between 0.63-0.74, respectively in our test set, and between 0.71-0.80 and 0.65-0.75 in the external validation cohort. We demonstrated the feasibility of ML for risk stratification in AML as a model disease for hematologic neoplasms, using a scalable and reusable ML framework. Our study illustrates the clinical applicability of ML as a decision support system in hematology.