Search

In spliceosome assembly, the 5 splice site is initially recognized by U1 snRNA. U1 leaves the spliceosome during the assembly process, therefore other factors contribute to the maintenance of 5 splice site identity as it is loaded into the catalytic site. Recent structural data suggest that human tri-snRNP 27K (SNRP27) M141 and SNU66 H734 interact to stabilize the U4/U6 quasi-pseudo knot at the base of the U6 snRNA ACAGAGA box in pre-B complex. Previously, we found that mutations in Caenorhabditis elegans at SNRP-27 M141 promote changes in alternative 5ss usage. We tested whether the potential interaction between SNRP-27 M141 and SNU-66 H765 (the C. elegans equivalent position to human SNU66 H734) contributes to maintaining 5 splice site identity during spliceosome assembly. We find that SNU-66 H765 mutants promote alternative 5 splice site usage. Many of the alternative 5 splicing events affected by SNU-66(H765G) overlap with those affected SNRP-27(M141T). Double mutants of snrp-27(M141T) and snu-66(H765G) are homozygous lethal. We hypothesize that mutations at either SNRP-27 M141 or SNU-66 H765 allow the spliceosome to load alternative 5 splice sites into the active site. Tests with mutant U1 snRNA and swapped 5 splice sites indicate that the ability of SNRP-27 M141 and SNU-66 H765 mutants to affect a particular 5 splice alternative splicing event is dependent on both the presence of a weaker consensus 5ss nearby and potentially nearby splicing factor binding sites. Our findings confirm a new role for the C terminus of SNU-66 in maintenance of 5 splice site identity during spliceosome assembly.

RNA helicases drive necessary rearrangements and ensure fidelity during the pre-mRNA splicing cycle. DEAD-box helicase DDX41 has been linked to human disease and has recently been shown to interact with DEAH-box helicase PRP22 in the spliceosomal C* complex, yet its function in splicing remains unknown. Depletion of DDX41 homolog SACY-1 from somatic cells has been previously shown to lead to changes in alternative 3 splice site (3ss) usage. Here, we show by transcriptomic analysis of published and novel data sets that SACY-1 perturbation causes a previously unreported pattern in alternative 3 splicing in introns with pairs of 3 splice sites ≤18 nt away from each other. We find that both SACY-1 depletion and the allele sacy-1(G533R) lead to a striking unidirectional increase in the usage of the proximal (upstream) 3ss. We previously discovered a similar alternative splicing pattern between germline tissue and somatic tissue, in which there is a unidirectional increase in proximal 3ss usage in the germline for ∼200 events; many of the somatic SACY-1 alternative 3 splicing events overlap with these developmentally regulated events. We generated targeted mutant alleles of the Caenorhabditis elegans homolog of PRP22, mog-5, in the region of MOG-5 that is predicted to interact with SACY-1 based on the human C* structure. These viable alleles, and a mimic of the myelodysplastic syndrome-associated allele DDX41(R525H), all promote the usage of proximal alternative adjacent 3 splice sites. We show that PRP22/MOG-5 and DDX41/SACY-1 have overlapping roles in proofreading the 3ss.

Branching is a critical step in RNA splicing that is essential for 5 splice site selection. Recent spliceosome structures have led to competing models for the recognition of the invariant adenosine at the branch point. However, there are no structures of any splicing complex with the adenosine nucleophile docked in the active site and positioned to attack the 5 splice site. Thus we lack a mechanistic understanding of adenosine selection and splice site recognition during RNA splicing. Here we present a cryo-electron microscopy structure of a group II intron that reveals that active site dynamics are coupled to the formation of a base triple within the branch-site helix that positions the 2-OH of the adenosine for nucleophilic attack on the 5 scissile phosphate. This structure, complemented with biochemistry and comparative analyses to splicing complexes, supports a base triple model of adenosine recognition for branching within group II introns and the evolutionarily related spliceosome.

The DYNC2H1 gene has been associated with short‐rib polydactyly syndrome (SRPS), among other skeletal ciliopathies. Two cases are presented of distinctive phenotypes resulting from splicing variants in DYNC2H1. The first is a 14‐week‐old fetus with enlarged nuchal translucency, oral hamartoma, malformed uvula, bifid epiglottis, short ribs, micromelia, long bone agenesis, polysyndactyly, heart defect, pancreatic cysts, multicystic dysplastic kidney, megabladder and trident acetabulum. A ciliopathies NGS panel revealed two compound heterozygous variants in DYNC2H1: c.7840‐18T>G r.7841_7964del p.Gly2614Aspfs*5 and c.11070G>A r.11044_11116del p.Ile3682Aspfs*2. Both variants were initially classified as variants of uncertain significance but were reclassified as likely pathogenic after PCR‐based RNA testing. The second is an 11‐year‐old overweight male with multiple accessory oral frenula, median cleft lip and alveolar ridge, polysyndactyly, brachydactyly, normal rib length, and hypogonadism. Exome sequencing revealed two compound heterozygous variants in DYNC2H1: c.6315del p.(Thr2106Glnfs*7), classified as likely pathogenic, and c.3303‐16A>G p.(?), classified as a variant of uncertain significance. PCR‐based RNA testing suggested that c.3303‐16A>G induces an in‐frame deletion: r.3303_3458del p.Asp1102_Arg1153del, although the normal transcript is still produced. These results are consistent with both SRPS type I/III in the first case and orofaciodigital syndrome in the second, an unprecedented description. This work thus improves the clinical and molecular knowledge of the phenotypes associated with splicing variants in the DYNC2H1 gene. [ABSTRACT FROM AUTHOR]

The fidelity of splice site selection is critical for proper gene expression. In particular, proper recognition of 3-splice site (3SS) sequences by the spliceosome is challenging considering the low complexity of the 3SS consensus sequence YAG. Here, we show that absence of the Prp18p splicing factor results in genome-wide activation of alternative 3SS in S. cerevisiae, including highly unusual non-YAG sequences. Usage of these non-canonical 3SS in the absence of Prp18p is enhanced by upstream poly(U) tracts and by their potential to interact with the first intronic nucleoside, allowing them to dock in the spliceosome active site instead of the normal 3SS. The role of Prp18p in 3SS fidelity is facilitated by interactions with Slu7p and Prp8p, but cannot be fulfilled by Slu7p, identifying a unique role for Prp18p in 3SS fidelity. This fidelity function is synergized by the downstream proofreading activity of the Prp22p helicase, but is independent from another late splicing helicase, Prp43p. Our results show that spliceosomes exhibit remarkably relaxed 3SS sequence usage in the absence of Prp18p and identify a network of spliceosomal interactions centered on Prp18p which are required to promote the fidelity of the recognition of consensus 3SS sequences.

