Magnetic resonance imaging of the erector spinae muscles in Duchenne muscular dystrophy: implication for scoliotic deformities
© Zoabli et al; licensee BioMed Central Ltd. 2008
Received: 12 September 2008
Accepted: 29 December 2008
Published: 29 December 2008
In Duchenne muscular dystrophy (DMD), the muscular degeneration often leads to the development of scoliosis. Our objective was to investigate how anatomical changes in back muscles can lead to scoliosis. Muscular volume and the level of fat infiltration in those muscles were thus evaluated, in non-scoliotic, pre-scoliotic and scoliotic patients. The overlying skin thickness over the apex level of scoliotic deformations was also measured to facilitate the interpretation of electromyographic signals when recorded on the skin surface.
In 8 DMD patients and two healthy controls with no known muscular deficiencies, magnetic resonance imaging (MRI) was used to measure continuously at 3 mm intervals the distribution of the erector spinae (ES) muscle in the T8-L4 region as well as fat infiltration in the muscle and overlying skin thickness: four patients were non-scoliotic (NS), two were pre-scoliotic (PS, Cobb angle < 15°) and two were scoliotic (S, Cobb angle ≥ 15°). For each subject, 63 images 3 mm thick of the ES muscle were obtained in the T8-L4 region on both sides of the spine. The pixel dimension was 0.39 × 0.39 mm. With a commercial software, on each 12 bits image, the ES contour on the left and on the right sides of the spine were manually determined as well as those of its constituents i.e., the iliocostalis (IL), the longissimus (LO) and the spinalis (SP) muscles. Following this segmentation, the surfaces within the contours were determined, the muscles volume were obtained, the amount of fat infiltration inside each muscle was evaluated and the overlying skin thickness measured.
The volume of the ES muscle of our S and PS patients was found smaller on the convex side relative to the concave one by 5.3 ± 0.7% and 2.8 ± 0.2% respectively. For the 4 NS patients, the volume difference of this muscle between right and left sides was 2.1 ± 1.5% and for the 2 controls, it was 1.4 ± 1.2%. Fat infiltration for the S and the PS patients was larger on the convex side than on the concave one (4.4 ± 1.6% and 4.5 ± 0.7% respectively) and the difference was more important near the apex. Infiltration was more important in the lateral IL muscle than in the medial SP and it was always larger near L2 than at any other spinal level. Fat infiltration was much more important in the ES for the DMD patients (49.9% ± 1.6%) than for the two controls (2.6 ± 0.8%). As for the overlying skin thickness measured near the deformity of the patients, it was larger on the concave than on the convex side: 14.8 ± 6.1 vs 13.5 ± 5.7 mm for the S and 10.3 ± 6.3 vs 9.8 ± 5.6 mm for the PS.
In DMD patients, our results indicate that a larger replacement of muscles fibers by fat infiltration on one side of the spine is a factor that can lead to the development of scoliosis. Efforts to slow such an infiltration on the most affected side of the spine could thus be beneficial to those patients by delaying the apparition of the scoliotic deformation. In addition to anatomical considerations, results obtained from the same patients but in experiments dealing with electromyography recordings, point to differences in the muscular contraction mechanisms and/or of the neural input to back muscles. This is similar to the adolescent idiopathic scoliosis (AIS) where a role of the nervous system in the development of the deformation has also been suggested.
Duchenne muscular dystrophy (DMD) is a severe genetic disease affecting one in 3500 boys . Due to lack of dystrophin, muscle fibers are susceptible to mechanical damage leading to their replacement with fat tissue . As muscular weakness progresses, most patients will experience falls resulting in disabling extremity fractures  and some of them will also develop a scoliosis . While many researches are related to the genetic aspects of the disease  and to muscle damage repair , investigations have also been carried on the identification of anatomical changes associated with the disease.
For example, it is possible with ultrasound (US) imaging to monitor DMD progression by observing the pathologic changes of skeletal muscle fibers, the level of fat infiltration and the effect of a therapy . With this imaging modality, the volume of superficial muscles can be estimated both at rest and under contraction). To assess fat infiltration in DMD back muscles, computerized tomography (CT) can also be used . For scoliotic DMD patients, Stern and Clark  reported that fat infiltration was more important on the concave side of their deformity. These authors found that the difference was proportional to the severity of the deformity and was correlated with the progression of the DMD within a 6-month period. This was not the case for the non-scoliotic patients where the amount of muscle fibers replaced by fat was found similar on both sides of the spine. Since CT imaging modality implies radiation to the patient, magnetic resonance imaging (MRI) is getting more often used . When body composition of DMD patients was evaluated with MRI , the mean fat mass represented 31.6 ± 17.0% of the body weight which is significantly higher than the 15.5 ± 12.1% obtained from growth charts. While MRI is capable of producing reliable results, processing of the collected images was lengthy few years ago but this has changed appreciably. For instances, with automatic segmentation methods, subcutaneous fat could be separated from muscle tissue within 2 s with an error < 3% .
Using MRI, the objective in the present study was to measure, along the spine, the volume and the level of fat infiltration in back muscles of DMD patients in order to better understand how the dystrophy affects those muscles. The overlying skin thickness over those muscles was also studied since it affects the characteristics of the surface electromyographic (EMG) signals recorded over those muscles.
Materials and methods
Information on the controls (C1 and C2) and on the DMD patients (P1 to P8) classified according to their Cobb angle.
Cobb angle (°)
Apex level, Convexity
DFZ duration (months)
The paraspinal muscles were imaged over the T8 to L4 region (9 cm above and 9 cm below T12-L1) where the apex of the scoliotic curvature in DMD patients is usually located. A 1.5 Tesla Symphony MRI Scanner (Siemens Medical Solutions, Malvern, PA, United States) was used with a 3D SE-T1 sequence (TR 550 ms, TE 14 ms, Nex 3) . Slice thickness was 3 mm, the pixel dimension was 0.39 × 0.39 mm and the 12 bits image size was 512 × 512 pixels. A set of 63 images were collected with a 7 slices interleave to reduce the effect of the breathing and motion during MRI acquisition. The acquisition time, including subject preparation, ranged from 18 to 24 min. When a set of images had been acquired, each image within the set was rapidly inspected to detect if blurring occurred due to motion artefacts of the subject since no sedation was used. Following such verification, a new series of images had to be acquired for two patients.
The manual segmentation procedure used for back muscles was initially developed for the upper arm where segmentation in axial, then sagittal, and coronal planes was developed to minimize error in the detection of muscles boundaries. Validity of this procedure was assessed with two observers performing independently the segmentation of the same images and using two different image processing software. To assess the accuracy of the segmentation procedure, a MRI plexiglass phantom (cylinder of ∅ = 187 mm, 60 mm thick with holes diameters from 11.5 down to 1.1 mm) was used as well as MRI images obtained from the upper limb of normal subjects and from the erector spinae (ES) of adolescent scoliotic patients. The accuracy on a cross section area determination varied with the size of the object: for the phantom, it was 0.2% for its external diameter of 187.0 mm, 2.9% for a hole of 11.0 mm and 4.9% for 7.0 mm. A look-up table of the estimated volume versus the diameter of each hole was obtained. The segmentation accuracy was then estimated to be = 1% for back muscles (∅~30 mm) and = 0.9% for biceps brachii (∅~50 mm). Accuracy was also found quite dependent on image contrast than on the human factors involved in the segmentation . In the present study with DMD patients, the data acquisition parameters offering the best contrast where chosen and the images were segmented by the observer who developed the technique (G.Z.). With a marching cube algorithm [17, 18]. a faceted volumic model of the ES muscle and of the infiltrated fat was obtained. With the look-up table, the accuracy in the determination of the volume of the ES and its 3 constituents was estimated to be ≤ 1%.
Following a manual delimitation of the outline of the left and right erector spinae (ES) and IL, LO and SP muscles with SliceOmatic software (Tomovision Inc., Montreal), a threshold was experimentally set at 512 over 4096 gray levels within each closed contour for the automatic detection of the infiltrated fat on each slice. For each muscle, the cross section area of 6 to 8 consecutives slices, approximately the thickness of a vertebra, were used to get a volume associated to each vertebral level. As for the amount of fat infiltration in a muscle, it was expressed in percentage of that muscle's volume. For the controls and the NS patients, % values on the right were averaged as well as those on the left. Similarly, for the PS and S patients, an average for the convex side and one for the concave side was obtained. No distinction was made whether the spine deformation was on the right or on the left
The overlying skin thickness was measured on both sides of the spine from T9 to T11 over IL, LO and SP muscles. The cutaneous folds caused by creases in the sheet over the MRI platform that could be associated to muscular material, eliminated by averaging 9 measures made on each side of the folds where contrast between fat and muscle was low.
Due to a reduced number of patients, statistical analysis only included mean, standard deviation and linear correlation calculations.
For the patients, the volume of the IL muscle ranged from 30 to 70 cm3 and it was larger for the controls (Fig. 2B). For the NS patients, a larger volume of IL (7.5 ± 3.7%) was found either on the right or on the left. For the PS and S patients, the IL volume was larger on the concave side respectively by 4.8 ± 1.9% and by 7.8 ± 4.0%. The volume of the LO muscle was slightly larger than the IL (Fig. 2C) and for the NS patients, a larger IL volume of 3.4 ± 1.0% was generally found on the left. As for the PS and S patients, volume of this muscle was larger on the concave side by 5.3 ± 3.9% and 5.7 ± 2.1% respectively. As for the SP muscle (Fig. 2D), the difference in volume was 3.1 ± 0.9% for the NS. A larger volume of SP was found on the concave side of one PS patient and the two S patients but the difference (1.8 ± 1.3%) was similar between the two groups. Volume difference between patients and the controls was less important for the SP than for the IL and the LO muscles.
The overlying skin thickness was measured at T10 on each side of the spine. From the 9 measures taken on each side, the NS patients had a similar thickness on both sides of the spine (an average of 11.3 ± 4.0 mm for the right side and 11.2 ± 4.2 mm for the left side). As for the PS and S patients, thickness was larger on the concave than on the convex side (10.3 ± 6.3 vs 9.8 ± 5.6 mm and 14.8 ± 6.1 vs 13.5 ± 5.7 mm respectively).
Using MRI to study the back muscles of DMD patients, fat infiltration was found larger on the convex side of the deviations of pre- and scoliotic patients while no appreciable difference was observed in non-scoliotic patients. A different result was obtained by Stern and Clark  when they studied the back muscles of 16 scoliotic DMD patients with CT scans. They reported that fat infiltration was more important on the concave side of the spinal deviation but when their density readings (Table 2 of their paper) are closely scrutinized, one can be lead to consider that fat infiltration was larger on the concave side for only half of their patients. Also, their density readings were taken only at T9 and L3 levels in the medial (i.e. SP) and lateral (i.e. IL) portions of the ES and fat infiltration in the LO was probably not taken into account. In contrast, our data were continuously collected at every 3 mm from T8 to L4 and all the constituents of the ES muscle were analyzed. Through the reliability offered by the use of our 3D segmentation procedure, the accuracy of volume determination can be considered similar to the one of the cross section area (i.e. ≤ 1%). For the ES muscle of the two controls, the smallest differences in muscle volume was 1.4 ± 1.2% and for fat infiltration it was 2.6 ± 0.8% while for the DMD patients, these differences were larger. All of our results are then outside the segmentation accuracy. Segmentation was easier to achieve for ES due to its size and for SP because fat infiltration was minimal. For more severely infiltrated LO and IL muscles, an additional segmentation cycle (axial, sagittal, then coronal plane) was done compared to the number of iterations used to delineate ES and SP.
Some of the differences between our results and those of Stern and Clark may be associated to the effects of the DFZ: all of our patients had been and some were still treated (n = 6) with DFZ while none of their patients received such a treatment. DFZ has an overall positive impact on DMD patients' quality of life: it improves cardiac function, prolongs walking, and seems to eliminate the need for spinal surgery, although vertebral fractures and stunted growth occur . Weight gain is a side effect of the medication but no correlation was found here between BMI and the duration of DFZ medication.
The presence of an uneven fat infiltration can induce a force unbalance on the spine which could lead to its buckling toward the weaker side. As the spine buckles, the skin becomes compressed on the concave side and stretched on the convex one. Our results indicate that the thickness difference was smaller for the PS than for the S patients (and absent for the NS patients). Due to the filtering effects of the overlying skin thickness, smaller EMG signals would be expected on the concave side. However, in another set of data collection with the same patients , a larger surface EMG signals were recorded on the concave side of the PS patients and on the convex side of a S patient. While smaller signals on the convex side of the PS may have result from the greater fat infiltration on that side, it cannot explain the larger signal on the convex of the S patient. It can thus be hypothesized that a more important neural input on one side of the spine can contribute to some extent to the presence of larger EMG signals on that side . In our pre-scoliotic DMD patients, a deficit in muscle fibers on one side of the spine and an imbalance in the neural input could thus be two factors leading to the deformation of the spine.
As the Cobb angle gets larger than 15°, it seems that compensatory mechanisms to prevent further deviation of the spine are set in action. For instance, for scoliotic DMD P7, a larger EMG activity was detected on his convex side  as it is frequently reported in AIS [22, 23]. The stretching of the skin on the convex side of a deviation causes a reduction in the overlying skin thickness, favoring the presence of a larger EMG signal but a greater fat infiltration on the convex side would counteract to some level such an increase. The presence of a larger EMG signal on the convex side may thus be attributed to a greater feedback from the more stretched muscle spindles and/or a larger neural input from the central nervous system. Imbalanced neural input was also considered to have a pathogenic importance in the etiology of AIS .
For Shimada , the trunk muscle imbalance would be one of the most important factors in the onset and progression of AIS since EMG amplitude is higher on the convex side when the deviation is progressive while no differences are observed in non progressive situation. According to the medical records of our PS and S patients, their spinal deformation did not progress over an 18 months period preceding our images acquisition. Either DFZ halted the progression of the deformation or their scoliosis was non progressive. Since higher EMG signals were still detected on the convex side, it seems that the scoliotic mechanism present in a DMD is somewhat different than in AIS.
To predict the progression of scoliosis in DMD patients, approaches such as Griffiths developmental scales, the Reynell language scales and the British picture vocabulary scales have been proposed [24, 25]. but their reliability seems to be limited as compared to the vital capacity at age 10 when ambulation usually ceases . As for Gibson and Wilkins , they proposed a clinical mobility index based on muscle weakness which is essentially caused by the progressive replacement of muscle tissue by fat. In our study, that mobility index was found correlated with fat infiltration in the ES (R2 = 0.70, p < 0.01) but not with BMI nor with the age of the patients. This lack of correlation with age may be specific to the ES and its constituents since for thigh muscles (hip, mid-thigh, and knee), older boys demonstrated more prominent fatty infiltration than younger ones . However, it has to be considered that despite a severe infiltration, NS P4, which had a very limited mobility, did not develop scoliosis since infiltration was similar on both sides of his spine. Therefore, an asymmetry in fat infiltration (Fig. 4) could help predict the development of a spinal deformation.
Only the ES muscle and its constituents were studied here since they can greatly influence the development of a scoliosis and their activity can be easily monitored from the skin surface. Investigation of less accessible trunk muscles, such as the abdominal and external oblique, psoas and quadratus lumborum, could be worth of investigation as they can also influence the stability of the spine.
DMD patients are frequently solicited to participate in a research protocol and we experienced difficulties with patients' recruitment. The small number of patients participating in our protocol constitutes the main limitation of the study. With fat infiltration appearing initially in the L2 region, images acquisition could be restricted in that zone thus reducing the time required for the execution of the MRI protocol. This could contribute to facilitate the recruitment of patients.
In our pre- and scoliotic DMD patients, the volume of muscle fibers of the erector spinae was found smaller and fat infiltration was more important on the convex than on the concave side. Anatomically, the presence of less muscle fibers on one side of the spine due to a larger fat infiltration is a factor that can contribute to the application of uneven forces on the spine leading to its deformation. Physiologically, uneven forces could also result from the action of the central nervous system on those muscles and this action seems to be different between the pre-scoliotic and scoliotic stages. Since infiltration along the spine took place initially in the IL muscle around the L2 level and spread upward and downward towards the LO and the SP, any effort to halt such a progression could help improve the condition of the DMD patients. As for the overlying skin thickness, it was larger on the concave side of the deviation and this may help explain in part why a smaller EMG signal is generally recorded on that side. While the scoliotic mechanism present in a DMD share similarities with AIS, it also has differences. These preliminary results deserve further experimental work for confirmation.
The authors thank Dr. Decarie and the radiology department of Ste-Justine Hospital for their collaboration. Special thanks to Mrs. Diane Choquette and Ginette Labrecque for their help in the programming of the acquisition sequence. Many thanks to Mr. Christian Bellefleur, and Mrs. Julie Joncas, the research nurse responsible for patients' recruitment and to Mrs. Linda Pelletier for reviewing the text. Work supported by the Canadian Institutes of Health Research, the Natural Sciences and Engineering Research Council of Canada (156144-04) and le Fonds Québécois de la Recherche sur la Nature et les Technologies.
- Voisin V, de la Porte S: Therapeutic strategies for Duchenne and Becker dystrophies. Int Rev Cytol. 2004, 240: 1-30.PubMedView Article
- Pedemonte M, Sandri C, Schiaffino S, Minetti C: Early decrease of IIx myosin heavy chain transcripts in Duchenne muscular dystrophy. Biochem Biophys Res Commun. 1999, 255: 466-9.PubMedView Article
- McDonald DG, Kinali M, Gallagher AC, Mercuri E, Muntoni F, Roper H, Jardine P, Jones DH, Pike MG: Fracture prevalence in Duchenne muscular dystrophy. Dev Med Child Neurol. 2002, 44: 695-8.PubMedView Article
- Suzuki S, Kasahara Y, Yamamoto S, Seto Y, Furukawa K, Nishino Y: Three-dimensional spinal deformity in scoliosis associated with cerebral palsy and with progressive muscular dystrophy. Spine. 1993, 18: 2290-4.PubMedView Article
- Becker K, Robb SA, Hatton Z, Yau SC, Abbs S, Roberts RG: Loss of a single amino acid from dystrophin resulting in Duchenne muscular dystrophy with retention of dystrophin protein. Hum Mutat. 2003, 21: 651-PubMedView Article
- Neri M, Torelli S, Brown S, Ugo I, Sabatelli P, Merlini L, Spitali P, Rimessi P, Gualandi F, Sewry C, Ferlini A, Muntoni F: Dystrophin levels as low as 30% are sufficient to avoid muscular dystrophy in the human. Neuromuscul Disord. 2007, 17: 913-8.PubMedView Article
- Hughes MS, Marsh JN, Wallace KD, Donahue TA, Connolly AM, Lanza GM, Wickline SA: Sensitive ultrasonic detection of dystrophic skeletal muscle in patients with duchenne muscular dystrophy using an entropy-based signal receiver. Ultrasound Med Biol. 2007, 33: 1236-43.PubMedPubMed CentralView Article
- Miyatani M, Kanehisa H, Ito M, Kawakami Y, Fukunaga T: The accuracy of volume estimates using ultrasound muscle thickness measurements in different muscle groups. Eur J Appl Physiol. 2004, 91: 264-72.PubMedView Article
- Stern L: New techniques in diagnosis, assessment of progression and research in muscular dystrophy. Aust Paediatr J. 1988, 24 (Suppl 1): 34-6.PubMed
- Stern LM, Clark BE: Investigation of scoliosis in Duchenne dystrophy using computerized tomography. Muscle Nerve. 1988, 11: 775-83.PubMedView Article
- Tingart MJ, Apreleva M, Lehtinen JT, Capell B, Palmer WE, Warner JJ: Magnetic resonance imaging in quantitative analysis of rotator cuff muscle volume. Clin Orthop. 2003, 415: 104-10.PubMedView Article
- Pichiecchio A, Uggetti C, Egitto MG, Berardinelli A, Orcesi S, Gorni KO, Zanardi C, Tagliabue A: Quantitative MR evaluation of body composition in patients with Duchenne muscular dystrophy. Eur Radiol. 2002, 12: 2704-9.PubMed
- Barra V, Boire JY: Segmentation of fat and muscle from MR images of the thigh by a possibilistic clustering algorithm. Comput Methods Programs Biomed. 2002, 68: 185-93.PubMedView Article
- Gibson DA, Wilkins KE: The management of spinal deformities in Duchenne muscular dystrophy. A new concept of spinal bracing. Clin Orthop. 1975, 108: 41-51.PubMedView Article
- Zoabli G, Mathieu PA, Aubin CE: Back muscles biometry in adolescent idiopathic scoliosis. Spine J. 2007, 7: 338-44.PubMedView Article
- Zoabli G, Mathieu PA, Aubin CE, Tinlot A, Beauséjour M, Feipel V, Malanda A: Assessment of manual segmentation of magnetic resonance images of skeletal muscles. IEEE Eng Med Biol Conf. Proceedings of the 23rd Annual International Conference, Istanbul. 2001, 3: 2685-7.
- Lorenson WE, Cline HE: Marching cubes: a high resolution 3D surface construction algorithm. Computer Graphics. 1987, 21: 163-9.View Article
- Jacq JJ, Roux C: Geodesic morphometry with applications to 3-D morpho-functional anatomy. Proc IEEE Emer Med Imag Tech. 2003, 91 (10): 1680-98.
- Houde S, Filiatrault M, Fournier A, Dube J, D'Arcy S, Berube D, Brousseau Y, Lapierre G, Vanasse M: Deflazacort use in duchenne muscular dystrophy: an 8-year follow-up. Pediatr Neurol. 2008, 38: 200-6.PubMedView Article
- Thouin JF, Mathieu PA, Ciolofan OC, Joncas J, Aubin CE: "EMG of back muscles in Duchenne Muscular Dystrophy patients,". International Research Society of Spinal Deformities, Symphosium 2004. 2004, Vancouver, BC, Canada, 293-6.
- Avikainen VJ, Rezasoltani A, Kauhanen HA: Asymmetry of paraspinal EMG-time characteristics in idiopathic scoliosis. J Spinal Disord. 1999, 12: 61-7.PubMed
- Shimada Y: [A study of trunk muscle in idiopathic scoliosis]. Nippon Seikeigeka Gakkai Zasshi. 1989, 63: 33-44.PubMed
- Odermatt D, Mathieu PA, Beausejour M, Labelle H, Aubin CE: Electromyography of scoliotic patients treated with a brace. J Orthop Res. 2003, 21: 931-6.PubMedView Article
- Smith RA, Sibert JR, Harper PS: Early development of boys with Duchenne muscular dystrophy. Dev Med Child Neurol. 1990, 32: 519-27.PubMedView Article
- Hiller LB, Wade CK: Upper extremity functional assessment scales in children with Duchenne muscular dystrophy: a comparison. Arch Phys Med Rehabil. 1992, 73: 527-34.PubMed
- Yamashita T, Kanaya K, Kawaguchi S, Murakami T, Yokogushi K: Prediction of progression of spinal deformity in Duchenne muscular dystrophy: a preliminary report. Spine. 2001, 26: E223-6.PubMedView Article
- Marden FA, Connolly AM, Siegel MJ, Rubin DA: Compositional analysis of muscle in boys with Duchenne muscular dystrophy using MR imaging. Skeletal Radiol. 2005, 34: 140-8.PubMedView Article
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.