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

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

Homozygous EPRS1 missense variant causing hypomyelinating leukodystrophy-15 alters variant-distal mRNA m6A site accessibility.

Hypomyelinating leukodystrophy (HLD) is an autosomal recessive disorder characterized by defective central nervous system myelination. Exome sequencing of two siblings with severe cognitive and motor impairment and progressive hypomyelination characteristic of HLD revealed homozygosity for a missense single-nucleotide variant (SNV) in EPRS1 (c.4444 C > A; p.Pro1482Thr), encoding glutamyl-prolyl-tRNA synthetase, consistent with HLD15. Patient lymphoblastoid cell lines express markedly reduced EPRS1 protein due to dual defects in nuclear export and cytoplasmic translation of variant EPRS1 mRNA. Variant mRNA exhibits reduced METTL3 methyltransferase-mediated writing of N6-methyladenosine (m6A) and reduced reading by YTHDC1 and YTHDF1/3 required for efficient mRNA nuclear export and translation, respectively. In contrast to current models, the variant does not alter the sequence of m6A target sites, but instead reduces their accessibility for modification. The defect was rescued by antisense morpholinos predicted to expose m6A sites on target EPRS1 mRNA, or by m6A modification of the mRNA by METTL3-dCas13b, a targeted RNA methylation editor. Our bioinformatic analysis predicts widespread occurrence of SNVs associated with human health and disease that similarly alter accessibility of distal mRNA m6A sites. These results reveal a new RNA-dependent etiologic mechanism by which SNVs can influence gene expression and disease, consequently generating opportunities for personalized, RNA-based therapeutics targeting these disorders.

m6A Reader YTHDC1 Impairs Respiratory Syncytial Virus Infection by Downregulating Membrane CX3CR1 Expression.

Respiratory syncytial virus (RSV) is the most prevalent cause of acute lower respiratory infection in young children. Currently, the first RSV vaccines are approved by the FDA. Recently, N6-methyladenosine (m6A) RNA methylation has been implicated in the regulation of the viral life cycle and replication of many viruses, including RSV. m6A methylation of RSV RNA has been demonstrated to promote replication and prevent anti-viral immune responses by the host. Whether m6A is also involved in viral entry and whether m6A can also affect RSV infection via different mechanisms than methylation of viral RNA is poorly understood. Here, we identify m6A reader YTH domain-containing protein 1 (YTHDC1) as a novel negative regulator of RSV infection. We demonstrate that YTHDC1 abrogates RSV infection by reducing the expression of RSV entry receptor CX3C motif chemokine receptor 1 (CX3CR1) on the cell surface of lung epithelial cells. Altogether, these data reveal a novel role for m6A methylation and YTHDC1 in the viral entry of RSV. These findings may contribute to the development of novel treatment options to control RSV infection.

KAP1 stabilizes MYCN mRNA and promotes neuroblastoma tumorigenicity by protecting the RNA m6A reader YTHDC1 protein degradation.

BACKGROUND: Neuroblastoma (NB) patients with amplified MYCN often face a grim prognosis and are resistant to existing therapies, yet MYCN protein is considered undruggable. KAP1 (also named TRIM28) plays a crucial role in multiple biological activities. This study aimed to investigate the relationship between KAP1 and MYCN in NB. METHODS: Transcriptome analyses and luciferase reporter assay identified that KAP1 was a downstream target of MYCN. The effects of KAP1 on cancer cell proliferation and colony formation were explored using the loss-of-function assays in vitro and in vivo. RNA stability detection was used to examine the influence of KAP1 on MYCN expression. The mechanisms of KAP1 to maintain MYCN mRNA stabilization were mainly investigated by mass spectrum, immunoprecipitation, RIP-qPCR, and western blotting. In addition, a xenograft mouse model was used to reveal the antitumor effect of STM2457 on NB. RESULTS: Here we identified KAP1 as a critical regulator of MYCN mRNA stability by protecting the RNA N6-methyladenosine (m6A) reader YTHDC1 protein degradation. KAP1 was highly expressed in clinical MYCN-amplified NB and was upregulated by MYCN. Reciprocally, KAP1 knockdown reduced MYCN mRNA stability and inhibited MYCN-amplified NB progression. Mechanistically, KAP1 regulated the stability of MYCN mRNA in an m6A-dependent manner. KAP1 formed a complex with YTHDC1 and RNA m6A writer METTL3 to regulate m6A-modified MYCN mRNA stability. KAP1 depletion decreased YTHDC1 protein stability and promoted MYCN mRNA degradation. Inhibiting MYCN mRNA m6A modification synergized with chemotherapy to restrain tumor progression in MYCN-amplified NB. CONCLUSIONS: Our research demonstrates that KAP1, transcriptionally activated by MYCN, forms a complex with YTHDC1 and METTL3, which in turn maintain the stabilization of MYCN mRNA in an m6A-dependent manner. Targeting m6A modification by STM2457, a small-molecule inhibitor of METTL3, could downregulate MYCN expression and attenuate tumor proliferation. This finding provides a new alternative putative therapeutic strategy for MYCN-amplified NB.

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.

