Perceptually relevant remapping of human somatotopy in 24 hours
Abstract
Experience-dependent reorganisation of functional maps in the cerebral cortex is well described in the primary sensory cortices. However, there is relatively little evidence for such cortical reorganisation over the short-term. Using human somatosensory cortex as a model, we investigated the effects of a 24 hr gluing manipulation in which the right index and right middle fingers (digits 2 and 3) were adjoined with surgical glue. Somatotopic representations, assessed with two 7 tesla fMRI protocols, revealed rapid off-target reorganisation in the non-manipulated fingers following gluing, with the representation of the ring finger (digit 4) shifted towards the little finger (digit 5) and away from the middle finger (digit 3). These shifts were also evident in two behavioural tasks conducted in an independent cohort, showing reduced sensitivity for discriminating the temporal order of stimuli to the ring and little fingers, and increased substitution errors across this pair on a speeded reaction time task.
https://doi.org/10.7554/eLife.17280.001eLife digest
The areas of the brain that receive inputs from our senses have a map-like structure. In an area called the visual cortex this map represents our field of vision; in the auditory cortex, it represents the range of different tones we can hear. The sense of touch is processed in the somatosensory cortex: an area of the brain that is organised around a map of the body, with adjacent regions of the cortex representing adjacent regions of the body. The clear structure of these brain regions makes them ideal for exploring how the organisation of the brain changes over time.
How quickly can changes to the touch inputs that the brain receives cause the map in the somatosensory cortex to reorganise? Can these effects be produced in just 24 hours? And would this remapping affect how we perceive touch? To investigate these questions, Kolasinski et al. glued together the right index and right middle fingers of healthy human volunteers. This separated the middle and ring fingers: a pair that usually move together due to the anatomical structure of the hand.
Functional magnetic resonance imaging of the brain’s activity revealed that within 24 hours of the gluing, the brain’s representation of the ring finger moved away from that of the middle finger, and towards the representation of the little finger. A perceptual judgment task mirrored this finding: after 24 hours of gluing, the participants became better at distinguishing between the middle and ring fingers and worse at distinguishing between the ring and little fingers. This is a powerful demonstration of the human brain’s potential to adapt and reorganise rapidly to changes to sensory inputs.
The sense of touch declines gradually with age and may also be reduced as a result of disease such as stroke. A long-term challenge is to understand how the sensory regions of the brain change during this loss of sensation. Further research could then investigate how to maintain the structure of the cortical map to prolong or restore high quality touch sensation.
https://doi.org/10.7554/eLife.17280.002Introduction
Evidence for experience-dependent plasticity in the adult mammalian brain exists across a variety of sensory modalities (Fu and Zuo, 2011). The ordered somatotopy of primary somatosensory cortex (SI) has long served as a model system for studies of cortical reorganisation, with a wealth of evidence from both the non-human primate and rodent literature linking both extreme and subtle peripheral changes in somatosensory inputs over months or years to alterations in cortical maps (Buonomano and Merzenich, 1998; Feldman and Brecht, 2005).
In the human brain, there has also been evidence of experience-dependent remapping in SI. Considerable emphasis has been placed upon cross-sectional differences in the cortical organisation of SI observed in specialist populations, such as musicians, or patients with sensorimotor dysfunction, such as focal dystonia (Elbert et al., 1995; Bara-Jimenez et al., 1998; Nelson et al., 2009; Kalisch et al., 2009). However, only limited longitudinal evidence exists that directly demonstrates remapping at the level of human SI, either in response to altered hand usage patterns (Stavrinou et al., 2007) or more intensive Hebbian co-activation paradigms delivering specific patterns of tactile stimulation to the fingertips (Pleger et al., 2001, 2003; Hodzic et al., 2004; Vidyasagar et al., 2014). There remains a limited understanding of the speed of SI plasticity and how cortical changes relate to behaviour. Here we address this gap in the literature, investigating the propensity for rapid experience-dependent cortical reorganisation and the behavioural relevance thereof.
Using a well-validated paradigm of single-subject fMRI mapping of human SI at 7 tesla (Sanchez-Panchuelo et al., 2010; Kolasinski et al., 2016), we asked two questions. First, can experience-dependent plastic remapping of SI somatotopy be elicited in the human cortex in just 24 hr? Second, are any observed cortical changes reflected in altered tactile function? A complementary combination of fMRI and behavioural psychophysics were used to investigate the effect of a 24 hr manipulation in which the right index (D2) and right middle digit (D3) were joined together using skin glue. Based on the primate literature (Clark et al., 1988; Allard et al., 1991), we hypothesised that forced co-use of the glued digit pair would result in increased tactile co-activation across the digits, and an increase in their shared cortical representation. An alternative hypothesis posited that compensatory behaviour during the 24 hr gluing manipulation would promote changes in the cortical representations of adjacent, but unaffected digits.
In the fMRI cohort, SI digit somatotopy of the right hand was mapped at 7 tesla (Kolasinski et al., 2016) after two periods of normal hand use (control 1 and control 2), and after the gluing manipulation. In each session, fMRI data were acquired during a block design task and phase-encoding task. We first asked whether the amount of cortical overlap between digit representations in SI changed after the gluing manipulation in comparison with two control sessions using the phase-encoding fMRI data. Measures of inter-digit overlap were calculated for adjacent digit pairs (D2-D3: index-middle, D3-D4: middle-ring, D4-D5: ring-little; Figure 1A). A two-way repeated measures ANOVA indicated a significant interaction between session and digit pair on the amount of cortical overlap (F(4,32) = 13.412, p<0.0005, η2:0.626) (Figure 1B). This was driven by a significant reduction in the cortical overlap between D3-D4 (Simple main effect: F(2,16) = 23.379, p<0.0005; Pairwise Sidak-corrected p<0.05) and a significant increase in the overlap of D4-D5 (Simple main effect: F(2,16) = 13.384, p<0.0005; Pairwise Sidak-corrected p<0.05) in the glued condition compared with both control sessions (Figure 1B). No significant change was found in the overlap between D2 and D3, the glued digits, where changes have been observed in similar but longer term studies in non-human primates (Clark et al., 1988; Allard et al., 1991). No shift in peak-to-peak distance between the digit representations (Supplementary file 1) or the overall surface area of each digit representation was observed (Figure 1—figure supplement 1). No systematic difference in the fit between the phase-encoding models and fMRI signal were observed across sessions (Supplementary file 2).
-
Figure 1—source data 1
- https://doi.org/10.7554/eLife.17280.004
To further explore the observed change in inter-digit overlap with no change in overall cortical surface area of each digit representation, the fMRI representations of digit 4 used to calculate the cortical overlap metrics (Dice coefficients) were visualised for each participant and session (Figure 2). This data revealed that the observed changes presented in Figure 1B were driven by an expansion in the representation of digit 4 adjacent to digit 5 and a corresponding contraction at the boundary with digit 3.
Representational similarity analysis conducted on the block design fMRI data also implicated changes in the representation of digit 4 (Figure 3). These data demonstrated a shift in the representation of D4 away from D3 and towards D5.
-
Figure 3—source data 1
- https://doi.org/10.7554/eLife.17280.010
Given the fMRI data were strongly suggestive of changes in cortical overlap outside of the glued digit pair (D2-D3), we sought to more fully investigate the observed pattern of short-term reorganisation by assessing the behavioural correlates of the gluing manipulation. The same experimental design and gluing manipulation were applied in a separate cohort. Instead of an fMRI scan, nine participants undertook behavioural psychophysics tasks during each session. The first task involved temporal order judgment (TOJ), where pairs of rapid vibrotactile stimuli were applied to adjacent digit pairs of the right hand at varying interstimulus intervals. Participants judged which digit was stimulated first. A psychometric function was fitted to the resulting accuracy data, from which a metric of tactile discrimination ability was calculated for each adjacent digit pair (Just Noticeable Difference: JND). Greater values of JND are associated with poorer discrimination across an adjacent digit pair (Figure 4A). We asked whether the ability to distinguish the order of two stimuli delivered in rapid succession, one each on two adjacent digits, differed in the session directly following the gluing manipulation in comparison with the control sessions.
-
Figure 4—source data 1
- https://doi.org/10.7554/eLife.17280.012
-
Figure 4—source data 2
- https://doi.org/10.7554/eLife.17280.013
We generated behavioural predictions based on our fMRI findings of altered cortical maps after the gluing manipulation (Figure 3). Specifically, we hypothesised that the gluing manipulation would result in improved tactile discrimination (reduced JND) across D3-D4, which demonstrates reduced cortical overlap, and diminished tactile discrimination (increased JND) across D4-D5, which demonstrates increased cortical overlap. Data from eight participants met the goodness of fit threshold for all sessions (R2 > 0.4). A two-way repeated measures ANOVA indicated a significant interaction between session and digit pair on tactile discrimination (JND) (F(4,28) = 14.613, p<0.0005, η2:0.676). This was driven by a significant reduction in JND across D3-D4 (Simple main effect: F(2,14) = 14.631, p<0.0005; Pairwise Sidak-corrected p<0.05), and a significant increase in the JND across D4-D5 (Simple main effect: F(2,14) = 10.578, p=0.002; Pairwise Sidak-corrected p<0.05) in the gluing condition compared with both control sessions (Figure 4B). No significant change in tactile discrimination was found for D2-D3, the glued digits. In other words, after the gluing manipulation, tactile discrimination was improved across D3-D4 and was diminished across D4-D5, with no change in the glued digits D2-D3. This finding was supported by similar results from a second task, assessing motor confusion, which involved rapidly cued button presses using individual digits. There was an increase in the level of confusion between digits 4 and 5 after the gluing manipulation compared with the control conditions (Figure 4—figure supplement 1).
The fMRI and psychophysics results present complementary evidence for functionally relevant reorganisation in human SI following just a 24 hr peripheral change in hand use. Previous longitudinal studies of detailed somatotopic remapping in humans using MEG or standard resolution fMRI have provided evidence for a general change in the distance between digit representations (Stavrinou et al., 2007; Mogilner et al., 1993), or a synchronisation in activity pattern of non-adjacent digits (Vidyasagar et al., 2014) after coupled stimulation across digits. Conversely, an intervention involving hand and arm immobilisation drove diminished tactile acuity and a corresponding reduction in BOLD activity in contralateral SI (Lissek et al., 2009). Our results build on this work, demonstrating that usage-dependent changes in the overlap of SI digit representations are associated with a corresponding change in the ability to differentiate tactile inputs to the digits. Specifically, an increase in the cortical overlap between D4 and D5 was accompanied by a reduction in the tactile discrimination ability across these two digits.
A consistent feature of our findings was that cortical reorganisation and behavioural change was not observed for the glued digit pair (D2-D3). Instead, we saw changes in cortical overlap, and corresponding change in behavioural performance, for the other digit pairs. Specifically, we saw a shift of D4 away from D3 and towards D5. These findings reject our hypothesis of increased tactile synchronisation driving an increase in cortical overlap. Instead, these data support the alternative hypothesis that compensatory behaviour in other digits during the 24 hr manipulation drives off-target effects.
In light of the considerable dexterous abilities specific to the human hand (Young, 2003) and the short duration of the gluing manipulation, it is feasible that the observed remapping reflects compensatory behavioural changes in the pattern of hand use during the manipulation. Rather than learning to coordinate and co-use the glued digits, participants may have adapted synergies with the non-manipulated digits. In this case, the increased physical separation of D3 and D4 during the gluing reduces the usually high degree of anatomical enslavement seen peripherally between these two digits (Kim et al., 2008), freeing D4 to function even more synergistically with D5. This argument is supported by the observed pattern of increased motor confusion between digits 4 and 5 after the gluing manipulation (Figure 4—figure supplement 1), demonstrating potential changes in movement synergies, that may explain the somatosensory changes observed. This finding complements recent observations that patterns of synergistic digit usage are strongly reflected in the generalisation of tactile learning from trained to adjacent untrained digits (Dempsey-Jones et al., 2016).
The two fMRI datasets showed independent evidence of a shift in the digit 4 representation during the gluing manipulation: away from digit 3, and towards digit 5 (Figures 1 and 3), raising questions regarding the underlying mechanism driving this reorganisation. Similar studies undertaken over a period of months in non-human primates have reported the emergence of overlapping receptive fields in response to the fusion of two adjacent digits (Clark et al., 1988; Allard et al., 1991). The spatial resolution of BOLD fMRI is, of course, poorer than invasive mapping techniques. Nonetheless, it is possible that the observed changes in cortical activation could result from shifts in population receptive fields at the boundary of cortical digit representations (Besle et al., 2014), or even changes in multidigit receptive fields (Thakur et al., 2012). As well as reorganisation in the cortex, strong evidence also suggests considerable plastic potential in subcortical structures, which is mirrored in the cortex (Jones, 1996; Rausell et al., 1998; Kambi et al., 2014). The remapping observed in SI here may therefore result from reorganisation in lower level nuclei of the ascending somatosensory pathway.
A shift in the cortical representation of digit 4 (Figure 2), rather than an overall enlargement of the representation, distinguishes the observed pattern of cortical change from previous studies of artificial tactile co-activation, wherein a generalised cortical magnification of the manipulated digit is observed (Pleger et al., 2001, 2003; Hodzic et al., 2004). The cortical shift perhaps reflects how naturalistic behavioural changes in hand use are reflected in the cortex, with remapping based on functional need; in this case, leading to perceptual improvements in certain digit pairs (D3/D4), and worsening in others (D4/D5).
What mechanism drives the observed change in cortical overlap? In vivo studies in rodent barrel cortex demonstrate that a 24 hr period of altered sensory input can induce increases in synaptic density in the corresponding cortical representation (Knott et al., 2002). This could unmask or potentiate pre-existing divergent or silenced connections between adjacent digit representations (DeFelipe et al., 1986; Recanzone et al., 1992; Huntley, 1997). The question remains: how could this cortical shift occur with no associated change in the surface area or location of peak activation of the digit 4 representation? Patterns of local horizontal connectivity in the cortex are established during development and moulded throughout life (Sur and Rubenstein, 2005). Homeostatic mechanisms also exist to maintain balance in cortical activity (Turrigiano and Nelson, 2000), and to preserve the scaling and pattern of cortical topography (Sharma et al., 2000; Vanderhaeghen et al., 2000). It is plausible that the gluing manipulation results in the potentiation of excitatory horizontal connectivity between the representations of digits 4 and 5, and the concurrent weakening of connectivity between the representations of digits 3 and 4. After 24 hr, we observe the resulting pattern of cortical reorganisation (Figures 2 and 3), wherein the relative positions of the representations have changed, but features of the maps firmly established in development, remain in place. A comparable process of reduced inhibition prior to remapping has been reported in studies of sensory deprivation in rodent barrel cortex. Reduced inputs after whisker clipping produce a local pattern of cortical disinhibition and broadening excitation, followed by subsequent cortical contraction of the representation (Albieri et al., 2015). In this case, changes in the relative inputs across different digits may prompt a transient disinhibitory cortical milieu, and the subsequent re-establishment of modified excitatory networks within the cortex, wherein the cortical representations have shifted to reflect new usage patterns. The unmasking of latent excitatory connections prior to any homeostatic increase in lateral inhibition to rebalance cortical excitability could explain the observation of more marked changes in tactile perception and representational surface area after very short and intense tactile training interventions, or rTMS protocols (Pleger et al., 2003; Tegenthoff et al., 2005; Ragert et al., 2008; Dinse and Tegenthoff, 2015).
The consistent cortical changes at the flanks of the digit 4 representation (Figure 2) in the absence of changes in the location of peak activation suggest that the latter may not be the most informative feature of rapid map reorganisation. The existence of such hard-wired elements in cortical somatotopy is supported by evidence of SI constraining the extent of reorganisation within the somatosensory system: acting as a blueprint for lower level nuclei (Zembrzycki et al., 2013), enforcing some limits on plasticity, and potentially explaining the persistence of SI somatotopic features long after marked peripheral changes, such as amputation (Kikkert et al., 2016).
Given the increasingly recognised contribution of motor efference information to activity patterns in somatosensory cortex (Lee et al., 2008, 2013), it is also possible that changes in motor outputs over the course of 24 hr may shape patterns of SI activity, promoting a top-down form of SI experience-dependent plasticity.
The rapidity of cortical changes shown herein supports further investigation of rapid and targeted therapeutic interventions in conditions such as focal dystonia, where maladaptive changes in somatotopy have previously been reported (Bara-Jimenez et al., 1998). The mechanistic link between short-term plasticity and long-term circuit changes remains unclear (Holtmaat and Svoboda, 2009; Johansen-Berg et al., 2012). Therefore, further questions remain as to how short-term peripheral manipulations could be used to induce sustained modification or refinement of fine grain functional cortical organisation, as a means to enhance or rehabilitate tactile function.
Materials and methods
All data were acquired in accordance with local central university research ethics committee approval (University of Oxford MSD-IDREC-C2-2013-05). Eighteen participants were recruited (age range 19–33 years, eight females; Supplementary file 3) and provided written informed consent. All participants were right handed according to the Edinburgh Handedness Inventory (Oldfield, 1971), had no history of neurological or psychiatric illness, and met local MRI safety criteria where appropriate.
fMRI experiment
Experimental design
Request a detailed protocolNine healthy adult volunteers participated in this experiment, in line with the closest previous studies of longitudinal remapping in human SI (Pleger et al., 2003; Young, 2003). Participants were scanned in four sessions. The order of experimental conditions is outlined in Figure 1—figure supplement 2. The first session chronologically (control 0) was treated as a habituation scan, to allow participants to get used to the unique aspects of being scanned in an ultra-high-field magnet, and the fMRI tasks. Included in the analysis were control 1 and control 2: two sessions that followed a period of normal hand use, and a glued session, which followed a 24 hr period of the gluing manipulation described above. The order in which the participants undertook the control and glued session was counterbalanced across the group.
During the gluing manipulation, the medial aspect of the right index digit and lateral aspect of the right middle digit were attached using a surgical grade cyanoacrylate skin glue (Derma+flex, Chemence Medical Products Inc., Alpharetta, GA), leaving the tips and pads of the digits exposed. The mixture octyl- and butyl-cyanoacrylates in this glue were selected for their combination of strength and flexibility. Participants were able to flex and extend the two attached digits and were instructed to go about the day-to-day routine as normal.
MRI data acquisition
Request a detailed protocolFunctional MRI data were acquired using a Siemens 7T Magnetom system with a 32-channel head coil. Blood oxygenation level dependent (BOLD) fMRI was acquired using a T2*-weighted multislice gradient echo planar imaging (EPI), with true axial slices centered on the left anatomical hand knob in the z-axis (TR 1500 ms, TE 25 ms, slice thickness 1.2 mm, in-plane resolution 1.2 × 1.2 mm, 22 axial slices, GRAPPA factor = 2). One single volume contrast matched EPI acquisition of the same resolution was also acquired for registration purposes. fMRI data acquisition, task and data preprocessing have been detailed previously (Kolasinski et al., 2016).
Structural MRI data were acquired using a 3T Siemens Trio system during one of the scan sessions with a multi-echo magnetisation prepared rapid acquisition gradient echo (MEMPRAGE) sequence (van der Kouwe et al., 2008) (TR 2530 ms, TE 1.69, 3.55, 5.41 and 7.27 ms, slice thickness 1.0 mm, in-plane resolution 1.0 × 1.0 mm, GRAPPA factor = 2).
fMRI tasks and stimuli
Request a detailed protocolfMRI data were acquired during an active motor task involving visually cued movements of individual digits, from digit 2 (D2: index) to digit 5 (D5: little) inclusive. The digits involved in the gluing manipulation were released prior to any scan. During the tasks, participants made individual movements of D2, D3, D4 and D5 in the form of button presses on an MRI-compatible button-box (in-house manufactured) placed on the right thigh of the participant. Instructions were delivered via a visual display projected into the scanner bore. Four white circles were presented, representing the four fingers of the right hand. The circles flashed individually at a frequency of 1 Hz, instructing movements of a specific digit at that rate. During each visit (Figure 1—figure supplement 2), a block design task and a phase-encoding task were undertaken during fMRI acquisition.
Phase-encoding fMRI task
Request a detailed protocolThe task followed a phase-encoding design, involving continuous button presses with no rest periods. The task involved 8 s blocks of movement of each digit. A phase-encoding forward task cycled through blocks of D2, D3, D4 and D5, with a total of eight repetitions of the cycle. A phase-encoding backward task cycled through blocks of D5, D4, D3 and D2, with a total of eight repetitions of the cycle. The total duration of task fMRI acquisition was 8 min 50 s.
Block design fMRI task
Request a detailed protocolA block design task was undertaken involving the same cued movements of a given digit in blocks of 12 s, contrasted with 12 s blocks of rest. Blocks of movement for each digit were repeated four times, giving a total of 16 movement blocks and 16 rest blocks. The total duration of the task was 6 min 30 s. Block order was counterbalanced, and different orders used across the different visits, counterbalanced across participants.
MRI data preprocessing
Request a detailed protocolfMRI data were processed using tools from the FMRIB Software Library (Jenkinson et al., 2012). All fMRI data were analysed at the level of individual participants and sessions; no group averaging was undertaken across session or across participants. fMRI data were subject to motion correction with FSL MCFLIRT (Jenkinson et al., 2002), the removal of non-brain tissue with the FSL BET (Smith, 2002), high-pass temporal filtering to remove slow drift in BOLD signal (100 s cut-off) and spatial smoothing (Gaussian kernel full-width half maximum 1.5 mm). fMRI data were subjected to exclusion in cases of visible spin history motion artifact as a result of sharp motion during one or more scan sessions (1 mm of absolute mean displacement in fewer than five volumes), as in previous studies (Kolasinski et al., 2016 ); however, no fMRI data were excluded from this study.
Phase-encoding fMRI analysis: inter-digit overlap
Request a detailed protocolThe phase-encoding fMRI data were used to map representations of individual digits in S1; using these representations, it is possible to derive a measure of overlap between the representations of adjacent digits. The phase-encoding task fMRI analysis has been outlined fully elsewhere (Kolasinski et al., 2016). The travelling wave approach described was applied to data derived from each individual participant from each scan session. In brief, this approach applies a cross-correlation with a number of iteratively time-shifted models to find a time point in the phase-encoding forward (D2-D5) and phase-encoding backward (D5-D2) cycle at which each voxel responded maximally. In this case, the reference model was a gamma-convolved boxcar: 8 s ‘on’ and 24 s ‘off’, repeated eight times. With each iterative shift of the model, a cross-correlation was calculated with the BOLD signal at each voxel, with sufficient shifts to cover the entire digit cycle from D2-D5 or D5-D2.
Each iteration of the model was assigned to a specific digit based on the position of the ‘on’ period in the digit cycle. The r-value maps for the lags of each digit were averaged to yield three-dimensional digit maps (D2, D3, D4 and D5) for the phase-encoding forward and backward tasks in each participant, in each session. These maps were resampled into the participant’s high-saturation EPI space from the corresponding session using FSL FLIRT. The forward and backward r-value maps were then averaged to yield a digit representation for each participant in each session. A corresponding z-statistic map was calculated for each of the four digit representations in each participant, in each session.
Digit maps from all sessions were resampled into the contrast matched EPI space of the control 0 habituation scan using FSL FLIRT (Jenkinson and Smith, 2001) (6 degrees of freedom, normalised correlation cost function) so as not to bias registration towards any individual session in the analysis. The control 0 contrast matched EPI image was registered to the structural MEMPRAGE image using boundary-based registration (Greve and Fischl, 2009) (BBR; Degrees of freedom: 6, FMRIB’s Automated Segmentation Tool (FAST) white matter segmentation, no search), refined using the white matter and pial surfaces using blink comparison in Freeview.
Cortical surface reconstructions were derived from MEMPRAGE images using FreeSurfer (Dale et al., 1999). The z-statistic digit maps from each session were projected from each participant’s raw BOLD fMRI EPI space onto the participant’s cortical surface using Connectome Workbench (Marcus et al., 2011). Once projected to the cortical surface, each digit representation for each session was thresholded using FDR (α = 0.01), and binarised for subsequent overlap analysis.
A measure of cortical overlap between adjacent digits (D2/D3, D3/D4 and D4/D5) was calculated from the binarised digit representations for each participant and session using the Dice coefficient (Dice, 1945). The metric varies from a value of 0, indicating no overlap between the representations, to a value of 1, indicating perfect overlap between the representations. The Dice coefficient was calculated within a FreeSurfer-derived SI region of interest, overlapping with the BOLD fMRI EPI acquisition volume. Where A and B are the area of the two digit representations, the Dice Coefficient is expressed as:
Patterns of peak activation for each digit during each session were calculated using Connectome Workbench, defined as the vertex with the peak z-statistic for each digit within the FreeSurfer-derived SI region of interest. The geodesic distance between adjacent digit peaks was subsequently calculated (Supplementary file 1).
Block design fMRI: representational similarity analysis
Request a detailed protocolMultivariate analysis was conducted using block design data to compare the representational structure of digit maps in the glued versus control sessions. fMRI data were analysed at the single-subject level using a standard GLM approach using FSL FEAT (Jenkinson et al., 2012), with an individual regressor for movement of each digit (D2, D3, D4, D5), using a gamma-HRF boxcar (four 12 s blocks of activity per digit during the task). The resulting parameter estimates and residuals were used for subsequent calculation of representational dissimilarity matrices (RDMs) using the MRC CBSU toolbox (http://www.mrc-cbu.cam.ac.uk/methods-and-resources/toolboxes/) (Nili et al., 2014). Parameter estimates were subjected to univariate noise normalisation by voxelwise covariance estimates. A region of interest was defined for each participant using the phase-encoding data, averaging across sessions, to yield a region of activation (FDR α = 0.01) in the anatomical hand knob. The noise-normalised parameter estimates underlying the ROI were then used to construct representational dissimilarity matrices using Euclidean distance for each session. The distance metrics for adjacent digit pairings were subjected to a two-way repeated measures ANOVA, with one factor: session (Control 1, Control 2, Glued) and a second factor: digit pair (D2/D3, D3/D4, D4/D5). For visualisation purposes, the individual RDMs for each session were subjected to classical multidimensional scaling (MDS). Within each session, the results of the MDS for each participant were aligned using Procrustes analysis. The standard errors were corrected for the reduction in variability induced by the alignment, as described previously (Ejaz et al., 2015).
Behavioural experiment
Experimental design
Request a detailed protocolNine healthy adult volunteers participated in this experiment, none of whom participated in experiment one. Participants attended four behavioural testing sessions in total. These sessions were scheduled in the same way as the scans the fMRI Experiment (Figure 1—figure supplement 2), with an initial habituation session (Control 0). Included in the analysis were control 1 and control 2: two sessions that followed a period of normal hand use, and a glue session, which followed a 24 hr period of the gluing manipulation (Figure 1—figure supplement 2). The gluing manipulation was identical to that applied in the fMRI experiment, and was not in place during any behavioural testing. The order in which the participants undertook the control and glued condition was counterbalanced across the group.
Temporal order judgment task
Request a detailed protocolParticipants performed a temporal order judgment (TOJ) task on vibrotactile stimuli presented on the distal pad of digits 2, 3, 4 and 5 of the right hand. This is a standard task used previously to study hand representation (Shore et al., 2005; Azañón et al., 2015). Participants sat in a testing room in a static chair placed a fixed distance from a table, with their arms resting on the table uncrossed. Their right hand was positioned in the midline, resting on a two-button response device (manufactured in-house). On each of the two buttons of the response device was mounted a vibrotactile stimulator delivering a fixed amplitude suprathreshold pulse (Two Oticon-A, 100 Ω bone conduction vibrators; 20 ms, 200 Hz sinusoidal; amplitude at least 10 times threshold), such that the ventral pads of adjacent digit pairs could be directly opposed to the flat surface of the stimulator. The relative position of the stimulators could be adjusted to accommodate different hand sizes and the relative lengths of different digit pairs. The view of the right hand and forearm was fully occluded and pink noise was used to mask any subtle auditory cues from the vibrotactile stimuli. Time of day was matched across conditions.
Participants performed a two alternative forced choice task following instructions on a laptop screen positioned at eye level at a distance of 56 cm. Participants performed three runs of this task, one for each digit pairing (D2-D3, D3-D4 and D4-D5; order counterbalanced across conditions and participants). Each run of the task was composed of 192 trials. During each trial, participants received one stimulus to each of the two digits and were asked to judge whether the left or right stimulus came first by pressing with the corresponding digit. Participants were cued on the screen with ‘Ready’ to alert of an oncoming stimulus, ‘Response’ to cue a choice and ‘No response’ in cases of failure to press respond within 2 s of the stimuli. Participants were instructed to fixate on a central fixation cross present on the screen at all times; monitoring was in place to ensure participants’ eyes remained open at all times. The inter stimulus intervals (ISIs) were +400,–400,+250,–250,+180,–180,+120,–120,+70,–70,+30 or 30 ms where, by convention, a positive ISI denotes that the right digit was stimulated first. For each digit pair, the 12 ISIs were presented 16 times in a randomly assigned order.
Temporal order judgment data analysis
Request a detailed protocolFor each run in each session, the proportion of right-first responses for each ISI was fitted to a logistic function (Figure 4A) using MATLAB and Statistics Toolbox Release 2014b (The MathWorks, Inc., Natick, MA). Data from one participant were excluded on the basis of a poor fit for digit pairs at one session (R2 < 0.4), threshold selected from earlier study (Azañón et al., 2015). The fit of the remaining data far exceeded this threshold (Range: 0.68–0.99). In TOJ tasks, the just noticeable difference (JND) is defined as half of the difference between two ISIs that yield the ‘right-most digit first’ judgment in 75% and 25% of the trials. Smaller values of JND are therefore associated with better performance in terms of temporal tactile resolution.
Motor confusion task
Request a detailed protocolA second behavioural test involved a motor confusion task. This task involved rapid, visually cued movements of individual digits, in the form of presses on a four-button button box (in-house design). The task was delivered in the same testing environment as the TOJ task. Movements were cued every 700 ms, in blocks of 210 button presses, separated by a 20-s rest period, with a total of four blocks acquired in each session. The order of button presses was pseudorandomised, and the same button press was not repeated twice consecutively; the overall sequence in each block included an equal number of all possible adjacent digit pairings, including both neighbouring digits (e.g. D2 D3) and non-neighbouring digits (e.g. D3 D5). Participants were instructed to respond as quickly and accurately as possible to each cue and to not correct any perceived errors through double presses.
Data were analysed including only responses within the 700-ms response window. Responses were categorised as correct or incorrect; for the latter, the digit used to make the response was recorded.
Statistics
Statistical analyses and graphing used JMP (Version 12.0, SAS Institute, Cary, NC) and Statistics Package for the Social Sciences (SPSS, Version 22.0, IBM Corporation, Armork, NY). The within-subjects data applied for the fMRI and psychophysics data were interrogated with independent two-way repeated measures ANOVAs, selected on the basis of data normality (Shapiro-Wilk confirmed). Significant interactions were further interrogated with analysis of simple main effects; full statistical outputs are provided in Supplementary file 4. All repeated measures ANOVAs were also run using a single value for the control condition by averaging the values from control A and control B; the same pattern of interactions and simple main effects were observed. Experimenters were not blinded to the condition (glued or control) during data acquisition for practical reasons.
References
-
Rapid bidirectional reorganization of cortical microcircuitsCerebral Cortex 25:3025–3035.https://doi.org/10.1093/cercor/bhu098
-
Reorganization of somatosensory area 3b representations in adult owl monkeys after digital syndactylyJournal of Neurophysiology 66:1048–1058.
-
Dynamic tuning of tactile localization to body postureCurrent Biology 25:512–517.https://doi.org/10.1016/j.cub.2014.12.038
-
Abnormal somatosensory homunculus in dystonia of the handAnnals of Neurology 44:828–831.https://doi.org/10.1002/ana.410440520
-
Cortical plasticity: from synapses to mapsAnnual Review of Neuroscience 21:149–186.https://doi.org/10.1146/annurev.neuro.21.1.149
-
Long-range focal collateralization of axons arising from corticocortical cells in monkey sensory-motor cortexJournal of Neuroscience 6:3749–3766.
-
Transfer of tactile perceptual learning to untrained neighboring fingers reflects natural use relationshipsJournal of Neurophysiology 115:1088–1097.https://doi.org/10.1152/jn.00181.2015
-
Hand use predicts the structure of representations in sensorimotor cortexNature Neuroscience 18:1034–1040.https://doi.org/10.1038/nn.4038
-
Experience-dependent structural plasticity in the cortexTrends in Neurosciences 34:177–187.https://doi.org/10.1016/j.tins.2011.02.001
-
Experience-dependent structural synaptic plasticity in the mammalian brainNature Reviews Neuroscience 10:647–658.https://doi.org/10.1038/nrn2699
-
A global optimisation method for robust affine registration of brain imagesMedical Image Analysis 5:143–156.https://doi.org/10.1016/S1361-8415(01)00036-6
-
BookDivergence of thalamocortical projections and limits on somatosensory cortical plasticityIn: Franzen O, Johannsson R, Terenius L, editors. Somesthesis and the Neurobiology of the Somatosensory Cortex. Birkhäuser Basel: Basel. pp. 1–15.https://doi.org/10.1007/978-3-0348-9016-8_1
-
Finger inter-dependence : linking the kinetic and kinematic variablesHuman Movement Science 27:408–422.https://doi.org/10.1016/j.humov.2007.08.005
-
Investigating the stability of Fine-Grain digit somatotopy in individual human participantsJournal of Neuroscience 36:1113–1127.https://doi.org/10.1523/JNEUROSCI.1742-15.2016
-
Motor modulation of afferent somatosensory circuitsNature Neuroscience 11:1430–1438.https://doi.org/10.1038/nn.2227
-
A disinhibitory circuit mediates motor integration in the somatosensory cortexNature Neuroscience 16:1662–1670.https://doi.org/10.1038/nn.3544
-
Informatics and data mining tools and strategies for the human connectome projectFrontiers in Neuroinformatics 5:4.https://doi.org/10.3389/fninf.2011.00004
-
Digit-specific aberrations in the primary somatosensory cortex in writer's crampAnnals of Neurology 66:146–154.https://doi.org/10.1002/ana.21626
-
A toolbox for representational similarity analysisPLoS Computational Biology 10:e1003553.https://doi.org/10.1371/journal.pcbi.1003553
-
Extensive divergence and convergence in the thalamocortical projection to monkey somatosensory cortexJournal of Neuroscience 18:4216–4232.
-
Mapping human somatosensory cortex in individual subjects with 7t functional MRIJournal of Neurophysiology 103:2544–2556.https://doi.org/10.1152/jn.01017.2009
-
Spatial modulation of tactile temporal-order judgmentsPerception 34:1251–1262.https://doi.org/10.1068/p3313
-
Fast robust automated brain extractionHuman Brain Mapping 17:143–155.https://doi.org/10.1002/hbm.10062
-
Second-order receptive fields reveal multidigit interactions in area 3b of the macaque monkeyJournal of Neurophysiology 108:243–262.https://doi.org/10.1152/jn.01022.2010
-
Hebb and homeostasis in neuronal plasticityCurrent Opinion in Neurobiology 10:358–364.https://doi.org/10.1016/S0959-4388(00)00091-X
-
A mapping label required for normal scale of body representation in the cortexNature Neuroscience 3:358–365.https://doi.org/10.1038/73929
-
Evolution of the human hand: the role of throwing and clubbingJournal of Anatomy 202:165–174.https://doi.org/10.1046/j.1469-7580.2003.00144.x
-
Sensory cortex limits cortical maps and drives top-down plasticity in thalamocortical circuitsNature Neuroscience 16:1060–1067.https://doi.org/10.1038/nn.3454
Article and author information
Author details
Funding
University College, Oxford
- James Kolasinski
Wellcome (104128/Z/14/Z)
- Tamar R Makin
Wellcome (102584/Z/13/Z)
- Charlotte J Stagg
Medical Research Council (MR/L009013/1)
- Saad Jbabdi
Wellcome (110027/Z/15/Z)
- Heidi Johansen-Berg
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
JK holds a Stevenson Junior Research Fellowship at University College, Oxford. TRM and CJS both hold Wellcome Trust/Royal Society Henry Dale Fellowships (TRM 104128/Z/14/Z; CJS: 102584/Z/13/Z). SJ holds an MRC Career Development Fellowship (MR/L009013/1). HJB holds a Wellcome Trust Principal Research Fellowship (110027/Z/15/Z). The work was additionally supported by the NIHR Oxford Biomedical Research Centre. Support for the 7T scanner was provided by the Medical Research Council. The authors thank Ellen Thomas, Maria Blöchl and Charles Spence for their assistance with piloting the behavioural paradigm, Harriet Dempsey-Jones and Lewis Gaul for their expertise in fitting psychophysics data, and Sanne Kikkert for her expertise in fMRI analysis.
Ethics
Human subjects: All data were acquired in accordance with local central university research ethics committee approval (University of Oxford MSD-IDREC-C2-2013-05). Eighteen participants were recruited, each providing written informed consent to take part in this study, and for the results of this study to be published.
Copyright
© 2016, Kolasinski 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
-
- 2,707
- views
-
- 499
- downloads
-
- 37
- 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
-
- Neuroscience
Motivation depends on dopamine, but might be modulated by acetylcholine which influences dopamine release in the striatum, and amplifies motivation in animal studies. A corresponding effect in humans would be important clinically, since anticholinergic drugs are frequently used in Parkinson’s disease, a condition that can also disrupt motivation. Reward and dopamine make us more ready to respond, as indexed by reaction times (RT), and move faster, sometimes termed vigour. These effects may be controlled by preparatory processes that can be tracked using electroencephalography (EEG). We measured vigour in a placebo-controlled, double-blinded study of trihexyphenidyl (THP), a muscarinic antagonist, with an incentivised eye movement task and EEG. Participants responded faster and with greater vigour when incentives were high, but THP blunted these motivational effects, suggesting that muscarinic receptors facilitate invigoration by reward. Preparatory EEG build-up (contingent negative variation [CNV]) was strengthened by high incentives and by muscarinic blockade, although THP reduced the incentive effect. The amplitude of preparatory activity predicted both vigour and RT, although over distinct scalp regions; frontal activity predicted vigour, whereas a larger, earlier, central component predicted RT. The incentivisation of RT was partly mediated by the CNV, though vigour was not. Moreover, the CNV mediated the drug’s effect on dampening incentives, suggesting that muscarinic receptors underlie the motivational influence on this preparatory activity. Taken together, these findings show that a muscarinic blocker impairs motivated action in healthy people, and that medial frontal preparatory neural activity mediates this for RT.
-
- Neuroscience
Complex macro-scale patterns of brain activity that emerge during periods of wakeful rest provide insight into the organisation of neural function, how these differentiate individuals based on their traits, and the neural basis of different types of self-generated thoughts. Although brain activity during wakeful rest is valuable for understanding important features of human cognition, its unconstrained nature makes it difficult to disentangle neural features related to personality traits from those related to the thoughts occurring at rest. Our study builds on recent perspectives from work on ongoing conscious thought that highlight the interactions between three brain networks – ventral and dorsal attention networks, as well as the default mode network. We combined measures of personality with state-of-the-art indices of ongoing thoughts at rest and brain imaging analysis and explored whether this ‘tri-partite’ view can provide a framework within which to understand the contribution of states and traits to observed patterns of neural activity at rest. To capture macro-scale relationships between different brain systems, we calculated cortical gradients to describe brain organisation in a low-dimensional space. Our analysis established that for more introverted individuals, regions of the ventral attention network were functionally more aligned to regions of the somatomotor system and the default mode network. At the same time, a pattern of detailed self-generated thought was associated with a decoupling of regions of dorsal attention from regions in the default mode network. Our study, therefore, establishes that interactions between attention systems and the default mode network are important influences on ongoing thought at rest and highlights the value of integrating contemporary perspectives on conscious experience when understanding patterns of brain activity at rest.