Abstract
Background:
Upper extremity impairment is a common disabling consequence of multiple sclerosis (MS). Little is known about factors influencing upper extremity functioning in MS.
Objective:
Identify associations with longitudinal changes in upper extremity performance.
Methods:
Using a real-world, retrospective cohort study (4140 persons with MS and 14,544 observations), multivariable repeated-measures analyses were completed for nine-hole peg test (9HPT) dominant and non-dominant hand times and their difference (asymmetry), with clinical, demographic, and social determinants as covariates. Interactions with time (disease duration) were examined.
Results:
Cross-sectionally, dominant, and non-dominant hand times were strongly correlated (ρ = 0.78), and moderately correlated with asymmetry, walking speeds, and depressive symptoms. Longitudinally, older age, males, Hispanics, Blacks, ever smokers, public insurance beneficiaries, greater lower-limb impairment, progressive disease, and greater depressive symptoms were independently associated with worsening dominant and non-dominant times. Public insurance beneficiaries, greater lower-limb impairment, progressive disease, and greater depressive symptoms were associated with greater asymmetry. Males accrued dominant hand dysfunction more rapidly and had greater asymmetry than females over time, while those with increasing lower-limb impairment more rapidly experienced asymmetry over time.
Discussion:
These findings identify high-risk groups and inform future research to guide strategies for preserving upper extremity function in multiple sclerosis.
Introduction
Upper extremity impairment is a common and disabling consequence of multiple sclerosis (MS) limiting personal and instrumental activities of daily functioning such as feeding, dressing, housework, and transportation. 1 It also impacts employment and overall quality of life.1–3 In a large North American Research Committee on MS (NARCOMS) cross-sectional study, hand dysfunction emerged early with only 40% of persons with MS (PwMS) reporting no impairment in the first year, and declining to 31%, 27%, and 19% in years 5, 10, and 15, respectively. 4 Moderate or severe impairment was reported in 8% in year 1, then in 12%, 18%, and 22% for years 5, 10, and 15, respectively. 4 Importantly, hand dysfunction may occur even in PwMS with preserved Expanded Disability Status Scale (EDSS). 5 These observations underscore the need to identify correlates of upper extremity dysfunction, considering both disease duration and impairment in other domains.
While little research has examined upper extremity symptoms and deficits in MS, multiple factors influencing other MS-related impairments may have overlapping and distinct effects. Neurological mechanisms link upper and lower extremity function—including shared corticospinal tract involvement and cerebellar or brainstem lesions—which can impair coordination, strength, and motor control.6,7 Thus, factors known to influence lower extremity function may also affect upper extremity function. Social determinants of health (i.e. individual- and neighborhood-level socioeconomic conditions) and modifiable behavioral/lifestyle factors (i.e. smoking and obesity) have also been linked to longitudinal changes in depressive symptoms, cognition, brain atrophy, and disability in PwMS.8–11 These factors exert broad systemic effects on inflammation and possibly on neurodegeneration, as well as access to care and rehabilitation, making them biologically and clinically relevant to upper extremity outcomes in MS.12–15 Emerging neuroimaging evidence reinforces this rationale, showing that both structural damage (e.g. regional atrophy and demyelination) and functional reorganization within sensorimotor networks are linked to upper-limb motor performance in PwMS. 16 Collectively, these findings suggest that upper- and lower-limb function share overlapping neural substrates, and that disruptions in these pathways may manifest asymmetrically, underscoring the importance of investigating clinical and social correlates of hand impairment alongside established measures of disability.
The nine-hole peg test (9HPT) is a widely used, reliable, and sensitive tool for measuring manual dexterity in PwMS, 17 and it is a key component of the MS Functional Composite. 18 9HPT correlates with greater disability, longer disease duration, EDSS, Timed 25-Foot Walk (T25FW), low-contrast letter acuity, Symbol Digit Modalities Test, and depressive symptom severity.19,20 Upper-limb asymmetry, defined as the absolute difference between dominant and non-dominant 9HPT times, has been reported at high prevalence among those with severe disability and correlates with EDSS.21,22 Asymmetry may also have functional implications, as greater reliance on the less impaired limb may foster compensatory mechanisms that reduce efficiency, increase fatigue, and impair bimanual tasks. Thus, monitoring inter-limb differences may offer complementary information to unilateral measures and may provide additional insight into disability progression.
Most existing MS studies of 9HPT are cross-sectional with limited scope (e.g. reporting correlations with sex and disease duration).23–25 Only one longitudinal study that examined dominant and non-dominant 9HPT impairment in ~650 young clinical trial participants, explored associations with time, age (18–25 vs 26–35 years), and sex. 26 Consequently, important knowledge gaps remain regarding the clinical and demographic factors influencing longitudinal 9HPT performance and asymmetry in diverse, real-world MS populations spanning the age spectrum. These relationships should also be evaluated for independent from lower-limb impairment and depressive symptoms, and whether they change over time. Therefore, our objectives were to: (1) identify demographic, clinical, and social determinants of health associated with longitudinal variation in 9HPT times and asymmetry, after accounting for lower-limb impairment, depressive symptoms, and disease course in PwMS and (2) assess whether any of these associations vary over time by testing interactions with disease duration.
Methods
Study population
We conducted a retrospective cohort study using real-world electronic health record (EHR) information from PwMS receiving care at the Cleveland Clinic’s Mellen Center for Multiple Sclerosis Treatment and Research (Cleveland, OH) and the Lou Ruvo Center for Brain Health (Las Vegas, NV)—tertiary MS referral centers with ~20,000 clinical encounters and ~1000 new patients per year. In 2007, the Knowledge Program © (KP) has electronically captured patient-reported and provider-entered outcomes, linked to the EHR. 27 Patients were included if they were US residents (eFigure 1) and had ⩾2 clinical encounters ⩾6 months apart between January 2008 and June 2016 with 9HPT recorded. 10 This study was approved by the institutional review boards of the Cleveland Clinic, Case Western Reserve University, and The MetroHealth System.
Dependent variables of interest
The 9HPT is a provider-administered task that objectively measures upper extremity function with exceptional reliability (>95%) in PwMS. 17 The average of two trials, time to place nine pegs in nine empty holes in a block followed by the removal of pegs, was recorded separately for self-reported dominant and non-dominant hands. The three outcomes were repeated measures of 9HPT times for dominant and non-dominant hands and asymmetry (absolute difference between hand times).
Independent variables of interest
We assessed clinical and demographic covariates, including disease duration (years), age (years), sex, ethnicity (Hispanic/Latinx, Not Hispanic), and race (White, Black, Other). We included area deprivation index (ADI; state-level ranking of census block socioeconomic disadvantage) of baseline residence. 28 Time-varying covariates were smoking status (never, current, former), insurance payer (private, self-pay, Medicare, and Medicaid as proxies for socioeconomic status (SES) in United States where self-pay, Medicare, and Medicaid are indicative of financial hardship, disability, or limited income-based coverage), marital status, assistive device use (none, unilateral, bilateral), T25FW (seconds; eFigure 2(a)), and depressive symptom severity (Patient Health Questionnaire 9; PHQ-9; eFigure 2(b)). Body mass index (BMI) was included if measured within 3 months of a 9HPT measure. EDSS was unavailable as it is not a part of routine care.
Statistical analysis
Descriptive statistics were reported as means (standard deviation (SD)) for continuous variables and percentages for categorical variables. At baseline, pairwise relationships between measures were evaluated using Spearman’s rank correlation, with 95% confidence intervals (CIs) and p-values derived from nonparametric bootstrapping (1000 repetitions).
Longitudinal associations were evaluated using generalized estimating equations with an identity link function and robust sandwich estimator using STATA v17.0 (StataCorp, TX; xtgee function). Given the unequal time intervals between encounters, we examined the quasi-likelihood for independence model criterion (QIC) for autoregressive (AR1) and exchangeable working correlation structures for the baseline models described below, for patients with ⩾2 encounters (a requirement for AR1). QIC was consistently lower for the exchangeable working correlation structure and, therefore, it was selected. Two longitudinal baseline models were generated per outcome: (1) all covariates except BMI (N = 4140 patients, 14,544 observations) and (2) including BMI (N = 2540 patients and 5090 observations). Fixed covariates were modeled as time-invariant, while smoking, insurance payer, marital status, assistive device use, T25FW, PHQ-9, and BMI were modeled as time-varying. Indicator variables retained observations with missing categorical data. Disease duration (years since diagnosis) was the primary time variable; interactions between disease duration and each covariate were tested to assess whether relationships changed over time, and significant interactions were retained in final models. Scatter, contour, and local polynomial smoothing plots with 95% CIs guided interpretations. Statistical significance was defined as two-sided α = 0.05.
Results
Patients averaged 3.5 encounters (SD = 2.8) over an average of 3.5 years (SD = 2.1; maximum 7.6 years). Baseline attributes are presented in Table 1 (N = 4140): mean age was 46.6 years (SD = 11.5), average disease duration 11.1 years (SD = 9.3), 74% were female, 95% non-Hispanic, 87% White, 11% Black, 45% never smokers, and 56% privately insured. Most patients had relapsing MS (77%) and did not require assistance ambulating (82%). Consistent with global trends, 90.4% were right-handed. On average, 9HPT dominant speeds were faster than non-dominant speeds (24.6 seconds (SD = 10.5) vs 26.6 seconds (SD = 12.3); p < 0.0001), and mean asymmetry was 4.7 seconds (SD = 10.1). Sixty-one percent had a BMI measure within 3 months of a 9HPT measure; these patients had less missing data for smoking, were less reliant on assisted devices, and less likely to have progressive disease (eTable 1).
Baseline attributes of the study population.
A subset of patients also had BMI measurements available within 3 months of baseline visit.
At baseline, dominant and non-dominant 9HPT times were strongly correlated (ρ = 0.78; 95% CI: 0.76, 0.79), both were similarly correlated with T25FW (ρ = 0.56–0.57) and PHQ-9 (ρ = 0.26–0.27) (Table 2). Asymmetry significantly correlated with all measures, with a stronger correlation with non-dominant (ρ = 0.51; 95% CI: 0.48, 0.53) than dominant (ρ = 0.28; 95% CI: 0.25, 0.31) hand times, suggesting the prior is a greater driver of asymmetry.
Spearman’s rank correlations (and 95% confidence intervals) among clinical measures at baseline. a
All correlations were statistically significant (p < 5 × 10−8).
Multivariable longitudinal associations for the three outcomes (Table 3; N = 4140 patients; 14,544 observations) showed that older age, male sex, Hispanic/Latinx ethnicity, Blacks race, ever smoking, public insurance or self-pay beneficiaries, assistive device users, progressive disease, slower T25FW times, and greater depressive symptom severity (PHQ-9) were independently associated with greater dominant hand impairment (longer/slower 9HPT times). A significant sex-disease duration interaction showed that males slowed an additional 1.5 second per 10 years (p = 0.0002) in dominant 9HPT times (Table 3, Figure 1(a)), beyond the average male–female difference of 1.6 second slower (p = 0.00078) dominant 9HPT times for males. In the raw data, treating all observations as independent, males had steeper slopes (eFigure 3(a)). However, after adjusting for covariates, the positive unadjusted slope in females (eFigure 3(a)) was no longer present (Figure 1(a)).
Longitudinal multivariable associations for the three hand impairment outcomes.
Two-sided p-values <0.05 are bolded.

