Alternative splicing of a chromatin modifier alters the transcriptional regulatory programs of stem cell maintenance and neuronal differentiation.

Development of embryonic stem cells (ESCs) into neurons requires intricate regulation of transcription, splicing, and translation, but how these processes interconnect is not understood. We found that polypyrimidine tract binding protein 1 (PTBP1) controls splicing of DPF2, a subunit of BRG1/BRM-associated factor (BAF) chromatin remodeling complexes. Dpf2 exon 7 splicing is inhibited by PTBP1 to produce the DPF2-S isoform early in development. During neuronal differentiation, loss of PTBP1 allows exon 7 inclusion and DPF2-L expression. Different cellular phenotypes and gene expression programs were induced by these alternative DPF2 isoforms. We identified chromatin binding sites enriched for each DPF2 isoform, as well as sites bound by both. In ESC, DPF2-S preferential sites were bound by pluripotency factors. In neuronal progenitors, DPF2-S sites were bound by nuclear factor I (NFI), while DPF2-L sites were bound by CCCTC-binding factor (CTCF). DPF2-S sites exhibited enhancer modifications, while DPF2-L sites showed promoter modifications. Thus, alternative splicing redirects BAF complex targeting to impact chromatin organization during neuronal development.

Alteration of RNA m6A methylation mediates aberrant RNA binding protein expression and alternative splicing in condyloma acuminatum.

Background: Condyloma acuminatum (CA) is caused by low-risk human papillomavirus, and is characterized by high recurrence after treatment. The RNA modification N6-methyladenosine (m6A) plays an important role during diverse viral infections, including high-risk HPV infection in cervical cancer. However, it is unclear whether low-risk HPV infection changes the RNA m6A methylation in CA. Methods: High-throughputm6A-sequencing was performed to profile the transcriptome-wide mRNA modifications of CA tissues infected by LR-HPVs and the paired normal tissues from CA patients. We further investigated the regulation of alternative splicing by RNA binding proteins (RBPs) with altered m6A modification and constructed a regulatory network among these RBPs, regulated alternative splicing events (RASEs) and regulated alternative splicing genes (RASGs) in CA. Results: The results show that the m6A level in CA tissues differed from that in the paired controls. Furthermore, cell cycle- and cell adhesion- associated genes with m6A modification were differentially expressed in CA tissues compared to the paired controls. In particular, seven RNA binding protein genes with specific m6A methylated sites, showed a higher or lower expression at the mRNA level in CA tissues than in the paired normal tissues. In addition, these differentially expressed RNA binding protein genes would regulate the alternative splicing pattern of apoptotic process genes in CA tissue. Conclusions: Our study reveals a sophisticated m6A modification profile in CA tissue that affects the response of host cells to HPV infection, and provides cues for the further exploration of the roles of m6A and the development of a novel treatment strategy for CA.

The Regulatory Network of hnRNPs Underlying Regulating PKM Alternative Splicing in Tumor Progression.

One of the hallmarks of cancer is metabolic reprogramming in tumor cells, and aerobic glycolysis is the primary mechanism by which glucose is quickly transformed into lactate. As one of the primary rate-limiting enzymes, pyruvate kinase (PK) M is engaged in the last phase of aerobic glycolysis. Alternative splicing is a crucial mechanism for protein diversity, and it promotes PKM precursor mRNA splicing to produce PKM2 dominance, resulting in low PKM1 expression. Specific splicing isoforms are produced in various tissues or illness situations, and the post-translational modifications are linked to numerous disorders, including cancers. hnRNPs are one of the main components of the splicing factor families. However, there have been no comprehensive studies on hnRNPs regulating PKM alternative splicing. Therefore, this review focuses on the regulatory network of hnRNPs on PKM pre-mRNA alternative splicing in tumors and clinical drug research. We elucidate the role of alternative splicing in tumor progression, prognosis, and the potential mechanism of abnormal RNA splicing. We also summarize the drug targets retarding tumorous splicing events, which may be critical to improving the specificity and effectiveness of current therapeutic interventions.

Assessing the Impact of Novel BRCA1 Exon 11 Variants on Pre-mRNA Splicing.

