
==== Front
J Neurophysiol
J Neurophysiol
JN
Journal of Neurophysiology
0022-3077
1522-1598
American Physiological Society Rockville, MD

38505916
JN-00407-2023
JN-00407-2023
10.1152/jn.00407.2023
Research Article
neuroimaging-meets-neurophysiologyNeuroimaging Meets NeurophysiologyActivity in the pontine reticular nuclei scales with handgrip force in humans
PONTINE RETICULAR NUCLEI SCALES WITH HANDGRIP FORCE
Danielson Tyler L. 1
Gould Layla A. 2
https://orcid.org/0000-0002-5108-8079
DeFreitas Jason M. 3
https://orcid.org/0000-0002-3144-8866
MacLennan Rob J. 4 5
Ekstrand Chelsea 6
https://orcid.org/0000-0003-0002-4021
Borowsky Ron 7
Farthing Jonathan P. 8
https://orcid.org/0000-0003-2258-1689
Andrushko Justin W. 9
1Applied Neuromuscular Physiology Laboratory, College of Education and Human Sciences, Oklahoma State University , Stillwater, Oklahoma, United States
2Division of Neurosurgery, Department of Surgery, College of Medicine, University of Saskatchewan , Saskatoon, Saskatchewan, Canada
3Department of Exercise Science, Falk College of Sport and Human Dynamics, Syracuse University , Syracuse, New York, United States
4Department of Neurology, College of Medicine, University of Florida , Gainesville, Florida, United States
5Brain Rehabilitation Research Center, Malcom Randall VA Medical Center , Gainesville, Florida, United States
6Department of Neuroscience, Faculty of Arts and Science, University of Lethbridge , Lethbridge, Alberta, Canada
7Department of Psychology and Health Studies, College of Arts and Science, University of Saskatchewan , Saskatoon, Saskatchewan, Canada
8College of Kinesiology, University of Saskatchewan , Saskatoon, Saskatchewan, Canada
9Department of Sport, Exercise and Rehabilitation, Faculty of Health and Life Sciences, https://ror.org/049e6bc10 Northumbria University , Newcastle upon Tyne, United Kingdom
Correspondence: J. W. Andrushko (Justin.andrushko@northumbria.ac.uk); J. P. Farthing (jon.farthing@usask.ca).
1 5 2024
20 3 2024
20 3 2024
131 5 807814
2 11 2023
21 2 2024
14 3 2024
Copyright © 2024 The Authors.
2024
The Authors
https://creativecommons.org/licenses/by/4.0/ Licensed under Creative Commons Attribution CC-BY 4.0. Published by the American Physiological Society.

The neural pathways that contribute to force production in humans are currently poorly understood, as the relative roles of the corticospinal tract and brainstem pathways, such as the reticulospinal tract (RST), vary substantially across species. Using functional magnetic resonance imaging (fMRI), we aimed to measure activation in the pontine reticular nuclei (PRN) during different submaximal handgrip contractions to determine the potential role of the PRN in force modulation. Thirteen neurologically intact participants (age: 28 ± 6 yr) performed unilateral handgrip contractions at 25%, 50%, 75% of maximum voluntary contraction during brain scans. We quantified the magnitude of PRN activation from the contralateral and ipsilateral sides during each of the three contraction intensities. A repeated-measures ANOVA demonstrated a significant main effect of force (P = 0.012, ηp2 = 0.307) for PRN activation, independent of side (i.e., activation increased with force for both contralateral and ipsilateral nuclei). Further analyses of these data involved calculating the linear slope between the magnitude of activation and handgrip force for each region of interest (ROI) at the individual-level. One-sample t tests on the slopes revealed significant group-level scaling for the PRN bilaterally, but only the ipsilateral PRN remained significant after correcting for multiple comparisons. We show evidence of task-dependent activation in the PRN that was positively related to handgrip force. These data build on a growing body of literature that highlights the RST as a functionally relevant motor pathway for force modulation in humans.

NEW & NOTEWORTHY In this study, we used a task-based functional magnetic resonance imaging (fMRI) paradigm to show that activity in the pontine reticular nuclei scales linearly with increasing force during a handgrip task. These findings directly support recently proposed hypotheses that the reticulospinal tract may play an important role in modulating force production in humans.

Listen to this article’s corresponding podcast at https://jneurophysiol.podbean.com/e/jnp-micro-podcast-pontine-reticular-nuclei-scales-with-handgrip-force/.

brain imaging
; descending motor pathways
; force modulation
motor control
; reticular formation
; Canadian Government | Natural Sciences and Engineering Research Council of Canada (NSERC) 10.13039/501100000038 2016-0529 Jonathan P. FarthingCanadian Government | Natural Sciences and Engineering Research Council of Canada (NSERC) 10.13039/501100000038 183968-2013-25 Ron BorowskyCall for PapersTrue
==== Body
pmcINTRODUCTION

Within the central nervous system, there are multiple descending pathways that control motor function, and the roles and relative contributions from each pathway may vary substantially across species. In nonmammals, such as reptiles and birds, the main motor pathways descending onto motoneurons originate in the brainstem (1). However, in humans, it is commonly thought that upper limb fine motor control and force modulation is primarily controlled through the corticospinal tract [CST (2, 3)]. An added difficulty in quantifying the role of the CST is how much its size and functionality differs between species. Within humans, the cortico-motoneuronal connections are substantially more developed than in nonhuman primates and other mammals (4), suggesting that the role of CST in humans likely cannot be inferred accurately from animal-model research. Increased development of the CST correlates positively with indices of dexterity (5), giving further evidence for the role of the CST in fine, distal motor production. Due to this specialization of the CST within humans, most research studies investigating the mechanisms of human motor control and muscle force have primarily focused on the CST, with equivocal and heterogeneous findings (6, 7).

