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 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.

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.

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.

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.

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.

Classification of introns, which is crucial to understanding their evolution and splicing, has historically been binary and has resulted in the naming of major and minor introns that are spliced by their namesake spliceosome. However, a broad range of intron consensus sequences exist, leading us to here reclassify introns as minor, minor-like, hybrid, major-like, major and non-canonical introns in 263 species across six eukaryotic supergroups. Through intron orthology analysis, we discovered that minor-like introns are a transitory node for intron conversion across evolution. Despite close resemblance of their consensus sequences to minor introns, these introns possess an AG dinucleotide at the -1 and -2 position of the 5' splice site, a salient feature of major introns. Through combined analysis of CoLa-seq, CLIP-seq for major and minor spliceosome components, and RNAseq from samples in which the minor spliceosome is inhibited we found that minor-like introns are also an intermediate class from a splicing mechanism perspective. Importantly, this analysis has provided insight into the sequence elements that have evolved to make minor-like introns amenable to recognition by both minor and major spliceosome components. We hope that this revised intron classification provides a new framework to study intron evolution and splicing., (© The Author(s) 2024. Published by Oxford University Press on behalf of Nucleic Acids Research.)

pre-mRNA splicing is a critical feature of eukaryotic gene expression. Both cis- and trans-splicing rely on accurately recognising splice site sequences by spliceosomal U snRNAs and associated proteins. Spliceosomal snRNAs carry multiple RNA modifications with the potential to affect different stages of pre-mRNA splicing. Here, we show that the conserved U6 snRNA m6A methyltransferase METT-10 is required for accurate and efficient cis- and trans-splicing of C. elegans pre-mRNAs. The absence of METT-10 in C. elegans and METTL16 in humans primarily leads to alternative splicing at 5' splice sites with an adenosine at +4 position. In addition, METT-10 is required for splicing of weak 3' cis- and trans-splice sites. We identified a significant overlap between METT-10 and the conserved splicing factor SNRNP27K in regulating 5' splice sites with +4A. Finally, we show that editing endogenous 5' splice site +4A positions to +4U restores splicing to wild-type positions in a mett-10 mutant background, supporting a direct role for U6 snRNA m6A modification in 5' splice site recognition. We conclude that the U6 snRNA m6A modification is important for accurate and efficient pre-mRNA splicing., (© The Author(s) 2024. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Alternative splicing plays a crucial role in protein diversity and gene expression regulation in higher eukaryotes, and mutations causing dysregulated splicing underlie a range of genetic diseases. Computational prediction of alternative splicing from genomic sequences not only provides insight into gene-regulatory mechanisms but also helps identify disease-causing mutations and drug targets. However, the current methods for the quantitative prediction of splice site usage still have limited accuracy. Here, we present DeltaSplice, a deep neural network model optimized to learn the impact of mutations on quantitative changes in alternative splicing from the comparative analysis of homologous genes. The model architecture enables DeltaSplice to perform "reference-informed prediction" by incorporating the known splice site usage of a reference gene sequence to improve its prediction on splicing-altering mutations. We benchmarked DeltaSplice and several other state-of-the-art methods on various prediction tasks, including evolutionary sequence divergence on lineage-specific splicing and splicing-altering mutations in human populations and neurodevelopmental disorders, and demonstrated that DeltaSplice outperformed consistently. DeltaSplice predicted ∼15% of splicing quantitative trait loci (sQTLs) in the human brain as causal splicing-altering variants. It also predicted splicing-altering de novo mutations outside the splice sites in a subset of patients affected by autism and other neurodevelopmental disorders (NDDs), including 19 genes with recurrent splicing-altering mutations. Integration of splicing-altering mutations with other types of de novo mutation burdens allowed the prediction of eight novel NDD-risk genes. Our work expanded the capacity of in silico splicing models with potential applications in genetic diagnosis and the development of splicing-based precision medicine., (© 2024 Xu et al.; Published by Cold Spring Harbor Laboratory Press.)

PARP inhibitor (PARPi) therapy has transformed outcomes for patients with homologous recombination DNA repair (HRR) deficient ovarian cancers, for example those with BRCA1 or BRCA2 gene defects. Unfortunately, PARPi resistance is common. Multiple resistance mechanisms have been described, including secondary mutations that restore the HR gene reading frame. BRCA1 splice isoforms △11 and △11q can contribute to PARPi resistance by splicing out the mutation-containing exon, producing truncated, partially functional proteins. However, the clinical impacts and underlying drivers of BRCA1 exon skipping are not fully understood.We analyzed nine ovarian and breast cancer patient derived xenografts (PDX) with BRCA1 exon 11 frameshift mutations for exon skipping and therapy response, including a matched PDX pair derived from a patient pre- and post-chemotherapy/PARPi. BRCA1 exon 11 skipping was elevated in PARPi resistant PDX tumors. Two independent PDX models acquired secondary BRCA1 splice site mutations (SSMs) that drive exon skipping, confirmed using qRT-PCR, RNA sequencing, immunoblotting and minigene modelling. CRISPR/Cas9-mediated disruption of splicing functionally validated exon skipping as a mechanism of PARPi resistance. SSMs were also enriched in post-PARPi ovarian cancer patient cohorts from the ARIEL2 and ARIEL4 clinical trials.Few PARPi resistance mechanisms have been confirmed in the clinical setting. While secondary/reversion mutations typically restore a gene's reading frame, we have identified secondary mutations in patient cohorts that hijack splice sites to enhance mutation-containing exon skipping, resulting in the overexpression of BRCA1 hypomorphs, which in turn promote PARPi resistance. Thus, BRCA1 SSMs can and should be clinically monitored, along with frame-restoring secondary mutations., (© 2024. The Author(s).)

