Drug trials in children engage with many ethical issues, from drug-related safety concerns to communication with patients and parents, and recruitment and informed consent procedures. This paper addresses the field of neuromuscular disorders where the possibility of genetic, mutation-specific treatments, has added new complexity. Not only must trial design address issues of equity of access, but researchers must also think through the implications of adopting a personalised medicine approach, which requires a precise molecular diagnosis, in addition to other implications of developing orphan drugs. It is against this background of change and complexity that the Project Ethics Council (PEC) was established within the TREAT-NMD EU Network of Excellence. The PEC is a high level advisory group that draws upon the expertise of its interdisciplinary membership which includes clinicians, lawyers, scientists, parents, representatives of patient organisations, social scientists and ethicists. In this paper we describe the establishment and terms of reference of the PEC, give an indication of the range and depth of its work and provide some analysis of the kinds of complex questions encountered. The paper describes how the PEC has responded to substantive ethical issues raised within the TREAT-NMD consortium and how it has provided a wider resource for any concerned parent, patient, or clinician to ask a question of ethical concern. Issues raised range from science related ethical issues, issues related to hereditary neuromuscular diseases and the new therapeutic approaches and questions concerning patients rights in the context of patient registries and bio-banks. We conclude by recommending the PEC as a model for similar research contexts in rare diseases.

Duchenne muscular dystrophy is caused by mutations in the DMD gene that disrupt the open reading frame and prevent the full translation of its protein product, dystrophin. Restoration of the open reading frame and dystrophin production can be achieved by exon skipping using antisense oligonucleotides targeted to splicing elements. This approach aims to transform the Duchenne muscular dystrophy phenotype to that of the milder disorder, Becker muscular dystrophy, typically caused by in-frame dystrophin deletions that allow the production of an internally deleted but partially functional dystrophin. There is ongoing debate regarding the functional properties of the different internally deleted dystrophins produced by exon skipping for different mutations; more insight would be valuable to improve and better predict the outcome of exon skipping clinical trials. To this end, we have characterized the clinical phenotype of 17 patients with Becker muscular dystrophy harbouring in-frame deletions relevant to on-going or planned exon skipping clinical trials for Duchenne muscular dystrophy and correlated it to the levels of dystrophin, and dystrophin-associated protein expression. The cohort of 17 patients, selected exclusively on the basis of their genotype, included 4 asymptomatic, 12 mild and 1 severe patient. All patients had dystrophin levels of >40% of control and significantly higher dystrophin (P = 0.013), β-dystroglycan (P = 0.025) and neuronal nitric oxide synthase (P = 0.034) expression was observed in asymptomatic individuals versus symptomatic patients with Becker muscular dystrophy. Furthermore, grouping the patients by deletion, patients with Becker muscular dystrophy with deletions with an end-point of exon 51 (the skipping of which could rescue the largest group of Duchenne muscular dystrophy deletions) showed significantly higher dystrophin levels (P = 0.034) than those with deletions ending with exon 53. This is the first quantitative study on both dystrophin and dystrophin-associated protein expression in patients with Becker muscular dystrophy with deletions relevant for on-going exon skipping trials in Duchenne muscular dystrophy. Taken together, our results indicate that all varieties of internally deleted dystrophin assessed in this study have the functional capability to provide a substantial clinical benefit to patients with Duchenne muscular dystrophy.

In Duchenne muscular dystrophy, the exon-skipping approach has obtained proof of concept in animal models, myogenic cell cultures, and following local and systemic administration in Duchenne patients. Indeed, we have previously demonstrated that low doses (7.5 mg/Kg/week) of 2′-O-methyl-phosphorothioate antisense oligoribonucleotides (AONs) adsorbed onto ZM2 nanoparticles provoke widespread dystrophin restoration 7 days after intraperitoneal treatment in mdx mice. In this study, we went on to test whether this dystrophin restoration was still measurable 90 days from the end of the same treatment. Interestingly, we found that both western blot and immunohistochemical analysis (up to 7% positive fibres) were still able to detect dystrophin protein in the skeletal muscles of ZM2-AON-treated mice at this time, and the level of exon-23 skipping could still be assessed by RT real-time PCR (up to 10% of skipping percentage). In contrast, the protein was undetectable by western blot analysis in the skeletal muscles of mdx mice treated with an identical dose of naked AON, and the percentage of dystrophin-positive fibres and exon-23 skipping were reminiscent of those of untreated mdx mice. Our data therefore demonstrate the long-term residual efficacy of this systemic low-dose treatment and confirm the protective effect nanoparticles exert on AON molecules.

