Activation of Cryptic 3' Splice-Sites by SRSF2 Contributes to Cassette Exon Skipping

Three exonic variants in the PHEX gene cause aberrant splicing in a minigene assay

Background: X-linked hypophosphatemia (XLH, OMIM 307800) is a rare phosphorus metabolism disorder caused by PHEX gene variants. Many variants simply classified as missense or nonsense variants were only analyzed at the DNA level. However, growing evidence indicates that some of these variants may alter pre-mRNA splicing, causing diseases. Therefore, this study aimed to use bioinformatics tools and a minigene assay to ascertain the effects of PHEX variations on pre-mRNA splicing. Methods: We analyzed 174 variants in the PHEX gene described as missense or nonsense variants. Finally, we selected eight candidate variants using bioinformatics tools to evaluate their effects on pre-mRNA splicing using a minigene assay system. The complementary DNA (cDNA) sequence for the PHEX gene (RefSeq NM_000444.6) serves as the basis for DNA variant numbering. Results: Of the eight candidate variants, three were found to cause abnormal splicing. Variants c.617T>G p.(Leu206Trp) and c.621T>A p.(Tyr207*) in exon 5 altered the splicing of pre-mRNA, owing to the activation of a cryptic splice site in exon 5, which produced an aberrant transcript lacking a part of exon 5, whereas variant c.1700G>C p.(Arg567Pro) in exon 16 led to the activation of a cryptic splice site in intron 16, resulting in a partial inclusion of intron 16. Conclusion: Our study employed a minigene system, which has a great degree of flexibility to assess abnormal splicing patterns under the circumstances of patient mRNA samples that are not available, to explore the impact of the exonic variants on pre-mRNA splicing. Based on the aforementioned experimental findings, we demonstrated the importance of analyzing exonic variants at the mRNA level.

SRSF2 plays an unexpected role as reader of m5C on mRNA, linking epitranscriptomics to cancer

A common mRNA modification is 5-methylcytosine (m5C), whose role in gene-transcript processing and cancer remains unclear. Here, we identify serine/arginine-rich splicing factor 2 (SRSF2) as a reader of m5C and impaired SRSF2 m5C binding as a potential contributor to leukemogenesis. Structurally, we identify residues involved in m5C recognition and the impact of the prevalent leukemia-associated mutation SRSF2P95H. We show that SRSF2 binding and m5C colocalize within transcripts. Furthermore, knocking down the m5C writer NSUN2 decreases mRNA m5C, reduces SRSF2 binding, and alters RNA splicing. We also show that the SRSF2P95H mutation impairs the ability of the protein to read m5C-marked mRNA, notably reducing its binding to key leukemia-related transcripts in leukemic cells. In leukemia patients, low NSUN2 expression leads to mRNA m5C hypomethylation and, combined with SRSF2P95H, predicts poor outcomes. Altogether, we highlight an unrecognized mechanistic link between epitranscriptomics and a key oncogenesis driver.

Identification and in-silico characterization of splice-site variants from a large cardiogenetic national registry

Splice-site variants in cardiac genes may predispose carriers to potentially lethal arrhythmias. To investigate, we screened 1315 probands and first-degree relatives enrolled in the Canadian Hearts in Rhythm Organization (HiRO) registry. 10% (134/1315) of patients in the HiRO registry carry variants within 10 base-pairs of the intron-exon boundary with 78% (104/134) otherwise genotype negative. These 134 probands were carriers of 57 unique variants. For each variant, American College of Medical Genetics and Genomics (ACMG) classification was revisited based on consensus between nine in silico tools. Due in part to the in silico algorithms, seven variants were reclassified from the original report, with the majority (6/7) downgraded. Our analyses predicted 53% (30/57) of variants to be likely/pathogenic. For the 57 variants, an average of 9 tools were able to score variants within splice sites, while 6.5 tools responded for variants outside these sites. With likely/pathogenic classification considered a positive outcome, the ACMG classification was used to calculate sensitivity/specificity of each tool. Among these, Combined Annotation Dependent Depletion (CADD) had good sensitivity (93%) and the highest response rate (131/134, 98%), dbscSNV was also sensitive (97%), and SpliceAI was the most specific (64%) tool. Splice variants remain an important consideration in gene elusive inherited arrhythmia syndromes. Screening for intronic variants, even when restricted to the ±10 positions as performed here may improve genetic testing yield. We compare 9 freely available in silico tools and provide recommendations regarding their predictive capabilities. Moreover, we highlight several novel cardiomyopathy-associated variants which merit further study.

