- Open Access
Dystrophin-deficient dogs with reduced myostatin have unequal muscle growth and greater joint contractures
Skeletal Musclevolume 6, Article number: 14 (2016)
Myostatin (Mstn) is a negative regulator of muscle growth whose inhibition promotes muscle growth and regeneration. Dystrophin-deficient mdx mice in which myostatin is knocked out or inhibited postnatally have a less severe phenotype with greater total mass and strength and less fibrosis and fatty replacement of muscles than mdx mice with wild-type myostatin expression. Dogs with golden retriever muscular dystrophy (GRMD) have previously been noted to have increased muscle mass and reduced fibrosis after systemic postnatal myostatin inhibition. Based partly on these results, myostatin inhibitors are in development for use in human muscular dystrophies. However, persisting concerns regarding the effects of long-term and profound myostatin inhibition will not be easily or imminently answered in clinical trials.
To address these concerns, we developed a canine (GRippet) model by crossbreeding dystrophin-deficient GRMD dogs with Mstn-heterozygous (Mstn +/−) whippets. A total of four GRippets (dystrophic and Mstn +/−), three GRMD (dystrophic and Mstn wild-type) dogs, and three non-dystrophic controls from two litters were evaluated.
Myostatin messenger ribonucleic acid (mRNA) and protein levels were downregulated in both GRMD and GRippet dogs. GRippets had more severe postural changes and larger (more restricted) maximal joint flexion angles, apparently due to further exaggeration of disproportionate effects on muscle size. Flexors such as the cranial sartorius were more hypertrophied on magnetic resonance imaging (MRI) in the GRippets, while extensors, including the quadriceps femoris, underwent greater atrophy. Myostatin protein levels negatively correlated with relative cranial sartorius muscle cross-sectional area on MRI, supporting a role in disproportionate muscle size. Activin receptor type IIB (ActRIIB) expression was higher in dystrophic versus control dogs, consistent with physiologic feedback between myostatin and ActRIIB. However, there was no differential expression between GRMD and GRippet dogs. Satellite cell exhaustion was not observed in GRippets up to 3 years of age.
Partial myostatin loss may exaggerate selective muscle hypertrophy or atrophy/hypoplasia in GRMD dogs and worsen contractures. While muscle imbalance is not a feature of myostatin inhibition in mdx mice, findings in a larger animal model could translate to human experience with myostatin inhibitors.
Duchenne muscular dystrophy (DMD) affects approximately 1 in 5000 newborn boys in whom DMD gene mutations lead to progressive degeneration of cardiac and skeletal muscle [1–4]. One strategy for promoting muscle regeneration involves inhibiting myostatin (Mstn; growth and differentiation factor 8 (GDF-8)), a negative regulator of muscle growth [5–8]. Humans , cattle , sheep , and dogs  with Mstn gene mutations have dramatic muscle hypertrophy. Dystrophin-deficient mdx mice in which myostatin is knocked out (Mstn −/−)  or inhibited postnatally [14, 15] also have a less severe phenotype with greater absolute force and less fibrosis of individual muscles.
Based on these findings, there has been increasing interest in treatments to inhibit myostatin and thus promote muscle growth [16, 17]. However, other studies have identified potential negative consequences to myostatin inhibition in skeletal muscle. Muscle tendons of Mstn −/− mice are hypocellular and more brittle . There are differential muscle effects in myostatin-null mice with, for example, the fast-twitch predominant extensor digital longus (EDL) muscle demonstrating reduced specific isometric force (force normalized by cross-sectional area (CSA)) and greater eccentric contraction decrement compared to the slow twitch soleus [19–21]. Elimination of myostatin in the dy(W) laminin alpha2-deficient murine model of congenital muscular dystrophy was associated with increased pre-weaning mortality, potentially due to reduced fat formation , while blockade of myostatin with transgenic expression of follistatin in the Dyf−/− model of LGMD2B and Myoshi myopathy exacerbated muscle degeneration with aging . Additional questions have been raised about potential exhaustion of the pool of muscle progenitor cells (i.e., satellite cells) undergoing multiple divisions in the absence of myostatin in muscular dystrophy .
While genetically engineered mice have provided an extremely powerful tool to study the molecular pathogenesis of disease [25, 26], results do not necessarily extrapolate to humans, presumably due to differences between murine and human size and physiology . These shortcomings are partially countered with canine models, which have been used extensively to study disease pathogenesis and treatment efficacy [28, 29]. This trend towards the use of dogs as models will likely accelerate with the recent sequencing of the canine genome . We, and others, have investigated the potential therapeutic role of myostatin inhibition in dogs. Adeno-associated virus (AAV8)-mediated over expression of the inhibitory myostatin propeptide was shown to enhance muscle growth in normal dogs . Analogous results were demonstrated in dystrophin-deficient golden retriever muscular dystrophy (GRMD) dogs . After 13 months, treated GRMD dogs had increased muscle weights, ranging from 49 % (tibialis cranialis) to 27 % (EDL), and a suggestion of reduced fibrosis by histochemistry.
To further study the effects of prolonged loss of myostatin in a large animal model of DMD, we developed a myostatin-deficient GRMD (GRippet) model by crossbreeding GRMD dogs with Mstn-heterozygous (Mstn +/−) whippets . GRippet and Mstn +/+ wild-type GRMD dogs, together with non-dystrophic control littermates, were assessed with functional tests, magnetic resonance imaging (MRI), and molecular/pathologic studies.
Dogs from a colony at the University of North Carolina at Chapel Hill (UNC-CH) received care and were assessed according to principles outlined in the National Research Council Guide for the Care and Use of Laboratory Animals. Studies were approved by the UNC-CH Institutional Animal Care and Use Committee (IACUC) through two protocols, UNC IACUC 08-103, Cross Breeding of Muscular Dystrophy and Myostatin-Null Dogs, and UNC IACUC 09-351, Standard Operating Procedures—Canine X-Linked Muscular Dystrophy. Independent of myostatin status, the dystrophic phenotype was determined based on elevation of serum creatine kinase (CK) and characteristic clinical signs. Genotype was confirmed by polymerase chain reaction (PCR) if CK results were ambiguous (Table 1).
All GRMD colonies worldwide have been derived from the same founder dog . This dog was initially used to establish a colony at Cornell University . The original dystrophic golden retriever was crossed with other breeds. Accordingly, GRMD is not a disease of purebred golden retriever dogs. We use the term GRMD to refer to dogs carrying the underlying splice site mutation [35, 36].
Whippet dogs homozygous for the myostatin-null allele (Mstn −/−) (so called bully whippets) have gross enlargement of muscle, while those that are heterozygous (Mstn +/−) have intermediate muscle mass and are significantly faster than their myostatin wild-type counterparts . Semen from a sire (Mstn +/−) of bully whippet dogs (Mstn −/−) was collected for subsequent artificial insemination of a GRMD carrier to produce a first litter. The second litter was produced by breeding a GRMD male to Speedy, a GRMD/Mstn +/− carrier from the first litter. A total of 10 dogs were studied. Buccal swabs from all dogs of the two litters were analyzed (DDC Veterinary, Fairfield, OH) to demonstrate the same Mstn gene deletion at nucleotides 939 and 940 described previously .
Prior to muscle biopsy and phenotypic tests, dogs were premedicated with acepromazine maleate (0.02 mg/kg), butorphanol (0.4 mg/kg), and atropine sulfate (0.04 mg/kg), masked, and then intubated and maintained with sevoflurane. Samples from the cranial sartorius (CS), vastus lateralis (VL), long digital extensor (LDE), and lateral head of the gastrocnemius (LHG) muscles were removed surgically via open biopsy at 8–9 months of age, snap frozen in isopentane cooled in liquid nitrogen, and stored at −80 °C.
Total cellular ribonucleic acid (RNA) was isolated from frozen skeletal muscle with Tripure reagent (Roche, Indianapolis, IN, USA) and DNase treated with deoxyribonucleic acid (DNA)-free kit (Applied Biosystems, Foster City, CA, USA). The RNA concentrations of the individual samples were measured using a Nanodrop 2000 spectrophotometer and assessed for quality using a 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA).
The canine myostatin gene was sequenced using primers that flank the 939–940 bp mutation site: F:GTGCTGTCGTTACCCTCTAA/R: GAGACATCTTTGTGGGAGTACAG (840–1040 bp). The 200 bp PCR product was cloned into plasmid cloning vector PCR2.1 with TA Cloning Kit (Invitrogen, Carlsbad, CA, USA), plasmid DNA was prepared using QIAprep Spin Miniprep Kit (Qiagen, Hilden, Germany), and the plasmid DNA was sequenced with M13 forward primer. Multiple colonies from each dog were selected, and the sequenced plasmids were assessed to demonstrate the 939–940 bp deletion in heterozygous dogs.
Myostatin gene expression in the CS and the VL from these 10 dogs was reported previously . Gene expression in these muscles plus the LDE and LHG was assessed for this current study. Quantitative realtime RT-PCR (qRT-PCR) primers were designed using Primer Express 3.0 software (Applied Biosystems, Foster City, CA, USA) for the normalizer gene hypoxanthine phosphoribosyltransferase 1 (HPRT1). The recommended ABI TaqMan Gene Expression Assay primer pair and probe were purchased for the canine Mstn gene (Applied Biosystems; catalog number cf02704228_m1). Samples of skeletal muscle RNA (100 ng) were reverse transcribed with oligo-dT and random octamer primers and Superscript II (Invitrogen). qRT-PCR was performed in duplicate reactions with Taqman Gene Expression PCR Master Mix on a 7900HT Fast Real-Time PCR System (Applied Biosystems).
Pelvic limb joint angles were measured at maximal flexion and extension, with additional calculation of the range of motion (ROM), using the method of Jaegger et al. . Our group had previously measured tibiotarsal joint (TTJ) angle in GRMD dogs with a different method in the context of serial peroneus longus force measurements [40, 41]. We continue to make this measurement to compare results to historical values. For this original TTJ angle measurement [40, 41], dogs were placed in dorsal recumbence, with the stifle held at a 90° angle and the tibia/fibula positioned parallel to the table. The angle formed by the flexor surface of the tarsus was measured using a goniometer centered over the lateral malleolus of the fibula. Values for normal dogs using this method approximate but are somewhat less than those recorded at maximal extension.
Tibiotarsal joint force and ECD
Force and eccentric contraction decrement (ECD) were assessed as previously described [38, 42, 43]. Briefly, TTJ flexion and extension torque was measured using a rapid-response servomotor/force transducer (model 310B LR, Aurora Scientific, Inc., Aurora, Ontario, Canada) controlled by a personal computer using custom LabView software. Supramaximal 150 V, 100-μs pulses were applied percutaneously (Model S48 Solid State Square Wave Stimulator, Grass Instruments, Quincy, MA, USA) in a 1.5-s tetanic run of 75 pulses (50/s) to the peroneal (flexion) and tibial (extension) nerves. The site of contact for the paw with the lever (moment arm) was estimated to be 75 % of the distance between the point of the hock and the distal digit. Torque (Newton (N)-meters (m)) was divided by the moment arm (m) to convert to force (N).
Eccentric contractions were induced by percutaneously stimulating the peroneal nerve using square wave pulses of 100-μs duration in a tetanic run for 1 s at a frequency of 50 Hz while simultaneously extending the TTJ with a servomotor (Aurora Scientific) [38, 43]. Thus, the muscles of the cranial tibial compartment were repeatedly stretched to induce mechanical damage. Three sets of 10 stretches for a total of 30, each set separated by 4 min, were performed. Contraction-induced injury was quantified by the force (torque) deficit (Fd) using the following equation: Fd = (Maximal isometric tetanic force [Po] before stretch − Po after stretch/Po before stretch) × 100.
Cranial sartorius circumference
We have previously shown that the CS muscle undergoes dramatic hypertrophy in GRMD dogs and that this hypertrophy tracks with postural abnormalities . Accordingly, as previously described , we used CS circumference measured at surgery during routine biopsy as a surrogate for muscle hypertrophy and associated postural changes in GRMD.
Studies were done on either a Siemens 3T Allegra Head-Only MRI scanner using a circular polarization (CP) head coil or on a whole body MAGNETOM Trio with Tim system with a 32-channel body coil (composed of a 16 element anterior part and 16 element posterior part) at the UNC-CH Biomedical Research Imaging Center [45, 46]. Dogs were placed on an MRI gantry in the sternal (prone) position with the pelvic limbs extended and scanned using a published protocol . T2-weighted image sequences with fat saturation (T2fs) and without fat saturation (T2w) were acquired using a variable-flip-angle turbo spin echo (TSE) sequence. A multi-spin-echo T2 (MSE-T2), with a ten-echo Carr-Purcell-Meiboom-Gill sequence, was acquired to calculate T2 mapping [45, 47]. The T2 mapping was calculated on a voxel-by-voxel basis for the proximal pelvic limbs of dogs in the transverse plane of MSE-T2 images by fitting an exponential decay curve to the signal intensity of the corresponding voxels using a linear-least-squares curve-fitting algorithm .
The major proximal pelvic limb muscles in T2w were registered and then manually segmented, while T2fs was used to identify the contour of each muscle. All proximal pelvic limb muscles were segmented but only five slices at the mid-femur were analyzed/quantified . The T2fs and T2 mapping were subsequently automatically aligned to the coordinate space of T2w using a rigid image registration technique . Thus, the muscle segmentation in T2w could be used in T2 mapping for the biomarker analysis. Texture analysis features were assessed as potential markers of patchy lesions such as necrosis [49, 50].
Femur length was used to correct absolute muscle volumes as previously described . To determine overall scores for T2 and the texture features in each group, the proportional muscle volume was considered, so as to calculate a weighted average.
Muscle histochemical and immunoblot studies
Morphometric analysis of muscle was performed in a blinded fashion. Muscle fibers were aligned for true cross sections (10 μm), stained with hematoxylin and eosin (H&E), and viewed with a light microscope (Axiocam System, Zeiss; GERMANY) at high power. Numbers per 100 myofibers were counted for degenerating fibers (hyaline and/or fragmenting fibers with/without evidence of myophagocytosis), regenerating fibers (small fibers with basophilic nuclei and prominent nuclei/nucleoli), and centrally nucleated fibers (fibers with nuclei located away from the sarcolemma within the cytoplasm; CNF). Muscle fiber CSAs were determined using Scio Image Software (National Institutes of Health (NIH)). A total of 200–600 fibers were measured for each muscle of each dog.
Immunoblots of cell lysates and tissue homogenates were performed as previously described . Hydroxyproline (HP) content of various muscles was measured using a modification of Woessner and colleagues  as described . Values are expressed as μg HP/mg of total muscle protein. For activin receptor type IIB (ActRIIB) Western blots, frozen muscle samples were homogenized in Tissue Extraction Reagent I (Invitrogen), followed by quantification with the BCA Protein Assay Kit (ThermoFisher Scientific, Grand Island, NY, USA). Five micrograms (μg) of total protein from each muscle were loaded onto Novex 10 % Tris-Glycine gels and transferred to polyvinylidene difluoride (PVDF) membranes. Antibodies used were ActRIIB (Abcam ab76940; Cambridge, MA, USA), at a concentration of 1:500, and secondary enhanced chemilumescent (ECL)-anti-mouse IgG, 1:5,000 (GE Healthcare NA931V; Pittsburgh, PA, USA). Following visualization with the Amersham ECL Prime Western Blotting Detection Reagent (GE Healthcare), quantification was performed using the NIH ImageJ program.
Myostatin (GDF-8) levels from the tissue lysates were determined using an immunoaffinity liquid chromatography-mass spectrometry (LC-MS) workflow. Tissue lysates were prepared by homogenization of each muscle sample using tissue protein extraction reagent (TPER, Pierce cat 78510) (ThermoFisher Scientific). Homogenate was centrifuged and 20 μL of the supernatant was used for myostatin analysis. Quantification of myostatin was performed against recombinantly produced myostatin protein (R&D Systems, Inc. Minneapolis, MN, USA; catalog number 788-G8). Recombinant protein standards were prepared in 1 % bovine serum albumin (BSA) in phosphate buffered saline containing 0.05 % Tween-20 (PBST). A seven-point standard curve was prepared from 0.1 to 12.8 ng/mL by serial dilution. Tissue levels were back calculated against this standard curve and normalized using tissue weight; levels are reported in pg of myostatin per mg of tissue. Samples were prepared by first adding 200 μL of 0.1 M Glycine buffer (pH 2.5) to 20 μL of sample and incubating for 1 h at room temperature. After 1 h, 50 μL of 1 M tris buffer (pH 8) was added to each sample, followed by 1 μg of biotin conjugated anti-GDF-8 antibody and incubated overnight at 4 °C. After incubation, streptavidin magnetic beads were used to isolate the myostatin/antibody complex. The magnetic beads were washed twice with PBST, followed by one PBS wash. Myostatin was then separated from the antibody and beads using 140 μL of 30 mM hydrochloric acid prepared in water. The acid eluate was then neutralized using 30 μL of 1 M tris buffer (pH 8). All samples were subsequently reduced using dithiothreitol (DTT), cysteines were alkylated with iodoacetamide, and the protein was digested using trypsin. Quantification occurred by nano liquid chromatography-mass spectrometry and liquid chromatography (LC-MS-MS) using an API 5500Q mass spectrometer against the tryptic peptide IPAMVVDR.
Satellite cell quantitation
For satellite cell quantification of muscle, Tabitha, Hagatha, Endora, and Esmerelda were anesthetized at 37 months of age and samples of the CS and VL were removed surgically and processed as described above. Muscle cryosections (10 μm) were immunostained using antibodies to Pax7 (Developmental Studies Hybridoma Bank (DSHB), Iowa City, IA; 1:10) followed by goat anti-mouse secondary antibody (Alexa Fluor 488, Invitrogen A11001; 1:500) and an antibody to laminin (2E8, DSHB, 1:100), and then goat anti-rabbit secondary antibody (Alexa Fluor 594, Invitrogen, A11012; 1:500). All nuclei were labeled with 4′,6-diamidino-2-phenylindole (DAPI) in mounting medium (P36930, Invitrogen). Image capture was performed using an upright microscope for bright field, differential interference contrast (DIC), and epifluorescence with an ApoTome structured illumination digital imaging system (Carl Zeiss Meditec, Inc, Dublin, CA). The image stacks consisted of six optical sections with 1.4 mm Z-steps. Nuclei staining positive for DAPI but outside laminin extracellular membranes were considered interstitial cells, while those inside laminin were considered myonuclei. Satellite cells were identified as Paired Box 7 (Pax7)+ nuclei within laminin staining as shown in Additional file 1: Figure S1 and as previously described .
The statistical package, R version 3.0.1 , was used for all statistical analyses. Analysis of variance (ANOVA) tests were used to compare results from individual dystrophic groups and controls. From the ANOVA tests, we examined whether at least one group was significantly different from the others. Once the null hypothesis was rejected in the ANOVA, we performed a post hoc Tukey’s honest significant differences (HSD) test to identify differences between a pair of groups. The model can be written as y ij = μ i + Ε ij , where the subscript i represents the group (1 = GRippets; 2 = GRMD dogs; 3 = non-dystrophic controls), the subscript j represents the jth subject, the y ij is the response variable of the jth subject of the ith group, the μ i is the mean of the ith group, and the error term Ε ij is assumed to have a standard normal distribution. This model was repeated to assess messenger ribonucleic acid (mRNA) and protein data, plus each MRI marker, texture feature, functional measurement, and histopathological lesion of each muscle. We computed Pearson correlation coefficients between MRI, histopathological, and mRNA/protein data. To account for multiple comparisons, we applied the false discovery rate (FDR) method to calculate the corrected p values . A value of p < 0.05 was considered significant. Trends were cited if p < 0.2.
Key molecular and phenotypic features of the dogs of this study are discussed below. Data are expressed as mean ± SD. Because of the small group sizes and stringency of the FDR statistical method, many of the differences between the two dystrophic groups did not reach significance. Significant (p < 0.05) findings and those that trended towards significance (p < 0.2) between GRMD and GRippet dogs are summarized in Table 2; some are illustrated in Fig. 1.
Posture, gait, and outcome
Objective measures of gait, such as accelerometry and video gait analysis, were not evaluated in these dogs. However, key postural changes and quantitative joint angle measurements (below) were compared between the two dystrophic groups. Overall, the dystrophic dogs had a more plantigrade posture and gait, with the limbs shifted forward under the torso. A plantigrade posture results from the carpi being hyperextended (dorsal flexion) and the tarsi being hyperflexed (palmar flexion). Postural abnormalities were more pronounced in the GRippets, equating to a more severe GRMD phenotype. These changes varied between the two litters. Dash (GRippet) from the first litter had a pronounced plantigrade stance, and his pelvic limbs were shifted forward (Fig. 2). He was no longer ambulatory at 9 months of age and was euthanized. Similar but milder postural changes were seen in the three GRippets of the second litter (Derwood, Abner, and Tabitha) (Fig. 3). While these dogs retained the ability to walk, they developed glossal hypertrophy, which interfered with their ability to eat (Fig. 4). Resulting weight loss necessitated euthanasia of Abner and Derwood at 9 and 11 months, respectively. The final GRippet, Tabitha, remained ambulatory but had pronounced postural abnormalities and was euthanized at 41 months.
Compared to control dogs, joint angles of GRMD dogs were similar, whereas those of GRippets differed. The GRippet hip angle with maximal flexion was larger (more restricted) (77.5° ± 4.93), and hip ROM was smaller (56.5 ± 12.5) versus controls (58.3° ± 5.51, p < 0.01 for flexion; 87.7 ± 15.9, p < 0.05 for ROM) (Fig. 1a). An unexpected and important finding was that GRippets had larger (more restricted) maximal TTJ (53.5° ± 7.85), stifle (42.0 ± 5.35), and hip (77.5° ± 4.93) (Fig. 1a) flexion angles compared to GRMD dogs (35.0 ± 7.00, 34.0 ± 3.61, and 63.0 ± 3.61, respectively) (p < 0.05 for TTJ and hip; p = 0.173 for stifle). The TTJ ROM motion trended towards being significantly decreased in the GRippets (p = 0.052) (Table 3).
Force and ECD
As with the joint angle data, GRippet body-mass-corrected force values (N/kg) varied markedly between the two litters. Dash from the first litter had dramatically higher TTJ tetanic flexion force (12.5) compared to the three GRMD dogs from the second litter (6.63 ± 1.59). A comparable but reversed differential was seen for extensor force, with Dash having a value of 11.1 versus 28.8 ± 2.75 for the other three GRippets. Flexion force values for GRippet (0.70 ± 0.41) and GRMD dogs (0.59 ± 0.32) were each lower than the controls (1.54 ± 0.04) (p < 0.05 for both). The GRMD (1.93 ± 0.41) and GRippet (1.92 ± 0.64) dogs also had lower extension force compared to controls (3.16 ± 0.48) (p < 0.05 for GRippets and p = 0.062 for GRMD). Values for ECD after either 10 or 30 contractions were higher in each dystrophic group versus controls (p < 0.05 for all). Neither force nor ECD values differed significantly between the two dystrophic groups (Table 4).
Cranial sartorius circumference
Cranial sartorius circumference corrected for body mass (mm/kg) trended towards being higher in the GRippets (5.05 ± 1.72; p = 0.053) but not GRMD dogs (3.39 ± 0.52, p = 0.617) when compared to controls (2.46 ± 0.25). Values for the two dystrophic groups did not differ (p = 0.219).
Effects of muscle atrophy or hypertrophy in dystrophic dogs were demonstrated by MRI, with parallel changes being seen with both absolute (mm3) and femur length-corrected (mm3/mm) volumes and the percent contribution that each muscle made to the overall CSA (Fig. 5 and Additional file 1: Figure S1; Additional file 2: Table S1).
Muscle volumes varied considerably among the dogs and were not clearly associated with Mstn gene status. As an example, overall mid-femur length-corrected muscle volume in the transverse plane ranged from 114 in Hagatha (GRMD) to 191 in Endora (GRMD carrier/Mstn +/+). Racer, the only dog that was normal for both genetic traits, had an intermediate value of 164. Esmeralda (GRMD normal/Mstn +/–) had a higher value (173), potentially in keeping with a hypertrophic effect of myostatin loss, but Endora’s value was even higher (191).
A differential effect of the dystrophy trait on individual muscles was demonstrated by assessing volumes of the quadriceps femoris (QF) and the combined cranial and caudal sartorius (total sartorius) muscles. The corrected QF volume ranged from 19.2 in Dash (GRippet) to 47.0 for Racer, while an opposite pattern was seen for the total sartorius, with values ranging from 4.3 for Racer to 13.0 for Dash.
GRippet and GRMD dogs versus non-dystrophic controls
When data from individual dystrophic groups were compared, controls had greater absolute (29,631 ± 274), and femur length-corrected (178 ± 13.4) volumes than either GRMD (19,682 ± 1255, absolute; 125 ± 8.87, corrected) or GRippet (21,791 ± 1409, absolute; 137 ± 7.76, corrected) dystrophic dogs (p < 0.01 for all). Differential effects on individual muscles paralleled those of the individual dogs. Control QF volumes were larger than those of either dystrophic group (p < 0.01 for both). Total sartorius volumes of the controls were smaller than those for either GRippet or GRMD dogs, with the GRippet value trending towards significance (p = 0.159 for absolute and 0.112 for corrected). In comparing either absolute or corrected volumes of the other muscles of the two dystrophic groups separately to controls, only the biceps femoris differed, being significantly smaller in both the GRMD and GRippet dogs (p < 0.01 for most).
With regard to CSA, the QF accounted for 25.4 ± 2.99 % in the controls, versus 21.6 ± 0.9 % in the GRMD dogs and 15.9 ± 0.61 % in the GRippets (Fig. 1b). Values for the GRippets differed significantly from controls (p < 0.01) but those for the GRMD dogs did not (p = 0.864). The CS/VL ratio showed a clear progression from the controls (0.24 ± 0.06), to the GRMD dogs (0.48 ± 0.04), to the GRippets (0.86 ± 0.19) (Fig. 1c), indicating that reduction in myostatin accentuated unequal muscle growth in dystrophic dogs. Differences for GRippets versus controls were significant (p < 0.05), but those for GRMD dogs were not (p = 0.864).
GRippet versus GRMD dogs
While GRippet absolute and femur length-corrected total muscle volumes were similar to those of the GRMD dogs, the pattern of differential muscle involvement between dystrophic and control dogs was exaggerated in the GRippets. For almost all muscles, the degree of atrophy or hypertrophy in GRMD dogs was more pronounced in the GRippets (Additional file 1: Figure S1). The differential effects on the quadriceps and sartorius muscles were highlighted by the CS/VL ratio, which trended towards being higher in the GRippets (p = 0.168).
Individual muscle volume on MRI and joint angles of the 10 dogs were correlated to clarify whether imbalanced muscle volume/strength contributed to postural changes. We expected that CS muscle volume would correlate negatively with TTJ angle measured by our original method  but neither the muscle volumes nor the percent CSA values correlated. This outcome may have been influenced by the age of the dogs. With our prior studies, the correlation was noted at 6 months of age; correlations have not been demonstrated at other ages. On the other hand, CS and total sartorius absolute volumes and percent CSA correlated positively with maximal hip flexion angle (r = 0.671 to 0.756; p < 0.05 for all but CS absolute volume for which p = 0.063). There was also a strong negative correlation between absolute and femur length-corrected QF volume (r = −0.829 and 0.809) and percent CSA (r = −0.935) versus maximal hip flexion (p < 0.05 for absolute and corrected volumes and 0.01 for percent CSA). These correlations could be in keeping with biomechanical relationships or represent independent variables tracking with the disease phenotype.
T2 mapping values and texture analysis
The overall T2 mapping value for the non-dystrophic control dogs (38.2 ± 5.00) was lower than that for the GRippets (52.3 ± 3.47) and GRMD dogs (51.9 ± 5.09) (p < 0.05 for both). T2 mapping values of individual muscles, with the exception of the sartorius heads, were also higher, although differences generally only trended towards significance. Values for texture analysis features did not differ significantly among the three groups (data not shown) (Additional file 3: Table S2).
Pathologic and molecular findings
Histopathologic changes (degenerating and regenerating fibers; centrally nucleated fibers)
Averaged scores for degenerating fibers and CNF were higher in GRippet and GRMD dogs versus controls (p < 0.01 for all but CNF in GRMD which was <0.05). In contrast, the difference for regenerating fiber scores in the two dystrophic groups only trended towards significance in GRMD dogs (p = 0.101) and was not significant in the GRippets (p = 0.523). Most scores for degenerating and regenerating fibers for individual muscles of separate dystrophic groups did not differ from control individual muscles. The one exception was the presence of greater numbers of degenerating fibers in the LDE of the GRMD dogs (p < 0.05). In contrast, CNF were increased in both dystrophic groups compared to controls in all muscles (p < 0.05 for most). Importantly, histopathologic lesions in GRippet and GRMD muscle did not differ. Specifically, there were no statistically significant differences in degenerating, regenerating, or CNF in GRippets compared to their GRMD littermates (Additional file 4: Table S3).
To investigate the potential that modulation of myostatin in muscular dystrophy might lead to exhaustion of the satellite cell pool, we quantified Pax7-positive nuclei in two dystrophic dogs, including one GRippet (Tabitha) and another Mstn +/+ GRMD dog (Hagatha), as well as two GRMD carriers that differed on myostatin status (Endora (Mstn +/+) and Esmerelda (Mstn +/−)). All dogs were adults, 37 months of age. There were no statistically significant differences in satellite cell numbers between muscles from dogs of different myostatin genotypes. However, somewhat unexpectedly, Pax7-positive cells accounted for a higher percentage of myonuclei in the two dystrophic dogs (4.3–5.7 %) versus the two GRMD carriers (0.6–2.8 %) for the four muscles studied (Table 5; Additional file 5: Figure S2).
CSA of myofibers
Consistent with the increase in muscle fiber size in the genetic absence or postnatal blockade of myostatin [6, 14], the CSA for the GRippet CS muscle (4133 ± 233) was larger than that of their GRMD littermates (2661 ± 316) (p < 0.01). In comparing the CSA of the CS fibers to controls (2,598 ± 346), GRippets had larger fibers (p < 0.01), while those of the GRMD dogs did not differ (p = 0.964). These findings of muscle fiber hypertrophy were not seen in other muscles. In fact, in the VL, the CSA was smaller in both dystrophic groups (GRippets, 1560 ± 155, p < 0.01; GRMD, 1850 ± 610; p < 0.05), compared to the controls (3,196 ± 326) (Fig. 1d) (Additional file 4: Table S3).
Hydroxyproline (HP) protein
Fibrosis is a major complicating factor of chronic myopathy. Studies have shown that myostatin stimulates muscle fibroblasts and that fibrosis is reduced in the absence of myostatin . Hydroxyproline content (μg HP/mg of total muscle protein) is a biomarker of muscle fibrosis and collagen content. Thus, we were interested in the relative HP content of the dystrophic and control groups. Overall HP content of GRMD dogs (26.5 ± 9.30) was higher than controls (9.53 ± 2.56) (p < 0.05), while that of the GRippet dogs (22.9 ± 4.79) trended towards being higher (p = 0.053). Although HP content for individual muscles from the dystrophic groups did not differ significantly when compared to controls, there were trends towards greater content in the GRippet VL (p = 0.074) and GRMD VL, CS, and LHG (all p < 0.2). Interestingly, the differential between the two dystrophic groups varied markedly among muscles, with GRippet values being lower in the CS and LHG but higher in the VL and LDE. None of these differences were significant (Additional file 4: Table S3).
Myostatin mRNA and protein expression
Myostatin mRNA data from the CS and VL of these same dogs have been reported previously . Levels in the CS of both dystrophic groups were significantly lower than controls, while those in the VL were lower but not significantly different. Gene expression in these muscles plus the LDE and LHG was assessed for this current study. Levels averaged for all four muscles were reduced to a comparable degree in GRippet and GRMD dogs at 8–9 months when compared to controls but differences were not significant. In assessing individual muscles, levels were reduced in the CS of GRMD (0.093 ± 0.059) and GRippet (0.041 ± 0.026) dogs, compared to controls (1.08 ± 0.458) (p < 0.01 for both) (Fig. 6). A less pronounced reduction was seen in the VL of GRMD (0.306 ± 0.20) and GRippet (0.335 ± 0.185) dogs versus controls (1.05 ± 0.388) (p < 0.05 for both) (Additional file 4: Table S3).
Similar to mRNA findings, average myostatin protein levels (pg/mg) for the three muscles assessed were lower in GRippet (2.94 ± 1.48) and GRMD (4.30 ± 0.98) dogs compared to controls (7.63 ± 3.51) but differences were not significant, nor did values differ between dystrophic groups (p = 0.770). Levels were significantly lower in the dystrophic versus control CS regardless of Mstn genotype (0.92 ± 0.73 in the GRippets and 1.63 ± 0.39 in GRMD) vs. 7.79 ± 2.18 in the two non-dystrophic dogs) (p < 0.05 for both). This is in keeping with the dramatic hypertrophy seen in the CS in even Mstn +/+ GRMD dogs [44, 56] and suggests that the hypertrophy may be more pronounced when myostatin is further down regulated. Consistent with this observation, Dash, the GRippet with the most pronounced postural changes and CS muscle hypertrophy, also had the lowest myostatin protein levels (0.21) by a wide margin. Values for the LDE and LGH did not differ between the dystrophic and control dogs.
Myostatin mRNA fold change and protein levels were correlated with MRI volumes and histopathologic changes among the dogs (all 10 for mRNA and 8 for protein) to help establish potential cause and effect links. Notably, there was a significant negative correlation between CS myostatin mRNA fold change (r = −0.7297) and protein (r = −0.7649) levels and CS percent CSA on MRI (both p < 0.05), suggesting myostatin down regulation may contribute to CS hypertrophy. A similar but positive correlation was seen between Mstn mRNA levels and VL volume corrected for femur length (r = 0.680; p < 0.05), providing a further potential link between down regulation of myostatin and VL atrophy.
With regard to histopathological findings, both myostatin mRNA and protein levels in the CS correlated negatively with the levels of degenerating fibers and CNF (r = −0.8178, p < 0.01 for mRNA; r = −0.7617, p < 0.05 for protein). A less pronounced effect was seen in the VL, with mRNA levels correlating negatively with degenerating fibers (r = −0.6433, p < 0.05) and a trend occurring for CNF (r = −0.5538, p = 0.097). These data indicate that the exaggerated lowering of myostatin could place muscle fibers at increased risk for injury (degenerating fibers) and subsequent regeneration (CNF).
Activin receptor type IIB (ActRIIB)
Expression of the putative receptor for myostatin, ActRIIB, varies among muscles in mice. For example, the fast-twitch EDL has higher levels than the slow twitch soleus . Thus, we were interested in whether the relative expression of ActRIIB protein among muscles of dogs in this study correlated with differential muscle size. ActRIIB protein expression normalized to beta-actin varied markedly among the four muscles. For Racer, the only completely normal dog, levels in the CS (0.077 arbitrary units) were considerably lower than those in the VL (0.299), LDE (0.408), and LHG (0.259). In contrast, for Endora (GRMD carrier, Mstn +/+), the CS ActRIIB level (0.190) was higher than that of either the LDE (0.090) or LHG (0.011) but lower than the VL (0.391). Thus, based on these limited data from the two non-dystrophic, Mstn +/+ dogs, there was no clear pattern of differential ActRIIB expression among the four muscles to account for the pattern of muscle atrophy/hypertrophy protein (Additional file 4: Table S3 and Fig. 7).
ActRIIB levels have not been assessed previously in GRMD dogs. Interestingly, levels for all muscles of the combined dystrophic dogs (0.98 ± 0.42) were higher compared to average overall values for Racer and Endora, the two non-dystrophic dogs with wild-type myostatin (0.22 ± 0.06) (p < 0.01). Moreover, ActRIIB levels in individual muscles were higher in dystrophic dogs (CS, p < 0.01; VL and LDE, p < 0.05) in all but the LHG (p = 0.585). However, neither overall mean (GRippet 0.92 ± 0.52; GRMD 1.06 ± 0.32) (p = 0.903) nor individual muscle ActRIIB expression differed between the Mstn +/+ and Mstn +/−dystrophic groups.
The level of feedback between myostatin and ActRIIB, or other growth factors and receptors, could differ among muscles, leading to variable atrophy or hypertrophy. We correlated levels of ActRIIB protein with myostatin mRNA and protein to further establish a potential relationship. We were particularly interested in the CS muscle, given its marked hypertrophy. For all dogs taken together, levels of CS ActRIIB trended to correlate negatively with those of myostatin mRNA (r = −0.726; p = 0.108) and protein (r = −0.6739; p = 0.142). This would be in keeping with an expected physiologic relationship between the two, with ActRIIB levels increasing as myostatin levels decline. The degree of correlation between myostatin mRNA and ActRIIB levels was assessed in the dystrophic groups alone to gain further insight into potential disease associations. No correlations were seen.
Development of myostatin inhibitors for muscle-wasting disorders is based on the assumption that enhanced muscle regeneration will improve muscle function. Because mdx mice, which lack myostatin [13, 14] or in which myostatin is blocked  are less severely affected, we hypothesized that the relative genetic loss of myostatin would improve signs in dystrophic dogs. Interpretation of data in this study is compromised by the small group sizes and the inherent phenotypic variation seen in GRMD dogs. With that said, phenotypic features of GRMD dogs heterozygous for myostatin (GRippets) and their GRMD littermates were similar. Thus, we were unable to demonstrate a beneficial effect of reduced myostatin in GRMD dogs. In fact, when compared to their GRMD dystrophic littermates, GRippets actually had more severe postural changes, evidenced by an exaggerated plantigrade stance, shifting of the pelvic limbs forward, and more restricted hip joint flexion. Strikingly, the relative absence of myostatin in the GRippets exaggerated pre-existing trends for either muscle atrophy or hypertrophy on MRI, most notably in the quadriceps and sartorius muscles, as shown by the progressive increase of the CS to VL volume ratio of the non-dystrophic controls (0.24 ± 0.06), to the GRMD dogs (0.48 ± 0.04), to the GRippets (0.86 ± 0.10). This differential pattern of muscle atrophy/hypertrophy appeared to have functional consequences, in that CS and VL volumes correlated, in opposite directions, with maximal hip flexion angle.
We have previously demonstrated unequal disease effects on TTJ flexors and extensors in natural history and preclinical studies in GRMD dogs [36, 42] and shown that proportional weakness of TTJ extensors or increased strength of flexors exaggerates TTJ joint flexor contractures . Moreover, the degree of CS hypertrophy in GRMD correlates negatively with the TTJ joint angle , suggesting that this muscle could play a role analogous to iliotibial band tightening in DMD . Differential muscle involvement also plays a role in the DMD disease course. Unbalanced forces acting on a joint, due to either disproportionate muscle weakening or shortening, lead to abnormal positioning. Several studies have assessed the proportional strength of joint agonist and antagonist muscles in DMD, with most concluding that muscle imbalance contributes to contractures [58–62]. Those noting a relationship identified a strong negative correlation between extensor muscle weakness and flexor contracture severity in DMD. As opposing extensor muscles weakened, flexor contractures worsened. To our knowledge, such contractures and postural changes are not seen in bully whippet dogs (Mstn −/− ) despite their gross muscle enlargement. Instead, a single bully whippet included in the original study was judged to run faster than myostatin wild-type dogs . This suggests that contractures are only precipitated when myostatin down regulation is coupled with the dystrophic state.
Key questions arising from our study relate to (1) why relative loss of myostatin would cause differential muscle effects that exaggerate the GRMD phenotype and (2) whether this untoward effect could potentially extend to postnatal treatment of humans with muscle-wasting disorders. Treatments directed at inhibiting or blocking myostatin are driven by a belief that myostatin is innately limiting the muscle regenerative response. Importantly, in considering potential therapies for DMD, one must take into account homeostatic compensatory mechanisms that are already in play. In cases of muscle loss, the body should naturally suppress myostatin expression to promote muscle regeneration. Indeed, microarray gene expression profiling in DMD , the mdx mouse , and GRMD dog  has shown that myostatin is down regulated in dystrophic muscle independent of modulating treatments. Consistent with findings reported here, the degree of myostatin mRNA down regulation in GRMD dogs also has been shown previously to vary among muscles, with levels being dramatically reduced in the CS, a hip flexor prone to hypertrophy (44, 56), and less so or not at all in selected extensor muscles . Not surprisingly, myostatin mRNA and protein levels in the CS of the GRippet dogs of our study were also markedly reduced when measured at 8–9 months of age. Interestingly, at this relatively late age, the GRippet levels were not significantly lower than those in their myostatin wild-type GRMD littermates. Had levels been measured earlier, particularly in advance of the rapid period of disease progression beginning at 3 months , we may have been able to distinguish a difference between GRMD and GRippet dogs. With this said, correlations between myostatin mRNA and protein levels with phenotypic features such as CS muscle volume on MRI and the degree of degenerating and regenerating fibers in both the CS and VL still suggested that the myostatin expression could contribute to variable muscle effects and the overall disease phenotype
The differential response of murine muscles to loss of myostatin is associated with varying levels of the myostatin receptor, ActRIIB . In comparing ActRIIB levels among the three groups of dogs in our study, dystrophic dogs taken together had significantly higher levels than non-dystrophic controls, suggesting a feedback mechanism whereby ActRIIB levels increase secondary to lowering of myostatin. Similarly, when all dogs were assessed collectively, myostatin mRNA and protein levels tended to correlate negatively with ActRIIB expression in the CS muscle and, to a lesser extent, the VL. This would be in keeping with an expected physiologic relationship between the two, with ActRIIB levels increasing as myostatin levels decline. On the other hand, when the dystrophic groups were assessed alone, no correlation was seen. This likely is simply a function of small group sizes and the stringency of the FDR statistical method but could point to dysregulation in dystrophic muscle.
Myostatin favors differentiation towards slow versus fast glycolytic fibers, due apparently to positive and negative regulation of MEF2 and MyoD, respectively . Consistent with this effect, cattle lacking Mstn  and Mstn-null mice  have increased numbers of fast glycolytic fibers. Therefore, in principle, fiber type distribution could influence effects of myostatin loss or blockade. Muscles with a predominantly fast glycolytic fiber type distribution could be more sensitive to the absence of myostatin and display greater hypertrophy. While most murine muscles respond similarly to genetic depletion of myostatin , the slow twitch soleus undergoes less pronounced hypertrophy than the fast-twitch EDL [20, 22, 24]. Based on the fast/slow fiber type distribution of canine muscles, the CS (49:51), VL (57:43), and LHG (50:50) are mixed muscles, while the LDE (71:29) is predominantly fast-twitch . Independent of myostatin status, the GRMD LDE undergoes dramatic atrophy, the VL and LHG are moderately atrophied, and the CS is markedly hypertrophied [44, 56]. Therefore, fiber type distribution, alone, does not appear to be a major factor in GRMD muscle atrophy/hypertrophy.
With regard to the second question posed above on the potential for our findings to translate to humans, it is not clear whether deleterious effects of genetic loss of myostatin will extend to postnatal treatments. During development, myostatin down regulates myogenic regulators such as MyoD to inhibit proliferation and differentiation of myoblasts [5, 69]. Genetic loss of myostatin in the mouse accelerates both primary (embryo) and secondary (fetal) myogenesis, with an increase in myofiber numbers (hyperplasia) and size (hypertrophy) . Myostatin-null mice have attained 87 % of their final adult number of myofibers by embryonic day 18.5 (near term) compared to only 73 % of wild-type mice . Further postnatal muscle enlargement relies more on an increase in protein synthesis rather than incorporation of more nuclei . Thus, one would logically suspect that effects of genetic versus postnatal loss or inhibition of myostatin would differ. However, the postural abnormalities seen in the GRippet dogs of our study did not become apparent until around 3 months of age. Accordingly, while the basis for differential muscle involvement may have been established in the embryo, phenotypic effects were not seen until after birth. Further, a previous study investigated the effect of postnatal stimulation of muscle growth in four GRMD dogs for 4 months beginning at 2.5 months of age. Results of this study were presented at the 16th International Congress of the World Muscle Society in 2011  but have not been published. This study used treatment with anti-ActRIIB antibody, an inhibitor of multiple BMP ligands including myostatin, to increase muscle growth. Overall muscle mass did not increase in treated dogs compared to GRMD controls. In keeping with our GRippet findings, several indices of motor function declined in treated dogs (Blot S, personal communication, 2015). While these findings are consistent with those of the GRippets, they may not have been specific to the effects of myostatin inhibition because a non-selective inhibitor was used.
We have previously shown that myostatin blockade with an AAV8-myostatin propeptide construct administered via regional limb delivery in 3-month-old normal dogs leads to muscle enlargement . As with the GRippets reported here, the effects were not uniform among muscles. In particular, the treated VL showed no increase in size, even though transgene levels were comparable to muscles that showed a treatment effect. In a separate study, a group of GRMD dogs were treated systemically with AAV8-myostatin dominant negative propeptide at 10 months of age and followed for 13 months with MRI and at necropsy to determine effects on distal pelvic limb muscle size . All treated GRMD muscles increased in size but the degree varied, ranging from 27 % in the LDE to 49 % in the cranial tibialis. Notably, none of the muscles in either of these studies decreased in size. This contrasts with our findings in the GRIppets, in which the relative absence of myostatin exaggerated pre-existing trends for either muscle atrophy or hypertrophy. Exaggeration of muscle atrophy, as with the quadriceps, is particularly puzzling and could have considerable potential clinical significance.
Treatments directed at inhibiting myostatin to increase the regenerative response of diseased muscle have been postulated to potentially exhaust satellite cell replicative capacity . Previous studies with dystrophic mouse models and genetic deletion or postnatal inhibition of myostatin are not adequate to address the question of satellite cell senescence since mice have considerably longer telomeres  and proportionally higher telomerase activity than humans . Like humans, dogs have relatively short telomeres and low telomerase activity [74, 75]. The GRippet is thus a good model to evaluate satellite cell populations in dystrophic muscle with chronic reduction in myostatin. Despite active degeneration and regeneration in GRippet muscles, satellite cells were the same proportion (approximately 5–6 %) of myonuclei as in GRMD muscles and did not show evidence of exhaustion.
Myostatin mRNA and protein levels were down regulated in GRMD dogs independent of myostatin status, suggesting an inherent feedback mechanism intended to promote muscle regeneration. The relative absence of myostatin in the GRippets exaggerated pre-existing trends for either muscle atrophy or hypertrophy on MRI. This differential pattern of muscle atrophy/hypertrophy appeared to have functional consequences, in that phenotypic features of the GRippets were, if anything, more severe than GRMD dogs with wild-type myostatin. Results from these GRippet dogs suggest that dystrophic muscles may be differentially affected by in utero myostatin loss. While data derived from dystrophic dogs lacking myostatin from inception will not necessarily extrapolate to postnatal treatments, this study reinforces the complexity of factors regulating muscle function and further emphasizes that murine models of muscular dystrophy may behave very differently than large animals such as dogs and humans. Although the lack of satellite cell exhaustion in GRippets is reassuring for future human experience with myostatin inhibitors, the muscle imbalance and severity of contractures requires continued careful evaluation. These findings may serve to inform future clinical trials of postnatal myostatin inhibition.
activin receptor type IIB
analysis of variance
bovine serum albumin
centrally nucleated fibers
differential interference contrast
Duchenne muscular dystrophy
Developmental Studies Hybridoma Bank
laminin alpha2-deficient mouse
- Dyf−/− :
mouse model of LGMD2B and Myoshi myopathy
eccentric contraction deficit
extensor digitorum longus
false discovery rate
growth and differentiation factor 8
golden retriever muscular dystrophy
hematoxylin and eosin
hypoxanthine phosphoribosyltransferase 1
honest significant differences
liquid chromatography-mass spectrometry
liquid chromatography-mass spectrometry and Liquid chromatography-tandem mass spectrometry
lateral head of gastrocnemius
magnetic resonance imaging
messenger ribonucleic acid
National Institutes of Health
paired box 7
phosphate buffered saline containing 0.05 % Tween-20
polymerase chain reaction
maximal isometric tetanic force
quantitative real-time polymerase chain reaction
range of motion
T2-weighted image sequences with fat saturation
T2-weighted image sequences without fat saturation
tissue protein extraction reagent
turbo spin echo
University of North Carolina-Chapel Hill
Hoffman EP, Brown RH, Kunkel LM. Dystrophin: the protein product of the Duchenne muscular dystrophy locus. Cell. 1987;51:919–28.
Grounds MD. Two-tiered hypotheses for Duchenne muscular dystrophy. Cell Mol Life Sci. 2008;65:1621–5.
Manzur AY, Muntoni F. Diagnosis and new treatments in muscular dystrophies. J Neurol Neurosurg Psy. 2009;80:706–14.
Mendell JR, Shilling C, Leslie ND, Flanigan KM, al-Dahhak R, Gastier-Foster J, et al. Evidence-based path to newborn screening for Duchenne muscular dystrophy. Ann Neurol. 2012;71:304–13.
Lee S-J. Regulation of muscle mass by myostatin. Annu Rev Cell Dev Biol. 2004;20:61–86.
McPherron AC, Lawler AM, Lee S-J. Regulation of skeletal muscle mass in mice by a new TGF-β superfamily member. Nature. 1997;387:83–90.
Lee S-J, McPherron AC. Myostatin and the control of skeletal muscle mass. Curr Opin Genet Dev. 1999;9:604–7.
McPherron AC, Lee S-J. Suppression of body fat accumulation in myostatin-deficient mice. J Clin Invest. 2002;109:595–601.
Schuelke M, Wagner KR, Stolz LE, Hubner C, Riebel T, Komen W, et al. Myostatin mutation associated with gross muscle hypertrophy in a child. N Engl J Med. 2007;350:2682–8.
McPherron AC, Lee S-J. Double muscling in cattle due to mutations in the myostatin gene. Proc Natl Acad Sci USA. 1997;94:12457–61.
Clop A, Marcq F, Takeda H, Pirottin D, Tordoir X, Bibé B, et al. A mutation creating a potential illegitimate microRNA target site in the myostatin gene affects muscularity in sheep. Nat Genet. 2006;38:813–8.
Mosher DS, Quignon P, Bustamante CD, Sutter NB, Mellersh CS, Parker HG, et al. A mutation in the myostatin gene increases muscle mass and enhances racing performance in heterozygote dogs. PLoS Genet. 2007;3:779–86.
Wagner KR, McPherron AC, Winik N, Lee S-J. Loss of myostatin attenuates severity of muscular dystrophy in mdx mice. Ann Neurol. 2002;52:832–6.
Wagner KR, Liu X, Chang X, Allen RE. Muscle regeneration in the prolonged absence of myostatin. Proc Natl Acad Sci USA. 2005;102:2519–24.
Bogdanovich S, Krag TO, Barton ER, Morris LD, Whittemore LA, Ahima RS, et al. Functional improvement of dystrophic muscle by myostatin blockade. Nature. 2002;420:418–21.
Bradley L, Yaworsky PJ, Walsh FS. Myostatin as a therapeutic target for musculoskeletal disease. Cell Mol Life Sci. 2008;65:2119–24.
Patel K, Amthor H. The function of myostatin and strategies of myostatin blockade—new hope for therapies aimed at promoting growth of skeletal muscle. Neuromuscul Disord. 2005;15:117–26.
Mendias CL, Bakhurin KI, Faulkner JA. Tendons of myostatin-deficient mice are small, brittle, and hypocellular. Proc Natl Acad Sci U S A. 2008;105:388–93.
Amthor H, Macharia R, Navarrete R, Schuelke M, Brown SC, Otto A, et al. Lack of myostatin results in excessive muscle growth but impaired force generation. Proc Natl Acad Sci U S A. 2007;104:1835–40.
Mendias CL, Marcin JE, Calerdon DR, Faulkner JA. Contractile properties of EDL and soleus muscles of myostatin-deficient mice. J Appl Physiol (1985). 2006;101:898–905.
Qaisar R, Renaud G, Morine K, Barton ER, Sweeney HL, Larsson L. Is functional hypertrophy and specific force coupled with the addition of myonuclei at the single muscle fiber level? FASEB J. 2012;26:1077–85.
Li ZF, Shelton GD, Engvall E. Elimination of myostatin does not combat muscular dystrophy in dy mice but increases postnatal lethality. Am J Pathol. 2005;166:491–7.
Lee YS, Lehar A, Sebald S, Liu M, Swaggart KA, Talbot Jr CC, et al. Muscle hypertrophy induced by myostatin inhibition accelerates degeneration in dysferlinopathy. Hum Mol Genet. 2015;24:5711–9.
Parsons SA, Millay DP, Sargent MA, McNally EM, Molkentin JD. Age-dependent effect of myostatin blockade on disease severity in a murine model of limb-girdle muscular dystrophy. Am J Pathol. 2006;168:1975–85.
Capecchi MR. Gene targeting in mice: functional analysis of the mammalian genome for the twenty-first century. Nat Rev Genet. 2005;6:507–12.
Sacca R, Engle SJ, Qin W, Stock JL, McNeish JD. Genetically engineered mouse models in drug discovery research. Methods Mol Biol. 2010;602:37–54.
Lin JH. Applications and limitations of genetically modified mouse models in drug discovery and development. Curr Drug Metabol. 2008;9:419–38.
Schneider MR, Wolf E, Braun J, Kolb HJ, Adler H. Canine embryo-derived stem cells and models for human diseases. Hum Mol Genet. 2007;17(R1):R42–47.
Tsai KL, Clark LA, Murphy KE. Understanding hereditary diseases using the dog and human as companion model systems. Mamm Genome. 2007;18:444–51.
Lindblad-Toh K, Wade CM, Mikkelsen TS, Karlsson EK, Jaffe DB, Kamal M, et al. Genome sequence, comparative analysis and haplotype structure of the domestic dog. Nature. 2005;438:803–19.
Qiao C, Li J, Zheng H, Bogan J, Li J, Yuan Z, et al. Hydrodynamic limb vein injection of AAV8 canine myostatin propeptide gene in normal dogs enhances muscle growth. Hum Gene Ther. 2009;20:1–10.
Bish LT, Sleeper MM, Forbes SC, Morine KJ, Reynolds C, Singletary GE, et al. Long-term systemic myostatin inhibition via liver-targeted gene transfer in golden retriever muscular dystrophy. Hum Gene Ther. 2011;12:1499–509.
Kornegay JN, Tuler SM, Miller DM, Levesque DC. Muscular dystrophy in a litter of golden retriever dogs. Muscle Nerve. 1988;11:1056–64.
Cooper BJ, Winand WN, Stedman H, Valentine BA, Hoffman EP, Kunkel LM, et al. The homologue of the Duchenne locus is defective in X-linked muscular dystrophy of dogs. Nature. 1988;334:154–6.
Sharp NJH, Kornegay JN, Van Camp SD, Herbstreith MH, Secore SL, Kettle S, et al. An error in dystrophin mRNA processing in golden retriever muscular dystrophy, an animal homologue of Duchenne muscular dystrophy. Genomics. 1991;13:115–21.
Kornegay JN, Bogan JR, Bogan DJ, Childers MK, Li J, Nghiem P, et al. Canine models of Duchenne muscular dystrophy and their use in therapeutic strategies. Mamm Genome. 2012;23:85–108.
Cotten SW, Kornegay JN, Bogan DJ, Wadosky KM, Patterson C, Willis MS. Genetic myostatin decrease in the golden retriever muscular dystrophy model does not significantly affect the ubiquitin proteasome system despite enhancing the severity of disease. Am J Transl Res. 2013;6:43–53.
Kornegay JN, Bogan JR, Bogan DJ, Childers MK, Grange RW. Golden retriever muscular dystrophy (GRMD): developing and maintaining a colony and physiological functional measurements. In: Duan D, editor. Muscle gene therapy: methods and protocols, Methods in Molecular Biology, vol. 709. New York: Humana Press; 2011. p. 105–23.
Jaegger G, Marcellin-Little DJ, Levine D. Reliability of goniometry in Labrador retrievers. Am J Vet Res. 2002;63:979–86.
Kornegay JN, Sharp NJH, Schueler RO, Betts CW. Tarsal joint contracture in dogs with golden retriever muscular dystrophy. Lab Anim Sci. 1994;44:331–3.
Kornegay JN, Sharp NJH, Bogan DJ, Van Camp SD, Metcalf JR, Schueler RO. Contraction tension and kinetics of the peroneus longus muscle in golden retriever muscular dystrophy. J Neurol Sci. 1994;123:100–7.
Kornegay JN, Bogan DJ, Bogan JR, Childers MK, Cundiff DD, Petroski GF, et al. Contraction force generated by tarsal joint flexion and extension in dogs with golden retriever muscular dystrophy. J Neurol Sci. 1999;166:115–21.
Childers MK, Grange RW, Kornegay JN. In vivo canine muscle function assay. J Vis Exp. 2011, Apr 5;(50). doi: 10.3791/2623.
Kornegay JN, Cundiff DD, Bogan DJ, Bogan JR, Okamura CS. The cranial sartorius muscle undergoes true hypertrophy in dogs with golden retriever muscular dystrophy. Neuromuscul Disord. 2003;13:493–500.
Wang J, Fan Z, Shiloh-Malawsky Y, An H, Kornegay JN, Styner M. A computerized MRI biomarker quantification scheme for a canine model of Duchenne muscular dystrophy. Int J Comput Assist Radiol Surg. 2013;8:763–74.
Fan Z, Wang J, Ahn M, Shiloh-Malawsky Y, Chahin N, Elmore S, et al. Characteristics of magnetic resonance imaging biomarkers in a natural history study of golden retriever muscular dystrophy. Neuromuscul Disord. 2014;24:178–91.
Hennig J, Speck O, Scheffler K. Optimization of signal behavior in the transition to driven equilibrium in steady-state free precession sequences. Magn Reson Med. 2002;48:801–9.
Wang J, Fan Z, Kornegay J, Styner M. MRI-based quantification of Duchenne muscular dystrophy in a canine model. Proc SPIE. 2011;8672:36720F.
Mahmoud-Ghoneim D, Cherel Y, Lemaire L, de Certaines JD, Maniere A. Texture analysis of magnetic resonance images of rat muscles during atrophy and regeneration. Magn Reson Imaging. 2006;24:167–71.
Wang J, Fan Z, Vandenborne K, Walter G, Shiloh-Malawsky Y, An H, et al. Statistical texture analysis based MRI quantification of Duchenne muscular dystrophy in a canine model. Proc SPIE. 2013;7965:79650G.
Li ZB, Kollias HD, Wagner KR. Myostatin directly regulates skeletal muscle fibrosis. J Biol Chem. 2008;283:19371–8.
Woessner Jr JF. The determination of hydroxyproline in tissue and protein samples containing small proportions of this imino acid. Arch Biochem Biophys. 1961;93:440–7.
Zhang Y, King OD, Rahimov F, Jones TI, Ward CW, Kerr JP, Liu N, Emerson CP Jr, Kunkel LM, Partridge TA, Wagner KR. Human skeletal muscle xenograft as a new preclinical model for muscle disorders. Hum Mol Genet. 2014;15:3180–8.
R Development Core Team: R. A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2012; ISBN 3-900051-07-0, URL http://www.R-project.org/.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Roy Stat Soc, Series B. 1995;57:289–300.
Nghiem PP, Hoffman EP, Mittal P, Brown K, Schatzberg SJ, Ghimbovschi S, et al. Sparing of the dystrophin-deficient cranial sartorius muscle is associated with classical and novel hypertrophy pathways in GRMD dogs. Am J Pathol. 2013;183:1411–24.
Rideau Y. Treatment of orthopedic deformity during the ambulatory stage of Duchenne muscular dystrophy. In: Serratrice G, Cros D, Desnuelle C, editors. Neuromuscular Disease. New York: Raven; 1984. p. 557–64.
Vignos Jr PJ, Spencer Jr GE, Archibald JC. Management of muscular dystrophy of childhood. JAMA. 1963;184:89–96.
Archibald KC, Vignos Jr PJ. A study of contractures in muscular dystrophy. Arch Phys Med Rehabil. 1959;40:150–7.
Siegel IM. Equinocavovarus in muscular dystrophy. Its treatment by percutaneous tarsal medullostomy and soft tissue release. Arch Surg. 1972;104:644–6.
Wagner MB, Vignos Jr PJ, Carlozzi C. Duchenne muscular dystrophy: a study of wrist and hand function. Muscle Nerve. 1989;12:236–44.
McDonald CM, Abresch RT, Carter GT, Fowler Jr WM, Johnson ER, Kilmer DD, et al. Profiles of neuromuscular diseases: Duchenne muscular dystrophy. Am J Phys Med Rehabil. 1995;74:S70–92.
Chen Y-W, Nagaraju K, Bakay M, McIntyre O, Rawat R, Shi R, et al. Early onset of inflammation and later involvement of TGFβ in Duchenne muscular dystrophy. Neurology. 2005;65:826–34.
Tseng BS, Zhao P, Pattison JS, Gordon SE, Granchelli JA, Madsen RW, et al. Regenerated mdx mouse skeletal muscle shows differential mRNA expression. J Appl Physiol. 2002;93:537–45.
Hennebry A, Berry C, Siriett V, O’Callaghan P, Chau L, Watson T, et al. Myostatin regulates fiber-type composition of skeletal muscle by regulating MEF2 and MyoD gene expression. Am J Physiol Cell Physiol. 2009;296:C525–34.
Stavaux D, Art T, McEntee K, Reznick M, Lekeux P. Muscle fibre type and size, and muscle capillary density in young double-muscled blue Belgian cattle. Zentralbl Veterinarmed A. 1994;41:229–36.
Girgenrath S, Song K, Whittemore LA. Loss of myostatin expression alters fiber-type distribution and expression of myosin heavy chain isoforms in slow- and fast-type skeletal muscle. Muscle Nerve. 2005;31:34–40.
Armstrong RB, Saubert IV CW, Seeherman HJ, Taylor CR. Distribution of fiber types in locomotory muscles of dogs. Am J Anat. 1982;163:87–98.
Manceau M, Gros J, Savage K, Thomé V, McPherron A, Paterson B, et al. Myostatin promotes the terminal differentiation of embryonic muscle progenitors. Genes Dev. 2008;22:668–81.
Matsakas A, Otto A, Elashry MI, Brown SC, Patel K. Altered primary and secondary myogenesis in the myostatin-null mouse. Rejuvenation Res. 2010;13:717–27.
Mouisel E, Barthélémy I, Ferry A, Garcia L, Blot S, Amthor H. The role of activin receptor IIB signaling on skeletal muscle and the possible therapeutic implication for Duchenne muscular dystrophy. Neuromuscul Disord. 2011;21:711.
Kipling D, Cooke HJ. Hypervariable ultra-long telomeres in mice. Nature. 1990;347:400–2.
Prowse KR, Greider CW. Developmental and tissue-specific regulation of mouse telomerase and telomere length. Proc Natl Acad Sci U S A. 1995;92:4818–22.
Argyle DJ, Nasir L. Telomerase: a potential diagnostic and therapeutic tool in canine oncology. Vet Pathol. 2003;40:1–7.
Cross SH, Allshire RC, McKay SJ, McGill NI, Cooke HJ. Cloning of human telomeres by complementation in yeast. Nature. 1989;338:771–4.
Kirberger RW, Steenkamp G, Spotswood TC, Boy SC, Miller DB, van Zyl M. Stenotic nasopharyngeal dysgenesis in the Dachshund: seven cases (2002-2004). J Am Anim Assoc. 2006;42:290–7.
These studies were supported, in part, by the Co-operative Program in Translational Research: Proposal for Establishment of the National Center for Canine Models of Duchenne Muscular Dystrophy (NCDMD) (1U24NS059696-01A1; NINDS/NIAMS) (J.N.K.).
The authors declare that they have no competing interests.
JNK conceived the study together with KRW, performed phenotypic studies, analyzed and interpreted the data, and wrote the manuscript. DJB, JRB, and JLD bred the dogs, provided animal care, contributed to phenotypic studies, and processed samples for analysis. JW and ZF analyzed MRI data. NL and LCW prepared muscle histochemical, immunoblot, and histopathological sections. RWG contributed to analysis of force and eccentric contraction decrement data. MA and HZ conducted the statistical analysis. CJB and CBL conducted confirmatory myostatin qRT-PCR and sequencing studies; SWC and MSW performed the initial myostatin qRT-PCR. JP and CAM quantified GDF8 protein. MAS devised the methods for MRI analysis and reviewed all MRI data. KRW conceived the study with JNK, analyzed the histochemical, immunoblot, and histopathological sections, and critically reviewed the manuscript. All authors read and approved the final manuscript.
Figure S1. Histograms depicting percent cross-sectional area contributed by each muscle on MRI at midthigh of control, GRMD, and GRippet dogs. Pre-existing atrophy or hypertrophy in GRMD muscles is generally more exaggerated in the GRippets. (TIF 664 kb)
Table S1. MRI volumes in non-dystrophic control, GRMD (Mstn +/+), and GRippet (Mstn +/−) Dogs (Mean ± SD). (DOCX 954 kb)
Table S2. T2 mapping values (Mean ± SD) in non-dystrophic control, GRMD (Mstn +/+), and GRippet (Mstn +/−) Dogs (mean ± SD). (DOCX 369 kb)
Table S3. Histopathologic, morphometric, and molecular findings in non-dystrophic control, GRMD (Mstn +/+), and GRippet (Mstn +/−) dogs. (DOCX 403 kb)
Figure S2. Pax7-positive cells in Mstn +/+ and Mstn +/− GRMD/GRippet and GRMD-carrier dogs. Laminin (A), Pax7 (B), DAPI (C), and merged (D) staining are seen. Satellite cells were defined as Pax7+ nuclei within the laminin + basal lamina (also see data in Additional file 4: Table S3). (TIF 915 kb)