Recent evidence, however, now suggests the reticulospinal tract (RST), a brainstem pathway thought to be limited in humans to postural and proximal muscle control, may be more involved in upper limb motor control than previously thought (8–10). In support of this, there has been evidence of distal RST innervation in nonhuman primate models (11) and in humans (8, 9). In addition to its contribution to gross movement control, recent papers have suggested that the RST may also have a role in force modulation (12–14). Some studies have used animal models to study the RST’s role in voluntary force production, but these models require the implantation of electrodes on the RST that is not possible in humans (10, 11, 15–17). For example, Glover and Baker (16) found that the RST activity appears to contribute to contraction force in macaque monkeys. Although, like our study, the evidence was only correlational. Others have examined this control in humans with stroke and spinal cord injury, allowing for models analogous to animal lesion studies (18, 19).

In humans, these studies primarily involve noninvasive stimulation to produce evoked contractions (19–21), whereas animal models have allowed for studies during voluntary movement (11, 15–17, 22). Examining the role of the RST during voluntary force production is important to elucidate how it contributes to the control of human movement. To our knowledge, no study to date has used functional magnetic resonance imaging (fMRI) to corroborate the role of the RST for controlling force in humans despite the benefits of its noninvasive nature. The purpose of this study was to quantify whether the pontine reticular nuclei (PRN), which give rise to the RST, are more active during handgrip contractions at higher compared with lower intensities.

MATERIALS AND METHODS

As a secondary analysis on a data set that has previously been published (23), we aimed to determine the relationship between fMRI blood oxygenation level dependent (BOLD) signal modulation in the pontine reticular formation nuclei and parametric increases in handgrip force during three separate fMRI runs (Fig. 1). Medullary reticular formation nuclei were not included in the analysis due to the medulla not being consistently in the field of view across participants during fMRI scans.

Figure 1. A: handgrip contractions at three different force levels were performed during three separate functional magnetic resonance imaging (fMRI) runs of the cerebrum and brainstem. The monitor in the top right shows our primary region of interest (ROI) in MNI-space, the pontine nuclei of the reticulospinal tract (RST). B: a simplified diagram of the hypothesized pathways involved in voluntary handgrip force production. For comparison, activity levels in the contralateral motor cortex (M1) were also assessed. Not shown, but also potential contributors, are the contralateral supplemental motor area (SMA) and its bilateral projections to the RST on each side, potential connections to the RST from the motor and premotor cortices, and crossed and uncrossed pathways from the ipsilateral M1. Figure was created in biorender.com.

Participants

Thirteen neurologically intact individuals volunteered to participate in this study (11 right-handed, 2 left-handed; sex: 7 females, 6 males; age: 28 ± 6 yr; height: 170.9 ± 9.8 cm; weight: 75.1 ± 16.7 kg). Participants were eligible if they were 18 yr or older and were MRI safe. The presence of intracranial metal clips, cardiac pacemaker, or any factor that would prevent a safe MRI scan disqualified individuals to participate.

Ethical Approval

All participants provided written informed consent before participation. The study adhered to the Declaration of Helsinki and received approval from the University of Saskatchewan Biomedical Research Ethics Board (Ethics No. Bio 01-125).

Experimental Outline

During a single fMRI session, participants engaged in three experimental tasks. These tasks were executed in a random sequence and consisted of submaximal unimanual right isometric handgrip contractions at 25%, 50%, and 75% of their maximum voluntary contraction (MVC). The session began with the acquisition of a structural MRI brain scan, followed by having the participants perform three MVCs. To prevent fatigue, a 1-min rest interval was allotted between each MVC attempt. Based on the MVC results, target force levels were set for the three submaximal experimental tasks. Participants then proceeded to complete the three experimental tasks during separate fMRI runs in random order.

Behavioral Motor Task

Each of the submaximal handgrip conditions were acquired in separate block design fMRI runs. For each run, participants performed 5 sets × 5 repetitions of either a 25%, 50%, or 75% MVC handgrip visuomotor target matching task with an MRI-compatible hand clench dynamometer (Biopac Systems Inc. Aero Camino Goleta, CA). Task blocks involved 1,650 ms [i.e., corresponding to the repetition time (TR) of the T2* fMRI scan] contractions alternating with 1,650 ms of rest (16,500 ms total task block). Each task block was separated by blocks of complete rest (16,500 ms total rest block). During each fMRI run, the participants wore MRI compatible goggles and viewed a computer screen projection, where they were able to view a custom-built LabView (v.8.6) interface. Participants saw clear target lines and go/no-go on screen flashing lights and were cued when to contract or relax. The LabView interface was triggered by the MRI through a fiber optic cable that sent a signal at the beginning of each TR. This signal went through a Siemens fMRI trigger converter and passed to the LabView interface via a BNC connection. This allowed for precise continuous synchronization between the MRI and the experimental paradigm computer throughout the scan (24). Target lines were presented relative to the individual’s peak MVC and force feedback was presented as a vertical force bar that was responsive to each participant’s handgrip contraction (i.e., harder contraction resulted in the bar rising vertically). Participants were visually instructed when to contract and when to relax with a green light that would turn on when the participant was to contract. A second red light was programmed to turn on during rest blocks to indicate a sustained rest. Participants were instructed to relax their nonactive arm throughout each of the fMRI runs.