Truncating the spliceosomal 'rope protein' Prp45 results in Htz1 dependent phenotypes.

Spliceosome assembly contributes an important but incompletely understood aspect of splicing regulation. Prp45 is a yeast splicing factor which runs as an extended fold through the spliceosome, and which may be important for bringing its components together. We performed a whole genome analysis of the genetic interaction network of the truncated allele of PRP45 (prp45(1-169)) using synthetic genetic array technology and found chromatin remodellers and modifiers as an enriched category. In agreement with related studies, H2A.Z-encoding HTZ1, and the components of SWR1, INO80, and SAGA complexes represented prominent interactors, with htz1 conferring the strongest growth defect. Because the truncation of Prp45 disproportionately affected low copy number transcripts of intron-containing genes, we prepared strains carrying intronless versions of SRB2, VPS75, or HRB1, the most affected cases with transcription-related function. Intron removal from SRB2, but not from the other genes, partly repaired some but not all the growth phenotypes identified in the genetic screen. The interaction of prp45(1-169) and htz1Δ was detectable even in cells with SRB2 intron deleted (srb2Δi). The less truncated variant, prp45(1-330), had a synthetic growth defect with htz1Δ at 16°C, which also persisted in the srb2Δi background. Moreover, htz1Δ enhanced prp45(1-330) dependent pre-mRNA hyper-accumulation of both high and low efficiency splicers, genes ECM33 and COF1, respectively. We conclude that while the expression defects of low expression intron-containing genes contribute to the genetic interactome of prp45(1-169), the genetic interactions between prp45 and htz1 alleles demonstrate the sensitivity of spliceosome assembly, delayed in prp45(1-169), to the chromatin environment.

Intramolecular autoinhibition regulates the selectivity of PRPF40A tandem WW domains for proline-rich motifs.

PRPF40A plays an important role in the regulation of pre-mRNA splicing by mediating protein-protein interactions in the early steps of spliceosome assembly. By binding to proteins at the 5´ and 3´ splice sites, PRPF40A promotes spliceosome assembly by bridging the recognition of the splices. The PRPF40A WW domains are expected to recognize proline-rich sequences in SF1 and SF3A1 in the early spliceosome complexes E and A, respectively. Here, we combine NMR, SAXS and ITC to determine the structure of the PRPF40A tandem WW domains in solution and characterize the binding specificity and mechanism for proline-rich motifs recognition. Our structure of the PRPF40A WW tandem in complex with a high-affinity SF1 peptide reveals contributions of both WW domains, which also enables tryptophan sandwiching by two proline residues in the ligand. Unexpectedly, a proline-rich motif in the N-terminal region of PRPF40A mediates intramolecular interactions with the WW tandem. Using NMR, ITC, mutational analysis in vitro, and immunoprecipitation experiments in cells, we show that the intramolecular interaction acts as an autoinhibitory filter for proof-reading of high-affinity proline-rich motifs in bona fide PRPF40A binding partners. We propose that similar autoinhibitory mechanisms are present in most WW tandem-containing proteins to enhance binding selectivity and regulation of WW/proline-rich peptide interaction networks.

Biological relevance of alternative splicing in hematologic malignancies.

