Finger somatotopy is preserved after tetraplegia but deteriorates over time
Abstract
Previous studies showed reorganised and/or altered activity in the primary sensorimotor cortex after a spinal cord injury (SCI), suggested to reflect abnormal processing. However, little is known about whether somatotopically specific representations can be activated despite reduced or absent afferent hand inputs. In this observational study, we used functional MRI and a (attempted) finger movement task in tetraplegic patients to characterise the somatotopic hand layout in primary somatosensory cortex. We further used structural MRI to assess spared spinal tissue bridges. We found that somatotopic hand representations can be activated through attempted finger movements in the absence of sensory and motor hand functioning, and no spared spinal tissue bridges. Such preserved hand somatotopy could be exploited by rehabilitation approaches that aim to establish new hand-brain functional connections after SCI (e.g. neuroprosthetics). However, over years since SCI the hand representation somatotopy deteriorated, suggesting that somatotopic hand representations are more easily targeted within the first years after SCI.
Introduction
A spinal cord injury (SCI) refers to damage of the spinal cord caused by a trauma, disease, or degeneration of vertebral discs. SCI mostly results in a partial or complete loss of motor control and sensation. Following a tetraplegia (or cervical SCI), individuals mostly experience a loss of muscle function and sensation in their limbs and torso (Curt et al., 1998; Kalsi-Ryan et al., 2014; Petersen et al., 2012). Accordingly, the primary somatosensory cortex (S1) mostly receives weakened, or is fully deprived of, afferent inputs and is exposed to altered motor behaviour (Ozdemir and Perez, 2018).
Research in non-human primate models of chronic and complete cervical SCI has shown that the S1 hand area becomes largely unresponsive to tactile hand stimulation after the injury (Jain et al., 2008; Kambi et al., 2014; Liao et al., 2021). The surviving finger-related activity became disorganised such that a few somatotopically appropriate sites but also other somatotopically non-matched sites were activated (Liao et al., 2021). Seminal non-human primate research has further demonstrated that SCI leads to extensive cortical reorganisation in S1, such that tactile stimulation of cortically adjacent body parts (e.g. of the face) activated the deprived brain territory (e.g. of the hand; Halder et al., 2018; Jain et al., 2008; Kambi et al., 2014). Although the physiological hand representation appears to largely be altered following a chronic cervical SCI in non-human primates, the anatomical isomorphs of individual fingers are unchanged (Jain et al., 1998). This suggests that while a hand representation can no longer be activated through tactile stimulation after the disruption of afferent spinal pathways, a latent and somatotopic hand representation may be preserved regardless of large-scale physiological reorganisation.
A similar pattern of results has been reported for human SCI patients. Transcranial magnetic stimulation (TMS) studies induced current in localised areas of SCI patients’ primary motor cortex (M1) to induce a peripheral muscle response. They found that representations of more impaired muscles retract or are absent while representations of less impaired muscles shift and expand (Fassett et al., 2018; Freund et al., 2011a; Levy et al., 1990; Streletz et al., 1995; Topka et al., 1991; Urbin et al., 2019). Similarly, human fMRI studies have shown that cortically neighbouring body part representations can shift towards, though do not invade, the deprived M1 and S1 cortex (Freund et al., 2011b; Henderson et al., 2011; Jutzeler et al., 2015; Wrigley et al., 2018; Wrigley et al., 2009). Other human fMRI studies hint at the possibility of latent somatotopic hand representations following SCI by showing that attempted movements with the paralysed and sensory-deprived body part can still evoke signals in the sensorimotor system (Cramer et al., 2005; Freund et al., 2011b; Kokotilo et al., 2009; Solstrand Dahlberg et al., 2018). This attempted ‘net’ movement activity was, however, shown to substantially differ from healthy controls: activity levels have been shown to be increased (Freund et al., 2011b; Kokotilo et al., 2009; Solstrand Dahlberg et al., 2018) or decreased (Hotz-Boendermaker et al., 2008), volumes of activation have been shown to be reduced (Cramer et al., 2005; Hotz-Boendermaker et al., 2008), activation was found in somatotopically non-matched cortical sites (Freund et al., 2011b), and activation was poorly modulated when patients switched from attempted to imagined movements (Cramer et al., 2005). These observations have therefore mostly been attributed to abnormal and/or disorganised processing induced by the SCI. It remains possible though that, despite certain aspects of sensorimotor activity being altered after SCI, somatotopically typical representations of the paralysed and sensory deprived body parts can be preserved (e.g. finger somatotopy of the affected hand(s)). Such preserved representations have the potential to be exploited in a functionally meaningful manner (e.g. via neuroprosthetics).
Case studies using intracortical stimulation in the S1 hand area to elicit finger sensations in tetraplegic patients hint at such preserved somatotopic representations (Fifer et al., 2020; Flesher et al., 2016), with one exception (Armenta Salas et al., 2018). Negative results were suggested to be due to a loss of hand somatotopy and/or reorganisation in S1 of the implanted tetraplegic patient, or due to potential misplacement of the implant (Armenta Salas et al., 2018). Whether fine-grained somatotopy is generally preserved in the tetraplegic patient population remains unknown. It is also unclear what clinical, behavioural, and structural spinal cord determinants may influence such representations to be maintained.
Here we used functional MRI (fMRI) and a visually cued (attempted) finger movement task in tetraplegic patients to examine whether hand somatotopy is preserved following a disconnect between the brain and the periphery. We instructed patients to perform the fMRI tasks with their most impaired upper limb and matched control participants’ tested hands to patients’ tested hands. If a patient was unable to make overt finger movements due to their injury, then we carefully instructed them to make attempted (i.e. not imagined) finger movements. To see whether patient’s maps exhibited characteristics of somatotopy, we visualised finger selectivity in S1 using a travelling wave approach. To investigate whether fine-grained hand somatotopy was preserved and could be activated in S1 following tetraplegia, we assessed inter-finger representational distance patterns using representational similarity analysis (RSA). These inter-finger distance patterns are thought to be shaped by daily life experience such that fingers used more frequently together in daily life have lower representational distances (Ejaz et al., 2015). RSA-based inter-finger distance patterns have been shown to depict the invariant representational structure of fingers in S1 and M1 better than the size, shape, and exact location of the areas activated by finger movements (Ejaz et al., 2015). Over the past years RSA has therefore regularly been used to investigate somatotopy of finger representations both in healthy (e.g. Akselrod et al., 2017; Ariani et al., 2020; Ejaz et al., 2015; Gooijers et al., 2021; Kieliba et al., 2021; Kolasinski et al., 2016b; Liu et al., 2021; Sanders et al., 2019) and patient populations (e.g. Dempsey-Jones et al., 2019; Ejaz et al., 2016; Kikkert et al., 2016; Wesselink et al., 2019). We closely followed procedures that have previously been used to map preserved and typical somatotopic finger selectivity and inter-finger representational distance patterns of amputees’ missing hands in S1 using volitional phantom finger movements (Kikkert et al., 2016; Wesselink et al., 2019). However, in amputees, these movements generally recruit the residual arm muscles that used to control the missing limb via intact connections between the brain and spinal cord. Whether similar preserved somatotopic mapping can be observed in SCI patients with diminished or no connections between the brain and the periphery is unclear. If finger somatotopy is preserved in tetraplegic patients, then we should find typical inter-finger representational distance patterns in the S1 hand area of these patients. By measuring a group of 14 chronic tetraplegic patients with varying amounts of spared spinal cord tissue at the lesion level (quantified by means of midsagittal tissue bridges based on sagittal T2w scans), we uniquely assessed whether preserved connections between the brain and periphery are necessary to preserve fine somatotopic mapping in S1 (Huber et al., 2017; Pfyffer et al., 2019). If spared connections between the periphery and the brain are not necessary for preserving hand somatotopy, then we should find typical inter-finger representational distance patterns even in patients without spared spinal tissue bridges. We also investigated what clinical and behavioural determinants may contribute to preserving S1 hand somatotopy after chronic tetraplegia. If spared sensorimotor hand function is not necessary for preserving hand somatotopy, then we should find typical inter-finger representational distance patterns even in patients who suffer from full sensory loss and paralysis of the hand(s).
Results
Patient impairments
We tested 14 chronic tetraplegic patients that were heterogenous in terms of completeness of the SCI (ranging from AIS-A to AIS-D), neurological level of the injury (ranging from C2 to C7), years since injury (ranging from 6 months to 33 years since SCI), and sensorimotor upper limb impairments (ranging from a Graded Redefined Assessment of Strength, Sensibility and Prehension test (GRASSP) score for both limbs of 21 to 220; healthy GRASSP score = 232). We further tested 18 age-, sex-, and handedness-matched healthy control participants.
Finger selectivity is preserved following tetraplegia
We used 3T fMRI and a travelling wave paradigm to investigate the somatotopic layout of finger selectivity on the S1 cortical surface (Besle et al., 2013; Kolasinski et al., 2016a). To visualise whether there was a consistent and somatotopic layout of finger selectivity across participants in the control and SCI patient groups, we created probability maps of finger selectivity. A characteristic hand map shows a gradient of finger preference, progressing from the thumb (laterally) to the little finger (medially). We found a characteristic progression of finger selective clusters in both the control and tetraplegic patient groups (Figure 1). Qualitative inspection suggests that inter-participant consistency was lowest for the little finger representation in both groups. It further appears that overall inter-participant consistency was reduced in the patient group compared to the control group.
Given the known the relatively high inter-participant variability in finger selectivity (Kolasinski et al., 2016a) and the clinical heterogeneity of our SCI patient group, we further visualised each individual participant’s finger selectivity map (Figure 2A and B). Overall, we found aspects of somatotopic finger selectivity in the maps of SCI patients’ hands, in which neighbouring clusters showed selectivity for neighbouring fingers in contralateral S1, similar to those observed in 18 age-, sex-, and handedness-matched healthy controls. Notably, a characteristic hand map was even found in a patient who suffered complete paralysis and sensory deprivation of the hands (Figure 2B, patient map 1; patient S01). Despite most maps (Figure 2, except patient map 3; patient S04) displaying aspects of characteristic finger selectivity, some finger representations were not visible in the thresholded patient and control maps.
To ensure that the observed finger selective clusters were not representing noise, but rather true finger selectivity, we calculated split-half consistency between two halves of the travelling wave dataset using the Dice overlap coefficient (DOC; Dice, 1945). Minimally thresholded finger maps were compared across the split-halves of the data within an S1 mask. Overall, split-half consistency was not significantly different between patients and controls, as tested using a robust mixed ANOVA (see Figure 2C; F(1,20.32) = 0.51, p=0.48). There was a significant difference in split-half consistency between pairs of same, neighbouring, and non-neighbouring fingers (F(2,18.59) = 159.69, p<0.001). This neighbourhood relationship was not significantly different between the control and patient groups (i.e. there was no significant interaction; F(2,18.59) = 2.44, p=0.11).
The DOC was highest for comparison of the same fingers between two halves of the dataset compared to neighbouring (controls: W = 171, p<0.001, BF10 = 2761.24; patients: W = 91, p<0.001, BF10 = 742.45) and non-neighbouring fingers (controls: W = 171, p<0.001, BF10 = 248.14; patients: W = 91, p<0.001, BF10 = 133.05). Moreover, neighbouring fingers showed greater overlap across the split-halves of the dataset than non-neighbouring fingers (controls: W = 171, p<0.001, BF10 = 835.29; patients: W = 91, p=0.001, BF10 = 65.32). This demonstrates that there was a somatotopic gradient in split-half consistency that was similar between the control and patient groups, demonstrating that the finger maps represented true finger selectivity.
Typical hand somatotopy is preserved following tetraplegia
Next, we assessed 3T univariate task-related activity during individual finger movements performed in a blocked design fashion. Task-related activity was quantified by extracting the percent signal change for finger movement (across all fingers) versus baseline within the contralateral S1 hand area (see Figure 3A). Overall, all patients were able to engage their S1 hand area by moving individual fingers (t(13) = 7.46, p<0.001; BF10 = 4.28 e + 3), as did controls (t(17) = 9.92, p<0.001; BF10 = 7.40 e + 5). Furthermore, patients’ task-related activity was not significantly different from controls (t(30) = –0.82, p=0.42; BF10 = 0.44), with the Bayes factor (BF) showing anecdotal evidence in favour of the null hypothesis. Similar results were found when exploring univariate task-related activity in the contralateral M1 hand ROI (see Figure 3—figure supplement 1).
While the travelling wave maps demonstrate finger selectivity, they provide little information about the overlap between finger representations. We therefore examined the intricate relationship between finger representations in the S1 hand area for all patients and controls using RSA (see Figure 3B and C). The resulting inter-finger distances were averaged across finger pairs within each participant to obtain an estimate for average inter-finger separability (see Figure 3D). We found that inter-finger separability in the S1 hand area was greater than 0 for patients (t(13) = 9.83, p<0.001; BF10 = 6.77 e + 4) and controls (t(17) = 11.70, p<0.001; BF10 = 6.92 e + 6), indicating that the S1 hand area in both groups contained information about individuated finger representations. Furthermore, for both controls (W = 171, p<0.001; BF10 = 4059) and patients (W = 105, p<0.001; BF10 = 279) there was significant greater separability (or representation strength) in the S1 hand area than in a control cerebral spinal fluid (CSF) ROI that would not be expected to contain information about individuated finger representations. We did not find a significant group difference in inter-finger separability of the S1 hand area (t(30) = 1.52, p=0.14; BF10 = 0.81), with the BF showing anecdotal evidence in favour of the null hypothesis.
We then tested whether the inter-finger distances were different across finger pairs between controls and tetraplegic patients using a robust mixed ANOVA with a within-participants factor for finger pair (10 levels) and a between-participants factor for group (two levels: controls and tetraplegic patients; Figure 3—figure supplement 2). We did not find a significant difference in inter-finger distances between patients and controls (F(1,21.66) = 1.50, p=0.23). The inter-finger distances were significantly different across finger pairs, as would be expected based on somatotopic mapping (F(9,15.38) = 27.22, p<0.001). This pattern of inter-finger distances was not significantly different between groups (i.e. no significant finger pair by group interaction; F(9,15.38) = 1.05, p=0.45). When testing for group differences per finger pair, the BF only revealed inconclusive evidence (BF >0.37 and < 1.11; note that we could not run a Bayesian ANOVA due to normality violations).
Although inter-finger separability was not significantly different between patients and controls, it is possible that the pattern of inter-finger distances was atypical in the patients. We therefore examined whether the inter-finger distance pattern was normal (or typical) in tetraplegic patients (see Figure 3E) by correlating each participant’s inter-finger distance pattern with a canonical inter-finger distance pattern. Tetraplegic patients’ typicality scores were compared to those of the controls and of a group congenital one-handers (data taken from an independent study; Wesselink et al., 2019). Congenital one-handers are born without a hand and therefore do not have a cortical representation of the missing hand (unlike amputees who develop the representation before losing the hand and hence have a ‘missing hand representation’; Wesselink et al., 2019). This group was therefore included as a control for absence of hand representation. We found a significant difference in typicality between tetraplegic patients, healthy controls, and congenital one-handers (H(2) = 26.64, p<0.001). As expected, post hoc tests revealed significantly higher typicality in controls compared to congenital one-handers (U = 0, p<0.001; BF10 = 113.60). Importantly, inter-finger distances typicality of the SCI patients was significantly higher than the typicality scores of the congenital one-handers (U = 4, p<0.001; BF10 = 90.33), but not significantly different from the typicality scores of the controls (U = 103, p=0.40; BF10 = 0.55). The BF for the comparison between SCI patients and controls showed anecdotal evidence for equivalence between both groups.
Typical hand somatotopy deteriorates over years after tetraplegia
Next, we aimed to understand which clinical, behavioural, and structural spinal cord determinants may allow hand representations in S1 to be maintained. We first explored correlations with S1 hand representation typicality. We found that the number of years since SCI significantly correlated with hand representation typicality (see Figure 4A; rs = –0.59, p=0.028), suggesting that S1 hand representation typicality may deteriorate over time after SCI. We further found that patients with more retained GRASSP motor function of the tested upper limb had more typical hand representations in S1 (see Figure 4B; rs = 0.60, p=0.02). We did not find a significant correlation between S1 hand representation typicality and GRASSP sensory function of the tested upper limb, spared midsagittal spinal tissue bridges at the lesion level, or cross-sectional spinal cord area (see Figure 4C–E; rs = 0.40, p=0.16, rs = 0.46, p=0.14, and rs = 0.33, p=0.25, respectively).
We further explored the hand representation typicality of patients S01 and S03 who did not have any spared midsagittal spinal tissue bridges at the lesion level, a complete (S01) or near complete (S03) hand paralysis, and a complete (S01) or near complete loss (S03) of hand sensory function (as assessed using the GRASSP test). Interestingly, both patients had a highly typical hand representation in S1 that was significantly different from congenital one-handers (i.e. who are born without a hand and do not have a missing hand representation; S01: t(12) = 3.20, p=0.008; S03: t(12) = 2.97, P = 0.01), but not controls (S01: t(17) = 0.95, p=0.36; S03: t(17) = 0.04, p=0.97). This suggests that retained connections between the periphery and the brain, retained motor functioning, and retained sensory functioning may not be necessary to maintain typical hand representations in S1.
We then ran an exploratory stepwise linear regression to investigate which of these clinical, behavioural, and structural spinal cord characteristics were predictive of hand representation typicality in S1. Years since SCI significantly predicted hand representation typicality in S1 with R2 = 0.40 (F(1,10) = 6.73, p=0.027). Motor function of the tested upper limb, sensory function of the tested hand, spared midsagittal spinal tissue bridges at the lesion level, and spinal cord area did not significantly add to the prediction (t = 1.43, p=0.19, t = 1.44, p=0.18, t = 1.19, p=0.26, and t = 0.41, p=0.69, respectively). This analysis suggests that while hand representations are preserved following tetraplegia, they may deteriorate over time.
To inspect this further, we bootstrapped the mean typicality of the able-bodied controls and congenital one-handers 10,000 times to infer their population means. While most tetraplegic patients’ typicality scores fell in, or very close to, the able-bodied controls’ distribution, we found that some SCI patients’ typicality scores fell in-between the able-bodied controls’ and congenital one-handers’ distributions (Figure 4F). This suggests that these patients may not be similar to congenital one-handers or to able-bodied controls. Interestingly, this included those patients for whom most years had passed since their SCI. This suggests that S1 hand representations might deteriorate after an SCI, but some weak hand information may be maintained in S1 even >30 years after an SCI.
Discussion
In this study, we investigated whether hand somatotopy is preserved and can be activated through attempted movements following tetraplegia. We tested a heterogenous group of SCI patients to examine what clinical, behavioural, and structural spinal cord determinants contribute to preserving S1 somatotopy. Our results revealed that detailed hand somatotopy can be preserved following tetraplegia, even in the absence of sensory and motor hand function and a lack of spared spinal tissue bridges. However, over time since SCI these finger maps deteriorated such that the hand somatotopy became less typical.
Spared spinal cord tissue bridges can be found in most patients with a clinically incomplete injury, their width being predictive of electrophysiological information flow, recovery of sensorimotor function, and neuropathic pain (Huber et al., 2017; Pfyffer et al., 2021; Pfyffer et al., 2019; Vallotton et al., 2019). However, in this study, spared midsagittal spinal tissue bridges at the lesion level and sensorimotor hand function did not seem necessary to maintain and activate a somatotopic hand representation in S1. We found a highly typical hand representation in two patients (S01 and S03) who did not have any spared spinal tissue bridges at the lesion level, a complete (S01) or near complete (S03) hand paralysis, and a complete (S01) or near complete loss (S03) of hand sensory function. Our predictive modelling results were in line with this notion and showed that these behavioural and structural spinal cord determinants were not predictive of hand representation typicality. Note however that our sample size was limited, and it is challenging to draw definite conclusions from non-significant predictive modelling results.
Time since injury was predictive of a deteriorated, or less typical, somatotopic S1 hand representation. In fact, both patients with a typical S1 hand representations but absent spinal tissue bridges suffered their SCI only 4 years ago (the group was on average 12 years since SCI). The hand representation typicality of patients who suffered their SCI further in the past were not similar to congenital one-handers’ nor to controls’ hand representation. Thus, S1 hand representations may deteriorate over time after an SCI, but some weak hand information appears to be maintained in S1 even >30 years after an SCI. This finding complements previous studies in amputees showing that missing hand somatotopy is preserved even decades after arm amputation and that years since injury was not related to missing hand representation typicality (Kikkert et al., 2016; Wesselink et al., 2019). While both amputees and SCI patients suffer from major sensory input loss and changed motor behaviour, their injuries are inherently different. The injured axons within amputees’ residual limb and the remaining part of the peripheral nerves mostly generate some spontaneous (ectopic) activity that is propagated to the brain and could contribute to maintaining S1 representations (Kikkert et al., 2019; Kikkert et al., 2018; Kikkert et al., 2016; Makin et al., 2013; Nyström and Hagbarth, 1981; Vaso et al., 2014). Furthermore, amputees mostly remain able to move and receive afferents from the residual limb muscles that used to control the missing hand as most of their motor system remains intact. Although their hand is missing, amputees’ residual arm muscles are often still used (either to move the residual limb and/or to control a prosthetic arm). Lastly, amputees’ vividness of kinaesthetic sensations during phantom finger movements was found to be predictive of the typicality of the S1 missing hand representation (Wesselink et al., 2019). A continued, though altered, experience relating to the missing hand in amputees may contribute to maintaining the somatotopic missing hand representation in S1. Contrarily, tetraplegic patients mostly have reduced or a complete loss of communication between the brain and periphery. They therefore have problems activating the adequate muscles and will lose orderly afferents from their muscles and the skin. It is possible that this continued disuse causes somatotopic S1 representations to deteriorate after SCI.
How may these representations be preserved over time and activated through attempted movements in the absence of peripheral information? S1 is reciprocally connected with various brain areas, for example,, M1, lateral parietal cortex, posterior parietal area 5, secondary somatosensory cortex, and supplementary motor cortex (Delhaye et al., 2018). After a loss of sensory inputs and paralysis through SCI, S1 representations may be activated and preserved through its interconnections with these areas. Firstly, it is possible that cortico-cortical efference copies may keep a representation ‘alive’ through occasional corollary discharge (London and Miller, 2013). While motor and sensory signals no longer pass through the spinal cord in the absence of spinal tissue bridges, S1 and M1 remain intact. When a motor command is initiated (e.g. in the form of an attempted hand movement), an efference copy is thought to be sent to S1 in the form of corollary discharge. This corollary discharge resembles the expected somatosensory feedback activity pattern and may drive somatotopic S1 activity even in the absence of ascending afferent signals from the hand (Adams et al., 2013; London and Miller, 2013). It is possible that our patients occasionally performed attempted movements which would result in corollary discharge in S1. Second, it is likely that attempting individual finger movements poses high attentional demands on tetraplegic patients. Accordingly, attentional processes might have contributed to eliciting somatotopic S1 activity. Evidence for this account comes from studies showing that it is possible to activate somatotopic S1 hand representations through attending to individual fingers (Puckett et al., 2017) or through touch observation (Kuehn et al., 2018). Attending to fingers during our attempted finger movement task may have been sufficient to elicit somatotopic S1 activity through top-down processes in the tetraplegic patients who lacked motor and sensory hand function. Furthermore, one might speculate that observing others’ or one’s own fingers being touched or directing attention to others’ hand movements or one’s own fingers may help preserve somatotopic representations. Third, it is possible that these somatotopic maps are relatively hardwired, and while they deteriorate over time, they never fully disappear. Indeed, somatotopic mapping of a sensory-deprived body part has been shown to be resilient after dystonia (Ejaz et al., 2016; though see Burman et al., 2008 and Elbert et al., 1998) and arm amputation (Bruurmijn et al., 2017; Kikkert et al., 2016; Wesselink et al., 2019). Fourth, it is possible that even though a patient is clinically assessed to be complete and is unable to perceive sensory stimuli on the deprived body part, there is still some ascending information flow that contributes to preserving somatotopy (Wrigley et al., 2018). A recent study found that although complete paraplegic SCI patients were unable to perceive a brushing stimulus on their toe, 48% of patients activated the location-appropriate S1 area (Wrigley et al., 2018). However, the authors of this study defined the completeness of patients’ injuries via behavioural testing, while we additionally assessed the retained connections passing through the SCI directly via quantification of spared spinal tissue bridges through structural MRI. It is unlikely that spinal tissue carrying somatotopically organised information would be missed by our assessment (Huber et al., 2017; Pfyffer et al., 2019). Our experiment did not allow us to tease apart these potential processes, and it is likely that various processes simultaneously influence the preservation of S1 somatotopy and elicited the observed somatotopic S1 activity.
Our finding of preserved S1 somatotopy may appear inconsistent with the wealth of evidence showing cortical reorganisation in S1 following SCI (Halder et al., 2018; Jain et al., 2008; Kambi et al., 2014). In these studies, experimenters indirectly probed the deprived S1 hand cortex via stimulation of cortically adjacent body parts. Human fMRI studies similarly probed the intact and cortically neighbouring body parts and suggested that their representations shift towards the deprived S1 cortex (Freund et al., 2011b; Freund et al., 2011a; Henderson et al., 2011; Jutzeler et al., 2015; Wrigley et al., 2018; Wrigley et al., 2009). TMS studies similarly induce current in localised areas of M1 to induce a peripheral muscle response of cortically neighbouring body parts. These studies demonstrated that the representations of less impaired muscles shift and expand following a complete or incomplete SCI (Fassett et al., 2018; Freund et al., 2011a; Levy et al., 1990; Streletz et al., 1995; Topka et al., 1991; Urbin et al., 2019). Our fMRI results showed that tetraplegic patients had a preserved somatotopic hand representation in S1, though this deteriorated over time. We did not probe body parts other than the hand and could therefore not investigate whether any remapping of other (neighbouring and/or intact) body part representations towards or into the deprived S1 hand cortex may have taken place. Whether reorganisation and preservation of the original function can simultaneously take place within the same cortical area therefore remains a topic for further investigation. It is possible that reorganisation and preservation of the original function could co-occur within cortical areas. Indeed, non-human primate studies demonstrated that remapping observed in S1 actually reflects reorganisation in subcortical areas of the somatosensory pathway, principally the brainstem (Chand and Jain, 2015; Kambi et al., 2014). As such, the deprived S1 area receives reorganised somatosensory inputs upon tactile stimulation of neighbouring intact body parts. This would simultaneously allow the original S1 representation of the deprived body part to be preserved, as observed in our results when we directly probed the deprived S1 hand area through attempted finger movements.
Together, our findings indicate that in the first years after a tetraplegia the somatotopic S1 hand representation is preserved and can be activated through attempted movements even in the absence of retained sensory hand function, motor hand function, and spared spinal tissue bridges. These preserved S1 finger maps could be exploited in a functionally meaningful manner by rehabilitation approaches that aim to establish new functional connections between the brain and the hand after a tetraplegia, for example, through neuroprosthetic limbs or advanced exoskeletons that are directly controlled by the brain (Ajiboye et al., 2017; Armenta Salas et al., 2018; Bouton et al., 2016; Lebedev and Nicolelis, 2017).
Materials and methods
Participants
15 chronic (i.e. > 6 months post injury) tetraplegic patients were recruited and 14 patients completed the measurements (mean age ± s.e.m. = 55 ± 3.6 years ; one female; six dominant left-handers; see Table 1 for demographic and clinical details). Patient inclusion criteria were as follows: aged 18–75 years, no MRI contraindications, at least 6 months post SCI, no neurological impairment or body function impairments not induced by SCI, and able to provide informed consent. 18 age-, sex-, and handedness-matched able-bodied control participants (age = 56 ± 3.6 years; one female; five dominant left-handers) also participated in this study. Control participant inclusion criteria were as follows: aged 18–75 years, no MRI contraindications, no impairment of body function induced by SCI, no neurological illness, no hand impairments, and able to provide informed consent.
Participants’ informed consent was obtained according to the Declaration of Helsinki prior to study onset. Ethical approval was granted by the Kantonale Ethikkommission Zürich (KEK-2018-00937). This study is registered on clinicaltrials.gov under NCT03772548. Two patients and one control participant were scanned twice due to excessive head motion during fMRI acquisition or suboptimal slice placement. One patient withdrew from the study prior to study completion. Data of one control participant were distorted and not usable for analysis. All data relating to these participants were discarded from all data analysis. Patient S02 did not complete the travelling wave measurements due to time constraints.
Clinical characterisation
Request a detailed protocolBehavioural testing was conducted in a separate session. We used the International Standards for Neurological Classification of Spinal Cord Injury (ISNCSCI) to neurologically classify patients’ completeness of injury and impairment level. We used the GRASSP assessment to define sensory and motor integrity of the upper limbs (Kalsi-Ryan et al., 2012). Each upper limb’s maximum score is 116 and refers to healthy conditions. We determined each patient’s most impaired upper limb according to the GRASSP. Note that GRASSP motor scores reflect overall upper-limb motor function (i.e. including arm and shoulder functioning; see Table 2 for muscle-specific GRASSP scores). GRASSP sensory scores are hand specific.
fMRI tasks
Request a detailed protocolWe employed two separate paradigms to uncover fine-grained somatotopic hand representations using fMRI: first, we used a travelling wave paradigm to investigate the somatotopic hand layout on the S1 cortical surface (Besle et al., 2013; Kolasinski et al., 2016a). Second, we employed a blocked design and RSA that takes into account the entire fine-grained activity pattern of each finger (i.e. including the representational inter-finger relationships; Ejaz et al., 2015; Kriegeskorte et al., 2008).
Participants were visually cued to perform individual finger movements while their palm was positioned up. Patients were instructed to perform the fMRI tasks with their most impaired upper limb (identified using the GRASSP, see previous section). Controls’ tested hands were matched to the patients. Due to their injury, not all patients were able to make overt finger movements. In these cases, patients were carefully instructed by the experimenter to make attempted (i.e. not imagined) finger movements. The experimenter explained that although it is not possible for the patient to perform an overt movement, attempting to perform the movement may still produce a motor command in the brain. Importantly, complete paraplegic patients are able to distinguish between attempted and imagined movements with their paralysed body part (Cramer et al., 2005; Hotz-Boendermaker et al., 2008; Sabbah et al., 2002). Furthermore, attempted foot movements activated SCI patients’ motor network similarly to controls performing overt foot movements (Cramer et al., 2005; Hotz-Boendermaker et al., 2008; Sabbah et al., 2002).
Participants saw five horizontally aligned white circles, corresponding to the five fingers, via a visual display viewed through a mirror mounted on the head coil. For participants moving their left hand, the leftmost and rightmost circles corresponded to the thumb and little finger, respectively. For participants moving their right hand, the leftmost circle corresponded to the little finger and the rightmost circle to the thumb. To cue a finger movement, the circle corresponding to this finger turned red. Participants performed self-paced flexion/extension movements with the cued finger for the duration of the colour change. Instructions were delivered using Psychtoolbox (v3) implemented in MATLAB (v2014). Head motion was minimised using over-ear MRI-safe headphones or padded cushions.
The travelling wave paradigm involved individuated finger movements in a set sequence. Each 10 s finger movement block was immediately followed by a movement block of a neighbouring finger. The forward sequence cycled through the fingers: thumb-index-middle-ring-little. To account for order-related biases due to the set movement cycle and sluggish haemodynamic response, we also collected data using a backward sequence: the backward sequence cycled through the movements in a reverse of the forward sequence: little-ring-middle-index-thumb fingers. The forward and backward sequences were employed in separate runs. A run lasted 6 min and 4 s, during which a sequence was repeated seven times. The forward and backward runs were repeated twice, with a total duration of 24 min and 16 s.
The blocked design consisted of six conditions: movement conditions for each of the five fingers and a rest condition. Finger movement instructions were as described above, and the word ‘Rest’ indicated the rest condition. A movement block lasted 8 s, and each condition was repeated five times per run in a counterbalanced order. Each run comprised a different block order and had a duration of 4 min and 14 s. We acquired four runs, with a total duration of 16 min and 56 s.
MRI acquisition
Request a detailed protocolMRI data were acquired using a Philips 3 Tesla Ingenia system (Best, The Netherlands) with a 17-channel HeadNeckSpine or, in case of participant discomfort due to the coil’s narrowness, a 15-channel HeadSpine coil. Anatomical T1-weighted images covering the brain and cervical spinal cord were acquired using the following acquisition parameters: 0.8 mm3 resolution, repetition time (TR) = 9.3 ms, echo time (TE) = 4.4 ms, and flip angle 8°. Anatomical T2-weighted images of the cervical spinal cord were acquired sagittally using the following acquisition parameters: 1 × 1 × 3 mm resolution, TR = 4500 ms, TE = 85 ms, flip angle = 90°, and slice gap = 0.3 mm, 15 slices. Task-fMRI data were acquired using an echo-planar-imaging (EPI) sequence with partial brain coverage: 22 sagittal slices were centred on the anatomical location of the hand knob with coverage over the thalamus and brainstem. We used the following acquisition parameters: 2 mm3 resolution, TR = 2000 ms, TE = 30 ms, flip angle = 82°, and SENSE factor = 2.2. We acquired 182 and 127 volumes for each of the travelling wave and blocked design runs, respectively.
fMRI analysis
Request a detailed protocolfMRI analysis was implemented using FSL v6.0 (https://fsl.fmrib.ox.ac.uk/fsl/fslwiki), Advanced Normalization Tools (ANTs) v2.3.1 (http://stnava.github.io/ANTs), the RSA toolbox (Nili et al., 2014; Wesselink and Maimon-Mor, 2017), and MATLAB (R2018a). Cortical surface visualisations were realised using FreeSurfer (https://surfer.nmr.mgh.harvard.edu/; Dale et al., 1999; Fischl et al., 2001) and Connectome Workbench (https://www.humanconnectome.org/software/connectome-workbench).
fMRI preprocessing
Request a detailed protocolCommon preprocessing steps were applied using FSL’s Expert Analysis Tool (FEAT). The following preprocessing steps were included: motion correction using MCFLIRT (Jenkinson et al., 2002), brain extraction using automated brain extraction tool BET (Smith, 2002), spatial smoothing using a 2 mm full-width-at-half-maximum (FWHM) Gaussian kernel, and high-pass temporal filtering with a 100 s (blocked design runs) or 90 s (travelling wave runs) cut-off.
Image registration
Request a detailed protocolImage co-registration was done in separate, visually inspected, steps. For each participant, a midspace was calculated between the four blocked design runs, that is, an average space in which images are minimally reoriented. We then transformed all fMRI data to this midspace using purely rigid probability mapping in ANTs. Next, we registered each participant’s midspace to the T1-weighted image, initially using 6 degrees of freedom and the mutual information cost function, and then optimised using boundary-based registration (BBR; Greve and Fischl, 2009) Each co-registration step was visually inspected and, if needed, manually optimised using blink comparison in Freeview.
Travelling wave analysis
Request a detailed protocolThe travelling wave approach is characterised by set finger movement cycles that are expected to result in neighbouring cortical activations. It is designed to capture voxels that show preferential activity to one condition, above and beyond all other conditions (i.e. winner-takes-all principle; testing for finger selectivity). The travelling wave approach is especially powerful to reveal the smooth progression of neighbouring representations that are specific for topographic maps. This technique is therefore frequently used to uncover retinotopic (DeYoe et al., 1996; Engel et al., 1997; Sereno et al., 1995), somatotopic (Besle et al., 2013; Kikkert et al., 2016; Kolasinski et al., 2016a; Mancini et al., 2012; Zeharia et al., 2015), and tonotopic representations (Da Costa et al., 2015; Talavage et al., 2004). Importantly, S1 finger movement somatotopy assessed using the travelling wave approach is highly consistent across scanning sessions (Kolasinski et al., 2016a).
Travelling wave analysis was conducted separately for each participant, closely following procedures previously described in Kikkert et al., 2016. A reference model was created using a gamma haemodynamic response function (HRF) convolved boxcar, using a 10 s ‘on’ (a single finger movement duration) and 40 s ‘off’ period (movement duration of all other fingers). This reference model was then shifted in time to model activity throughout the full movement cycle. Since we had a 2 s TR and a 50 s movement cycle, the reference model was shifted 25 times.
Within each individual run, each voxel’s preprocessed BOLD signal time course was cross-correlated with each of the 25 reference models. This resulted in 25 r-values per voxel per run that were normalised using a Fisher r-to-z transformation. To create finger-specific (i.e. hard-edged) maps, we assigned individual lags to specific fingers and averaged the z-values across these lags. This resulted in five z-values, one for each finger per voxel per run. These z-values were averaged across runs per voxel and per finger assignment. As a result, each individual voxel now contained five averaged z-values, one for each finger. Next, a winner-take-all approach was used to assign each voxel to one finger exclusively based on the maximum z-value, providing us with finger specificity.
To visualise the smooth gradient of progression across fingers, we produced lag-specific maps. The backward runs’ standardised cross-correlation z-values were lag-reversed and averaged with the forward runs per voxel and per lag. This resulted in 25 z-values per voxel (one per lag). Next, we used a winner-take-all principle to find the maximum z-value across lags for each voxel, providing us with lag specificity.
Cortical surface projections were constructed from participant’s T1-weighted images. The winner-take-all finger-specific and lag-specific gradient maps were projected onto the cortical surface using cortical-ribbon mapping. Thresholding was applied to the winner-take-all finger-specific maps and lag-specific maps on the cortical surface using a false discovery criterion q < 0.05 based on the native (3D) values. The false discovery rate (FDR) thresholded finger-specific maps were combined to form a hand map. Within this hand map, the lag-specific map was used to visualise the smooth gradient of progression across fingers. We were unable to find a characteristic hand map in one patient who anecdotally reported post hoc that he was unable to performed attempted finger movements.
Inter-participant probability of finger selectivity
Request a detailed protocolTo visualise inter-participant consistency of somatotopic finger-selective representations, we calculated cortical activation probability maps. To ensure that the tested hemisphere was consistently aligned for all participants, we first flipped the acquired T1-weighted images and the travelling wave winner-take-all finger-specific maps along the x-axis for the left-hand tested participants. For these left-hand tested participants, we created new cortical surface projections using their flipped T1-weighted images. Each participant’s cortical surface was then inflated into a sphere and aligned to the FreeSurfer 2D average atlas using sulcal depth and curvature information. The travelling wave winner-take-all finger-specific maps were resampled to the FreeSurfer 2D average atlas and thresholded using a false discovery criterion q < 0.05 based on the native (3D) values. We then calculated finger-specific probability maps for the control and SCI patient groups, separately.
Spatial correspondence of finger maps over time: DOC analysis
Request a detailed protocolTo confirm that the travelling wave finger-specific maps did not represent random noise, we quantified spatial consistency of finger preference between two halves of the data using the Dice overlap coefficient (DOC; Dice, 1945). The DOC calculates the spatial overlap between two representations relative to the total area of these representations. The DOC ranges from 0 (no spatial overlap) to 1 (perfect spatial overlap). If A and B represent the areas of two representations, then the DOC is expressed as
We followed previously described procedures for calculating the DOC between two halves of the travelling wave data (Kikkert et al., 2016; Kolasinski et al., 2016a; Sanders et al., 2019). The averaged finger-specific maps of the first forward and backward runs formed the first data half. The averaged finger-specific maps of the second forward and backward runs formed the second data half. The finger-specific clusters were minimally thresholded (Z > 2) on the cortical surface and masked using an S1 ROI that was created based on Brodmann area parcellation using FreeSurfer. We used minimally thresholded finger-specific clusters for DOC analysis to ensure that we were sensitive to overlaps that would be missed when using high thresholds (see Figure 2—figure supplement 1 for a visualisation of the minimally thresholded split-half hard-edged finger maps used to calculate the DOC). Note that the same results were found when thresholding the finger-specific clusters using an FDR q < 0.05 criterion (see Figure 2—figure supplement 2). The DOC was calculated between same, neighbouring, and non-neighbouring fingers across the two data halves (see Figure 2C).
If the finger maps would be spatially consistent and represent true finger selectivity, then one would expect a higher DOC between pairs of ‘same’ fingers across the two data halves compared to neighbouring and non-neighbouring finger pairs. One would further expect to find a somatotopic relationship in the DOCs: that is, a higher DOC between neighbouring compared to non-neighbouring finger pairs. We tested whether the somatotopic relationship in the DOCs was different in controls and patients using a robust mixed ANOVA with a within-participants factor for finger pair (three levels: same, neighbouring and non-neighbouring finger pairs) and a between-participants factor for group (two levels: controls and SCI patients).
Univariate analysis
Request a detailed protocolTo assess univariate task-related activity of the blocked design data, time-series statistical analysis was carried out per run using FMRIB’s Improved Linear Model (FILM) with local autocorrelation correction, as implemented in FEAT. We obtained activity estimates using a general linear modelling (GLM) based on the double-gamma HRF and its temporal derivative. Each finger movement condition was contrasted with rest. A further contrast was defined for overall task-related activity by contrasting all movement conditions with rest. A fixed effects higher-level analysis was ran for each participant to average across runs.
We defined an S1 hand ROI by converting the S1 ROI used to calculate split-half consistency to volumetric space. Any holes were filled and non-zero voxels were mean dilated. Next, the axial slices spanning 2 cm medial/lateral to the hand knob (Yousry et al., 1997) were identified on the 2 mm MNI standard brain (min-max MNI z-coordinates = 40–62). This mask was non-linearly transformed to each participant’s native structural space. Finally, we used this mask to restrict the S1 ROI and extracted an S1 hand area ROI. The percent signal change for overall task-related activity was then extracted for voxels underlying this S1 hand ROI per participant. A similar analysis was used to investigate overall task-related activity in an M1 hand ROI (see Figure 3—figure supplement 1). We further compared activity levels in finger-specific ROIs in S1 between groups and conducted a geodesic distance analysis to assess whether the finger representations of SCI patients were aligned differently and/or shifted compared to the control group (see Appendix 1).
Representational similarity analysis
Request a detailed protocolWhile the traditional travelling wave approach is powerful to uncover the somatotopic finger arrangement, a fuller description of hand representation can be obtained by taking into account the entire fine-grained activity pattern of all fingers. RSA-based inter-finger overlap patterns have been shown to depict the invariant representational structure of fingers better than the size, shape, and exact location of the areas activated by finger movements (Ejaz et al., 2015). RSA-based measures are furthermore not prone to some of the problems of measurements of finger selectivity (e.g. dependence on map thresholds). We estimated inter-finger overlap using RSA. Note that it is also possible to estimate somatotopic overlap from travelling wave data using an iterated Multigrid Priors (iMGP) method and population-receptive field modelling (Da Rocha Amaral et al., 2020; Puckett et al., 2020). We computed the distance between the activity patterns measured for each finger pair within the S1 hand ROI using the cross-validated squared Mahalanobis distance (or crossnobis distance; Nili et al., 2014). We extracted the blocked design voxel-wise parameter estimates (betas) for each finger movement condition versus rest (identified in the univariate analysis) and the model fit residuals under the S1 hand ROI. We prewhitened the extracted betas using the model fit residuals. We then calculated the cross-validated squared Mahalanobis distances between each possible finger pair, using our four runs as independent cross-validation folds, and averaged the resulting distances across the folds. If it is impossible to statistically differentiate between conditions (i.e. when this parameter is not represented in the ROI), the expected value of the distance estimate would be 0. If it is possible to distinguish between activity patterns, this value will be larger than 0.
The distance values for all finger pairs were assembled in a representational dissimilarity matrix (RDM), with a width and height corresponding to the five finger movement conditions. Since the RDM is mirrored across the diagonal with meaningless zeros on the diagonal, all statistical analyses were conducted on the 10 unique off-diagonal values of the RDM. We first estimated the strength of the finger representation or ‘finger separability’ by averaging the 10 unique off-diagonal values of the RDM. If there is no information in the ROI that can statistically distinguish between the finger conditions, then due to cross-validation the expected separability would be 0. If there is differentiation between the finger conditions, the separability would be larger than 0 (Nili et al., 2014). Note that this does not directly indicate that this region contains topographic information, but rather that this ROI contains information that can distinguish between the finger conditions. To further ensure that our S1 hand ROI was activated distinctly for different fingers, we created a CSF ROI that would not contain finger-specific information. We repeated our RSA analysis in this ROI and statistically compared the separability of the CSF and S1 hand area ROIs. Second, we tested whether the inter-finger distances were different between controls and patients using a robust mixed ANOVA with a within-participants factor for finger pair (10 levels) and a between-participants factor for group (two levels: controls and SCI patients). Third, we estimated the somatotopic typicality (or normality) of each participant’s RDM by calculating a Spearman correlation with a canonical RDM. We followed previously described procedures for calculating the typicality score (Dempsey-Jones et al., 2019; Ejaz et al., 2015; Kieliba et al., 2021; Wesselink et al., 2019). The canonical RDM was based on 7T finger movement fMRI data in an independently acquired cohort of healthy controls (n = 8). The S1 hand ROI used to calculate this canonical RDM was defined similarly as in the current study (see Wesselink et al., 2021 for details). Note that results were unchanged when calculating typicality scores using a canonical RDM based on the averaged RDM of the age-, sex-, and handedness-matched control group tested in this study (see Appendix 1). The typicality scores were Fisher r-to-z transformed prior to statistical analysis (the rs typicality scores are used solely for visualisation). Controls’ and SCI patients’ typicality scores were compared to each other and to those of a group of individuals with congenital hand loss (n = 13), hereafter one-handers, obtained in another study (data publicly available on https://osf.io/gmvua/; Wesselink et al., 2019). Congenital one-handers are born without a hand and do not have an S1 hand representation contralateral to the missing hand.
Finally, we performed multidimensional scaling (MDS) to visualise the distance structure of the RDM in an intuitive manner. MDS projects the higher-dimensional RDM into a lower-dimensional space, while preserving the inter-finger distance values as well as possible (Borg and Groenen, 2005). MDS was performed for each individual participant and then averaged per group after Procrustes alignment to remove arbitrary rotation induced by MDS.
Structural MRI analysis
Midsagittal tissue bridges analysis
Request a detailed protocolWe used sagittal T2w structural images of the cervical spinal cord at the lesion level to quantify spared tissue bridges. It has been shown that already in the sub-acute stage after SCI oedema and haemorrhage have largely resolved and hyperintense signal changes reliably reflect intramedullary neural damage (Huber et al., 2017; Pfyffer et al., 2019; Vallotton et al., 2019). We closely followed previously described procedures (Huber et al., 2017; Pfyffer et al., 2019; Vallotton et al., 2019). We used Jim 7.0 software (Xinapse Systems, Aldwincle, UK) for manual lesion segmentation at the lesion level, for which high intra- and interobserver reliability has previously been reported (Huber et al., 2017; Pfyffer et al., 2019). The experimenter conducting the manual segmentation was blinded to patient identity. We only included patients’ T2w scans if the lesion (i.e. hyperintense, CSF-filled cystic cavity) was clearly visible on the midsagittal slice. We excluded images of two patients with metal artefacts or insufficient data quality which would not allow a reliable quantification of lesion measures. Tissue bridges were defined as the relatively hypointense intramedullary region between the hyperintense CSF on one side and the cystic cavity on the other side. We assessed the width of ventral and dorsal tissue bridges on the midsagittal slice and summed these to get the total width of tissue bridges.
Cervical cross-sectional spinal cord area analysis
Request a detailed protocolWe used Jim 7.0 software (Xinapse Systems) to extract the cross-sectional spinal cord area at cervical level C2/3 of the spinal cord from the sagittal T2w scans. We used multi-planar reconstruction of sagittal images, resulting in 10 contiguous axial slices at C2/3 with a thickness of 2 mm (Losseff et al., 1996). Using the active-surface model from Horsfield et al., 2010, the cross-sectional spinal cord area was calculated semi-automatically for every slice and averaged over all 10 slices.
Statistical analysis
Request a detailed protocolStatistical analysis was carried out using SPSS (v25). Standard approaches were used for statistical analysis, as mentioned in the Results section. If normality was violated (assessed using the Shapiro–Wilk test), non-parametric statistical testing or robust ANOVAs (in RStudio v1.4; WRS2 package; Mair and Wilcox, 2020) were used. We used a Crawford–Howell t-test to compare single patients to the congenital and control groups (Corballis, 2009). All testing was two-tailed, and corrected p-values were calculated using the Benjamini–Hochberg procedure to control the FDR with q < 0.05. The correlational analysis was considered exploratory and we did not correct for multiple comparisons in this analysis.
Bayesian analysis was carried out using JASP (v0.12.2) for the main comparisons to investigate support for the null hypothesis with the Cauchy prior width set at 0.707 (i.e. JASP’s default). Following the conventional cut-offs, a BF smaller than 1/3 is considered substantial evidence in favour of the null hypothesis. A BF greater than 3 is considered substantial evidence, and a BF greater than 10 is considered strong evidence in favour of the alternative hypothesis. A BF between 1/3 and 3 is considered weak or anecdotal evidence (Dienes, 2014; Kass and Raftery, 1995).
Appendix 1
Percent signal change in finger-specific clusters
To assess whether finger movement activity levels were different between patients and controls, we created finger-specific ROIs and extracted the activity level of the corresponding finger movement for each participant. To create the finger-specific ROIs, we thresholded the probability finger surface maps that were created from the travelling wave data of the control group (see main text) at 25% (i.e. meaning that at least 5 out of 18 control participants needed to significantly activate a vertex for this vertex to be included in the ROI) and binarised. We then used the separately acquired blocked design data to extract the finger movement activity levels underlying these finger-specific ROIs. We first flipped the contrast images resulting from the fixed effects analysis (i.e. that was ran to average across the four blocked design runs) along the x-axis for the left-hand tested participants. Each participant’s contrast maps were then resampled to the FreeSurfer 2D average atlas, and the averaged z-standardised activity level was extracted for each finger movement vs. rest contrast underlying the finger-specific ROIs.
We compared the activity levels for each finger movement in the corresponding finger ROI (i.e. thumb movement activity in the thumb ROI, index finger movement activity in the index finger ROI, etc.) between groups. After correction for multiple comparisons, there was no significant difference between groups for the thumb (U = 93, p=0.37), index (t(30) = –0.003, p=0.99), middle (t(30) = 1.11, p=0.35), ring (t(30) = 2.02, p=0.13), or little finger (t(30) = 2.14, p=0.20).
Geodesic distance analysis
To assess whether the finger representations of the SCI patients were aligned differently and/or shifted compared to the control participants, the cortical distance was calculated between each participant’s peak finger activation and a reference anchor in the S1 foot cortex. We used the blocked design data for this analysis as the travelling wave paradigm does not lend itself to accurately extract peak finger activation locations due to the consecutive movement task.
Each participant’s fixed effect contrast maps that were already resampled to the FreeSurfer 2D average atlas were masked using the S1 hand ROI. The geodesic distance was then calculated between the peak activated vertex in the S1 hand cortex and a reference anchor in the S1 foot cortex (MNI coordinates: −16.93, –32.06, 73.80). Note that geodesic distance measures calculated onto average cortical surfaces are not confounded by inter-subject variability in gyrification, follow the anatomical constraints of the brain (i.e. do not cross white matter boundaries), and allow comparison between different subjects. This is a key advantage of this approach over using Euclidean distance measures (Mancini et al., 2019). A robust mixed ANOVA with a within-participants factor for finger (five levels: thumb, index, middle, ring, and little finger) and a between-participants factor for group (two levels: controls and patients) showed that, as expected, the geodesic distances were significantly different across fingers (i.e. a main effect for fingers; F(1,18.81) = 4.29, p=0.05). While the peak finger vertexes were located more lateral for the patient compared to the control group (i.e. a main effect for group; F(4,13.93) = 85.66, p<0.001), the geodesic distance pattern across the fingers was not significantly different across groups (i.e. no significant finger by group interaction; F(4,13.93) = 0.91, p=0.48). When testing for group differences per finger, the BF only revealed inconclusive evidence (BF >0.34 and < 1.88; note that we could not run a Bayesian ANOVA due to normality violations).
Typicality analysis using a canonical RDM based on the controls’ average RDM
To ensure that our typicality results did not change when using a canonical inter-finger RDM based on the age-, sex-, and handedness-matched subjects tested in this study, we recalculated the typicality scores of all participants using the averaged inter-finger RDM of our control sample as the canonical RDM. We found a strong and highly significant correlation between the typicality scores calculated using the canonical inter-finger RDM from the independent dataset (reported in the main text) and the typicality scores calculated using our controls’ average RDM. This was true for both the SCI patient (rs = 0.92, p<0.001) and control groups (rs = 0.78, p<0.001).
We then repeated all typicality analysis reported in the main text. As expected, we found the same results using the typicality scores calculated using our controls’ average RDM as when using the canonical inter-finger RDM from the independent dataset: There was a significant difference in typicality between tetraplegic patients, healthy controls, and congenital one-handers (H(2) = 27.61, p<0.001). We further found significantly higher typicality in controls compared to congenital one-handers (U = 0, p<0.001; BF10 = 76.11). Importantly, the typicality scores of the SCI patients were significantly higher than the congenital one-handers (U = 2, p<0.001; BF10 = 50.98), but not significantly different from the controls (U = 94, p=0.24; BF10 = 0.55). Number of years since SCI significantly correlated with hand representation typicality (rs = −0.54, p=0.05) and patients with more retained GRASSP motor function of the tested upper limb had more typical hand representations in S1 (rs = 0.58, p=0.03). There was no significant correlation between S1 hand representation typicality and GRASSP sensory function of the tested upper limb, spared midsagittal spinal tissue bridges at the lesion level, or cross-sectional spinal cord area (rs = 0.40, p=0.15, rs = 0.50, p=0.10, and rs = 0.48, p=0.08, respectively). An exploratory stepwise linear regression analysis revealed that years since SCI significantly predicted hand representation typicality in S1 with R2 = 0.33 (F(1,10) = 4.98, p=0.05). Motor function, sensory function, spared midsagittal spinal tissue bridges at the lesion level, and spinal cord area did not significantly add to the prediction (t = 1.31, p=0.22, t = 1.62, p=0.14, t = 1.70, p=0.12, and t = 1.09, p=0.30, respectively).
Data availability
Full details of the experimental protocol are available on clinicaltrials.gov under the number NCT03772548. Data is shared on https://osf.io/e8u95/.
References
-
Predictions not commands: Active inference in the motor systemBrain Structure & Function 218:611–643.https://doi.org/10.1007/s00429-012-0475-5
-
Brain motor system function after chronic, complete spinal cord injuryBrain: A Journal of Neurology 128:2941–2950.https://doi.org/10.1093/brain/awh648
-
Functional outcome following spinal cord injury: significance of motor-evoked potentials and ASIA scoresArchives of Physical Medicine and Rehabilitation 79:81–86.https://doi.org/10.1016/S0003-9993(98)90213-1
-
Neural basis of touch and proprioception in primate cortexComprehensive Physiology 8:1575–1602.https://doi.org/10.1002/cphy.c170033
-
Organized Toe Maps in Extreme Foot UsersCell Reports 28:2748–2756.https://doi.org/10.1016/j.celrep.2019.08.027
-
Using Bayes to get the most out of non-significant resultsFrontiers in Psychology 5:1–17.https://doi.org/10.3389/fpsyg.2014.00781
-
Hand use predicts the structure of representations in sensorimotor cortexNature Neuroscience 18:1034–1040.https://doi.org/10.1038/nn.4038
-
BookFinger Representations in Sensorimotor Cortex Are Not Disrupted in Musicians’ DystoniaSociety for Neuroscience Annual MeetingSan Diego, USA: Society for Neuroscience Annual Meeting.
-
Automated manifold surgery: Constructing geometrically accurate and topologically correct models of the human cerebral cortexIEEE Transactions on Medical Imaging 20:70–80.https://doi.org/10.1109/42.906426
-
Intracortical microstimulation of human somatosensory cortexScience Translational Medicine 8:1–11.https://doi.org/10.1126/scitranslmed.aaf8083
-
Corticomotor representation to a human forearm muscle changes following cervical spinal cord injuryEuropean Journal of Neuroscience 34:1839–1846.https://doi.org/10.1111/j.1460-9568.2011.07895.x
-
Disability, atrophy and cortical reorganization following spinal cord injuryBrain: A Journal of Neurology 134:1610–1622.https://doi.org/10.1093/brain/awr093
-
Are midsagittal tissue bridges predictive of outcome after cervical spinal cord injury?Annals of Neurology 81:740–748.https://doi.org/10.1002/ana.24932
-
Large-scale reorganization in the somatosensory cortex and thalamus after sensory loss in macaque monkeysThe Journal of Neuroscience 28:11042–11060.https://doi.org/10.1523/JNEUROSCI.2334-08.2008
-
Outcome of the upper limb in cervical spinal cord injury: Profiles of recovery and insights for clinical studiesThe Journal of Spinal Cord Medicine 37:503–510.https://doi.org/10.1179/2045772314Y.0000000252
-
Bayes FactorsJournal of the American Statistical Association 90:773–795.https://doi.org/10.1080/01621459.1995.10476572
-
Reaffirming the link between chronic phantom limb pain and maintained missing hand representationCortex; a Journal Devoted to the Study of the Nervous System and Behavior 106:174–184.https://doi.org/10.1016/j.cortex.2018.05.013
-
Neural basis of induced phantom limb pain reliefAnnals of Neurology 85:59–73.https://doi.org/10.1002/ana.25371
-
Investigating the stability of fine-grain digit somatotopy in individual human participantsThe Journal of Neuroscience 36:1113–1127.https://doi.org/10.1523/JNEUROSCI.1742-15.2016
-
Representational similarity analysis - connecting the branches of systems neuroscienceFrontiers in Systems Neuroscience 2:1–28.https://doi.org/10.3389/neuro.06.004.2008
-
Visually-driven maps in area 3bThe Journal of Neuroscience 38:1295–1310.https://doi.org/10.1523/JNEUROSCI.0491-17.2017
-
Corticocuneate projections are altered after spinal cord dorsal column lesions in New World monkeysThe Journal of Comparative Neurology 529:1669–1702.https://doi.org/10.1002/cne.25050
-
Responses of somatosensory area 2 neurons to actively and passively generated limb movementsJournal of Neurophysiology 109:1505–1513.https://doi.org/10.1152/jn.00372.2012
-
Robust statistical methods in R using the WRS2 packageBehavior Research Methods 52:464–488.https://doi.org/10.3758/s13428-019-01246-w
-
Fine-grained nociceptive maps in primary somatosensory cortexThe Journal of Neuroscience 32:17155–17162.https://doi.org/10.1523/JNEUROSCI.3059-12.2012
-
Fine-Grained Mapping of Cortical Somatotopies in Chronic Complex Regional Pain SyndromeThe Journal of Neuroscience 39:9185–9196.https://doi.org/10.1523/JNEUROSCI.2005-18.2019
-
A Toolbox for Representational Similarity AnalysisPLOS Computational Biology 10:e1003553.https://doi.org/10.1371/journal.pcbi.1003553
-
Afferent Input and Sensory Function after Human Spinal Cord InjuryJournal of Neurophysiology 119:134–144.https://doi.org/10.1152/jn.00354.2017
-
Chronic Cervical Spinal Cord Injury: DTI Correlates with Clinical and Electrophysiological MeasuresJournal of Neurotrauma 29:1556–1566.https://doi.org/10.1089/neu.2011.2027
-
Predictive Value of Midsagittal Tissue Bridges on Functional Recovery After Spinal Cord InjuryNeurorehabilitation and Neural Repair 35:33–43.https://doi.org/10.1177/1545968320971787
-
Fast robust automated brain extractionHuman Brain Mapping 17:143–155.https://doi.org/10.1002/hbm.10062
-
Brain changes after spinal cord injury, a quantitative meta-analysis and reviewNeuroscience and Biobehavioral Reviews 90:272–293.https://doi.org/10.1016/j.neubiorev.2018.04.018
-
Tonotopic Organization in Human Auditory Cortex Revealed by Progressions of Frequency SensitivityJournal of Neurophysiology 91:1282–1296.https://doi.org/10.1152/jn.01125.2002
-
Localization of the motor hand area to a knob on the precentral gyrus A new landmarkBrain: A Journal of Neurology 120:141–157.https://doi.org/10.1093/brain/120.1.141
Article and author information
Author details
Funding
Swiss National Science Foundation (320030_175616)
- Nicole Wenderoth
ETH Zurich Postdoctoral Fellowship Program
- Sanne Kikkert
Swiss National Science Foundation (PCEFP3_181362 / 1)
- Patrick Freund
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank our participants for taking part in the study. We thank Lydia Kämpf, Nicolin Gauler, and Silvia Hofer for their assistance with data collection, Daan Wesselink for help with data analysis, and Tamar Makin for comments on the manuscript.
Ethics
Participants' informed consent was obtained according to the Declaration of Helsinki prior to study onset. Ethical approval was granted by the Kantonale Ethikkommission Zürich (KEK-2018-00937).
Copyright
© 2021, Kikkert et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 983
- views
-
- 110
- downloads
-
- 21
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Developmental Biology
- Neuroscience
Two major ligand-receptor signaling axes – endothelin Edn3 and its receptor Ednrb, and glial-derived neurotrophic factor (GDNF) and its receptor Ret – are required for migration of enteric nervous system (ENS) progenitors to the hindgut. Mutations in either component cause colonic aganglionosis, also called Hirschsprung disease. Here, we have used Wnt1Cre and Pax2Cre in mice to show that these driver lines label distinct ENS lineages during progenitor migration and in their terminal hindgut fates. Both Cre lines result in Hirschsprung disease when combined with conditional Ednrb or conditional Ret alleles. In vitro explant assays and analysis of lineage-labeled mutant embryos show that GDNF but not Edn3 is a migration cue for cells of both lineages. Instead, Edn3-Ednrb function is required in both for GDNF responsiveness albeit in different ways: by expanding the Ret+ population in the Pax2Cre lineage, and by supporting Ret function in Wnt1Cre-derived cells. Our results demonstrate that two distinct lineages of progenitors give rise to the ENS, and that these divergently utilize endothelin signaling to support migration to the hindgut.
-
- Neuroscience
The slow-intrinsic-pacemaker dopaminergic (DA) neurons originating in the ventral tegmental area (VTA) are implicated in various mood- and emotion-related disorders, such as anxiety, fear, stress and depression. Abnormal activity of projection-specific VTA DA neurons is the key factor in the development of these disorders. Here, we describe the crucial role of the NALCN and TRPC6, non-selective cation channels in mediating the subthreshold inward depolarizing current and driving the firing of action potentials of VTA DA neurons in physiological conditions. Furthermore, we demonstrate that down-regulation of TRPC6 protein expression in the VTA DA neurons likely contributes to the reduced activity of projection-specific VTA DA neurons in chronic mild unpredictable stress (CMUS) depressive mice. In consistent with these, selective knockdown of TRPC6 channels in the VTA DA neurons conferred mice with depression-like behavior. This current study suggests down-regulation of TRPC6 expression/function is involved in reduced VTA DA neuron firing and chronic stress-induced depression-like behavior of mice.