Lowering expression of prion protein (PrP) is a well-validated therapeutic strategy in prion disease, but additional modalities are urgently needed. In other diseases, small molecules have proven capable of modulating pre-mRNA splicing, sometimes by forcing inclusion of cryptic exons that reduce gene expression. Here, we characterize a cryptic exon located in human PRNP's sole intron and evaluate its potential to reduce PrP expression through incorporation into the 5' untranslated region. This exon is homologous to exon 2 in nonprimate species but contains a start codon that would yield an upstream open reading frame with a stop codon prior to a splice site if included in PRNP mRNA, potentially downregulating PrP expression through translational repression or nonsense-mediated decay. We establish a minigene transfection system and test a panel of splice site alterations, identifying mutants that reduce PrP expression by as much as 78%. Our findings nominate a new therapeutic target for lowering PrP., Competing Interests: Conflict of interest S. M. V. acknowledges speaking fees from Ultragenyx, Illumina, Biogen, and Eli Lilly, consulting fees from Invitae and Alnylam, and research support from Ionis, Gate, and Sangamo. E.V.M. acknowledges speaking fees from Eli Lilly, consulting fees from Deerfield and Alnylam, and research support from Ionis, Gate, and Sangamo., (Copyright © 2024 The Authors. Published by Elsevier Inc. All rights reserved.)

Congenital contractural arachnodactyly (CCA) is a rare connective tissue disorder characterized by arachnodactyly, multiple joint contractures, progressive kyphoscoliosis, pectus deformity and abnormal crumpled ears. FBN2 is the only gene currently known to be associated with CCA. In this study, we report on a prenatal case presented with skeletal, cardiac and spinal malformations. And his father had elongated limbs, contractures of the proximal interphalangeal joints, high myopia and scoliosis. We conducted whole exome sequencing (WES) on the fetus-parental trio and a heterozygous variant (hg19 chr5:127,673,685, c.3598 + 4A > G, NM_001999.4) in intron 27 of the FBN2 gene was successfully identified, inherited from the father. Reverse transcriptase-polymerase chain reaction (RT-PCR) was performed to evaluate the potential splicing effect of this variant, which confirmed that the variant caused a deletion of exon 27 (126 bp) by disrupting the splice-donor site and destroyed the 17th calcium-binding epidermal growth factor-like (cbEGF) domain. Our research not only finds the etiology of the disease in affected individuals and expands the mutation spectrum of FBN2 gene, but also provides genetic counseling and fertility guidance for this family., (© 2023. The Author(s).)

Intronic deletions that critically shorten donor-to-branchpoint (D-BP) distance of a precursor mRNA impose biophysical space constraint on assembly of the U1/U2 spliceosomal complex, leading to canonical splicing failure. Here we use a series of β-globin (HBB) gene constructs with intron 1 deletions to define D-BP lengths that present low/no risk of mis-splicing and lengths which are critically short and likely elicit clinically relevant mis-splicing. We extend our previous observation in EMD intron 5 of 46 nt as the minimal productive D-BP length, demonstrating spliceosome assembly constraint persists at D-BP lengths of 47-56 nt. We exploit the common HBB exon 1 β-thalassemia variant that strengthens a cryptic donor (NM&#95;000518.5(HBB):c.79G > A) to provide a simple barometer for the earliest signs of space constraint, via cryptic donor activation. For clinical evaluation of intronic deletions, we assert D-BP lengths > 60 nt present low mis-splicing risk while space constraint increases exponentially with D-BP lengths < 55 nt, with critical risk and profound splicing abnormalities with D-BP lengths < 50 nt., (© 2024. The Author(s).)

Identification of splice sites is a critical step in pre-messenger RNA (pre-mRNA) splicing because the definition of the exon/intron boundaries controls what nucleotides are incorporated into mature mRNAs. The intron boundary with the upstream exon is initially identified through interactions with the U1 small nuclear ribonucleoprotein (snRNP). This involves both base-pairing between the U1 snRNA and the pre-mRNA as well as snRNP proteins interacting with the 5' splice site (5'ss)/snRNA duplex. In yeast, this duplex is buttressed by two conserved protein factors, Yhc1 and Luc7. Luc7 has three human paralogs (LUC7L, LUC7L2, and LUC7L3), which play roles in alternative splicing. What domains of these paralogs promote splicing at particular sites is not yet clear. Here, we humanized the zinc finger (ZnF) domains of the yeast Luc7 protein in order to understand their roles in splice site selection using reporter assays, transcriptome analysis, and genetic interactions. Although we were unable to determine a function for the first ZnF domain, humanization of the second ZnF domain to mirror that found in LUC7L or LUC7L2 resulted in altered usage of nonconsensus 5'ss. In contrast, the corresponding ZnF domain of LUC7L3 could not support yeast viability. Further, humanization of Luc7 can suppress mutation of the ATPase Prp28, which is involved in U1 release and exchange for U6 at the 5'ss. Our work reveals a role for the second ZnF of Luc7 in splice site selection and suggests that different ZnF domains may have different ATPase requirements for release by Prp28., (© 2024 Carrocci et al.; Published by Cold Spring Harbor Laboratory Press for the RNA Society.)

