585357

research-article2015

NNRXXX10.1177/1545968315585357Neurorehabilitation and Neural RepairNeurorehabilitation and Neural RepairCheung et al

Original Article

Lesion Characteristics of Individuals With Upper Limb Spasticity After Stroke

Neurorehabilitation and Neural Repair 1­–8 © The Author(s) 2015 Reprints and permissions: sagepub.com/journalsPermissions.nav DOI: 10.1177/1545968315585357 nnr.sagepub.com

Daniel K. Cheung1,2,3, Seth A. Climans4, Sandra E. Black, MD, FRCP1,2,3,5, Fuqiang Gao, MD2, Gregory M. Szilagyi2,3, and George Mochizuki, PhD1,2,3

Abstract This study explores the relationship between lesion location and volume and upper limb spasticity after stroke. Ninetyseven stroke patients (51 with spasticity) were included in the analysis (age = 67.5 ± 13.3 years, 57 males). Lesions were traced from computed tomography and magnetic resonance images and coregistered to a symmetrical brain template. Lesion overlays from the nonspastic group were subtracted from the spastic group to determine the regions of the brain more commonly lesioned in spastic patients. Similar analysis was performed across groups of participants whose upper limb (elbow or wrist) Modified Ashworth Scale (MAS) score ranged from 1 (mild) to 4 (severe). Following subtraction analysis and Fisher’s exact test, the putamen was identified as the area most frequently lesioned in individuals with spasticity. More severe spasticity was associated with a higher lesion volume. This study establishes the neuroanatomical correlates of poststroke spasticity and describes the relationship between lesion characteristics and the severity of spasticity using mixed brain imaging modalities, including computed tomography imaging, which is more readily available to clinicians. Understanding the association between lesion location and volume with the development and severity of spasticity is an important first step toward predicting the development of spasticity after stroke. Such information could inform the implementation of intervention strategies during the recovery process to minimize the extent of impairment. Keywords stroke, spasticity, lesion–symptom mapping

Introduction One common consequence of stroke is spasticity—a velocity-dependent increase in tonic stretch reflexes.1 Upper limb spasticity is associated with poorer function,2 in part due to restrictions in range of motion and speed of limb movements.3 As a positive symptom of the upper motor neuron syndrome, poststroke spasticity occurs due to changes in the net balance of excitatory and inhibitory inputs on spinallevel circuitry secondary to injury in higher order centers or in their descending pathways (ie, pyramidal or parapyramidal fibers).4,5 While the physiological mechanisms that lead to spasticity are generally understood, the link between the neuroanatomical correlates and clinical observations are less established. Little is known about the characteristics of stroke lesions (ie, location and volume) that modulate muscle tone and stretch sensitivity; although, recent reports associate severe spasticity (Modified Ashworth Scale [MAS] ≥ 2) with injury to subcortical structures.6 In spite of this, there is variability in brain regions identified as being associated with the observed presence of spasticity. For example, Ryu and colleagues reported that 87/104 (84%) of

the individuals with spasticity presented with ischemic or hemorrhagic strokes involving the anterior or middle cerebral arteries.7 Other authors have identified the basal ganglia, thalamus, and capsule as specific neuroanatomical substrates of spasticity based on lesion overlays in 4 subjects with spasticity.8 This apparent variability in the anatomic correlates of spasticity is a potential contributor to the variability in clinical presentation, assessment, and efficacy of spasticity management protocols. Brain lesion profiles, which combine both lesion location and lesion volume, are important predictors of the severity of motor impairments. These profiles reveal stronger 1

Heart and Stroke Foundation Canadian Partnership for Stroke Recovery, Toronto, Ontario, Canada 2 Sunnybrook Research Institute, Toronto, Ontario, Canada 3 University of Toronto, Toronto, Ontario, Canada 4 Western University, London, Ontario, Canada 5 Sunnybrook Health Sciences Centre, Toronto, Ontario, Canada Corresponding Author: George Mochizuki, Sunnybrook Research Institute, 2075 Bayview Avenue, Rm M6-178, Toronto, Ontario, Canada M4N3M5. Email: [email protected]

Downloaded from nnr.sagepub.com at WESTERN OREGON UNIVERSITY on May 25, 2015

2

Neurorehabilitation and Neural Repair 

