Abstract
Background Duchenne muscular dystrophy (DMD) is caused by DMD mutations leading to dystrophin loss. Full length Dp427 is the primary dystrophin isoform expressed in skeletal muscle and is also expressed in the central nervous system (CNS). Two shorter isoforms, Dp140 and Dp71, are highly expressed in the CNS. While a role for Dp140 and Dp71 on DMD CNS co-morbidities is well known, relationships between lack of Dp140 and Dp71 and DMD motor outcomes are not. We have conducted a series of investigations addressing this.
Methods Functional outcome data from 387 DMD boys aged 4.0-15.4 years was subdivided by DMD mutation expected effect on isoform expression; Group 1 (Dp427 absent, Dp140/Dp71 present, n=201); group 2 (Dp427/Dp140 absent, Dp71 present, n=152); and group 3 (Dp427/Dp140/Dp71 absent, n=34). Relationships between isoform group and North Star ambulatory assessment (NSAA) scores, 10m walk/run and rise times were explored using regression analysis. We used Capillary Western immunoassay (Wes) analysis to study Dp427, Dp140 and Dp71 production in wild-type and DMD skeletal muscle and myogenic cultures. Grip strength was studied in wild-type, mdx (Dp427 absent, Dp140/Dp71 present), mdx52 (Dp427/Dp140 absent, Dp71 present) and DMD-null (lacking all isoforms) mice.
Results In DMD boys, we found a strong association between isoform group and motor function. In DMD boys, mean NSAA scores at 5 years of age were 6.1 points lower in group 3 than group 1 (p<0.01) and 4.9 points lower in group 3 than group 2 (p=0.05). Mean peak NSAA scores were 4.0 points lower in group 3 than group 1 (p<0.01), 2.4 points lower in group 3 than group 2 (p=0.09) and 1.6 points lower in group 2 than group 1 (p=0.04).
Average grip strength in peak force at 3 months of age was higher in mdx than mdx52 mice (p=0.01).
Dp427, but not Dp71, was produced in normal skeletal muscle; low levels of Dp71 were detected in DMD skeletal muscle. High Dp71 levels were present in wild-type and DMD myogenic cultures.
Conclusions DMD boys lacking Dp140 and Dp140/Dp71 displayed worse motor function with a cumulative effect of isoform loss. DMD mouse models lacking Dp427 and Dp140 had lower grip strength than those lacking Dp427 but not Dp140. Our results highlight the importance of considering the effects of dystrophin isoform loss on DMD motor impairment, with important implications for understanding the complex relationship between brain and muscle function in DMD and patient stratification for clinical trials.
Introduction
Duchenne muscular dystrophy (DMD) is an X-linked recessive condition caused by DMD mutations leading to deficiency of dystrophin[1,2]
Boys with DMD typically present with delayed motor milestones, frequent falls and speech delay[2]. Neurobehavioural comorbidities, including intellectual disability, attention deficit hyperactivity disorder and/or autism spectrum disorder, occur in approximately one third of patients[2]. Loss of ambulation typically occurs by 12 years of age, followed by progression to cardiomyopathy and respiratory insufficiency[2]. Mean age of survival is in the late twenties[2].
DMD shows significant clinical heterogeneity. Some heterogeneity is related to DMD genotypes allowing production of low levels of dystrophin[3]. In addition, changes in genes other than DMD, that modify the clinical course, have recently been identified[4–6]. However, each of these modifiers contributes only modestly to DMD trajectory, suggesting the existence of other factors.
The DMD locus encodes for multiple isoforms[7]. Full-length Dp427 exists in 3 almost identical isoforms. Dp427m, is the isoform most highly expressed in skeletal muscle. Dp427c and Dp427p are expressed in the cortex, amygdala, hippocampus and cerebellum[7]. Additional shorter dystrophin isoforms are driven by promoters in introns further downstream the gene. Two of which, Dp140 and Dp71, are abundantly expressed in the brain[7,8]. Whilst in DMD, all mutations affect Dp427 production, the DMD mutation can also result in disruption of one or both of Dp140 and Dp71. DMD mutations involving the genomic region upstream of intron 44 only affect Dp427; mutations involving the genomic region from exon 51 to exon 62 inclusive, but not involving the genomic region of exon 63 and/or the genomic region downstream of exon 63, also affect Dp140, but not Dp71, and mutations involving the genomic region of exon 63 and/or the genomic region downstream of exon 63 affect Dp427, Dp140 and Dp71[7,8]. The effect of mutations involving the genomic region from exon 45 to exon 50 inclusive and not involving the genomic region of exon 51 and/or the genomic region downstream of exon 51 on Dp140 production cannot be predicted precisely, as the Dp140 promoter is located in intron 44 and its translation start site is located in exon 51[8,9]. Other isoforms also exist (Dp260 and Dp116), with patterns of expression limited to the retina and peripheral nerves respectively[7].
Several studies have suggested a role of Dp140 and Dp71 in DMD central nervous system (CNS) involvement[10–13]. In a European DMD cohort, 15% of boys lacking only Dp427 had intellectual disability, compared with 25% of boys lacking Dp427 and Dp140 and 64% of boys lacking Dp427, Dp140 and Dp71[13].
We hypothesised that DMD mutations affecting different dystrophin isoforms may have a different impact on DMD motor outcomes. We conducted a series of studies to address this. We evaluated the relationship between mutations leading to loss of Dp427/Dp140 and Dp427/Dp140/Dp71 compared to loss of Dp427 alone on motor function in a large cohort of DMD boys. We focused on four motor outcomes commonly used in clinical trials: (1) North Star ambulatory assessment (NSAA) score, (2) rise from supine time and velocity, (3) 10 metre walk/run time and velocity, and (4) age at loss of ambulation (LOA) in boys with DMD. We evaluated grip strength in wild-type and 3 DMD mouse models with mutations associated with different patterns of dystrophin isoform expression. Finally, we evaluated dystrophin isoform protein production in wild type and DMD skeletal muscle and myogenic cells in culture.
Methods
Standard protocol, approvals, registrations, and patient consents
The Dubowitz Neuromuscular Centre coordinates the North Star network of 23 UK neuromuscular centres looking after DMD patients[13–15]. Clinical data is prospectively collected in routine clinical appointments every six months. There is a national training programme and standard operating procedures for each site. Clinical data is stored in an electronic database managed by CertusLtd[13]. Participants meeting inclusion criteria were recruited from this network.
Study design and participants
387 participants aged 4.0-15.4 years, with a genetically confirmed DMD diagnosis or a diagnosis of intermediate muscular dystrophy (IMD) and an out-of-frame DMD deletion or duplication were included. Participants with in-frame DMD deletions/duplications and an IMD phenotype, BMD, manifesting carriers, clinical trial participants and boys under 4 years were excluded.
NorthStar natural history study data
Clinical data reviewed from the NorthStar natural history study included DMD mutation, date of DMD diagnosis, NSAA scores, rise from supine time, 10 metre walk/run time, age of LOA and presence or absence of learning disabilities. The NSAA is a 17-item DMD-specific motor function scale with a maximum possible total score of 34[16].
For the functional data, we analysed data at 5 years of age, as this is the standard lower end of age at recruitment in most clinical trials; the age of peak functional achievement, and the age at loss of ambulation.
DMD mutation data
Participants were grouped into 3 groups based on predicted DMD mutation effects on dystrophin isoform expression: Group 1 (Dp427 absent, Dp140/Dp71 present, n=201); group 2 (Dp427/Dp140 absent, Dp71 present, n=152); and group 3 (Dp427/Dp140/Dp71 absent, n=34). Based on the genomic organisation of the DMD locus[7,8], patients with mutations only involving the genomic region upstream of intron 44 were considered to be Dp427 negative and Dp140/Dp71 positive (group 1), patients with mutations involving the genomic region from exon 51 to exon 62 inclusive, but not involving the genomic region of exon 63 and/or the genomic region downstream of exon 63 were considered Dp427/Dp140 negative and Dp71 positive (group 2) and patients with mutations involving the genomic region of exon 63 and/or the genomic region downstream of exon 63 were considered to be Dp427/Dp140/Dp71 negative (group 3)[8,9]. Patients with mutations involving the genomic region from exon 45 to exon 50 inclusive and not involving the genomic region of exon 51 and/or the genomic region downstream of exon 51 were excluded from analysis due to the difficulty in predicting the effects of these mutations on Dp140 expression, as the Dp140 promoter is located in intron 44 and its translation start site is located in exon 51[8,9].
Cognitive status
Cognitive status was evaluated in a subset of boys by both formal testing of intelligence quotient (IQ) and presence or absence of learning difficulties. Presence or absence of learning difficulties was collected for the NorthStar natural history study and determined by parental and/or educational report.
Formal IQ testing was carried out in a subset of 40 boys who have participated in other studies undertaken at our centre - Ricotti et al[10] and Maresh et al (manuscript in preparation). Methods of IQ assessments are reported in our previous publication[10].
Participants were categorised into two groups based on cognitive status. The cognitive impairment group (n=56) had learning disabilities (as reported by parent/guardian(s) and/or educational services) and/or an IQ of less than 85 (>1 standard deviation below the mean population IQ)[17,18]. The normal cognition group (n=71) had no learning disability (as reported by parent/guardian(s) and/or educational services) and/or an IQ of 85 or above.
Statistical analysis for analysis in DMD boys
Patient characteristics were summarised using mean and standard deviation for continuous data and frequency and proportion for categorical data. GC regimen was summarised as the regime (Daily, Intermittent/Other or None) taken for the longest duration over both the full longitudinal period (majority regimen) and the duration prior to peak NSAA (early regime). Where age of GC initiation was not available (139/340 of boys on GC), this was estimated based on mean age of GC initiation in the cohort. Ambulatory function at 5 years was described for patients for whom a visit was recoded between 4.5-5.5 years of age.
Comparisons between isoform groups and neurocognitive impairment groups were made using one-way analysis of variance and Chi-squared tests. For the 3 main longitudinal outcomes, NSAA total score, 10m walk/run velocity and rise time velocity, relationships with age were explored for each isoform group; additionally, the relationship for NSAA total score with age was explored in the 2 cognition groups. Fractional polynomial regression, accounting for the longitudinal data, was used to find best fitting models, by comparing model deviances22. Noting the different observed mean peak functional scores for the dystrophin isoform groups and neurocognitive groups, we focused subsequent analysis on the age range 5.5–8 years, where we expected the majority of participants to reach peak motor function. Within this age interval, we calculated maximum observed motor function scores and age of maximum attainment for the 3 main outcomes.
Multivariable regression analysis was used to compare maximum scores achieved between dystrophin isoform groups and neurocognitive groups. We assessed whether differences in the majority GC regimen used up to the time of maximum score accounted for differences between dystrophin isoform groups. Kaplan-Meier survival estimation was used to estimate median time to LOA for each dystrophin isoform group. P value less than or equal to 0.05 was considered statistically significant.
Statistical analysis and creation of figures was conducted by Dr Deborah Ridout in Stata v 15 StataCorp. 2017. Stata Statistical Software: Release 15. College Station, TX: StataCorp LLC
Mouse DMD models
All mice used in this study were maintained at the National Center of Neurology and Psychiatry (NCNP)[19]. 3 different dystrophic mouse models were studied, with differential dystrophin isoform involvement; mdx mice (Dp427 absent, Dp140/Dp71 present), mdx-52 mice (Dp427/Dp140 absent, Dp71 present) and DMD-null mice lacking all dystrophin isoforms[20–22]. Mdx52 mice were kindly provided by Dr T. Sasaoka (Brain Research Institute, Niigata University, Niigata, Japan)[23]. DMD-null mice were generated by Dr Kazunori Hanaoka[22]. C57BL/6 mice were used as controls to match the background of mdx, mdx52 and DMD-null mice. Genotyping was performed using previously described PCR method[21–23]. Animal care was provided by the Small Animal Research Facility at the NCNP. Mice were allowed ad libitum access to food and drinking water. All behavioural experiments were performed between 9 am and 1 pm in strict accordance with the regulations of the National Institute of Neuroscience and the National Center of Neurology and Psychiatry (Japan) for animal experiments and were approved by the Animal Investigation Committee of the Institute. Statistical analysis and figures were carried out using GraphPad Prism8 (GraphPad Software Inc., La Jolla, CA).
Mouse forelimb grip strength test
The grip strength test was conducted on male dystrophic mouse models, including mdx (n=9), mdx52 (n=10), DMD-null (n=3) mice at 3 months of age. Mice were tested to determine peak paw grip strength, positioned horizontally from a grip bar using a grip strength meter (MK-380 M; Muromachi Kikai Co., Ltd., Japan) as previously described, and pulled back slowly and steadily until these mice released their grip[24]. This was repeated 6 times, and peak force for the forelimb paws was measured. These mice weighed on the date of their grip strength testing to be normalized by their body weight. The degree of grip fatigue (Force-decline rate) is calculated by comparing the first two pull to the last two pull. In the formula (5th +6th) / (1st +2nd) *100 (%) gives a measure of fatigue.
Dystrophin isoform protein production in wild type and DMD skeletal muscle and myogenic cells
All patients and controls provided written informed consent for skin and muscle biopsy samples. Fibroblasts and muscle tissues were supplied by the MRC Centre for Neuromuscular Disease Biobank London.
Tissue sampling was approved by the NHS national Research Ethics Service as follows under the studies ‘Setting up of a rare diseases biological samples bank (biobank) for research to facilitate pharmacological, gene and cell therapy trials in neuromuscular disorders (NMD)’ (REC reference number: 06/Q0406/33 - Hammersmith and Queen Charlotte’s and Chelsea Research Ethics Committee, ‘The use of cells as a model system to study pathogenesis and therapeutic strategies for Neuromuscular Disorders’ (REC reference 13/LO/1826 - London - Stanmore Research Ethics Committee) and ‘Genes and Proteins in Neuromuscular Disorders’ (REC reference 13/LO/1894 - London - Camberwell St Giles Research Ethics Committee).
Proteins were extracted from human and DMD dermal fibroblasts differentiated into myotubes after transfection by a lentiviral-mediated myoD construct carrying a puromycin selection cassette (for transduced cell enrichment) and a DsRed cassette (for assessing transduction efficiency). Protein lysates from control and DMD muscle were prepared from snap frozen tissues. Protein concentrations were measured using the Pierce BCA Protein assay kit (Thermo Scientific 23225), according to the manufacturer’s instructions.
Capillary Western immunoassay (Wes) analysis was performed on a Wes system (ProteinSimple) according to the manufacturer’s instructions using a 66–440 kDa Separation Module (ProteinSimple). For dystrophin a rabbit polyclonal anti-dystrophin antibody (ab15277, Abcam, dilution 1/50) and an anti-rabbit secondary antibody (042-206, Protein Simple) were used. The following amount of protein lysate (μg per well/capillary) were loaded for each sample: 1µg control fibroblasts MyoD-transfected; 1µg DMD fibroblasts MyoD-transfected; 0.125µg control muscle; 0.250μg DMD muscle.
Results
Patient characteristics
Participants lacking both the Dp427 and Dp140 isoforms, but not Dp71 (group 2) were on average diagnosed with DMD 6 months earlier than those who only lacked Dp427 (group 1, p=0.05, Table 1). There was no difference between groups 1 and 3 (p=0.99) and 2 and 3 (p=0.99).
Patient characteristics: DMD=Duchenne muscular dystrophy, sd=standard deviation, NSAA=North Star ambulatory assessment, G1=group 1, G2=group 2, G3=group 3, GC=glucocorticoid, na = not estimable.
We did not find a statistically significant difference in age of GC initiation (p=0.22), GC regime (p=0.09) or GC use recorded at any time between dystrophin isoform group, table 1. There was a non-significant trend for those in group 3 to be more likely to be on no GC (p=0.09, table 1).
Ambulatory function in DMD boys at 5 years of age
Mean NSAA score at 5 years of age was over 6 points lower in participants lacking all isoforms (Dp427/Dp140/Dp71; group 3) compared to those only lacking Dp427 (group 1; p<0.01) and 4.9 points lower in group 3 than those lacking both Dp427 and Dp140, but able to express Dp71 (group 2) (p=0.05, Table 1). There was no difference between groups 1 and 2 (p=0.99).
Mean rise from supine times at 5 years of age were 4.6 seconds slower in group 3 than group 1 (p<0.01) and 4.6 seconds slower in group 3 than group 2 (p=0.01, Table 1). There was no difference between groups 1 and 2 (p=0.99).
There was no difference in mean 10m walk/run time at 5 years of age between dystrophin isoform groups (p=0.99 for all group-wise comparisons). 10m walk/run time can be considered to be a less complex assessment than rise from supine time and the NSAA.
These results suggest that at 5 years of age, lack of Dp71 in addition to Dp427 and Dp140 was associated with poorer motor function in the two more complex motor assessments.
Peak motor function in DMD boys in the dystrophin isoform groups
Previous studies have reported a mean peak in NSAA score at 6-7 years of age[14,15]. In our NSAA score trajectory models, peak NSAA scores in the different isoform groups occurred within the age range of 5.5-8.0 years (Fig 1). Therefore, to assess if the isoform groups differed in maximal functional achievements, we looked at differences in peak mean NSAA score, 10m walk/run velocity and rise from supine time in the observed data for those aged 5.5-8.0 years in the dystrophin isoform groups (Table 2).
Mean peak NSAA scores for the dystrophin isoform groups overall and stratified by GC regime, mean peak 10m walk/run velocity and mean peak rise from supine time velocities not stratified by GC regimen. G1 = group 1, G2 = group 2, G3 = group 3 and GC=glucocorticoid. This table considers the dataset for a subset (n=263) of boys, aged 5.5 – 8.0 years, for whom majority GC regime data and peak NSAA scores were available for all 262 boys. This subset of 5.5-8.0 years of age was used as this is the age in which the majority of participants reach peak motor function[14,15]
Estimated mean NSAA score trajectory with age models in the dystrophin isoform groups. Each line represents estimated mean NSAA score plotted against age for the dystrophin isoform group. Group 1 = blue, group 2 = red, group 3 = green.
Mean peak NSAA scores were progressively lower for groups lacking more dystrophin isoforms; 4.0 points lower in group 3 than group 1 (p<0.01), 2.4 points lower in group 3 than group 2 (p=0.09) and 1.6 points lower in group 2 than group 1 (p=0.04; Table 2), after adjusting for GC regime. The age at attainment of peak score was not significantly different between the isoform groups (p=0.10).
Similarly, mean peak 10m walk/run velocity was 0.4 m/s (95%CI 0.2-0.7) slower in group 3 than group 1 (p<0.001) and 0.2 m/s (95%CI (0.05-0.4) slower in group 2 than group 1 (p<0.001), and non-significantly lower in group 3 compared to group 2 (p=0.10), after adjusting for GC regime. There was no difference in age of mean peak 10m walk/run velocity between isoform groups (p=0.91). Peak rise from supine time velocity was 0.04 m/s slower in group 3 than group 1 and 0.03 m/s slower in group 3 than group 2, however these differences were not significant (p=0.61).
Age of loss of ambulation (LOA) in DMD boys
The median ages of loss of ambulation were lower in groups 2 and 3 than group 1 (15.7 years in group 1, 13.0 years in group 2 and 13.6 years in group 3), however these differences were not significant (Table 1, p=0.67).
Effects of dystrophin isoform group on functional trajectories
NorthStar ambulatory assessment trajectory stratified by dystrophin isoform group
In NSAA score trajectory models, mean NSAA scores were lower in group 2 than group 1 (Fig 1). From 4.5-8.5 years of age, mean NSAA scores were lower in group 3 than in groups 1 and 2 with a cumulative effect of loss of isoforms (Fig 1). Peak estimated NSAA scores were lower in groups 2 and 3 with a cumulative effect of loss of isoforms (Fig 1).
10m walk/run velocity trajectory model stratified by dystrophin isoform group
In models of 10m walk/run velocity trajectory with age, velocities were higher in group 1 than groups 2 and 3 with a cumulative effect of loss of isoforms from approximately 6-8 years of age (Fig 2).
Estimated mean 10m walk/run velocity trajectory models in the dystrophin isoform groups. Each line represents estimated mean NSAA score trajectory plotted against age for the dystrophin isoform group. Group 1 = blue, group 2 = red, group 3 = green.
Rise from supine time velocity trajectories stratified by dystrophin isoform group
In mean rise from supine time velocity (m/s) trajectory models, group 3 had a slower estimated mean rise from supine time velocity than groups 1 and 2 from approximately 5-9 years of age (Fig 3).
Estimated mean rise from supine time velocity trajectory models in the dystrophin isoform groups. Each line represents estimated mean NSAA score trajectory plotted against age for the dystrophin isoform group. Group 1 = blue, group 2 = red, group 3 = green.
Influence of cognitive impairment on functional outcomes
We evaluated the influence of cognitive impairment on peak NSAA scores in a subset of 127 boys for whom cognition group, GC regimen and peak NSAA scores were available.
We first confirmed a strong cumulative association between lack of brain dystrophin isoforms and cognitive impairment (p<0.001, Supplementary Table 1). Cognitive impairment was found in 32% (20/63) of those in group 1, 53% (29/55) of those in group 2 and 78% (7/9) of those in group 3 (Supplementary Table 1), in keeping with previously published studies[10–13].
Mean peak NSAA scores were 2.2 points lower in those with cognitive impairment than those with normal cognition (p=0.04), after adjusting for GC regimen (Supplementary table 1).
In NSAA trajectory models stratified by cognitive group, mean NSAA scores were lower in those with cognitive impairment than those without cognitive impairment at all ages (Fig 4)
Estimated mean NSAA score trajectory with age models in those with and without cognitive impairment. Each line represents estimated mean NSAA score trajectory plotted against age for the cognition group. Black = normal cognition and red = impaired cognition.
Forelimb grip strength measurements in mdx, mdx52 and DMD-null mice at the age of 3 months
Average grip strength in peak force at 3 months of age was significantly higher for mdx mice than that of mdx52 mice (p = 0.01), with statistical differences between WT and mdx52 mice, and between WT and DMD-null mice (Fig. 5a). There were no statistically significant differences in the degree of grip fatigue at the age of 3 months between the three DMD mouse strains.
Repeat testing of grip strength in WT (n=4), mdx (n=9), mdx52 (n=10) and DMD-null (n=3) mice at the age of 3 months. Mean grip strength at the age of 3 months; n = 3-10, Mean ± SEM, One-way ANOVA with Dunnett p = 0.0005 (WT vs mdx52), p = 0.0052 (WT vs DMD-null), p = 0.01 (mdx vs mdx52). (5a). The degree of grip fatigue at the age of 3 months. Mean ± SEM, One-way ANOVA with Dunnett (5b). There was no statistically significant difference in grip fatigue at the age of 3 months between the 3 DMD mouse strains.
Dystrophin isoform protein production in wild type and DMD skeletal muscle and myogenic cells
A band corresponding to Dp71 was seen in control and DMD MyoD transfected fibroblasts and in DMD skeletal muscle, but not in control skeletal muscle (Fig 6). The band corresponding to full length dystrophin was only seen in MyoD transfected control fibroblasts and control skeletal muscle, but not DMD MyoD transfected fibroblasts or skeletal muscle (Fig 6).
Virtual blot of dystrophin signal detected by Wes using a C-terminus polyclonal anti dystrophin antibody (Abcam 15277) in: 1) control fibroblasts MyoD transfected; 2) DMD fibroblasts MyoD transfected; 3) control muscle; 4) DMD muscle. The amount of proteins used for each sample were: 1) 1µg; 2) 1µg; 3) 0.2501µg; 4) 0.1251µg. As previously reported, the full length dystrophin signal detected by Wes displays a molecular weight of approximately 300 kDa, which is lower than the predicted 427 kDa for the full-length (muscle) dystrophin isoform[25].
Discussion
DMD is characterised by significant variability in clinical progression[15,26]. Latent class trajectory analysis identified different classes of progression, but did not necessarily provide information on the drivers of the observed variability[15]. The reduced ability to identify drivers of this variability complicates assessments of DMD boys in clinical practice and clinical trials.
We hypothesised that DMD mutations expected to differentially impact on expression of different dystrophin isoforms could differentially affect motor function in DMD in both humans and mouse models. We focused on dystrophin isoforms previously demonstrated to play an important role in brain function. We performed a series of investigations characterising dystrophin isoform production in wild-type and DMD myogenic and skeletal muscle and characterising patterns of motor impairment in 3 DMD mouse models, wild-type mice and a large cohort of DMD boys with differential dystrophin isoform involvement.
Our results clearly demonstrate marked reductions in motor function achievement at 5 years of age (reduced mean NSAA scores and rise from supine times) and peak motor function (reduced mean peak NSAA scores and 10m walk/run velocities) in DMD boys lacking Dp140 and Dp71, with a clear cumulative effect of loss of isoforms. The differences in NSAA score observed between the isoform groups was considerable, substantially exceeding the minimally clinically important difference in NSAA score of approximately 3 points, highlighting the clear clinical significance of our findings[27].
We assessed the impact of loss of different dystrophin isoforms on grip strength in mdx (Dp427 absent, Dp140/Dp71 present); mdx52 (Dp427/Dp140 absent, Dp71 present); and DMD-null mice (Dp427/Dp140/Dp71 absent). As expected, mdx mice had significantly reduced grip strength compared to wild type mice. When assessing the impact of isoform loss, average grip strength in mdx mice at 3 months of age was significantly higher than in mdx52 mice (p = 0.01), but no further decline was observed in the DMD-null mice lacking all 3 isoforms.
When we studied dystrophin isoform protein production, we could not detect Dp140 in normal or DMD skeletal muscle or myogenic cultures. With the highly sensitive Wes technique, we detected low levels of Dp71 production in skeletal muscle of DMD boys, but not controls. This contrasted with the high Dp71 levels detected in both normal and DMD myogenic cultures.
This study directly correlates motor impairment in DMD humans and DMD mouse models to different patterns of dystrophin isoform expression.
A previous clinical study found DMD children with intellectual disability often have worse motor outcomes[28]. However, when genotype grouping was done, patients were subdivided between those with mutations proximal or distal to exon 30. Hence not allowing to distinguish the role of Dp140 and Dp71 in motor outcomes[28].
While our data unequivocally demonstrate differences in motor performance between these isoform groups, they do not necessarily provide a definitive mechanistic explanation for this. We considered a number of potential factors: differences in glucocorticoid use; a direct role of Dp140 and Dp71 in skeletal muscle function, differences in comprehension of motor outcome instructions; or impact of lack of Dp140 and Dp71 on higher aspects of motor coordination and planning.
With respect to glucocorticoid use, we found no difference in glucocorticoid regime between isoform groups and the pattern of cumulative marked reductions in motor outcomes in DMD boys lacking Dp140 and Dp71 remained present after adjusting for glucocorticoid regime.
Regarding the possibility that DMD boys lacking Dp140 and Dp71 were less able to understand instructions for performing assessments, we consider this unlikely as the components of the assessments are part of the children’s daily routine, the assessments are carried out every 6 months in hospital and the specialist physiotherapists have a standard protocol and regular training in carrying out the assessments.
With regards to the mdx mice study, the mdx-52 mice had reduced grip strength compared to the mdx mice. However, we did not detect a significant difference in grip strength between the DMD-null and mdx52 mice.
Our studies in DMD boys and dystrophic mouse models provide clear evidence of a relationship between patterns of DMD isoform involvement and motor performance and open the field to various possible explanations to account for these findings. The fact that the pattern of Dp140 expression is confined to the CNS makes the possibility of a skeletal muscle role for Dp140 not plausible. However, a possible, and previously not appreciated, role of Dp140 on motor neuron function cannot be excluded and could have contributed to the observed differences between Dp427 only and Dp427 and Dp140 deficient boys and mice.
Regarding Dp71, we confirm previous findings of Dp71 production in DMD muscle, but not control muscle[29]. In DMD muscle, Dp71 is not localised at the sarcolemma, C-terminal dystrophin staining is the standard diagnostic test and DMD patients universally demonstrate absent C-terminal dystrophin at the sarcolemma[30]. This observation, together with the absence of an actin binding domain in Dp71, argues against a peripheral effect of this isoform in functional motor outcomes, which is also confirmed by the dystrophic mouse studies. We confirm that myotubes in culture produce high levels of Dp71, and that this capacity is retained in DMD patients with mutations upstream of the Dp71 promoter. This, and previous findings, strongly, indicate a role for Dp71 in myoblast cell proliferation and satellite cell activation[31–33]. This suggests that the low level of Dp71 production in DMD muscle is a secondary phenomenon related to the muscle degeneration and regeneration characteristic of DMD.
One hypothesis we consider plausible is that the observed differences in motor function could be related to a direct impact of CNS involvement on motor function[11–13,28]. Whilst the concept of an impact of brain function on motor performance is novel in the DMD field, this is a well-recognised concept in other conditions. Westendorp et al. compared gross motor skills of children with learning disabilities (LD) with age-matched typically developing children[34]. Children with LD had lower locomotor scores than typically developing children[34]. Several studies have found that children with ASD without skeletal muscle involvement have both gross and fine motor delay[35,36]. These observations are relevant as the prevalence of LD is higher in DMD boys, particularly in those lacking Dp140 and Dp71[10].
Previous studies have hypothesised a possible role of the cerebellum and/or the cerebellar thalamic cortical connectivity in DMD cognitive difficulties[37]. Battini et al found a deficits of executive function, specifically in tasks requiring planning and directing goal oriented behaviour, in DMD boys[37]. Hellebrekers et al found Dp140 negative males had slower information processing speeds than Dp140 positive males[38].
With reference to the mdx mice, Chaussenot et al demonstrated executive dysfunction in a Dp71 deficient mouse model, whilst no deep CNS phenotyping is available for mdx52 mice[39].
Dp427m, Dp140 and Dp71 are all expressed in the adult human cerebellum[7]. The cerebellum plays a crucial role in the control of goal-directed movements and timing of coordinated movement[40]. Features of cerebellar movement dysfunction in children include gross motor delay and poor coordination with complex movements (with relative sparing of simple motor tasks)[41]. Interestingly, of the 3 main ambulatory motor outcomes we considered, the NSAA score was the only outcome showing impairment in those lacking the shorter isoforms at both 5 years and peak motor function. Of these 3 outcomes, the NSAA is the most complex, requiring the highest level of motor coordination and planning.
Taken together, we consider a plausible contributor for the patterns of more severe motor impairment in boys and mdx mice lacking Dp140 and Dp140/Dp71 to be the associated deficits in the higher centres of motor function, control and coordination, with a possible role of deficits in the cerebellum and/or the cerebellar thalamic cortical connectivity.
The strengths of this study include the large number of patients, data collected in a ‘real world’ clinical setting and the longitudinal nature of the data. Limitations included some missing data due to the real-world data collection in routine clinical appointments. However, all missing data was accounted for statistically.
In summary, we found that mutations resulting in lack of Dp140 or Dp140/Dp71 in addition to Dp427 are associated with poorer motor function performance in DMD boys and dystrophic mouse models, with a cumulative effect of loss of isoforms in DMD boys. Given the major role of these isoforms in brain function, we hypothesise this could be at least partly related to deficits in the higher centres of motor function, control and coordination in those lacking Dp140 and Dp140/Dp71.
Irrespective of the possible pathophysiological explanation, our novel findings provide evidence for a relationship between the site of the DMD mutations and the effect on shorter dystrophin isoform production on not only cognitive, but also motor outcomes. This has crucial implications for clinical practice and clinical trial design involving patients with DMD.
Data Availability
The data that support the findings of this study are available from the corresponding author, [FM], upon reasonable request.
Author’s financial relationships and potential conflicts of interest
Mary Chesshyre, Deborah Ridout, Kate Maresh, Lianne Abbott, Vandana Ayyar Gupta, Marion Main, Yoko Ookubo, Yasumasa Hashimoto and Silvia Torelli report no financial disclosures or potential conflicts of interest.
Yoshitsugu Aoki is the grant holder of the grant Grants in Aid for Research on Nervous and Mental Disorders (grant number 2-6)
Kate Maresh has received payments to her institution from Great Ormond Street Hospital Children’s Charity and MRC Centre for Neuromuscular Diseases, Queen Square, London in the last 36 months.
Valeria Ricotti – Co-founder, EVP and CMO of DiNAQOR, shareholder of Solid Biosciences
Adnan Manzur is the clinical lead for the North Star clinical network and is one of the co-grant holders from the MDUK for maintenance of the network.
Mariacristina Scoto has received speaker and consultancy honoraria For Roche, Avexis, Santhera and Biogen
Giovanni Baranello has received consultancy honoraria from AveXis, Roche, Biogen, PTC, and Sarepta Therapeutics. Dr Giovanni Baranello has received speaker honoraria from AveXis, Roche and PTC
Francesco Muntoni is supported by the NIHR Great Ormond Street Hospital Biomedical Research Centre and has received speaker and consultancy honoraria from Sarepta Therapeutics, Avexis, PTC Therapeutics, Roche and Pfizer.
Ethical guidelines and consent
The North Star clinical network study was reviewed by the North Sheffield Research Ethics Committee and it was felt that it didn’t need to go to a research ethics committee. Written informed consent was obtained for the collection of all clinical data and the NorthStar clinical network project has Caldicott Guardian approval. All clinical assessments are conducted according to the principles of the Declaration of Helsinki (2000) and its later amendments and the Principles of Good Clinical Practice.
Tissue sampling was approved by the NHS national Research Ethics Service as follows under the studies ‘Setting up of a rare diseases biological samples bank (biobank) for research to facilitate pharmacological, gene and cell therapy trials in neuromuscular disorders (NMD)’ (REC reference number: 06/Q0406/33 - Hammersmith and Queen Charlotte’s and Chelsea Research Ethics Committee, ‘The use of cells as a model system to study pathogenesis and therapeutic strategies for Neuromuscular Disorders’ (REC reference 13/LO/1826 - London - Stanmore Research Ethics Committee) and ‘Genes and Proteins in Neuromuscular Disorders’ (REC reference 13/LO/1894 - London - Camberwell St Giles Research Ethics Committee).
Animal studies were performed in accordance with the appropriate ethics committee and have therefore been performed with the ethical standards laid down in the 1964 Declaration of Helsinki and its later amendments.
Acknowledgements
We are grateful to the DMD patients and their families, the North Star clinical network study group and its senior clinico-academic coordinator Dr Vandana Ayyar Gupta, which is co-led by Francesco Muntoni (f.muntoni{at}ucl.ac.uk) and Adnan Y Manzur (adnan.manzur{at}gosh.nhs.uk), Muscular Dystrophy UK (MDUK) for funding the North Star network and Certus Technology Associates Limited for hosting the database. Note that the lead authors for the North Star Group are Dr Francesco Muntoni, f.muntoni{at}ucl.ac.uk; and Dr Adnan Manzur, a.manzur{at}gosh.nhs.uk. The North Star DMD Network is supported by a grant from Muscular Dystrophy UK to Dr Adnan Manzur and Prof Francesco Muntoni, at UCL. We are very grateful to Cyrille Vaillend (Neuroscience Paris-Saclay Institute (Neuro-PSI), UMR 9197, Université Paris Sud, CNRS, Université Paris Saclay, 91190 Orsay, France) for his helpful review of and feedback on the manuscript.
The support of L’Association Française contre les Myopathies for the ‘Outcome measures in Duchenne Muscular Dystrophy: A Natural History Study: Imaging and Neuropsychology Sub-study’ is gratefully acknowledged. The support of Great Ormond Street Hospital Children’s Charity for ‘A Study of Emotional Function in Duchenne Muscular Dystrophy (EmoDe Study)’ is gratefully acknowledged. Kate Maresh was supported by the Medical Research Council (MRC) via the MRC Centre for Neuromuscular Diseases, Institute of Neurology, University College London.
The support of Great Ormond Street Hospital Children’s Charity and the NIHR GOSH biomedical research centre (BRC) for the ‘Neurodevelopmental, emotional, and behavioural problems in Duchenne muscular dystrophy in relation to underlying dystrophin gene mutations’ study is gratefully acknowledged19.
This work is supported by the NIHR GOSH BRC. The views expressed are participants of the author(s) and not necessarily participants of the NHS, the NIHR or the Department of Health.
The support of EUH2020 grant 83245 Brain Involvement iN Dystrophinopathies is also gratefully acknowledged.
The support of the grant Grants in Aid for Research on Nervous and Mental Disorders (grant number 2-6 to Yoshitsugu Aoki) for the functional testing using dystrophic and wild-type mice is gratefully acknowledged.
References
