DMD antisense oligonucleotide mediated exon skipping efficiency correlates with flanking intron retention time and target position within the exon.

Mutations in the DMD gene are causative for Duchenne muscular dystrophy (DMD). Antisense oligonucleotide (AON) mediated exon skipping to restore disrupted dystrophin reading frame is a therapeutic approach that allows production of a shorter but functional protein. As DMD causing mutations can affect most of the 79 exons encoding dystrophin, a wide variety of AONs are needed to treat the patient population. Design of AONs is largely guided by trial-and-error, and it is yet unclear what defines the skippability of an exon. Here, we use a library of phosphorodiamidate morpholino oligomer (PMOs) AONs of similar physical properties to test the skippability of a large number of DMD exons. The DMD transcript is non-sequentially spliced, meaning that certain introns are retained longer in the transcript than downstream introns. We tested whether the relative intron retention time has a significant effect on AON efficiency, and found that targeting an out-of-frame exon flanked at its 5'-end by an intron that is retained in the transcript longer ('slow' intron) leads to overall higher exon skipping efficiency than when the 5'-end flanking intron is 'fast'. Regardless of splicing speed of flanking introns, we find that positioning an AON closer to the 5'-end of the target exon leads to higher exon skipping efficiency opposed to targeting an exons 3'-end. The data enclosed herein can be of use to guide future target selection and preferential AON binding sites for both DMD and other disease amenable by exon skipping therapies.

SplicePie: a novel analytical approach for the detection of alternative, non-sequential and recursive splicing.

Alternative splicing is a powerful mechanism present in eukaryotic cells to obtain a wide range of transcripts and protein isoforms from a relatively small number of genes. The mechanisms regulating (alternative) splicing and the paradigm of consecutive splicing have recently been challenged, especially for genes with a large number of introns. RNA-Seq, a powerful technology using deep sequencing in order to determine transcript structure and expression levels, is usually performed on mature mRNA, therefore not allowing detailed analysis of splicing progression. Sequencing pre-mRNA at different stages of splicing potentially provides insight into mRNA maturation. Although the number of tools that analyze total and cytoplasmic RNA in order to elucidate the transcriptome composition is rapidly growing, there are no tools specifically designed for the analysis of nuclear RNA (which contains mixtures of pre- and mature mRNA). We developed dedicated algorithms to investigate the splicing process. In this paper, we present a new classification of RNA-Seq reads based on three major stages of splicing: pre-, intermediate- and post-splicing. Applying this novel classification we demonstrate the possibility to analyze the order of splicing. Furthermore, we uncover the potential to investigate the multi-step nature of splicing, assessing various types of recursive splicing events. We provide the data that gives biological insight into the order of splicing, show that non-sequential splicing of certain introns is reproducible and coinciding in multiple cell lines. We validated our observations with independent experimental technologies and showed the reliability of our method. The pipeline, named SplicePie, is freely available at: https://github.com/pulyakhina/splicing_analysis_pipeline. The example data can be found at: https://barmsijs.lumc.nl/HG/irina/example_data.tar.gz.

Non-sequential and multi-step splicing of the dystrophin transcript.

The dystrophin protein encoding DMD gene is the longest human gene. The 2.2 Mb long human dystrophin transcript takes 16 hours to be transcribed and is co-transcriptionally spliced. It contains long introns (24 over 10kb long, 5 over 100kb long) and the heterogeneity in intron size makes it an ideal transcript to study different aspects of the human splicing process. Splicing is a complex process and much is unknown regarding the splicing of long introns in human genes. Here, we used ultra-deep transcript sequencing to characterize splicing of the dystrophin transcripts in 3 different human skeletal muscle cell lines, and explored the order of intron removal and multi-step splicing. Coverage and read pair analyses showed that around 40% of the introns were not always removed sequentially. Additionally, for the first time, we report that non-consecutive intron removal resulted in 3 or more joined exons which are flanked by unspliced introns and we defined these joined exons as an exon block. Lastly, computational and experimental data revealed that, for the majority of dystrophin introns, multistep splicing events are used to splice out a single intron. Overall, our data show for the first time in a human transcript, that multi-step intron removal is a general feature of mRNA splicing.

A novel feed-forward loop between ARIH2 E3-ligase and PABPN1 regulates aging-associated muscle degeneration.

Alanine expansion mutations in poly(A)-binding protein nuclear 1 (PABPN1) cause muscle weakness in the late-onset disorder oculopharyngeal muscular dystrophy. In affected muscles, expanded PABPN1 forms nuclear aggregates, depleting levels of soluble PABPN1 and inducing a genome-wide shift from distal to proximal polyadenylation site usage. PABPN1 protein accumulation is regulated by the ubiquitin proteasome system, which is highly dysregulated in oculopharyngeal muscular dystrophy. We show that ARIH2 E3-ligase regulates PABPN1 protein accumulation and aggregation. Levels of ARIH2 mRNA are regulated by PABPN1 via proximal polyadenylation site usage. We demonstrate that masking the proximal polyadenylation site in ARIH2 3' untranslated region by antisense oligonucleotides elevates the expression of ARIH2 and PABPN1 and restores myogenic defects that are induced by ARIH2 or PABPN1 down-regulation in cell culture. In vivo ARIH2 mRNA levels significantly decrease from midlife in vastus lateralis muscles and highly correlate with muscle degeneration. We suggest that the expression of both genes is maintained by a feed-forward loop between mRNA stability regulated by PABPN1 and protein turnover regulated by ARIH2.

Longitudinal Assessment of Creatine Kinase, Creatine/Creatinineratio, and Myostatin as Monitoring Biomarkers in Becker Muscular Dystrophy.