HBV Enhances Sorafenib Resistance in Hepatocellular Carcinoma by Reducing Ferroptosis via SRSF2-Mediated Abnormal PCLAF Splicing

Hepatocellular carcinoma (HCC) is one of the most lethal human cancers. Hepatitis B virus (HBV) infection accounts for nearly 50% of HCC cases. Recent studies indicate that HBV infection induces resistance to sorafenib, the first-line systemic treatment for advanced HCC for more than a decade, from 2007 to 2020. Our previous research shows that variant 1 (tv1) of proliferating cell nuclear antigen clamp-associated factor (PCLAF), overexpressed in HCC, protects against doxorubicin-induced apoptosis. However, there are no reports on the relevance of PCLAF in sorafenib resistance in HBV-related HCC. In this article, we found that PCLAF levels were higher in HBV-related HCC than in non-virus-related HCC using bioinformatics analysis. Immunohistochemistry (IHC) staining of clinical samples and the splicing reporter minigene assay using HCC cells revealed that PCLAF tv1 was elevated by HBV. Furthermore, HBV promoted the splicing of PCLAF tv1 by downregulating serine/arginine-rich splicing factor 2 (SRSF2), which hindered the inclusion of PCLAF exon 3 through a putative cis-element (116-123), "GATTCCTG". The CCK-8 assay showed that HBV decreased cell susceptibility to sorafenib through SRSF2/PCLAF tv1. HBV reduced ferroptosis by decreasing intracellular Fe²⁺ levels and activating GPX4 expression via the SRSF2/PCLAF tv1 axis, according to a mechanism study. Suppressed ferroptosis, on the other hand, contributed to HBV-mediated sorafenib resistance through SRSF2/PCLAF tv1. These data suggested that HBV regulated PCLAF abnormal alternative splicing by suppressing SRSF2. HBV caused sorafenib resistance by reducing ferroptosis via the SRSF2/PCLAF tv1 axis. As a result, the SRSF2/PCLAF tv1 axis may be a prospective molecular therapeutic target in HBV-related HCC, as well as a predictor of sorafenib resistance. The inhibition of the SRSF2/PCLAF tv1 axis may be crucial in the emergence of systemic chemotherapy resistance in HBV-associated HCC.

Identification of novel deep intronic PAH gene variants in patients diagnosed with phenylketonuria

Phenylketonuria (PKU) is caused by phenylalanine hydroxylase (PAH) gene variants. Previously, 94.21% of variants were identified using Sanger sequencing and multiplex ligation-dependent probe amplification. To investigate the remaining variants, we performed whole-genome sequencing for four patients with PKU and unknown genotypes to identify deep intronic or structural variants. We identified three novel heterozygous variants (c.706+368T>C, c.1065+241C>A, and c.1199+502A>T) in a deep PAH gene intron. We detected a c.1199+502A>T variant in 60% (6/10) of PKU patients with genetically undetermined PKU. In silico predictions indicated that the three deep variants may impact splice site selection and result in the inclusion of a pseudo-exon. A c.1199+502A>T PAH minigene and reverse transcription PCR (RT-PCR) on blood RNA from a PKU patient with biallelic variants c.1199+502A>T and c.1199G>A confirmed that the c.1199+502A>T variant may strengthen the predicted branch point and leads to the inclusion of a 25-nt pseudo-exon in the PAH mRNA. Reverse transcription polymerase chain reaction (RT-PCR) on the minigene revealed that c.706+368T>C may create an SRSF2 (SC35) binding site via a 313-nt pseudo-exon, whereas c.1065+241C>A may produce an 81-nt pseudo-exon that strengthens the predicted SRSF1 (SF2/ASF) binding site. These results augment current knowledge of PAH genotypes and show that deep intronic analysis of PAH can genetically diagnose PKU.

The Impact of Pro-Inflammatory Cytokines on Alternative Splicing Patterns in Human Islets

Alternative splicing (AS) within the β cell has been proposed as one potential pathway that may exacerbate autoimmunity and unveil novel immunogenic epitopes in type 1 diabetes (T1D). We employed a computational strategy to prioritize pathogenic splicing events in human islets treated with IL-1β + IFN-γ as an ex vivo model of T1D and coupled this analysis with a k-mer based approach to predict RNA binding proteins involved in AS. In total, 969 AS events were identified in cytokine-treated islets, with the majority (44.8%) involving a skipped exon. ExonImpact identified 129 events predicted to impact protein structure. AS occurred with high frequency in MHC Class II-related mRNAs, and targeted qPCR validated reduced inclusion of Exon5 in the MHC Class II gene HLA-DMB. Single molecule RNA FISH confirmed increased HLA-DMB splicing in pancreatic sections from human donors with established T1D and autoantibody positivity. Serine and Arginine Rich Splicing Factor 2 was implicated in 37.2% of potentially pathogenic events, including Exon5 exclusion in HLA-DMB. Together, these data suggest that dynamic control of AS plays a role in the β cell response to inflammatory signals during T1D evolution.