Alternative splicing (AS) is a strictly regulated process that generates multiple mRNA variants from a single gene, thus contributing to proteome diversity. Transcriptome-wide sequencing studies revealed networks of functionally coordinated splicing events, which produce isoforms with distinct or even opposing functions. To date, several mechanisms of AS are deregulated in leukemic cells, mainly due to mutations in splicing and/or epigenetic regulators and altered expression of splicing factors (SFs). In this review, we discuss aberrant splicing events induced by mutations affecting SFs (SF3B1, U2AF1, SRSR2, and ZRSR2), spliceosome components (PRPF8, LUC7L2, DDX41, and HNRNPH1), and epigenetic modulators (IDH1 and IDH2). Finally, we provide an extensive overview of the biological relevance of aberrant isoforms of genes involved in the regulation of apoptosis (e. g. BCL-X, MCL-1, FAS, and c-FLIP), activation of key cellular signaling pathways (CASP8, MAP3K7, and NOTCH2), and cell metabolism (PKM).

Predicting survival in patients with myelodysplastic/myeloproliferative neoplasms with SF3B1 mutation and thrombocytosis.

We investigated data from 180 consecutive patients with myelodysplastic/myeloproliferative neoplasms with SF3B1 mutation and thrombocytosis (MDS/MPN-SF3B1-T) who were diagnosed according to the 2022 World Health Organization (WHO) classification of myeloid neoplasms to identify covariates associated with survival. At a median follow-up of 48 months (95% confidence interval [CI] 35-61 months), the median survival was 69 months (95% CI 59-79 months). Patients with bone marrow ring sideroblasts (RS) < 15% had shorter median overall survival (OS) than did those with bone marrow RS ≥ 15% (41 months [95% CI 32-50 months] versus 76 months [95% CI 59-93 months]; P < 0.001). According to the univariable analyses of OS, age ≥ 65 years (P < 0.001), hemoglobin concentration (Hb) < 80 g/L (P = 0.090), platelet count (PLT) ≥ 800 × 10E + 9/L (P = 0.087), bone marrow RS < 15% (P < 0.001), the Revised International Prognostic Scoring System (IPSS-R) cytogenetic category intermediate/poor/very poor (P = 0.005), SETBP1 mutation (P = 0.061) and SRSF2 mutation (P < 0.001) were associated with poor survival. Based on variables selected from univariable analyses, two separate survival prediction models, a clinical survival model, and a clinical-molecular survival model, were developed using multivariable analyses with the minimum value of the Akaike information criterion (AIC) to specifically predict outcomes in patients with MDS/MPN-SF3B1-T according to the 2022 WHO classification.

YTHDC1 promotes the malignant progression of gastric cancer by promoting ROD1 translocation to the nucleus.

RNA-binding proteins (RBPs) make vital impacts on tumor progression and are important potential targets for tumor treatment. Previous studies have shown that RBP regulator of differentiation 1 (ROD1), enriched in the nucleus, is abnormally expressed and functions as a splicing factor in tumors; however, the mechanism underlying its involvement in gastric cancer (GC) is unknown. In this study, ROD1 is found to stimulate GC cell proliferation and metastasis and is related to poor patient prognosis. In vitro experiments showed that ROD1 influences GC proliferation and metastasis through modulating the imbalance of the level of the oncogenic gene OIP5 and the tumor suppressor gene GPD1L. Further studies showed that the N6-methyladenosine (m6A) "reader" protein YTHDC1 can interact with ROD1 and regulate the balance of the expression of the downstream molecules OIP5/GPD1L by promoting the nuclear enrichment of ROD1. Therefore, YTHDC1 stimulates GC development and progression through modulating nuclear enrichment of the splicing factor ROD1.

RBM22 regulates RNA polymerase II 5' pausing, elongation rate, and termination by coordinating 7SK-P-TEFb complex and SPT5.

BACKGROUND: Splicing factors are vital for the regulation of RNA splicing, but some have also been implicated in regulating transcription. The underlying molecular mechanisms of their involvement in transcriptional processes remain poorly understood. RESULTS: Here, we describe a direct role of splicing factor RBM22 in coordinating multiple steps of RNA Polymerase II (RNAPII) transcription in human cells. The RBM22 protein widely occupies the RNAPII-transcribed gene locus in the nucleus. Loss of RBM22 promotes RNAPII pause release, reduces elongation velocity, and provokes transcriptional readthrough genome-wide, coupled with production of transcripts containing sequences from downstream of the gene. RBM22 preferentially binds to the hyperphosphorylated, transcriptionally engaged RNAPII and coordinates its dynamics by regulating the homeostasis of the 7SK-P-TEFb complex and the association between RNAPII and SPT5 at the chromatin level. CONCLUSIONS: Our results uncover the multifaceted role of RBM22 in orchestrating the transcriptional program of RNAPII and provide evidence implicating a splicing factor in both RNAPII elongation kinetics and termination control.

