Splicing is controlled by a large set of regulatory elements (SREs) including splicing enhancers and silencers, which are involved in exon recognition. Variants at these motifs may dysregulate splicing and trigger loss-of-function transcripts associated with disease. Our goal here was to study the alternatively spliced exons 8 and 10 of the breast cancer susceptibility gene CHEK2. For this purpose, we used a previously published minigene with exons 6-10 that produced the expected minigene full-length transcript and replicated the naturally occurring events of exon 8 [Δ(E8)] and exon 10 [Δ(E10)] skipping. We then introduced 12 internal microdeletions of exons 8 and 10 by mutagenesis in order to map SRE-rich intervals by splicing assays in MCF-7 cells. We identified three minimal (10-, 11-, 15-nt) regions essential for exon recognition: c.863_877del [ex8, Δ(E8): 75%] and c.1073_1083del and c.1083_1092del [ex10, Δ(E10): 97% and 62%, respectively]. Then 87 variants found within these intervals were introduced into the wild-type minigene and tested functionally. Thirty-eight of them (44%) impaired splicing, four of which (c.883G>A, c.883G>T, c.884A>T, and c.1080G>T) induced negligible amounts (<5%) of the minigene full-length transcript. Another six variants (c.886G>A, c.886G>T, c.1075G>A, c.1075G>T, c.1076A>T, and c.1078G>T) showed significantly strong impacts (20-50% of the minigene full-length transcript). Thirty-three of the 38 spliceogenic variants were annotated as missense, three as nonsense, and two as synonymous, underlying the fact that any exonic change is capable of disrupting splicing. Moreover, c.883G>A, c.883G>T, and c.884A>T were classified as pathogenic/likely pathogenic variants according to ACMG/AMP (American College of Medical Genetics and Genomics/Association for Molecular Pathology)-based criteria. © 2024 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Fukutin (FKTN) c.647+2084G>T creates a pseudo-exon with a premature stop codon, which causes Fukuyama congenital muscular dystrophy (FCMD). We aimed to ameliorate aberrant splicing of FKTN caused by this variant. We screened compounds focusing on splicing regulation using the c.647+2084G>T splicing reporter and discovered that the branchpoint, which is essential for splicing reactions, could be a potential therapeutic target. To confirm the effectiveness of branchpoints as targets for exon skipping, we designed branchpoint-targeted antisense oligonucleotides (BP-AONs). This restored normal FKTN mRNA and protein production in FCMD patient myotubes. We identified a functional BP by detecting splicing intermediates and creating BP mutations in the FKTN reporter gene; this BP was non-redundant and sufficiently blocked by BP-AONs. Next, a BP-AON was designed for a different FCMD-causing variant, which induces pathogenic exon trapping by a common SINE-VNTR-Alu-type retrotransposon. Notably, this BP-AON also restored normal FKTN mRNA and protein production in FCMD patient myotubes. Our findings suggest that BPs could be potential targets in exon-skipping therapeutic strategies for genetic disorders.

The HLA genes are amongst the most polymorphic in the human genome. Alternative splicing could add an extra layer of complexity, but has not been studied extensively. Here, we applied an RNA based approach to study the influence of allele polymorphism on alternative splicing of HLA-C in peripheral blood. RNA was isolated from these peripheral cells, converted into cDNA and amplified specifically for 12 common HLA-C allele groups. Through subsequent sequencing of HLA-C, we observed alternative splicing variants of HLA-C*04 and *16 that resulted in exon 5 skipping and were co-expressed with the mature transcript. Investigation of intron 4 sequences of HLA-C*04 and *16 compared with other HLA-C alleles demonstrated no effect on predicted splice sites and branch point. To further investigate if the unique polymorphic positions in exon 5 of HLA-C*04 or *16 may facilitate alternative splicing by acting on splicing regulatory elements (SRE), in-silico splicing analysis was performed. While the HLA-C*04 specific SNP in exon 5 had no effect on predicted exonic SRE, the HLA-C*16 specific exon 5 SNP did alter exonic SRE. Our findings provide experimental and theoretical support for the concept that polymorphisms within the HLA-C alleles influence the alternative splicing of HLA-C.

