Experimental evidence in animal models indicates that the brainstem plays a major role in sensory modulation. However, mapping functional activity within the human brainstem presents many methodological challenges. These constraints have deterred essential research into human sensory brainstem processing. Here, using a 3T functional Magnetic Resonance Imaging (fMRI) sequence optimised for the brainstem, combined with uni- and multivariate analysis approaches, we investigated the extent to which functional activity of neighbouring somatosensory nuclei can be delineated in the brainstem, thalamus, and primary somatosensory cortex (S1). Whilst traditional univariate approaches offered limited differentiation between adjacent hand and face activation in the brainstem, multivariate classification enabled above-chance decoding of these activity patterns across S1, the thalamus, and the brainstem. Our findings establish a robust methodological approach to explore signal processing within the brainstem and across the entire somatosensory stream. This is a fundamental step towards broadening our understanding of somatosensory processing within humans and determining what changes in sensory integration may occur in clinical populations following sensory deprivation.

Accurately localising touch and proprioception is vital for daily functioning. If our brain were to mislocate a hand touch as being on the face, it would cause disoriented movements, making simple tasks like eating or dressing nearly impossible. However, despite its crucial role, our understanding of somatosensory processing in humans remains limited. Importantly, the relationship between somatosensory organisation in the brainstem and the plasticity of cortical representations remains unexplored, as it is extremely challenging to study human brainstem nuclei in vivo. Here, as a fundamental step towards imaging functionally distinct somatosensory processing pathways, we investigate the extent to which we can delineate the activity of neighbouring somatosensory nuclei across the brainstem, thalamus, and primary somatosensory cortex (S1) using 3T functional Magnetic Resonance Imaging (fMRI).

The processing of afferent inputs is organised in a somatotopic fashion along the somatosensory processing stream, which facilitates the accurate interpretation of afferent signals and coordination of movements. This map-like organisation showcases a point-to-point correspondence between spatial activity patterns in the central nervous system (CNS) and body parts. For example, tactile and proprioceptive signals from the hands enter the CNS via the dorsal horn and then travel through the spinal cord before reaching the first CNS relay station in the ipsilateral cuneate nucleus of the medulla (Fig. 1). From the cuneate nucleus, the signals are projected to the contralateral ventroposterior lateral (VPL) nucleus of the thalamus and finally to the hand area of the primary somatosensory cortex and higher-order cortical areas (Purves, 2019). Conversely, afferent signals from the face synapse adjacently in the ipsilateral trigeminal nuclei before reaching the contralateral ventroposterior medial (VPM) nucleus of the thalamus and, finally, the face area of the primary somatosensory cortex and higher-order cortical regions (Craven, 2011; Purves, 2019).

Fig. 1.

Processing of tactile and proprioceptive inputs across the central nervous system. Tactile and proprioceptive sensations from the hand (purple) ascend via the ipsilateral dorsal column of the spinal cord, reaching the cuneate (Cu) nucleus in the brainstem (1; purple). From there, signals are relayed to the contralateral ventroposterior lateral (VPL) nucleus of the thalamus (2; purple) and then to the primary somatosensory cortex (3; purple). Conversely, inputs from the face (yellow) are relayed to the ipsilateral spinal trigeminal (SpV) nucleus in the brainstem (1; yellow) before projecting to the contralateral ventroposterior medial (VPM) nucleus of the thalamus (2; yellow) and then to the primary somatosensory cortex (3; yellow). Figure created with BioRender.com.

Fig. 1.

Processing of tactile and proprioceptive inputs across the central nervous system. Tactile and proprioceptive sensations from the hand (purple) ascend via the ipsilateral dorsal column of the spinal cord, reaching the cuneate (Cu) nucleus in the brainstem (1; purple). From there, signals are relayed to the contralateral ventroposterior lateral (VPL) nucleus of the thalamus (2; purple) and then to the primary somatosensory cortex (3; purple). Conversely, inputs from the face (yellow) are relayed to the ipsilateral spinal trigeminal (SpV) nucleus in the brainstem (1; yellow) before projecting to the contralateral ventroposterior medial (VPM) nucleus of the thalamus (2; yellow) and then to the primary somatosensory cortex (3; yellow). Figure created with BioRender.com.

Close modal

The brainstem somatosensory nuclei have been termed the somatosensory gateway to the brain (Versteeg, Chowdhury et al., 2021). Studies in animal models have indicated that they play a critical role in sensory modulation by attenuating, enhancing, or rerouting incoming signals (Jain et al., 2000; Keine et al., 2016; Kubota et al., 2024; Versteeg, Chowdhury et al., 2021; Versteeg, Rosenow et al., 2021). For instance, a recent study in non-human primates demonstrated that somatosensory signals are attenuated in the cuneate nucleus during voluntary movement as compared to passive movement (Kubota et al., 2024). Additionally, non-human primate models of spinal cord injury have identified these brainstem nuclei as key sites for driving upstream plasticity, where reorganisation between inputs from the adjacent cuneate and trigeminal nuclei leads to the cortical expansion of the face representations in S1 (Halder et al., 2018; Jain et al., 2000; Kaas et al., 1999; Kambi et al., 2014). This plasticity has been linked to both recovery of function and the development of maladaptive circuitry that relates to neuropathic pain (Jutzeler et al., 2015; Zaaimi et al., 2012; Zörner et al., 2014). Therefore, investigating the brainstem is not only crucial for a complete systems-level understanding of somatosensory processing but may also provide insight into the driving forces of adaptive plasticity in clinical cohorts with impaired somatosensation.

fMRI offers a non-invasive solution to probe brainstem function in vivo. However, while it is regularly used to assess somatotopic organisation and processing in cortical areas, specifically S1 (Besle et al., 2013; Ejaz et al., 2015; Gozzi et al., 2024; Kikkert et al., 2016, 2021, 2023; Kolasinski et al., 2016; Kuehn et al., 2014, 2018; Martuzzi et al., 2014; Odermatt et al., 2024; Puckett et al., 2017; Rabe et al., 2023; Sanchez-Panchuelo et al., 2010, 2012; Sanders et al., 2023; Wesselink et al., 2019), several challenges arise when using fMRI to assess somatosensory processing at the level of the brainstem (Beissner, 2015; Brooks et al., 2013; Sclocco et al., 2018). First, the brainstem is highly vulnerable to motion and physiological noise. Specifically, respiration and the pulsation of nearby arteries and cerebrospinal fluid (CSF) can introduce fMRI signal inhomogeneities, which severely impact signal quality (Brooks, Faull, Pattinson, & Jenkinson, 2013). Second, the somatosensory nuclei are small, measuring only a few millimetres in width and spanning just a few voxels in typical fMRI images. Moreover, the somatosensory nuclei,such as the cuneate and trigeminal, lay directly adjacent to each other, making it challenging to differentiate their activity given the spatial resolution limits of fMRI, which is typically 2–4 mm (Glover, 2011; Menon & Goodyear, 2001). Finally, the brainstem serves as a primary hub for many diverse physiological processes, with nuclei responsible for autonomic control, pain modulation, and emotion processing closely surrounding those involved in sensorimotor signal transmission. This can increase uncertainty when attributing signals to individual nuclei (Harrison et al., 2023; Paxinos & Mai, 2004). Despite these limitations, recent advancements in fMRI acquisition and pre-processing techniques have significantly mitigated these challenges, enabling researchers to characterise brainstem activity in humans (Beissner et al., 2014; Harrison et al., 2023; Harvey et al., 2008; Matt et al., 2019; Reddy et al., 2024). However, very few fMRI investigations have explored human somatosensory brainstem processing.

To date, human imaging studies involving the cuneate and trigeminal nuclei have primarily demonstrated activation in these regions as secondary outcomes whilst investigating alternate physiological processes (Faull et al., 2015; Pattinson et al., 2009) or while optimising brainstem pre-processing methods (Beissner et al., 2011; Harvey et al., 2008). Crucially, researchers have yet to establish the extent to which fMRI can be used to delineate the functional activity of these neighbouring nuclei. Notably, previous studies have demonstrated the capability of fMRI to localise univariate activity in adjacent sensory and motor brainstem nuclei (Komisaruk et al., 2002; Reddy et al., 2024). However, these studies did not directly delineate the extent to which these overlapping activations could be differentiated. This distinction between localisation and discrimination is crucial for determining the somatotopic organisation of somatosensory representations and assessing if activity is reorganised in response to environmental stimuli or injury.

In this study, we combined fMRI with a hand-and-face movement paradigm to map activation across the somatosensory nuclei of the brainstem, thalamus, and S1. In comparison to purely sensory paradigms, motor tasks typically yield stronger and more robust univariate somatosensory activity across the somatosensory stream (Charyasz et al., 2023; Komisaruk et al., 2002; Kündig et al., 2024; Sanders et al., 2023). We used both traditional univariate (i.e., activity magnitude measured via a general linear model) and state-of-the-art multivariate analysis approaches (i.e., generalisation of voxel representations assessed via a linear classifier) to evaluate the extent to which we can delineate hand and face activity across these somatosensory nuclei. If 3T fMRI can resolve activity in adjacent brainstem nuclei, then we should observe a significant difference in the magnitude of hand- and face-related activation across anatomical masks of the cuneate and spinal trigeminal nuclei. However, given the limitations of univariate approaches in detecting subthreshold or overlapping representations, we anticipate that multivariate classification will provide greater sensitivity in distinguishing hand and face activity, particularly in the brainstem (Ejaz et al., 2015; Haufe et al., 2014; Hebart & Baker, 2018; Muret et al., 2022; Wesselink et al., 2022). Therefore, we expect that multivariate classification accuracy will be significantly above chance within the somatosensory nuclei of interest, reflecting the presence of distinct but distributed hand and face activity patterns. By integrating both univariate and multivariate analyses, we aim to establish a comprehensive methodological approach to explore signal processing in humans across the somatosensory stream.

2.1 Participants

We recruited 21 healthy adults between 18 and 75 years old. Given the novelty of this approach, there is no well-established method for calculating sample size. We conducted an a priori power analysis based on our previous study demonstrating significant S1 activity during hand movement (Kikkert et al., 2021), aiming for a power level of 80% for a two-tailed one-sample t-test (alpha levl = 0.05), which resulted in a necessary sample size of 12 participants (calculated with G*Power; Faul et al., 2007). However, group sizes larger than 20 are recommended for fMRI experiments and parametric statistical analysis (Pajula & Tohka, 2016), which informed our decision to recruit a larger sample size than required by the power analysis. We further conducted a post-hoc sensitivity analysis to assess the minimum detectable effect size as detailed in the Supplementary Materials.

Exclusion criteria were as follows: no contraindications for MRI, absence of neurological disorders, no hand impairments, and the capacity to provide informed consent. Notably, one participant was rescanned due to abnormal fMRI noise components identified through independent component analysis (ICA). Another participant was excluded due to corruption of their fMRI data during data acquisition, resulting in a final sample of 20 participants (mean age ± s.e.m. = 50.8 ± 3.5 years; 2 females; 3 dominant left-handers). The data used in this study were part of a larger investigation into somatosensory processing in spinal cord injury (SCI) patients. As such, the healthy control participants reported in the current manuscript were age-, sex-, and handedness-matched to the SCI cohort. Handedness was determined using the Edinburgh Handedness Inventory, in which participants self-reported their preferred hand for a range of activities of daily living (Oldfield, 1971). Note that the exclusion of left-handed participants did not alter the main results of the current manuscript (see post hoc sensitivity and generalisation analyses in the Supplementary Materials for details).

The study received ethical approval from the Ethics Committee of the Canton of Zürich (KEK-2018-00937), and written informed consent was obtained from all participants prior to inclusion. The study is registered on clinicaltrials.gov under NCT03772548.

2.2 Experimental design

We designed a blocked fMRI paradigm to investigate the activation of adjacent medullary nuclei involved in processing sensory input from the face and hands. The blocked design included three movement conditions: self-paced left-hand open-and-closing, right-hand open-and-closing, and face movement (i.e., continuous jaw clenching and self-paced lip pursing, as previously implemented by Beissner et al., 2011), as well as a rest condition. While tactile paradigms are an intuitive solution to localise sensory nuclei (Komisaruk et al., 2002; Reddy et al., 2024), motor tasks often evoke stronger and more robust activation across the somatosensory stream (Charyasz et al., 2023; Komisaruk et al., 2002; Sanders et al., 2023). This is particularly advantageous for imaging subcortical regions like the brainstem, which exhibit a lower signal-to-noise ratio (Brooks et al., 2013).