Handgrip Motor Task Processing

The handgrip force data were processed in Matlab using custom scripts (25). Data were processed with a fourth-order 100 Hz low-pass Butterworth filter and full-wave rectified. The mean from within the full width half maximum of each contraction was then normalized to the mean kg-force of the highest MVC and expressed as a percentage of MVC. The mean of each task block (five contractions) was calculated and used in subsequent analyses.

fMRI Parameters

All scans were acquired in a Siemens 3 T MAGNETOM Skyra MRI (Siemens Healthcare, Erlangen, Germany). At the start of the session, a whole brain T1-weighted (T1w) structural scan was acquired using a high-resolution magnetization prepared rapid acquisition gradient echo (MPRAGE) sequence consisting of 192 T1w echo-planar imaging slices (1-mm slice thickness with no gap), with an in-plane resolution of 1 × 1 mm [field of view = 256 × 256; TR = 1,900 ms; echo time (TE) = 2.08 ms]. Each of the experimental fMRI runs was acquired using T2*-weighted single-shot gradient-echo planar imaging with an interleaved ascending sequence, consisting of 105 volumes (TR = 1,650 ms; TE = 30 ms) of 25 axial slices of 4 mm thickness (1-mm gap) with an in-plane resolution of 2.7 mm × 2.7 mm (field of view = 250) using a flip angle of 90°. The top two coil sets (16 channels) of a 20-channel Siemens head coil (Siemens Healthcare) were used. Runs consisted of a 10-volume alternating block design beginning with five volumes for stabilization (task, rest; 105 volumes total).

fMRI Preprocessing

Functional MRI data processing was carried out using FMRI Expert Analysis Tool (FEAT) Version 6.00, as part of FSL [FMRIB’s Software Library (26)]. Boundary-based registration was used to register the functional image to the high-resolution T1w image. Registration of the functional images to the T1w image was carried out using FMRIB’s Linear Image Registration Tool [FLIRT (27, 28)], and then nonlinearly (rigid + affine + deformable syn) registered to the MNI152_1mm brain template using antsRegistrationSyN and antsApplyTransforms from Advanced Normalization Tools [ANTs (29, 30)]. The following preprocessing was applied: motion correction using Motion Correction FMRIB’s Linear Image Registration Tool [MCFLIRT (28)]; non-brain removal using FNIRT-based brain extraction as implemented in the fsl_anat anatomical processing script (31); spatial smoothing using a Gaussian kernel of FWHM 6 mm; grand-mean intensity normalization of the entire four-dimension (4-D) data set by a single multiplicative factor (32). Next, Independent Component Analysis Automatic Removal of Motion Artifacts (ICA-AROMA) was used remove motion-related noise from each functional run (32). Following ICA-AROMA data clean up, data were high-pass temporal filtered with a 0.01-Hz cutoff. Time-series statistical analyses were carried out using FMRIB’s Improved Linear Model (FILM) with local autocorrelation correction (33). Z (Gaussianised T/F) statistic images were constructed nonparametrically using Gaussian Random Field theory-based maximum height thresholding with a corrected significance threshold of P = 0.05 (34).

Region of Interest Signal Change