Our study focused on assessing the effects of three newly identified BRCA1 exon 11 variants (c.1019T>C, c.2363T>G, and c.3192T>C) on breast cancer susceptibility. Using computational predictions and experimental splicing assays, we evaluated their potential as pathogenic mutations. Our in silico analyses suggested that the c.2363T>G and c.3192T>C variants could impact both splicing and protein function, resulting in the V340A and V788G mutations, respectively. We further examined their splicing effects using minigene assays in MCF7 and SKBR3 breast cancer cell lines. Interestingly, we found that the c.2363T>G variant significantly altered splicing patterns in MCF7 cells but not in SKBR3 cells. This finding suggests a potential influence of cellular context on the variant's effects. While attempts to correlate in silico predictions with RNA binding factors were inconclusive, this observation underscores the complexity of splicing regulation. Splicing is governed by various factors, including cellular contexts and protein interactions, making it challenging to predict outcomes accurately. Further research is needed to fully understand the functional consequences of the c.2363T>G variant in breast cancer pathogenesis. Integrating computational predictions with experimental data will provide valuable insights into the role of alternative splicing regulation in different breast cancer types and stages.

SR proteins in cancer: function, regulation, and small inhibitor.

Alternative splicing of pre-mRNAs is a fundamental step in RNA processing required for gene expression in most metazoans. Serine and arginine-rich proteins (SR proteins) comprise a family of multifunctional proteins that contain an RNA recognition motif (RRM) and the ultra-conserved arginine/serine-rich (RS) domain, and play an important role in precise alternative splicing. Increasing research supports SR proteins as also functioning in other RNA-processing-related mechanisms, such as polyadenylation, degradation, and translation. In addition, SR proteins interact with N6-methyladenosine (m6A) regulators to modulate the methylation of ncRNA and mRNA. Dysregulation of SR proteins causes the disruption of cell differentiation and contributes to cancer progression. Here, we review the distinct biological characteristics of SR proteins and their known functional mechanisms during carcinogenesis. We also summarize the current inhibitors that directly target SR proteins and could ultimately turn SR proteins into actionable therapeutic targets in cancer therapy.

Genomic deletions explain the generation of alternative BRAF isoforms conferring resistance to MAPK inhibitors in melanoma.

Resistance to MAPK inhibitors (MAPKi), the main cause of relapse in BRAF-mutant melanoma, is associated with the production of alternative BRAF mRNA isoforms (altBRAFs) in up to 30% of patients receiving BRAF inhibitor monotherapy. These altBRAFs have been described as being generated by alternative pre-mRNA splicing, and splicing modulation has been proposed as a therapeutic strategy to overcome resistance. In contrast, we report that altBRAFs are generated through genomic deletions. Using different in vitro models of altBRAF-mediated melanoma resistance, we demonstrate the production of altBRAFs exclusively from the BRAF V600E allele, correlating with corresponding genomic deletions. Genomic deletions are also detected in tumor samples from melanoma and breast cancer patients expressing altBRAFs. Along with the identification of altBRAFs in BRAF wild-type and in MAPKi-naive melanoma samples, our results represent a major shift in our understanding of mechanisms leading to the generation of BRAF transcripts variants associated with resistance in melanoma.

Antisense oligonucleotide therapeutic approach for Timothy syndrome.

Timothy syndrome (TS) is a severe, multisystem disorder characterized by autism, epilepsy, long-QT syndrome and other neuropsychiatric conditions1. TS type 1 (TS1) is caused by a gain-of-function variant in the alternatively spliced and developmentally enriched CACNA1C exon 8A, as opposed to its counterpart exon 8. We previously uncovered several phenotypes in neurons derived from patients with TS1, including delayed channel inactivation, prolonged depolarization-induced calcium rise, impaired interneuron migration, activity-dependent dendrite retraction and an unanticipated persistent expression of exon 8A2-6. We reasoned that switching CACNA1C exon utilization from 8A to 8 would represent a potential therapeutic strategy. Here we developed antisense oligonucleotides (ASOs) to effectively decrease the inclusion of exon 8A in human cells both in vitro and, following transplantation, in vivo. We discovered that the ASO-mediated switch from exon 8A to 8 robustly rescued defects in patient-derived cortical organoids and migration in forebrain assembloids. Leveraging a transplantation platform previously developed7, we found that a single intrathecal ASO administration rescued calcium changes and in vivo dendrite retraction of patient neurons, suggesting that suppression of CACNA1C exon 8A expression is a potential treatment for TS1. Broadly, these experiments illustrate how a multilevel, in vivo and in vitro stem cell model-based approach can identify strategies to reverse disease-relevant neural pathophysiology.

Alternative splicing in prostate cancer progression and therapeutic resistance.