Accelerated DNA replication fork speed due to loss of R-loops in myelodysplastic syndromes with SF3B1 mutation.

Myelodysplastic syndromes (MDS) with mutated SF3B1 gene present features including a favourable outcome distinct from MDS with mutations in other splicing factor genes SRSF2 or U2AF1. Molecular bases of these divergences are poorly understood. Here we find that SF3B1-mutated MDS show reduced R-loop formation predominating in gene bodies associated with intron retention reduction, not found in U2AF1- or SRSF2-mutated MDS. Compared to erythroblasts from SRSF2- or U2AF1-mutated patients, SF3B1-mutated erythroblasts exhibit augmented DNA synthesis, accelerated replication forks, and single-stranded DNA exposure upon differentiation. Importantly, histone deacetylase inhibition using vorinostat restores R-loop formation, slows down DNA replication forks and improves SF3B1-mutated erythroblast differentiation. In conclusion, loss of R-loops with associated DNA replication stress represents a hallmark of SF3B1-mutated MDS ineffective erythropoiesis, which could be used as a therapeutic target.

Myelodysplastic Neoplasms (MDS) with Ring Sideroblasts or SF3B1 Mutations: The Improved Clinical Utility of World Health Organization and International Consensus Classification 2022 Definitions, a Single-Centre Retrospective Chart Review.

Myelodysplastic neoplasms (MDS) with ring sideroblasts (RS) are diagnosed via bone marrow aspiration in the presence of either (i) ≥15% RS or (ii) 5-14% RS and an SF3B1 mutation. In the MEDALIST trial and in an interim analysis of the COMMANDS trial, lower-risk MDS-RS patients had decreased transfusion dependency with luspatercept treatment. A total of 6817 patients with suspected hematologic malignancies underwent molecular testing using a next-generation-sequencing-based genetic assay and 395 MDS patients, seen at our centre from 1 January 2018 to 31 May 2023, were reviewed. Of these, we identified 39 evaluable patients as having lower-risk MDS with SF3B1 mutations: there were 20 (51.3%) males and 19 (48.7%) females, with a median age of 77 years (range of 57 to 92). Nineteen (48.7%) patients had an isolated SF3B1 mutation with a mean variant allele frequency of 35.2% +/- 8.1%, ranging from 7.4% to 46.0%. There were 29 (74.4%) patients with ≥15% RS, 6 (15.4%) with 5 to 14% RS, one (2.6%) with 1% RS, and 3 (7.7%) with no RS. Our study suggests that a quarter of patients would be missed based on the morphologic criterion of only using RS greater than 15% and supports the revised 2022 definitions of the World Health Organization (WHO) and International Consensus Classification (ICC), which shift toward molecularly defined subtypes of MDS and appropriate testing.

Association of poly(rC)-binding protein-2 with sideroflexin-3 through TOM20 as an iron entry pathway to mitochondria.

Iron is essential for all the lives and mitochondria integrate iron into heme and Fe-S clusters for diverse use as cofactors. Here, we screened mitochondrial proteins in KU812 human chronic myelogenous leukemia cells by glutathione S-transferase pulldown assay with PCBP2 to identify mitochondrial receptors for PCBP2, a major cytosolic Fe(II) chaperone. LC-MS analyses identified TOM20, sideroflexin-3 (SFXN3), SFXN1 and TOM70 in the affinity-score sequence. Stimulated emission depletion microscopy and proteinase-K digestion of mitochondria in HeLa cells revealed that TOM20 is located in the outer membrane of mitochondria whereas SFXN3 is located in the inner membrane. Although direct association was not observed between PCBP2 and SFXN3 with co-immunoprecipitation, proximity ligation assay demonstrated proximal localization of PCBP2 with TOM20 and there was a direct binding between TOM20 and SFXN3. Single knockdown either of PCBP2 and SFXN3 in K562 leukemia cells significantly decreased mitochondrial catalytic Fe(II) and mitochondrial maximal respiration. SFXN3 but not MFRN1 knockout (KO) in mouse embryonic fibroblasts decreased FBXL5 and heme oxygenase-1 (HO-1) but increased transferrin uptake and induced ferritin, indicating that mitochondrial iron entry through SFXN3 is distinct. MFRN1 KO revealed more intense mitochondrial Fe(II) deficiency than SFXN3 KO. Insufficient mitochondrial heme synthesis was evident under iron overload both with SFXN3 and MFRN KO, which was partially reversed by HO-1 inhibitor. Conversely, SFXN3 overexpression caused cytosolic iron deficiency with mitochondrial excess Fe(II), which further sensitized HeLa cells to RSL3-induced ferroptosis. In conclusion, we discovered a novel pathway of iron entry into mitochondria from cytosol through PCBP2-TOM20-SFXN3 axis.