Loss of nuclear TDP-43 is a hallmark of neurodegeneration in TDP-43 proteinopathies, including amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). TDP-43 mislocalization results in cryptic splicing and polyadenylation of pre-messenger RNAs (pre-mRNAs) encoding stathmin-2 (also known as SCG10), a protein that is required for axonal regeneration. We found that TDP-43 binding to a GU-rich region sterically blocked recognition of the cryptic 3 splice site in STMN2 pre-mRNA. Targeting dCasRx or antisense oligonucleotides (ASOs) suppressed cryptic splicing, which restored axonal regeneration and stathmin-2-dependent lysosome trafficking in TDP-43-deficient human motor neurons. In mice that were gene-edited to contain human STMN2 cryptic splice-polyadenylation sequences, ASO injection into cerebral spinal fluid successfully corrected Stmn2 pre-mRNA misprocessing and restored stathmin-2 expression levels independently of TDP-43 binding.

Alternative splicing enables higher eukaryotes to expand mRNA diversity from a finite number of genes through highly combinatorial splice site selection mechanisms that are influenced by the sequence of competing splice sites, cis-regulatory elements binding trans-acting factors, the length of exons and introns harbouring alternative splice sites and RNA secondary structures at putative splice junctions. To test the hypothesis that the intron definition or exon definition modes of splice site recognition direct the selection of alternative splice patterns, we created a database of alternative splice site usage (ALTssDB). When alternative splice sites are embedded within short introns (intron definition), the 5' and 3' splice sites closest to each other across the intron preferentially pair, consistent with previous observations. However, when alternative splice sites are embedded within large flanking introns (exon definition), the 5' and 3' splice sites closest to each other across the exon are preferentially selected. Thus, alternative splicing decisions are influenced by the intron and exon definition modes of splice site recognition. The results demonstrate that the spliceosome pairs splice sites that are closest in proximity within the unit of initial splice site selection.

The spliceosome undergoes extensive rearrangements as it assembles onto precursor messenger RNAs. In the earliest assembly step, U1snRNA identifies the 5 splice site. However, U1snRNA leaves the spliceosome relatively early in assembly, and 5 splice site identity is subsequently maintained through interactions with U6snRNA, protein factor PRP8, and other components during the rearrangements that build the catalytic site. Using a forward genetic screen in Caenorhabditis elegans, we have identified suppressors of a locomotion defect caused by a 5ss mutation. Here we report three new suppressor alleles from this screen, two in PRP8 and one in SNRNP200/BRR2. mRNASeq studies of these suppressor strains indicate that they also affect specific native alternative 5ss, especially for suppressor PRP8 D1549N. A strong suppressor at the unstructured N-terminus of SNRNP200, N18K, indicates a novel role for this region. By examining distinct changes in the splicing of native genes, examining double mutants between suppressors, comparing these new suppressors to previously identified splicing suppressors from yeast, and mapping conserved suppressor residues onto cryoEM structural models of assembling human spliceosomes, we conclude that there are multiple interactions at multiple stages in spliceosome assembly responsible for maintaining the initial 5ss identified by U1snRNA for entry into the catalytic core.

Objectives: Methylmalonic acidemia (MMA) is a rare inborn genetic disorder that is characterized by increased levels of methylmalonic acid in blood plasma and urine. Isolated methylmalonic acidemia is one of the most common types of MMA and is caused by mutations in the gene encoding methyl-malonyl coenzyme A mutase (MMUT). In this study, we investigated the possible mechanisms underlying the symptoms of isolated MMA in a patient by molecular analysis. Methods: PCR amplification and Sanger sequencing analysis was performed to identify variants in the MMUT gene in the proband and his family. Furthermore, minigene constructs were generated to validate the splicing defects in the MMUT gene variant identified in the proband. Results: The 3-year-old patient was admitted to the hospital with symptoms of MMA, including fever, convulsions, and vomiting. He showed metabolic acidosis, high levels of methylmalonic acid in blood and urine, and normal blood homocysteine levels. Genetic analysis demonstrated that the patient was a compound heterozygous carrier of two variants in the MMUT gene: a missense c.278G > A variant that has already been reported in a patient with the severe mut⁰ phenotype; and a novel splice site variant c.2125-2A > G. RT-PCR analysis showed that, while the novel variant clearly alters splicing, a minor amount of a full-length transcript is generated, suggesting that a wild-type protein may be produced although at a lower quantitative level. The patient's condition improved after treatment with vitamin B12. Serious complications were not reported during follow-up at age 5. Conclusions: We identified a novel splice site variant that partially disrupts normal splicing of the MMUT pre-mRNA. Production of a reduced amount of full-length transcript is responsible for the mild clinical phenotype observed in this patient. Functional studies have proven useful in exploring the genotype‐phenotype association and in providing guidance for the genetic diagnosis of MMA.

The major spliceosome mediates pre-mRNA splicing by recognizing the highly conserved sequences at the 5' and 3' splice sites and the branch point. More than 150 proteins participate in the splicing process and are organized in the spliceosomal A, B, and C complexes. FRA10AC1 is a peripheral protein of the spliceosomal C complex and its ortholog in the green alga facilitates recognition or interaction with splice sites. We identified biallelic pathogenic variants in FRA10AC1 in five individuals from three consanguineous families. The two unrelated Patients 1 and 2 with loss-of-function variants showed developmental delay, intellectual disability, and no speech, while three siblings with the c.494_496delAAG (p.Glu165del) variant had borderline to mild intellectual disability. All patients had microcephaly, hypoplasia or agenesis of the corpus callosum, growth retardation, and craniofacial dysmorphism. FRA10AC1 transcripts and proteins were drastically reduced or absent in fibroblasts of Patients 1 and 2. In a heterologous expression system, the p.Glu165del variant impacts intrinsic stability of FRA10AC1 but does not affect its nuclear localization. By co-immunoprecipitation, we found ectopically expressed HA-FRA10AC1 in complex with endogenous DGCR14, another component of the spliceosomal C complex, while the splice factors CHERP, NKAP, RED, and SF3B2 could not be co-immunoprecipitated. Using an in vitro splicing reporter assay, we did not obtain evidence for FRA10AC1 deficiency to suppress missplicing events caused by mutations in the highly conserved dinucleotides of 5' and 3' splice sites in an in vitro splicing assay in patient-derived fibroblasts. Our data highlight the importance of specific peripheral spliceosomal C complex proteins for neurodevelopment. It remains possible that FRA10AC1 may have other and/or additional cellular functions, such as coupling of transcription and splicing reactions.