Prostate cancer (CaP) remains the second leading cause of cancer deaths in western men. CaP mortality results from diverse molecular mechanisms that mediate resistance to the standard of care treatments for metastatic disease. Recently, alternative splicing has been recognized as a hallmark of CaP aggressiveness. Alternative splicing events cause treatment resistance and aggressive CaP behavior and are determinants of the emergence of the two major types of late-stage treatment-resistant CaP, namely castration-resistant CaP (CRPC) and neuroendocrine CaP (NEPC). Here, we review recent multi-omics data that are uncovering the complicated landscape of alternative splicing events during CaP progression and the impact that different gene transcript isoforms can have on CaP cell biology and behavior. We discuss renewed insights in the molecular machinery by which alternative splicing occurs and contributes to the failure of systemic CaP therapies. The potential for alternative splicing events to serve as diagnostic markers and/or therapeutic targets is explored. We conclude by considering current challenges and promises associated with splicing-modulating therapies, and their potential for clinical translation into CaP patient care.

Activation of Cryptic Donor Splice Sites within the UDP-Glucuronosyltransferase (UGT)1A First-Exon Region Generates Variant Transcripts That Encode UGT1A Proteins with Truncated Aglycone-Binding Domains.

The human UDP-glucuronosyltransferases (UGTs) have crucial roles in metabolizing and clearing numerous small lipophilic compounds. The UGT1A locus generates nine UGT1A mRNAs, 65 spliced transcripts, and 34 circular RNAs. In this study, our analysis of published UGT-RNA capture sequencing (CaptureSeq) datasets identified novel splice junctions that predict 24 variant UGT1A transcripts derived from ligation of exon 2 to unique sequences within the UGT1A first-exon region using cryptic donor splice sites. Of these variants, seven (1A1_n1, 1A3_n3, 1A4_n4, 1A5_n1, 1A8_n2, 1A9_n2, 1A10_n7) are predicted to encode UGT1A proteins with truncated aglycone-binding domains. We assessed their expression profiles and deregulation in cancer using four RNA sequencing (RNA-Seq) datasets of paired normal and cancerous drug-metabolizing tissues from large patient cohorts. Variants were generally coexpressed with their canonical counterparts with a higher relative abundance in tumor than in normal tissues. Variants showed tissue-specific expression with high interindividual variability but overall low abundance. However, 1A8_n2 showed high abundance in normal and cancerous colorectal tissues, with levels that approached or surpassed canonical 1A8 mRNA levels in many samples. We cloned 1A8_n2 and showed expression of the predicted protein (1A8_i3) in human embryonic kidney (HEK)293T cells. Glucuronidation assays with 4-methylumbelliferone (4MU) showed that 1A8_i3 had no activity and was unable to inhibit the activity of 1A8_i1 protein. In summary, the activation of cryptic donor splice sites within the UGT1A first-exon region expands the UGT1A transcriptome and proteome. The 1A8_n2 cryptic donor splice site is highly active in colorectal tissues, representing an important cis-regulatory element that negatively regulates the function of the UGT1A8 gene through pre-mRNA splicing. SIGNIFICANT STATEMENT: The UGT1A locus generates nine canonical mRNAs, 65 alternately spliced transcripts, and 34 different circular RNAs. The present study reports a series of novel UDP-glucuronosyltransferase (UGT)1A variants resulting from use of cryptic donor splice sites in both normal and cancerous tissues, several of which are predicted to encode variant UGT1A proteins with truncated aglycone-binding domains. Of these, 1A8_n2 shows exceptionally high abundance in colorectal tissues, highlighting its potential role in the first-pass metabolism in gut through the glucuronidation pathway.

The splicing factor WBP11 mediates MCM7 intron retention to promote the malignant progression of ovarian cancer.

Accumulating studies suggest that splicing factors play important roles in many diseases including human cancers. Our study revealed that WBP11, a core splicing factor, is highly expressed in ovarian cancer (OC) tissues and associated with a poor prognosis. WBP11 inhibition significantly impaired the proliferation and mobility of ovarian cancer cells in vitro and in vivo. Furthermore, FOXM1 transcriptionally activated WBP11 expression by directly binding to its promoter in OC cells. Importantly, RNA-seq and alternative splicing event analysis revealed that WBP11 silencing decreased the expression of MCM7 by regulating intron 4 retention. MCM7 inhibition attenuated the increase in malignant behaviors of WBP11-overexpressing OC cells. Overall, WBP11 was identified as an oncogenic splicing factor that contributes to malignant progression by repressing intron 4 retention of MCM7 in OC cells. Thus, WBP11 is an oncogenic splicing factor with potential therapeutic and prognostic implications in OC.