FSL’s featquery tool was used to extract the BOLD signal percent change from three regions of interest (ROI) in participants’ native-space. For each ROI, ANTs were used to transform ROIs from MNI-space to native-space prior to signal change extraction. The percent signal change was extracted from the left contralateral and right ipsilateral pontine reticular formation nuclei (left: PnO_PnC_l; right: PnO_PnC_r) from the Brainstem navigator toolkit [https://www.nitrc.org/projects/brainstemnavig (35–41)]. The Navigator Toolkit atlas was developed in vivo by acquiring multiple contrasts (T1w, T2w, T2*, diffusion-weighted imaging) at an ultra-high field strength (7 T) with high spatial resolution (1.1 mm isotropic), and distortion matched diffusion and T2w echo planar images. This approach facilitated the accurate delineation of white and gray matter boundaries on a subject-by-subject level after data were registered to a 1-mm stereotaxic space (MNI standard space). Anatomically defined subnuclei were able to be identified in stereotaxic space based on their microstructural diffusion properties when assessed in conjunction with the other image contrasts. This approach avoided misregistration confounds across the anatomical images. Neighboring nuclei and landmarks were used to initially identify structures, and labels were averaged across subjects, thereby making a probabilistic atlas. Finally, this was validated by inter-rater agreement, internal consistency, and nuclei volume (35–41). In addition, the percent signal change was extracted from left contralateral motor cortex (M1) using an ROI mask from the Brainnetome atlas [left M1: A4ul_l (42)]. The extracted values were then used in subsequent statistical analyses.

Statistical Analysis

An omnibus ROI [contralateral pontine reticular nuclei (PRN), ipsilateral PRN] × condition (25%, 50%, 75% MVC) repeated-measures analysis of variance (RM-ANOVA) was run using the magnitude of brain activation from each of the ROI’s (i.e., BOLD signal percent change from rest) to determine whether PRN activation changed with handgrip force modulation. Data were assessed for normality using Shapiro–Wilks in addition to Zskewness and Zkurtosis. ANOVA tests are robust against violations to normality (43, 44) and thus minor violations to normality as determined by Zskewness and Zkurtosis were considered acceptable for running the parametric model. Sphericity was assessed with Mauchley’s test of sphericity and where violated the Greenhouse-Geisser correction was used. Hedges g effect sizes were calculated and reported to adjust for the small sample size. A Bonferroni adjustment was used to correct for multiple comparisons where appropriate post hoc testing was conducted. The α level for significance was set to α = 0.05.

In addition, for each ROI the slope coefficient of brain activation across the three submaximal handgrip force conditions (percent change of BOLD signal activation/percent target force) was calculated on a subject-by-subject basis using linear regression. The slope coefficients were then analyzed using separate one-sample t tests to determine if each ROI scales its magnitude of activation across the varying submaximal force outputs. A positive slope indicates activation in that ROI increases with handgrip force. A negative slope would represent a decrease in activation with an increase in handgrip force and a nonsignificant slope (i.e., not different from zero) would represent that activation did not vary as a function of handgrip force. For this analysis, the α level was set to α = 0.017 (0.05/3 ROIs; contralateral PRN, ipsilateral PRN, contralateral M1) and P values ≤0.017 were considered significant.

RESULTS

Magnitude of BOLD Activation across Contraction Intensities

RM-ANOVA was performed to test the effect of force (25%, 50%, and 75% MVC) and side (ipsilateral and contralateral) on PRN activation. Mauchly’s test of sphericity was violated for Force (P = 0.023) and Force × Side (P = 0.006). There was a significant main effect for force [F(1.337,16.039) = 5.306, P = 0.012, ηp2 = 0.307], but no main effect for side [F(1,12) = 0.413, P = 0.533, ηp2 = 0.033]. The interaction effect (Force × Side) was also not significant [F(1.243,14.921) = 0.957, P = 0.364, ηp2 = 0.074]. A Bonferroni post hoc analysis showed that activation of pontine nuclei at 75% MVC was significantly greater than at 25% MVC [mean difference (MD) = −0.2, SE = 0.068, P = 0.02] as well as at 50% MVC (MD = −0.179, SE = 0.068, P = 0.042). The difference in activation between 25% and 50% MVC was not significant (MD = −0.021, SE = 0.068, P = 1.000; Fig. 2A).

Figure 2. A: the activation in the pontine nuclei of the reticular formation during contractions of varying intensities. The contralateral and ipsilateral nuclei are collapsed together because there was neither a side × force intensity interaction nor a main effect for side as a factor. Pontine activation did show a main effect for force intensity, with post hoc pairwise comparisons shown above the box plots. B: the relationship between % activation and % force from a single region of interest (ROI) from a participant with a linear fit to clearly demonstrate the slope analysis. The relationship was fit with a linear model, and the slope interpreted as the degree to which activation for that region changes with increases in force. C: the slope coefficients from each subject for each ROI. The values above the scatter plot demonstrate the percent of subjects that exhibited a positive slope (i.e., positive scaling in activation with increasing force) for each ROI. The P values above the box plots are from one-sample t tests, showing that the mean of the slopes was significantly different than zero in the pons before correcting for multiple comparisons, but not in motor cortex (M1). Only the ipsilateral pontine reticular nuclei (PRN) significantly scaled after statistical correction (corrected α = 0.017). The thick gray horizontal lines in A and C extend the zero point for easy comparison.

ROI BOLD Activation Scaling with Force

Linear regression analyses were conducted on each ROI for each subject to determine the slope coefficient, which quantifies the rate of increase in activation with increasing contraction intensity (Fig. 2B). Next, using separate one sample t tests of the individual slopes, we showed at the group-level pontine reticular activation was significantly different from zero on both the ipsilateral [t(12) = 2.672, P = 0.01, g = 0.712] and contralateral [t(12) = 1.792, P = 0.049, g = 0.478] sides before correction; however, after a manual Bonferroni correction for multiple comparisons, only the ipsilateral side remained significant. The contralateral M1 was not statistically different from zero [t(12) = 0.245, P = 0.405, g = 0.065], indicating it did not reliably scale with force (Fig. 2C).

DISCUSSION

The present study examined the relationship between unimanual handgrip force production and BOLD signal activation of the contralateral and ipsilateral PRN of the RST, as well as the contralateral M1 for comparison. Our primary finding is that activation within the PRN consistently scaled linearly with increased force production. The traditional understanding is that the RST primarily innervates axial and proximal muscle groups in humans (45). However, over the past 15 yr a small body of literature has begun to challenge this belief (8, 9, 11, 46) with research suggesting the RST may reach distal musculature of the arms and hands and modulate force (12–14). With this task, we have shown, for the first time, PRN BOLD signal activation significantly scales with handgrip force, where a positive relationship between the magnitude of activation and force output are observed. Specifically, we observe significant scaling in the ipsilateral PRN (with a large effect size) after correcting for multiple comparisons. There was also scaling with a moderate effect size in the contralateral PRN. In addition, the RM-ANOVA showed that side (ipsilateral vs. contralateral) had no effect on activation across the force conditions. These findings, paired with the lack of scaling in contralateral M1 provide some support for the notion of the PRN and possibly the RST pathway contribute to high force outputs, rather than the traditional belief that human distal motor control is exclusively governed by contralateral cortical regions (47). Firing rates within populations of M1 cells are highly variable with force modulation (6, 16, 47), which is supported by a lack of M1 BOLD signal scaling with varying handgrip contraction forces (23). Glover and Baker (16) also found that cells in the reticular formation almost exclusively had a positive correlation between firing rate and force. Our data in humans support their conclusion, as we found no robust support for the role of M1 in force modulation in our previous publication (23), and add to it with evidence that force modulation may be driven by the PRN, and possibly the RST.

During isometric contractions, the recruitment of motor units tend to occur in an orderly fashion from low-threshold to higher-threshold (48). As force demand increases, larger and higher-threshold motor units are recruited (49–52), and the firing rates of smaller, lower-threshold motor units that were already active increase. These recruitment properties may also affect neuromuscular plasticity, as the repeated activation of the larger motor units with high-intensity resistance training often leads to greater gains in strength compared with lower-intensity training (53). In addition to the impact of high-intensity contractions on motoneuron recruitment, and strength, a recent review by Hortobágyi et al. (54) reported high-intensity resistance training also produces greater motor-related neuroplasticity in the neurologically intact young adult brain. However, mechanistic studies into the specific neural pathways and adaptations to resistance training that affords this motoneuronal recruitment remains incompletely understood. This incomplete mechanistic understanding is due in part to research focusing on the CST, which has resulted in equivocal insight. Here, we add to the growing body of literature that pinpoints the PRN, and as a potentially relevant nuclei for force control with a handgrip task, which, until recently (12–14), was thought to be exclusively controlled by the CST (55). Specifically, we show evidence for intensity-dependent activation of the PRN that is correlated with force output, which suggests this is a task relevant pathway for handgrip force. Furthermore, based on the principles of use-dependent plasticity, the repeated preferential activation of the PRN with high-intensity contractions over the course of a resistance training intervention may lead to the strengthening of the RST pathway’s integrity in a manner that facilitates strength adaptations. Recent findings have led to proposals that the RST may have a role in modulating motor unit discharge rates, which would in turn alter force and rate of force production (14). Škarabot et al. (14) proposed hypothesis about the role of RST in modulating force that is supported by the findings of this study.

Evaluating the RST in applied settings presents its own set of challenges, such as the need for external stimuli paired with auditory perturbations (e.g., transcranial magnetic stimulation or cervical medullary electrical stimulation paired with loud acoustic startle sounds ≥ 110 dB), in addition to a high variance in responses. Here, using fMRI, we show a noninvasive evaluation of the PRN in a controlled environment during voluntary motor behaviors. We argue this is an important contribution to the field investigating the role of the RST, as these data provide support to previous findings without the need for external perturbations. Furthermore, when paired with previous research (12, 13, 15), our findings indicate that the RST may be a functionally relevant pathway contributing to muscle force and, by extension, possibly chronic strength training adaptations.

Limitations and Future Directions

Investigating patterns of brain activation during high-intensity contractions with fMRI is not a trivial task, as head motion commonly occurs during effortful movements, and can result in the loss of data due to motion-related contamination. With that, the investigation of fMRI-related BOLD signal activation during different types of motor tasks may not be possible, and thus these experiments are limited to simpler isolated hand movement patterns that are less likely to produce high amounts of motion.

Though we present evidence for the role of the PRN in gross distal hand movement, future research may consider testing PRN activation and RST excitability during finer distal movements. By focusing on dexterous movements such as flexion, extension, abduction, adduction, and circumduction, RST activation bias toward certain types of movement may be determined. Future research into this area should investigate the contributions of the RST to different tasks and their adaptations to training interventions with varying manipulations. A deeper understanding of how the RST contributes to motor function could lead to better rehabilitative modalities for those with dysfunction of the CST, such as individuals with stroke (8, 12). In addition, we focused only on the pontine nuclei of the RST. The medullary nuclei had to be excluded because the medulla was cut off from the scan for many of the participants in favor of capturing complete cortex. Future work should examine whether the medullary nuclei also play a potential role in force modulation. Finally, future research would benefit from replication of this experiment but with a larger sample size to build confidence in data interpretation and robustness of findings. Two final limitations or considerations of this work, that is inherent with all fMRI research, that should be noted when interpreting the data in relation to single-cell invasive measures of neurophysiology are 1) fMRI is not directly measuring neural activity, but it is modeling the hemodynamic response which is related to the local field potentials at the cellular level (56), and 2) the spatial resolution of fMRI is not capable of measuring single-cell neural activity. Within a typical fMRI voxel, there can be a few million individual neurons, and even more synapses that make up the measured signal (57). These two considerations make it difficult to directly compare these population-level hemodynamic measures with findings from single-cell recordings in invasive studies.