Instructions were presented as white text on a black screen (‘left hand’, ‘right hand’, ‘face’, or ‘rest’) and viewed by participants through a mirror mounted on the head coil. We further instructed participants to focus on a green fixation cross throughout the scan to minimise eye movements and avoid excessive activation of the abducens nucleus, located caudal to the principle sensory division of the trigeminal nuclei (Paxinos & Mai, 2004). Each movement block lasted 12 seconds, followed by a 6-second rest block, with each condition repeated eight times per run in a counterbalanced order (total run duration = 7 minutes 12 seconds). In total, we acquired four runs, each with a different block order, resulting in 32 trials per condition. Instructions were delivered using Psychtoolbox (Psychophysics Toolbox V.3), which was implemented in Matlab v9.5 (R2018b). Participant head motion was minimised by padded cushions and/or over-ear MRI-safe headphones.

2.3 MRI acquisition

We collected MRI data with a 3 Tesla Philips Ingenia system (Best, The Netherlands) and a 15-channel HeadSpine coil. We acquired a T1-weighted structural scan covering the brain and brainstem to facilitate the co-registration of functional images with anatomical regions of interest (sequence parameters: echo time (TE): 4.4 ms; repetition time (TR): 9.3 ms; flip angle: 8˚; resolution: 0.75 mm3). For task-fMRI, we implemented an echo-planar-imaging (EPI) sequence with partial brain coverage. We recorded blood-oxygenation-level dependent (BOLD) signals from 32 sagittal slices aligned parallel to the floor of the fourth ventricle, covering the brainstem, thalamus, and S1 (sequence parameters: single-shot gradient echo, TE: 26 ms; TR: 2500 ms; flip angle: 82˚; spatial resolution: 1.8 x 1.8 x 1.8 mm3; field of view (FOV): 210 x 187.6 x 70.2 mm3; SENSE factor: 2.1; with 179 volumes per run). Across four runs, we collected 160 volumes per movement condition and 236 volumes for the rest condition. Additionally, we acquired a whole-brain, single-volume EPI image (with parameters identical to the task-fMRI protocol) to enhance co-registration between the fMRI data and anatomical scans. To facilitate post hoc physiological noise correction of the fMRI data, we recorded chest movement data using a respiratory bellow throughout the scans (sampling rate: 496 Hz). A pulse oximeter was attached to the second digit of each participant’s foot to record oxygen saturation and cardiac pulse. All physiological measurement devices were connected to a data acquisition device (Invivo; Phillips Medical Systems, Orlando, Florida) coupled with a desktop computer running the corresponding recording software.

2.4 MRI analysis