The 2.2 Mb long dystrophin (DMD) gene, the largest gene in the human genome, corresponds to roughly 0.1% of the entire human DNA sequence. Mutations in this gene cause Duchenne muscular dystrophy and other milder X-linked, recessive dystrophinopathies. Using a custom-made tiling array, specifically designed for the DMD locus, we identified a variety of novel long non-coding RNAs (lncRNAs), both sense and antisense oriented, whose expression profiles mirror that of DMD gene. Importantly, these transcripts are intronic in origin and specifically localized to the nucleus and are transcribed contextually with dystrophin isoforms or primed by MyoD-induced myogenic differentiation. Furthermore, their forced ectopic expression in both human muscle and neuronal cells causes a specific and negative regulation of endogenous dystrophin full length isoforms and significantly down-regulate the activity of a luciferase reporter construct carrying the minimal promoter regions of the muscle dystrophin isoform. Consistent with this apparently repressive role, we found that, in muscle samples of dystrophinopathic female carriers, lncRNAs expression levels inversely correlate with those of muscle full length DMD isoforms. Overall these findings unveil an unprecedented complexity of the transcriptional pattern of the DMD locus and reveal that DMD lncRNAs may contribute to the orchestration and homeostasis of the muscle dystrophin expression pattern by either selective targeting and down-modulating the dystrophin promoter transcriptional activity.

Background

Although Duchenne and Becker muscular dystrophies, X-linked recessive myopathies, predominantly affect males, a clinically significant proportion of females manifesting symptoms have also been reported. They represent an heterogeneous group characterized by variable degrees of muscle weakness and/or cardiac involvement. Though preferential inactivation of the normal X chromosome has long been considered the principal mechanism behind disease manifestation in these females, supporting evidence is controversial.

Methods

Eighteen females showing a mosaic pattern of dystrophin expression on muscle biopsy were recruited and classified as symptomatic (7) or asymptomatic (11), based on the presence or absence of muscle weakness. The causative DMD gene mutations were identified in all cases, and the X-inactivation pattern was assessed in muscle DNA. Transcriptional analysis in muscles was performed in all females, and relative quantification of wild-type and mutated transcripts was also performed in 9 carriers. Dystrophin protein was quantified by immunoblotting in 2 females.

Results

The study highlighted a lack of relationship between dystrophic phenotype and X-inactivation pattern in females; skewed X-inactivation was found in 2 out of 6 symptomatic carriers and in 5 out of 11 asymptomatic carriers. All females were characterized by biallelic transcription, but no association was found between X-inactivation pattern and allele transcriptional balancing. Either a prevalence of wild-type transcript or equal proportions of wild-type and mutated RNAs was observed in both symptomatic and asymptomatic females. Moreover, very similar levels of total and wild-type transcripts were identified in the two groups of carriers.

Conclusions

This is the first study deeply exploring the DMD transcriptional behaviour in a cohort of female carriers. Notably, no relationship between X-inactivation pattern and transcriptional behaviour of DMD gene was observed, suggesting that the two mechanisms are regulated independently. Moreover, neither the total DMD transcript level, nor the relative proportion of the wild-type transcript do correlate with the symptomatic phenotype.

Collagen VI is a major extracellular matrix (ECM) protein with a critical role in maintaining skeletal muscle functional integrity. Mutations in COL6A1, COL6A2 and COL6A3 genes cause Ullrich Congenital Muscular Dystrophy (UCMD), Bethlem Myopathy, and Myosclerosis. Moreover, Col6a1−/− mice and collagen VI deficient zebrafish display a myopathic phenotype. Recently, two additional collagen VI chains were identified in humans, the α5 and α6 chains, however their distribution patterns and functions in human skeletal muscle have not been thoroughly investigated yet. By means of immunofluorescence analysis, the α6 chain was detected in the endomysium and perimysium, while the α5 chain labeling was restricted to the myotendinous junctions. In normal muscle cultures, the α6 chain was present in traces in the ECM, while the α5 chain was not detected. In the absence of ascorbic acid, the α6 chain was mainly accumulated into the cytoplasm of a sub-set of desmin negative cells, likely of interstitial origin, which can be considered myofibroblasts as they expressed α-smooth muscle actin. TGF-β1 treatment, a pro-fibrotic factor which induces trans-differentiation of fibroblasts into myofibroblasts, increased the α6 chain deposition in the extracellular matrix after addition of ascorbic acid. In order to define the involvement of the α6 chain in muscle fibrosis we studied biopsies of patients affected by Duchenne Muscular Dystrophy (DMD). We found that the α6 chain was dramatically up-regulated in fibrotic areas where, in contrast, the α5 chain was undetectable. Our results show a restricted and differential distribution of the novel α6 and α5 chains in skeletal muscle when compared to the widely distributed, homologous α3 chain, suggesting that these new chains may play specific roles in specialized ECM structures. While the α5 chain may have a specialized function in tissue areas subjected to tensile stress, the α6 chain appears implicated in ECM remodeling during muscle fibrosis.