Splicing factor SRSF2-centric gene regulation

Serine/arginine-rich splicing factor 2 (SRSF2) is a splicing factor that is widely expressed in a variety of mammalian cell types. Increasing evidence has confirmed that SRSF2 plays vital roles in a number of biological and pathological processes. Therefore, it is important to understand how its expression is regulated, and how it regulates the expression of its target genes. Recently, we found that SRSF2 expression could be upregulated by herpes simplex virus-1 (HSV-1) infection, and that altered SRSF2 expression, in turn, epigenetically regulates the transcription of HSV-1 genes. Further studies on T cell exhaustion demonstrated that upregulated SRSF2 in exhausted T cells elevated the levels of multiple immune checkpoint molecules by associating with the acyl-transferases, P300 and CBP, and by altering histone modification near the transcription start sites of these genes, thereby influencing signal transducer and activator of transcription 3 binding to these gene promoters. These findings suggest that SRSF2 acts as an important sensor and effector during disease progression. Here, we discuss the molecules that regulate SRSF2 gene expression and their associated mechanisms, and the mechanisms via which SRSF2 regulates the expression of target genes, thus providing novel insights into the central role of SRSF2 in gene regulation.

SRSF9 Regulates Cassette Exon Splicing of Caspase-2 by Interacting with Its Downstream Exon

Alternative splicing (AS) is an important posttranscriptional regulatory process. Damaged or unnecessary cells need to be removed though apoptosis to maintain physiological processes. Caspase-2 pre-mRNA produces pro-apoptotic long mRNA and anti-apoptotic short mRNA isoforms through AS. How AS of Caspase-2 is regulated remains unclear. In the present study, we identified a novel regulatory protein SRSF9 for AS of Caspase-2 cassette exon 9. Knock-down (KD) of SRSF9 increased inclusion of cassette exon and on the other hand, overexpression of SRSF9 decreased inclusion of this exon. Deletion mutagenesis demonstrated that exon 9, parts of intron 9, exon 8 and exon 10 were not required for the role of SRSF9 in Caspase-2 AS. However, deletion and substitution mutation analysis revealed that AGGAG sequence located at exon 10 provided functional target for SRSF9. In addition, RNA-pulldown mediated immunoblotting analysis showed that SRSF9 interacted with this sequence. Gene ontology analysis of RNA-seq from SRSF9 KD cells demonstrates that SRSF9 could regulate AS of a subset of apoptosis related genes. Collectively, our results reveal a basis for regulation of Caspase-2 AS.

U2AF65-Dependent SF3B1 Function in SMN Alternative Splicing

Splicing factor 3b subunit 1 (SF3B1) is an essential protein in spliceosomes and mutated frequently in many cancers. While roles of SF3B1 in single intron splicing and roles of its cancer-linked mutant in aberrant splicing have been identified to some extent, regulatory functions of wild-type SF3B1 in alternative splicing (AS) are not well-understood yet. Here, we applied RNA sequencing (RNA-seq) to analyze genome-wide AS in SF3B1 knockdown (KD) cells and to identify a large number of skipped exons (SEs), with a considerable number of alternative 5′ splice-site selection, alternative 3′ splice-site selection, mutually exclusive exons (MXE), and retention of introns (RI). Among altered SEs by SF3B1 KD, survival motor neuron 2 (SMN2) pre-mRNA exon 7 splicing was a regulatory target of SF3B1. RT-PCR analysis of SMN exon 7 splicing in SF3B1 KD or overexpressed HCT116, SH-SY5Y, HEK293T, and spinal muscular atrophy (SMA) patient cells validated the results. A deletion mutation demonstrated that the U2 snRNP auxiliary factor 65 kDa (U2AF65) interaction domain of SF3B1 was required for its function in SMN exon 7 splicing. In addition, mutations to lower the score of the polypyrimidine tract (PPT) of exon 7, resulting in lower affinity for U2AF65, were not able to support SF3B1 function, suggesting the importance of U2AF65 in SF3B1 function. Furthermore, the PPT of exon 7 with higher affinity to U2AF65 than exon 8 showed significantly stronger interactions with SF3B1. Collectively, our results revealed SF3B1 function in SMN alternative splicing.