Abstract
Introduction
Hemiparesis after stroke is particularly common, with up to 70% of individuals experiencing mild to severe upper extremity (UE) dysfunction (Nakayama, Jorgensen, Raaschou, & Olsen, 1994). In individuals with stroke, it has been suggested that the extent of structural damage to cortical areas likely underlies the potential for recovery (Schaechter, Moore, Connell, Rosen, & Dijkhuizen, 2006). Therefore, assessing morphological differences within the motor cortex may provide valuable information of the neural events that underlie UE hemiparesis.
After stroke, morphological and structural changes are triggered spontaneously at the onset of infarction and may persist into the chronic phase post stroke (Dimyan & Cohen, 2011). Spontaneous recovery of motor function appears to be related to lesion location (Shelton & Reding, 2001), however there is a weak relationship between lesion location and UE functional recovery after stroke (Gauthier et al., 2014; Gauthier, Taub, Mark, Perkins, & Uswatte, 2009). Though the survival of peri-infarct tissue is an important component of recovery (Cramer, Shah, Juranek, Crafton, & Le, 2006), it is becoming clear that regions distant from the lesion are also important for recovery of motor function. For instance, after subcortical stroke, recovery of arm function has been linked with greater functional activity within cortical areas remote from the lesion (Cramer, 2008; Cramer et al., 1997; Weiller, Ramsay, Wise, Friston, & Frackowiak, 1993). Further, functional imaging of individuals with chronic stroke shows that activation of the ipsilesional hemisphere is associated with improved motor outcome (Calautti et al., 2010; Ward, Brown, Thompson, & Frackowiak, 2003) and structural integrity of cortical areas (Schaechteret al., 2006). It appears that morphological changes and the structural integrity of surviving cortical areas contribute to the overall magnitude and pattern of recovery of UE function in individuals with chronic stroke.
Metabolic changes in the cortex after stroke have previously been investigated using biological markers of membrane integrity and neuronal damage, and have been described in individuals with both acute and chronic stroke (Craciunas et al., 2013; Federico et al., 1998). One metabolite NAA, is considered to reflect neuronal structural and functional integrity (Demougeot, Marie, Giroud, & Beley, 2004). Magnetic resonance spectroscopic (MRS) studies of individuals in the acute stage post stroke have linked decreases in NAA concentration to poor functional recovery of the hemiparetic UE(Federico et al., 1998). NAA changes have also been observed in individuals with chronic stroke within ipsilesional normal appearing grey matter (Munoz Maniega et al., 2008) and have been associated with morphological changes, including neuronal death in stroke (Demougeot et al., 2004). In individuals with chronic stroke, lower NAA concentrations have been reported within both ipsi- and contra-lesional primary motor cortices (Cirstea et al., 2011; Craciunas et al.,2013).
Another metabolite that has been investigated in individuals with chronic stroke is glutamate (and its readily available precursor glutamine, collectively referred to as Glx). In a study by Cirstea et al. (2011), ipsilesional concentrations of NAA and Glx in normal appearing grey matter were both positively correlated with arm motor impairment in individuals with chronic stroke.
While MRS has provided information about neuronal integrity and metabolism after stroke, it remains unclear how regional structural differences, including cortical thickness, may be associated with metabolic function after stroke. Increases in cortical thickness in the sensorimotor cortical areas have been reported in the contralesional hemisphere of individuals with sub-acute (Brodtmann et al., 2012) and chronic stroke (Schaechter et al., 2006). Voxel-based morphometric investigation in individuals with chronic subcortical stroke revealed grey matter atrophy in motor areas, which may be associated with hemiparetic arm impairment (Gauthier, Taub, Mark, Barghi, & Uswatte, 2012). In the current study, we performed metabolic and structural assessment of the primary motor cortex using H1 magnetic resonance spectroscopy and cortical thickness measurement to investigate whether morphological differences are associated with UE motor performance.
To our knowledge, no study has investigated how differences in cortical thickness relate to metabolic changes in individuals with chronic subcortical stroke. Total NAA (NAA and N-Acetylaspartylglutamate (NAAG)) and Glx were selected as metabolites of interest a priori because: 1) both are altered in individuals with chronic stroke, 2) the magnitude of change in these has been associated with the level of UE impairment (Cirstea et al., 2011, 2012), and 3) an observed decrease in cortical thickness may be explained by some form of neuronal cell degradation, which could be indexed using tNAA (Munoz Maniega et al.,2008).
The primary aim of this study was to characterize metabolic and morphometric properties of the precentral gyrus. In addition we considered how these properties related to UE motor impairment in individuals with chronic subcortical ischemic stroke. Our hypothesis was that individuals with stroke would exhibit reduced precentral gyrus cortical thickness and lower tNAA and Glx concentrations in the ipsilesional hemisphere. Further, we hypothesized that measures of ipsilesional cortical thickness, and tNAA and Glx concentrations, would be associated with UE motor function after stroke.
Materials and methods
Participants
Seventeen individuals (12 male, 5 female; age 68 +/–9.6 years) with subcortical ischemic stroke in the chronic phase (>6 months) of recovery, and 11 neurologically healthy controls (4 male, 7 female; age 60.4 +/–6.0 years) were recruited for this study. The research ethics board at the University of British Columbia approved all aspects of this work and consent was obtained according to the Declaration of Helsinki. Individuals were excluded if they: 1) were in the acute (0–3 months) or sub-acute (3–6 months) phase of recovery, 2) had contraindications to MRI, 3) had previous history of stroke, epilepsy, neurodegenerative disorder or head trauma, or 4) showed a visible lesion extending in to the precentral gyrus grey matter. Upon enrolment, each participant underwent neuroimaging assessment on a Philips Achieva 3.0T MRI scanner at the University of British Columbia, which was followed by functional assessment in the Brain Behaviour Lab.
Magnetic resonance imaging acquisition
A high-resolution T1-weighted scan (TR/TE: 7425/3.64 ms, FOV = 256 × 256 mm, 160 sagittal slices, flip angle 6°, voxel size = 1 mm3) was performed on each enrolled participant. In addition, single voxel 1H-MRS (PRESS, TR/TE = 2000/35 ms, sampling frequency: 2000 Hz, data points = 1024, signal averages = 128) was obtained to measure absolute metabolite concentrations. A single 30 mm × 22 mm × 15 mm voxel was centered over the hand area of the primary motor cortex in both ipsilesional and contralesional hemispheres. Projection based autoshimming was carried out. Unsuppressed water spectra were acquired to enable metabolite concentration estimations. Localization of hand representation for each individual and corresponding voxel placement was determined using known anatomical coordinates (Yousry et al., 1997) (Fig. 1).
Cortical morphometry analysis
Cortical reconstruction and volumetric segmentation was performed with the Freesurfer image analysis suite, which is documented and freely available for download online (http://surfer.nmr.mgh.havard.edu). The technical details of these procedures are described in prior publications (A. M. Dale, Fischl, &Sereno, 1999; Fischl & Dale, 2000). A number of steps for processing are required, and include: 1) motion correction and averaging (Reuter, Rosas, & Fischl, 2010) of multiple volumetric T1 weighted images (when more than one is available), 2) removal of non-brain tissue using a hybrid watershed/surface deformation procedure (Segonne et al., 2004), 3) segmentation of the subcortical white matter and deep grey matter volumetric structures (Fischl et al., 2002; Fischl et al., 2004), 4) intensity normalization (Sled, Zijdenbos, & Evans, 1998), 5) tessellation of the grey matter white matter boundary, and 6) automated topology correction (Fischl, Liu, & Dale, 2001; Segonne, Pacheco, & Fischl, 2007) and surface deformation following intensity gradients. This final step is done to optimally place the grey/white and grey/cerebrospinal fluid borders at the location where the greatest shift in intensity defines the transition to the other tissue class (A. M. Dale et al., 1999; Anders M. Dale & Sereno, 1993; Fischl & Dale, 2000). MRI-derived measures of human cerebral cortical thickness have been shown to be highly reliable (Han et al., 2006), and have been described in a number of pathological conditions such as Huntington’s disease (Rosas et al., 2002), schizophrenia (Kuperberg et al.,2003) and multiple sclerosis (Sailer et al., 2003). Recently the technique has been used to investigate regional differences in cortical thickness in sub-acute stroke (Brodtmann et al., 2012) and chronic stroke populations (Borich, Neva, & Boyd, 2015). We compared average cortical thickness in the precentral gyrus in both Stroke (ipsilesional/contralesional) and Control (non-dominant/dominant) groups.
Magnetic resonance spectroscopy analysis
Metabolic measures were quantified using LCModel version 6.2-2B which provides automatic quantification of in vivo proton magnetic resonance spectra (Provencher, 1993). Water scaling was used to estimate absolute metabolite concentrations using the reference area of the unsuppressed water signal along with the assumed water concentration of 35880 mM within the largely grey matter voxel. Concentrations of tNAA and Glx were extracted as the primary metabolites of interest. Technical details can be found in Supplementary data.
Upper-extremity motor assessment
UE motor performance was assessed using the WMFT, which measures how rapidly individuals are able to perform a series of 15 functional tasks in the laboratory. A trained physiotherapist carried out all components of the assessment. Recently, Hodics et al. (2012) validated an approach that characterizes the WMFT by converting the time to complete each task into task rate data using the following equation:
This method of calculation has shown to be more sensitive in assessing the functional abilities of more severely affected individuals with stroke. Thus, participants in the current study were assessed using the standard 15-task WMFT; subsequently we calculated task-rate data (Hodics et al. (2012).
Statistical analysis
To examine differences in tNAA and Glx concentration and precentral gyrus thickness a between group (Stroke, Control) multivariate analysis of variance (MANOVA) was conducted. Follow-up univariate analyses of variance (ANOVAs) for each dependent variable were conducted. Parametric correlation analysis using Pearson’s r was conducted to examine the relationship between tNAA and Glx and precentral gyrus thickness in both hemispheres. Multiple linear regression analyses were performed to assess the amount of variance in UE motor performance (WMFT rate) explained by demographic information (age and post-stroke duration), metabolite concentration (tNAA and Glx), and cortical thickness. In the Stroke group, age and post stroke duration were entered first to account for post-stroke (Graham et al., 1993) and age related changes in metabolites (Angelie et al., 2001; Kaiser, Schuff, Cashdollar, & Weiner, 2005) and cerebral tissue (Resnick, Pham, Kraut, Zonderman, & Davatzikos, 2003). Concentration of ipsilesional cerebral metabolites tNAA and Glx were the next predictors entered, based on previous reports suggesting their association with UE motor impairment (Cirstea et al., 2011, 2012). Ipsilesional cortical thickness was the final predictor entered into the model. When comparing between the Stroke and Control groups, data from the ipsilesional cortex in individuals with stroke were compared to data extracted from the non-dominant cortex in our Control group.
Results
Participants
Participant demographics from the Stroke group (Table 1) and the Control group (Table 2) are listed below. The location of each lesion was determined from T1-weighted MRI.
Anatomic and metabolic differences in chronic stroke
There was an overall main effect of Group across the measures (F = 4.972; p = 0.004). Concentrations of tNAA (F = 6.47; p = 0.017), Glx (F = 10.37; p = 0.001), and cortical thickness (F = 6.12; p = 0.021) were all significantly lower bilaterally in the Stroke group compared to the Control group. Within the Stroke group, significantly less ipsilesional tNAA concentration (F = 4.14; p = 0.05), Glx (F = 11.72; p = 0.002) and cortical thickness (F = 6.37; p = 0.017) was also observed compared to the contralesional hemisphere (Fig. 2). No significant difference between hemispheres was found in the Control group for any of the variables (Fig. 3).
Relationship between cortical thickness and metabolite measures
A significant positive correlation was found between tNAA and precentral gyrus thickness in both ipsilesional (r = 0.782; p < 0.001) and contralesional (r = 0.517; p = 0.033) hemispheres of the Stroke group (Table 3, Fig. 4). No correlations between Glx concentration and precentral gyrus thickness was observed. No correlations were observed in theControl group.
Relationship between metabolite concentration, cortical thickness and upper extremity motor performance
Regression analysis results are summarized in Table 4 and Fig. 5. In the Stroke group, after adjusting for age and post-stroke duration, ipsilesional metabolite concentrations (tNAA and Glx) accounted for 53% of the total variance of hemiparetic UE WMFT score (R2 = 0.535, p = 0.043). tNAA concentration explained a larger amount of variance in WMFT motor performance than Glx. Adding ipsilesional precentral gyrus cortical thickness to the model did not explain an additional amount of unique variance in WMFT score, and the overall model was not significant (R2 change = 0.005, p = 0.088).
Discussion
We observed thinner cortex in the hand area of the precentral gyrus, as well as lower tNAA and Glx concentration in individuals with stroke, particularly in the ipsilesional hemisphere. Lower bilateral tNAA concentration was associated with a less cortical thickness bilaterally in the Stroke group. Our results suggest that ipsilesional precentral gyrus thickness and tNAA concentration were associated with UE motor performance.
Metabolite concentrations
To our knowledge, this is the first study to report significantly lower Glx concentration in the precentral gyrus of individuals in the chronic stage post-stroke. We also observed a lower concentration of the cerebral metabolite tNAA, which is consistent with previous findings in individuals with stroke (Cirstea et al., 2011; Craciunas et al., 2013). Because NAA and NAAG are believed to exist almost exclusively in intact neurons, a possible explanation for our finding of abnormal Glx and tNAA concentrations includes neuronal death (Pereira et al., 1999).Neuronal death following ischemic stroke is often attributed to improper homeostatic balance of glutamate. The release of glutamate following ischemia results in a cascade of cellular events such as intracellular calcium increase and mitochondrial dysfunction that ultimately leads to cellular death (Hazell, 2007). Acutely after stroke, glutamate concentration increases for approximately 6 hours after the onset of infarction (Davalos, Castillo, Serena, &Noya, 1997); this may cause permanent damage. Because lower ipsilesional tNAA and Glx was accompanied by a thinner precentral gyrus, neuronal death may, in part, explain our results. Supporting this possibility, age-related brain shrinkage has been associated with decreased NAA and Glx concentration in the motor cortex (Kaiser et al., 2005). Similarly, in the current study we observed a positive correlation between tNAA concentration and precentral gyrus thickness.
Cortical thickness
Precentral gyrus thickness was also significantly lower in the Stroke group compared to the Control group, and ipsilesional thickness in the Stroke group was significantly associated with UE motor performance. These results are consistent with previous reports suggesting that bilateral diffuse tissue loss may occur in individuals with chronic stroke (Gauthier et al., 2012; Kraemer et al., 2004).Further, decreases in grey matter density in non-infarcted motor regions in individuals with chronic stroke correlates with WMFT of UE motor function (Gauthier et al., 2012). Interestingly, these results differ from previous reports of increased contralesional sensorimotor cortical thickness in individuals with chronic stroke (Schaechter et al., 2006), and greater contralesional paracentral and superior frontal cortical thickness in individuals with sub-acute stroke(Brodtmann et al., 2012). It is possible that these differences could be explained by compensatory mechanisms within the contralesional hemisphere, such as structural plasticity within cortical regions in response to brain injury. Functional reorganization of contralesional motor cortical areas has been described previously in animal models of stroke undergoing rehabilitation training (Nudo, Wise, SiFuentes, & Milliken, 1996), and has been linked with improved hand function. There is evidence to suggest that contralesional compensatory changes are dependent upon increased use of the less-impaired limb (Nudo, Plautz, & Frost, 2001). This idea is supported by previous work that demonstrated a relationship between grey matter changes and arm motor impairment in chronic stroke patients undergoing specific rehabilitation programs (Gauthier et al., 2008). Indeed, individuals engaged in Constraint-Induced Movement Therapy (CIMT), a rehabilitation therapy that combines restraint of the non-hemiparetic limb with intensive use of the hemiparetic limb, show significant increases in grey matter in sensory and motor areas both contralateral and ipsilateral to the hemiparetic arm (Gauthier et al., 2008). Unfortunately, it is common in stroke survivors for daily use of the hemiparetic arm to be severely limited (Rand & Eng, 2015), therefore, it is possible that non-use of the hemiparetic arm might explain the observed differences in ipsilesional precentral gyrusthickness.
Upper extremity motor assessment
In our regression analyses, age and time post stroke alone were not significantly associated with UE motor performance. When these variables were combined with tNAA and Glx concentration the model was able to explain a significant amount of variance in hemiparetic UE motor function. tNAA explained the greatest amount of variance in motor function as shown by WMFT scores. These results suggest that tNAA concentration, a marker for neuronal integrity, in spared ipsilesional motor cortical areas may contribute to UE motor function in individuals in the chronic stage post-stroke. These data support previous work noting decreased grey matter structural density in individuals with chronic stroke, suggesting that grey matter atrophy may contribute to arm motor function (Gauthier et al., 2012). The addition of precentral gyrus thickness to our model only resulted in a small increase in R2, but the overall model lost significance, suggesting that the two measures (tNAA concentration and cortical thickness) may beredundant.
Contrary to our hypothesis, we observed no association between Glx and UE motor function. Our original hypothesis that ipsilesional Glxconcentrations would be associated with UE motor performance was based on: 1) reported intracortical excitability differences after stroke (Carmichael, 2012; Tarkka, Kononen, Pitkanen, Sivenius, & Mervaalat, 2008), and 2) past findings showing that Glx concentration in motor areas correlated with arm motor impairment in individuals with chronic stroke (Cirstea et al., 2011). Further, Glx concentration in the motor cortex has been correlated with global motor cortical excitability (assessed by transcranial magnetic stimulation) (Stagg et al., 2011), which was a predictor of motor recovery after stroke (Jung, Kim, Kim, & Paik, 2012). Further investigation of the relationship between glutamate and glutamine concentration and motor function in individuals after stroke will be required to resolve these issues.
Limitations
There were limitations to this study. First, our sample size was small, however previous studies using magnetic resonance spectroscopy (Cirstea et al., 2012; Craciunas et al., 2013) and other imaging approaches (Mang et al., 2015; Pavlakis et al., 1999) in individuals with chronic stroke have had comparable sample sizes. Our design was cross-sectional and our findings could be used to provide preliminary evidence that may inform future longitudinal or interventional studies. Larger longitudinal studies should be performed in order to: 1) further examine relationships between lesion volume and location; 2) understand how these measures may be associated with real world arm function after stroke; and 3) determine how physical rehabilitation may be related to metabolites or cortical thickness in chronicstroke.
The mean age of the Stroke group was greater than that of the neurologically healthy Control group. Because normal aging is associated with cerebral tissue loss (Resnick et al., 2003) and changes in NAA and glutamate concentration (Angelie et al., 2001; Kaiser et al., 2005), it is possible that the group differences observed in these measures can be explained by age differences. Yet the age difference between the groups in our current study was relatively small (Stroke group age 68 +/–9.7 years; Control group age 60.4 +/–6.0 years). There was also an unequal distribution of men and women between the Stroke group (12 male, 5 female) and the Control group (4 male, 7 female). This is important to note because some evidence suggests women have poorer functional outcome after ischemic stroke than men (Gibson, 2013). While it is possible this gender imbalance in our sample affected the results, a study published from the Registry of the Canada Stroke Network found that women had only slightly worse functional status at 6 months post stroke, and that there were no striking sex differences in stroke severity, 6-month mortality or quality of life (Kapral et al., 2005). Additionally, pre-stroke hemisphere dominance information was not collected, however, we compared data from the ipsilesional cortex in individuals with stroke to data extracted from the non-dominant cortex in ourControl group.
For H1MRS data collection, we used a technique known as single voxel spectroscopy (SVS), which obtains a signal from a single volume of interest. An alternative technique to SVS, known as chemical shift imaging (CSI) obtains a chemical signal simultaneously from multiple voxels in brain regions of interest. A limitation with SVS is that our voxel size (30 × 22 × 15 mm) extended beyond the hand area of the primary motor cortex, and therefore captured data from the cortex and underlying white matter. We selected SVS because it is the preferred method for obtaining consistent, high quality in vivo spectra (Drost, Riddle, & Clarke, 2002). Further, the integrity of white matter tracts underlying the motor cortex is associated with motor output in chronic stroke (Schulz et al., 2015), and is therefore relevant to the overall objectives of the study.
Last our MRS acquisition did not include saturation bands for suppression of lipid signal from sub-cutaneous regions outside the brain. Since the voxel placement was relatively close to sub-cutaneous regions outside the brain, saturations bands might have reduced contributions to the spectrum from lipids near the voxel. A small number (2-3) of both control and stroke subject spectra included lipid peaks, which overlapped slightly with the NAA spectral peak. However, in all cases the lipid peak was very well fit by LCModel and the NAA fits had Cramer-Rao lower bounds of 3% or less.
Conclusions
These findings suggest a relationship between arm motor function, metabolite concentrations and cortical thickness. Taken together our results advance knowledge regarding the interactions amongst these factors, and suggest biochemical differences may be important markers of post-stroke motor dysfunction.
Footnotes
Acknowledgments
This work was supported by funding from the Canadian Institutes for Health Research (MOP- 106651; PI: L.A.B).