Deciphering alternative splicing events and their therapeutic implications in colorectal Cancer.

Colorectal cancer (CRC) is one of the most common malignant tumors with complex molecular regulatory mechanisms. Alternative splicing (AS), a fundamental regulatory process of gene expression, plays an important role in the occurrence and development of CRC. This study analyzed AS Percent Spliced In (PSI) values from 49 pairs of CRC and normal samples in the TCGA SpliceSeq database. Using Lasso and SVM, AS features that can differentiate colorectal cancer from normal were screened. Univariate COX regression analysis identified prognosis-related AS events. A risk model was constructed and validated using machine learning, Kaplan-Meier analysis, and Decision Curve Analysis. The regulatory effect of protein arginine methyltransferase 5 (PRMT5) on poly(RC) binding protein 1 (PCBP1) was verified by immunoprecipitation experiments, and the effect of PCBP1 on the AS of Obscurin (OBSCN) was verified by PCR. Five AS events, including HNF4A.59461.AP and HNF4A.59462.AP, were identified, which can distinguish CRC from normal tissue. A machine learning model using 21 key AS events accurately predicted CRC prognosis. High-risk patients had significantly shorter survival times. PRMT5 was found to regulate PCBP1 function and then influence OBSCN AS, which may drive CRC progression. The study concluded that some AS events is significantly different in CRC and normal tissues, and some of these AS events are related to the prognosis of CRC. In addition, PRMT family-driven arginine modifications play an important role in CRC-specific AS events.

Alternative splicing of CsWRKY21 positively regulates cold response in tea plant.

Alternative splicing (AS) was an important post-transcriptional mechanism that involved in plant resistance to adversity stress. WRKY transcription factors function as transcriptional activators or repressors to modulate plant growth, development and stress response. However, the role of alternate splicing of WRKY in cold tolerance is poorly understood in tea plants. In this study, we found that the CsWRKY21 transcription factor, a member of the WRKY IId subfamily, was induced by low temperature. Subcellular localization and transcriptional activity assays showed that CsWRKY21 localized to the nucleus and had no transcriptional activation activity. Y1H and dual-luciferase reporter assays showed that CsWRKY21 suppressed expression of CsABA8H and CsUGT by binding with their promoters. Transient overexpression of CsABA8H and CsUGT reduced abscisic acid (ABA) content in tobacco leaves. Furthermore, we discovered that CsWRKY21 undergoes AS in the 5'UTR region. The AS transcript CsWRKY21-b was induced at low temperature, up to 6 folds compared to the control, while the full-length CsWRKY21-a transcript did not significantly change. Western blot analysis showed that the retention of introns in the 5'UTR region of CsWRKY21-b led to higher CsWRKY21 protein content. These results revealed that alternative splicing of CsWRKY21 involved in cold tolerance of tea plant by regulating the protein expression level and then regulating the content of ABA, and provide insights into molecular mechanisms of low temperature defense mediated by AS in tea plant.

RNA splicing regulator EIF3D regulates the tumor microenvironment through immunogene-related alternative splicing in head and neck squamous cell carcinoma.

Study finds that eukaryotic translation initiation factor 3 subunit D (EIF3D) may play an important role in aberrant alternative splicing (AS) events in tumors. AS possesses a pivotal role in both tumour progression and the constitution of the tumour microenvironment (TME). Regrettably, our current understanding of AS remains circumscribed especially in the context of immunogene-related alternative splicing (IGAS) profiles within Head and Neck Squamous Cell Carcinoma (HNSC). In this study, we comprehensively analyzed the function and mechanism of action of EIF3D by bioinformatics analysis combined with in vitro cellular experiments, and found that high expression of EIF3D in HNSC was associated with poor prognosis of overall survival (OS) and progression-free survival (PFS). The EIF3D low expression group had a higher degree of immune infiltration and better efficacy against PD1 and CTLA4 immunotherapy compared to the EIF3D high expression group. TCGA SpliceSeq analysis illustrated that EIF3D influenced differentially spliced alternative splicing (DSAS) events involving 105 differentially expressed immunogenes (DEIGs). We observed an induction of apoptosis and a suppression of cell proliferation, migration, and invasion in EIF3D knock-down FaDu cells. RNA-seq analysis unveiled that 531 genes exhibited differential expression following EIF3D knockdown in FaDu cells. These include 52 DEIGs. Furthermore, EIF3D knockdown influenced the patterns of 1923 alternative splicing events (ASEs), encompassing 129 IGASs. This study identified an RNA splicing regulator and revealed its regulatory role in IGAS and the TME of HNSC, suggesting that EIF3D may be a potential target for predicting HNSC prognosis and immunotherapeutic response.