DMD pathogenic variants for Duchenne and Becker muscular dystrophy are detectable with high sensitivity by standard clinical exome analyses of genomic DNA. However, up to 7% of DMD mutations are deep intronic and analysis of muscle-derived RNA is an important diagnostic step for patients who have negative genomic testing but abnormal dystrophin expression in muscle. In this study, muscle biopsies were evaluated from 19 patients with clinical features of a dystrophinopathy, but negative clinical DMD mutation analysis. Reverse transcription-polymerase chain reaction or high-throughput RNA sequencing methods identified 19 mutations with one of three pathogenic pseudoexon types: deep intronic point mutations, deletions or insertions, and translocations. In association with point mutations creating intronic splice acceptor sites, we observed the first examples of DMD pseudo 3'-terminal exon mutations causing high efficiency transcription termination within introns. This connection between splicing and premature transcription termination is reminiscent of U1 snRNP-mediating telescripting in sustaining RNA polymerase II elongation across large genes, such as DMD. We propose a novel classification of three distinct types of mutations identifiable by muscle RNA analysis, each of which differ in potential treatment approaches. Recognition and appropriate characterization may lead to therapies directed toward full-length dystrophin expression for some patients.

MotivationCircular RNA (circRNA) is a novel class of long non-coding RNAs that have been broadly discovered in the eukaryotic transcriptome. The circular structure arises from a non-canonical splicing process, where the donor site backspliced to an upstream acceptor site. These circRNA sequences are conserved across species. More importantly, rising evidence suggests their vital roles in gene regulation and association with diseases. As the fundamental effort toward elucidating their functions and mechanisms, several computational methods have been proposed to predict the circular structure from the primary sequence. Recently, advanced computational methods leverage deep learning to capture the relevant patterns from RNA sequences and model their interactions to facilitate the prediction. However, these methods fail to fully explore positional information of splice junctions and their deep interaction.ResultsWe present a robust end-to-end framework, Junction Encoder with Deep Interaction (JEDI), for circRNA prediction using only nucleotide sequences. JEDI first leverages the attention mechanism to encode each junction site based on deep bidirectional recurrent neural networks and then presents the novel cross-attention layer to model deep interaction among these sites for backsplicing. Finally, JEDI can not only predict circRNAs but also interpret relationships among splice sites to discover backsplicing hotspots within a gene region. Experiments demonstrate JEDI significantly outperforms state-of-the-art approaches in circRNA prediction on both isoform level and gene level. Moreover, JEDI also shows promising results on zero-shot backsplicing discovery, where none of the existing approaches can achieve.Availability and implementationThe implementation of our framework is available at https://github.com/hallogameboy/JEDI.Supplementary informationSupplementary data are available at Bioinformatics online.

Hereditary Diffuse Gastric Cancer (HDGC) syndrome is associated with CDH1 germline likely pathogenic/pathogenic variants. Carriers of CDH1 germline likely pathogenic/pathogenic variants are predisposed to diffuse gastric cancer and lobular breast cancer. This study aims to classify the CDH1 c.[715G>A] missense variant identified in a diffuse gastric cancer prone family by performing splicing studies. RT-PCR and subsequent cloning experiments were performed to investigate whether this variant completely disrupts normal splicing. This variant preferentially abolishes normal splicing through activation of a cryptic 3' acceptor splice site within exon 6 of CDH1, presumably leading to a premature protein truncation within first extracellular domain repeat of E-cadherin protein. Our results contributed to evidence necessary to resolve pathogenicity classification of this variant, indicating that this variant is to be classified as pathogenic.

Trafficking protein particle (TRAPP) complexes, which include the TRAPPC4 protein, regulate membrane trafficking between lipid organelles in a process termed vesicular tethering. TRAPPC4 was recently implicated in a recessive neurodevelopmental condition in four unrelated families due to a shared c.454+3A>G splice variant. Here, we report 23 patients from 17 independent families with an early-infantile-onset neurodegenerative presentation, where we also identified the homozygous variant hg38:11:119020256 A>G (NM_016146.5:c.454+3A>G) in TRAPPC4 through exome or genome sequencing. No other clinically relevant TRAPPC4 variants were identified among any of over 10,000 patients with neurodevelopmental conditions. We found the carrier frequency of TRAPPC4 c.454+3A>G was 2.4-5.4 per 10,000 healthy individuals. Affected individuals with the homozygous TRAPPC4 c.454+3A>G variant showed profound psychomotor delay, developmental regression, early-onset epilepsy, microcephaly and progressive spastic tetraplegia. Based upon RNA sequencing, the variant resulted in partial exon 3 skipping and generation of an aberrant transcript owing to use of a downstream cryptic splice donor site, predicting a premature stop codon and nonsense mediated decay. These data confirm the pathogenicity of the TRAPPC4 c.454+3A>G variant, and refine the clinical presentation of TRAPPC4-related encephalopathy.

Metabolic labeling is a widely used tool to investigate different aspects of pre-mRNA splicing and RNA turnover. The labeling technology takes advantage of native cellular machineries where a nucleotide analog is readily taken up and incorporated into nascent RNA. One such analog is 4-thiouridine (4sU). Previous studies demonstrated that the uptake of 4sU at elevated concentrations (>50μM) and extended exposure led to inhibition of rRNA synthesis and processing, presumably induced by changes in RNA secondary structure. Thus, it is possible that 4sU incorporation may also interfere with splicing efficiency. To test this hypothesis, we carried out splicing analyses of pre-mRNA substrates with varying levels of 4sU incorporation (0-100%). We demonstrate that increased incorporation of 4sU into pre-mRNAs decreased splicing efficiency. The overall impact of 4sU labeling on pre-mRNA splicing efficiency negatively correlates with the strength of splice site signals such as the 3' and the 5' splice sites. Introns with weaker splice sites are more affected by the presence of 4sU. We also show that transcription by T7 polymerase and pre-mRNA degradation kinetics were impacted at the highest levels of 4sU incorporation. Increased incorporation of 4sU caused elevated levels of abortive transcripts, and fully labeled pre-mRNA is more stable than its uridine-only counterpart. Cell culture experiments show that a small number of alternative splicing events were modestly, but statistically significantly influenced by metabolic labeling with 4sU at concentrations considered to be tolerable (40 μM). We conclude that at high 4sU incorporation rates small, but noticeable changes in pre-mRNA splicing can be detected when splice sites deviate from consensus. Given these potential 4sU artifacts, we suggest that appropriate controls for metabolic labeling experiments need to be included in future labeling experiments.

Alternative splicing is responsible for much of the transcriptomic and proteomic diversity observed in eukaryotes and involves combinatorial regulation by many cis-acting elements and trans-acting factors. SR and hnRNP splicing regulatory proteins often have opposing effects on splicing efficiency depending on where they bind the pre-mRNA relative to the splice site. Position-dependent splicing repression occurs at spliceosomal E-complex, suggesting that U1 snRNP binds but cannot facilitate higher order spliceosomal assembly. To test the hypothesis that the structure of U1 snRNA changes during activation or repression, we developed a method to structure-probe native U1 snRNP in enriched conformations that mimic activated or repressed spliceosomal E-complexes. While the core of U1 snRNA is highly structured, the 5' end of U1 snRNA shows different SHAPE reactivities and psoralen crosslinking efficiencies depending on where splicing regulatory elements are located relative to the 5' splice site. A motif within the 5' splice site binding region of U1 snRNA is more reactive toward SHAPE electrophiles when repressors are bound, suggesting U1 snRNA is bound, but less base-paired. These observations demonstrate that splicing regulators modulate splice site selection allosterically.