The recognition of the 5' splice site (5' ss) is one of the earliest steps of pre-mRNA splicing. To better understand, the mechanism and regulation of 5' ss recognition, we selectively humanized components of the yeast U1 (yU1) snRNP to reveal the function of these components in 5' ss recognition and splicing. We targeted U1C and Luc7, two proteins that interact with and stabilize the yU1 snRNA and the 5' ss RNA duplex. We replaced the zinc-finger (ZnF) domain of yeast U1C (yU1C) with its human counterpart, which resulted in a cold-sensitive growth phenotype and moderate splicing defects. We next added an auxin-inducible degron to yeast Luc7 (yLuc7) protein (to mimic the lack of Luc7Ls in human U1 snRNP). We found that Luc7-depleted yU1 snRNP resulted in the concomitant loss of Prp40 and Snu71 (two other essential yU1 snRNP proteins), and further biochemical analyses suggest a model of how these three proteins interact with each other in the U1 snRNP. The loss of these proteins resulted in a significant growth retardation accompanied by a global suppression of pre-mRNA splicing. The splicing suppression led to mitochondrial dysfunction as revealed by a release of Fe 2+ into the growth medium and an induction of mitochondrial reactive oxygen species. Together, these observations indicate that the human U1C ZnF can substitute that of yeast, Luc7 is essential for the incorporation of the Luc7-Prp40-Snu71 trimer into yU1 snRNP, and splicing plays a major role in the regulation of mitochondrial function in yeast., (© 2024 Chalivendra et al.; Published by Cold Spring Harbor Laboratory Press for the RNA Society.)

Bioinformatics is a contemporary interdisciplinary area focused on analyzing the growing number of genome sequences. Gene variants are differences in DNA sequences among individuals within a population. Splice site recognition is a crucial step in the process of gene expression, where the coding sequences of genes are joined together to form mature messenger RNA (mRNA). These genetic variants that disrupt genes are believed to be the primary reason for neuro-developmental disorders like ASD (Autism Spectrum Disorder) is a neuro-developmental disorder that is diagnosed in individuals, families, and society and occurs as the developmental delay in one among the hundred genes that are associated with these disorders. Missense variants, premature stop codons, or deletions alter both the quality and quantity of encoded proteins. Predicting genes within exons and introns presents main challenges, such as dealing with sequencing errors, short reads, incomplete genes, overlapping, and more. Although many traditional techniques have been utilized in creating an exon prediction system, the primary challenge lies in accurately identifying the length and spliced strand location classification of exons in conjunction with introns. From now on, the suggested approach utilizes a Deep Learning algorithm to analyze intricate and extensive genomic datasets. M-LSTM is utilized to categorize three binary combinations (EI as 1, IE as 2, and none as 3) using spliced DNA strands. The M-LSTM system is able to sequence extensive datasets, ensuring that long information can be stored without any impact on the current input or output. This enables it to recognize and address long-term connections and problems with rapidly increasing gradients. The proposed model is compared internally with Naïve Bayes and Random Forest to assess its efficacy. Additionally, the proposed model's performance is forecasted by utilizing probabilistic parameters like recall, F1-score, precision, and accuracy to assess the effectiveness of the proposed system., 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 Elsevier B.V. All rights reserved.)

CRISPR-Cas technology has transformed functional genomics, yet understanding of how individual exons differentially shape cellular phenotypes remains limited. Here, we optimized and conducted massively parallel exon deletion and splice-site mutation screens in human cell lines to identify exons that regulate cellular fitness. Fitness-promoting exons are prevalent in essential and highly expressed genes and commonly overlap with protein domains and interaction interfaces. Conversely, fitness-suppressing exons are enriched in nonessential genes, exhibiting lower inclusion levels, and overlap with intrinsically disordered regions and disease-associated mutations. In-depth mechanistic investigation of the screen-hit TAF5 alternative exon-8 revealed that its inclusion is required for assembly of the TFIID general transcription initiation complex, thereby regulating global gene expression output. Collectively, our orthogonal exon perturbation screens established a comprehensive repository of phenotypically important exons and uncovered regulatory mechanisms governing cellular fitness and gene expression., Competing Interests: Declaration of interests The authors declare no competing interests., (Published by Elsevier Inc.)

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.