Multifaceted roles of PDCD6 both within and outside the cell.

Programmed cell death protein 6 (PDCD6) is an evolutionarily conserved Ca2+-binding protein. PDCD6 is involved in regulating multifaceted and pleiotropic cellular processes in different cellular compartments. For instance, nuclear PDCD6 regulates apoptosis and alternative splicing. PDCD6 is required for coat protein complex II-dependent endoplasmic reticulum-to-Golgi apparatus vesicular transport in the cytoplasm. Recent advances suggest that cytoplasmic PDCD6 is involved in the regulation of cytoskeletal dynamics and innate immune responses. Additionally, membranous PDCD6 participates in membrane repair through endosomal sorting complex required for transport complex-dependent membrane budding. Interestingly, extracellular vesicles are rich in PDCD6. Moreover, abnormal expression of PDCD6 is closely associated with many diseases, especially cancer. PDCD6 is therefore a multifaceted but pivotal protein in vivo. To gain a more comprehensive understanding of PDCD6 functions and to focus and stimulate PDCD6 research, this review summarizes key developments in its role in different subcellular compartments, processes, and pathologies.

A common cellular response to broad splicing perturbations is characterized by metabolic transcript downregulation driven by the Mdm2-p53 axis.

Disruptions in core cellular processes elicit stress responses that drive cell-state changes leading to organismal phenotypes. Perturbations in the splicing machinery cause widespread mis-splicing, resulting in p53-dependent cell-state changes that give rise to cell-type-specific phenotypes and disease. However, a unified framework for how cells respond to splicing perturbations, and how this response manifests itself in nuanced disease phenotypes, has yet to be established. Here, we show that a p53-stabilizing Mdm2 alternative splicing event and the resulting widespread downregulation of metabolic transcripts are common events that arise in response to various splicing perturbations in both cellular and organismal models. Together, our results classify a common cellular response to splicing perturbations, put forth a new mechanism behind the cell-type-specific phenotypes that arise when splicing is broadly disrupted, and lend insight into the pleiotropic nature of the effects of p53 stabilization in disease.

Identification of Nonfunctional Alternatively Spliced Isoforms of STING in Human Acute Myeloid Leukemia.

Lack of robust activation of Stimulator of Interferon Genes (STING) pathway and subsequent induction of type I IFN responses is considered a barrier to antitumor immunity in acute myeloid leukemia (AML). Using common human AML cell lines as in vitro tools to evaluate the efficacy of novel STING agonists, we found most AML lines to be poor producers of IFNs upon exposure to extremely potent agonists, suggesting cell-intrinsic suppression of STING signaling may occur. We observed unexpected patterns of response that did not correlate with levels of STING pathway components or of known enzymes associated with resistance. To identify a genetic basis for these observations, we cloned and sequenced STING from the cDNA of human AML cell lines and found both frequent mutations and deviations from normal RNA splicing. We identified two novel spliced isoforms of STING in these lines and validated their expression in primary human AML samples. When transduced into reporter cells, these novel STING isoforms exhibited complete insensitivity to agonist stimulation. These observations identify alternative splicing as a mechanism of STING pathway suppression and suggest that most AML silences the STING pathway through direct modification rather than through engagement of external inhibitory factors. SIGNIFICANCE: We find that AML acquires resistance to innate immune activation via the STING pathway through aberrant splicing of the STING transcript including two novel forms described herein that act as dominant negatives. These data broaden understanding of how cancers evolve STING resistance, and suggest that the AML tumor microenvironment, not the cancer cell, should be the target of therapeutic interventions to activate STING.

Interrogations of single-cell RNA splicing landscapes with SCASL define new cell identities with physiological relevance.