GPATCH8 modulates mutant SF3B1 mis-splicing and pathogenicity in hematologic malignancies.

Mutations in the RNA splicing factor gene SF3B1 are common across hematologic and solid cancers and result in widespread alterations in splicing, yet there is currently no therapeutic means to correct this mis-splicing. Here, we utilize synthetic introns uniquely responsive to mutant SF3B1 to identify trans factors required for aberrant mutant SF3B1 splicing activity. This revealed the G-patch domain-containing protein GPATCH8 as required for mutant SF3B1-induced splicing alterations and impaired hematopoiesis. GPATCH8 is involved in quality control of branchpoint selection, interacts with the RNA helicase DHX15, and functionally opposes SURP and G-patch domain containing 1 (SUGP1), a G-patch protein recently implicated in SF3B1-mutant diseases. Silencing of GPATCH8 corrected one-third of mutant SF3B1-dependent splicing defects and was sufficient to improve dysfunctional hematopoiesis in SF3B1-mutant mice and primary human progenitors. These data identify GPATCH8 as a novel splicing factor required for mis-splicing by mutant SF3B1 and highlight the therapeutic impact of correcting aberrant splicing in SF3B1-mutant cancers.

A phase I study to evaluate the safety, pharmacokinetics, and pharmacodynamics of PF-06939999 (PRMT5 inhibitor) in patients with selected advanced or metastatic tumors with high incidence of splicing factor gene mutations.

BACKGROUND: Protein arginine methyltransferase 5 (PRMT5) methylates multiple substrates dysregulated in cancer, including spliceosome machinery components. PF-06939999 is a selective small-molecule PRMT5 inhibitor. PATIENTS AND METHODS: This phase I dose-escalation and -expansion trial (NCT03854227) enrolled patients with selected solid tumors. PF-06939999 was administered orally once or twice a day (q.d./b.i.d.) in 28-day cycles. The objectives were to evaluate PF-06939999 safety and tolerability to identify maximum tolerated dose (MTD) and recommended part 2 dose (RP2D), and assess pharmacokinetics (PK), pharmacodynamics [changes in plasma symmetric dimethylarginine (SDMA) levels], and antitumor activities. RESULTS: In part 1 dose escalation, 28 patients received PF-06939999 (0.5 mg q.d. to 6 mg b.i.d.). Four of 24 (17%) patients reported dose-limiting toxicities: thrombocytopenia (n = 2, 6 mg b.i.d.), anemia (n = 1, 8 mg q.d.), and neutropenia (n = 1, 6 mg q.d.). PF-06939999 exposure increased with dose. Steady-state PK was achieved by day 15. Plasma SDMA was reduced at steady state (58%-88%). Modulation of plasma SDMA was dose dependent. No MTD was determined. In part 2 dose expansion, 26 patients received PF-06939999 6 mg q.d. (RP2D). Overall (part 1 + part 2), the most common grade ≥3 treatment-related adverse events included anemia (28%), thrombocytopenia/platelet count decreased (22%), fatigue (6%), and neutropenia (4%). Three patients (6.8%) had confirmed partial response (head and neck squamous cell carcinoma, n = 1; non-small-cell lung cancer, n = 2), and 19 (43.2%) had stable disease. No predictive biomarkers were identified. CONCLUSIONS: PF-06939999 demonstrated a tolerable safety profile and objective clinical responses in a subset of patients, suggesting that PRMT5 is an interesting cancer target with clinical validation. However, no predictive biomarker was identified. The role of PRMT5 in cancer biology is complex and requires further preclinical, mechanistic investigation to identify predictive biomarkers for patient selection.

SF3B3-regulated mTOR alternative splicing promotes colorectal cancer progression and metastasis.