Highlights

► We study collagen VI alpha 5 and alpha 6 chains expression in human skeletal muscle. ► ► We show that alpha 5 chain is exclusively present at the myotendinous junctions. ► We show that alpha 6 chain is present in the ECM but not at the basement membranes. ► We show that alpha 6 chain increases in Duchenne muscular dystrophy muscle fibrosis. ► We hypothesize the possible formation of alpha1-2-6 or alpha1-2-5 chains trimers.

Background

In X-linked dilated cardiomyopathy due to dystrophin mutations which abolish the expression of the M isoform (5'-XLDC), the skeletal muscle is spared through the up-regulation of the Brain (B) isoform, a compensatory mechanism that does not appear to occur in the heart of affected individuals.

Methods

We quantitatively studied the expression topography of both B and M isoforms in various human heart regions through in-situ RNA hybridization, Reverse-Transcriptase and Real-Time PCR experiments. We also investigated the methylation profile of the B promoter region in the heart and quantified the B isoform up regulation in the skeletal muscle of two 5'-XLDC patients.

Results

Unlike the M isoform, consistently detectable in all the heart regions, the B isoform was selectively expressed in atrial cardiomyocytes, but absent in ventricles and in conduction system structures. Although the level of B isoform messenger in the skeletal muscle of 5'-XLDC patients was lower that of the M messenger present in control muscle, it seems sufficient to avoid an overt muscle pathology. This result is consistent with the protein level in XLDC patients muscles we previously quantified. Methylation studies revealed that the B promoter shows an overall low level of methylation at the CG dinucleotides in both atria and ventricles, suggesting a methylation-independent regulation of the B promoter activity.

Conclusions

The ventricular dilatation seen in 5'-XLDC patients appears to be functionally related to loss of the M isoform, the only isoform transcribed in human ventricles; in contrast, the B isoform is well expressed in heart but confined to the atria. Since the B isoform can functionally replace the M isoform in the skeletal muscle, its expression in the heart could potentially exert the same rescue function. Methylation status does not seem to play a role in the differential B promoter activity in atria and ventricles, which may be governed by other regulatory mechanisms. If these mechanisms could be deduced, de-silencing of the B isoform may represent a therapeutic strategy in 5'-XLDC patients.

Elucidation of new biomarkers and potential drug targets from high-throughput profiling data is a challenging task due to a limited number of available biological samples and questionable reproducibility of differential changes in cross-dataset comparisons. In this paper we propose a novel computational approach for drug and biomarkers discovery using comprehensive analysis of multiple expression profiling datasets.

The new method relies on aggregation of individual profiling experiments combined with leave-one-dataset-out validation approach. Aggregated datasets were studied using Sub-Network Enrichment Analysis algorithm (SNEA) to find consistent statistically significant key regulators within the global literature-extracted expression regulation network. These regulators were linked to the consistent differentially expressed genes.

We have applied our approach to several publicly available human muscle gene expression profiling datasets related to Duchenne muscular dystrophy (DMD). In order to detect both enhanced and repressed processes we considered up- and down-regulated genes separately. Applying the proposed approach to the regulators search we discovered the disturbance in the activity of several muscle-related transcription factors (e.g. MYOG and MYOD1), regulators of inflammation, regeneration, and fibrosis. Almost all SNEA-derived regulators of down-regulated genes (e.g. AMPK, TORC2, PPARGC1A) correspond to a single common pathway important for fast-to-slow twitch fiber type transition. We hypothesize that this process can affect the severity of DMD symptoms, making corresponding regulators and downstream genes valuable candidates for being potential drug targets and exploratory biomarkers.

Author Summary

Comparison of gene expression in diseased and normal tissue is a powerful tool of studying processes involved in pathogenesis and searching for potential drug targets and biomarkers of the disease's progression and treatment outcome. We have developed a novel approach for systematic knowledge-driven analysis of gene expression profiling data, which can suggest the underlying cause of the observed differential expression by identifying which expression regulators might be involved. These regulators can not only be the promising subjects of further investigation, but also potential drug targets, as normalization of their activity might alleviate some of the disease's symptoms. The targets downstream of suggested regulators can be proposed as exploratory biomarkers in disease treatment and prognosis. We used our approach to analyze public gene expression datasets of Duchenne muscular dystrophy – a progressive inherited disease in males. Some of the regulators and biomarkers that we found were already investigated in the context of DMD, while some of them were not yet studied and may be of interest for biological and clinical studies.

Six individuals with Ullrich congenital muscular dystrophy (UCMD) and mutations in the genes-encoding collagen VI, aging 5–9, received 3–5 mg/kg of cyclosporine A (CsA) daily for 1 to 3.2 years. The primary outcome measure was the muscle strength evaluated with a myometer and expressed as megalimbs. The megalimbs score showed significant improvement (P = 0.01) in 5 of the 6 patients. Motor function did not change. Respiratory function deteriorated in all. CsA treatment corrected mitochondrial dysfunction, increased muscle regeneration, and decreased the number of apoptotic nuclei. Results from this study demonstrate that long-term treatment with CsA ameliorates performance in the limbs, but not in the respiratory muscles of UCMD patients, and that it is well tolerated. These results suggest considering a trial of CsA or nonimmunosuppressive cyclosporins, that retains the PTP-desensitizing properties of CsA, as early as possible in UCMD patients when diaphragm is less compromised.

Antisense oligonucleotides (AONs) are being developed as RNA therapeutic molecules for Duchenne muscular dystrophy. For oligonucleotides with the 2′-O-methyl-phosphorothioate (2OMePS) RNA chemistry, proof of concept has been obtained in patient-specific muscle cell cultures, the mouse and dog disease models, and recently by local administration in Duchenne patients. To further explore the pharmacokinetic (PK)/pharmacodynamic (PD) properties of this chemical class of oligonucleotides, we performed a series of preclinical studies in mice. The results demonstrate that the levels of oligonucleotides in dystrophin-deficient muscle fibers are much higher than in healthy fibers, leading to higher exon-skipping levels. Oligonucleotide levels and half-life differed for specific muscle groups, with heart muscle showing the lowest levels but longest half-life (~46 days). Intravenous (i.v.), subcutaneous (s.c.), and intraperitoneal (i.p.) delivery methods were directly compared. For each method, exon-skipping and novel dystrophin expression were observed in all muscles, including arrector pili smooth muscle in skin biopsies. After i.v. administration, the oligonucleotide peak levels in plasma, liver, and kidney were higher than after s.c. or i.p. injections. However, as the bioavailability was similar, and the levels of oligonucleotide, exon-skipping, and dystrophin steadily accumulated overtime after s.c. administration, we selected this patient-convenient delivery method for future clinical study protocols.

Background

Multiple sclerosis (MS) is a complex disorder thought to result from an interaction between environmental and genetic predisposing factors which have not yet been characterised, although it is known to be associated with the HLA region on 6p21.32. Recently, a picture of chronic cerebrospinal venous insufficiency (CCSVI), consequent to stenosing venous malformation of the main extra-cranial outflow routes (VM), has been described in patients affected with MS, introducing an additional phenotype with possible pathogenic significance.

Methods

In order to explore the presence of copy number variations (CNVs) within the HLA locus, a custom CGH array was designed to cover 7 Mb of the HLA locus region (6,899,999 bp; chr6:29,900,001-36,800,000). Genomic DNA of the 15 patients with CCSVI/VM and MS was hybridised in duplicate.

Results

In total, 322 CNVs, of which 225 were extragenic and 97 intragenic, were identified in 15 patients. 234 known polymorphic CNVs were detected, the majority of these being situated in non-coding or extragenic regions. The overall number of CNVs (both extra- and intragenic) showed a robust and significant correlation with the number of stenosing VMs (Spearman: r = 0.6590, p = 0.0104; linear regression analysis r = 0.6577, p = 0.0106).

The region we analysed contains 211 known genes. By using pathway analysis focused on angiogenesis and venous development, MS, and immunity, we tentatively highlight several genes as possible susceptibility factor candidates involved in this peculiar phenotype.

Conclusions

The CNVs contained in the HLA locus region in patients with the novel phenotype of CCSVI/VM and MS were mapped in detail, demonstrating a significant correlation between the number of known CNVs found in the HLA region and the number of CCSVI-VMs identified in patients. Pathway analysis revealed common routes of interaction of several of the genes involved in angiogenesis and immunity contained within this region. Despite the small sample size in this pilot study, it does suggest that the number of multiple polymorphic CNVs in the HLA locus deserves further study, owing to their possible involvement in susceptibility to this novel MS/VM plus phenotype, and perhaps even other types of the disease.

Background

Molecular characterization of collagen-VI related myopathies currently relies on standard sequencing, which yields a detection rate approximating 75-79% in Ullrich congenital muscular dystrophy (UCMD) and 60-65% in Bethlem myopathy (BM) patients as PCR-based techniques tend to miss gross genomic rearrangements as well as copy number variations (CNVs) in both the coding sequence and intronic regions.

Methods

We have designed a custom oligonucleotide CGH array in order to investigate the presence of CNVs in the coding and non-coding regions of COL6A1, A2, A3, A5 and A6 genes and a group of genes functionally related to collagen VI. A cohort of 12 patients with UCMD/BM negative at sequencing analysis and 2 subjects carrying a single COL6 mutation whose clinical phenotype was not explicable by inheritance were selected and the occurrence of allelic and genetic heterogeneity explored.

Results

A deletion within intron 1A of the COL6A2 gene, occurring in compound heterozygosity with a small deletion in exon 28, previously detected by routine sequencing, was identified in a BM patient. RNA studies showed monoallelic transcription of the COL6A2 gene, thus elucidating the functional effect of the intronic deletion. No pathogenic mutations were identified in the remaining analyzed patients, either within COL6A genes, or in genes functionally related to collagen VI.

Conclusions

Our custom CGH array may represent a useful complementary diagnostic tool, especially in recessive forms of the disease, when only one mutant allele is detected by standard sequencing. The intronic deletion we identified represents the first example of a pure intronic mutation in COL6A genes.

Background

The commonest pathogenic DMD changes are intragenic deletions/duplications which make up to 78% of all cases and point mutations (roughly 20%) detectable through direct sequencing. The remaining mutations (about 2%) are thought to be pure intronic rearrangements/mutations or 5'-3' UTR changes. In order to screen the huge DMD gene for all types of copy number variation mutations we designed a novel custom high density comparative genomic hybridisation array which contains the full genomic region of the DMD gene and spans from 100 kb upstream to 100 kb downstream of the 2.2 Mb DMD gene.

Results

We studied 12 DMD/BMD patients who either had no detectable mutations or carried previously identified quantitative pathogenic changes in the DMD gene. We validated the array on patients with previously known mutations as well as unaffected controls, we identified three novel pure intronic rearrangements and we defined all the mutation breakpoints both in the introns and in the 3' UTR region. We also detected a novel polymorphic intron 2 deletion/duplication variation. Despite the high resolution of this approach, RNA studies were required to confirm the functional significance of the intronic mutations identified by CGH. In addition, RNA analysis identified three intronic pathogenic variations affecting splicing which had not been detected by the CGH analysis.

Conclusion

This novel technology represents an effective high throughput tool to identify both common and rarer DMD rearrangements. RNA studies are required in order to validate the significance of the CGH array findings. The combination of these tools will fully cover the identification of causative DMD rearrangements in both coding and non-coding regions, particularly in patients in whom standard although extensive techniques are unable to detect a mutation.

Dystrophin is a subsarcolemmal protein that, by linking the actin cytoskeleton to the extracellular matrix via dystroglycans, is critical for the integrity of muscle fibers. Here, we report that epidermal melanocytes, obtained from conventional skin biopsy, express dystrophin with a restricted localization to the plasma membrane facing the dermal–epidermal junction. In addition the full-length muscle isoform mDp427 was clearly detectable in melanocyte cultures as assessed by immunohistochemistry, RNA, and Western blot analysis. Melanocytes of Duchenne muscular dystrophy (DMD) patients did not express dystrophin, and the ultrastructural analysis revealed typical mitochondrial alterations similar to those occurring in myoblasts from the same patients. Mitochondria of melanocytes from DMD patients readily accumulated tetramethylrhodamine methyl ester, indicating that they are energized irrespective of the presence of dystrophin but, at variance from mitochondria of control donors, depolarized upon the addition of oligomycin, suggesting that they are affected by a latent dysfunction unmasked by inhibition of the ATP synthase. Pure melanocyte cultures can be readily obtained by conventional skin biopsies and may be a feasible and reliable tool alternative to muscle biopsy for functional studies in dystrophinopathies. The mitochondrial dysfunction occurring in DMD melanocytes could represent a promising cellular biomarker for monitoring dystrophinopathies also in response to pharmacological treatments. J. Cell. Physiol. 228: 1323–1331, 2013. © 2012 Wiley Periodicals, Inc.