RNA splicing shapes the gene regulatory programs that underlie various physiological and disease processes. Here, we present the SCASL (single-cell clustering based on alternative splicing landscapes) method for interrogating the heterogeneity of RNA splicing with single-cell RNA-seq data. SCASL resolves the issue of biased and sparse data coverage on single-cell RNA splicing and provides a new scheme for classifications of cell identities. With previously published datasets as examples, SCASL identifies new cell clusters indicating potentially precancerous and early-tumor stages in triple-negative breast cancer, illustrates cell lineages of embryonic liver development, and provides fine clusters of highly heterogeneous tumor-associated CD4 and CD8 T cells with functional and physiological relevance. Most of these findings are not readily available via conventional cell clustering based on single-cell gene expression data. Our study shows the potential of SCASL in revealing the intrinsic RNA splicing heterogeneity and generating biological insights into the dynamic and functional cell landscapes in complex tissues.

Detecting and understanding meaningful cancerous mutations based on computational models of mRNA splicing.

Cancer research has long relied on non-silent mutations. Yet, it has become overwhelmingly clear that silent mutations can affect gene expression and cancer cell fitness. One fundamental mechanism that apparently silent mutations can severely disrupt is alternative splicing. Here we introduce Oncosplice, a tool that scores mutations based on models of proteomes generated using aberrant splicing predictions. Oncosplice leverages a highly accurate neural network that predicts splice sites within arbitrary mRNA sequences, a greedy transcript constructor that considers alternate arrangements of splicing blueprints, and an algorithm that grades the functional divergence between proteins based on evolutionary conservation. By applying this tool to 12M somatic mutations we identify 8K deleterious variants that are significantly depleted within the healthy population; we demonstrate the tool's ability to identify clinically validated pathogenic variants with a positive predictive value of 94%; we show strong enrichment of predicted deleterious mutations across pan-cancer drivers. We also achieve improved patient survival estimation using a proposed set of novel cancer-involved genes. Ultimately, this pipeline enables accelerated insight-gathering of sequence-specific consequences for a class of understudied mutations and provides an efficient way of filtering through massive variant datasets - functionalities with immediate experimental and clinical applications.

Satellite double-stranded RNA induces mesenchymal transition in pancreatic cancer by regulating alternative splicing.

Human satellite II (HSATII), composed of tandem repeats in pericentromeric regions, is aberrantly transcribed in epithelial cancers, particularly pancreatic cancer. Dysregulation of repetitive elements in cancer tissues can facilitate incidental dsRNA formation; however, it remains controversial whether dsRNAs play tumor-promoting or tumor-suppressing roles during cancer progression. Therefore, we focused on the double-stranded formation of HSATII RNA and explored its molecular function. The overexpression of double-stranded HSATII (dsHSATII) RNA promoted mesenchymal-like morphological changes and enhanced the invasiveness of pancreatic cancer cells. We identified an RNA-binding protein, spermatid perinuclear RNA-binding protein (STRBP), which preferentially binds to dsHSATII RNA rather than single-stranded HSATII RNA. The mesenchymal transition of dsHSATII-expressing cells was rescued by STRBP overexpression. Mechanistically, STRBP is involved in the alternative splicing of genes associated with epithelial-mesenchymal transition (EMT). We also confirmed that isoform switching of CLSTN1, driven by dsHSATII overexpression or STRBP depletion, induced EMT-like morphological changes. These findings reveal a novel tumor-promoting function of dsHSATII RNA, inducing EMT-like changes and cell invasiveness, thus enhancing our understanding of the biological significance of aberrant expression of satellite arrays in malignant tumors.

Exogenous hydrogen sulphide promotes plant flowering through the Arabidopsis splicing factor AtU2AF65a.

Alternative splicing (AS) is an important regulatory mode at the post-transcriptional level, through which many flowering genes regulate floral transition by producing multiple transcripts, and splicing factors have essential roles in this process. Hydrogen sulphide (H2S) is a newly found gasotransmitter that has critical physiological roles in plants, and one of its potential modes of action is via persulfidation of target proteins at specific cysteine sites. Previously, it has been shown that both the splicing factor AtU2AF65a and H2S are involved in the regulation of plant flowering. This study found that, in Arabidopsis, the promoting effect of H2S on flowering was abolished in atu2af65a-4 mutants. Transcriptome analyses showed that when AtU2AF65a contained mutations, the regulatory function of H2S during the AS of many flowering genes (including SPA1, LUH, LUG and MAF3) was inhibited. The persulfidation assay showed that AtU2AF65a can be persulfidated by H2S, and the RNA immunoprecipitation data indicated that H2S could alter the binding affinity of AtU2AF65a to the precursor messenger RNA of the above-mentioned flowering genes. Overall, our results suggest that H2S may regulate the AS of flowering-related genes through persulfidation of splicing factor AtU2AF65a and thus lead to early flowering in plants.