The contribution of deep intronic splice-altering variants to hereditary breast and ovarian cancer (HBOC) is unknown. Current computational in silico tools to predict spliceogenic variants leading to pseudoexons have limited efficiency. We assessed the performance of the SpliceAI tool combined with ESRseq scores to identify spliceogenic deep intronic variants by affecting cryptic sites or splicing regulatory elements (SREs) using literature and experimental datasets. Our results with 233 published deep intronic variants showed that SpliceAI, with a 0.05 threshold, predicts spliceogenic deep intronic variants affecting cryptic splice sites, but is less effective in detecting those affecting SREs. Next, we characterized the SRE profiles using ESRseq, showing that pseudoexons are significantly enriched in SRE-enhancers compared to adjacent intronic regions. Although the combination of SpliceAI with ESRseq scores (considering ∆ESRseq and SRE landscape) showed higher sensitivity, the global performance did not improve because of the higher number of false positives. The combination of both tools was tested in a tumor RNA dataset with 207 intronic variants disrupting splicing, showing a sensitivity of 86%. Following the pipeline, five spliceogenic deep intronic variants were experimentally identified from 33 variants in HBOC genes. Overall, our results provide a framework to detect deep intronic variants disrupting splicing.

Discriminating which nucleotide variants cause disease or contribute to phenotypic traits remains a major challenge in human genetics. In theory, any intragenic variant can potentially affect RNA splicing by altering splicing regulatory elements (SREs). However, these alterations are often ignored mainly because pioneer SRE predictors have proved inefficient. Here, we report the first large-scale comparative evaluation of four user-friendly SRE-dedicated algorithms (QUEPASA, HEXplorer, SPANR, and HAL) tested both as standalone tools and in multiple combined ways based on two independent benchmark datasets adding up to >1,300 exonic variants studied at the messenger RNA level and mapping to 89 different disease-causing genes. These methods display good predictive power, based on decision thresholds derived from the receiver operating characteristics curve analyses, with QUEPASA and HAL having the best accuracies either as standalone or in combination. Still, overall there was a tight race between the four predictors, suggesting that all methods may be of use. Additionally, QUEPASA and HEXplorer may be beneficial as well for predicting variant-induced creation of pseudoexons deep within introns. Our study highlights the potential of SRE predictors as filtering tools for identifying disease-causing candidates among the plethora of variants detected by high-throughput DNA sequencing and provides guidance for their use in genomic medicine settings.

It is possible to estimate the prior probability of pathogenicity for germline disease gene variants based on bioinformatic prediction of variant effect/s. However, routinely used approaches have likely led to the underestimation and underreporting of variants located outside donor and acceptor splice site motifs that affect messenger RNA (mRNA) processing. This review presents information about hereditary cancer gene germline variants, outside native splice sites, with experimentally validated splicing effects. We list 95 exonic variants that impact splicing regulatory elements (SREs) in BRCA1, BRCA2, MLH1, MSH2, MSH6, and PMS2. We utilized a pre-existing large-scale BRCA1 functional data set to map functional SREs, and assess the relative performance of different tools to predict effects of 283 variants on such elements. We also describe rare examples of intronic variants that impact branchpoint (BP) sites and create pseudoexons. We discuss the challenges in predicting variant effect on BP site usage and pseudoexonization, and suggest strategies to improve the bioinformatic prioritization of such variants for experimental validation. Importantly, our review and analysis highlights the value of considering impact of variants outside donor and acceptor motifs on mRNA splicing and disease causation.