(a) Plot of predicted dominant hand times by sex and disease duration and (b) Plot of predicted non-dominant hand times by sex and disease duration.
For non-dominant times, associations were similar (Table 3). Greater ADI was only significantly associated with worsening non-dominant 9HPT, and Hispanics showed a trend toward slower times (p = 0.06). The effect sizes for all attributes were larger in the non-dominant hand model compared with the dominant hand model even with comparable intercepts. There was no evidence for any interactions with disease duration. For non-dominant hand times, the raw data (eFigure 3(b)) and model-predicted estimates (Figure 1(b)) the slopes for males and females only differed by the intercept.
For asymmetry, public insurance or self-pay beneficiaries, assistive device use, progressive disease, slower T25FW, and greater PHQ-9 were significantly associated with worsening (Table 3). Sex and T25FW interacted with disease duration. In males, asymmetry increased by 1 second per decade (Table 3 and Figure 2), driven by worsening non-dominant hand times while dominant hand performance plateaued after 75 seconds (eFigure 4). The T25FW, a 1 second increase in T25FW corresponded to an additional 0.05 second increase in asymmetry per decade (Table 3 and Figure 3).

Plot of sex by disease duration interaction from asymmetry model.

Plot of timed 25-foot walk by disease duration interaction from asymmetry model.
Including BMI as a covariate (N = 2540 patients and 5090 observations) did not change results, and BMI was not associated with 9HPT outcomes (data not shown).
Discussion
This is the first repeated-measures longitudinal study examining dominant and non-dominant 9HPT times, as well as asymmetry, conducted in a large, diverse, real-world MS patient population. We observed independent associations between multiple clinical, demographic, and social determinants of health and worsening dominant 9HPT times, with more pronounced effects on non-dominant hands. Males exhibited a faster rate of decline in dominant hand function with increasing disease duration. While fewer attributes were associated with increases in asymmetry, two key interactions emerged: males and individuals with slower walking speeds experienced more pronounced discordance between hand times as a function of disease duration.
In general, these 9HPT and asymmetry associations parallel longitudinal findings for T25FW, PHQ-9, perceived global disability, and cognitive impairment, except for the null findings for BMI.8–10,29 Sexual dimorphism in MS risk, presentation, and progression is well documented, 30 and our results for sex are consistent with a secondary analysis of clinical trials data, 26 and we also demonstrated greater asymmetry in males. Not only do males on average have worse upper extremity functioning, they more rapidly accrued dominant hand dysfunction and greater discordance between hands over the disease course evidenced by statistical interactions. These findings are biologically plausible: upper- and lower-limb function share overlapping neural substrates, and lesions affecting the corticospinal tracts, cerebellum, or brainstem can impair coordination and motor control across both limbs.6,7 Neuroimaging studies further support this rationale, showing that structural damage and functional reorganization within sensorimotor networks are directly linked to upper- and lower-limb performance in PwMS.6,16,31 In parallel, the observed effects of social determinants and modifiable lifestyle factors align with their systemic influence on inflammation, neurodegeneration, and access to rehabilitation. Together, these biological and contextual mechanisms provide a coherent framework for why both clinical and social factors shape the trajectory of upper extremity dysfunction in MS.
Disease duration, assistive device use, slower walking speed, and progressive disease were independently associated with worse dominant and non-dominant hand times and greater asymmetry over time, as expected given our understanding of MS. Older age was associated with slower 9HPT times but not asymmetry, suggesting symmetric decline in upper extremity functioning with aging. Depressive symptom severity negatively impacted all three outcomes, paralleling longitudinal associations for T25FW and perceived global disability. 8 While these relationships may not reflect causal relationships, we hypothesize they may reflect shared biological processes (i.e. diminished repair capacity and accelerated age-related decline).
We observed greater hand impairment for Hispanic/Latinx and Blacks compared with non-Hispanics and Whites. A prior cross-sectional study observed Black PwMS had 2.1 second slower 9HPT (averaged across hands) compared with White PwMS. 32 We observed similar magnitudes longitudinally for the dominant and non-dominant hand. However, the underlying race difference is likely larger than what we report as potential mediators (factors downstream of race in causal paths to upper extremity impairment, for example, insurance payer, ADI) were included as covariates. In post hoc analyses of our models excluding insurance payer and ADI, the difference between Black versus White PwMS increased by 15% (2.22 seconds; 95% CI = 1.34, 3.09) and 26% (2.84s; 95% CI = 1.73, 3.94) for dominant and non-dominant 9HPT times, respectively. For Hispanic/Latinx versus non-Hispanic PwMS, dominant and non-dominant 9HPT times increased by 3% (1.91s; 95% CI = 0.14, 3.69) and 5% (3.14s; 95% = CI 0.004, 6.28), respectively. The null results for race and ethnicity in the asymmetry model were unchanged. These post hoc analyses suggest social inequities may mediate portions of the ethnoracial differences in upper extremity dysfunction. However, the relative sample sizes of these ethnoracial subgroups were modest and may be subject to sampling variability. Thorough mediation and moderation analyses are warranted to generate clearer insights.
Two of the modifiable behavioral and lifestyle factors (tobacco smoke exposure, obesity), social determinants of health (individual SES, care access reflected by insurance payer, neighborhood-level SES reflected by ADI), and other social attributes (marital status) investigated had profound associations with upper extremity function. Smoking has been consistently associated with worsening MS outcomes, including disease progression, walking speed, perceived disability, depressive symptom, processing speed, and gray matter loss, with greater deficits in current smokers.8,9,33 Our results builds on a cross-sectional study where active smokers with MS self-reported greater hand impairment. 34 Here, current and former smokers had significant worsen of 9HPT times, with greater deficits experienced by active smokers independent of disease course and lower-limb functioning. Public insurance beneficiaries also had greater deficits across the three outcomes, even with adjustment for covariates often considered as confounders (i.e. age and race/ethnicity). Since patients were receiving care at a tertiary referral facility, these associations more likely reflect adverse SES conditions rather than variation in care continuity/access, and parallels associations observed for other MS outcomes.8,9,32 A noteworthy finding was that ADI was only associated with worsening non-dominant 9HPT times and not dominant 9HPT times. We speculate this finding may reflect compensatory prioritization of dominant hand function for PwMS living in low-resourced settings. This aligns with growing evidence that adverse neighborhood SES independently contributes to worse MS outcomes. 8
There are several strengths to this study, including the opportunity to simultaneously examine independent relationships for multiple clinical, demographic, and social determinants of health in a large, diverse, real-world MS patient population. We were able to examine three outcomes, as well as interactions with time. A primary limitation is that this study was comprised of patients seen within a single tertiary healthcare enterprise, although they were residents spanning the US with a preponderance of Ohio residents (eFigure 1). Second, the modest sample sizes for ethnoracial minorities indicate that those findings should be interpreted with caution. Third, BMI data were only available for a subset of encounters for patients less reliant on assisted devices and who had a lower prevalence of progressive disease. The absence of detailed disease-modifying treatment information is another limitation. DMT use, type, and treatment history may influence the rate decline, and accounting for these factors could refine the models. However, as our cohort was drawn from long-term patients managed at tertiary MS centers, it is likely that the majority were receiving guideline-based care; nonetheless, variation in DMT efficacy, timing of initiation, and access, potentially influenced by insurance payer, may partially underlie some observed associations. Similarly, physical activity is known to affect functional outcomes in MS, yet it was not available and may have led to residual confounding. Another limitation is that this study was restricted to those with 9HPT for both sides; thus, generalizability to those with full impairment in one arm may be limited. Nonetheless, recent analyses of progressive MS trials reinforce the importance of the 9HPT, showing that it remains one of the few outcomes capable of detecting meaningful change even in non-ambulatory patients and that alternative scoring approaches (e.g. either-hand 9HPT) can increase sensitivity to worsening. 35 These complementary findings support the broader validity of 9HPT-based measures as indicators of upper-limb disability across the full spectrum of MS.
In summary, we identified multiple prominent and independent associations for longitudinal changes in upper extremity performance in MS, which may inform prognostication efforts and highlight the multifaceted nature of MS disability accrual and progression. These findings contribute to the growing literature on sex dimorphisms, ethnoracial differences, and highlight the need for quantifying health disparities and understanding the broader impact of social determinants of health in MS. Importantly, future research should also investigate other aspects of hand functioning, such as grip strength, as well as the concordance with patient-reported outcomes, which are likely to capture different components of upper extremity function that are most relevant to daily activities. Also, the influence of various disease-modifying treatments needs to be examined. Finally, it is essential to characterize the drivers of upper extremity performance in vulnerable populations, including structural barriers to care, access to rehabilitation services, and biological contributors to differential disease progression. Identifying potentially modifiable risk factors, such as active smoking, may help guide equitable care strategies aimed at preserving upper extremity function in MS patients.
Supplemental Material
sj-docx-1-msj-10.1177_13524585251384060 – Supplemental material for Longitudinal modeling of upper extremity function in multiple sclerosis: Associations for clinical and sociodemographic factors
Supplemental material, sj-docx-1-msj-10.1177_13524585251384060 for Longitudinal modeling of upper extremity function in multiple sclerosis: Associations for clinical and sociodemographic factors by Farren BS Briggs, Jiayue Yang, Karlo Toljan, Devon S Conway, Carrie M Hersh, Marisa P McGinley, Robert A Bermel, Deborah M Miller, Alessandro S De Nadai, Douglas D Gunzler and Daniel Ontaneda in Multiple Sclerosis Journal
Footnotes
Acknowledgements
We thank Mr. Scott Husak for his assistance with data curation and management.
Author Contributions
FB was responsible for conception and design of the study; all authors provided critical feedback on the study design. DO, DC, CH, MG, RB, and DM contributed to data acquisition. FB, AD, and DG guided the statistical analysis plan. JY conducted exploratory analyses. FB conducted final analyses. JY and KT drafted portions of the manuscript; FB completed the first draft and generated all figures. All authors guided interpretations and edited and approved the final draft. The authors have full access to the data used in the analyses, and they have the right to publish all data, separate and apart from the guidance of any sponsor. FB takes full responsibility for the data, the analyses and interpretation, and the conduct of the research.
Data Availability
Data are available from the Cleveland Clinic to qualified investigators.
Declaration of Conflicting Interests
The authors declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The authors received no financial support for the research, authorship, and/or publication of this article.
ORCID iDs
Supplemental material
Supplemental material for this article is available online.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