Conclusions

In a small cohort of neurologically intact young adults, we demonstrate that activation of the PRN of the RST scales in magnitude with handgrip force. These data provide preliminary insight into the potential role of the PRN and possibly the RST during distal muscular force modulation, a concept most often considered under exclusive control by the CST in humans. These data add to the growing body of literature suggesting a potentially important role of the RST for neuromuscular strength adaptations in humans and lays the foundation for which future research can further investigate RST function to disentangle the neural underpinnings of force modulation and the neural adaptations governing increases in strength.

DATA AVAILABILITY

Data are available upon reasonable request.

GRANTS

This study was funded by a Natural Sciences and Engineering Research Council of Canada (NSERC) Discovery Grant 2016-0529 (to Jonathan P. Farthing), in addition to an NSERC Discovery Grant: 183968-2013-25 (to Ron Borowsky).

DISCLOSURES

No conflicts of interest, financial or otherwise, are declared by the authors.

AUTHOR CONTRIBUTIONS

L.A.G., R.B., J.P.F., and J.W.A. conceived and designed research; L.A.G., C.E., R.B., J.P.F., and J.W.A. performed experiments; T.L.D., C.E., and J.W.A. analyzed data; J.M.D., R.B., J.P.F., and J.W.A. interpreted results of experiments; J.M.D. and J.W.A. prepared figures; T.L.D. and J.W.A. drafted manuscript; J.M.D., R.J.M., R.B., J.P.F., and J.W.A. edited and revised manuscript; T.L.D., L.A.G., J.M.D., R.J.M., C.E., R.B., J.P.F., and J.W.A. approved final version of manuscript.

