Abstract

The purpose of this study was to determine if variables calculated from diffusion tensor imaging (DTI) would serve as a reliable marker of damage after a muscle strain injury in dystrophic (mdx) and wild type (WT) mice. Unilateral injury to the tibialis anterior muscle (TA) was induced in vivo by 10 maximal lengthening contractions. High resolution T1- and T2-weighted structural MRI, including T2 mapping and spin echo DTI was acquired on a 7T small animal MRI system. Injury was confirmed by a significant loss of isometric torque (85% in mdx versus 42% in WT). Greater increases in apparent diffusion coefficient (ADC), axial, and radial diffusivity (AD and RD) of the injured muscle were present in the mdx mice versus controls. These changes were paralleled by decreases in fractional anisotropy (FA). Additionally, T2 was increased in the mdx mice, but the spatial extent of the changes was less than those in the DTI parameters. The data suggest that DTI is an accurate indicator of muscle injury, even at early time points where the MR signal changes are dominated by local edema.

1. Introduction

The muscular dystrophies (MDs) are a heterogeneous group of inherited disorders characterized by progressive weakness and degeneration of skeletal muscles. Duchenne muscular dystrophy (DMD), the most common form of MD, is an X-linked disorder that was first described over a century ago [1]. DMD is caused by the absence of dystrophin, a 427 kDa protein found on the cytoplasmic surface of the plasma membrane of muscle fibers (the sarcolemma) in skeletal and cardiac muscle. It is well established that the absence of dystrophin leads to the impaired linkage between the cytoskeleton and the extracellular matrix protein laminin, which is thought to result in fragility of the sarcolemma and a disruption in the proper transmission of force during contraction.

Muscle strains are one of the most common complaints treated by physicians [2]. When an activated muscle lengthens because the external load exceeds the tension generated by the muscle contraction, this is termed a lengthening (“eccentric”) contraction. The force generated during a maximal lengthening contraction is at least twofold the force developed during a maximal isometric contraction; therefore lengthening contractions are more likely to produce damage than either isometric or concentric contractions [3]. Eccentric contractions are especially harmful to dystrophic muscle [4, 5].

Plain films, or X-rays, are not very useful for imaging muscle pathology, unless heterotopic bone formation has occurred within the muscle. Diagnosis of acute muscle strains is still typically made based on physical examination and patient history, but muscle injuries can be detected with MR imaging methods [6]. As noninvasive technology continues to improve and imaging, such as magnetic resonance imaging (MRI), becomes more commonplace, these tools will play a greater role in diagnosis, prognosis, and in rehabilitation planning [7]. Unlike X-rays, MRI offers superb tissue contrast and has high sensitivity to the hemorrhage and edema that follow muscle injuries. This, together with the capability to evaluate multiple arbitrary anatomic planes, make it the ideal technique to evaluate muscle injures. Because of the similarities between muscle injury and muscle disease, MRI may yield useful information for both of these conditions.

Conventionally, muscle strains are revealed best by T2-weighted MRI images, which optimize contrast between injured muscles with edema (increased signal intensity) and normal uninjured muscles. More recently, diffusion tensor imaging (DTI) has been explored as a more accurate marker for muscle damage compared to T2-weighted MRI [8]. The variables obtained with DTI, at least in other tissues such as the brain [9], show a strong and rapid response to damage, whereas the T2 signal can take a prolonged period to change. DTI is based on measurement of the apparent diffusion of water in tissues. Self-diffusion of water in tissue is restricted by membranes, resulting in an apparent diffusion coefficient (ADC) that is lower than the free diffusion coefficient and is orientation dependent for elongated structures. In short, DTI can be used to determine the three-dimensional architecture of skeletal muscle [10, 11] as well as provide useful information regarding muscle damage [12].

The purpose of this study was to determine if measurements obtained from DTI would serve as a reliable marker of damage after a muscle strain injury in dystrophic (mdx) and control mice. We hypothesize that DTI biomarkers would provide a more informative assessment of muscle injury than T2 and that these measures will further elucidate the increased susceptibility to injury of the mdx model.

2. Materials and Methods

2.1. Injury

