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.

DOI: http://dx.doi.org/10.7554/eLife.17280.001

eLife 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.

DOI: http://dx.doi.org/10.7554/eLife.17280.002

Main text

Introduction

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).

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.

Figure 2.
Download figureOpen in new tabFigure 2. Pattern of shift in the cortical representation of digit 4 from inter-digit overlap analysis.

Data displayed for all nine participants (AI) showing the outline of the digit 4 representation mapped during the three sessions (Control 1, Control 2, Glued) overlaid on the individual participants’ cortical surface reconstruction (sulci: dark grey, gyri: light grey; zoomed panel showing anatomical hand knob). The location of the peak vertex is shown with a coloured circle. These results demonstrated a consistent shift in the representation of digit 4 at the level of individual participants, such that while there is a minimal change in peak activation and area of the representation, the flank adjacent to the representation of digit 5 expands, and the flank adjacent to digit 3 contracts, consistent with the observed changes in cortical overlap (Figure 1B).

DOI: http://dx.doi.org/10.7554/eLife.17280.008

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.
Download figureOpen in new tabFigure 3. Representational similarity analysis of block design data yields complementary evidence of shift in the S1 representation of digit 4 away from digit 3 and towards digit 5.

(A) Noise normalised parameter estimates from a standard GLM for each digit were used to construct representational dissimilarity matrices (RDMs) using Euclidean distance within a hand knob ROI derived individually for each subject from phase-encoding data from all sessions. Average raw distance values are shown for each session (B) Multidimensional scaling and Procrustes analysis of individual participants’ distance matrices at each time point demonstrate schematically the observed shift in the representation of digit 4, away from digit 3, and towards digit 5, consistent with the observed pattern of cortical overlap and tactile discrimination changes. Two-way repeated measures ANOVA indicated a significant interaction between session and digit pair on the amount of cortical overlap (F(2.0,16.2) = 4.430, p=0.029, η2:0. 356), driven by a shift in the representation of digit 4 away from digit 3 (Simple Main Effect: F(2,16) = 16.076; pairwise comparisons glued vs. control 1 and control 2: p<0.01). Multidimensional scaling yields the spatial relationship of representations in arbitrary units (a.u.).

DOI: http://dx.doi.org/10.7554/eLife.17280.009

Figure 3—source data 1.Data presented Figure 3A.

Representational similarity analysis distance metrics for adjacent digit pairs (D2–D3, D3–D4 and D4–D5) across three conditions (Control 1, Control 2 and Glued).

DOI: http://dx.doi.org/10.7554/eLife.17280.010

Download source data [figure-3—source-data-1.media-2.csv]

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.

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

Nine 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

Functional 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

fMRI 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

The 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

A 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

fMRI 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

The 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:2|AB||A|+|B|

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

Multivariate 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

Nine 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

Participants 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

For 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

A 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

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.

Decision letter

Richard Ivry, Reviewing editor, University of California, Berkeley, United States

In the interests of transparency, eLife includes the editorial decision letter and accompanying author responses. A lightly edited version of the letter sent to the authors after peer review is shown, indicating the most substantive concerns; minor comments are not usually included.

Thank you for submitting your article "Rapid remapping of human somatosensory cortex" for consideration by eLife. Your article has been reviewed by three peer reviewers, and the evaluation has been overseen by Rich Ivry, the Reviewing Editor, and Sabine Kastner as the Senior Editor. The following individual involved in review of your submission has agreed to reveal their identity: Hubert Dinse (Reviewer #1).

The reviewers have discussed the reviews with one another and the Reviewing Editor has drafted this decision to help you prepare a revised submission.

Summary:

Three reviewers have evaluated the paper. We all find the paper intriguing with the potential to make a nice addition to the literature on short-term plasticity in the human somatosensory system. The converging evidence from fMRI and psychophysics has the potential to make a compelling story. However, there are a number of substantive issues that must be addressed in a revision to make the case more convincing.

Essential Revisions:

1) The main issues concern the analyses.

1A) We take it that your key findings, namely that the main reorganization effects are observed in the untreated fingers, were unexpected – that you had anticipated observing changes in the representation of the treated fingers (perhaps in addition to the untreated fingers). Granted, one has to follow the data. However, the situation is puzzling. You report no change in the peak locations, representational areas of all fingers, and overlap of D2-D3 (treated fingers). You do find changes in overlap of D3-D4 (smaller overlap) and D4-D5 (greater overlap). If we put this all together, the results indicate that there is a change in the flank of D4 such that the region near D3 shifts to the other side (near D5). That is, with no change in peak of D4 or its overall area, there should be some sort of marked asymmetry in representation of D4, with much larger flank on side towards D5 (with corresponding reduction in flank on D3 side to keep area constant). Is this correct? Or does the significant overlap shift emerge as a consequence of non-significant changes in finger peaks, areas, etc.? We would like to see a graphical depiction to better understand this puzzling pattern, perhaps some pre-post variant on the graph you show in Figure 1A. Ideally, this would be evident in maps of individuals (as you did in your previous J Neurosci paper) as well as in some sort of group analysis. Be sure to clarify how you measured the peaks.