Although immune checkpoint blockade (ICB) therapy has revolutionized cancer treatment, many patients do not respond or develop resistance to ICB. N6 -methylation of adenosine (m6A) in RNA regulates many pathophysiological processes. Here, we show that deletion of the m6A demethylase Alkbh5 sensitized tumors to cancer immunotherapy. Alkbh5 has effects on m6A density and splicing events in tumors during ICB. Alkbh5 modulates Mct4/Slc16a3 expression and lactate content of the tumor microenvironment and the composition of tumor-infiltrating Treg and myeloid-derived suppressor cells. Importantly, a small-molecule Alkbh5 inhibitor enhanced the efficacy of cancer immunotherapy. Notably, the ALKBH5 gene mutation and expression status of melanoma patients correlate with their response to immunotherapy. Our results suggest that m6A demethylases in tumor cells contribute to the efficacy of immunotherapy and identify ALKBH5 as a potential therapeutic target to enhance immunotherapy outcome in melanoma, colorectal, and potentially other cancers.

The decoy exon model has been proposed to regulate a subset of intron retention (IR) events involving predominantly larger introns (>1 kb). Splicing reporter studies have shown that decoy splice sites are essential for activity, suggesting that decoys act by engaging intron-terminal splice sites and competing with cross-intron interactions required for intron excision. The decoy model predicts that antisense oligonucleotides may be able to block decoy splice sites in endogenous pre-mRNA, thereby reducing IR and increasing productive gene expression. Indeed, we now demonstrate that targeting a decoy 5' splice site in the O-GlcNAc transferase (OGT) gene reduced IR from ∼80% to ∼20% in primary human erythroblasts, accompanied by increases in spliced OGT RNA and OGT protein expression. The remaining OGT IR was refractory to antisense treatment and might be mediated by independent mechanism(s). In contrast, other retained introns were strongly dependent on decoy function, since antisense targeting of decoy 5' splice sites greatly reduced (SNRNP70) or nearly eliminated (SF3B1) IR in two widely expressed splicing factors, and also greatly reduced IR in transcripts encoding the erythroid-specific structural protein, α-spectrin (SPTA1). These results show that modulating decoy exon function can dramatically alter IR and suggest that dynamic regulation of decoy exons could be a mechanism to fine-tune gene expression post-transcriptionally in many cell types.

While splicing changes caused by somatic mutations in SF3B1 are known, identifying full-length isoform changes may better elucidate the functional consequences of these mutations. We report nanopore sequencing of full-length cDNA from CLL samples with and without SF3B1 mutation, as well as normal B cell samples, giving a total of 149 million pass reads. We present FLAIR (Full-Length Alternative Isoform analysis of RNA), a computational workflow to identify high-confidence transcripts, perform differential splicing event analysis, and differential isoform analysis. Using nanopore reads, we demonstrate differential 3' splice site changes associated with SF3B1 mutation, agreeing with previous studies. We also observe a strong downregulation of intron retention events associated with SF3B1 mutation. Full-length transcript analysis links multiple alternative splicing events together and allows for better estimates of the abundance of productive versus unproductive isoforms. Our work demonstrates the potential utility of nanopore sequencing for cancer and splicing research.

RNA splicing and spliceosome assembly in eukaryotes occur mainly during transcription. However, co-transcriptional splicing has not yet been explored in plants. Here, we built transcriptomes of nascent chromatin RNAs in Arabidopsis thaliana and showed that nearly all introns undergo co-transcriptional splicing, which occurs with higher efficiency for introns in protein-coding genes than for those in noncoding RNAs. Total intron number and intron position are two predominant features that correlate with co-transcriptional splicing efficiency, and introns with alternative 5' or 3' splice sites are less efficiently spliced. Furthermore, we found that mutations in genes encoding trans-acting proteins lead to more introns with increased splicing defects in nascent RNAs than in mature RNAs, and that introns with increased splicing defects in mature RNAs are inefficiently spliced at the co-transcriptional level. Collectively, our results not only uncovered widespread co-transcriptional splicing in Arabidopsis but also identified features that may affect or be affected by co-transcriptional splicing efficiency.

Pathogenic variants in the ATP-binding cassette transporter A4 (ABCA4) gene cause a continuum of retinal disease phenotypes, including Stargardt disease. Noncanonical splice site (NCSS) and deep-intronic variants constitute a large fraction of disease-causing alleles, defining the functional consequences of which remains a challenge. We aimed to determine the effect on splicing of nine previously reported or unpublished NCSS variants, one near exon splice variant and nine deep-intronic variants in ABCA4, using in vitro splice assays in human embryonic kidney 293T cells. Reverse transcription-polymerase chain reaction and Sanger sequence analysis revealed splicing defects for 12 out of 19 variants. Four deep-intronic variants create pseudoexons or elongate the upstream exon. Furthermore, eight NCSS variants cause a partial deletion or skipping of one or more exons in messenger RNAs. Among the 12 variants, nine lead to premature stop codons and predicted truncated ABCA4 proteins. At least two deep-intronic variants affect splice enhancer and silencer motifs and, therefore, these conserved sequences should be carefully evaluated when predicting the outcome of NCSS and deep-intronic variants.

The generation of protein coding mRNAs from pre-mRNA is a fundamental biological process that is required for gene expression. Alternative pre-mRNA splicing is responsible for much of the transcriptomic and proteomic diversity observed in higher order eukaryotes. Aberrations that disrupt regular alternative splicing patterns are known to cause human diseases, including various cancers. Alternative splicing is a combinatorial process, meaning many factors affect which two splice sites are ligated together. The features that dictate exon inclusion are comprised of splice site strength, intron-exon architecture, RNA secondary structure, splicing regulatory elements, promoter use and transcription speed by RNA polymerase and the presence of post-transcriptional nucleotide modifications. A comprehensive view of all of the factors that influence alternative splicing decisions is necessary to predict splicing outcomes and to understand the molecular basis of disease. This article is part of a Special Issue entitled: RNA structure and splicing regulation edited by Francisco Baralle, Ravindra Singh and Stefan Stamm.

