Categorization of 77 dystrophin exons into 5 groups by a decision tree using indexes of splicing regulatory factors as decision markers

Ensemble-Learning and Feature Selection Techniques for Enhanced Antisense Oligonucleotide Efficacy Prediction in Exon Skipping

Antisense oligonucleotide (ASO)-mediated exon skipping has become a valuable tool for investigating gene function and developing gene therapy. Machine-learning-based computational methods, such as eSkip-Finder, have been developed to predict the efficacy of ASOs via exon skipping. However, these methods are computationally demanding, and the accuracy of predictions remains suboptimal. In this study, we propose a new approach to reduce the computational burden and improve the prediction performance by using feature selection within machine-learning algorithms and ensemble-learning techniques. We evaluated our approach using a dataset of experimentally validated exon-skipping events, dividing it into training and testing sets. Our results demonstrate that using a three-way-voting approach with random forest, gradient boosting, and XGBoost can significantly reduce the computation time to under ten seconds while improving prediction performance, as measured by R² for both 2'-O-methyl nucleotides (2OMe) and phosphorodiamidate morpholino oligomers (PMOs). Additionally, the feature importance ranking derived from our approach is in good agreement with previously published results. Our findings suggest that our approach has the potential to enhance the accuracy and efficiency of predicting ASO efficacy via exon skipping. It could also facilitate the development of novel therapeutic strategies. This study could contribute to the ongoing efforts to improve ASO design and optimize gene therapy approaches.

Analysis of Pathogenic Pseudoexons Reveals Novel Mechanisms Driving Cryptic Splicing

Understanding pre-mRNA splicing is crucial to accurately diagnosing and treating genetic diseases. However, mutations that alter splicing can exert highly diverse effects. Of all the known types of splicing mutations, perhaps the rarest and most difficult to predict are those that activate pseudoexons, sometimes also called cryptic exons. Unlike other splicing mutations that either destroy or redirect existing splice events, pseudoexon mutations appear to create entirely new exons within introns. Since exon definition in vertebrates requires coordinated arrangements of numerous RNA motifs, one might expect that pseudoexons would only arise when rearrangements of intronic DNA create novel exons by chance. Surprisingly, although such mutations do occur, a far more common cause of pseudoexons is deep-intronic single nucleotide variants, raising the question of why these latent exon-like tracts near the mutation sites have not already been purged from the genome by the evolutionary advantage of more efficient splicing. Possible answers may lie in deep intronic splicing processes such as recursive splicing or poison exon splicing. Because these processes utilize intronic motifs that benignly engage with the spliceosome, the regions involved may be more susceptible to exonization than other intronic regions would be. We speculated that a comprehensive study of reported pseudoexons might detect alignments with known deep intronic splice sites and could also permit the characterisation of novel pseudoexon categories. In this report, we present and analyse a catalogue of over 400 published pseudoexon splice events. In addition to confirming prior observations of the most common pseudoexon mutation types, the size of this catalogue also enabled us to suggest new categories for some of the rarer types of pseudoexon mutation. By comparing our catalogue against published datasets of non-canonical splice events, we also found that 15.7% of pseudoexons exhibit some splicing activity at one or both of their splice sites in non-mutant cells. Importantly, this included seven examples of experimentally confirmed recursive splice sites, confirming for the first time a long-suspected link between these two splicing phenomena. These findings have the potential to improve the fidelity of genetic diagnostics and reveal new targets for splice-modulating therapies.

Pseudoexons of the DMD Gene

The DMD gene is the largest in the human genome, with a total intron content exceeding 2.2Mb. In the decades since DMD was discovered there have been numerous reported cases of pseudoexons (PEs) arising in the mature DMD transcripts of some individuals, either as the result of mutations or as low-frequency errors of the spliceosome. In this review, I collate from the literature 58 examples of DMD PEs and examine the diversity and commonalities of their features. In particular, I note the high frequency of PEs that arise from deep intronic SNVs and discuss a possible link between PEs induced by distal mutations and the regulation of recursive splicing.

Neuronal SH-SY5Y cells use the C-dystrophin promoter coupled with exon 78 skipping and display multiple patterns of alternative splicing including two intronic insertion events

Duchenne muscular dystrophy (DMD) is a progressive muscle wasting disease caused by mutations in the dystrophin gene. One-third of DMD cases are complicated by mental retardation. Here, we used reverse transcription PCR to analyze the pattern of dystrophin transcripts in neuronal SH-SY5Y cells. Among the three alternative promoters/first exons at the 5'-end, only transcripts containing the brain cortex-specific C1 exon could be amplified. The C-transcript appeared as two products: a major product of the expected size and a minor larger product that contained the cryptic exon 1a between exons C1 and 2. At the 3'-end there was complete exon 78 skipping. Together, these findings indicate that SH-SY5Y cells have neuron-specific characteristics with regard to both promoter activation and alternative splicing. We also revealed partial skipping of exons 9 and 71. Four amplified products were obtained from a fragment covering exons 36-41: a strong expected product, two weak products lacking either exon 37 or exon 38, and a second strong larger product with a 568-bp insertion between exons 40 and 41. The inserted sequence matched the 3'-end of intron 40 perfectly. We concluded that a cryptic splice site was activated in SH-SY5Y cells to create the novel, unusually large, exon 41e (751 bp). In total, we identified seven alternative splicing events in neuronal SH-SY5Y cells, and calculated that 32 dystrophin transcripts could be produced. Our results may provide clues in the analysis of transcriptype-phenotype correlations as regards mental retardation in DMD.

Tissue- and case-specific retention of intron 40 in mature dystrophin mRNA