1B) We would like a sense of the variability of individual somatotopic maps. Two measures to consider here. First, a measure of goodness-of-fit between the phase-encoding model and brain activity (either the average R2 or r-value for each digit). Second, a quantitative evaluation of inter-subject stability, perhaps performed by calculating Cronback's α (Cronbach, 1951; Grinband et al., Neuroimage 2008) to measure the correlation between the spatial patterns of digit representations of each subject.

1C) Other authors have shown that the representational structure of the fingers may depict the invariant organization of the somatosensory cortex better than the exact spatial distribution of the finger clusters (Ejaz et al., 2015). Such measures are not prone to some of the problems of the measurements of cortical overlap used here (e.g., dependence on map thresholds. Representational spaces (Rss) can be computed and compared between sessions, allowing an assessment of the distance between the brain activity patterns related to each digit before and after gluing. This would allow you to confirm that the representational structure of the digits is modified in the post-gluing session only, and not in the control conditions.

1D) Please clarify how binary images were obtained to calculate the Dice coefficient. Did the binary maps obtained through different thresholds show the same patterns of overlap across conditions (i.e., those reported in Figure 1B)?

1E) Methods suggest that the individual maps of correlation coefficients were averaged and the z-score was calculated on the group level map. Would it be better to transform the single-subject maps to z-scores before averaging them?

1F) Not clear if you have averaged across positive and negative BOLD signals. Would analyzing negative and positive BOLD separately reveal a different picture of the changes?

2) While the psychophysics is impressive, the temporal order task seems to provide an indirect way to assess change in overlap representation. A spatial localization task would provide a more direct test (less confusing between D3 and D4, more confusions between D4 and D5), as well as add to null results of no change in confusions for D2 and D3. This experiment would further bolster your conclusions, especially important here given that the imaging results are unexpected.

3) The methods section is difficult to sort through. A few of us came away from text assuming you had done before/after contrasts, with one contrast Baseline vs. Control and the other Baseline vs. Gluing. Your design figure suggests that you may have planned things this way for decided there were concerns with the first session (familiarity with scanner environment?) and thus developed a different set of pairings. As such, one contrast is separated by 24 hours and the other by 4 weeks. We recognize that, by counterbalancing things, this is not a confound. Indeed, one could say it is even more impressive that there is little change between Control 1 and 2 given that the time between scans is much greater for this contrast (for half the participants) compared to between (one of the) control scans and Gluing scan. However, we do think it important to report if there is an effect of delay. That is, do you get more stability when the control scans are separated by 24 hours vs. 4 weeks (and similarly, more difference with Gluing when there is 4 weeks vs. 24 hours). The n will be quite small (4 vs. 5) in an analysis of order but we think the order issue should be addressed. Given that the small n may lack power in a statistical test, we would also like to see a supplemental figure with the fMRI data for the two orders.

4) While appreciating the novelty of your design, we do believe that there is a substantial body of literature, both human and non-human, on short-term plasticity. You cite the Pleger et al. '03 paper which is most similar to the current work. There are other uncited papers that seem relevant to the discussion, some using fMRI and some using other methods. Please consider relevance of Pleger et al. 2001 (PNAS), Dinse et al. 2003 (Science), Hodzik et al. 2004 (J Neurosci), Tegenthoff et al. 2005 (PLoS Biol), Lissek et al. 2009, Kalisch et al. 2009, and Mogilner et al. 1993 (PNAS).

5) Your method is similar to that used by Merzenich and colleagues in the late 80's. They, of course, looked at long-term changes and observed dramatically different results. A (slightly) expanded discussion of these differences and possible mechanisms would better situate your paper in the somatosensory reorganization literature.

DOI: http://dx.doi.org/10.7554/eLife.17280.019

Author response

Comments

If your username is different from your full name, we require you to identify yourself within the comment itself. Comments are checked by a moderator (and/or an eLife editor) before they appear. Comments should be constructive, relevant to the article, conform to our terms and conditions, and include any pertinent competing interests.