ACKNOWLEDGMENTS

The authors thank Doug Renshaw, Sul Ross State University, for assistance with data collection, and Shawn Reinink, University of Saskatchewan, for contribution to the LabView software interface used in this experiment.
==== Refs
REFERENCES

1. Kalman M. Evolution of the brain: at the reptile-bird transition. In: Encyclopedia of Neuroscience, edited by Binder MD, Hirokawa N, Windhorst U. Berlin, Heidelberg: Springer Berlin Heidelberg, 2009, p. 1305–1312.
2. Evarts EV. Activity of pyramidal tract neurons during postural fixation. J Neurophysiol 32 : 375–385, 1969. doi:10.1152/JN.1969.32.3.375. 4977837
3. Evarts EV. Relation of pyramidal tract activity to force exerted during voluntary movement. J Neurophysiol 31 : 14–27, 1968. doi:10.1152/JN.1968.31.1.14. 4966614
4. Lemon RN. Descending pathways in motor control. Annu Rev Neurosci 31 : 195–218, 2008. doi:10.1146/annurev.neuro.31.060407.125547. 18558853
5. Heffner R, Masterton B. Variation in form of the pyramidal tract and its relationship to digital dexterity. Brain Behav Evol 12 : 161–200, 1975. doi:10.1159/000124401. 1212616
6. Maier MA, Bennett KMB, Hepp-Reymond MC, Lemon RN. Contribution of the monkey corticomotoneuronal system to the control of force in precision grip. J Neurophysiol 69 : 772–785, 1993. doi:10.1152/jn.1993.69.3.772. 8463818
7. Kidgell DJ, Bonanno DR, Frazer AK, Howatson G, Pearce AJ. Corticospinal responses following strength training: a systematic review and meta-analysis. Eur J Neurosci 46 : 2648–2661, 2017. doi:10.1111/EJN.13710. 28921683
8. Baker SN, Perez MA. Reticulospinal contributions to gross hand function after human spinal cord injury. J Neurosci 37 : 9778–9784, 2017. doi:10.1523/JNEUROSCI.3368-16.2017. 28871033
9. Honeycutt CF, Kharouta M, Perreault EJ. Evidence for reticulospinal contributions to coordinated finger movements in humans. J Neurophysiol 110 : 1476–1483, 2013. doi:10.1152/jn.00866.2012. 23825395
10. Maitland S, Baker SN. Ipsilateral motor evoked potentials as a measure of the reticulospinal tract in age-related strength changes. Front Aging Neurosci 13 : 612352, 2021. doi:10.3389/FNAGI.2021.612352. 33746734
11. Riddle CN, Edgley SA, Baker SN. Direct and indirect connections with upper limb motoneurons from the primate reticulospinal tract. J Neurosci 29 : 4993–4999, 2009. doi:10.1523/JNEUROSCI.3720-08.2009. 19369568
12. Akalu Y, Frazer AK, Howatson G, Pearce AJ, Siddique U, Rostami M, Tallent J, Kidgell DJ. Identifying the role of the reticulospinal tract for strength and motor recovery: a scoping review of nonhuman and human studies. Physiol Rep 11 : e15765, 2023. doi:10.14814/PHY2.15765. 37474275
13. Atkinson E, Škarabot J, Ansdell P, Goodall S, Howatson G, Thomas K. Does the reticulospinal tract mediate adaptation to resistance training in humans? J Appl Physiol (1985) 133 : 689–696, 2022. doi:10.1152/japplphysiol.00264.2021. 35834623
14. Škarabot J, Folland JP, Holobar A, Baker SN, Del Vecchio A. Startling stimuli increase maximal motor unit discharge rate and rate of force development in humans. J Neurophysiol 128 : 455–469, 2022. doi:10.1152/jn.00115.2022. 35829632
15. Glover IS, Baker SN. Cortical, corticospinal, and reticulospinal contributions to strength training. J Neurosci 40 : 5820–5832, 2020. doi:10.1523/JNEUROSCI.1923-19.2020. 32601242
16. Glover IS, Baker SN. Both corticospinal and reticulospinal tracts control force of contraction. J Neurosci 42 : 3150–3164, 2022. doi:10.1523/JNEUROSCI.0627-21.2022. 35241490
17. Tapia JA, Tohyama T, Poll A, Baker SN. The existence of the StartReact effect implies reticulospinal, not corticospinal, inputs dominate drive to motoneurons during voluntary movement. J Neurosci 42 : 7634–7647, 2022. doi:10.1523/JNEUROSCI.2473-21.2022. 36658461
18. Hammerbeck U, Tyson SF, Samraj P, Hollands K, Krakauer JW, Rothwell J. The strength of the corticospinal tract not the reticulospinal tract determines upper-limb impairment level and capacity for skill-acquisition in the sub-acute post-stroke period. Neurorehabil Neural Repair 35 : 812–822, 2021. doi:10.1177/15459683211028243. 34219510
19. Xu J, Ejaz N, Hertler B, Branscheidt M, Widmer M, Faria AV, Harran MD, Cortes JC, Kim N, Celnik PA, Kitago T, Luft AR, Krakauer JW, Diedrichsen J. Separable systems for recovery of finger strength and control after stroke. J Neurophysiol 118 : 1151–1163, 2017. doi:10.1152/JN.00123.2017. 28566461
20. Alagona G, Delvaux V, Gérard P, De Pasqua V, Pennisi G, Delwaide PJ, Nicoletti F, Maertens De Noordhout A. Ipsilateral motor responses to focal transcranial magnetic stimulation in healthy subjects and acute-stroke patients. Stroke 32 : 1304–1309, 2001. doi:10.1161/01.STR.32.6.1304. 11387491
21. Ziemann U, Ishii K, Borgheresi A, Yaseen Z, Battaglia F, Hallett M, Cincotta M, Wassermann EM. Dissociation of the pathways mediating ipsilateral and contralateral motor-evoked potentials in human hand and arm muscles. J Physiol 518 : 895–906, 1999. doi:10.1111/J.1469-7793.1999.0895P.X. 10420023
22. Baker SN. The primate reticulospinal tract, hand function and functional recovery. J Physiol 589 : 5603–5612, 2011. doi:10.1113/jphysiol.2011.215160. 21878519
23. Andrushko JW, Gould LA, Renshaw DW, Ekstrand C, Hortobágyi T, Borowsky R, Farthing JP. High force unimanual handgrip contractions increase ipsilateral sensorimotor activation and functional connectivity. Neuroscience 452 : 111–125, 2021. doi:10.1016/j.neuroscience.2020.10.031. 33197497
24. Ekstrand C, Neudorf J, Kress S, Borowsky R. Structural connectivity predicts functional activation during lexical and sublexical reading. Neuroimage 218 : 117008, 2020. doi:10.1016/J.NEUROIMAGE.2020.117008. 32485306
25. The MathWorks Inc. R2018a. Natick, MA: The MathWorks Inc.
26. Jenkinson M, Beckmann CF, Behrens TEJJ, Woolrich MW, Smith SM. FSL. Neuroimage 62 : 782–790, 2012. doi:10.1016/j.neuroimage.2011.09.015. 21979382
27. Jenkinson M, Smith S. A global optimisation method for robust affine registration of brain images. Med Image Anal 5 : 143–156, 2001. doi:10.1016/S1361-8415(01)00036-6. 11516708
28. Jenkinson M, Bannister P, Brady M, Smith S. Improved optimization for the robust and accurate linear registration and motion correction of brain images. Neuroimage 17 : 825–841, 2002. doi:10.1016/S1053-8119(02)91132-8. 12377157
29. Avants BB, Epstein CL, Grossman M, Gee JC. Symmetric diffeomorphic image registration with cross-correlation: evaluating automated labeling of elderly and neurodegenerative brain. Med Image Anal 12 : 26–41, 2008. doi:10.1016/J.MEDIA.2007.06.004. 17659998
30. Avants B, Tustison NJ, Song G. Advanced normalization tools: V1.0. Insight J 2009. doi:10.54294/uvnhin.
31. Smith SM, Jenkinson M, Woolrich MW, Beckmann CF, Behrens TEJ, Johansen-Berg H, Bannister PR, De Luca M, Drobnjak I, Flitney DE, Niazy RK, Saunders J, Vickers J, Zhang Y, De Stefano N, Brady JM, Matthews PM. Advances in functional and structural MR image analysis and implementation as FSL. NeuroImage 23 , Suppl 1 : S208–S219, 2004. doi:10.1016/j.neuroimage.2004.07.051. 15501092
32. Pruim RHR, Mennes M, van Rooij D, Llera A, Buitelaar JK, Beckmann CF. ICA-AROMA: a robust ICA-based strategy for removing motion artifacts from fMRI data. Neuroimage 112 : 267–277, 2015. doi:10.1016/j.neuroimage.2015.02.064. 25770991
33. Woolrich MW, Ripley BD, Brady M, Smith SM. Temporal autocorrelation in univariate linear modeling of FMRI data. Neuroimage 14 : 1370–1386, 2001. doi:10.1006/nimg.2001.0931. 11707093
34. Worsley K. Statistical analysis of activation images. In: Functional MRI: An Introduction to Methods, edited by Jezzard P, Matthews PM, Smith SM. Oxford, UK: Oxford University Press, 2001, p. 251–270.
35. Bianciardi M, Toschi N, Edlow BL, Eichner C, Setsompop K, Polimeni JR, Brown EN, Kinney HC, Rosen BR, Wald LL. Toward an in vivo neuroimaging template of human brainstem nuclei of the ascending arousal, autonomic, and motor systems. Brain Connect 5 : 597–607, 2015. doi:10.1089/brain.2015.0347. 26066023
36. García-Gomar MG, Strong C, Toschi N, Singh K, Rosen BR, Wald LL, Bianciardi M. In vivo probabilistic structural atlas of the inferior and superior colliculi, medial and lateral geniculate nuclei and superior olivary complex in humans based on 7 tesla MRI. Front Neurosci 13 : 764, 2019. doi:10.3389/fnins.2019.00764. 31440122
37. Singh K, Indovina I, Augustinack JC, Nestor K, García-Gomar MG, Staab JP, Bianciardi M. Probabilistic template of the lateral parabrachial nucleus, medial parabrachial nucleus, vestibular nuclei complex, and medullary viscero-sensory-motor nuclei complex in living humans from 7 tesla MRI. Front Neurosci 13 : 1425, 2019. doi:10.3389/fnins.2019.01425. 32038134
38. Singh K, García-Gomar MG, Bianciardi M. Probabilistic atlas of the mesencephalic reticular formation, isthmic reticular formation, microcellular tegmental nucleus, ventral tegmental area nucleus complex, and caudal-rostral linear raphe nucleus complex in living humans from 7 tesla magnetic resonance imaging. Brain Connect 11 : 613–623, 2021. doi:10.1089/brain.2020.0975. 33926237
39. Bianciardi M, Strong C, Toschi N, Edlow BL, Fischl B, Brown EN, Rosen BR, Wald LL. A probabilistic template of human mesopontine tegmental nuclei from in vivo 7 T MRI. Neuroimage 170 : 222–230, 2018. doi:10.1016/j.neuroimage.2017.04.070. 28476663
40. García-Gomar MG, Videnovic A, Singh K, Stauder M, Lewis LD, Wald LL, Rosen BR, Bianciardi M. Disruption of brainstem structural connectivity in REM sleep behavior disorder using 7 tesla magnetic resonance imaging. Mov Disord 37 : 847–853, 2022. doi:10.1002/mds.28895. 34964520
41. Bianciardi M. Brainstem Navigator. NITRC, 2021. https://www.nitrc.org/projects/brainstemnavig/.
42. Fan L, Li H, Zhuo J, Zhang Y, Wang J, Chen L, Yang Z, Chu C, Xie S, Laird AR, Fox PT, Eickhoff SB, Yu C, Jiang T. The human brainnetome atlas: a new brain atlas based on connectional architecture. Cereb Cortex 26 : 3508–3526, 2016. doi:10.1093/cercor/bhw157. 27230218
43. Schmider E, Ziegler M, Danay E, Beyer L, Bühner M. Is it really robust?: reinvestigating the robustness of ANOVA against violations of the normal distribution assumption. Methodology 6 : 147–151, 2010. doi:10.1027/1614-2241/a000016.
44. Blanca MJ, Alarcón R, Arnau J, Bono R, Bendayan R. Non-normal data: is ANOVA still a valid option? Psicothema 29 : 552–557, 2017. doi:10.7334/psicothema2016.383. 29048317
45. Kuypers HGJM. Anatomy of the descending pathways. In: Comprehensive Physiology, edited by Terjung R. Wiley, 2011, p. 597–666. doi:10.1002/cphy.cp010213.
46. Zaaimi B, Edgley SA, Soteropoulos DS, Baker SN. Changes in descending motor pathway connectivity after corticospinal tract lesion in macaque monkey. Brain 135 : 2277–2289, 2012. doi:10.1093/BRAIN/AWS115. 22581799
47. Schmidt EM, Jost RG, Davis KK. Reexamination of the force relationship of cortical cell discharge patterns with conditioned wrist movements. Brain Res 83 : 213–223, 1975. doi:10.1016/0006-8993(75)90931-2. 1109294
48. De Luca CJ, Hostage EC. Relationship between firing rate and recruitment threshold of motoneurons in voluntary isometric contractions. J Neurophysiol 104 : 1034–1046, 2010 [E rratum in J Neurophysiol 107: 1544, 2012]. doi:10.1152/JN.01018.2009. 20554838
49. Milner-Brown HS, Stein RB, Yemm R. The orderly recruitment of human motor units during voluntary isometric contractions. J Physiol 230 : 359–370, 1973. doi:10.1113/jphysiol.1973.sp010192. 4350770
50. Milner-Brown HS, Stein RB, Yemm R. Changes in firing rate of human motor units during linearly changing voluntary contractions. J Physiol 230 : 371–390, 1973. doi:10.1113/JPHYSIOL.1973.SP010193. 4708898
51. de Luca CJ, LeFever RS, McCue MP, Xenakis AP. Behaviour of human motor units in different muscles during linearly varying contractions. J Physiol 329 : 113–128, 1982. doi:10.1113/JPHYSIOL.1982.SP014293. 7143246
52. Milner-Brown HS, Stein RB, Yemm R. The contractile properties of human motor units during voluntary isometric contractions. J Physiol 228 : 285–306, 1973. doi:10.1113/JPHYSIOL.1973.SP010087. 4687100
53. Fry AC. The role of resistance exercise intensity on muscle fibre adaptations. Sports Med 34 : 663–679, 2004. doi:10.2165/00007256-200434100-00004. 15335243
54. Hortobágyi T, Vetrovsky T, Balbim GM, Sorte Silva NCB, Manca A, Deriu F, Kolmos M, Kruuse C, Liu-Ambrose T, Radák Z, Váczi M, Johansson H, dos Santos PCR, Franzén E, Granacher U. The impact of aerobic and resistance training intensity on markers of neuroplasticity in health and disease. Ageing Res Rev 80 : 101698, 2022. doi:10.1016/J.ARR.2022.101698. 35853549
55. Natali AL, Reddy V, Bordoni B. Neuroanatomy, Corticospinal Cord Tract (Online). Treasure Island, FL: StatPearls Publishing, 2023.
56. Logothetis NK, Pauls J, Augath M, Trinath T, Oeltermann A. Neurophysiological investigation of the basis of the fMRI signal. Nature 412 : 150–157, 2001. doi:10.1038/35084005. 11449264
57. Huettel SA, Song AW, McCarthy G. Functional Magnetic Resonance Imaging (2nd ed.). Sunderland, MA: Sinauer Associates Inc., 2004.
