Original Research ARTICLE
Identification of Candidate Protein Markers in Skeletal Muscle of Laminin-211-Deficient CMD Type 1A-Patients
- 1Department of Pediatric Neurology, Developmental Neurology and Social Pediatrics, University of Duisburg-Essen, Essen, Germany
- 2Leibniz-Institut für Analytische Wissenschaften -ISAS- e.V., Dortmund, Germany
- 3Department of Clinical Neuroscience, University of Cambridge, Cambridge, United Kingdom
Laminin-211 deficiency leads to the most common form of congenital muscular dystrophy in childhood, MDC1A. The clinical picture is characterized by severe muscle weakness, brain abnormalities and delayed motor milestones defining MDC1A as one of the most severe forms of congenital muscular diseases. Although the molecular genetic basis of this neurological disease is well-known and molecular studies of mouse muscle and human cultured muscle cells allowed first insights into the underlying pathophysiology, the definition of marker proteins in human vulnerable tissue such as skeletal muscle is still lacking. To systematically address this need, we analyzed the proteomic signature of laminin-211-deficient vastus muscle derived from four patients and identified 86 proteins (35 were increased and 51 decreased) as skeletal muscle markers and verified paradigmatic findings in a total of two further MDC1A muscle biopsies. Functions of proteins suggests fibrosis but also hints at altered synaptic transmission and accords with central nervous system alterations as part of the clinical spectrum of MDC1A. In addition, a profound mitochondrial vulnerability of the laminin-211-deficient muscle is indicated and also altered abundances of other proteins support the concept that metabolic alterations could be novel mechanisms that underline MDC1A and might constitute therapeutic targets. Intersection of our data with the proteomic signature of murine laminin-211-deficient gastrocnemius and diaphragm allowed the definition of nine common vulnerable proteins representing potential tissue markers.
Congenital muscular dystrophies (CMD) comprise a heterogenous group of genetically caused neuromuscular diseases with muscle weakness apparent at birth or in the first 6 months of life. Most of the different subtypes are of autosomal recessive inheritance (1). Laminin-211 (formerly merosin) -deficient CMD type 1A (MDC1A) is caused by recessive mutations in the LAMA2 gene (encoding for the α2 subunit of laminin-211) and constitutes approximately 10–30% of total CMD cases in the European population. Laminin-211 is expressed in the brain vasculature, the skeletal muscle basal lamina as well as in the myotendinous and neuromuscular junctions (2). White matter T2 signal hyperintensity as reflection of increased interstitial water content occurs in almost all patients after 6 months of age (3, 4). MDC1A-patients suffer from muscular weakness associated with elevated serum creatine kinase (CK) levels, poor suck and cry, multiple joint contractures and delayed motor development. Most of the MDC1A-patients never achieve independent ambulation (5–7). Extramuscular manifestations include seizures in 30% of patients, demyelinating neuropathy and CNS abnormalities such as polymicrogyria and cortical bandlike heterotopia. Mental retardation rarely occurs (8). Whereas patients with a complete deficiency of laminin-211 present with a severe clinical spectrum of the disease, a partial deficiency of the protein leads to milder phenotypes (6). In 2011, Gawlik and Durbeej speculated that the pathogenicity of LAMA2 mutations, which disrupt the assembly of the corresponding laminin-211 protein with other basal lamina components, explains the full penetrance of the phenotype (9).
Laminin is a cell-adhesion molecule localized to the basement membrane of skeletal muscle. The biological functions of laminins such as modulation of cytoskeleton and intracellular signaling pathways are accomplished via the interaction with transmembrane receptors which—in skeletal muscle—are represented by dystroglycan and integrin α7β1 as the two major receptors for laminin-211 (10, 11). In this context, laminin has also been postulated to protect the muscle fibers from damage under the constant stress of contractions (9). Notably, while laminin-211 (composed of α2, β1, and γ1 chains) was first isolated from placenta and originally called merosin (12), it is now well established that laminin-211 is the main laminin isoform in skeletal muscle (12–14). In this context, it is important to note that laminin-211 function has been linked to muscle development and through the formation of laminin networks also to cytoskeleton and intracellular signaling pathways. Moreover, it is believed that laminin-211 influences signal transmission events and muscle innervation via modulation of NMJ-integrity and function (9).
Although clinical features of LAMA2-patients, especially regarding the manifestation in skeletal muscle are well described, molecular signatures of muscle pathology defining protein markers of the effect of loss of functional laminin-211 remain elusive. Utilizing label free mass spectrometry-based protein quantification, we here systematically address laminin-211-deficiency-related protein changes to define tissue biomarkers in skeletal muscle of MDC1A-patients from the severe disease spectrum. This in turn can provide potential insights into the underlying pathophysiology and thus represent valuable outcome measures for therapeutic intervention concepts such as the emerging gene therapy (15). Hence, data presented in this article provides important molecular insights into MDC1A-related pathophysiology due to alterations in protein composition beyond the extra-cellular matrix, an important aspect for the definition of further or alternative therapeutic intervention concepts.
Materials and Methods
The muscle biopsy specimen investigated in this study have been initially collected for diagnostic purposes. Written informed consent was obtained from the participants (or rather their legal guardians) for the participation into the subsequent research as well as for publication of the findings (including any potentially-identifying information). Ethical approval was not required as per the local legislation.
Immunofluorescence-Based Studies in Muscle Biopsies
Five micrometer cryosections of muscle biopsy specimen were generated, dried and afterwards re-hydrated in phosphate-buffered saline (PBS) followed by exposure to primary antibodies (Table 1) diluted in 1% bovine serum albumin (BSA). After incubation of the sections in a wet-chamber for one hour at 25°C (incubator), those were washed three times in PBS. Next, the secondary antibodies (Alexa 488-conjugated) 1:500 diluted in 1% BSA were added and sections were again incubated in a wet-chamber for one hour at 25°C. Prior sections were finally covered in mounting medium (Dako), three washing steps in PBS were performed.
Immunoblot-Based Studies on Whole Muscle Protein Extracts
Twenty micrograms of patient skeletal muscle protein lysates were prepared using NuPAGE LDS sample buffer (Thermo) and NuPAGE reducing agent, denatured at room temperature for 10 minutes and loaded into 4–12% gradient SDS-PAGE gels (Thermo, NP0322). Proteins were separated and transferred to a PVDF membrane using an iBlot2 dry transfer system (Thermo). PVDF membranes were blocked using 5% milk in Tris buffer-saline with 1% tween-20 (TBS-T) at room temperature for 1 hour before probing using monoclonal antibodies targeting NDUFB8, SDHB, UQCRC2, COX2, and ATP5A (Abcam, ab110411) as well as antibodies targeting GAPDH (Abcam, ab8245) and VDAC1 (Abcam, ab14734). Membranes were washed in TBS-T and probed again for 1 h with HRP-conjugated anti-Mouse IgG (Abcam, ab97023) prior to chemiluminescent imaging. Quantification was performed using ImageJ.
Proteomic Profiling in Human Skeletal Muscle
Ammonium hydrogen carbonate (NH4HCO3), anhydrous magnesium chloride (MgCl2), guanidine hydrochloride (GuHCl), iodoacetamide (IAA), and urea were purchased from Sigma-Aldrich, Steinheim, Germany. Tris base was obtained from Applichem Biochemica, Darmstadt, Germany and Sodium dodecyl sulfate (SDS) was purchased from Carl Roth, Karlsruhe, Germany. Dithiothreitol (DTT), EDTA-free protease inhibitor (Complete Mini) tablets were obtained from Roche Diagnostics, Mannheim, Germany. Sodium chloride (NaCl) and calcium chloride (CaCl2) were from Merck, Darmstadt. Sequencing grade modified trypsin was from Promega, Madison, WI USA. Benzonase® Nuclease was purchased from Novagen. Bicinchoninic acid assay (BCA) kit was acquired from Thermo Fisher Scientific, Dreieich, Germany. All chemicals for ultra-pure HPLC solvents such as formic acid (FA), trifluoroacetic acid (TFA) and acetonitrile (ACN) were obtained from Biosolve, Valkenswaard, The Netherlands.
Cell Lysis, Sample Preparation and Trypsin Digestion
In total eight samples derived from four healthy controls (gender- and age-matched) and four MDC1A-patients were processes independently. All muscle samples were collected from the mid portion of vastus lateralis. Approximately 10 slices of 10 μm of muscle were lysed in 50 μL of lysis buffer (50 mM Tris-HCl (pH 7.8) 150 mM NaCl, 1% SDS, and Complete Mini) using a manual glass grinder. Then samples were centrifuged for 5 min at 4°C and 5.000 g. Protein concentration of the supernatant was determined by BCA assay (according to the manufacturer's protocol). In order to reduce the cysteines 10 mM of DTT were added to the samples followed by incubation at 56°C for 30 min. Next, the free thiol groups were alkylated with 30 mM IAA at room temperature (RT) in the dark for 30.
Sample digestion and cleanup were performed using filter-aided sample preparation (FASP) as described previously (16, 17) with some minor changes. Briefly, 100 μg of protein lysate was diluted 10-fold with freshly prepared 8 M urea/100 mM Tris-HCl (pH 8.5) buffer (18) and placed on a Microcon centrifugal device (30 KDa cutoff). Afterwards, the filter was centrifuged at 13,500 g at RT for 15 min (all the following centrifugation steps were performed under the same conditions). Three washing steps were carried out with 100 μL of 8 M urea/100 mM Tris-HCl (pH 8.5). For buffer exchange, the device was washed thrice with 100 μL of 50 mM NH4HCO3 (pH 7.8). Next, 100 μL of the digestion buffer (trypsin (Promega), 1:25 w/w, protease to substrate, (0.2 M GuHCl and 2 mM CaCl2 in 50 mM NH4HCO3 pH 7.8), was added to the filter which contains the bound proteins and the samples were incubated at 37°C for 14 h. Resulting tryptic peptides were recovered by centrifugation with 50 μL of 50 mM NH4HCO3 followed by 50 μL of ultra-pure water. In the final step, tryptic peptides were acidified by adding 5 μl of 10 % TFA (v/v). The digests were quality controlled as described previously (19).
Samples (technical duplicates) were measured using an Ultimate 3000 nano RSLC system coupled to an Orbitrap Fusion Lumos mass spectrometer (both Thermo Scientific) and analyzed in a randomized order to minimize systematic errors. Firstly, peptides were preconcentrated on a 100 μm × 2 cm C18 trapping column for 10 min using 0.1 % TFA (v/v) at a flow rate of 20 μL/min. Next the separation of the peptides was performed on a 75 μm × 50 cm C18 main column (both Pepmap, Thermo Scientific) with a 120 min LC gradient ranging from 3 to 35 % of 84 % ACN, 0.1 % FA (v/v) at a flow rate of 230 nL/min. MS1 spectra was acquired in the Orbitrap from 300 to 1,500 m/z at a resolution of 120,000 using the polysiloxane ion at m/z 445.12003 as lock mass (20), with maximum injection times of 50 ms. Next, top 10 most intense signals were selected for fragmentation by HCD with a collision energy of 30%. MS2 spectra were acquired in the ion trap at a resolution of 120,000, with maximum injection times of 300 ms and a dynamic exclusion of 15 s. The ACG target was set at 2.0 × 103 ions for MS2 and 2.0 × 105 for MS2.
Label Free Data Analysis
Data analysis of the acquired label free quantitative MS data was performed using the Progenesis LC-MS software from Non-linear Dynamics (Newcastle upon Tyne, U.K.). Alignment of MS raw data was conducted by Progenesis, which automatically selected one of the LC-MS files as reference. Next, peak picking was performed and only features within retention time and m/z windows from 0 to 120 min and 300–1,500 m/z, with charge states +2, +3, and +4 were considered for peptide statistics, analysis of variance (ANOVA). The MS/MS spectra were exported as peak lists which were searched against a concatenated target/decoy version of the mouse Uniprot database (Homo sapiens with 20273 entries, downloaded 22.07.2017) using Mascot 2.4 (Matrix Science), MS-GF+, and X!Tandem Jackhammer (2013.06.15) with the help of searchGUI 1.14.4 (21). Trypsin with a maximum of two missed cleavages was selected as enzyme. Carbamidomethylation of Cys was set as fixed and oxidation of Met was selected as variable modification. MS and MS/MS tolerances were set to 10 ppm and 0.5 Da, respectively. Combined search results were filtered at a false discovery rate (FDR) of 1 % on the protein level and exported using PeptideShaker 0.28.0 (http://code.google.com/p/peptide-shaker/). Data was reimported into Progenesis and peptide sequences containing oxidized Met were excluded for further analysis. Only proteins that were quantified with unique peptides were exported. For each protein, the average of the normalized abundances (obtained from Progenesis) was calculated and used to determine the ratios between patients and control. Only proteins which were (i) commonly quantified in all the replicates with (ii) one unique peptides, (iii) an ANOVA p-value of ≤0.05 (Progenesis) and (iv) an average ratio ≤ log2 0.99or ≥ log2−0.9 were considered as up respectively down regulated.
Data Plotting and Pathway Analysis
All data was plotted using Origin 6.0 and Adobe Illustrator. For pathway analysis, the GO ontology, KEGG and Reactome were used and data manually filtered for relevant pathways. The Proteomap was the online tool available (https://www.proteomaps.net/). The annotation of these proteomaps is based on the KEGG database platform, each protein is shown by a polygon, and functionally relevant proteins are arranged as neighbors. Additionally, polygon areas represent protein abundances weighted by protein size.
Clinical and Genetic Findings
We followed six patients diagnosed with CMD type 1A in the Department of Neuropediatrics of the University Children's Hospital, University Duisburg-Essen, in a tertiary care setting. Mutations in LAMA2 were found in all patients; however, one patient, who was severely affected, was lost of follow-up at the age of 2 years (patient 1).
Table 2 summarizes the molecular genetic, clinical and neuroradiologic findings. All six patients presented with symptoms within the first 4 months of life: poor head control, generalized muscular hypotonia and muscular weakness, poor spontaneous movements, delayed motor milestones, and high levels of creatine kinase (CK). The highest measured CK levels are listed in Table 2. The six children never achieved ambulation. Patient 5 was able to walk with support, but after a febrile virus infection with rhabdomyolysis at the age of 20 months she had a dramatic loss of motor function. Currently, at the age of 35 months, she showed a recovery with the ability to stand with support. All patients had dysphagia, recurrent chest infections and need varying levels of pulmonal support (i.e., cough assist). Additionally, joint contractures were present in all six patients. In patients 2–6, the underlying genetic mutations segregate with the phenotypes.
Four of our six patients underwent brain MRI, and three of the four showed white matter lesions (WMLs) on T2-weighted brain MRI images. The WMLs were characterized by a diffuse symmetrical distribution in cerebral areas that are normally myelinated in the developing brain (Figure 1) without involvement of thalamus or brainstem in both cases. Lissencephaly could not be detected in our patients. Patient 1 had a brain MRI performed at the age 4 months, which might be too early to detect WML in laminin-211-deficiency. All patients had normal cognitive development with no seizures being reported.
Figure 1. MRI findings in MDC1A-patients included in the study. Brain MRI TIRM (Turbo-Inversion Recovery-Magnitude, axial) show diffuse symmetrical white matter lesions predominantly in the periventricular regions and no involvement of thalamus or brainstem in both cases. (A) Patient 3 with expanded supratentorial leukencephalopathy at two years of age. (B) Patient 4 diffuse leukencephalopathy.
Muscle Biopsy Findings
Frozen biopsied muscle tissue samples from all six patients were analyzed and showed severe features of muscular dystrophy including increased extracellular connective tissue with fibrosis, cell necrosis, numerous central nuclei, and high fiber size variability (Figure 2A). All patients presented total absence of laminin-211 in our immunofluorescence-based microscopic examinations (Figure 2A). Moreover, laminin-511 was mild do moderately increased in the muscle biopsy specimen of the six patients as shown for two representative cases in Figure 2B. As controls, we selected subjects who had underwent a muscle biopsy, but who were found not to have any signs of a skeletal muscle disease.
Figure 2. Muscle biopsy findings in MDC1A-patients included in the study. (A) Muscle biopsy cryostat sections stained with hematoxylin and eosin (HE) and immunofluorescence staining with anti-laminin-211 antibodies reveals absent immunoreactivity in the patient-derived sections. (B) Immunofluorescence staining of muscle cryostat sections with anti-laminin-511 antibody reveals in increased immunoreactivity in muscle biopsies of MDC1A-patients as indicated for two representative cases shown here. Ctrl, control muscle; numbers 1, 2, 3, 4, 5, and 6 refer to muscle biopsies derived from MDC1A-patients with the corresponding patient numbers listed in Table 2.
Loss of LAMA2 Induces Changes in General Protein Composition in Human Skeletal Muscle
Proteomic profiling is a useful approach to obtain unbiased insights into the molecular etiology of diseases such as muscular disorders (24, 25). To identify, molecular marker proteins for muscle fiber vulnerability due to the loss of functional laminin-211 (protein encoded by LAMA2), the proteomic signature of muscle biopsy specimen of four MDC1A-patients (patients 1–4; see Table 2) from the severe spectrum of the disease was compared to the signature in four controls. Using liquid chromatography coupled to tandem mass spectrometry (LC-MS/MS) (label free protein quantification) (Figure 3), we quantified 1977 proteins out of which 86 proteins presented (4.3%) with statistically significant altered abundances (35 were increased and 51 decreased; Table 3). The considerable decrease of laminin-211 (−3.84; log2) in the patient-derived samples hereby reflects the sensitivity of our proteomic profiling approach as well as the reliability of the data. To obtain insights into the molecular etiopathology of MDC1A, pathway analyses were performed utilizing DAVID (https://david.ncifcrf.gov/), KEGG (www.genome.jp/kegg/pathway.html), Reactome (www.reactome.org) and Proteomap (www.proteomaps.net). In addition, information concerning subcellular localization and function has been extracted from uniport (www.uniprot.org) for each protein altered in abundance to provide a complete picture of the molecular basis of MDC1A manifestation in skeletal muscle (Table 3). Data analysis via diverse pathway analysis tools revealed potential changes in axon guidance, cell cycle and DNA-repair, extracellular matrix organization, metabolism of fatty acids, sugar, lipids and proteins. In addition, protein changes seem to impact on vesicular transport machinery, TCA and respiratory electron transport as well as proper muscle contraction (Figure 3C). Data analysis via Proteomap confirmed the impact of detected protein changes on cell cycle, cellular metabolism, muscle contraction as well as on the composition of the extracellular matrix. Based on the functional properties of some of the proteoglycans with altered abundances such as COL1A1 and -A2, LAMA4 and ITGA6, changes in PI3K/Akt-signaling as a pathway involved in a variety of cellular functions has been highlighted by the Proteomaps-representation (Figure 3D). This molecular observation accords with a postulated function of laminin-211 in signal transmission (9).
Figure 3. Proteomic profiling of four MDC1A-patient derived muscles. (A) Methodological workflow applied in the study. (B) Results of our label-free proteomic profiling are shown as a volcano plot. All points (each one represents a protein) over the horizontal line have a statistically significant p-ANOVA of ≤0.05. In total, 51 proteins are decreased (red points) and 35 proteins are increased in abundance (green points) in laminin-211-deficient skeletal muscle. (C) Results of an in silico pathway analysis of proteins altered in abundance utilizing DAVID, KEGG, and Reactome platforms indicating alterations of ECM protein composition, cellular metabolism, axon guidance, cell cycle, and DNA repair, muscle contraction and of vesicular protein transport. (D) Results of an in silico pathway analysis of proteins altered in abundance utilizing Proteomap platform confirming the perturbed cellular functions indicated by the results of the previous pathway analysis and moreover indicating perturbed function of cellular signaling cascades such as PI3K-Akt signaling.
Perturbed proteins represent tissue markers for LAMA2-related congenital myopathy and—based on information available in uniport—mainly localize to the sarcoplasm, to nuclei, the SR and to mitochondria thus suggesting a global unified organelle vulnerability in LAMA2-mutant skeletal muscle with a predominance to mitochondria. In addition, functions of several of the accord proteins accord with the known involvement of laminin-211 in ECM-composition, signal transduction, cytoskeleton, mitochondrial homeostasis and muscle cell development. Interestingly, functions of a variety of proteins increased in abundance hint toward the activation of neuroprotective mechanisms including neuromuscular transmission.
Altered Abundance of Paradigmatic Proteins and Mitochondrial Vulnerability Can Be Confirmed in Independent Muscle Biopsies Derived From MDC1A-Patients
To study the potential of altered proteins as tissue markers and reliability of our proteomic findings further immunofluorescence studies were performed on the muscle biopsies derived from the four patients included in the proteomic profiling (patients 1–4) as well as on biopsies derived from two further MDC1A-cases (patients 5 and 6). Doing so, we focused on abundances and distribution of agrin, glypican-1, and glucose-6-phosphate isomerase as paradigmatic proteins. Immunofluorescence-based studies of agrin showed occasionally focal sarcoplasmic accumulations most likely leading to the detected increase of overall agrin protein level (Table 3) as shown for two representative cases in Figure 4A. In comparison to the controls which show an enrichment of glypican-1 at the sarcolemma, in MDC1A-patient derived biopsies, small sarcoplasmic dots immunoreactive for the protein as well as an enrichment of these dots at the (sub-)sarcolemmal region could be observed. This irregular protein-distribution is shown for two representative cases in Figure 4A and most likely causes the 3.25-fold (log2 ratio) increase detected by proteomic profiling. Moreover, immunofluorescence studies confirmed the decreased expression of glucose-6-phosphate isomerase in muscle biopsies derived from four MDC1A-cases (shown for two representative cases in Figure 4A).
Figure 4. Confirmation study of proteomic findings in independent MDC1A-patient derived muscle biopsies. (A) Immunofluorescence studies of spectrin have been carried out to visualize sarcolemma in muscle cryostat sections of control and patient-derived muscles and to show similar abundance and localization of the protein. Further immunofluorescence studies of glucose-6-phosphate isomerase reveal a focal dot-like sarcoplasmic staining/ enrichment in muscle biopsy specimen derived of controls (one representative control is shown here) but show a decrease in the biopsies of two representative MDC1A-patients. Immunofluorescence studies of glypican-1 showed focal enrichment at the sarcolemma of control muscle (one representative control is shown). However, in the LAMA2-diseased muscle fibers, small dots sarcoplasmic, and (sub-)sarcolemmal dots immunoreactive for glypican-1 were detected. Immunofluorescence studies of agrin revealed of a sarcoplasmic enrichment/ aggregation in MDC1A-patient derived muscle fibers (white arrows) compared to control muscle fibers. (B) Histological NADH-TR staining revealed “sarcoplasmic gaps” (white arrows in patient 2; second and third column) and reduced enzyme activity (white arrows in patient 2, 3, 5, and 6) as well as focal enrichment in some myofibers (black arrows in patient 2, 3, and 6). Modified Gomori Trichrome staining revealed sarcoplasmic areas with minor staining (white arrows in patient 2, 3, 5, and 6) as well as areas of sarcoplasmic and subsarcolemmal enrichment (black arrows in patient 2,3, and 5). (C) Immunoblot analysis of proteins representative for the respiratory chain complexes I-V and of VDAC1 (porin) revealed a statistically significant decrease of complexes III and V as well as of VDAC1 in MDC1A-patient derived protein muscle extracts compared to controls (diagram). GAPDH has been used as loading control and proteins have been quantified against GAPDH (ns, not significant; *statistically significant, **statistically very significant).
Given that our proteomic data revealed a profound vulnerability of NADH dehydrogenases, NADH muscle histology (NADH-TR to highlight the oxidative enzyme activity) has been investigated in biopsies of six MDC1A-cases and revealed no structural abnormalities like cores, whorles or lobulated fibers but reduced staining and even gaps in some muscle fibers indicative for reduced enzyme activity (shown for four representative cases in Figure 4B). Modified Gomori Trichrome staining has additionally been performed in these patients to study mitochondrial abnormalities and occasionally revealed both, focal accumulations and reduced staining suggesting irregular mitochondrial distribution (shown for four representative cases in Figure 4B). Immunofluorescence studies of Glucose-6-phosphate isomerase (GPI), glypican-1 and agrin have been carried out twice with similar results. Prompted by the proteomic results suggestive for vulnerability of the respiratory chain, respective complexes have been investigated by immunoblotting utilizing antibodies targeting proteins localized to the respective complexes and muscle protein extracts derived from two controls and three MDC1A-patients. Results of these studies confirm a vulnerability of the respiratory chain complexes (Figure 4C) and thus support our proteomic findings which were indicative for mitochondrial vulnerability. Immunoblot studies have been carried out twice with similar results.
A complete loss of laminin-211, which is encoded by LAMA2, causes a severe congenital muscular dystrophy with onset of symptoms in the first few months of life. Partial laminin-211 deficiency can be caused not only by primarily LAMA2 mutations, but also secondarily by other muscular dystrophies, including dystroglycanopathy (26, 27). This in turn suggests the need to define tissue marker proteins for patients with classical MDC1A based on LAMA2 mutations leading to laminin-211 deficiency. To systematically address this need, proteomic profiling has been carried out using muscle biopsy specimens derived from MDC1A patients with myopathy and brain malformations. Results not only confirm laminin-211 deficiency but also accord with elevated CK level in the patients as intramuscular CK was significantly decreased in the patient-derived muscle. In addition, decrease of Mth938 domain-containing protein (involved in preadipocyte differentiation and adipogenesis) as well as of prostaglandin reductase 2 (involved in inhibition of adipocyte differentiation) might provide a molecular link to a replacement of degenerating muscle fibers by fatty tissue. Based on the considerable increase of proteins localized to the extracellular matrix (Table 3) our data moreover confirm proliferation of connective tissue as a pathophysiological hallmark of fibrosis observed in the muscle biopsy specimen of our patients (Figure 1). The increase of LAMA4 (laminin subunit α4/ Lm-411) as well as of LAMA5 (laminin subunit α5/ Lm-511; identified by our routine diagnostic staining (Figure 2B) and 1.7-fold (log2 ratio) increased with a p-Anova of 0.19 in our proteomic findings; data not shown) most likely reflects a cellular attempt to (partially) compensate the loss of laminin-211 and to avoid a complete breakdown of muscle fibers. This assumption is not only supported by the results of a previous study showing that laminin-111 protein therapy reduces muscle pathology and improves viability of a MDC1A mouse model (28) but also by the known important role of LAMA4 for NMJ-integrity (29), the concomitant increase of integrin alpha-6/beta-1, a receptor for laminin and the increase of the basal cell adhesion molecule/ laminin alpha-5 receptor (Table 3) to optimize laminin binding to laminin-211-deficient muscle fibers. However, despite the increase of LAMA4 in MDC1A, the resultant laminin-protein complex is known to bind merely poorly to integrins and α-dystroglycan (30).
Given that subtle neuromuscular junction (NMJ) defects have been reported in laminin α2 chain-deficient mice (31), the detected abnormal sarcoplasmic accumulation and thus decreased release to the synaptic cleft of muscle agrin, a basal lamina glycoprotein crucial for the formation and the maintenance of NMJs, might contribute to impaired neuromuscular transmission as part of the MDC1A-pathophysiology. Along this line, increase of glypican-1, might indicate activation of a compensatory mechanism to avoid NMJ-breakdown resulting from profound de-innervation. Notably, mini-agrin has been shown to bind to the basement membrane and the DGC via α-dystroglycan and thus ameliorate muscle pathology in vivo (32) and transgenic expression of mini-agrin (contains binding sites for laminins and α-dystroglycan) and αLNNd (recombinant protein linking LAMA4 to α-dystroglycan) in a mouse model for MDC1A fully restored basement membrane stability. This effect resulted in recovery of muscle force and size leading to increased overall body weight, and extended life span (33). Moreover, glypican-1 has recently been linked to the pathophysiology of a muscular dystrophy complicated by a myasthenic syndrome (34) as well as of laminin α4, which has an important role in NMJ-integrity (29). However, the conclusion of endogenous activation of compensatory mechanisms in muscle of MDC1A-patients is further supported by the increase of C4b-binding protein alpha chain belonging to the complement system that deposits its activation products on innervated motor end-plates in ALS-patients (35).These combined findings in turn underline the significance our proteomic profiling data aiding the identification of marker proteins with pathophysiological relevance as well as impacts for attempts to develop new treatment strategies for MDC1A.
However, as the NMJ represents a paradigmatic synapse, altered abundance of a variety of other proteins including 14-3-3F, AKR7A2, protein FAM162A, GTPC1, NIBAN, PRVA, SRC8, Succinate-semialdehyde dehydrogenase, TTHY (Table 3), might correlate with the broad activation of neuroprotective mechanisms in the etiopathology of MDC1A and thus also link to the CNS alterations of the children. NudC domain-containing protein 2 modulates the LIS1/dynein pathway via stabilization of lissencephaly protein 1 (LIS1) with Hsp90, a cellular chaperone. Importantly, the p.L279P of NudC domain-containing protein 2 influences LIS1 stability (36). Given that lissencephaly has been described in MDC1A-patients (37) increase of this protein not only accords with the concept of activation of proteins involved in neuroprotection (in none of our patients lissencephaly was found via MRI) but also provides a first molecular hint to the manifestation of lissencephaly in the etiopathology of MDC1A. Neuronal vulnerability associated with altered synaptic transmission—as one pathophysiological cascade among others—is indicated by the decrease of glucose-6-phosphate isomerase, acting as a neurotrophic factor for spinal and sensory neurons (38), of phosphatidylethanolamine-binding protein 1 involved in the function of the presynaptic cholinergic neurons of the central nervous system, of dihydropyrimidinase-related protein 2 involved in neuronal development and polarity and in axon growth and guidance (Table 3). Activation of protective mechanisms toward maintenance of muscle contraction is for instance indicated by increase of ADDG, CNN3 and myosin light chain 4 as well as MYL6B (Table 3). Increased abundance of sarcolemmal membrane-associated protein involved in myoblast fusion (Table 3) might hint to activation or muscle fiber regeneration.
Interestingly, the results of our proteomic profiling indicate a profound vulnerability of mitochondria in laminin-211-deficient muscle based on the downregulation of a variety of mitochondrial proteins (Table 3 and Figure 4C) and reduced mitochondrial respiration and ATP production in MDC1A-patient derived muscle cells has been reported as the result of changes in abundances of transcripts encoding for mitochondrial key players (39). Among the dysregulated expression of genes related to energy production in myotubes, Fontes-Oliveira and colleagues described reduced level of NDUFA8 and decrease of NADH dehydrogenase [ubiquinone] 1 alpha subcomplex subunit 8 encoded by this transcript has been identified in our study (Table 3). Thus, our proteomic data confirm an impairment of the mitochondrial bioenergetic status on the protein level in human MDC1A muscle cells also in complex laminin-211-deficient muscle progressed in the etiopathology. Clincally, one of our patients (patient 5) presented with rhabdomylosis and loss of the ability to stand and walk with support after a febrile infection. One might speculate that this unusual disease course also hints toward severe metabolic impairment.
As our study is based on the utilization of material derived from patients with progressed status of the disease, it is not possible to clearly differentiate between primary pathophysiological events and molecular changes occurring as a phenomenon secondary to the primary cascades. However, intersection with data obtained from skeletal muscle derived from murine animal models with early (and unified) manifestation of the disease might harbor some limitations regrading general validity of data comparison but still allows to define candidate marker proteins of general significance involved in the initiate pathophysiology. De Oliveira and co-workers studied the proteomic signature of diaphragm and gastrocnemius muscle derived from dy3K/dy3K mice (MDC1A mouse model). Out of the approximately 700 identified proteins, 113 and 101 respectively, were differentially expressed in the investigated tissues (40). Notably, there was no overlap between the proteins increased in the muscle of the mouse model and the patients although few proteins with elevated levels in murine laminin-211-deficient muscle are related to extracellular matrix composition and might hint to early stages of fibrosis. Interestingly, parvalbumin alpha has been identified as decreased in their study but as increased in our proteomic profiling of MDC1A-patient derived quadriceps muscle. As in skeletal muscle, parvalbumin is thought to be involved in relaxation after contraction, its increase might accord with the activation of compensatory mechanisms upon disease progression. However, on a general note, one might speculate that either the differences in the investigated muscles (diaphragm and gastrocnemius in mice and vastus in human) or the progressed disease in our patients compared to the young-aged animals (four weeks) might explain the missing overlap between increased proteins in mice and human. Nevertheless, when focussing on the decreased proteins, a variety of different ones are downregulated in both, human and mice thus representing promising tissue markers for MDC1A across species. These proteins include creatine kinase M-type, voltage-dependent anion-selective channel protein 1, glucose-6-phosphate isomerase, cytoplasmic malate dehydrogenase, cytochrome C oxidase subunit 5A, mitochondrial NADH-dehydrogenase [ubiquinone] flavoprotein 1 and dihydrolipollysine-residue succinyltransferase component of 2-oxoglutarate dehydrogenase complex as well as NADH dehydrogenase [ubiquinone] 1β subcomplex subunit 10 and 2-oxoglutarate dehydrogenase.
The intersected data is indeed indicative of defective metabolism—especially mitochondrial function—in laminin-211-deficient muscle cells of MDC1A patients and hence confirm the findings of Fontes-Oliveira and colleagues in mice (39). This is moreover underlined by the results of our histological NADH-TR and Gomori Trichome staining supporting the concept of impaired mitochondrial activity. Along this line, the findings of both studies support the hypothesis that skeletal muscle metabolism might be a promising pharmacological target to improve muscle function, energy efficiency and tissue maintenance of MDC1A (39). Notably, in a pre-clinical study utilizing Lama2-deficient mice, mitochondria have been proven to represent a promising therapeutic target (41).
Analysis of the proteomic signature of proximal muscle derived from four children suffering from MDC1A allowed the identification of 86 proteins with altered abundance and potential pathophysiological impact. Focussing on molecular changes in skeletal muscle, a variety of affected proteins are linked to the vulnerability of the central nervous system, especially to altered synaptic transmission in the disease caused by MDC1A. Moreover, this suggests that therapeutic intervention targeting the synaptic dysfunction might represent a promising element of the overall concept in the treatment of MDC1A as already shown in mice treated with agrin. Although our patients do not present with lissencephaly, increased NudC domain-containing protein 2 might represent a potential protein modifier of the central nervous system phenotype of MDC1A. Proteomic signature of laminin-211-deficient muscle moreover indicated a profound mitochondrial vulnerability with predominant decrease of proteins belonging to complex 1. A comparison of our data with the proteomic signature of laminin-211-deficient distal muscle and diaphragm from mice allowed the identification of nine decreased proteins representing potential markers for MDC1A. However, more comprehensive studies utilizing skeletal muscle derived from MDC1A-patients (sub-cohorts at different stages of the disease) are needed to ultimately define vulnerable proteins as definite MDC1A tissue markers.
Contribution to the Field Statement
MDC1A is the most common form of congenital muscular dystrophies. The disease is caused by mutations in the LAMA2 gene and molecular studies of mouse muscle and human cultured muscle cells already allowed first insights into the underlying pathophysiology. However, the definition of (candidate) marker proteins in human skeletal muscle is still lacking. To address this gap of knowledge, we conducted a study to investigate the proteomic signature of laminin-211-deficient vastus muscle obtained from four MDC1A-patients. Results of our unbiased screening allowed the identification of potential tissue marker proteins which might also be involved in the pathophysiology of the disease.
Ethics commitee of Duisburg-Essen University.
HK conceptualized and designed the study, drafted the first version of the manuscript, interpreted results, and was principally responsible for the final content. DH carried out the initial analysis and interpretation of results, contributed to the discussion. AR carried out the initial analysis, conceptualized, and designed the study, reviewed and revised the manuscript for important intellectual content. RH and MJ performed the analysis of mitochondrial proteins. US reviewed and revised the manuscript for important intellectual content.
This work was supported by a grant from the French Muscular Dystrophy Association (AFM-Téléthon; #21644; grant to AR). Financial support was moreover obtained from the Ministerium für Innovation, Wissenschaft und Forschung des Landes Nordrhein-Westfalen, from the Senatsverwaltung für Wirtschaft, Technologie und Forschung des Landes Berlin and from the Bundesministerium für Bildung und Forschung.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
1. Philpot J, Sewry C, Pennock J, Dubowitz V. Clinical phenotype in congenital muscular dystrophy: correlation with expression of merosin in skeletal muscle. Neuromuscul Disord. (1995) 5:301–5. doi: 10.1016/0960-8966(94)00069-L
3. Tan E, Topaloglu H, Sewry C, Zorlu Y, Naom I, Erdem S, et al. Late onset muscular dystrophy with cerebral white matter changes due to partial merosin deficiency. Neuromuscul Disord. (1997) 7:85–9. doi: 10.1016/S0960-8966(96)00421-X
5. Herrmann R, Straub V, Meyer K, Kahn T, Wagner M, Voit T. Congenital muscular dystrophy with laminin alpha 2 chain deficiency: identification of a new intermediate phenotype and correlation of clinical findings to muscle immunohistochemistry. Eur J Pediatr. (1996) 155:968–76. doi: 10.1007/BF02282889
6. Cohn RD, Herrmann R, Sorokin L, Wewer UM, Voit T. Laminin alpha2 chain-deficient congenital muscular dystrophy: variable epitope expression in severe and mild cases. Neurology. (1998) 51:94–100. doi: 10.1212/WNL.51.1.94
7. Pegoraro E, Marks H, Garcia CA, Crawford T, Mancias P, Connolly AM, et al. Laminin alpha2 muscular dystrophy: genotype/phenotype studies of 22 patients. Neurology. (1998) 51:101–10. doi: 10.1212/WNL.51.1.101
8. Xiong H, Tan D, Wang S, Song S, Yang H, Gao K, et al. Genotype/phenotype analysis in Chinese laminin-alpha2 deficient congenital muscular dystrophy patients. Clin Genet. (2015) 87:233–43. doi: 10.1111/cge.12366
12. Ehrig K, Leivo I, Argraves WS, Ruoslahti E, Engvall E. Merosin, a tissue-specific basement membrane protein, is a laminin-like protein. Proc Natl Acad Sci USA. (1990) 87:3264–8. doi: 10.1073/pnas.87.9.3264
13. Leivo I, Engvall E. Merosin, a protein specific for basement membranes of Schwann cells, striated muscle, and trophoblast, is expressed late in nerve and muscle development. Proc Natl Acad Sci USA. (1988) 85:1544–8. doi: 10.1073/pnas.85.5.1544
14. Patton BL, Miner JH, Chiu AY, Sanes JR. Distribution and function of laminins in the neuromuscular system of developing, adult, and mutant mice. J Cell Biol. (1997) 139:1507–21. doi: 10.1083/jcb.139.6.1507
15. Kemaladewi DU, Maino E, Hyatt E, Hou H, Ding M, Place KM, et al. Correction of a splicing defect in a mouse model of congenital muscular dystrophy type 1A using a homology-directed-repair-independent mechanism. Nat Med. (2017) 23:984–9. doi: 10.1038/nm.4367
19. Burkhart JM, Schumbrutzki C, Wortelkamp S, Sickmann A, Zahedi RP. Systematic and quantitative comparison of digest efficiency and specificity reveals the impact of trypsin quality on MS-based proteomics. J Proteomics. (2012) 75:1454–62. doi: 10.1016/j.jprot.2011.11.016
20. Olsen JV, De Godoy LMF, Li G, Macek B, Mortensen P, Pesch R, et al. Parts per million mass accuracy on an orbitrap mass spectrometer via lock mass injection into a C-trap. Mol Cell Proteomics. (2005) 4:2010–21. doi: 10.1074/mcp.T500030-MCP200
21. Vaudel M, Barsnes H, Berven FS, Sickmann A, Martens L. SearchGUI: an open-source graphical user interface for simultaneous OMSSA and X!Tandem searches. Proteomics. (2011) 11:996–9. doi: 10.1002/pmic.201000595
22. Geranmayeh F, Clement E, Feng LH, Sewry C, Pagan J, Mein R, et al. Genotype-phenotype correlation in a large population of muscular dystrophy patients with LAMA2 mutations. Neuromuscul Disord. (2010) 20:241–50. doi: 10.1016/j.nmd.2010.02.001
23. Siala O, Louhichi N, Triki C, Moriniere M, Rebai A, Richard P, et al. Severe MDC1A congenital muscular dystrophy due to a splicing mutation in the LAMA2 gene resulting in exon skipping and significant decrease of mRNA level. Genet Test. (2007) 11:199–207. doi: 10.1089/gte.2006.0517
24. Gonzalez Coraspe JA, Weis J, Anderson ME, Munchberg U, Lorenz K, Buchkremer S, et al. Biochemical and pathological changes result from mutated Caveolin-3 in muscle. Skelet Muscle. (2018) 8:28. doi: 10.1186/s13395-018-0173-y
25. Roos A, Thompson R, Horvath R, Lochmuller H, Sickmann A. Intersection of proteomics and genomics to “solve the unsolved” in rare disorders such as neurodegenerative and neuromuscular diseases. Proteomics Clin Appl. (2018) 12. doi: 10.1002/prca.201700073
26. Hayashi YK, Engvall E, Arikawa-Hirasawa E, Goto K, Koga R, Nonaka I, et al. Abnormal localization of laminin subunits in muscular dystrophies. J Neurol Sci. (1993) 119:53–64. doi: 10.1016/0022-510X(93)90191-Z
27. Chan SH, Foley AR, Phadke R, Mathew AA, Pitt M, Sewry C, et al. Limb girdle muscular dystrophy due to LAMA2 mutations: diagnostic difficulties due to associated peripheral neuropathy. Neuromuscul Disord. (2014) 24:677–83. doi: 10.1016/j.nmd.2014.05.008
28. Rooney JE, Knapp JR, Hodges BL, Wuebbles RD, Burkin DJ. Laminin-111 protein therapy reduces muscle pathology and improves viability of a mouse model of merosin-deficient congenital muscular dystrophy. Am J Pathol. (2012) 180:1593–602. doi: 10.1016/j.ajpath.2011.12.019
29. Patton BL, Cunningham JM, Thyboll J, Kortesmaa J, Westerblad H, Edstrom L, et al. Properly formed but improperly localized synaptic specializations in the absence of laminin alpha4. Nat Neurosci. (2001) 4:597–604. doi: 10.1038/88414
30. Talts JF, Sasaki T, Miosge N, Gohring W, Mann K, Mayne R, et al. Structural and functional analysis of the recombinant G domain of the laminin alpha4 chain and its proteolytic processing in tissues. J Biol Chem. (2000) 275:35192–9. doi: 10.1074/jbc.M003261200
31. Desaki J, Matsuda S, Sakanaka M. Morphological changes of neuromuscular junctions in the dystrophic (dy) mouse: a scanning and transmission electron microscopic study. J Electron Microsc. (1995) 44:59–65.
32. Moll J, Barzaghi P, Lin S, Bezakova G, Lochmuller H, Engvall E, et al. An agrin minigene rescues dystrophic symptoms in a mouse model for congenital muscular dystrophy. Nature. (2001) 413:302–7. doi: 10.1038/35095054
33. Reinhard JR, Lin S, Mckee KK, Meinen S, Crosson SC, Sury M, et al. Linker proteins restore basement membrane and correct LAMA2-related muscular dystrophy in mice. Sci Transl Med. (2017) 9:eaal4649. doi: 10.1126/scitranslmed.aal4649
35. Bahia El Idrissi N, Bosch S, Ramaglia V, Aronica E, Baas F, Troost D. Complement activation at the motor end-plates in amyotrophic lateral sclerosis. J Neuroinflammation. (2016) 13:72. doi: 10.1186/s12974-016-0538-2
36. Zhu XJ, Liu X, Jin Q, Cai Y, Yang Y, Zhou T. The L279P mutation of nuclear distribution gene C (NudC) influences its chaperone activity and lissencephaly protein 1 (LIS1) stability. J Biol Chem. (2010) 285:29903–10. doi: 10.1074/jbc.M110.105494
37. Kubota A, Ishiura H, Mitsui J, Sakuishi K, Iwata A, Yamamoto T, et al. A homozygous LAMA2 mutation of c.818G>A caused partial merosin deficiency in a Japanese patient. Intern Med. (2018) 57:877–82. doi: 10.2169/internalmedicine.9588-17
38. Haga A, Niinaka Y, Raz A. Phosphohexose isomerase/autocrine motility factor/neuroleukin/maturation factor is a multifunctional phosphoprotein. Biochim Biophys Acta. (2000) 1480:235–44. doi: 10.1016/S0167-4838(00)00075-3
39. Fontes-Oliveira CC, Steinz M, Schneiderat P, Mulder H, Durbeej M. Bioenergetic impairment in congenital muscular dystrophy type 1A and leigh syndrome muscle cells. Sci Rep. (2017) 7:45272. doi: 10.1038/srep45272
40. De Oliveira BM, Matsumura CY, Fontes-Oliveira CC, Gawlik KI, Acosta H, Wernhoff P, et al. Quantitative proteomic analysis reveals metabolic alterations, calcium dysregulation, and increased expression of extracellular matrix proteins in laminin alpha2 chain-deficient muscle. Mol Cell Proteomics. (2014) 13:3001–13. doi: 10.1074/mcp.M113.032276
Keywords: laminin-211, laminin-α2, congenital muscular dystrophy, agrin, NudC domain-containing protein 2, muscle proteomics
Citation: Kölbel H, Hathazi D, Jennings M, Horvath R, Roos A and Schara U (2019) Identification of Candidate Protein Markers in Skeletal Muscle of Laminin-211-Deficient CMD Type 1A-Patients. Front. Neurol. 10:470. doi: 10.3389/fneur.2019.00470
Received: 11 October 2018; Accepted: 17 April 2019;
Published: 07 May 2019.
Edited by:Stefan Bittner, Johannes Gutenberg University Mainz, Germany
Reviewed by:Michele H. Jacob, Tufts University, United States
Dione Kobayashi, Independent Researcher, MA, Cambridge, United States
Copyright © 2019 Kölbel, Hathazi, Jennings, Horvath, Roos and Schara. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.