All protocols were approved by the University of Maryland Institutional Animal Care & Use Committee (IACUC). The injury model results in a significant and reproducible injury and has been described in detail previously [1316]. Briefly, 5 adult healthy (C57BL/10ScSn) and 5 dystrophic (C57BLScSn-DMDmdx) male mice were purchased (The Jackson Laboratory, Bar Harbor, Me) and were approximately 8 weeks old at the start of experiments. Before each injury protocol, mice were anesthetized (~4-5% isoflurane in an induction chamber, then ~2% isoflurane via a nosecone for maintenance) using a precision vaporizer (cat no. 91103, Vet Equip, Inc, Pleasanton, Calif). Sterile ophthalmic cream (Paralube Vet Ointment, PharmaDerm, Floham Park, NJ) was applied to each eye to protect the corneas from drying. During the procedure, the animal was kept warm by use of a heat lamp.

With the animal supine, the hindlimb was stabilized, and the foot was secured onto a plate, the axis of which was attached to a stepper motor (model T8904, NMB Technologies, Chatsworth, Calif) and a torque sensor (QWFK-8M, Sensotec) to measure torque. The fibular nerve was stimulated via subcutaneous needle electrodes (J05 Needle Electrode Needles, 36BTP, Jari Electrode Supply, Gilroy, Calif), and proper electrode position was determined by a series of isometric twitches. Impulses generated by an S48 square pulse stimulator (Grass Instruments, West Warwick, RI) were 1 ms in duration and passed through a PSIU6 stimulator isolation unit (Grass Instruments, West Warwick, RI). In addition to visual confirmation of isolated dorsiflexion, an increase in twitch torque in response to increasing voltage indicated that opposing plantarflexor muscles were not being simultaneously stimulated [17].

The majority of torque produced by the dorsiflexors is from the tibialis anterior muscle (TA) [18], and we have shown previously that this model results in injury to this muscle [15, 16, 19, 20]. A custom program based on commercial software (Labview version 8.5, National Instruments, Austin, Tex) was used to synchronize contractile activation and the onset of ankle rotation. Injury resulted from 10 forced lengthening (plantarflexion) contractions through a 0°–70° arc of motion (with the foot orthogonal to the tibia, considered 0°). A sham procedure was performed whereby the identical protocol was performed (including stimulation), but without lengthening.

2.2. Outcome Measures

A maximal isometric contraction (200 ms duration) of the dorsiflexors was used to measure maximal torque before injury. For each animal, maximal isometric torque was also measured 5 minutes after injury (to measure force lost due to injury). All isometric contractions were performed with the ankle at 20° into plantarflexion, a position that results in optimal force production.

2.3. In Vivo Imaging

High resolution MRI (magnetic resonance imaging) was performed on a Bruker Biospin (Billerica, Mass) 7.0 Tesla MR system equipped with a 12 cm gradient insert (660 mT/m maximum gradient, 4570 T/m/s maximum slew rate) running Paravision 5.0 software to assess muscle damage on the day of injury. An MR-compatible small-animal monitoring and gating system (SA Instruments, Inc.) was used to monitor respiration rate and body temperature. Body temperature was maintained at 36-37°C using a warm water circulator. A custom-made holder was used to position the mouse in the supine position with both legs parallel to the bore of the magnet from knee to foot. Two hours after injury [21], animals were placed inside the scanner, with a four-channel receive-only surface coil placed anterior to the TA. Anesthesia was maintained using 2% isofluorane and was adjusted to maintain a safe respiration rate.

After imaging localizers, the following MR scans were performed: T1-weighted rapid acquisition with relaxation enhancement (RARE) with the following parameters: TE = 9.52 ms, TR = 1800, echo train length = 4, NA = 8, in-plane resolution m, and slice thickness = 750 μm; dual echo PD/T2 RARE: TE = 19.0/57.1 ms, TR = 5000 ms, echo train length = 4, NA = 1, in-plane resolution m, and slice thickness = 750 μm. Spin echo (SE) diffusion tensor image data was acquired using 12 noncolinear directions: -value = 350 s/m , TE = 26 ms, TR = 4500 ms, NA = 1, in-plane resolution m, and slice thickness = 750 μm. Multislice multiecho (MSME) T2 mapping image data was acquired in the same slice positions as the DTI data using 16 TEs = 11.4 ms to 182.5 ms with TE = 11.4 ms, TR = 10000 ms, NA = 1, in-plane resolution m, and slice thickness = 750 μm.