Analysis of the MRI data was performed using tools from FSL v5.0.9 (http://fsl.fmrib.ox.ac.uk/fsl), Matlab v9.5 (R2018b), FreeSurfer v6.0 (https://surfer.nmr.mgh.harvard.edu/), SPM12 (http://www.fil.ion.ucl.ac.uk/spm/), the PhysIO toolbox (Kasper et al., 2017), the Decoding Toolbox (TDT) v3.999F (https://sites.google.com/site/tdtdecodingtoolbox/), and Connectome Workbench v1.5 (https://www.humanconnectome.org/software/connectome-workbench).

2.4.1 Preprocessing

As brainstem fMRI signals are highly susceptible to corruption by physiological and other noise sources, we applied an optimised brainstem fMRI preprocessing pipeline based on previous research (Beissner, 2015; Beissner et al., 2011; Faull et al., 2015; Matt et al., 2019; Meissner et al., 2024). This pipeline included brain extraction using FSL’s automated brain extraction tool (BET; Smith, 2002), motion correction using FMRIB’s Linear Image Registration Tool (MCFLIRT; Jenkinson et al., 2002), a 90s high-pass temporal filter implemented via FSL’s Expert Analysis Tool (FEAT; Woolrich et al., 2001), and spatial smoothing. We opted for a 2 mm full width at half maximum (FWHM) Gaussian kernel, as excessive smoothing of brainstem activity is not recommended due to a potential reduction in signal-to-noise ratio caused by averaging signal of interest with surrounding tissue. This kernel size is comparable to the diameter of the target nuclei and aligns with previous recommendations by Beissner et al. (2011).

Further denoising was completed using FSL’s Physiological Noise Modelling (PNM; Brooks et al., 2008). In this step, cardiac and respiratory phases were assigned slicewise to each volume of the EPI image time series. The physiological noise regression model comprised a total of 34 regressors: 4th-order harmonics for both the cardiac and respiratory cycles, 2nd-order harmonics for their interactions, and individual regressors for heart rate and respiration volume. We visually inspected the waveforms to ensure accurate identification of the peaks of the cardiac and respiratory cycles, making adjustments as needed. The resulting voxelwise confound regressors were then incorporated into the first-level general linear model (GLM).

For the multivariate analysis, physiological nuisance regressors were modelled as low-order Fourier series using the RETROICOR method (Glover et al., 2000), implemented in the PhysIO toolbox (Kasper et al., 2017) to allow compatibility with the first-level GLM in SPM12 (http://www.fil.ion.ucl.ac.uk/spm/). The regressors included 3rd-order cardiac harmonics, 4th-order respiratory harmonics, and their 1st-order interaction terms, along with individual regressors for respiratory volume per time (RVT) and heart rate variability, as recommended by Harvey et al. (2008).

Finally, independent component analysis (ICA) was conducted using the Multivariate Exploratory Linear Optimized Decomposition interface (MELODIC; Smith et al., 2004), which enables the decomposition of fMRI data into temporal and spatial components. The resulting components were classified as either “Noise” or “Signal” based on visual inspection of their spatial maps, time series, and power spectra information, following the guidelines published by Griffanti et al. (2017). Two independent researchers labelled the components to assess interrater reliability. Cohen’s kappa indicated strong agreement between raters, with an estimated score of 0.85 (95% CI: 0.82 to 0.88). The ICA-derived noise regressors were subsequently added to the GLM design matrix. The mean temporal signal-to-noise ratio (tSNR) in the brainstem after data preprocessing was 46.53 (s.e.m = 1.30; Supplementary Fig. 1). This is comparable to that reported by Brooks et al. (2013) for 3T brainstem data after noise correction (mean brainstem tSNR = 42).

2.4.2 Image co-registration

Image co-registration was systematically completed in separate, visually inspected steps. First, we aligned each participant’s fMRI data to their whole-brain EPI image using FMRIB’s Linear Image Registration Tool (FLIRT), employing a mutual information cost function and 6 degrees of freedom (Jenkinson & Smith, 2001). Next, each participant’s EPI images were registered to the T1-weighted image, initially using 6 degrees of freedom and the mutual information cost function, and then optimised using boundary-based registration (BBR; Greve & Fischl, 2009). Subsequently, each T1-weighted structural image was co-registered to the MNI152 (1 mm3) template using FMRIB’s nonlinear registration tool (FNIRT) with 12 degrees of freedom (Andersson et al., 2007). Following each step, the image co-registration was visually inspected using blink comparison in Freeview (FreeSurfer) to ensure precise alignment. While standard quality assessments primarily focus on cortical anatomical landmarks (e.g., sulci and gyri), our study additionally prioritised alignment of brainstem structures, particularly the floor of the IVth ventricle, due to our focus on the adjacent dorsal medullary nuclei. Suboptimal alignment was identified as an anterior translation of the floor of the IVth ventricle in the functional image relative to the structural scan. In this instance, we repeated the co-registration without adding the whole-brain EPI image and reinspected the alignment. If the adjustment improved alignment, evidenced by a visually diminished overlap with the cerebrospinal fluid-filled space, we proceeded with the modified registration.

2.4.3 Region of interest definition

To assess hand and face activity in our regions of interest (ROIs) along the somatosensory pathway, we defined masks of the brainstem, thalamus, and postcentral gyrus in MNI 1 mm3 space, as well as surface S1 hand and face area masks.

Whole masks of the right and left postcentral gyrus, bilateral thalamus, and brainstem were derived from the Harvard-Oxford cortical and subcortical structural atlases (Desikan et al., 2006; Frazier et al., 2005). The superior limit of the postcentral gyrus mask was restricted to 1 cm above the anatomical hand knob in S1 (z = 63), and the inferior limit was set to the border of the secondary somatosensory cortex. The thalamus mask was thresholded at 10% to prevent overlap with the ventricles in the MNI 1 mm3 brain template. Additionally, the brainstem mask was thresholded at 50%, and its caudal boundary was manually extended to include the full extent of the cuneate and spinal trigeminal nuclei (z = -71). To reduce the number of voxel-wise comparisons and increase the sensitivity of the univariate non-parametric group analysis, the brainstem mask was further restricted to include only the medulla by removing axial slices above the inferior pontine sulcus (z = -43). Notably, the whole brainstem mask was used for the multivariate classification analysis to avoid a strict a priori hypothesis regarding the location of decoding information and to exploit the covariation of data across the brainstem for noise suppression. Furthermore, the thalamus masks and S1 masks were combined across hemispheres for the multivariate analysis.

For the univariate ROI analysis, specific brainstem ROIs, including the spinal trigeminal (bilateral) and left and right cuneate nuclei, were anatomically defined on the MNI 1 mm3 standard template brain. Notably, we selected the spinal nucleus of the trigeminal nerve as it runs immediately lateral to the cuneate nuclei, providing an ideal model for assessing our research question. We first determined the locations of the right cuneate and spinal trigeminal nuclei using anatomical coordinates extracted from a probabilistic stereotactic brainstem atlas published by Afshar (1978) and further detailed by Beissner et al. (2011). These masks were drawn relative to three anatomical planes: the fastigial floor line, the floor of the IVth ventricle, and the midsagittal plane. For each coordinate, the 90% confidence bounds were calculated and rounded to whole millimetres before drawing. To avoid overlapping ROIs, we excluded voxels that fell within the 90% confidence bounds of the spinal trigeminal and cuneate nuclei. The final volumes of the spinal trigeminal and the cuneate ROIs (unilateral) were 448 mm3 and 296 mm3, respectively. Each ROI was then mirrored along the midsagittal plane to create the contralateral ROI. The cuneate masks were further restricted to below the level of the pyramidal decussation to limit potential signal overlap from both the left and right hands. Additionally, anatomic thalamic masks of the VPL and VPM nuclei were obtained from the Thalamus Morel Atlas (Krauth et al., 2010).

Finally, both S1 hand and face area ROIs were defined on the fs_LR template surface. The hand ROI was defined by selecting the nodes approximately 1 cm below and above the anatomical hand knob within S1, as defined by the Glasser et al. (2016) atlas. The face ROI was defined as starting 1 cm lateral to the S1 hand area and included all S1 lateral to this boundary. A 1 cm gap was maintained between the S1 hand and face ROIs to exclude the portion of S1 where significant overlap between face and hand representations would be expected.

2.5 Univariate analysis

Time-series statistical analysis was performed per run using FEAT with FMRIB’s Improved Linear Model (FILM) with local autocorrelation correction (Woolrich et al., 2001). Activity estimates were derived from a general linear model (GLM), in which regressors were convolved with three basic functions generated from FMRIB’s Linear Optimal Basis Sets (FLOBS) toolkit (Woolrich, Behrens, & Smith et al., 2004). This approach was selected to accommodate the variability of BOLD responses in brainstem regions, which may not be adequately represented by the canonical HRF alone (Wall et al., 2009). While this method convolves the GLM regressors with both the canonical HRF and its temporal and dispersion derivatives, only the canonical HRF regression parameter estimates were used for subsequent higher-level analyses to prevent overestimating statistical effects (as implemented in Meissner et al. (2024) and Faull et al. (2015)). As detailed in Section 2.4.1, voxel-wise confound lists from PNM and ICA noise components were added to the GLM as nuisance regressors. Three main contrasts were defined for the ROI analysis: each movement condition (left hand, right hand, face) versus rest. Additionally, the contrasts left-hand > all, right-hand > all, and face > all were defined for exploring whole-brain activity maps.

A second-level analysis was run for each participant to average across runs using FMRIB’s Local Analysis of Fixed Effects (Woolrich, Behrens, Beckmann, et al., 2004). To minimise the risk of confounding results due to motion (e.g., noise-induced variability, partial volume effects, and motion-induced overlap of border voxels across ROIs), we excluded runs from the first-level analysis with an absolute mean displacement >1 mm. In total, 4 out of the 80 acquired runs were excluded, with no more than one run being omitted per participant.

To visualise group-level activation, we conducted region-specific third-level statistical analyses utilising defined ROIs for the brainstem and thalamus as previously implemented in fMRI analyses of the subcortical regions (Brooks et al., 2017; Oliva et al., 2021, 2022; Reddy et al., 2024). Notably, parametric, cluster-based correction methods, which are typically employed and suitable for widespread cortical activations, lack sensitivity to detect the smaller, more variable clusters in subcortical areas (Friston et al., 1994; Woo et al., 2014). Additionally, brainstem spatial smoothing constraints (as discussed in Section 2.4.1) are incompatible with parametric Gaussian field-based inferences, resulting in overly conservative significance estimates (Beissner et al., 2011; Worsley & Friston, 1995). Therefore, to increase sensitivity to activation in the brainstem and thalamus, we implemented FSL’s nonparametric permutation test RANDOMISE with threshold-free cluster enhancement (TFCE) and family-wise error (FWE) correction (Smith & Nichols, 2009; Winkler et al., 2014). Second-level COPE maps in 1 mm3 MNI standard space were merged across participants for one-sample t-tests within anatomical masks of the medulla and thalamus (left and right combined) separately. The default parameters for variance smoothing (5 mm; as in Reddy et al. (2024) and Kündig et al. (2024)) and the number of permutations (5000) were utilised (Winkler et al., 2014). The resulting TFCE statistical p-maps were thresholded using an FWE correction threshold of p < .05. However, as the complete activation patterns inform our multivariate approach, we utilise the ‘highlighting’ method proposed by Taylor et al. (2023) and additionally display subthreshold activity with a reduced opacity. For statistical inference of cortical activations, second-level statistical analyses were repeated in individual participant space. For optimal cortical alignment, the subsequent COPE maps were transformed to the cortical surface and concatenated for group-level analysis. The cortical surface projections were constructed from each participant’s native T1-weighted images. For each contrast of interest, we implemented FSL’s permutation analysis of linear models (PALM) with sign flipping and 5000 permutations of a two-sided one-sample t-test. The resulting group z-statistic images were thresholded at z > 3.1 and FWE corrected at p < .01.

Finally, to explore the extent to which univariate activity can be distinguished within precise ROIs of the somatosensory relay nuclei, we extracted the task-related activity. This was quantified as the mean z-statistic for each movement condition compared to rest within the anatomical masks of the cuneate and spinal trigeminal nuclei (brainstem), the VPL and VPM (thalamus), and the hand and face cortical areas (S1) (as described in Section 2.4.3). Using repeated-measures ANOVA, we tested the hypothesis that each nucleus would exhibit the highest univariate activity for the somatotopically appropriate body part, for example, activity in the left cuneate nucleus would be highest during left-hand movement compared to right-hand or face movement.

2.6 Multivariate classification analysis

Multivariate classification is a machine-learning technique that decodes different experimental conditions based on their neuronal activation patterns (Mur et al., 2009; Norman et al., 2006; Pereira et al., 2009). Unlike univariate ROI analyses, which provide insight into the strength of voxel-wise activations that are spatially averaged across a specified region, multivariate approaches evaluate simultaneous activity across multiple voxels. In the instance of overlapping activity, as anticipated in the brainstem, this approach may offer greater sensitivity for discriminating between conditions (Ejaz et al., 2015; Muret et al., 2022; Wesselink et al., 2022). Moreover, classification approaches incorporate subthreshold information to provide valuable insight into the underlying representational overlap (Haxby, 2012; Haxby et al., 2001; Haynes & Rees, 2006; Mur et al., 2009; Norman et al., 2006). While univariate approaches may demonstrate that two tasks activate the same voxel, they do not statistically evaluate this shared spatial representation, which is of key importance in assessing changes in somatotopic organisation. Importantly for our research question, multivariate analysis may provide a more comprehensive look at the extent to which we can separate between neural responses of adjacent brainstem nuclei. Therefore, we explored whether we could train a support vector machine (SVM) classifier to decode between left-hand, right-hand, and face activity within the brainstem, thalamus, and S1.

For the classification analysis, we generated beta images for the three main contrasts (left-hand, right-hand, and face vs. rest) using a first-level GLM in SPM12. The settings were matched to those used in the univariate first-level GLM. We set the global masking threshold to zero to ensure that all voxels were retained for the classification analysis regardless of their intensity. As in the univariate analysis, we modelled both the temporal and dispersion derivatives and extracted the canonical HRF regression parameter estimates for classification. This resulted in a total of 240 beta images in 1 mm3 MNI standard space, that is, one beta per condition and per run for each participant. Transforming the first-level beta images to the standard template not only increases the number of samples to train the classifier but also allows us to classify across participants. These beta weights were then used to allow for a multiclass SVM in which we compared the three experimental conditions with each other. To do so, we used the default model parameters of the decoding toolbox (TDT; Hebart, Görgen, & Haynes, 2014) to implement a leave-one-participant-out cross-validation scheme (k = 20, L2-norm multiclass support vector machine with a linear kernel, C = 1). Notably, in contrast to traditional machine learning, the objective here is not to optimise the classifier’s performance but instead to demonstrate that the underlying brain signals contain condition-specific information relevant to classification. We calculated the decoding accuracy of each condition across all the voxels within masks of the brainstem, thalamus, and S1 (see Section 2.4.3 for ROI details). We then extracted the average test classification accuracies and a confusion matrix of the predicted and true labels. Finally, we generated a null distribution from 1000 random permutations of the condition labels. We then calculated the empirical p-value of the average test classification accuracy based on the proportion of permutation accuracies equal to or greater than this score. Furthermore, as decoding of body-part selective activity has not been conducted before within the brainstem, we explored the generalisation of our SVM decoder as a function of sample size and ROI (Supplementary Fig. 2).

Importantly, the above classification analysis does not require a strict apriori hypothesis regarding the location of neural activity, unlike univariate ROI analyses, which are constrained by the need for precise ROI masks to capture activity accurately. Conversely, obtaining significant classification accuracy across the brainstem does not provide insight into whether the classification analysis is sensitive to information encoded in the areas surrounding the target nuclei. Therefore, we also conducted a searchlight analysis to investigate which voxels within the brainstem were associated with classification accuracies above chance. We used a leave-one-participant-out cross-validation scheme (k = 20, L2-norm support vector machine with a linear kernel, C = 1). This analysis produces a searchlight map that assigns multiclass classification accuracy to each voxel within the brainstem ROI. Notably, across-participant searchlight analyses assume that the spatial variation among participants is smaller than the searchlight radius and may otherwise be prone to spatial distortions or exaggeration (Etzel et al., 2013). For optimal sampling, Kriegeskorte et al. (2006) recommend a searchlight radius comparable to the size and compactness of the region. To balance these factors, we implemented a 3 mm searchlight radius, which is comparable to the diameter of the targeted nuclei and may reduce the risk of spatial inaccuracies if the information is present in clusters within these nuclei (Afshar, 1978; Etzel et al., 2013; Knights et al., 2022).

To identify voxels with classification accuracies significantly above chance, we employed the statistical approach recommended by Stelzer et al. (2013). First, we ran 1000 permutations with random shuffling of the condition labels to create a null distribution of classification accuracies across all voxels in the brainstem. We then determined the classification accuracy corresponding to a critical cut-off of 0.05, thresholding both the initial searchlight map and the 1000 permutation maps by this value. Next, we conducted cluster correction by extracting the maximum cluster size within the thresholded permutation maps to create a null distribution of maximum cluster sizes. Finally, we calculated the cluster size corresponding to a p-value of 0.05 and masked the thresholded searchlight map to only include clusters larger than this value.

2.7 Statistical analysis

Statistical analysis for the univariate ROI analysis was conducted using RStudio (v4.0.3) and GraphPad Prism (v10.3.0). Outlier values ± 2.5 standard deviations from the mean were replaced with the mean. All statistical testing was conducted with and without outliers to confirm that this adjustment did not change the results. We used standard approaches for statistical analysis, as detailed in the Results section. Parametric or non-parametric testing was used depending on normality, which was confirmed with the Shapiro-Wilk test. In addition, Mauchly’s test was implemented to assess sphericity across each ROI, and Bonferroni or Dunn’s (Friedmans ANOVA) corrections were applied for post-hoc comparisons. For the univariate ROI analysis, planned contrasts were implemented to compare the mean of the target condition with each non-target condition.

We investigated whether 3T fMRI, with an optimised acquisition and preprocessing pipeline, could effectively distinguish between neural responses of neighbouring somatosensory nuclei in the brainstem, thalamus, and S1. We delineated left-hand, right-hand, and face movement-related activity utilising two approaches. First, using traditional univariate voxel-wise and ROI analyses, we assessed the extent to which we could distinguish between conditions based on their overall activation strength. Second, using multivariate classification analysis, we explored the degree to which somatosensory nuclei can be distinguished based on their underlying activation patterns for the different conditions, thereby providing greater insight into the representational overlap of these signals. In addition, we conducted a multiclass searchlight analysis to determine the specific regions of the brainstem encoding relevant information for the classification of hand and face movement.

3.1 Univariate analysis reveals a partial differentiation of hand and face activation within somatosensory relay nuclei

We first explored body-part-specific activity across the somatosensory processing stream. To do so, we created group activity maps for each condition contrasted with all other conditions. We found significant group-level activation in the anticipated nuclei and areas of the brainstem, thalamus, and S1 (Fig. 2a).

Fig. 2.

Univariate analyses can partially differentiate hand and face activity within neighbouring somatosensory nuclei. (A) Group-level activity maps (n = 20) during left-hand, right-hand, and face movements. To obtain body-part-specific activity, we contrasted the conditions as follows: left-hand > right-hand and face (purple), right-hand > left-hand and face (turquoise), and face > left and right hands (yellow). Anatomical regions of interest for the somatosensory processing stream are displayed in the left column, demonstrating anticipated areas of activation (left cuneate nucleus = purple; right cuneate nucleus = turquoise; spinal trigeminal nucleus = yellow). The corresponding family-wise error-corrected (PFWE) activity maps are shown in the right column. Cortical activations are projected onto a cortical surface, thresholded at Z > 3.1 and PFWE < 0.01; subthreshold cortical activity is omitted for clarity. For the brainstem and thalamus, suprathreshold activity clusters are outlined in black (threshold-free cluster enhancement, PFWE < 0.05). Subthreshold activity is shown with reduced opacity. PFWE values are indicated alongside their respective colour bars. (B) Region of interest (ROI) analysis shows the mean z-statistic for each task (contrasted with rest) within anatomical masks of the left cuneate, right cuneate, and spinal trigeminal nuclei (visualised in A, left panel). Individual participants are plotted as dots. Statistical significance is denoted by: ****p < .0001, ***p < .001, **p < .01. p-Values are corrected for multiple comparisons within each ROI.

Fig. 2.

Univariate analyses can partially differentiate hand and face activity within neighbouring somatosensory nuclei. (A) Group-level activity maps (n = 20) during left-hand, right-hand, and face movements. To obtain body-part-specific activity, we contrasted the conditions as follows: left-hand > right-hand and face (purple), right-hand > left-hand and face (turquoise), and face > left and right hands (yellow). Anatomical regions of interest for the somatosensory processing stream are displayed in the left column, demonstrating anticipated areas of activation (left cuneate nucleus = purple; right cuneate nucleus = turquoise; spinal trigeminal nucleus = yellow). The corresponding family-wise error-corrected (PFWE) activity maps are shown in the right column. Cortical activations are projected onto a cortical surface, thresholded at Z > 3.1 and PFWE < 0.01; subthreshold cortical activity is omitted for clarity. For the brainstem and thalamus, suprathreshold activity clusters are outlined in black (threshold-free cluster enhancement, PFWE < 0.05). Subthreshold activity is shown with reduced opacity. PFWE values are indicated alongside their respective colour bars. (B) Region of interest (ROI) analysis shows the mean z-statistic for each task (contrasted with rest) within anatomical masks of the left cuneate, right cuneate, and spinal trigeminal nuclei (visualised in A, left panel). Individual participants are plotted as dots. Statistical significance is denoted by: ****p < .0001, ***p < .001, **p < .01. p-Values are corrected for multiple comparisons within each ROI.

Close modal

At the cortical level, we found significant task-selective brain activity during left- and right-hand movement in the contralateral S1 and primary motor (M1) hand cortices (Fig. 2a; top panel). We further detected significant bilateral activity elicited by face movement laterally to the hand representations, consistent with the established somatotopic organisation of these body parts in the cortex (Penfield & Boldrey, 1937; Penfield & Rasmussen, 1950; Roux et al., 2018). To enhance sensitivity for detecting subcortical activation, we conducted region-based group-level analyses within anatomically defined masks of the medulla and thalamus. As expected, left- and right-hand movements significantly activated the contralateral VPL nucleus of the thalamus (Fig 2a; middle panel). In contrast, face movement did not yield significant activation in the right thalamus. However, uncorrected TFCE p-maps revealed subthreshold face-task-selective activation medial to the hand clusters, overlapping with the VPM nuclei. In the medulla, left-hand movement elicited significant activation clusters overlapping with the left cuneate nucleus. While the right-hand movement did not produce significant clusters, subthreshold activation was localised within the right cuneate nucleus (Fig. 2a; bottom panel). Face movements significantly activated the spinal trigeminal nuclei bilaterally, as expected. Notably, this activation was widespread, with one significant cluster located between the cuneate nuclei.

Next, we aimed to explore the extent to which univariate activity can be distinguished within precise ROIs of the somatosensory relay nuclei. To do so, we extracted the task-related activity, quantified as the mean z-statistic, for each movement condition compared to rest within anatomical masks of the cuneate and spinal trigeminal nuclei (brainstem), the VPL and VPM (thalamus), hand and face cortical areas (S1) (Fig. 2b). During movement of the hands, we anticipated higher activation in the ipsilateral cuneate, contralateral VPL, and contralateral S1 hand area, as compared to the opposite hand and face movement. Conversely, during movement of the face, we anticipated higher bilateral activation within the spinal trigeminal, VPM and S1 face areas, as compared to left or right-hand movement.

The statistical results of the ROI-specific comparisons and the number of outlier values are summarised in Table 1. In S1 (Fig. 2b; top panel), we reliably distinguished between the hand and face movement conditions across all ROIs. Similarly, in the thalamus (Fig. 2b; middle panel), we were able to distinguish between hand and face movement within the left and right VPL ROIs. However, in the VPM ROI, we only found a trend towards higher face than left-hand activity. Finally, in the brainstem (Fig. 2b; bottom panel), while we could distinguish between ipsilateral and contralateral hand activity in the cuneate ROIs, we could not distinguish between the ipsilateral hand and adjacent face activity. These findings suggest a limited ability to differentiate between hand and facial activity within the anatomical confines of our cuneate ROIs using univariate analysis. However, in support of our hypothesis, we found that face activity was significantly higher than both left- and right-hand activity within the spinal trigeminal ROI. Due to the apparent left-side dominance in the group activity maps for the brainstem hand activity, we compared the mean left-hand activity in the left cuneate ROI to the mean right-hand activity in the right cuneate ROI. However, we found there was no significant difference between the activity levels for the two hands (t(19) = 1.08, p = .30, Cohen’s d = 0.24). Additionally, to assess the power of our study to detect observed effects in the brainstem, we conducted a post hoc sensitivity analysis (see Supplementary Materials for details). This analysis suggested that the sample size was adequate for detecting the observed effect sizes in the spinal trigeminal and left cuneate ROIs, though a larger sample may have been beneficial for the right cuneate ROI analysis.

Table 1.

Summary of statistical effects for the univariate ROI analyses

ROIComparisonTest statisticp-ValueEffect sizeOutliers
S1 
Left S1 Hand Area Main effect F(2, 38) = 303.4 <.0001 η² = 0.941 
Right vs. Left hand t(38) = 19.74 <.0001  
Right hand vs. Face t(38) = 22.63 <.0001  
Right S1 Hand Area Main effect χ²(2) = 30.40 <.0001 W = 0.760 
Left vs. Right hand Z = 4.427 <.0001  
Left hand vs. Face Z = 5.060 <.0001  
Bilateral S1 Face Area Main effect F(2, 38) = 54.64 <.0001 η² = 0.742 
Face vs. Left hand t(38) = 7.459 <.0001  
Face vs. Right hand t(38) = 7.446 <.0001  
Thalamus 
Left VPL Main effect F(2, 38) = 52.63 <.0001 η² = 0.735 
Right vs. Left hand t(38) = 10.08 <.0001  
Right hand vs. Face t(38) = 6.687 <.0001  
Right VPL Main effect F(2, 38) = 36.13 <.0001 η² = 0.655 
Left vs. Right hand t(38) = 7.841 <.0001  
Left hand vs. Face t(38) = 6.763 <.0001  
Bilateral VPM Main effect F(2, 38) = 10.54 <.001 η² = 0.357 
Face vs. Right hand t(38) = 4.582 <.0001  
Face vs. Left hand t(38) = 2.052 0.094  
Brainstem 
Left Cuneate Main effect F(2, 38) = 24.40 <.0001 η² = 0.562 
Left vs. Right hand t(38) = 6.079 <.0001  
Left hand vs. Face t(38) = 0.061 >.9999  
Right Cuneate Main effect χ²(2) = 9.70 .008 W = 0.243 
Right vs. Left hand Z = 3.004 .0053  
Right hand vs. Face Z = 0.791 .858  
Spinal Trigeminal Main effect F(2, 38) = 32.25 <.0001 η² = 0.629 
Face vs. Left hand t(38) = 6.708 <.0001  
Face vs. Right hand t(38) = 7.179 <.0001  
ROIComparisonTest statisticp-ValueEffect sizeOutliers
S1 
Left S1 Hand Area Main effect F(2, 38) = 303.4 <.0001 η² = 0.941 
Right vs. Left hand t(38) = 19.74 <.0001  
Right hand vs. Face t(38) = 22.63 <.0001  
Right S1 Hand Area Main effect χ²(2) = 30.40 <.0001 W = 0.760 
Left vs. Right hand Z = 4.427 <.0001  
Left hand vs. Face Z = 5.060 <.0001  
Bilateral S1 Face Area Main effect F(2, 38) = 54.64 <.0001 η² = 0.742 
Face vs. Left hand t(38) = 7.459 <.0001  
Face vs. Right hand t(38) = 7.446 <.0001  
Thalamus 
Left VPL Main effect F(2, 38) = 52.63 <.0001 η² = 0.735 
Right vs. Left hand t(38) = 10.08 <.0001  
Right hand vs. Face t(38) = 6.687 <.0001  
Right VPL Main effect F(2, 38) = 36.13 <.0001 η² = 0.655 
Left vs. Right hand t(38) = 7.841 <.0001  
Left hand vs. Face t(38) = 6.763 <.0001  
Bilateral VPM Main effect F(2, 38) = 10.54 <.001 η² = 0.357 
Face vs. Right hand t(38) = 4.582 <.0001  
Face vs. Left hand t(38) = 2.052 0.094  
Brainstem 
Left Cuneate Main effect F(2, 38) = 24.40 <.0001 η² = 0.562 
Left vs. Right hand t(38) = 6.079 <.0001  
Left hand vs. Face t(38) = 0.061 >.9999  
Right Cuneate Main effect χ²(2) = 9.70 .008 W = 0.243 
Right vs. Left hand Z = 3.004 .0053  
Right hand vs. Face Z = 0.791 .858  
Spinal Trigeminal Main effect F(2, 38) = 32.25 <.0001 η² = 0.629 
Face vs. Left hand t(38) = 6.708 <.0001  
Face vs. Right hand t(38) = 7.179 <.0001  

Analyses were conducted using repeated-measures ANOVA or Friedman’s ANOVA.

Overall, these results indicate that activation among somatosensory nuclei can be distinguished to some extent using univariate analysis of 3T fMRI data. However, this distinction varies across conditions and locations. Specifically, our brainstem ROI analysis of the cuneate nuclei did not clearly differentiate between hand and facial activation. The anatomical definition of our ROIs was designed so that the cuneate and spinal trigeminal nuclei were adjacent, with no intervening gap. This approach is constrained by the need for precise ROI masks to accurately capture activity. Therefore, to examine the extent to which we can distinguish between conditions without being limited by this potential confound, we next conducted a multivariate classification analysis.

3.2 Multivariate classification can differentiate between the activation of neighbouring brainstem nuclei involved in somatosensory processing

To investigate the extent to which left-hand, right-hand, and face movement can be distinguished through their neuronal activity patterns, we conducted a classification analysis. Specifically, we tested whether we could train a multiclass support vector machine (SVM) to decode left-hand, right-hand, and face activation across participants within bilateral anatomical masks of the brainstem, thalamus, and S1. As anticipated, we found a high classification accuracy in S1, where fMRI signals are robust and less prone to noise corruption, with an average multiclass classification accuracy across participants of 98.3% (s.e.m. = 1.30, p = <.001; empirical chance level: 33.3%). Moreover, a confusion matrix comparing the frequency of predicted to true labels did not display a bias towards an individual movement condition (Fig. 3). Importantly, we found that condition-specific activity could also be distinguished significantly above chance level in both the thalamus (accuracy = 60.4%, s.e.m. = 2.76, p = <.001) and the brainstem (accuracy = 50.8%, s.e.m. = 2.25, p = <.001) (Fig. 3). The confusion matrices revealed that we could distinguish between all conditions in both these subcortical areas. Post hoc analysis showed that brainstem classification accuracy plateaued at n = 10 participants (Supplementary Fig. 2), indicating that the current sample size was sufficient to achieve stable classification performance. Additionally, to confirm that ipsilateral hand and face classification was not solely driven by contralateral face activity, we replicated these findings using masks of the left and right brainstem, divided along the midsagittal plane (x = 0) (Supplementary Fig. 3). These results indicate that there is sufficient information present in 3T fMRI data subcortically to decode between neighbouring somatosensory representations of the face and hands.

Fig. 3.

Hand and face activation patterns can be distinguished across the somatosensory processing stream using multivariate classification analysis. We trained linear support vector machines to classify voxel-wise activity patterns across participants in three anatomical regions of interest (ROIs): (A) primary somatosensory cortex (S1), (B) thalamus, and (C) brainstem. Each confusion matrix depicts the percentage of correct and incorrect classifications for the left-hand, right-hand, and face conditions (rows = true labels, columns = predicted labels).

Fig. 3.

Hand and face activation patterns can be distinguished across the somatosensory processing stream using multivariate classification analysis. We trained linear support vector machines to classify voxel-wise activity patterns across participants in three anatomical regions of interest (ROIs): (A) primary somatosensory cortex (S1), (B) thalamus, and (C) brainstem. Each confusion matrix depicts the percentage of correct and incorrect classifications for the left-hand, right-hand, and face conditions (rows = true labels, columns = predicted labels).

Close modal

While the classification analysis described above shows that we can reliably decode between conditions based on brainstem activity patterns, it does not specify which regions of the brainstem enable our classifier to differentiate between conditions. We, therefore, conducted a searchlight analysis to investigate which voxels within the brainstem were associated with classification accuracies above chance level. The resulting accuracy map (Fig. 4) was thresholded by a critical cut-off of 39.7% accuracy (p < .05) and cluster-corrected (k = 639 voxels; p < .05). As anticipated, the voxels with significant classification accuracies were localised in the dorsal medulla, overlapping with the anatomical masks of our target nuclei: the cuneate and spinal trigeminal nuclei. More extensive significant clusters were observed in the left dorsal medulla. Significant voxels were also located in the dorsal pons surrounding the trigeminal motor and principal sensory nuclei. This result demonstrates that using multivariate classification analysis, 3T fMRI, and a body movement task, it is possible to distinguish between the neuronal activity patterns of adjacent brainstem nuclei involved in somatosensory processing.

Fig. 4.

Searchlight analysis across the brainstem demonstrates that voxels differentiating between face and hand movement are localised in the dorsal medulla. We conducted a searchlight analysis to identify the contribution of individual brainstem voxels to the classification of hand and face movements. The left panel demonstrates the brainstem region of interest (white) implemented in the analysis. The expected areas of neural activation for somatosensory input from the face (spinal trigeminal nucleus; yellow) and hands (cuneate nucleus; purple) are overlaid. The right panel displays the searchlight accuracy map illustrating voxels with classification accuracies significantly above chance (thresholded and cluster-corrected at a critical cut-off of 39.7%, p = < .05; k = 639 voxels). The results indicate that the activation patterns for hand and face movement are predominately located in the dorsal medulla, overlapping with the anticipated nuclei involved in the somatosensory processing of face and hand input.

Fig. 4.

Searchlight analysis across the brainstem demonstrates that voxels differentiating between face and hand movement are localised in the dorsal medulla. We conducted a searchlight analysis to identify the contribution of individual brainstem voxels to the classification of hand and face movements. The left panel demonstrates the brainstem region of interest (white) implemented in the analysis. The expected areas of neural activation for somatosensory input from the face (spinal trigeminal nucleus; yellow) and hands (cuneate nucleus; purple) are overlaid. The right panel displays the searchlight accuracy map illustrating voxels with classification accuracies significantly above chance (thresholded and cluster-corrected at a critical cut-off of 39.7%, p = < .05; k = 639 voxels). The results indicate that the activation patterns for hand and face movement are predominately located in the dorsal medulla, overlapping with the anticipated nuclei involved in the somatosensory processing of face and hand input.

Close modal

Functional mapping of the human brainstem presents many challenges. These methodological constraints have deterred essential investigations into human sensory brainstem processing. Here, we delineate the functional activity of two adjacent somatosensory nuclei, namely the spinal trigeminal and cuneate nuclei. Our results demonstrate that functional activity can be distinguished in humans across the somatosensory stream through fMRI and optimised data acquisition and analysis strategies. This is a fundamental step towards broadening our understanding of somatosensory processing within humans and determining what changes in organisation may occur following sensory deprivation.

We utilised both univariate and multivariate analyses to delineate the activity of neighbouring somatosensory nuclei. Notably, in the univariate group analysis, neither right-hand activity clusters in the brainstem nor face activity clusters in the right thalamus reached the significance threshold. A previous study by Reddy et al. (2024) demonstrated significant activation in the brainstem during both hand and foot stimulation using fMRI. However, while Reddy et al. contrasted each body part with rest, we compared each body part to all other conditions. Therefore, instead of simply identifying clusters associated with each body part, we sought body-part selective clusters to assess how effectively we could differentiate between tasks. Importantly, the uncorrected group activity maps revealed that, although weak, face and right-hand selective clusters spatially overlapped with the anticipated nuclei involved in somatosensory processing. This finding highlights a limitation of univariate analyses, which are constrained by the chosen significance threshold, whereas multivariate analyses can combine subthreshold information across voxels to enhance signals that may otherwise go undetected (Hebart & Baker, 2018).

Similarly, in the univariate ROI analysis, while we were able to distinguish between all conditions in S1, we were only able to partially differentiate between conditions in the brainstem and thalamus. Given that atlas-based masks do not draw a distinct border between the hand and face nuclei in these regions, even slight anatomical differences between individuals could cause inaccuracies in estimating these borders, leaving no margin for error. Moreover, the normalisation process during co-registration may have caused a slight misalignment of the anatomical borders in individual participants, causing an artificial overlap in activity. As with the S1 analysis, it would be beneficial to leave a border of empty voxels between ROIs to minimise this potential overlap. However, the small size of the brainstem nuclei and the resolution limits of 3T imaging make such an approach unfeasible. Higher-field imaging, such as 7T fMRI, may provide the spatial resolution needed to define a precise boundary between nuclei. However, it may also increase susceptibility to artefacts and signal inhomogeneities (Smith et al., 2014). Alternatively, the inability to separate hand and face activity may reflect a (partially) shared somatotopic overlap between these representations within the cuneate nucleus, as previously observed in the human primary sensorimotor cortex (Gordon et al., 2023; Muret & Makin, 2021; Muret et al., 2022; Schellekens et al., 2018). Indeed, separation between the two nuclei has proved unsuccessful in animal research, where direct injection of tracers has led to the labelling of the neighbouring nucleus (Jain et al., 2000; Liao et al., 2021; Marfurt & Rajchert, 1991). Therefore, analysis techniques that do not rely on such a strict a priori region specification may be more appropriate for differentiating between these potentially overlapping somatosensory representations.

Multivariate classification analyses provide a solution to evaluating neuronal representations across conditions within the entire somatosensory system. Notably, this technique offers greater sensitivity and specificity than univariate analysis by targeting distributed patterns of activity and exploiting the covariance of data across voxels (Haynes & Rees, 2005; Hebart & Baker, 2018; Jimura & Poldrack, 2012; Kamitani & Tong, 2005). More importantly, this method may provide a more comprehensive understanding of how well we can separate neural responses of adjacent brainstem nuclei, even when the signal-to-noise ratio is low (Haufe et al., 2014; Hebart & Baker, 2018). Indeed, through this approach, we found that we could differentiate between the neural representations of left-hand, right-hand, and face movement across the brainstem, thalamus, and S1. To infer which areas of the brainstem encoded information for distinguishing these tasks, we also conducted a searchlight analysis. We found that the hand and face movement conditions were specifically encoded in the dorsal medulla, overlapping with the anticipated nuclei involved in somatosensory processing. Notably, we observed a left-sided dominance in the brainstem for both the univariate and searchlight analyses. However, quantitatively, we did not find a statistically significant difference in activity between the hands. Therefore, we are cautious not to overinterpret this pattern. Future studies are needed to explore whether larger and more heterogeneous cohorts may provide more robust evidence of lateralisation. Notably, our study cohort was relatively homogeneous in age. Age-related changes in brain structure and function, such as reduced neuroplasticity and altered functional connectivity, may influence sensorimotor processing and neural activation patterns (Bardouille & Bailey, 2019; Cassady et al., 2019; Seidler et al., 2010). Therefore, further research specifically targeting different age groups is needed to assess whether similar patterns are observed across the lifespan.

To our knowledge, multivariate classification has not yet been applied to distinguish body-part selective activity within the brainstem. Our findings demonstrate a novel solution for exploring signal processing within the brainstem and across the entire somatosensory system. Moreover, given the brainstem’s diverse roles in pain modulation (Youssef et al., 2016), auditory processing (De Martino et al., 2013), visual perception (Limbrick-Oldfield et al., 2012), respiration (Faull et al., 2015), emotion processing (Satpute et al., 2013), behavioural arousal (Aston-Jones et al., 1999), and autonomic control (Samuels & Szabadi, 2008), we suggest that multivariate approaches may allow crucial insights into a broad range of physiological processes. Importantly, in our classification analysis, we utilised an across-participant decoding design as opposed to the approach that is commonly used for cortical classification analysis, where inferences are made by combining within-participant accuracy scores (Gallivan et al., 2011; Haxby et al., 2001; Kamitani & Tong, 2005; Nambu et al., 2015; Odermatt et al., 2024). In regions with a low SNR, like the brainstem, aggregating data across participants may provide more robust results and has been demonstrated to detect smaller effects if the inter-individual variability remains moderate (Wang et al., 2020). Critically, this method enables us to make inferences about somatosensory processing differences between patient and control populations. Therefore, future studies may extend this analysis to investigate changes in somatotopic organisation between healthy and patient cohorts that experience sensory deprivation, for example, spinal cord injury, limb loss, or peripheral neuropathy. Understanding the extent to which neural representations of sensory stimuli are preserved or altered following major sensory deprivation may assist in the development of therapeutic interventions that aim to promote adaptive plasticity and prevent the development of maladaptive circuitry.

Together, our results demonstrate that sensory activation in adjacent brainstem nuclei can be distinguished in healthy individuals using a movement paradigm and 3T fMRI with optimised data acquisition and analysis strategies. This provides the methodological foundation for future investigations into a wide range of brainstem-mediated functions, including the plasticity or stability of somatosensory processing in clinical populations.

Data analysis was performed using publicly available toolboxes and script templates, with modifications as described in the Methods section of the manuscript. Researchers wishing to replicate or build on this work are encouraged to contact the corresponding author for further guidance. The raw data analysed in this study cannot be shared due to confidentiality restrictions outlined in the ethics agreement approved by the Ethics Committee of the Canton of Zürich (KEK-2018-00937). However, aggregated and anonymised group-level data will be made available to support reproducibility.

P.H., N.W., S.M., and S.K. conceptualised and designed the study. P.H., S.K., and I.O. programmed the task and analysis scripts. P.H., S.K., F.R., and I.O. acquired the data. P.H., S.K., O.H., P.F., and N.W. planned the analysis. P.H. analysed the data. P.H., S.K., P.F., O.H., and N.W. interpreted the data. P.H. and S.K. drafted the manuscript, and all authors substantively revised it.

We thank the participants for their participation in the study. We thank Sijamini Baskaralingam and Charlotte Meninghin for their help with data collection. We thank Roger Luchinger and the Institute for Biomedical Imaging for scanning support. During this project, P.H. and S.K. were funded by the Swiss National Science Foundation Ambizione Grant PZ00P3_208996 and Grant 32003B_207719. N.W. is additionally supported by the National Research Foundation, Prime Minister’s Office, Singapore under its Campus for Research Excellence and Technological Enterprise (CREATE) program (FHT). O.H. is supported by the Royal Society of New Zealand and New Zealand Lottery Health Research. S.M. is supported by the ETH Career Seed Award.

The authors declare no competing interest.

Supplementary material for this article is available with the online version here: https://doi.org/10.1162/imag_a_00581.

Afshar
,
F
.
,
Watkins
,
E. S
.
, &
Yap
,
J. C
. (
1978
).
Stereotaxic atlas of the human brainstem and cerebellar nuclei: A variability study
.
Raven Press
.
Andersson
,
J. L. R.
,
Jenkinson
,
M.
,
&
Smith
,
S
.
(
2007
).
Non-linear registration, aka spatial normalisation
(FMRIB Technical Report TR07JA2). FMRIB Centre, University of Oxford. https://www.fmrib.ox.ac.uk/analysis/techrep
Aston-Jones
,
G.
,
Rajkowski
,
J.
, &
Cohen
,
J.
(
1999
).
Role of locus coeruleus in attention and behavioral flexibility
.
Biol Psychiatry
,
46
(
9
),
1309
1320
. https://doi.org/10.1016/s0006-3223(99)00140-7
Bardouille
,
T.
, &
Bailey
,
L.
(
2019
).
Evidence for age-related changes in sensorimotor neuromagnetic responses during cued button pressing in a large open-access dataset
.
Neuroimage
,
193
,
25
34
. https://doi.org/10.1016/j.neuroimage.2019.02.065
Beissner
,
F.
(
2015
).
Functional MRI of the brainstem: Common problems and their solutions
.
Clin Neuroradiol
,
25
(
Suppl 2
),
251
257
. https://doi.org/10.1007/s00062-015-0404-0
Beissner
,
F.
,
Deichmann
,
R.
, &
Baudrexel
,
S.
(
2011
).
fMRI of the brainstem using dual-echo EPI
.
Neuroimage
,
55
(
4
),
1593
1599
. https://doi.org/10.1016/j.neuroimage.2011.01.042
Beissner
,
F.
,
Schumann
,
A.
,
Brunn
,
F.
,
Eisenträger
,
D.
, &
Bär
,
K. J.
(
2014
).
Advances in functional magnetic resonance imaging of the human brainstem
.
Neuroimage
,
86
,
91
98
. https://doi.org/10.1016/j.neuroimage.2013.07.081
Besle
,
J.
,
Sánchez-Panchuelo
,
R.-M.
,
Bowtell
,
R.
,
Francis
,
S.
, &
Schluppeck
,
D.
(
2013
).
Single-subject fMRI mapping at 7 T of the representation of fingertips in S1: A comparison of event-related and phase-encoding designs
.
J Neurophysiol
,
109
(
9
),
2293
2305
. https://doi.org/10.1152/jn.00499.2012
Brooks
,
J. C.
,
Beckmann
,
C. F.
,
Miller
,
K. L.
,
Wise
,
R. G.
,
Porro
,
C. A.
,
Tracey
,
I.
, &
Jenkinson
,
M.
(
2008
).
Physiological noise modelling for spinal functional magnetic resonance imaging studies
.
Neuroimage
,
39
(
2
),
680
692
. https://doi.org/10.1016/j.neuroimage.2007.09.018
Brooks
,
J. C.
,
Faull
,
O. K.
,
Pattinson
,
K. T.
, &
Jenkinson
,
M.
(
2013
).
Physiological noise in brainstem FMRI
.
Front Hum Neurosci
,
7
,
623
. https://doi.org/10.3389/fnhum.2013.00623
Brooks
,
J. C. W.
,
Davies
,
W.-E.
, &
Pickering
,
A. E.
(
2017
).
Resolving the Brainstem contributions to attentional analgesia
.
J Neurosci
,
37
(
9
),
2279
2291
. https://doi.org/10.1523/jneurosci.2193-16.2016
Cassady
,
K.
,
Gagnon
,
H.
,
Lalwani
,
P.
,
Simmonite
,
M.
,
Foerster
,
B.
,
Park
,
D.
,
Peltier
,
S. J.
,
Petrou
,
M.
,
Taylor
,
S. F.
,
Weissman
,
D. H.
,
Seidler
,
R. D.
, &
Polk
,
T. A.
(
2019
).
Sensorimotor network segregation declines with age and is linked to GABA and to sensorimotor performance
.
Neuroimage
,
186
,
234
244
. https://doi.org/10.1016/j.neuroimage.2018.11.008
Charyasz
,
E.
,
Heule
,
R.
,
Molla
,
F.
,
Erb
,
M.
,
Kumar
,
V. J.
,
Grodd
,
W.
,
Scheffler
,
K.
, &
Bause
,
J.
(
2023
).
Functional mapping of sensorimotor activation in the human thalamus at 9.4 Tesla
.
Front Neurosci
,
17
,
1116002
. https://doi.org/10.3389/fnins.2023.1116002
Craven
,
J
.
(
2011
).
Major ascending and descending tracts in the spinal cord
.
Anaesthesia & Intensive Care Medicine
,
12
(
1
),
26
27
. https://doi.org/10.1016/j.mpaic.2010.10.001
De Martino
,
F.
,
Moerel
,
M.
,
van de Moortele
,
P. F.
,
Ugurbil
,
K.
,
Goebel
,
R.
,
Yacoub
,
E.
, &
Formisano
,
E.
(
2013
).
Spatial organization of frequency preference and selectivity in the human inferior colliculus
.
Nat Commun
,
4
,
1386
. https://doi.org/10.1038/ncomms2379
Desikan
,
R. S.
,
Ségonne
,
F.
,
Fischl
,
B.
,
Quinn
,
B. T.
,
Dickerson
,
B. C.
,
Blacker
,
D.
,
Buckner
,
R. L.
,
Dale
,
A. M.
,
Maguire
,
R. P.
,
Hyman
,
B. T.
,
Albert
,
M. S.
, &
Killiany
,
R. J.
(
2006
).
An automated labeling system for subdividing the human cerebral cortex on MRI scans into gyral based regions of interest
.
Neuroimage
,
31
(
3
),
968
980
. https://doi.org/10.1016/j.neuroimage.2006.01.021
Ejaz
,
N.
,
Hamada
,
M.
, &
Diedrichsen
,
J.
(
2015
).
Hand use predicts the structure of representations in sensorimotor cortex
.
Nat Neurosci
,
18
(
7
),
1034
1040
. https://doi.org/10.1038/nn.4038
Etzel
,
J. A.
,
Zacks
,
J. M.
, &
Braver
,
T. S.
(
2013
).
Searchlight analysis: Promise, pitfalls, and potential
.
Neuroimage
,
78
,
261
269
. https://doi.org/10.1016/j.neuroimage.2013.03.041
Faul
,
F.
,
Erdfelder
,
E.
,
Lang
,
A.-G.
, &
Buchner
,
A.
(
2007
).
G*Power 3: A flexible statistical power analysis program for the social, behavioral, and biomedical sciences
.
Behav Res Methods
,
39
(
2
),
175
191
. https://doi.org/10.3758/BF03193146
Faull
,
O. K.
,
Jenkinson
,
M.
,
Clare
,
S.
, &
Pattinson
,
K. T. S.
(
2015
).
Functional subdivision of the human periaqueductal grey in respiratory control using 7tesla fMRI
.
Neuroimage
, 113,
356
364
. https://doi.org/10.1016/j.neuroimage.2015.02.026
Frazier
,
J. A.
,
Chiu
,
S.
,
Breeze
,
J. L.
,
Makris
,
N.
,
Lange
,
N.
,
Kennedy
,
D. N.
,
Herbert
,
M. R.
,
Bent
,
E. K.
,
Koneru
,
V. K.
,
Dieterich
,
M. E.
,
Hodge
,
S. M.
,
Rauch
,
S. L.
,
Grant
,
P. E.
,
Cohen
,
B. M.
,
Seidman
,
L. J.
,
Caviness
,
V. S.
, &
Biederman
,
J.
(
2005
).
Structural brain magnetic resonance imaging of limbic and thalamic volumes in pediatric bipolar disorder
.
Am J Psychiatry
,
162
(
7
),
1256
1265
. https://doi.org/10.1176/appi.ajp.162.7.1256
Friston
,
K. J
.
,
Worsley
,
K. J
.
,
Frackowiak
,
R. S
.
,
Mazziotta
,
J. C
.
, &
Evans
,
A. C
.
(
1994
).
Assessing the significance of focal activations using their spatial extent
.
Hum Brain Mapp
,
1
(
3
),
210
220
. https://doi.org/10.1002/hbm.460010306
Gallivan
,
J. P.
,
McLean
,
D. A.
,
Valyear
,
K. F.
,
Pettypiece
,
C. E.
, &
Culham
,
J. C.
(
2011
).
Decoding action intentions from preparatory brain activity in human parieto-frontal networks
.
J Neurosci
,
31
(
26
),
9599
9610
. https://doi.org/10.1523/jneurosci.0080-11.2011
Glasser
,
M. F.
,
Coalson
,
T. S
.
,
Robinson
,
E. C.
,
Hacker
,
C. D.
,
Harwell
,
J.
,
Yacoub
,
E
.
,
Ugurbil
,
K
.
,
Andersson
,
J.
,
Beckmann
,
C. F.
,
Jenkinson
,
M.
,
Smith
,
S. M.
, &
Van Essen
,
D. C.
(
2016
).
A multi-modal parcellation of human cerebral cortex
.
Nature
,
536
(
7615
),
171
178
. https://doi.org/10.1038/nature18933
Glover
,
G. H.
(
2011
).
Overview of functional magnetic resonance imaging
.
Neurosurg Clin N Am
,
22
(
2
),
133
139
, vii. https://doi.org/10.1016/j.nec.2010.11.001
Glover
,
G. H.
,
Li
,
T. Q.
, &
Ress
,
D.
(
2000
).
Image-based method for retrospective correction of physiological motion effects in fMRI: RETROICOR
.
Magn Reson Med
,
44
(
1
),
162
167
. https://doi.org/10.1002/1522-2594(200007)44:1<162::aid-mrm23>3.0.co;2-e
Gordon
,
E. M.
,
Chauvin
,
R. J.
,
Van
,
A. N.
,
Rajesh
,
A.
,
Nielsen
,
A.
,
Newbold
,
D. J.
,
Lynch
,
C. J.
,
Seider
,
N. A.
,
Krimmel
,
S. R.
,
Scheidter
,
K. M.
,
Monk
,
J.
,
Miller
,
R. L.
,
Metoki
,
A.
,
Montez
,
D. F.
,
Zheng
,
A.
,
Elbau
,
I.
,
Madison
,
T.
,
Nishino
,
T.
,
Myers
,
M. J.
, …
Dosenbach
,
N. U. F
. (
2023
).
A somato-cognitive action network alternates with effector regions in motor cortex
.
Nature
,
617
(
7960
),
351
359
. https://doi.org/10.1038/s41586-023-05964-2
Gozzi
,
N.
,
Chee
,
L.
,
Odermatt
,
I.
,
Kikkert
,
S.
,
Preatoni
,
G.
,
Valle
,
G.
,
Pfender
,
N.
,
Beuschlein
,
F.
,
Wenderoth
,
N.
,
Zipser
,
C.
, &
Raspopovic
,
S.
(
2024
).
Wearable neuroprosthesis improves mobility and reduces pain in neuropathic participants
.
medRxiv
, 2024.2005.2008.24306164. https://doi.org/10.1101/2024.05.08.24306164
Greve
,
D. N.
, &
Fischl
,
B.
(
2009
).
Accurate and robust brain image alignment using boundary-based registration
.
Neuroimage
,
48
(
1
),
63
72
. https://doi.org/10.1016/j.neuroimage.2009.06.060
Griffanti
,
L.
,
Douaud
,
G.
,
Bijsterbosch
,
J.
,
Evangelisti
,
S.
,
Alfaro-Almagro
,
F.
,
Glasser
,
M. F.
,
Duff
,
E. P.
,
Fitzgibbon
,
S.
,
Westphal
,
R.
,
Carone
,
D.
,
Beckmann
,
C. F.
, &
Smith
,
S. M.
(
2017
).
Hand classification of fMRI ICA noise components
.
Neuroimage
,
154
,
188
205
. https://doi.org/10.1016/j.neuroimage.2016.12.036
Halder
,
P.
,
Kambi
,
N.
,
Chand
,
P.
, &
Jain
,
N.
(
2018
).
Altered expression of reorganized inputs as they ascend from the cuneate nucleus to cortical area 3b in monkeys with long-term spinal cord injuries
.
Cereb Cortex
,
28
(
11
),
3922
3938
. https://doi.org/10.1093/cercor/bhx256
Harrison
,
O. K.
,
Brooks
,
J. C. W.
, &
Clare
,
S.
(
2023
).
Chapter 18—Brainstem imaging
. In
K.
Markenroth Bloch
,
M.
Guye
, &
B. A.
Poser
(Eds.),
Advances in magnetic resonance technology and applications
(Vol.
10
, pp.
287
299
).
Academic Press
. https://doi.org/10.1016/B978-0-323-99898-7.00017-1
Harvey
,
A. K.
,
Pattinson
,
K. T. S.
,
Brooks
,
J. C. W.
,
Mayhew
,
S. D.
,
Jenkinson
,
M.
, &
Wise
,
R. G.
(
2008
).
Brainstem functional magnetic resonance imaging: Disentangling signal from physiological noise
.
J Magn Reson Imaging
,
28
(
6
),
1337
1344
. https://doi.org/10.1002/jmri.21623
Haufe
,
S.
,
Meinecke
,
F.
,
Görgen
,
K.
,
Dähne
,
S.
,
Haynes
,
J.-D.
,
Blankertz
,
B.
, &
Bießmann
,
F.
(
2014
).
On the interpretation of weight vectors of linear models in multivariate neuroimaging
.
Neuroimage
,
87
,
96
110
. https://doi.org/10.1016/j.neuroimage.2013.10.067
Haxby
,
J. V.
(
2012
).
Multivariate pattern analysis of fMRI: The early beginnings
.
Neuroimage
,
62
(
2
),
852
855
. https://doi.org/10.1016/j.neuroimage.2012.03.016
Haxby
,
J. V.
,
Gobbini
,
M. I.
,
Furey
,
M. L.
,
Ishai
,
A.
,
Schouten
,
J. L.
, &
Pietrini
,
P.
(
2001
).
Distributed and overlapping representations of faces and objects in ventral temporal cortex
.
Science
,
293
(
5539
),
2425
2430
. https://doi.org/10.1126/science.1063736
Haynes
,
J.-D.
, &
Rees
,
G.
(
2005
).
Predicting the orientation of invisible stimuli from activity in human primary visual cortex
.
Nature neuroscience
,
8
(
5
),
686
691
. https://doi.org/10.1038/nn1445
Haynes
,
J.-D.
, &
Rees
,
G.
(
2006
).
Decoding mental states from brain activity in humans
.
Nat Rev Neurosci
,
7
(
7
),
523
534
. https://doi.org/10.1038/nrn1931
Hebart
,
M. N.
, &
Baker
,
C. I.
(
2018
).
Deconstructing multivariate decoding for the study of brain function
.
Neuroimage
,
180
,
4
18
. https://doi.org/10.1016/j.neuroimage.2017.08.005
Hebart
,
M. N.
,
Görgen
,
K.
, &
Haynes
,
J. D
.
(
2014
).
The Decoding Toolbox (TDT): A versatile software package for multivariate analyses of functional imaging data
.
Front Neuroinform
,
8
,
88
. https://doi.org/10.3389/fninf.2014.00088
Jain
,
N.
,
Florence
,
S. L.
,
Qi
,
H. X.
, &
Kaas
,
J. H.
(
2000
).
Growth of new brainstem connections in adult monkeys with massive sensory loss
.
Proc Natl Acad Sci U S A
,
97
(
10
),
5546
5550
. https://doi.org/10.1073/pnas.090572597
Jenkinson
,
M.
,
Bannister
,
P.
,
Brady
,
M.
, &
Smith
,
S.
(
2002
).
Improved optimization for the robust and accurate linear registration and motion correction of brain images
.
Neuroimage
,
17
(
2
),
825
841
. https://doi.org/10.1016/s1053-8119(02)91132-8
Jenkinson
,
M.
, &
Smith
,
S.
(
2001
).
A global optimisation method for robust affine registration of brain images
.
Med Image Analysis
,
5
(
2
),
143
156
. https://doi.org/10.1016/S1361-8415(01)00036-6
Jimura
,
K.
, &
Poldrack
,
R. A.
(
2012
).
Analyses of regional-average activation and multivoxel pattern information tell complementary stories
.
Neuropsychologia
,
50
(
4
),
544
552
. https://doi.org/10.1016/j.neuropsychologia.2011.11.007
Jutzeler
,
C. R
.
,
Curt
,
A
.
, &
Kramer
,
J. L
.
(
2015
).
Relationship between chronic pain and brain reorganization after deafferentation: A systematic review of functional MRI findings
.
NeuroImage Clin
,
9
,
599
606
. https://doi.org/10.1016/j.nicl.2015.09.018
Kaas
,
J. H.
,
Florence
,
S. L.
, &
Jain
,
N.
(
1999
).
Subcortical contributions to massive cortical reorganizations
.
Neuron
,
22
(
4
),
657
660
. https://doi.org/10.1016/S0896-6273(00)80725-4
Kambi
,
N.
,
Halder
,
P.
,
Rajan
,
R.
,
Arora
,
V.
,
Chand
,
P.
,
Arora
,
M.
, &
Jain
,
N.
(
2014
).
Large-scale reorganization of the somatosensory cortex following spinal cord injuries is due to brainstem plasticity
.
Nat Commun
,
5
(
1
),
3602
. https://doi.org/10.1038/ncomms4602
Kamitani
,
Y.
, &
Tong
,
F.
(
2005
).
Decoding the visual and subjective contents of the human brain
.
Nat Neurosci
,
8
(
5
),
679
685
. https://doi.org/10.1038/nn1444
Kasper
,
L.
,
Bollmann
,
S.
,
Diaconescu
,
A. O.
,
Hutton
,
C.
,
Heinzle
,
J.
,
Iglesias
,
S.
,
Hauser
,
T. U.
,
Sebold
,
M.
,
Manjaly
,
Z. M.
,
Pruessmann
,
K. P.
, &
Stephan
,
K. E.
(
2017
).
The PhysIO toolbox for modeling physiological noise in fMRI data
.
J Neurosci Methods
,
276
,
56
72
. https://doi.org/10.1016/j.jneumeth.2016.10.019
Keine
,
C
.
,
Rübsamen
,
R.
, &
Englitz
,
B
. (
2016
).
Inhibition in the auditory brainstem enhances signal representation and regulates gain in complex acoustic environments
.
eLife
,
5
,
e19295
. https://doi.org/10.7554/eLife.19295
Kikkert
,
S.
,
Kolasinski
,
J.
,
Jbabdi
,
S.
,
Tracey
,
I.
,
Beckmann
,
C. F.
,
Johansen-Berg
,
H.
, &
Makin
,
T. R.
(
2016
).
Revealing the neural fingerprints of a missing hand
.
eLife
,
5
,
e15292
. https://doi.org/10.7554/eLife.15292
Kikkert
,
S.
,
Pfyffer
,
D.
,
Verling
,
M.
,
Freund
,
P.
, &
Wenderoth
,
N.
(
2021
).
Finger somatotopy is preserved after tetraplegia but deteriorates over time
.
eLife
,
10
,
e67713
. https://doi.org/10.7554/elife.67713
Kikkert
,
S.
,
Sonar
,
H. A.
,
Freund
,
P.
,
Paik
,
J.
, &
Wenderoth
,
N.
(
2023
).
Hand and face somatotopy shown using MRI-safe vibrotactile stimulation with a novel soft pneumatic actuator (SPA)-skin interface
.
Neuroimage
,
269
,
119932
. https://doi.org/10.1016/j.neuroimage.2023.119932
Knights
,
E.
,
Smith
,
F. W.
, &
Rossit
,
S.
(
2022
).
The role of the anterior temporal cortex in action: Evidence from fMRI multivariate searchlight analysis during real object grasping
.
Sci Rep
,
12
(
1
),
9042
. https://doi.org/10.1038/s41598-022-12174-9
Kolasinski
,
J.
,
Makin
,
T. R.
,
Jbabdi
,
S.
,
Clare
,
S.
,
Stagg
,
C. J.
, &
Johansen-Berg
,
H.
(
2016
).
Investigating the stability of fine-grain digit somatotopy in individual human participants
.
J Neurosci
,
36
(
4
),
1113
1127
. https://doi.org/10.1523/jneurosci.1742-15.2016
Komisaruk
,
B. R.
,
Mosier
,
K. M.
,
Liu
,
W.-C.
,
Criminale
,
C.
,
Zaborszky
,
L.
,
Whipple
,
B.
, &
Kalnin
,
A.
(
2002
).
Functional localization of brainstem and cervical spinal cord nuclei in humans with fMRI
.
Am J Neuroradiol
,
23
(
4
),
609
617
. https://www.ajnr.org/content/ajnr/23/4/609.full.pdf
Krauth
,
A.
,
Blanc
,
R.
,
Poveda
,
A.
,
Jeanmonod
,
D.
,
Morel
,
A.
, &
Székely
,
G.
(
2010
).
A mean three-dimensional atlas of the human thalamus: Generation from multiple histological data
.
Neuroimage
,
49
(
3
),
2053
2062
. https://doi.org/10.1016/j.neuroimage.2009.10.042
Kriegeskorte
,
N.
,
Goebel
,
R.
, &
Bandettini
,
P.
(
2006
).
Information-based functional brain mapping
.
Proc Natl Acad Sci USA
,
103
(
10
),
3863
3868
. https://doi.org/10.1073/pnas.0600244103
Kubota
,
S.
,
Sasaki
,
C.
,
Kikuta
,
S.
,
Yoshida
,
J.
,
Ito
,
S.
,
Gomi
,
H.
,
Oya
,
T.
, &
Seki
,
K.
(
2024
).
Modulation of somatosensory signal transmission in the primate cuneate nucleus during voluntary hand movement
.
Cell Reports
, 43(3), 113884. https://doi.org/10.1016/j.celrep.2024.113884
Kuehn
,
E.
,
Haggard
,
P.
,
Villringer
,
A.
,
Pleger
,
B.
, &
Sereno
,
M. I.
(
2018
).
Visually-driven maps in area 3b
.
J Neurosci
,
38
(
5
),
1295
1310
. https://doi.org/10.1523/jneurosci.0491-17.2017
Kuehn
,
E.
,
Mueller
,
K.
,
Turner
,
R.
, &
Schütz-Bosbach
,
S.
(
2014
).
The functional architecture of S1 during touch observation described with 7 T fMRI
.
Brain Struct Funct
,
219
(
1
),
119
140
. https://doi.org/10.1007/s00429-012-0489-z
Kündig
,
C. W.
,
Finsterbusch
,
J.
,
Freund
,
P.
, &
David
,
G.
(
2024
).
Functional magnetic resonance imaging of the lumbosacral cord during a lower extremity motor task
.
Imaging Neurosci
,
2
,
1
19
. https://doi.org/10.1162/imag_a_00227
Liao
,
C. C.
,
Qi
,
H. X.
,
Reed
,
J. L.
,
Jeoung
,
H. S.
, &
Kaas
,
J. H.
(
2021
).
Corticocuneate projections are altered after spinal cord dorsal column lesions in New World monkeys
.
J Comp Neurol
,
529
(
7
),
1669
1702
. https://doi.org/10.1002/cne.25050
Limbrick-Oldfield
,
E. H.
,
Brooks
,
J. C.
,
Wise
,
R. J.
,
Padormo
,
F.
,
Hajnal
,
J. V.
,
Beckmann
,
C. F.
, &
Ungless
,
M. A.
(
2012
).
Identification and characterisation of midbrain nuclei using optimised functional magnetic resonance imaging
.
Neuroimage
,
59
(
2
),
1230
1238
. https://doi.org/10.1016/j.neuroimage.2011.08.016
Marfurt
,
C. F.
, &
Rajchert
,
D. M.
(
1991
).
Trigeminal primary afferent projections to “non-trigeminal” areas of the rat central nervous system
.
J Comp Neurol
,
303
(
3
),
489
511
. https://doi.org/10.1002/cne.903030313
Martuzzi
,
R.
,
van der Zwaag
,
W.
,
Farthouat
,
J.
,
Gruetter
,
R.
, &
Blanke
,
O.
(
2014
).
Human finger somatotopy in areas 3b, 1, and 2: A 7T fMRI study using a natural stimulus
.
Hum Brain Mapp
,
35
(
1
),
213
226
. https://doi.org/10.1002/hbm.22172
Matt
,
E.
,
Fischmeister
,
F. P. S.
,
Amini
,
A.
,
Robinson
,
S. D.
,
Weber
,
A.
,
Foki
,
T.
,
Gizewski
,
E. R.
, &
Beisteiner
,
R.
(
2019
).
Improving sensitivity, specificity, and reproducibility of individual brainstem activation
.
Brain Struct Funct
,
224
(
8
),
2823
2838
. https://doi.org/10.1007/s00429-019-01936-3
Meissner
,
S. N.
,
Bächinger
,
M.
,
Kikkert
,
S.
,
Imhof
,
J.
,
Missura
,
S.
,
Carro Dominguez
,
M.
, &
Wenderoth
,
N.
(
2024
).
Self-regulating arousal via pupil-based biofeedback
.
Nat Hum Behav
,
8
(
1
),
43
62
. https://doi.org/10.1038/s41562-023-01729-z
Menon
,
R. S.
, &
Goodyear
,
B. G.
(
2001
).
Spatial and temporal resolution in fMRI
. In
P.
Jezzard
,
P. M.
Matthews
, &
S. M.
Smith
(Eds.),
Functional magnetic resonance imaging: An introduction to methods
(pp.
146
158
).
Oxford University Press
. https://doi.org/10.1093/acprof:oso/9780192630711.003.0007
Mur
,
M.
,
Bandettini
,
P. A.
, &
Kriegeskorte
,
N.
(
2009
).
Revealing representational content with pattern-information fMRI—An introductory guide
.
Soc Cogn Affect Neurosci
,
4
(
1
),
101
109
. https://doi.org/10.1093/scan/nsn044
Muret
,
D.
, &
Makin
,
T. R.
(
2021
).
The homeostatic homunculus: Rethinking deprivation-triggered reorganisation
.
Curr Opin Neurobiol
,
67
,
115
122
. https://doi.org/10.1016/j.conb.2020.08.008
Muret
,
D.
,
Root
,
V.
,
Kieliba
,
P.
,
Clode
,
D.
, &
Makin
,
T. R.
(
2022
).
Beyond body maps: Information content of specific body parts is distributed across the somatosensory homunculus
.
Cell Reports
,
38
(
11
),
110523
. https://doi.org/10.1016/j.celrep.2022.110523
Nambu
,
I.
,
Hagura
,
N.
,
Hirose
,
S.
,
Wada
,
Y.
,
Kawato
,
M.
, &
Naito
,
E.
(
2015
).
Decoding sequential finger movements from preparatory activity in higher-order motor regions: A functional magnetic resonance imaging multi-voxel pattern analysis
.
Eur J Neurosci
,
42
(
10
),
2851
2859
. https://doi.org/10.1111/ejn.13063
Norman
,
K. A.
,
Polyn
,
S. M.
,
Detre
,
G. J.
, &
Haxby
,
J. V.
(
2006
).
Beyond mind-reading: Multi-voxel pattern analysis of fMRI data
.
Trends Cogn Sci
,
10
(
9
),
424
430
. https://doi.org/10.1016/j.tics.2006.07.005
Odermatt
,
I. A.
,
Schulthess-Lutz
,
M.
,
Mihelj
,
E.
,
Howell
,
P.
,
Heimhofer
,
C.
,
McMackin
,
R.
,
Ruddy
,
K.
,
Freund
,
P.
,
Kikkert
,
S.
, &
Wenderoth
,
N.
(
2024
).
TMS-based neurofeedback training of mental finger individuation induces neuroplastic changes in the sensorimotor cortex
.
bioRxiv
, 2024.2005.2016.594100. https://doi.org/10.1101/2024.05.16.594100
Oldfield
,
R. C.
(
1971
).
The assessment and analysis of handedness: The Edinburgh inventory
.
Neuropsychologia
,
9
(
1
),
97
113
. https://doi.org/10.1016/0028-3932(71)90067-4
Oliva
,
V.
,
Gregory
,
R.
,
Davies
,
W.-E.
,
Harrison
,
L.
,
Moran
,
R.
,
Pickering
,
A. E.
, &
Brooks
,
J. C. W.
(
2021
).
Parallel cortical-brainstem pathways to attentional analgesia
.
Neuroimage
, 226, 117548. https://doi.org/10.1016/j.neuroimage.2020.117548
Oliva
,
V.
,
Hartley-Davies
,
R.
,
Moran
,
R.
,
Pickering
,
A. E.
, &
Brooks
,
J. C. W.
(
2022
).
Simultaneous brain, brainstem, and spinal cord pharmacological-fMRI reveals involvement of an endogenous opioid network in attentional analgesia
.
eLife
,
11
,
e71877
. https://doi.org/10.7554/eLife.71877
Pajula
,
J.
, &
Tohka
,
J.
(
2016
).
How many is enough? Effect of sample size in inter-subject correlation analysis of fMRI
.
Comput Intell Neurosci
,
2016
,
2094601
. https://doi.org/10.1155/2016/2094601
Pattinson
,
K. T.
,
Mitsis
,
G. D
.
,
Harvey
,
A. K.
,
Jbabdi
,
S.
,
Dirckx
,
S.
,
Mayhew
,
S. D.
,
Rogers
,
R.
,
Tracey
,
I.
, &
Wise
,
R. G
.
(
2009
).
Determination of the human brainstem respiratory control network and its cortical connections in vivo using functional and structural imaging
.
NeuroImage
,
44
(
2
),
295
305
. https://doi.org/10.1016/j.neuroimage.2008.09.007
Paxinos
,
G.
, &
Mai
,
J. K.
(
2004
).
The human nervous system
. Elsevier. https://doi.org/10.1016/b978-012547626-3/50001-6
Penfield
,
W.
, &
Boldrey
,
E.
(
1937
).
Somatic motor and sensory representation in the cerebral cortex of man as studied by electrical stimulation
.
Brain
,
60
(
4
),
389
443
. https://doi.org/10.1093/brain/60.4.389
Penfield
,
W.
, &
Rasmussen
,
T.
(
1950
).
The cerebral cortex of man: A clinical study of localization of function
.
Journal of the American Medical Association
,
144
(
16
),
1412
1412
. https://doi.org/10.1001/jama.1950.02920160086033
Pereira
,
F.
,
Mitchell
,
T.
, &
Botvinick
,
M
.
(
2009
).
Machine learning classifiers and fMRI: A tutorial overview
.
NeuroImage
,
45
(
Suppl 1
),
S199
S209
. https://doi.org/10.1016/j.neuroimage.2008.11.007
Puckett
,
A. M.
,
Bollmann
,
S.
,
Barth
,
M.
, &
Cunnington
,
R.
(
2017
).
Measuring the effects of attention to individual fingertips in somatosensory cortex using ultra-high field (7T) fMRI
.
Neuroimage
,
161
,
179
187
. https://doi.org/10.1016/j.neuroimage.2017.08.014
Purves
,
D
.
(
2019
).
Neuroscience
(6th international ed.). Oxford University Press.
Rabe
,
F.
,
Kikkert
,
S.
, &
Wenderoth
,
N.
(
2023
).
Performing a vibrotactile discrimination task modulates finger representations in primary somatosensory cortex
.
J Neurophysiol
,
130
(
4
),
1015
1027
. https://doi.org/10.1152/jn.00428.2022
Reddy
,
N. A.
,
Clements
,
R. G.
,
Brooks
,
J. C. W.
, &
Bright
,
M. G.
(
2024
).
Simultaneous cortical, subcortical, and brainstem mapping of sensory activation
.
Cerebr Cortex
,
34
(
6
). https://doi.org/10.1093/cercor/bhae273
Roux
,
F. E.
,
Djidjeli
,
I.
, &
Durand
,
J. B.
(
2018
).
Functional architecture of the somatosensory homunculus detected by electrostimulation
.
J Physiol
,
596
(
5
),
941
956
. https://doi.org/10.1113/jp275243
Samuels
,
E. R.
, &
Szabadi
,
E.
(
2008
).
Functional neuroanatomy of the noradrenergic locus coeruleus: Its roles in the regulation of arousal and autonomic function part II: Physiological and pharmacological manipulations and pathological alterations of locus coeruleus activity in humans
.
Curr Neuropharmacol
,
6
(
3
),
254
285
. https://doi.org/10.2174/157015908785777193
Sanchez-Panchuelo
,
R. M.
,
Besle
,
J.
,
Beckett
,
A.
,
Bowtell
,
R.
,
Schluppeck
,
D.
, &
Francis
,
S.
(
2012
).
Within-digit functional parcellation of Brodmann areas of the human primary somatosensory cortex using functional magnetic resonance imaging at 7 tesla
.
J Neurosci
,
32
(
45
),
15815
15822
. https://doi.org/10.1523/jneurosci.2501-12.2012
Sanchez-Panchuelo
,
R. M.
,
Francis
,
S.
,
Bowtell
,
R.
, &
Schluppeck
,
D.
(
2010
).
Mapping human somatosensory cortex in individual subjects with 7T functional MRI
.
J Neurophysiol
,
103
(
5
),
2544
2556
. https://doi.org/10.1152/jn.01017.2009
Sanders
,
Z.-B.
,
Dempsey-Jones
,
H.
,
Wesselink
,
D. B.
,
Edmondson
,
L. R.
,
Puckett
,
A. M.
,
Saal
,
H. P.
, &
Makin
,
T. R.
(
2023
).
Similar somatotopy for active and passive digit representation in primary somatosensory cortex
.
Hum Brain Mapp
,
44
(
9
),
3568
3585
. https://doi.org/10.1002/hbm.26298
Satpute
,
A. B.
,
Wager
,
T. D.
,
Cohen-Adad
,
J.
,
Bianciardi
,
M.
,
Choi
,
J. K.
,
Buhle
,
J. T.
,
Wald
,
L. L.
, &
Barrett
,
L. F.
(
2013
).
Identification of discrete functional subregions of the human periaqueductal gray
.
Proc Natl Acad Sci U S A
,
110
(
42
),
17101
17106
. https://doi.org/10.1073/pnas.1306095110
Schellekens
,
W.
,
Petridou
,
N.
, &
Ramsey
,
N. F.
(
2018
).
Detailed somatotopy in primary motor and somatosensory cortex revealed by Gaussian population receptive fields
.
Neuroimage
,
179
,
337
347
. https://doi.org/10.1016/j.neuroimage.2018.06.062
Sclocco
,
R.
,
Beissner
,
F.
,
Bianciardi
,
M.
,
Polimeni
,
J. R.
, &
Napadow
,
V.
(
2018
).
Challenges and opportunities for brainstem neuroimaging with ultrahigh field MRI
.
Neuroimage
,
168
,
412
426
. https://doi.org/10.1016/j.neuroimage.2017.02.052
Seidler
,
R. D.
,
Bernard
,
J. A.
,
Burutolu
,
T. B.
,
Fling
,
B. W.
,
Gordon
,
M. T.
,
Gwin
,
J. T.
,
Kwak
,
Y.
, &
Lipps
,
D. B.
(
2010
).
Motor control and aging: Links to age-related brain structural, functional, and biochemical effects
.
Neurosci Biobehav Rev
,
34
(
5
),
721
733
. https://doi.org/10.1016/j.neubiorev.2009.10.005
Smith
,
S. A.
,
Dortch
,
R. D.
,
Barry
,
R. L.
, &
Gore
,
J. C.
(
2014
).
Chapter 2.4—Ultra-high field spinal cord imaging
. In
J.
Cohen-Adad
&
C. A. M.
Wheeler-Kingshott
(Eds.),
Quantitative MRI of the spinal cord
(pp.
106
119
).
Academic Press
. https://doi.org/10.1016/B978-0-12-396973-6.00008-3
Smith
,
S. M.
(
2002
).
Fast robust automated brain extraction
.
Hum Brain Mapp
,
17
(
3
),
143
155
. https://doi.org/10.1002/hbm.10062
Smith
,
S. M.
,
Jenkinson
,
M.
,
Woolrich
,
M. W.
,
Beckmann
,
C. F.
,
Behrens
,
T. E.
,
Johansen-Berg
,
H.
,
Bannister
,
P. R.
,
De Luca
,
M.
,
Drobnjak
,
I.
,
Flitney
,
D. E.
,
Niazy
,
R. K.
,
Saunders
,
J.
,
Vickers
,
J.
,
Zhang
,
Y.
,
De Stefano
,
N.
,
Brady
,
J. M.
, &
Matthews
,
P. M.
(
2004
).
Advances in functional and structural MR image analysis and implementat
.
Neuroimage
,
23
(
Suppl. 1
),
S208
S219
. https://doi.org/10.1016/j.neuroimage.2004.07.051
Smith
,
S. M.
, &
Nichols
,
T. E.
(
2009
).
Threshold-free cluster enhancement: Addressing problems of smoothing, threshold dependence and localisation in cluster inference
.
Neuroimage
,
44
(
1
),
83
98
. https://doi.org/10.1016/j.neuroimage.2008.03.061
Stelzer
,
J.
,
Chen
,
Y.
, &
Turner
,
R
.
(
2013
).
Statistical inference and multiple testing correction in classification-based multi-voxel pattern analysis (MVPA): Random permutations and cluster size control
.
NeuroImage
,
65
,
69
82
. https://doi.org/10.1016/j.neuroimage.2012.09.063
Taylor
,
P. A.
,
Reynolds
,
R. C.
,
Calhoun
,
V.
,
Gonzalez-Castillo
,
J.
,
Handwerker
,
D. A.
,
Bandettini
,
P. A.
,
Mejia
,
A. F.
, &
Chen
,
G.
(
2023
).
Highlight results, don’t hide them: Enhance interpretation, reduce biases and improve reproducibility
.
NeuroImage
,
274
,
120138
. https://doi.org/10.1016/j.neuroimage.2023.120138
Versteeg
,
C.
,
Chowdhury
,
R. H.
, &
Miller
,
L. E.
(
2021
).
Cuneate nucleus: The somatosensory gateway to the brain
.
Curr Opin Physiol
,
20
,
206
215
. https://doi.org/10.1016/j.cophys.2021.02.004
Versteeg
,
C.
,
Rosenow
,
J. M.
,
Bensmaia
,
S. J.
, &
Miller
,
L. E.
(
2021
).
Encoding of limb state by single neurons in the cuneate nucleus of awake monkeys
.
J Neurophysiol
,
126
(
2
),
693
706
. https://doi.org/10.1152/jn.00568.2020
Wall
,
M. B.
,
Walker
,
R.
, &
Smith
,
A. T.
(
2009
).
Functional imaging of the human superior colliculus: An optimised approach
.
Neuroimage
,
47
(
4
),
1620
1627
. https://doi.org/10.1016/j.neuroimage.2009.05.094
Wang
,
Q.
,
Cagna
,
B.
,
Chaminade
,
T.
, &
Takerkart
,
S.
(
2020
).
Inter-subject pattern analysis: A straightforward and powerful scheme for group-level MVPA
.
Neuroimage
, 204, 116205. https://doi.org/10.1016/j.neuroimage.2019.116205
Wesselink
,
D. B.
,
Sanders
,
Z. B.
,
Edmondson
,
L. R.
,
Dempsey-Jones
,
H.
,
Kieliba
,
P.
,
Kikkert
,
S.
,
Themistocleous
,
A. C.
,
Emir
,
U.
,
Diedrichsen
,
J.
,
Saal
,
H. P.
, &
Makin
,
T. R.
(
2022
).
Malleability of the cortical hand map following a finger nerve block
.
Sci Adv
,
8
(
16
),
eabk2393
. https://doi.org/10.1126/sciadv.abk2393
Wesselink
,
D. B.
,
van den Heiligenberg
,
F. M. Z.
,
Ejaz
,
N.
,
Dempsey-Jones
,
H.
,
Cardinali
,
L.
,
Tarall-Jozwiak
,
A.
,
Diedrichsen
,
J.
, &
Makin
,
T. R.
(
2019
).
Obtaining and maintaining cortical hand representation as evidenced from acquired and congenital handlessness
.
eLife
,
8
,
e37227
. https://doi.org/10.7554/eLife.37227
Winkler
,
A. M.
,
Ridgway
,
G. R.
,
Webster
,
M. A.
,
Smith
,
S. M.
, &
Nichols
,
T. E.
(
2014
).
Permutation inference for the general linear model
.
Neuroimage
, 92,
381
397
. https://doi.org/10.1016/j.neuroimage.2014.01.060
Woo
,
C. W.
,
Krishnan
,
A.
, &
Wager
,
T. D.
(
2014
).
Cluster-extent based thresholding in fMRI analyses: Pitfalls and recommendations
.
NeuroImage
,
91
,
412
419
. https://doi.org/10.1016/j.neuroimage.2013.12.058
Woolrich
,
M. W.
,
Behrens
,
T. E.
, &
Smith
,
S. M.
(
2004
).
Constrained linear basis sets for HRF modelling using Variational Bayes
.
Neuroimage
,
21
(
4
),
1748
1761
. https://doi.org/10.1016/j.neuroimage.2003.12.024
Woolrich
,
M. W.
,
Behrens
,
T. E. J.
,
Beckmann
,
C. F.
,
Jenkinson
,
M.
, &
Smith
,
S. M.
(
2004
).
Multilevel linear modelling for FMRI group analysis using Bayesian inference
.
Neuroimage
,
21
(
4
),
1732
1747
. https://doi.org/10.1016/j.neuroimage.2003.12.023
Woolrich
,
M. W.
,
Ripley
,
B. D.
,
Brady
,
M.
, &
Smith
,
S. M.
(
2001
).
Temporal autocorrelation in univariate linear modeling of FMRI data
.
Neuroimage
,
14
(
6
),
1370
1386
. https://doi.org/10.1006/nimg.2001.0931
Worsley
,
K. J.
, &
Friston
,
K. J
.
(
1995
).
Analysis of fMRI time-series revisited—again
.
NeuroImage
,
2
(
3
),
173
181
. https://doi.org/10.1006/nimg.1995.1023
Youssef
,
A. M.
,
Macefield
,
V. G.
, &
Henderson
,
L. A.
(
2016
).
Pain inhibits pain; human brainstem mechanisms
.
Neuroimage
,
124
(
Pt A
),
54
62
. https://doi.org/10.1016/j.neuroimage.2015.08.060
Zaaimi
,
B.
,
Edgley
,
S. A.
,
Soteropoulos
,
D. S.
, &
Baker
,
S. N.
(
2012
).
Changes in descending motor pathway connectivity after corticospinal tract lesion in macaque monkey
.
Brain
,
135
(
7
),
2277
2289
. https://doi.org/10.1093/brain/aws115
Zörner
,
B.
,
Bachmann
,
L. C.
,
Filli
,
L.
,
Kapitza
,
S.
,
Gullo
,
M.
,
Bolliger
,
M.
,
Starkey
,
M. L.
,
Röthlisberger
,
M.
,
Gonzenbach
,
R. R.
, &
Schwab
,
M. E.
(
2014
).
Chasing central nervous system plasticity: The brainstem’s contribution to locomotor recovery in rats with spinal cord injury
.
Brain
,
137
(
6
),
1716
1732
. https://doi.org/10.1093/brain/awu078
This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International (CC BY 4.0) license, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. For a full description of the license, please visit https://creativecommons.org/licenses/by/4.0/legalcode.

Supplementary data