Myelodysplastic syndromes (MDS) with ring sideroblasts are hematopoietic stem cell disorders with erythroid dysplasia and mutations in the SF3B1 splicing factor gene. Patients with MDS with SF3B1 mutations often accumulate excessive tissue iron, even in the absence of transfusions, but the mechanisms that are responsible for their parenchymal iron overload are unknown. Body iron content, tissue distribution, and the supply of iron for erythropoiesis are controlled by the hormone hepcidin, which is regulated by erythroblasts through secretion of the erythroid hormone erythroferrone (ERFE). Here, we identified an alternative ERFE transcript in patients with MDS with the SF3B1 mutation. Induction of this ERFE transcript in primary SF3B1-mutated bone marrow erythroblasts generated a variant protein that maintained the capacity to suppress hepcidin transcription. Plasma concentrations of ERFE were higher in patients with MDS with an SF3B1 gene mutation than in patients with SF3B1 wild-type MDS. Thus, hepcidin suppression by a variant ERFE is likely responsible for the increased iron loading in patients with SF3B1-mutated MDS, suggesting that ERFE could be targeted to prevent iron-mediated toxicity. The expression of the variant ERFE transcript that was restricted to SF3B1-mutated erythroblasts decreased in lenalidomide-responsive anemic patients, identifying variant ERFE as a specific biomarker of clonal erythropoiesis.

Deletions in the 16.6 kb mitochondrial genome have been implicated in numerous disorders that often display muscular and/or neurological symptoms due to the high-energy demands of these tissues. We describe a catalogue of 4489 putative mitochondrial DNA (mtDNA) deletions, including their frequency and relative read rate, using a combinatorial approach of mitochondria-targeted PCR, next-generation sequencing, bioinformatics, post-hoc filtering, annotation, and validation steps. Our bioinformatics pipeline uses MapSplice, an RNA-seq splice junction detection algorithm, to detect and quantify mtDNA deletion breakpoints rather than mRNA splices. Analyses of 93 samples from postmortem brain and blood found (i) the 4977 bp 'common deletion' was neither the most frequent deletion nor the most abundant; (ii) brain contained significantly more deletions than blood; (iii) many high frequency deletions were previously reported in MitoBreak, suggesting they are present at low levels in metabolically active tissues and are not exclusive to individuals with diagnosed mitochondrial pathologies; (iv) many individual deletions (and cumulative metrics) had significant and positive correlations with age and (v) the highest deletion burdens were observed in major depressive disorder brain, at levels greater than Kearns-Sayre Syndrome muscle. Collectively, these data suggest the Splice-Break pipeline can detect and quantify mtDNA deletions at a high level of resolution.

Many viruses shut off host gene expression to inhibit antiviral responses. Viral proteins and host proteins required for viral replication are typically spared in this process, but the mechanisms of target selectivity during host shutoff remain poorly understood. Using transcriptome-wide and targeted reporter experiments, we demonstrate that the influenza A virus endoribonuclease PA-X usurps RNA splicing to selectively target host RNAs for destruction. Proximity-labeling proteomics reveals that PA-X interacts with cellular RNA processing proteins, some of which are partially required for host shutoff. Thus, PA-X taps into host nuclear pre-mRNA processing mechanisms to destroy nascent mRNAs shortly after their synthesis. This mechanism sets PA-X apart from other viral host shutoff proteins that target actively translating mRNAs in the cytoplasm. Our study reveals a unique mechanism of host shutoff that helps us understand how influenza viruses suppress host gene expression.

Pre-mRNA splicing must occur with extremely high fidelity. Spliceosomes assemble onto pre-mRNA guided by specific sequences (5' splice site, 3' splice site, and branchpoint). When splice sites are mutated, as in many hereditary diseases, the spliceosome can aberrantly select nearby pseudo- or "cryptic" splice sites, often resulting in nonfunctional protein. How the spliceosome distinguishes authentic splice sites from cryptic splice sites is poorly understood. We performed a Caenorhabditis elegans genetic screen to find cellular factors that affect the frequency with which the spliceosome uses cryptic splice sites and identified two alleles in core spliceosome component Prp8 that alter cryptic splicing frequency. Subsequent complementary genetic and structural analyses in yeast implicate these alleles in the stability of the spliceosome's catalytic core. However, despite a clear effect on cryptic splicing, high-throughput mRNA sequencing of these prp-8 mutant C. elegans reveals that overall alternative splicing patterns are relatively unchanged. Our data suggest the spliceosome evolved intrinsic mechanisms to reduce the occurrence of cryptic splicing and that these mechanisms are distinct from those that impact alternative splicing.

The epithelial cell adhesion molecule gene (EPCAM, previously known as TACSTD1 or TROP1) encodes a membrane-bound protein that is localized to the basolateral membrane of epithelial cells and is overexpressed in some tumors. Biallelic mutations in EPCAM cause congenital tufting enteropathy (CTE), which is a rare chronic diarrheal disorder presenting in infancy. Monoallelic deletions of the 3' end of EPCAM that silence the downstream gene, MSH2, cause a form of Lynch syndrome, which is a cancer predisposition syndrome associated with loss of DNA mismatch repair. Here, we report 13 novel EPCAM mutations from 17 CTE patients from two separate centers, review EPCAM mutations associated with CTE and Lynch syndrome, and structurally model pathogenic missense mutations. Statistical analyses indicate that the c.499dupC (previously reported as c.498insC) frameshift mutation was associated with more severe treatment regimens and greater mortality in CTE, whereas the c.556-14A>G and c.491+1G>A splice site mutations were not correlated with treatments or outcomes significantly different than random simulation. These findings suggest that genotype-phenotype correlations may be useful in contributing to management decisions of CTE patients. Depending on the type and nature of EPCAM mutation, one of two unrelated diseases may occur, CTE or Lynch syndrome.

Alternative splicing diversifies mRNA transcripts in human cells. While the spliceosome pairs exons with a high degree of accuracy, the rates of rare aberrant and non-canonical pre-mRNA splicing have not been evaluated at the nucleotide level to determine the quantity and identity of these events across splice junctions. Using ultra-deep sequencing the frequency of aberrant and non-canonical splicing events for three splice junctions flanking exon 7 of SMN1 were determined at single nucleotide resolution. After correction for background noise introduced by PCR amplification and sequencing steps, pre-mRNA splicing was shown to maintain a low overall rate of aberrant and non-canonically spliced events. Several previously unannotated splicing events across 3 exon|intron junctions in SMN1 were identified. Mutations within SMN exon 7 were shown to affect splicing fidelity by modulating RNA secondary structures, by altering the binding site of regulatory proteins and by changing the 5' splice site strength. Mutations also create a truncated SMN1 exon 7 through the introduction of a de novo non-canonical 5' splice site. The results from the ultra-deep sequencing approach highlight the impressive fidelity of pre-mRNA splicing and demonstrate that the immediate sequence context around splice sites is the main driving force behind non-canonical splice site pairing.

In a recent issue of Molecular Cell, Boehm et al. (2018), Blazquez et al. (2018), and Gonatopoulos-Pournatzis et al. (2018) uncover novel mechanisms by which the cell regulates splicing of cryptic splice sites and microexons.