correlations to functional outcomes than lesion size only.9 Lesion location (independent of volume) is also associated with the extent of motor recovery. This is apparent from reports linking recovery with lesions sparing the motor cortex10 or those reporting greater motor impairment with lesions involving the corticospinal tracts.11 While lesions have functional and location-dependent consequences in relation to paresis, a similar relationship with spasticity has not been clearly identified. Observing differences in lesion overlaps between stroke survivors with and without spasticity is one approach to probing this relationship. Analysis using lesion–symptom mapping can reveal relationships between sites of brain injury and behavioral outcomes and can be used to associate stroke infarct characteristics with specific behavioral phenotypes.12 Using computed tomography (CT) or magnetic resonance imaging (MRI), lesion–symptom mapping has revealed relationships between brain lesion characteristics and functional outcomes, motor recovery,9 and gait,12 though by the correlational nature of such studies, causality cannot be claimed. Nevertheless, evaluating these relationships can expand our understanding of what characteristics of a stroke lesion may correlate with spasticity and could help in clinical management. In an effort to systematically examine the relationship between the presence and severity of spasticity with stroke lesion size and location, the current study aimed to investigate the neuroanatomical correlates of upper limb spasticity after stroke through the use of clinical scans and manual lesion tracings. A better understanding of the association between brain lesion profile and severity of spasticity could complement predictive models for spasticity,13 which, in turn, could affect the implementation of spasticity-management protocols in the early stages after stroke when the potential for recovery may be the highest. Given the pathways implicated in the pathophysiology of spasticity14 and neuroanatomical correlates of severe spasticity,6 we hypothesized that subcortical lesions would better correlate with spasticity after stroke than cortical lesions.

Methods Subjects Study data were retrospectively drawn from the Heart and Stroke Foundation Canadian Partnership for Stroke Recovery Longitudinal Database. This database is used to prospectively track stroke recovery longitudinally from 5 partnering rehabilitation hospitals using a battery of standardized clinical measures. Participants in the database provided written informed consent to participate in the protocol, which was approved by the research ethics board at each institution. Data from 297 participants were available at the time of analysis. The database was queried for participants

between 3 months and 2 years poststroke at the time of assessment with both (a) a documented MAS score ranging from 0 to 4 of the upper limb (wrist or elbow) as measured by a physiotherapist and (b) an available axial, non-contrast CT or MRI scan. Though mixing the CT and MRI scans may contribute to variability in lesion volumes and locations depending on the imaging modality used, methods for normalizing brain images to templates based on mixed modalities (CT and MRI) have been established.15 We believed it was important to include CT scans as they are readily available through routine clinical care. If an individual had differing spasticity scores at the elbow and wrist, the higher of the 2 scores was used. Subjects with a known stroke prior to the baseline stroke or transient ischemic attack were excluded. Subjects taking antispastic medication were also excluded so as to avoid biasing the MAS scores toward the lower end of severity. Available clinical history was assessed and subjects with confounding nonstroke pathologies or unusable scans were eliminated from the subject pool (see Results).

Scan Selection Scan selection was based on scan timing, lesion visibility, scan sensitivity, and scan modality. Preference was given to scans taken closest to the assessment date, and MRI was preferred over CT if available. Within MR imaging, T1-weighted imaging was preferred to Fluid Attenuated Inversion Recovery (FLAIR) or T2-weighted imaging, which in turn was preferred over diffusion-weighted imaging.

Lesion Tracing Lesions were manually traced, slice-by-slice, using an inhouse modified version of ITK-SNAP22 called ITKSNAP-SB. Lesion tracers had access to radiology reports, but were blind to clinical status and spasticity scores. Lesion core, perinecrotic tissue, lacunes, and contralateral lesions were included in the tracings, consistent with procedures previously used by our group.12 Lacunar infarction was distinguished from enlarged Virchow–Robin spaces using criteria set out by Bokura and colleagues.16 Tracings were reviewed and confirmed by an experienced research radiologist (FG).

Template Alignment Skulls were digitally stripped away from each scan to normalize brains onto a common template. For CT scans, skull-stripping was accomplished with slice-by-slice semi-automated seeding and manual limit-setting within Analyze10 (Biomedical Imaging Resource, Mayo Clinic, Rochester, MN). For high-resolution T1-weighted scans, the automated algorithm Brain Extraction Tool17 was used to remove the

Downloaded from nnr.sagepub.com at WESTERN OREGON UNIVERSITY on May 25, 2015

3

Cheung et al skull from the images. Manual editing of the resultant brain mask was then performed in ITK-SNAP-SB. For FLAIR and T2-weighted scans, skull-stripping was accomplished sliceby-slice using a manual outline tracing within Analyze10. All brains were aligned using an affine transformation to the symmetric T1-weighted template brain available through the International Consortium for Brain Mapping in 2009.18 This was accomplished by using Automated Image Registration19 on the skull-stripped brains. The traced lesions were then transformed into template space. The resultant transformations were reviewed by the research radiologist (FG). Gaussian smoothing (1 mm radius) was performed on the lesions. Using MRIcro software (Rorden, Columbia, SC), region of interest (ROI) images were created for each lesion and flipped onto the right hemisphere. Images with similarly sized lesions in both hemispheres remained unflipped. Primary hemisphere (if any) of stroke was determined from radiologic reports.