T2 mapping was performed using custom software written in MATLAB (The Mathworks; Natick, Mass) using nonlinear least squares to fit the measured data at each pixel to the canonical T2 signal equation. Diffusion tensor reconstruction and tractography was performed using TrackVis (Martinos Center for Biomedical Imaging; Massachusetts General Hospital; Boston, Mass) and MATLAB to calculate mean diffusivity (MD), fractional anisotropy (FA), radial diffusivity (RD), and axial diffusivity (AD) images. A region of interest (ROI) in the TA and tractography was used to guide bilateral segmentation of the TA. Diffusion tensor tractography was calculated using the fiber assignment by continuous tracking (FACT) method [22] with a termination criteria of an angle greater than 35°, followed by subsequent spline filtering. Tracts were restricted to those traveling through several transverse slices of manually traced regions of interest within each of the left and right TAs. This was used to create an image mask for each muscle, which was then divided into proximal, middle, and distal sections of approximately equal length. Finally, the masks were used to calculate average measurements of MD, FA, RD, AD, and T2 within each section of the injured and uninjured muscle. These measurements were compared between normal and dystrophic mice for both the uninjured and the injured side using a Wilcoxon Rank Sum test.

2.4. Histology

After MR and functional data were collected, anesthetized animals were fixed with 4% paraformaldehyde via perfusion through the left ventricle to preserve tissue morphology. TAs were harvested, weighed, snap frozen in liquid nitrogen, and then stored at −80°C. Animals not perfused fixed were euthanized by carbon dioxide inhalation followed by decapitation. To study the integrity of the muscle fiber membrane, additional mice ( each group) received an intraperitoneal injection of 1% Evans blue dye (EBD; Sigma, St. Louis, Mo) (wt/vol) in phosphate-buffered saline (PBS, pH 7.4) at a volume of 1% body mass (BM) (1 mg EBD/0.1 mL PBS/10 g BM). This solution was sterilized by passage through a Millex-GP 0.22 μm filter (Millipore, Bedford, Mass) and administered 24 h before death of the animal to assure a good signal [23]. EBD binds to albumin and is detected by fluorescence microscopy (at 568 nm) in the extracellular space. Presence of the protein-bound dye inside the muscle fiber indicates damage to the sarcolemma. Here, the sections were assessed blindly, and the myofibers were judged in a binary fashion, as positive or negative for intracellular EBD.

Transverse sections were cut on a cryostat (10 μm thickness) and collected onto glass slides (Superfrost Plus; VWR, West Chester, Pa). For EBD-injected tissue, myofibers with dye-labeled sarcoplasm were quantified in cross-sections under fluorescence optics [16, 24]. At least 400 fibers in 10 optical fields were assessed in 3 mice; results are expressed as percent labeled fibers. In sections from muscles without EBD, sections were stained with hematoxylin and eosin (H&E) for tissue evaluation. Sections were randomized and viewed at 100x magnification in a Zeiss Axioskop light microscope, and pictures were taken with a digital camera (AxioCam HR using AxioVision 3.0). Each optical field contained an average of fibers, and more than 45 fields were counted per muscle.

3. Results

To induce injury in the tibialis anterior muscle (TA), we superimposed a lengthening contraction onto a maximal isometric contraction (Figure 1(a)). To assess the amount of injury, we measured maximal torque before and after injury (Figure 1(b)). While the mechanism of injury was identical between animals, the injury was much more severe in dystrophic mice, with an average force loss of 85%, compared to 42% ( ) in normal mice.

Histological changes paralleled functional changes to some extent. For H&E staining, there was only minimal evidence of perivascular inflammation in the wild-type tissue after injury (Figure 2(b)) compared to undisturbed healthy muscle (Figure 2(a)). Even without injury, centrally nucleated fibers, fiber size variation, mild inflammation, and slight increases in endomysial connective tissue in the uninjured mdx mice are all consistent with ongoing muscle damage and regeneration (Figure 2(c)). In contrast, the myonecrosis, myophagocytosis, and foci of inflammation surrounding individual muscle fibers seen in the injured mdx muscle (Figure 2(d)) are all consistent with acute muscle injury. The amount of membrane damage, as evidenced by the number of fibers with intracellular EBD, also reflected the magnitude of functional injury (Figures 3(a)–3(e)). In wild-type tissue sections (Figure 3(a)), there were few fibers that could be found to have intracellular EBD (1%); however after the relatively mild injury protocol (Figure 3(b)), the number of fibers with membrane damage did not increase significantly ( %, ). The uninjured mdx tissue (Figure 3(c)) already had a small amount of EBD-positive fibers ( %, ), but there was almost a 5-fold increase in the number of fibers with membrane damage after injury ( %, , Figure 3(d)).