The spliceosome catalyzes the removal of introns from pre-messenger RNA (mRNA) and subsequent pairing of exons with remarkable fidelity. Some exons are known to be skipped or included in the mature mRNA in a cell type- or context-dependent manner (cassette exons), thereby contributing to the diversification of the human proteome. Interestingly, splicing is initiated (and sometimes completed) co-transcriptionally. Here, we develop a kinetic mathematical modeling framework to investigate alternative co-transcriptional splicing (CTS) and, specifically, the control of cassette exons' inclusion. We show that when splicing is co-transcriptional, default splice patterns of exon inclusion are more likely than when splicing is post-transcriptional, and that certain exons are more likely to be regulatable (i.e. cassette exons) than others, based on the exon-intron structure context. For such regulatable exons, transcriptional elongation rates may affect splicing outcomes. Within the CTS paradigm, we examine previously described hypotheses of co-operativity between splice sites of short introns (i.e. 'intron definition') or across short exons (i.e. 'exon definition'), and find that models encoding these faithfully recapitulate observations in the fly and human genomes, respectively.

The group II intron and the spliceosome share a common active site architecture and are thought to be evolutionarily related. Here we report the 3.7 Å crystal structure of a eukaryotic group II intron in the lariat-3' exon form, immediately preceding the second step of splicing, analogous to the spliceosomal P complex. This structure reveals the location of the intact 3' splice site within the catalytic core of the group II intron. The 3'-OH of the 5' exon is positioned in close proximity to the 3' splice site for nucleophilic attack and exon ligation. The active site undergoes conformational rearrangements with the catalytic triplex having different configurations before and after the second step of splicing. We describe a complete model for the second step of group II intron splicing that incorporates a dynamic catalytic triplex being responsible for creating the binding pocket for 3' splice site capture.

The tri-snRNP 27K protein is a component of the human U4/U6-U5 tri-snRNP and contains an N-terminal phosphorylated RS domain. In a forward genetic screen in C. elegans, we previously identified a dominant mutation, M141T, in the highly-conserved C-terminal region of this protein. The mutant allele promotes changes in cryptic 5 splice site choice. To better understand the function of this poorly characterized splicing factor, we performed high-throughput mRNA sequencing analysis on worms containing this dominant mutation. Comparison of alternative splice site usage between the mutant and wild-type strains led to the identification of 26 native genes whose splicing changes in the presence of the snrp-27 mutation. The changes in splicing are specific to alternative 5 splice sites. Analysis of new alleles suggests that snrp-27 is an essential gene for worm viability. We performed a novel directed-mutation experiment in which we used the CRISPR-cas9 system to randomly generate mutations specifically at M141 of SNRP-27. We identified eight amino acid substitutions at this position that are viable, and three that are homozygous lethal. All viable substitutions at M141 led to varying degrees of changes in alternative 5 splicing of native targets. We hypothesize a role for this SR-related factor in maintaining the position of the 5 splice site as U1snRNA trades interactions at the 5 end of the intron with U6snRNA and PRP8 as the catalytic site is assembled.

During terminal erythropoiesis, the splicing machinery in differentiating erythroblasts executes a robust intron retention (IR) program that impacts expression of hundreds of genes. We studied IR mechanisms in the SF3B1 splicing factor gene, which expresses ∼50% of its transcripts in late erythroblasts as a nuclear isoform that retains intron 4. RNA-seq analysis of nonsense-mediated decay (NMD)-inhibited cells revealed previously undescribed splice junctions, rare or not detected in normal cells, that connect constitutive exons 4 and 5 to highly conserved cryptic cassette exons within the intron. Minigene splicing reporter assays showed that these cassettes promote IR. Genome-wide analysis of splice junction reads demonstrated that cryptic noncoding cassettes are much more common in large (>1 kb) retained introns than they are in small retained introns or in nonretained introns. Functional assays showed that heterologous cassettes can promote retention of intron 4 in the SF3B1 splicing reporter. Although many of these cryptic exons were spliced inefficiently, they exhibited substantial binding of U2AF1 and U2AF2 adjacent to their splice acceptor sites. We propose that these exons function as decoys that engage the intron-terminal splice sites, thereby blocking cross-intron interactions required for excision. Developmental regulation of decoy function underlies a major component of the erythroblast IR program.

We previously investigated the transcriptome and proteome profiles of the murine ocular lens at six developmental time points including two embryonic (E15 and E18) and four postnatal time points (P0, P3, P6, and P9). Here, we extend our analyses to identify novel transcripts and peptides in developing  mouse lens. We identified a total of 9,707 novel transcripts and 325 novel fusion genes in developing mouse lens. Additionally, we identified 13,281 novel alternative splicing (AS) events in mouse lens including 6,990 exon skipping (ES), 2,447 alternative 3' splice site (A3SS), 1,900 alternative 5' splice site (A5SS), 1,771 mutually exclusive exons (MXE), and 173 intron retention (IR). Finally, we integrated our OMIC (Transcriptome and Proteome) datasets identifying 20 novel peptides in mouse lens. All 20 peptides were validated through matching MS/MS spectra of synthetic peptides. To the best of our knowledge, this is the first report integrating OMIC datasets to identify novel peptides in developing murine lens.

Prp8 is an essential protein that regulates spliceosome assembly and conformation during pre-mRNA splicing. Recent cryo-EM structures of the spliceosome model Prp8 as a scaffold for the spliceosome's catalytic U snRNA components. Using a new amino acid probing strategy, we identified a dynamic region in human Prp8 that is positioned to stabilize the pre-mRNA in the spliceosome active site through interactions with U5 snRNA. Mutagenesis of the identified Prp8 residues in yeast indicates a role in 5' splice site recognition. Genetic interactions with spliceosome proteins Isy1, which buttresses the intron branch point, and Snu114, a regulatory GTPase that directly contacts Prp8, further corroborate a role for the same Prp8 residues in substrate positioning and activation. Together the data suggest that adjustments in interactions between Prp8 and U5 snRNA help establish proper positioning of the pre-mRNA into the active site to enhance 5' splice site fidelity.

Alternative pre-mRNA splicing plays a major role in expanding the transcript output of human genes. This process is regulated, in part, by the interplay of trans-acting RNA binding proteins (RBPs) with myriad cis-regulatory elements scattered throughout pre-mRNAs. These molecular recognition events are critical for defining the protein-coding sequences (exons) within pre-mRNAs and directing spliceosome assembly on noncoding regions (introns). One of the earliest events in this process is recognition of the 3' splice site (3'ss) by U2 small nuclear RNA auxiliary factor 2 (U2AF2). Splicing regulators, such as the heterogeneous nuclear ribonucleoprotein A1 (HNRNPA1), influence spliceosome assembly both in vitro and in vivo, but their mechanisms of action remain poorly described on a global scale. HNRNPA1 also promotes proofreading of 3'ss sequences though a direct interaction with the U2AF heterodimer. To determine how HNRNPA1 regulates U2AF-RNA interactions in vivo, we analyzed U2AF2 RNA binding specificity using individual-nucleotide resolution crosslinking immunoprecipitation (iCLIP) in control and HNRNPA1 overexpression cells. We observed changes in the distribution of U2AF2 crosslinking sites relative to the 3'ss of alternative cassette exons but not constitutive exons upon HNRNPA1 overexpression. A subset of these events shows a concomitant increase of U2AF2 crosslinking at distal intronic regions, suggesting a shift of U2AF2 to "decoy" binding sites. Of the many noncanonical U2AF2 binding sites, Alu-derived RNA sequences represented one of the most abundant classes of HNRNPA1-dependent decoys. We propose that one way HNRNPA1 regulates exon definition is to modulate the interaction of U2AF2 with decoy or bona fide 3'ss.