Accurate interpretation of genomic variants that alter RNA splicing is critical to precision medicine. We present a computational framework, Prediction of variant Effect on Percent Spliced In (PEPSI), that predicts the splicing impact of coding and noncoding variants for the Fifth Critical Assessment of Genome Interpretation (CAGI5) \"Vex-seq\" challenge. PEPSI is a random forest regression model trained on multiple layers of features associated with sequence conservation and regulatory sequence elements. Compared to other splicing defect prediction tools from the literature, our framework integrates secondary structure information in predicting variants that disrupt splicing regulatory elements (SREs). We applied our model to classify splice-disrupting variants among 2,094 single-nucleotide polymorphisms from the Exome Aggregation Consortium using model-predicted changes in percent spliced in (ΔPSI) associated with tested variants. Benchmarking our model against widely used state-of-the-art tools, we demonstrate that PEPSI achieves comparable performance in terms of sensitivity and precision. Moreover, we also show that using secondary structure context can help resolve several cases where changes in the counts of SREs do not correspond with the directionality of ΔPSI measured for tested variants.

Transcription of the HIV-1 provirus generates a viral pre-mRNA, which is alternatively spliced into more than 50 HIV-1 mRNAs encoding all viral proteins. Regulation of viral alternative splice site usage includes the presence of splicing regulatory elements (SREs) which can dramatically impact RNA expression and HIV-1 replication when mutated. Recently, we were able to show that two viral SREs, GI3-2 and ESEtat, are important players in the generation of viral vif, vpr and tat mRNAs. Furthermore, we demonstrated that masking these SREs by transfected locked nucleic acid (LNA) mixmers affect the viral splicing pattern and viral particle production. With regard to the development of future therapeutic LNA mixmer-based antiretroviral approaches, we delivered the GI3-2 and the ESEtat LNA mixmers \"nakedly\", without the use of transfection reagents (gymnosis) into HIV-1 infected cells. Surprisingly, we observed that gymnotically-delivered LNA mixmers accumulated in the cytoplasm, and seemed to co-localize with GW bodies and induced degradation of mRNAs containing their LNA target sequence. The GI3-2 and the ESEtat LNA-mediated RNA degradation resulted in abrogation of viral replication in HIV-1 infected Jurkat and PM1 cells as well as in PBMCs.

Alternative splicing plays a key role in the HIV-1 life cycle and is essential to maintain an equilibrium of mRNAs that encode viral proteins and polyprotein-isoforms. In particular, since all early HIV-1 proteins are expressed from spliced intronless and late enzymatic and structural proteins from intron containing, i.e. splicing repressed viral mRNAs, cellular splicing factors and splicing regulatory proteins are crucial for the replication capacity. In this review, we will describe the complex network of cis-acting splicing regulatory elements (SREs), which are mainly localized in the neighbourhoods of all HIV-1 splice sites and warrant the proper ratio of individual transcript isoforms. Since SREs represent binding sites for trans-acting cellular splicing factors interacting with the cellular spliceosomal apparatus we will review the current knowledge of interactions between viral RNA and cellular proteins as well as their impact on viral replication. Finally, we will discuss potential therapeutic approaches targeting HIV-1 alternative splicing.

This issue dedicated to the code of life tackles very challenging and open questions in Biology. The genetic code, brilliantly uncovered over 50 years ago is an example of a univocal biological code. In fact, except for very few and marginal variations, it is the same from bacteria to man, the RNA stretch: 5\' GUGUUC 3\' reads as the dipeptide: Val-Phe in bacteria, in yeast, in Arabidopsis, in zebra fish, in mouse and in human. A degree of ambiguity is possible if mutations are introduced in the tRNAs in a way that the anticodon reads one amino acid but the aminoacyl-transferase attaches a different one onto the tRNA. These were the very useful suppressor genes that aided greatly the study of bacterial genetics. Other biological codes however, are more akin to social codes and are less amenable to an unambiguous deciphering. Legal and ethical codes, weather we like it or not, are flexible and depend on the structure and history of the society that has produced them, as well as a specific point in time. The codes that govern RNA splicing have similar characteristics. In fact, the splicing code depends on a myriad of different factors that in part are influenced by the background in which they are read such as different cells, tissues or developmental stages. Given the complexity of the splicing process, the construction of an algorithm that can define exons or their fate with certainty has not yet been achieved. However a substantial amount of information towards the deciphering of the splicing code has been gathered and in this manuscript we summarize the point reached.