Lesion Analysis A linear sum of the binary tracings was calculated based on the lesion overlap at each voxel. This resulted in a lesion density map. The Talaraich coordinates acquired within MRIcro were exported to determine the specific brain structures associated with the lesion overlays, specifically the voxel with the most overlaps. A variant of voxel-based Lesion-Symptom Mapping20 adapted for categorical lesion status and spasticity data using a combination of in-house lesion-mapping software and statistical analyses performed in version 3.1.0 of the R statistical computing environment (R Core Team, 2014) was used. Compiled across the template-space lesion masks of all subjects, the list of binary class labels (0 = intact, 1 = lesioned) at a given voxel coordinate formed a “patient– lesion pattern,”21 dividing subjects into 2 groups with respect to the presence or absence of a lesion in the same approximate anatomical location. Because discretely labeled lesions (“1s”) form contiguous regions within each patient, the spatial nonindependence of voxel-wise lesion status is also reflected in the tendency of neighboring voxels to share patient–lesion patterns in the aggregate data. After identifying all unique patient–lesion patterns in template space, R was used to perform Fisher’s exact test on contingency tables combining each of the corresponding distinct groupings of observations with binarized spasticity scores (0 = MAS = 0, 1 = MAS ≥ 1). All voxels with a lesion in at least 20% of all subjects, spastic or nonspastic, were included in the analysis. This threshold was chosen to achieve a minimum detectable effect size in accordance with the prevalence of spasticity at varying time points of recovery (1 day to 6 months) reported to range between 19% and 27%.22,23 A statistical map of P values was generated by assigning the result from each patient–lesion partitioning of the data to all corresponding voxels. The final maps were filtered using a P value

threshold, Bonferroni-corrected with respect to the total number of unique patient–lesion patterns in the data set.21 Using SPSS (IBM SPSS Statistics 21), Pearson’s χ2 was used to determine whether the spastic and nonspastic groups differed in the type of stroke causing the lesion (ischemic or hemorrhagic) or in the imaging modality used to characterize the lesion (CT or MR).

Results Patient Demographic and Clinical Characteristics In total, 102 subjects met the initial inclusion criteria. Following initial examination of available brain scans, 5 subjects were excluded—3 because the available brain imaging showed insufficient evidence of stroke and 1 due to severe herpesviral encephalitis. The other subject was excluded because the scans had large portions of skull missing due to a decompressive craniectomy. Ninety-seven subjects (57 male; 46 with MAS of 0, 51 with MAS of 1, 1+, 2, 3, or 4) were included in the analysis (Table 1). Overall, the average age of participants included in the analysis was 67.5 ± 13.3 years. The median time between date of stroke and date of scan (64 CT, 5 T1, 10 T2, 18 FLAIR) for all participants was 4 days (range = 0-608), and the median time between date of stroke and assessment date was 187 days (range = 122-627). The median National Institutes of Health Stroke Scale total score for all participants was 3 (range = 0-12), and the median Chedoke-McMaster Stroke Assessment arm and hand score was 4 (range = 1-7) and 5 (range = 1-7), respectively. The majority of participants experienced ischemic stroke (11 cases of lacunar stroke) and 19 had primary hemorrhagic stroke. Neither the type of lesion (ischemic or hemorrhagic) nor the type of imaging (CT or MR) was significantly associated with the presence of spasticity, χ2(1) = 1.83, P = .176, and χ2(1) = 0.37, P = .544, respectively. Demographic data were comparable between the spastic and nonspastic groups. Overall, the spastic and nonspastic groups differed in clinical presentation. Individuals with spasticity were characterized as having more severe strokes and greater impairment in the hand and arm. Clinical scores are presented in Table 1.

Lesion Analysis: Differences Between the Spasticity and No Spasticity Groups When the brain lesion overlays of nonspastic group were subtracted from the spastic participant group, the resulting lesion map revealed the basal ganglia to be the brain region with the most overlaps. Following Fisher’s exact test of each individual subdivision, the resulting P values were plotted using MRIcro, and the basal ganglia—more specifically, the putamen—was confirmed as the area of peak

Downloaded from nnr.sagepub.com at WESTERN OREGON UNIVERSITY on May 25, 2015

4

Neurorehabilitation and Neural Repair 

Table 1.  Demographic and Clinical Data of Stroke Patients Based on Modified Ashworth Score of the Upper Limba.

N Sex, male/female Age, years Lesion side, right/left/both Scan type, MRI/CT Scan date, days poststroke*,e Assessment date, days poststroke*,e MAS assessment site, wrist/elbow/same NIHSS total* CMSA arm* CMSA hand*

Nonspastic

Spastic

46 27/19 70.0 (13.1) 18/26/2 16/30 3 (0-448) 190 (126-372) — 1 (0-12) 6 (2-7) 6 (2-7)

51 30/21 65.3 (13.2) 33/17/1 17/34 4 (0-608) 184.5 (122-627) 20/18/13 4 (0-11) 3 (1-7) 3 (1-6)

Pb .990c .086d

0.273f 0.312f

Lesion Characteristics of Individuals With Upper Limb Spasticity After Stroke.

This study explores the relationship between lesion location and volume and upper limb spasticity after stroke. Ninety-seven stroke patients (51 with ...
410KB Sizes 2 Downloads 7 Views