For the past decade, cancer genomic studies have focused on mutations leading to splice-site disruption, overlooking those having splice-creating potential. Here, we applied a bioinformatic tool, MiSplice, for the large-scale discovery of splice-site-creating mutations (SCMs) across 8,656 TCGA tumors. We report 1,964 originally mis-annotated mutations having clear evidence of creating alternative splice junctions. TP53 and GATA3 have 26 and 18 SCMs, respectively, and ATRX has 5 from lower-grade gliomas. Mutations in 11 genes, including PARP1, BRCA1, and BAP1, were experimentally validated for splice-site-creating function. Notably, we found that neoantigens induced by SCMs are likely several folds more immunogenic compared to missense mutations, exemplified by the recurrent GATA3 SCM. Further, high expression of PD-1 and PD-L1 was observed in tumors with SCMs, suggesting candidates for immune blockade therapy. Our work highlights the importance of integrating DNA and RNA data for understanding the functional and the clinical implications of mutations in human diseases.

The spliceosome undergoes dramatic changes in a splicing cycle. Structures of B, Bact, C, C*, and intron lariat spliceosome complexes revealed mechanisms of 5'-splice site (ss) recognition, branching, and intron release, but lacked information on 3'-ss recognition, exon ligation, and exon release. Here we report a cryo-electron microscopy structure of the postcatalytic P complex at 3.3-angstrom resolution, revealing that the 3' ss is mainly recognized through non-Watson-Crick base pairing with the 5' ss and branch point. Furthermore, one or more unidentified proteins become stably associated with the P complex, securing the 3' exon and potentially regulating activity of the helicase Prp22. Prp22 binds nucleotides 15 to 21 in the 3' exon, enabling it to pull the intron-exon or ligated exons in a 3' to 5' direction to achieve 3'-ss proofreading or exon release, respectively.

The spliceosome is the complex molecular machinery that sequentially assembles on eukaryotic messenger RNA precursors to remove introns (pre-mRNA splicing), a physiologically regulated process altered in numerous pathologies. We report transcriptome-wide analyses upon systematic knock down of 305 spliceosome components and regulators in human cancer cells and the reconstruction of functional splicing factor networks that govern different classes of alternative splicing decisions. The results disentangle intricate circuits of splicing factor cross-regulation, reveal that the precise architecture of late-assembling U4/U6.U5 tri-small nuclear ribonucleoprotein (snRNP) complexes regulates splice site pairing, and discover an unprecedented division of labor among protein components of U1 snRNP for regulating exon definition and alternative 5' splice site selection. Thus, we provide a resource to explore physiological and pathological mechanisms of splicing regulation.

The two transesterification reactions of pre-mRNA splicing require highly complex yet well-controlled rearrangements of small nuclear RNAs and proteins (snRNP) in the spliceosome. The efficiency and accuracy of these reactions are critical for gene expression, as almost all human genes pass through pre-mRNA splicing. Key parameters that determine the splicing outcome are the length of the intron, the strengths of its splicing signals and gaps between them, and the presence of splicing controlling elements. In particular, the gap between the branchpoint (BP) and the 3' splice site (ss) of introns is a major determinant of the splicing efficiency. This distance falls within a small range across the introns of an organism. The constraints exist possibly because BP and 3'ss are recognized by BP-binding proteins, U2 snRNP and U2 accessory factors (U2AF) in a coordinated manner. Furthermore, varying distances between the two signals may also affect the second transesterification reaction since the intervening RNA needs to be accurately positioned within the complex RNP machinery. Splicing such pre-mRNAs requires cis-acting elements in the RNA and many trans-acting splicing regulators. Regulated pre-mRNA splicing with BP-distant 3'ss adds another layer of control to gene expression and promotes alternative splicing., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024. Published by Elsevier B.V.)

The mRNAs encoding the human papillomavirus type 16 (HPV16) E6 and E7 oncogene mRNAs are subjected to extensive alternative RNA splicing at multiple regulated splice sites. One of the most extensively used 5'-splice sites in the HPV16 genome is named SD880 and is located immediately downstream of the E7 open reading frame. Here, we show that a cluster of three GGG-motifs adjacent to HPV16 SD880 interacts with heterogeneous nuclear ribonucleoprotein (hnRNP) H that cooperates with SD880 to stimulate splicing to the upstream HPV16 3'-splice site SA742. This splice site is located in the E7 coding region and is required for the production of the HPV16 226^742 mRNA that encodes the E6^E7 fusion protein. Enhancement of HPV16 E6^E7 mRNA production by hnRNP H occurred at the expense of the intron-retained E6 mRNAs and the spliced E7 mRNAs, demonstrating that hnRNP H controls the relative levels of E6, E7, and E6^E7 proteins. Unexpectedly, overexpression of hnRNP H also promoted retention of the downstream E1 encoding intron and enhanced E1 protein production. We concluded that hnRNP H plays an important role in the HPV16 gene expression program.IMPORTANCEHere, we show that hnRNP H binds to multiple GGG-motifs downstream of human papillomavirus type 16 (HPV16) splice site SD880 and acts in concert with SD880 to promote expression of the HPV16 E6^E7 mRNA. The E6^E7 protein has been shown previously to stabilize the HPV16 E6 and E7 oncoproteins and may as such contribute to the carcinogenic properties of HPV16. In its capacity of major regulator of HPV16 oncogene expression, hnRNP H may be exploited as a target for antiviral drugs to HPV16., Competing Interests: The authors declare no conflict of interest.

The spliceosome is a eukaryotic multimegadalton RNA-protein complex that removes introns from transcripts. The spliceosome ensures the selection of each exon-intron boundary through multiple recognition events. Initially, the 5' splice site (5' SS) and branch site (BS) are bound by the U1 small nuclear ribonucleoprotein (snRNP) and the U2 snRNP, respectively, while the 3' SS is mostly determined by proximity to the branch site. A large number of splicing factors recognize the splice sites and recruit the snRNPs before the stable binding of the snRNPs occurs by base-pairing the snRNA to the transcript. Fidelity of this process is crucial, as mutations in splicing factors and U2 snRNP components are associated with many diseases. In recent years, major advances have been made in understanding how splice sites are selected in Saccharomyces cerevisiae and humans. Here, I review and discuss the current understanding of the recognition of splice sites by the spliceosome with a focus on recognition and binding of the branch site by the U2 snRNP in humans., (© 2024 Tholen; Published by Cold Spring Harbor Laboratory Press for the RNA Society.)