MRI revealed a clear increase on the T2 signal after injury (Figures 4(a) and 4(b)). Quantitative changes in T2 are apparent in Figures 4(c) and 4(d), where the expected T2 signal was fit across multiple echo times in T2-weighted images [25]. There was no significant difference in parameters between the injured and uninjured leg in the normal mice. Further, there were no significant differences in any of the measured parameters for the TA between normal and dystrophic mice in the uninjured side. However, there were marked differences in parameters in the injured side between normal and dystrophic mice and between the injured and uninjured leg in the dystrophic mice (Table 1). As in prior studies [14], the sham animals did not have an increase in T2 signal by the time imaging was performed.

Figure 5 shows a representative picture of 3D tractography. When comparing parameter of diffusion on the injured side, dystrophic mice showed significantly increased MD and AD and decreased FA ( ) in the proximal and middle TA compared to wild-type uninjured mice. These differences were much more marked than the more limited differences we measured with T2. While not significant, RD trended toward increased values. These data are plotted against force loss in Figure 6. DTI tractography (see Figure 5, e.g.) was used to visualize the TAs.

4. Discussion

Since X-rays reveal little regarding skeletal muscle anatomy or function, muscle injury is still typically diagnosed by taking a history and performing a physical examination. However, modalities such as MRI are becoming more affordable and more commonplace, and there is a surge of interest in using them to learn more about both muscle injury and muscle disease. Acute muscle strain injury is a stochastic event, and that would not be ethical or desirable to induce in human subjects; therefore most human studies involve imaging after “eccentric exercise” or “exertional muscle injury” (i.e., overuse type injuries). By using an established animal model to induce an acute contraction-induced injury, we were able to compare changes in muscle function, in vivo imaging, and histology shortly after injury. Our injury method clearly resulted in muscle damage, as indicated by the measurable loss of force, changes in H&E staining, and an uptake of EBD by injured muscle fibers, as well as marked changes in MRI parameters, especially DTI.

There is a growing number of animal studies that use MRI to assess muscle injury; however, these include such methods as overuse (running downhill), myotoxins, denervation, hindlimb suspension reloading, and ischemia-reperfusion [8, 12, 2631], rather than by acute muscle strain injury. This can make reproducibility and even relevance to functional outcome difficult to determine. Our in vivo animal model of injury provides several advantages, such as the ability to study exclusively one type of contraction (eccentric), to measure force directly from an individual muscle, and to characterize the biochemical and morphological state of the muscle at defined times following injury and during recovery. Because the animal is anesthetized, we can reliably induce damage to a single muscle group with a known velocity, arc of motion, and contractile level. Effort, such as in a maximal voluntary contraction (MVC), and pain are not confounding factors when assessing changes in function. Downhill treadmill running is an exercise that is sometimes used as an injury model, but there are problems with compliance [32]. The in vivo injury model provides a known dose of injury to a specific muscle group in all animals.

Regardless of the method used to induce damage, the T2 signal in MRI is typically increased significantly after muscle injury, but after some protocols there is a long-term persistence of this signal [33, 34], with a time scale that is not in line with the temporal changes of inflammation and edema [28]. Although some studies have shown that T2 values gradually peak hours or days after injury [26, 30], these employed a myotoxin to injure the muscle. Myotoxins provide a model to study necrosis, inflammation, and massive degeneration and regeneration but do not provide a physiological model of injury. In a recent paper, Mathur et al. [32] used downhill running to study the changes in T2 after eccentric exercise of mdx mice. They used older animals (5–15 months) but still found that mdx mice were more susceptible to damage, with T2 returning to preexercise levels 10 days after exertion. Interestingly, they noted a strong relationship between the percent of EBD-positive area and the percentage of pixels with elevated T2 ( ). EBD binds to albumin, a blood plasma protein, and is detected by fluorescence microscopy in the extracellular space. Presence of the protein-bound dye inside a muscle fiber indicates damage to the sarcolemma. Since albumin is a relatively large protein (~67 kDa), their findings suggest that downhill running results in relatively large sarcolemmal tears of damaged myofibers.