BACKGROUND: Aberrant alternative splicing (AS) is a pervasive event during colorectal cancer (CRC) development. SF3B3 is a splicing factor component of U2 small nuclear ribonucleoproteins which are crucial for early stages of spliceosome assembly. The role of SF3B3 in CRC remains unknown. METHODS: SF3B3 expression in human CRCs was analyzed using publicly available CRC datasets, immunohistochemistry, qRT-PCR, and western blot. RNA-seq, RNA immunoprecipitation, and lipidomics were performed in SF3B3 knockdown or overexpressing CRC cell lines. CRC cell xenografts, patient-derived xenografts, patient-derived organoids, and orthotopic metastasis mouse models were utilized to determine the in vivo role of SF3B3 in CRC progression and metastasis. RESULTS: SF3B3 was upregulated in CRC samples and associated with poor survival. Inhibition of SF3B3 by RNA silencing suppressed the proliferation and metastasis of CRC cells in vitro and in vivo, characterized by mitochondria injury, increased reactive oxygen species (ROS), and apoptosis. Mechanistically, silencing of SF3B3 increased mTOR exon-skipped splicing, leading to the suppression of lipogenesis via mTOR-SREBF1-FASN signaling. The combination of SF3B3 shRNAs and mTOR inhibitors showed synergistic antitumor activity in patient-derived CRC organoids and xenografts. Importantly, we identified SF3B3 as a critical regulator of mTOR splicing and autophagy in multiple cancers. CONCLUSIONS: Our findings revealed that SF3B3 promoted CRC progression and metastasis by regulating mTOR alternative splicing and SREBF1-FASN-mediated lipogenesis, providing strong evidence to support SF3B3 as a druggable target for CRC therapy.

Development and validation of AI/ML derived splice-switching oligonucleotides.

Splice-switching oligonucleotides (SSOs) are antisense compounds that act directly on pre-mRNA to modulate alternative splicing (AS). This study demonstrates the value that artificial intelligence/machine learning (AI/ML) provides for the identification of functional, verifiable, and therapeutic SSOs. We trained XGboost tree models using splicing factor (SF) pre-mRNA binding profiles and spliceosome assembly information to identify modulatory SSO binding sites on pre-mRNA. Using Shapley and out-of-bag analyses we also predicted the identity of specific SFs whose binding to pre-mRNA is blocked by SSOs. This step adds considerable transparency to AI/ML-driven drug discovery and informs biological insights useful in further validation steps. We applied this approach to previously established functional SSOs to retrospectively identify the SFs likely to regulate those events. We then took a prospective validation approach using a novel target in triple negative breast cancer (TNBC), NEDD4L exon 13 (NEDD4Le13). Targeting NEDD4Le13 with an AI/ML-designed SSO decreased the proliferative and migratory behavior of TNBC cells via downregulation of the TGFß pathway. Overall, this study illustrates the ability of AI/ML to extract actionable insights from RNA-seq data.

YTHDC1-dependent m6A modification modulated FOXM1 promotes glycolysis and tumor progression through CENPA in triple-negative breast cancer.

Triple-negative breast cancer (TNBC) exhibits heightened aggressiveness compared with other breast cancer (BC) subtypes, with earlier relapse, a higher risk of distant metastasis, and a worse prognosis. Transcription factors play a pivotal role in various cancers. Here, we found that factor forkhead box M1 (FOXM1) expression was significantly higher in TNBC than in other BC subtypes and normal tissues. Combining the findings of Gene Ontology (GO) enrichment analysis and a series of experiments, we found that knockdown of the FOXM1 gene attenuated the ability of TNBC cells to proliferate and metastasize both in vivo and in vitro. In addition, Spearman's test showed that FOXM1 significantly correlated with glycolysis-related genes, especially centromere protein A (CENPA) in datasets (GSE76250, GSE76124, GSE206912, and GSE103091). The effect of silencing FOXM1 on the inhibition of CENPA expression, TNBC proliferation, migration, and glycolysis could be recovered by overexpression of CENPA. According to MeRIP, the level of m6A modification on FOMX1 decreased in cells treated with cycloleucine (a m6A inhibitor) compared with that in the control group. The increase in FOXM1 expression caused by YTHDC1 overexpression could be reversed by the m6A inhibitor, which indicated that YTHDC1 enhanced FOXM1 expression depending on m6A modification. Therefore, we concluded that the YTHDC1-m6A modification/FOXM1/CENPA axis plays an important role in TNBC progression and glycolysis.