The dystrophin gene, which is mutated in Duchenne muscular dystrophy (DMD), comprises 79 exons that show multiple alternative splicing events. Intron retention, a type of alternative splicing, may control gene expression. We examined intron retention in dystrophin introns by reverse-transcription PCR from skeletal muscle, focusing on the nine shortest (all <1000 bp), because these are more likely to be retained. Only one, intron 40, was retained in mRNA; sequencing revealed insertion of a complete intron 40 (851 nt) between exons 40 and 41. The intron 40 retention product accounted for 1.2% of the total product but had a premature stop codon at the fifth intronic codon. Intron 40 retention was most strongly observed in the kidney (36.6%) and was not obtained from the fetal liver, lung, spleen or placenta. This indicated that intron retention is a tissue-specific event whose level varies among tissues. In two DMD patients, intron 40 retention was observed in one patient but not in the other. Examination of splicing regulatory factors revealed that intron 40 had the highest guanine-cytosine content of all examined introns in a 30-nt segment at its 3' end. Further studies are needed to clarify the biological role of intron 40-retained dystrophin mRNA.

Flexibility in the N-terminal actin-binding domain: clues from in silico mutations and molecular dynamics

Dystrophin is a long, rod-shaped cytoskeleton protein implicated in muscular dystrophy (MDys). Utrophin is the closest autosomal homolog of dystrophin. Both proteins have N-terminal actin-binding domain (N-ABD), a central rod domain and C-terminal region. N-ABD, composed of two calponin homology (CH) subdomains joined by a helical linker, harbors a few disease causing missense mutations. Although the two proteins share considerable homology (>72%) in N-ABD, recent structural and biochemical studies have shown that there are significant differences (including stability, mode of actin-binding) and their functions are not completely interchangeable. In this investigation, we have used extensive molecular dynamics simulations to understand the differences and the similarities of these two proteins, along with another actin-binding protein, fimbrin. In silico mutations were performed to identify two key residues that might be responsible for the dynamical difference between the molecules. Simulation points to the inherent flexibility of the linker region, which adapts different conformations in the wild type dystrophin. Mutations T220V and G130D in dystrophin constrain the flexibility of the central helical region, while in the two known disease-causing mutants, K18N and L54R, the helicity of the region is compromised. Phylogenetic tree and sequence analysis revealed that dystrophin and utrophin genes have probably originated from the same ancestor. The investigation would provide insight into the functional diversity of two closely related proteins and fimbrin, and contribute to our understanding of the mechanism of MDys.

In silico screening based on predictive algorithms as a design tool for exon skipping oligonucleotides in Duchenne muscular dystrophy

The use of antisense 'splice-switching' oligonucleotides to induce exon skipping represents a potential therapeutic approach to various human genetic diseases. It has achieved greatest maturity in exon skipping of the dystrophin transcript in Duchenne muscular dystrophy (DMD), for which several clinical trials are completed or ongoing, and a large body of data exists describing tested oligonucleotides and their efficacy. The rational design of an exon skipping oligonucleotide involves the choice of an antisense sequence, usually between 15 and 32 nucleotides, targeting the exon that is to be skipped. Although parameters describing the target site can be computationally estimated and several have been identified to correlate with efficacy, methods to predict efficacy are limited. Here, an in silico pre-screening approach is proposed, based on predictive statistical modelling. Previous DMD data were compiled together and, for each oligonucleotide, some 60 descriptors were considered. Statistical modelling approaches were applied to derive algorithms that predict exon skipping for a given target site. We confirmed (1) the binding energetics of the oligonucleotide to the RNA, and (2) the distance in bases of the target site from the splice acceptor site, as the two most predictive parameters, and we included these and several other parameters (while discounting many) into an in silico screening process, based on their capacity to predict high or low efficacy in either phosphorodiamidate morpholino oligomers (89% correctly predicted) and/or 2'O Methyl RNA oligonucleotides (76% correctly predicted). Predictions correlated strongly with in vitro testing for sixteen de novo PMO sequences targeting various positions on DMD exons 44 (R² 0.89) and 53 (R² 0.89), one of which represents a potential novel candidate for clinical trials. We provide these algorithms together with a computational tool that facilitates screening to predict exon skipping efficacy at each position of a target exon.

A novel splicing silencer generated by DMD exon 45 deletion junction could explain upstream exon 44 skipping that modifies dystrophinopathy

Duchenne muscular dystrophy (DMD), a progressive muscle-wasting disease, is mostly caused by exon deletion mutations in the DMD gene. The reading frame rule explains that out-of-frame deletions lead to muscle dystrophin deficiency in DMD. In outliers to this rule, deletion junction sequences have never previously been explored as splicing modulators. In a Japanese case, we identified a single exon 45 deletion in the patient's DMD gene, indicating out-of-frame mutation. However, immunohistochemical examination disclosed weak dystrophin signals in his muscle. Reverse transcription-PCR amplification of DMD exons 42 to 47 revealed a major normally spliced product with exon 45 deletion and an additional in-frame product with deletion of both exons 44 and 45, indicating upstream exon 44 skipping. We considered the latter to underlie the observed dystrophin expression. Remarkably, the junction sequence cloned by PCR walking abolished the splicing enhancer activity of the upstream intron in a chimeric doublesex gene pre-mRNA in vitro splicing. Furthermore, antisense oligonucleotides directed against the junction site counteracted this effect. These indicated that the junction sequence was a splicing silencer that induced upstream exon 44 skipping. It was strongly suggested that creation of splicing regulator is a modifier of dystrophinopathy.

Hidden dangers: a cryptic exon disrupts BRCA2 mRNA

The first mutation that disrupts BRCA2 mRNA by including a novel, cryptic exon is reported in this issue. The mutation lies deep within an intron and would not have been detected by conventional screening methods. In the future, more mutations may be discovered by direct mRNA analysis.