McIntosh et al. [35] provided one of the first studies to use MRI to assess skeletal muscle in mdx mice. These authors noted heterogeneous signal intensity on T2-weighted images in mdx hindlimb muscles. These foci of high intensity were subsequently confirmed by Walter et al. [36], who used more advanced techniques to obtain higher resolution images as well to quantify the changes in the T2 signal that occurred in the mdx (and sg ) mice compared to controls. Walter et al. also showed that delivery of sg by AAV (adenoassociated virus) into hindlimb muscles was not only able to rescue to phenotype (i.e., return of normal histology) but also restored the T2 properties of skeletal muscle.

DTI parameters have been compared to longitudinal sections of the rat TA, the results of which show that DTI directions actually represent the local orientation of fibers in the rat TA muscle [37]. Damon et al. have used DTI fiber tracking successfully to measure pennation angles of myofibers in human skeletal muscle [10, 11, 38, 39], and subsequent studies have also used DTI in mice, both for three-dimensional reconstruction of fiber tracks [40] and as a measure of damage induced by ischemia reperfusion [8]. Heemskerk et al. provide an excellent description of measured changes in DTI parameters after ischemia-reperfusion type of injury in mouse skeletal muscle [12]. They found that DTI indices changed in response to ischemia-reperfusion and that the indices correlated with histopathology, although they caution against overinterpretation since their study was not designed to achieve exact correlation. Such correlation requires optimal MRI protocols, an exact spatial correspondence between histology and MRI slices, and specific immunohistochemical labeling to identify specific cell types.

One can use DTI to compare muscle architecture of dystrophic and healthy hindlimb muscles of mice in vivo. This could be extremely valuable to muscle physiologists who need to obtain the cross-sectional area (CSA) of a muscle to calculate the “specific force” (force normalized to CSA) without removing the muscles. In addition to the computer-generated models of fiber tracks, we were able to collect quantifiable data. Using our injury model, the diffusion parameters (decrease in FA and increases in MD and AD) and T2 findings are consistent with increased edema. However, the lack of significant changes in RD may suggest increased diffusion along the sarcolemma as a result of cell swelling. Interestingly, significant changes in DTI parameters were evident in the middle and proximal sections of the TA, and significant T2 changes were evident in the proximal TA. Such findings support earlier work suggesting that damage to the myofibers is not limited to the muscle-tendon junction, where soreness often occurs [14].

There are certain limitations to any study of muscle injury and disease. In addition to technical issues involved with imaging [41], position of the limb [42], and the lack of any single biological marker to explain the loss of force after injury [43, 44], one must still use caution when interpreting T2 intensity or DTI findings. The exact meaning of T2 changes are still to be elucidated [28, 32, 45]. While tractography corresponds to normal anatomy, there is a tendency after injury for the more superficial fascicles to appear abnormal (Figure 5, arrows). This three-dimensional reconstruction can appear to have fibers that “punch through” the strong connective tissue (epimysium) that surrounds that muscle and enter into the space occupied by another muscle, when indeed this is not the case. This is likely to be even more of a problem when using DTI on patients with DMD, as they have significant progressive damage over time and a gradual replacement of muscle tissue by fat, something that does not occur in the mdx mouse model of DMD. Such issues are still being worked on and can hopefully be resolved. However, this change in normal direction of myofibers is clearly indicative of damage, and tractography in healthy muscle is still a useful tool to obtain volume in order to calculate “specific force” (force per unit cross-sectional area).

In summary, we have used MRI to study injured muscle by observing edema and muscle fiber diffusion changes with T2 and DTI sequences, respectively. DTI exhibited significant changes in both the middle and proximal sections of the TA, whereas T2 changes were only significant in the proximal section. This is, to the best of our knowledge, the first animal study to use DTI in a reliable and valid animal model of eccentric injury. The results suggest that DTI is a reliable marker in the assessment of acute muscle injury, even at early time points where the MR signal changes are dominated by local edema.

Acknowledgment

This work was supported by Grants to RML from the National Institutes of Health (K01AR053235 and 1R01AR059179).