BACKGROUND AND OBJECTIVES: The slow and variable disease progression of Becker muscular dystrophy (BMD) urges the development of biomarkers to facilitate clinical trials. We explored changes in 3 muscle-enriched biomarkers in serum of patients with BMD over 4-year time and studied associations with disease severity, disease progression, and dystrophin levels in BMD. METHODS: We quantitatively measured creatine kinase (CK) using the International Federation of Clinical Chemistry reference method, creatine/creatinineratio (Cr/Crn) using liquid chromatography-tandem mass spectrometry, and myostatin with ELISA in serum and assessed functional performance using the North Star Ambulatory Assessment (NSAA), 10-meter run velocity (TMRv), 6-Minute Walking Test (6MWT), and forced vital capacity in a 4-year prospective natural history study. Dystrophin levels were quantified in the tibialis anterior muscle using capillary Western immunoassay. The correlation between biomarkers, age, functional performance, mean annual change, and prediction of concurrent functional performance was analyzed using linear mixed models. RESULTS: Thirty-four patients with 106 visits were included. Eight patients were nonambulant at baseline. Cr/Crn and myostatin were highly patient specific (intraclass correlation coefficient for both = 0.960). Cr/Crn was strongly negatively correlated, whereas myostatin was strongly positively correlated with the NSAA, TMRv, and 6MWT (Cr/Crn rho = -0.869 to -0.801 and myostatin rho = 0.792 to 0.842, all p < 0.001). CK showed a negative association with age (p = 0.0002) but was not associated with patients' performance. Cr/Crn and myostatin correlated moderately with the average annual change of the 6MWT (rho = -0.532 and 0.555, p = 0.02). Dystrophin levels did not correlate with the selected biomarkers nor with performance. Cr/Crn, myostatin, and age could explain up to 75% of the variance of concurrent functional performance of the NSAA, TMRv, and 6MWT. DISCUSSION: Both Cr/Crn and myostatin could potentially serve as monitoring biomarkers in BMD, as higher Cr/Crn and lower myostatin were associated with lower motor performance and predictive of concurrent functional performance when combined with age. Future studies are needed to more precisely determine the context of use of these biomarkers.

Peripheral blood transcriptome profiling enables monitoring disease progression in dystrophic mice and patients.

DMD is a rare disorder characterized by progressive muscle degeneration and premature death. Therapy development is delayed by difficulties to monitor efficacy non-invasively in clinical trials. In this study, we used RNA-sequencing to describe the pathophysiological changes in skeletal muscle of 3 dystrophic mouse models. We show how dystrophic changes in muscle are reflected in blood by analyzing paired muscle and blood samples. Analysis of repeated blood measurements followed the dystrophic signature at five equally spaced time points over a period of seven months. Treatment with two antisense drugs harboring different levels of dystrophin recovery identified genes associated with safety and efficacy. Evaluation of the blood gene expression in a cohort of DMD patients enabled the comparison between preclinical models and patients, and the identification of genes associated with physical performance, treatment with corticosteroids and body measures. The presented results provide evidence that blood RNA-sequencing can serve as a tool to evaluate disease progression in dystrophic mice and patients, as well as to monitor response to (dystrophin-restoring) therapies in preclinical drug development and in clinical trials.

Antisense-Mediated Skipping of Dysferlin Exons in Control and Dysferlinopathy Patient-Derived Cells.

Dysferlinopathies encompass a spectrum of progressive muscular dystrophies caused by the lack of dysferlin due to missense mutations in the dysferlin gene or mutations causing premature truncation of protein translation. Dysferlin is a modular protein, and dysferlins lacking one or more repetitive domains have been shown to retain functionality. As such, antisense-mediated exon skipping has been proposed as a therapy for dysferlinopathy. By skipping the mutated exon, the reading frame would be maintained, while the mutation would be bypassed, thus allowing production of an internally deleted, but partially functional, dysferlin. We previously showed that dysferlin exon skipping is feasible in control cell lines. We here evaluated exon skipping and dysferlin protein restoration in patient-derived cells requiring the skipping of exon 9, 29, 30, or 34. Exon 30 skipping was possible at high efficiency, but did not result in increased dysferlin. We discovered that the alleged exon 30 mutation was in fact a polymorphism and identified a splicing mutation in intron 28 as the disease-causing mutation. While exon skipping was feasible for each of the other cell lines, no increases in dysferlin protein could be detected by western blotting.

Molecular signatures of age-associated chronic degeneration of shoulder muscles.

Chronic muscle diseases are highly prevalent in the elderly causing severe mobility limitations, pain and frailty. The intrinsic molecular mechanisms are poorly understood due to multifactorial causes, slow progression with age and variations between individuals. Understanding the underlying molecular mechanisms could lead to new treatment options which are currently limited. Shoulder complaints are highly common in the elderly, and therefore, muscles of the shoulder's rotator cuff could be considered as a model for chronic age-associated muscle degeneration. Diseased shoulder muscles were characterized by muscle atrophy and fatty infiltration compared with unaffected shoulder muscles. We confirmed fatty infiltration using histochemical analysis. Additionally, fibrosis and loss of contractile myosin expression were found in diseased muscles. Most cellular features, including proliferation rate, apoptosis and cell senescence, remained unchanged and genome-wide molecular signatures were predominantly similar between diseased and intact muscles. However, we found down-regulation of a small subset of muscle function genes, and up-regulation of extracellular region genes. Myogenesis was defected in muscle cell culture from diseased muscles but was restored by elevating MyoD levels. We suggest that impaired muscle functionality in a specific environment of thickened extra-cellular matrix is crucial for the development of chronic age-associated muscle degeneration.