Splice site recognition is essential for defining the transcriptome. Drugs like risdiplam and branaplam change how human U1 snRNP recognizes particular 5' splice sites (5'SS) and promote U1 snRNP binding and splicing at these locations. Despite the therapeutic potential of 5'SS modulators, the complexity of their interactions and snRNP substrates have precluded defining a mechanism for 5'SS modulation. We have determined a sequential binding mechanism for modulation of -1A bulged 5'SS by branaplam using a combination of ensemble kinetic measurements and colocalization single molecule spectroscopy (CoSMoS). Our mechanism establishes that U1-C protein binds reversibly to U1 snRNP, and branaplam binds to the U1 snRNP/U1-C complex only after it has engaged with a -1A bulged 5'SS. Obligate orders of binding and unbinding explain how reversible branaplam interactions cause formation of long-lived U1 snRNP/5'SS complexes. Branaplam targets a ribonucleoprotein, not only an RNA duplex, and its action depends on fundamental properties of 5'SS recognition., (© 2024. The Author(s).)

Loss of function of the RNA-binding protein TDP-43 (TDP-LOF) is a hallmark of amyotrophic lateral sclerosis (ALS) and other neurodegenerative disorders. Here we describe TDP-REG, which exploits the specificity of cryptic splicing induced by TDP-LOF to drive protein expression when and where the disease process occurs. The SpliceNouveau algorithm combines deep learning with rational design to generate customizable cryptic splicing events within protein-coding sequences. We demonstrate that expression of TDP-REG reporters is tightly coupled to TDP-LOF in vitro and in vivo. TDP-REG enables genomic prime editing to ablate the UNC13A cryptic donor splice site specifically upon TDP-LOF. Finally, we design TDP-REG vectors encoding a TDP-43/Raver1 fusion protein that rescues key pathological cryptic splicing events, paving the way for the development of precision therapies for TDP43-related disorders.

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., Competing Interests: Declaration of interests O.A.-W. is a founder and scientific advisor of Codify Therapeutics, holds equity in this company, and receives research funding from this company. O.A.-W. has served as a consultant for Foundation Medicine Inc., Merck, Prelude Therapeutics, Amphista Therapeutics, MagnetBio, and Janssen, and is on the Scientific Advisory Board of Envisagenics Inc., Harmonic Discovery Inc., and Pfizer Boulder; O.A.-W. has received prior research funding from H3B Biomedicine, Nurix Therapeutics, Minovia Therapeutics, and LOXO Oncology unrelated to the current manuscript., (Copyright © 2024 The Authors. Published by Elsevier Inc. All rights reserved.)

The process of splicing messenger RNA to remove introns plays a central role in creating genes and gene variants. We describe Splam, a novel method for predicting splice junctions in DNA using deep residual convolutional neural networks. Unlike previous models, Splam looks at a 400-base-pair window flanking each splice site, reflecting the biological splicing process that relies primarily on signals within this window. Splam also trains on donor and acceptor pairs together, mirroring how the splicing machinery recognizes both ends of each intron. Compared to SpliceAI, Splam is consistently more accurate, achieving 96% accuracy in predicting human splice junctions., (© 2024. The Author(s).)

Splicing is an important step of gene expression in all eukaryotes. Splice sites might be used with different efficiency, giving rise to alternative splicing products. At the same time, splice sites might be used at a variable rate. We used 5-ethynyl uridine labeling to sequence a nascent transcriptome of HeLa cells and deduced the rate of splicing for each donor and acceptor splice site. The following correlation analysis showed a correspondence of primary transcript features with the rate of splicing. Some dependencies we revealed were anticipated, such as a splicing rate decrease with a decreased complementarity of the donor splice site to U1 and acceptor sites to U2 snRNAs. Other dependencies were more surprising, like a negative influence of a distance to the 5' end on the rate of the acceptor splicing site utilization, or the differences in splicing rate between long, short, and RBM17-dependent introns. We also observed a deceleration of last intron splicing with an increase of the distance to the poly(A) site, which might be explained by the cooperativity of the splicing and polyadenylation. Additional analysis of splicing kinetics of SF3B4 knockdown cells suggested the impairment of a U2 snRNA recognition step. As a result, we deconvoluted the effects of several examined features on the splicing rate into a single regression model. The data obtained here are useful for further studies in the field, as they provide general splicing rate dependencies as well as help to justify the existence of slowly removed splice sites., (© 2024 Bolikhova et al.; Published by Cold Spring Harbor Laboratory Press for the RNA Society.)

U1 snRNP plays a critical role in 5'-splice site recognition and is a frequent target of alternative splicing factors. These factors transiently associate with human U1 snRNP and are not amenable for structural studies, while their Saccharomyces cerevisiae (yeast) homologs are stable components of U1 snRNP. Here, we report the cryoEM structure of yeast U1 snRNP at 3.6 Å resolution with atomic models for ten core proteins, nearly all essential domains of its RNA, and five stably associated auxiliary proteins. The foot-shaped yeast U1 snRNP contains a core in the "ball-and-toes" region architecturally similar to the human U1 snRNP. All auxiliary proteins are in the "arch-and-heel" region and connected to the core through the Prp42/Prp39 paralogs. Our demonstration that homodimeric human PrpF39 directly interacts with U1C-CTD, mirroring yeast Prp42/Prp39, supports yeast U1 snRNP as a model for understanding how transiently associated auxiliary proteins recruit human U1 snRNP in alternative splicing.

We report the clinical, biochemical, and molecular findings in two brothers with encephalopathy and multi-systemic disease. Abnormal transferrin glycoforms were suggestive of a type I congenital disorder of glycosylation (CDG). While exome sequencing was negative for CDG related candidate genes, the testing revealed compound heterozygous mutations in the mitochondrial elongation factor G gene (GFM1). One of the mutations had been reported previously while the second, novel variant was found deep in intron 6, activating a cryptic splice site. Functional studies demonstrated decreased GFM1 protein levels, suggested disrupted assembly of mitochondrial complexes III and V and decreased activities of mitochondrial complexes I and IV, all indicating combined OXPHOS deficiency.

A major goal of biomedicine is to understand the function of every gene in the human genome. Loss-of-function mutations can disrupt both copies of a given gene in humans and phenotypic analysis of such 'human knockouts' can provide insight into gene function. Consanguineous unions are more likely to result in offspring carrying homozygous loss-of-function mutations. In Pakistan, consanguinity rates are notably high. Here we sequence the protein-coding regions of 10,503 adult participants in the Pakistan Risk of Myocardial Infarction Study (PROMIS), designed to understand the determinants of cardiometabolic diseases in individuals from South Asia. We identified individuals carrying homozygous predicted loss-of-function (pLoF) mutations, and performed phenotypic analysis involving more than 200 biochemical and disease traits. We enumerated 49,138 rare (