Published on in Vol 13 (2025)

Preprints (earlier versions) of this paper are available at https://preprints.jmir.org/preprint/57599, first published .
Accelerometry-Assessed Physical Activity and Circadian Rhythm to Detect Clinical Disability Status in Multiple Sclerosis: Cross-Sectional Study

Accelerometry-Assessed Physical Activity and Circadian Rhythm to Detect Clinical Disability Status in Multiple Sclerosis: Cross-Sectional Study

Accelerometry-Assessed Physical Activity and Circadian Rhythm to Detect Clinical Disability Status in Multiple Sclerosis: Cross-Sectional Study

1Department of Neurology, Johns Hopkins University School of Medicine, 600 North Wolfe Street, Pathology 627, Baltimore, MD, United States

2Department of Biostatistics, Johns Hopkins University Bloomberg School of Public Health, Baltimore, MD, United States

Corresponding Author:

Nicole Bou Rjeily, MD


Background: Tools for measuring clinical disability status in people with multiple sclerosis (MS) are limited. Accelerometry objectively assesses physical activity and circadian rhythmicity profiles in the real-world environment and may potentially distinguish levels of disability in MS.

Objective: This study aims to determine if accelerometry can detect differences in physical activity and circadian rhythms between relapsing-remitting multiple sclerosis (RRMS) and progressive multiple sclerosis (PMS) and to assess the interplay within person between the 2 domains of physical activity (PA) and circadian rhythm (CR) in relation to MS type.

Methods: This study represents an analysis of the baseline data from the prospective HEAL-MS (home-based evaluation of actigraphy to predict longitudinal function in multiple sclerosis) study. Participants were divided into 3 groups based on the Expanded Disability Status Scale (EDSS) criteria for sustained disability progression: RRMS-Stable, RRMS-Suspected progression, and PMS. Baseline visits occurred between January 2021 and March 2023. Clinical outcome measures were collected by masked examiners. Participants wore the GT9X Link ActiGraph on their nondominant wrists for 2 weeks. After adjusting for age, sex, and BMI, a logistic regression model was fitted to evaluate the association of each accelerometry metric with odds of PMS versus RRMS. We also evaluated the association of accelerometry metrics in differentiating the 2 RRMS subtypes. The Joint and Individual Variation Explained (JIVE) model was used to assess the codependencies between the PA and CR domains and their joint and individual association with MS subtype.

Results: A total of 253 participants were included: 86 with RRMS-Stable, 82 with RRMS-Suspected progression, and 85 with PMS. Compared to RRMS, participants with PMS had lower total activity counts (β=−0.32, 95% CI −0.61 to −0.03), lower time spent in moderate to vigorous physical activity (β=−0.01, 95% CI −0.02 to −0.004), higher active-to-sedentary transition probability (β=5.68, 95% CI 1.86-9.5), lower amplitude (β=−0.0004, 95% CI −0.0008 to −0.0001), higher intradaily variability (β=4.64, 95% CI 1.45-7.84), and lower interdaily stability (β=−4.43, 95% CI −8.77 to −0.10). Using the JIVE model for PA and CR domains, PMS had higher first joint component (β=0.367, 95% CI 0.088-0.656), lower PA-1 component (β=−0.441, 95% CI −0.740 to −0.159), and lower PA-2 component (β=−0.415, 95% CI −0.717 to −0.126) compared to RRMS. No significant differences were detected between the 2 RRMS subtypes except for lower relative amplitude in those with suspected progression (β=−5.26, 95% CI −10.80 to −0.20).

Conclusions: Accelerometry detected differences in physical activity patterns between RRMS and PMS. More advanced analytic techniques may help discern differences between the 2 RRMS subgroups. Longitudinal follow-up is underway to assess the potential for accelerometry to detect or predict disability progression.

JMIR Mhealth Uhealth 2025;13:e57599

doi:10.2196/57599

Keywords



Multiple sclerosis (MS) is a chronic disease affecting the central nervous system in which those affected typically have intermittent neurologic symptoms and signs early in the course (relapsing-remitting multiple sclerosis [RRMS]) but often, subsequently (or, less commonly, from onset), slowly accumulate disability (progressive multiple sclerosis [PMS]). Currently available tools to measure disability in MS are limited. The Expanded Disability Status Scale (EDSS) is the most widely used clinical outcome measure and is considered by the Food and Drug Administration to be the gold standard for phase 3 trials in PMS. The EDSS, however, is semiquantitative, has limited reliability, and only captures a person’s state at one short point in time during a clinical visit [Inojosa H, Schriefer D, Ziemssen T. Clinical outcome measures in multiple sclerosis: a review. Autoimmun Rev. May 2020;19(5):102512. [CrossRef] [Medline]1]. People with MS may have different disability states at different times throughout the day, as symptoms can fluctuate with temperature, fatigue, stress, and other factors [Albrecht H, Wötzel C, Erasmus LP, Kleinpeter M, König N, Pöllmann W. Day-to-day variability of maximum walking distance in MS patients can mislead to relevant changes in the Expanded Disability Status Scale (EDSS): average walking speed is a more constant parameter. Mult Scler. Apr 2001;7(2):105-109. [CrossRef] [Medline]2,Christogianni A, Bibb R, Davis SL, et al. Temperature sensitivity in multiple sclerosis: an overview of its impact on sensory and cognitive symptoms. Temperature (Austin). Jul 3, 2018;5(3):208-223. [CrossRef]3]. The limitations of the EDSS lead to long delays in confirming that a person has transitioned from RRMS to PMS and also inflate the sample size and follow-up time required for MS trials for which disability progression is the primary end point [Katz Sand I, Krieger S, Farrell C, Miller AE. Diagnostic uncertainty during the transition to secondary progressive multiple sclerosis. Mult Scler. Oct 2014;20(12):1654-1657. [CrossRef]4,Inojosa H, Proschmann U, Akgün K, Ziemssen T. A focus on secondary progressive multiple sclerosis (SPMS): challenges in diagnosis and definition. J Neurol. Apr 2021;268(4):1210-1221. [CrossRef] [Medline]5].

Triaxial accelerometry is a safe and relatively inexpensive tool that may offer an objective and sensitive measure of disability in people with MS. With the use of an accelerometer worn on the wrist, real-time information about physical activity and circadian rhythmicity patterns can be collected in a person’s natural environment. Such data may allow detection of variation in activity that may be missed during clinical visits [Block VJ, Lizée A, Crabtree-Hartman E, et al. Continuous daily assessment of multiple sclerosis disability using remote step count monitoring. J Neurol. Feb 2017;264(2):316-326. [CrossRef] [Medline]6,Dalla-Costa G, Radaelli M, Maida S, et al. Smart watch, smarter EDSS: improving disability assessment in multiple sclerosis clinical practice. J Neurol Sci. Dec 15, 2017;383:166-168. [CrossRef] [Medline]7]. Several studies have looked at the use of accelerometry in MS to identify associations with sleep, fatigue, depression, or disability outcomes (EDSS and Patient-Determined Disease Steps); however, whether wrist-worn accelerometry can differentiate between people with RRMS and PMS remains an existing question [Klaren RE, Motl RW, Dlugonski D, Sandroff BM, Pilutti LA. Objectively quantified physical activity in persons with multiple sclerosis. Arch Phys Med Rehabil. Dec 2013;94(12):2342-2348. [CrossRef] [Medline]8-Macdonald E, Buchan D, Cerexhe L, Renfrew L, Sculthorpe N. Accelerometer measured physical activity and sedentary time in individuals with multiple sclerosis versus age matched controls: a systematic review and meta-analysis. Mult Scler Relat Disord. Jan 2023;69:104462. [CrossRef] [Medline]16]. We hypothesized that people with RRMS and PMS have different physical activity and circadian rhythmicity patterns. Due to the interplay between physical activity and circadian rhythm domains within person, we also aimed to study joint and individual variations in these 2 domains and their associations with MS subtype.


Participant Selection and Study Design

This study includes an analysis of the baseline data from the longitudinal observational study HEAL-MS (home-based evaluation of actigraphy to predict longitudinal function in multiple sclerosis; Figure S1 in

Multimedia Appendix 1

Illustrations of the study framework along with additional statistical analyses and sensitivity checks.

DOCX File, 796 KBMultimedia Appendix 1). Eligible participants followed at the Johns Hopkins Multiple Sclerosis Precision Medicine Center of Excellence were recruited between January 2021 and March 2023. Participants had to have a diagnosis of MS, were aged ≥40 years, had no apparent comorbidities that may limit physical activity (eg, heart failure and end-stage renal disease), had no MS relapse within the last 6 months prior to enrollment, and had baseline EDSS score ≤6.5. These criteria were chosen so that the participants with RRMS were of the age at which transition to PMS might begin within the next few years, the changes in accelerometry measures were more likely to be due to MS and not comorbid conditions, and the baseline measures were not affected by recovery from a recent relapse. Medical records were reviewed by an MS-trained neurologist to ensure eligibility.

Participants were divided into 3 groups with a target of 85 participants each: those with stable RRMS who had no suspected or confirmed progression (RRMS-Stable), those with PMS who had confirmed disability worsening on EDSS (PMS), and those with RRMS who were suspected to be progressing clinically (inferred by an MS-trained neurologist from medical record notes or personal knowledge of the participant) but did not have sustained disability worsening on exam (RRMS-Suspected progression). Confirmed disability worsening on exam was defined as an EDSS change of ≥1.0 point if baseline EDSS was ≤5.5, or of ≥0.5 points if baseline EDSS was ≥6.0; this change should be sustained for ≥24 weeks, not in the context of an explanatory relapse. The 2 RRMS groups were matched based on age (target ±2 years with prespecified ability to relax the criterion if necessary), sex at birth, race or ethnicity, and efficacy class of current disease-modifying therapy (no treatment, first-line treatment [injectables and oral therapies (except for cladribine or ofatumumab)], or higher-efficacy treatment [infusion therapies and cladribine]). The PMS group was expected to be a little older and to be prescribed different (or even no) medications due to the likely lack of neuroinflammation and relapse activity in this group. However, enrollment of a comprehensive cohort of participants with respect to sex and race or ethnicity was attempted.

Clinical Measures

The EDSS exam was conducted after the fitting of the accelerometer and was performed by a masked EDSS-trained physician. A modified multiple sclerosis functional composite (MSFC), which includes the 9-hole peg test, the timed 25-foot walk test, the high- and low-contrast letter acuity (binocular, 2.5% contrast Sloan charts), and the Symbol Digit Modalities Test, was performed by a masked, trained study team member. Participants also completed patient-reported outcomes such as the International Physical Activity Questionnaire (IPAQ).

Accelerometry Measures

Accelerometry metrics were derived from accelerometry data collected with the GT9X Link ActiGraph using a built-in triaxial accelerometer [User Guide: ActiGraph GT9X Link + ActiLife, Revision 6. Amazon S3. Feb 18, 2020. URL: https:/​/s3.​amazonaws.com/​actigraphcorp.com/​wp-content/​uploads/​2020/​03/​05155628/​ActiGraph_Link_UserGuide_E.​200.​6001_Revision6_FINAL.​pdf [Accessed 2023-09-03] 17]. All participants were instructed to wear the ActiGraph device on the wrist of their nondominant hand for 24 hours a day over a duration of 2 weeks. Accelerometers were set to capture 3D acceleration at 30 Hz with the acceleration range of ±8 G. The raw acceleration data (.gt3x) were downloaded from the device using ActiLife v6.134 Lite Edition. Binary raw activity data (Hz-level accelerometry data) were read by read.gt3x package into an R data frame (R Foundation for Statistical Computing) and transformed into 60-second epochs activity count data in 1440 minutes per day (12 AM to 11:59 PM) analytic format. The activity counts are vector magnitude-based activity counts.

The following criteria were applied to define valid days: (1) intervals of 90 minutes or longer with all minute-level activity counts equal to 0 were defined as nonwear intervals [Choi L, Liu Z, Matthews CE, Buchowski MS. Validation of accelerometer wear and nonwear time classification algorithm. Med Sci Sports Exerc. Feb 2011;43(2):357-364. [CrossRef] [Medline]18], (2) valid days were defined as those with total wear-time ≥90% of the day (≥1296 minutes of wear), and (3) each participant should have ≥3 valid days of accelerometry data.

When a wear period overlapped with the transition to or from daylight saving time, the following adjustments were made. If the clock was set 1 hour forward, 1 hour of data were missing, so imputation was used by averaging the same hour from other valid days for this participant. If the clock was set 1 hour back, duplicate data were generated for that hour; in this instance, we used the average of the duplicated hour, and the loss of 1 hour on the final day was imputed by averaging the same hour from the other valid days for this participant. The wear period overlapped with the transition to or from daylight saving time for only 5 participants.

Accelerometry measures in this study are presented as averages over all valid days and included (1) measures of volume and intensity: total activity count (TAC; daily and 2-hour specific), total log-transformed activity count (TLAC), the total daily number of nonactive minutes, the number of minutes spent in moderate to vigorous physical activity (MVPA), and the number of minutes spent in light intensity physical activity (LIPA); (2) measures of composition: MVPA/LIPA (the ratio of the time spent in MVPA over the time spent in LIPA) and MVPA/nonactivity (the ratio of the time spent in MVPA over the time spent nonactive); (3) measures of fragmentation: sedentary-to-active transition probability (SATP) and active-to-sedentary transition probability (ASTP) [Di J, Leroux A, Urbanek J, et al. Patterns of sedentary and active time accumulation are associated with mortality in US adults: the NHANES study. bioRxiv. Preprint posted online on Aug 31, 2017. [CrossRef]19,Schrack JA, Kuo PL, Wanigatunga AA, et al. Active-to-sedentary behavior transitions, fatigability, and physical functioning in older adults. J Gerontol A Biol Sci Med Sci. Mar 14, 2019;74(4):560-567. [CrossRef] [Medline]20]; and (4) diurnal landmarks of rest-activity rhythms: average log acceleration during the most active 10 hours of the day (M10), midpoint of M10, average log acceleration during the least active 5 hours of the day (L5), midpoint of L5, daytime activity ratio estimate (DARE), amplitude, relative amplitude (RA), midline estimating statistic of rhythm (MESOR), acrophase, intradaily variability (IV), and interdaily stability (IS) [Leroux A, Xu S, Kundu P, et al. Quantifying the predictive performance of objectively measured physical activity on mortality in the UK Biobank. J Gerontol A Biol Sci Med Sci. Jul 13, 2021;76(8):1486-1494. [CrossRef] [Medline]21].

The following sample-specific cut points were developed for discriminating between nonactive, LIPA, and MVPA minutes: (1) nonactive minute was defined as a minute with activity counts ≤2000, (2) LIPA minute was defined as activity counts >2000 and ≤6750, and (3) MVPA minute was defined as activity counts >6750. These cut points were specific to our cohort and defined based on an internal calibration procedure that maximized concordance between accelerometry-based estimates and the participants’ self-reported physical activity outcomes in the IPAQ (Table S1 and Figure S2 in

Multimedia Appendix 1

Illustrations of the study framework along with additional statistical analyses and sensitivity checks.

DOCX File, 796 KBMultimedia Appendix 1 [Welk GJ, Lamoureux NR, Zeng C, et al. Equating NHANES monitor-based physical activity to self-reported methods to enhance ongoing surveillance efforts. Med Sci Sports Exerc. Jun 1, 2023;55(6):1034-1043. [CrossRef] [Medline]22]). Our derived cut points were similar to those defined in other population cohorts [Welk GJ, Lamoureux NR, Zeng C, et al. Equating NHANES monitor-based physical activity to self-reported methods to enhance ongoing surveillance efforts. Med Sci Sports Exerc. Jun 1, 2023;55(6):1034-1043. [CrossRef] [Medline]22-Rhudy MB, Dreisbach SB, Moran MD, Ruggiero MJ, Veerabhadrappa P. Cut points of the Actigraph GT9X for moderate and vigorous intensity physical activity at four different wear locations. J Sports Sci. Mar 2020;38(5):503-510. [CrossRef] [Medline]24]. Rather than using predetermined cutoffs reported in the literature, we believe this method was more suitable for this study due to different cohort demographics and clinical conditions that may alter average physical activity intensity.

Statistical Analyses

All statistical analyses were conducted using R (version 4.3.1; R Foundation for Statistical Computing). Descriptive statistics of the demographics, clinical data, and accelerometry measures for each of the 3 groups were summarized as mean (SD) or median (IQR). Differences between the combined RRMS group (RRMS-Stable and RRMS-Suspected progression) and PMS group were compared using t tests and Pearson χ² test as appropriate. After adjusting for age, sex, and BMI, a logistic regression model was fitted to evaluate the association of each accelerometry metric with odds of PMS (vs RRMS). We also evaluated the association of accelerometry metrics in differentiating the RRMS subtypes (RRMS-Suspected progression vs RRMS-Stable). To account for differences in duration and timing of sleep, sensitivity analyses were performed by re-exploring the same accelerometry-derived measures within the most active 10 hours of each day (M10) rather than the 24-hour daily total.

Functional Principal Component Analysis

Functional principal component (fPC) analysis, applied to diurnal rest-activity rhythms (minute-level activity profiles), captures main patterns of temporal allocation of activity and provides important diurnal landmarks. The choice of the first 5 fPCs as measures of circadian rhythmicity is based on 2 factors. First, they capture most of the diurnal variability (76%) in our data. Second, the use of fPC analysis in accelerometry studies of diurnal rest-activity rhythms in large national cohorts such as the National Health and Nutrition Examination Survey (NHANES) and the UK Biobank mostly focuses on the first 4 fPCs since functional principal components beyond the 4th typically capture a relatively low proportion of variance in the data, and their shapes tend to be very difficult to interpret due to their very dynamic shape [Xiao Q, Lu J, Zeitzer JM, Matthews CE, Saint-Maurice PF, Bauer C. Rest-activity profiles among U.S. adults in a nationally representative sample: a functional principal component analysis. Int J Behav Nutr Phys Act. Mar 24, 2022;19(1):32. [CrossRef] [Medline]25-Stein MJ, Baurecht H, Sedlmeier AM, et al. Association between circadian physical activity patterns and mortality in the UK Biobank. Int J Behav Nutr Phys Act. Sep 1, 2023;20(1):102. [CrossRef] [Medline]27]. Nevertheless, in our case, the fifth fPC explained a significant proportion of the total variability (4%) and was also included. Figure S3 in

Multimedia Appendix 1

Illustrations of the study framework along with additional statistical analyses and sensitivity checks.

DOCX File, 796 KBMultimedia Appendix 1 shows that the additional fPCs beyond the 5th did not provide meaningful explanations of variability.

Joint and Individual Variation Explained

The Joint and Individual Variation Explained (JIVE) is an integrative dimension reduction technique that can be applied to multiple features grouped within several domains [Di J, Spira A, Bai J, et al. Joint and individual representation of domains of physical activity, sleep, and circadian rhythmicity. Stat Biosci. Jul 2019;11(2):371-402. [CrossRef] [Medline]28]. The physical activity (PA) domain was characterized by measures of total volume of physical activity, times spent in different intensities of activity and composition ratios of those times, measures of activity fragmentation/continuity, and temporal-local measures of total activity presented in 2-hour bins. The circadian rhythmicity (CR) domain was characterized by measures of the strength of circadian rhythms (RA and related components such as M10 and timing of M10), measures of variability including IS and IV, parametric measures extracted from the cosinor model (MESOR, amplitude, acrophase), and the first 5 fPCs. There is a significant amount of interdependence between the PA and CR domains since the measures belong to the same participants over the same time frame. Variations in PA and CR domains that are independent from each other are also expected. JIVE is a model developed to explore such situations and to decompose the joint and individual effects [Lock EF, Hoadley KA, Marron JS, Nobel AB. Joint and individual variation explained (JIVE) for integrated analysis of multiple data types. Ann Appl Stat. Mar 1, 2013;7(1):523-542. [CrossRef] [Medline]29]. Figure 1 demonstrates the conceptual representation of the benefits of JIVE to better understand the association between MS subtype and the joint and individual variation of PA and CR domains. We applied JIVE to the domains of PA and CR to create three groups of latent variables that explain (1) joint variation shared across PA and CR, (2) individual variation specific to each of the domains, and (3) remaining unexplained variation.

Figure 1. Conceptual representation of associations between multiple sclerosis (MS) subtype and the physical activity (PA) and circadian rhythm (CR) domains performed in standard regression modeling (left) versus the associations between MS subtype and 3 sets of independent (uncorrelated) latent variables representing joint-PA-CR, individual-PA, and individual-CR information after JIVE (Joint and Individual Variation Explained) decomposition (right). Regression analysis with JIVE components as predictors can reveal and distinguish joint and individual associations between the PA and CR domains and MS subtype.

Ethical Considerations

The study was approved by the institutional review board at Johns Hopkins University (IRB00243681). All participants provided written informed consent and had the ability to opt out at any point. All data are stored on REDCap (Research Electronic Data Capture), and only study members who are approved by the institutional review board have access to the data. The data are deidentified to the greatest extent possible. All participants received US $40 as compensation for the baseline visit to offset study-related travel and parking costs. All study procedure costs, including shipping of the ActiGraph, were covered as part of the study.


Demographics and Clinical Outcomes

A total of 275 people with MS completed a baseline study visit as part of the HEAL-MS study; 253 had retrievable accelerometry data and were included in this baseline analysis (16 never returned the device, while 6 had insufficient valid days). Participants were divided into three groups: 86 with RRMS-Stable, 82 with RRMS-Suspected progression, and 85 with PMS. Baseline visits occurred between January 2021 and March 2023. Table 1 shows the demographic and clinical data for the 3 groups.

Table 1. Demographics and clinical data collected at the baseline visit for HEAL-MSa participants.
RRMSbPMSc (n=85)P valued
RRMS-Stable (n=86)RRMS-Suspected progression (n=82)
Demographics
 Age (years), mean (SD)53 (7)54 (8)57 (9)<.001e
 Sex (female), n (%)61 (71)61 (74)57 (67).36
 Race, n (%).67
  White64 (74)63 (77)67 (79)
  Black19 (22)14 (17)13 (15)
  Other3 (4)5 (6)5 (6)
 Hispanic or Latino, n (%)4 (5)3 (4)3 (4).79
 BMI, mean (SD)29.4 (6.4)28.2 (5.6)28 (6.2).35
 Employment status, n (%).18
  Full time45 (58)37 (51)32 (43)
  Part time4 (5)6 (8)7 (9)
  Homemaker7 (9)2 (3)4 (5)
  Retired13 (17)13 (18)18 (24)
  Unemployed3 (4)3 (4)0 (0)
  Disability6 (8)12 (16)14 (19)
Clinical data
 EDSSf score, median (IQR)2.0 (1.5-3.0)3.0 (2.0-3.5)5.0 (3.5-6.0)<.001e
 MSFCg, median (IQR)
  25-foot walk4.5 (4.0-5.1)5.3 (4.3-6.2)7.2 (5.8-9.5)<.001e
  Dh-9hpgi21.3 (19.4-24.0)23.8 (20.9-27.3)26.9 (22.1-34.4)<.001e
  NDj-9hpg22.7 (20.8-25.2)25.5 (21.9-28.9)27.4 (23.5-34.6)<.001e
  SDMTk51 (47-59)48 (38-57)42 (34-50)<.001e
 Tremorl, median (IQR)0 (0-0.5)0 (0-0.5)0 (0-1).20

aHEAL-MS: home-based evaluation of actigraphy to predict longitudinal function in multiple sclerosis.

bRRMS: relapsing-remitting multiple sclerosis.

cPMS: progressive multiple sclerosis.

dP value from t test or Pearson χ² test comparing RRMS (RRMS-Stable and RRMS-Suspected progression) vs PMS.

eP<.05

fEDSS: Expanded Disability Status Scale.

gMSFC: multiple sclerosis functional composite.

hD: dominant hand.

i9hpg: 9-hole peg test.

jND: nondominant hand.

kSDMT: Symbol Digit Modalities Test.

lBain Score for Tremor Severity (BSTS) was used for tremor evaluation.

Accelerometry Measures Over 24 Hours

Figure S4 in

Multimedia Appendix 1

Illustrations of the study framework along with additional statistical analyses and sensitivity checks.

DOCX File, 796 KBMultimedia Appendix 1 shows the raw accelerometry data for 1 participant over the course of the week. The fPCs captured important temporal or diurnal landmarks (Figure 2) and are comparable to those captured in large national surveys such as UK Biobank and NHANES [Xiao Q, Lu J, Zeitzer JM, Matthews CE, Saint-Maurice PF, Bauer C. Rest-activity profiles among U.S. adults in a nationally representative sample: a functional principal component analysis. Int J Behav Nutr Phys Act. Mar 24, 2022;19(1):32. [CrossRef] [Medline]25,Stein MJ, Baurecht H, Sedlmeier AM, et al. Association between circadian physical activity patterns and mortality in the UK Biobank. Int J Behav Nutr Phys Act. Sep 1, 2023;20(1):102. [CrossRef] [Medline]27]. Accelerometry metrics for the 3 groups are shown in Table 2. TAC and MVPA, the 2 main measures of daily physical activity volume, as well as M10, amplitude, and MESOR, measures of circadian rhythmicity, were significantly different between the RRMS and PMS groups, with greater intensity of activity in the RRMS groups compared to the PMS group. ASTP was also significantly different, with higher values (more frequent transitions to sedentary time) in the PMS group. When TAC was divided into 2-hour windows, it was significantly lower in the PMS group at specific times of the day (6 AM to 12 PM and 2 PM to 10 PM; Figure 3). Figure S5 in

Multimedia Appendix 1

Illustrations of the study framework along with additional statistical analyses and sensitivity checks.

DOCX File, 796 KB
Multimedia Appendix 1
demonstrates the correlation plot for all accelerometry measures; the green box demonstrates the temporal profiles of the circadian measures associated with TAC over 24 hours. MESOR, amplitude, and fPC1 are highly correlated with TAC between 9 AM and 8 PM.

Figure 2. First 5 functional principal components (fPCs). fPC 2, 3, and 4 values are multiplied by −1.
Table 2. Mean (SD) baseline accelerometry metrics over 24 hours for HEAL-MSa participants.
RRMSbPMSc (n=85)P valued
RRMS-Stable (n=86)RRMS-Suspected progression (n=82)
Metrics of physical activity volume and fragmentation
TACf (x103)2113 (566)2181 (724)1909 (717).008e
TLACg6277 (868)6278 (933)6196 (999).51
Nonactive minutes1067 (102)1052 (114)1087 (130).08
MVPAh59 (34)61 (52)41 (43)<.001e
LIPAi314 (84)327 (83)312 (105).52
SATPj0.09 (0.03)0.09 (0.03)0.09 (0.03).21
ASTPk0.27 (0.06)0.26 (0.06)0.30 (0.11).003e
Metrics of circadian rhythm
M10l2482 (665)2540 (844)2161 (806)<.001e
Midpoint of M10m2:20 PM (85)2:15 PM (86)2:14 PM (77).72
L5n111 (67)135 (109)146 (124).09
Midpoint of L5m3:58 AM (68)3:58 AM (80)4:05 AM (76).57
RAo0.91 (0.05)0.89 (0.07)0.87 (0.09)<.001e
DAREp0.70 (0.05)0.70 (0.05)0.69 (0.05).049e
Amplitude2290 (1081)2218 (847)1923 (839).008e
MESORq1313 (623)1264 (390)1120 (429).01e
Acrophase14.56 (1.92)14.61 (1.35)14.76 (1.31).40
IVr0.52 (0.11)0.49 (0.09)0.49 (0.08).02e
ISs0.26 (0.07)0.27 (0.07)0.25 (0.06).08
fPC1t1912 (17,574)4242 (22,747)−5206 (21,441).003e
fPC2933 (11,249)−1596 (11,922)507 (10,181).59
fPC3235 (9505)−166 (8702)−108 (6432).89
fPC4−1489 (6339)−585 (8694)451 (4989).10
fPC549 (6677)−489 (5772)677 (5744).27

aHEAL-MS: home-based evaluation of actigraphy to predict longitudinal function in multiple sclerosis.

bRRMS: relapsing-remitting multiple sclerosis.

cPMS: progressive multiple sclerosis.

dP value from t test comparing RRMS (RRMS-Stable and RRMS-Suspected progression) vs PMS.

eP<.05

fTAC: total activity count.

gTLAC: total log-transformed activity count.

hMVPA: moderate to vigorous physical activity.

iLIPA: light intensity physical activity.

jSATP: sedentary-to-active transition probability.

kASTP: active-to-sedentary transition probability.

lM10: average log acceleration during the most active 10 hours of the day.

mMidpoints of M10 and L5 are represented as time of the day in hour:minutes (±SD in minutes).

nL5: average log acceleration during the least active 5 hours of the day.

oRA: relative amplitude.

pDARE: daytime activity ratio estimate.

qMESOR: midline estimating statistic of rhythm.

rIV: intradaily variability.

sIS: interdaily stability.

tfPC: functional principal component.

Figure 3. Total activity counts (TAC) for each 2-hour interval over the course of 24 hours for PMS versus RRMS. The 2-sample t test P values show the 2-hour intervals when TAC was significantly different between the 2 groups. As an additional reference, the logistic regression P values are shown with the TAC 2-hour intervals used as predictors, after adjusting for age, sex, and BMI. MS: multiple sclerosis; PMS: progressive multiple sclerosis; RRMS: relapsing-remitting multiple sclerosis.

In multivariable-adjusted models (Table 3), participants with PMS had lower TAC (β=−0.32, 95% CI −0.61 to −0.03; P=.03), lower MVPA (β=−.01, 95% CI −0.02 to −0.004; P=.004), lower M10 (β=−0.0006, 95% CI −0.001 to −0.0002; P=.003), and higher ASTP (β=5.68, 95% CI 1.86-9.50; P=.004) when compared to those with RRMS. PMS (vs combined RRMS) was also associated with lower TAC from 8 AM to 12 PM and from 2 PM to 8 PM (Figure 3). Lower values of compositional ratios MVPA/LIPA and MVPA/nonactive were also associated with PMS (vs RRMS), indicating that individuals with PMS tend to spend relatively less time in moderate/vigorous activity compared to the time in light intensity activity or nonactivity (Table 3). Participants with PMS also had lower amplitude (β=−0.0004, 95% CI −0.0008 to −0.0001; P=.01), lower MESOR (β=−0.0009, 95% CI −0.002 to −0.0002; P=.01), higher IV (β=4.64, 95% CI 1.45-7.84; P=.004), lower IS (β=−4.43, 95% CI −8.77 to −0.10; P=.045), and lower fPC1 (β=−0.39, 95% CI −0.69 to −0.09; P=.01) (Table 3). There were no significant associations of accelerometry measures between RRMS subgroups except for a difference in RA (β=−5.26, 95% CI −10.80 to −0.20; P=.049), with lower RA in RRMS-Suspected progression (Table 3; Figure S6 in

Multimedia Appendix 1

Illustrations of the study framework along with additional statistical analyses and sensitivity checks.

DOCX File, 796 KBMultimedia Appendix 1).

Table 3. Multivariable logistic regression model to evaluate the association of each accelerometry metric with multiple sclerosis subtype. The multivariable logistic regression models were adjusted for age, sex, and BMI.
PMSa vs RRMSb (reference group)RRMS-Suspected progression vs RRMS-Stable (reference group)
EstimateP value95% CIEstimateP value95% CI
TACd,e−0.32.03c−0.61 to −0.030.14.42−0.20 to 0.50
TLACd,f−0.04.79−0.31 to 0.230.05.77−0.28 to 0.38
Nonactive minutes0.002.16−0.0007 to 0.004−0.002.29−0.005 to 0.001
MVPAg−0.01.004c−0.02 to −0.0040.002.70−0.006 to 0.009
LIPAh−0.0006.69−0.004 to 0.0020.002.25−0.002 to 0.006
SATPi−3.004.59−13.84 to 7.831.60.81−11.16 to 14.47
ASTPj5.68.004c1.86 to 9.50−4.44.13−10.47 to 1.05
MVPA/LIPA−4.25.002c−6.99 to −1.51−0.31.81−2.91 to 2.27
MVPA/nonactive−7.90.02c−14.53 to −1.272.35.46−3.73 to 9.04
LIPA/nonactive−0.36.74−2.50 to 1.781.74.22−0.98 to 4.58
M10k−0.0006.003c−0.001 to −0.00020.0001.63−0.0003 to 0.0005
Midpoint of M10−0.0004.82−0.004 to 0.003−0.0006.75−0.004 to 0.003
L5l0.002.11−0.0005 to 0.0050.004.070 to 0.008
Midpoint of L50.0009.71−0.004 to 0.005−0.004.16−0.01 to 0.002
RAm−6.00.002c−9.73 to −2.26−5.26.049c−10.80 to −0.20
DAREn−5.80.04c−11.25 to −0.34−0.34.91−6.42 to 5.74
Amplitude−0.0004.01c−0.0008 to −0.0001−0.0001.59−0.0004 to −0.0002
MESORo−0.0009.01c−0.002 to −0.0002−0.0002.48−0.0009 to 0.0004
Acrophase0.08.37−0.10 to 0.270.02.84−0.17 to 0.20
IVp4.64.004c1.45 to 7.84−0.30.88−4.21 to 3.62
ISq−4.43.045c−8.77 to −0.100.53.82−3.97 to 5.03
fPC1r−0.39.01c−0.69 to −0.090.14.39−0.19 to 0.47
fPC20.04.79−0.23 to 0.31−0.25.12−0.58 to 0.07
fPC30.08.57−0.19 to 0.350.009.96−0.30 to 0.32
fPC40.19.18−0.08 to 0.460.11.51−0.21 to 0.42
fPC50.15.27−0.12 to 0.42−0.05.74−0.36 to 0.26

aPMS: progressive multiple sclerosis.

bRRMS: relapsing-remitting multiple sclerosis.

cP<.05.

dRegression inputs were scaled by dividing TAC and TLAC values by 1 SD.

eTAC: total activity count.

fTLAC: total log-transformed activity count.

gMVPA: moderate to vigorous physical activity.

hLIPA: light intensity physical activity.

iSATP: sedentary-to-active transition probability.

jASTP: active-to-sedentary transition probability.

kM10: average log acceleration during the most active 10 hours of the day.

lL5: average log acceleration during the least active 5 hours of the day.

mRA: relative amplitude.

nDARE: daytime activity ratio estimate.

oMESOR: midline estimating statistic of rhythm.

pIV: intradaily variability.

qIS: interdaily stability.

rfPC: functional principal component.

Accelerometry Measures Within M10

To account for differences in duration and timing of sleep, the accelerometry metrics were explored within the M10 period. The results are summarized in Table S2 in

Multimedia Appendix 1

Illustrations of the study framework along with additional statistical analyses and sensitivity checks.

DOCX File, 796 KBMultimedia Appendix 1. TAC, MVPA, and ASTP remained significantly different between the RRMS and PMS groups. Additionally, the PMS group displayed more nonactive time compared to the RRMS group. Logistic regression identified significant associations between the MS subtype (PMS vs RRMS) and TAC, MVPA, nonactive minutes, ASTP, and MVPA/LIPA (Table S3 in

Multimedia Appendix 1

Illustrations of the study framework along with additional statistical analyses and sensitivity checks.

DOCX File, 796 KB
Multimedia Appendix 1
). No significant associations were identified when comparing the 2 RRMS subgroups.

The JIVE Components

We also applied JIVE, a novel integrative dimension reduction technique. Figure 1 presents the conceptualized diagram of JIVE decomposition and how JIVE can be used to better understand codependencies between PA and CR domains and their joint and individual association with MS subtype. Figure 4 shows the estimated JIVE components obtained from the accelerometry-derived domains of PA and CR. The directions of loadings (+ or −) and squares of loadings for the estimated JIVE joint component and individual PA and CR components are shown in Tables 4 and 5, and Table S4 in

Multimedia Appendix 1

Illustrations of the study framework along with additional statistical analyses and sensitivity checks.

DOCX File, 796 KBMultimedia Appendix 1, respectively. Higher positive scores of the first joint JIVE component primarily represent lower volume of physical activity (captured by TAC) and weaker circadian rhythms (captured by M10, MESOR, amplitude, and fPC1) (Table 4). Higher positive individual PA-1 scores capture less frequent transitions to active behavior (SATP), less time spent in LIPA (LIPA, LIPA/nonactive, TLAC), and more time spent nonactive. Higher positive individual PA-2 scores capture more time spent in MVPA (MVPA, MVPA/LIPA, MVPA/nonactive) and more activity in the morning and evening hours (TAC 8 to 10 AM and TAC 6 PM to 12 AM) (Table 5). Joint JIVE component and individual PA-1 and PA-2 components were found to be significantly associated with MS subtype (Table 6). No significant associations were identified when comparing the 2 RRMS subgroups.

Figure 4. Joint and Individual Variation Explained (JIVE) by the 2 JIVE accelerometry-derived domains of circadian rhythm (CR) and physical activity (PA).
Table 4. JIVEa joint component with variables with more than 5% of proportional variation.
VariableJoint component loadingRelative magnitude of joint component loading
M10b−0.26653850.0710428
fPC1c−0.26082760.0680310
MESORd−0.25195080.0634792
Amplitude−0.24089920.0580324
TACe−0.22565100.0509184

aJIVE: Joint and Individual Variation Explained.

bM10: average log acceleration during the most active 10 hours of the day.

cfPC: functional principal component.

dMESOR: midline estimating statistic of rhythm.

eTAC: total activity count.

Table 5. Individual JIVEa components of physical activity with variables with more than 5% of proportional variation.
VariablePA component loadingRelative magnitude of PA component loading
Physical activity component 1
SATPb−0.44937030.2019337
LIPAc−0.41411490.1714912
LIPA/nonactive−0.36649750.1343204
TLACd−0.33659270.1132946
MVPA/LIPAe0.32931820.1084505
Nonactive minutes0.24481850.0599361
Physical activity component 2
MVPA/LIPA0.47289960.2236340
MVPA0.39701610.1576218
MVPA/nonactive0.36734620.1349432
TACf 8 PM to 10 PM0.33415790.1116615
TAC 8 AM to 10 AM0.30994020.0960629
TAC 10 PM to 12 AM0.24771340.0613620
TAC 6 PM to 8 PM0.24302020.0590588
Physical activity component 3
TAC 12 AM to 2 AM0.60356900.3642955
TAC 10 PM to 12 AM0.45509180.2071086
TAC 4 AM to 6 AM−0.35006180.1225433
TAC 2 AM to 4 AM−0.30850890.0951777
Physical activity component 4
ASTPg−0.51348490.2636667
TAC 8 PM to 10 PM0.38035490.1446698
TAC 2 PM to 4 PM−0.35885930.1287800
SATP−0.32833070.1078011
TLAC−0.27391820.0750312
TAC 10 AM to 12 PM0.25166890.0633372

aJIVE: Joint and Individual Variation Explained.

bSATP: sedentary-to-active transition probability.

cLIPA: light intensity physical activity.

dTLAC: total log-transformed activity count.

eMVPA: moderate to vigorous physical activity.

fTAC: total activity count.

gASTP: active-to-sedentary transition probability.

Table 6. Multivariable logistic regression model to evaluate the association of each JIVEa component with multiple sclerosis subtype. The multivariable logistic regression models were adjusted for age, sex, and BMI.
PMSb vs RRMSc (reference group)RRMS-Suspected progression vs RRMS-Stable (reference group)
EstimateP value95% CIEstimateP value95% CI
Joint component0.3665.01d0.0881 to 0.6558−0.1266.44−0.4525 to 0.1897
PA-1e−0.441.003d−0.7397 to −0.159−0.0183.91−0.3404 to 0.3034
PA-2−0.4147.006d−0.7174 to −0.1255−0.3066.07−0.6502 to 0.0207
PA-3−0.0527.71−0.3282 to 0.2204−0.2622.11−0.5887 to 0.0529
PA-40.235.11−0.046 to 0.52720.1849.25−0.129 to 0.5098
CR-1f0.2588.07−0.0162 to 0.54110.0107.95−0.3044 to 0.3275
CR-2−0.2732.052−0.5536 to 0.0001−0.2991.07−0.6344 to 0.019
CR-3−0.142.32−0.4266 to 0.1350.0815.61−0.2345 to 0.4031

aJIVE: Joint and Individual Variation Explained.

bPMS: progressive multiple sclerosis.

cRRMS: relapsing-remitting multiple sclerosis.

dP<.05.

ePA: individual physical activity component.

fCR: individual circadian rhythm component.


Principal Findings

This cross-sectional analysis was conducted to determine if a set of detailed accelerometry-derived summary measures of physical activity and circadian rhythmicity distinguishes between MS subtypes. As an extension of the work we and others had previously published, in which physical activity counts were on average lower in those with MS who were more disabled [Keller JL, Tian F, Fitzgerald KC, et al. Using real-world accelerometry-derived diurnal patterns of physical activity to evaluate disability in multiple sclerosis. J Rehabil Assist Technol Eng. 2022;9. [CrossRef] [Medline]9-Ezeugwu V, Klaren RE, Hubbard EA, Manns PT, Motl RW. Mobility disability and the pattern of accelerometer-derived sedentary and physical activity behaviors in people with multiple sclerosis. Prev Med Rep. 2015;2:241-246. [CrossRef] [Medline]11], we show here that participants with PMS exerted less total daily volume of physical activity than those with RRMS, and that these differences seemed more pronounced during specific times of the day. Sensitivity analyses to account for differences in timing and duration of sleep using the accelerometry metrics during the M10 period demonstrated very similar findings, which, if confirmed in subsequent studies, might help improve the feasibility of data collection (eg, if participants are not eager to wear the device overnight).

This work newly demonstrates that beyond average overall activity measures, specific PA indices had stronger discriminative properties than others in distinguishing PMS from RRMS. For example, those with PMS were more likely to transition from active to nonactive behavior (higher ASTP) and had less time spent in MVPA, instead favoring a greater amount of time in which they are inactive or only lightly active. ASTP may be a complementary measure of functionality to the measures of volume of daily activity such as TAC and MVPA. In the general population, higher ASTP is associated with worse health and functional status and higher fatiguability in older adults [Schrack JA, Kuo PL, Wanigatunga AA, et al. Active-to-sedentary behavior transitions, fatigability, and physical functioning in older adults. J Gerontol A Biol Sci Med Sci. Mar 14, 2019;74(4):560-567. [CrossRef] [Medline]20]. Our findings are consistent, as people with PMS have overall greater disability and symptoms. Therefore, it may be useful for future studies to explore novel measures of fragmentation of physical activity to determine indices that could better distinguish MS subtype.

Although differences in patterns of activity were demonstrated between the PMS and RRMS subgroups, exploring accelerometry-measured sleep and circadian/diurnal rhythmicity may discern even more significant differences. For example, we previously demonstrated that among people with MS, diurnal variability is greater in those with higher EDSS [Keller JL, Tian F, Fitzgerald KC, et al. Using real-world accelerometry-derived diurnal patterns of physical activity to evaluate disability in multiple sclerosis. J Rehabil Assist Technol Eng. 2022;9. [CrossRef] [Medline]9]. Herein, we confirmed that differences in TAC within specific times of the day (8 AM to 12 PM; 2 PM to 8 PM) appeared to drive the daily average activity measures. Participants with PMS also had higher IV and lower IS, further demonstrating the importance of detecting differences in diurnal rhythmicity between the 2 groups. The observed differences in our cohort herein serve as motivation for further exploration with more advanced analytical techniques beyond the use of averages of daily CR measures. Our next steps will include further investigation of more complex measures of fragmentation of physical activity, distributional and temporally local distributional characteristics of activity, and parameters of sleep and circadian disruptions [Di J, Leroux A, Urbanek J, et al. Patterns of sedentary and active time accumulation are associated with mortality in US adults: the NHANES study. bioRxiv. Preprint posted online on Aug 31, 2017. [CrossRef]19,Di J, Spira A, Bai J, et al. Joint and individual representation of domains of physical activity, sleep, and circadian rhythmicity. Stat Biosci. Jul 2019;11(2):371-402. [CrossRef] [Medline]28,Salomon A, Galperin I, Buzaglo D, et al. Fragmentation, circadian amplitude, and fractal pattern of daily-living physical activity in people with multiple sclerosis: is there relevant information beyond the total amount of physical activity? Mult Scler Relat Disord. Dec 2022;68:104108. [CrossRef] [Medline]30-Ghosal R, Ghosh SK, Schrack JA, Zipunnikov V. Distributional outcome regression and its application to modelling continuously monitored heart rate and physical activity. arXiv. Preprint posted online on Jan 26, 2023. [CrossRef]32].

An application of JIVE demonstrated a significant overlap between variation in physical activity and circadian rhythmicity domains. Joint component 1 was significant and was primarily driven by the total volume of PA along with measures of strength of CR (M10, fPC1, MESOR, amplitude), suggesting that a joint study of these 2 domains may provide a more accurate picture and capture the interdependence between the domains and their combined association with MS subtype. The significant association of MS subtype with the individual PA-1 and PA-2 components confirmed the results from the standalone associations of individual measures with MS subtype. Specifically, more fragmented physical activity profiles with lower times spent in MVPA and LIPA can be used to differentiate PMS and RRMS.

Accelerometry did not distinguish between the RRMS subgroups cross-sectionally. Correctly classifying people with MS as having transitioned from RRMS to PMS often takes several years and is usually done retrospectively [Ziemssen T, Bhan V, Chataway J, et al. Secondary progressive multiple sclerosis: a review of clinical characteristics, definition, prognostic tools, and disease-modifying therapies. Neurol Neuroimmunol Neuroinflamm. Jan 2023;10(1):381-384. [CrossRef] [Medline]33]. In designing this cohort, we hypothesized that data acquired in the early phase of cohort follow-up would predict subsequent worsening as measured by EDSS and that accelerometry may help to reduce the time of diagnostic uncertainty. Since we report only cross-sectional data, we did not expect substantial differences between the RRMS groups, especially as these groups were matched for key characteristics (age, sex, race or ethnicity, and disease modifying therapy class). Longitudinal tracking of the participants will allow for tracking within-person changes in disability trajectories. Whether including additional features of accelerometry data, such as more advanced measures of sleep or circadian rhythm, will better distinguish RRMS subpopulations even at a single time point will be an important next step.

Finally, although our findings are currently not immediately applicable in clinical practice, the study suggests it may be possible to alert clinicians about certain patterns in 24-hour behaviors that may affect the lives of people with MS. Though it can be difficult to clinically distinguish PMS from RRMS, accelerometry measures appear to detect distinct patterns of activity. People with PMS seem to have increased variability and fragmentation of their activity patterns, and intermittently reduced function is commonly reported, subjectively, by people with more advanced MS during clinical visits. This increases the confidence that accelerometry can in fact measure clinically meaningful patterns in people with MS. It is important for future efforts to recapitulate these findings with commercial accelerometers in order to support scalability and to translate the findings into results that clinicians and patients can easily understand and monitor over time. For instance, the increased variability in people with PMS indicates that their experience of their disease may be very different from the clinical picture evident during their clinical visits. A discussion with patients regarding their physical activity, supplemented by accelerometry data, could shed light on other aspects of their lives that should be acknowledged. Understanding when patients are likely to be at their peak activity levels may help plan their day better (when to plan effortful tasks or schedule breaks) and may inform lifestyle modifications such as exercise or even physical therapy interventions to evaluate for improvements.

Limitations

There are limitations to the study. The classification of RRMS as stable versus suspected progression is subjective, and it is likely that some of these individuals were misclassified, where a stable person was deemed as having suspected progression or vice versa, which could contribute to the lack of ability to detect differences in accelerometry measures between the 2 RRMS groups. We chose to create this subclassification to mirror the “clinical intuition” whereby physicians and, often, their patients worry they are progressing even if they do not meet the criteria for sustained EDSS worsening, and we wanted to be sure that sufficient numbers of such people were enrolled in the study. The longitudinal trajectories of the participants with RRMS will be informative; since all are at risk, based on age, of transitioning to a PMS phenotype, we can evaluate this risk as it relates to accelerometry measures not only using these subjective subtypes but also using a more agnostic, data-driven approach within the participants with RRMS subgroup as a whole. Another potential bias related to our study is the Hawthorne effect, where participants could potentially alter their physical activity because they were being monitored. Although this risk cannot be eliminated, it was minimized by disabling the features that display activity levels on the ActiGraph monitors (ie, the number of steps); only the time of day was displayed. It is possible that relevant comorbidities were missed during chart review although this is unlikely to affect a substantial enough proportion of individuals to meaningfully have influenced results. Finally, this is a cross-sectional study, and we will be able to better explore the dynamic changes in activity measured by accelerometry as well as fluctuations in the EDSS with the longitudinal follow-up, which will provide deeper insight into the evolution of MS progression measures.

Conclusions

In conclusion, accelerometry-derived measures captured differences in physical activity and circadian rhythm profiles between people with PMS and RRMS. Further analytic approaches that consider more sophisticated approaches to time (such as circadian patterns) may help further refine the ability to detect differences between these subgroups at a single time point, whereas longitudinal follow-up is most likely to allow for the identification of accelerometry patterns that portend subsequent disability worsening as defined by EDSS. We specifically hypothesize that the measures that distinguish PMS from RRMS at baseline will be those to change first in the participants with RRMS who transition to PMS in the longitudinal follow-up of this cohort. Enabling the identification of MS progression or risk may help with clinical decision-making (such as choosing an appropriate disease-modifying therapy for a given person) or may enable advancements in telemedicine by tracking functions remotely. More importantly, these technologies will likely facilitate therapeutic development for PMS by improving disability detection and worsening, alone or in combination with other metrics such as neuroimaging, which will enable faster and smaller trials of promising interventions.

Acknowledgments

This study was funded by the National Institute of Nursing Research, National Institutes of Health (grant 5R01NR018851). The funder played no role in the study design, data collection, analysis and interpretation of data, or the writing of this manuscript.

We would also like to thank MS4MS for their support of our wearables program.

Data Availability

The datasets generated or analyzed during the study are available from the corresponding author on reasonable request and with appropriate approvals in place. The source code associated with this study is available upon request and will be made publicly accessible on GitHub upon the publication of this paper.

Conflicts of Interest

EMM discloses research funding from Biogen and Roche/Genentech, consulting for BeCareLink, LLC, and royalties for editorial duties for UpToDate. PAC is/has been Principal Investigator on grants to Johns Hopkins University from Genentech and Principia. He has served as a paid consultant for Lilly, Avidea Technologies, Idorsia, Nervgen, and Biogen. All other authors have no disclosures.

Multimedia Appendix 1

Illustrations of the study framework along with additional statistical analyses and sensitivity checks.

DOCX File, 796 KB

  1. Inojosa H, Schriefer D, Ziemssen T. Clinical outcome measures in multiple sclerosis: a review. Autoimmun Rev. May 2020;19(5):102512. [CrossRef] [Medline]
  2. Albrecht H, Wötzel C, Erasmus LP, Kleinpeter M, König N, Pöllmann W. Day-to-day variability of maximum walking distance in MS patients can mislead to relevant changes in the Expanded Disability Status Scale (EDSS): average walking speed is a more constant parameter. Mult Scler. Apr 2001;7(2):105-109. [CrossRef] [Medline]
  3. Christogianni A, Bibb R, Davis SL, et al. Temperature sensitivity in multiple sclerosis: an overview of its impact on sensory and cognitive symptoms. Temperature (Austin). Jul 3, 2018;5(3):208-223. [CrossRef]
  4. Katz Sand I, Krieger S, Farrell C, Miller AE. Diagnostic uncertainty during the transition to secondary progressive multiple sclerosis. Mult Scler. Oct 2014;20(12):1654-1657. [CrossRef]
  5. Inojosa H, Proschmann U, Akgün K, Ziemssen T. A focus on secondary progressive multiple sclerosis (SPMS): challenges in diagnosis and definition. J Neurol. Apr 2021;268(4):1210-1221. [CrossRef] [Medline]
  6. Block VJ, Lizée A, Crabtree-Hartman E, et al. Continuous daily assessment of multiple sclerosis disability using remote step count monitoring. J Neurol. Feb 2017;264(2):316-326. [CrossRef] [Medline]
  7. Dalla-Costa G, Radaelli M, Maida S, et al. Smart watch, smarter EDSS: improving disability assessment in multiple sclerosis clinical practice. J Neurol Sci. Dec 15, 2017;383:166-168. [CrossRef] [Medline]
  8. Klaren RE, Motl RW, Dlugonski D, Sandroff BM, Pilutti LA. Objectively quantified physical activity in persons with multiple sclerosis. Arch Phys Med Rehabil. Dec 2013;94(12):2342-2348. [CrossRef] [Medline]
  9. Keller JL, Tian F, Fitzgerald KC, et al. Using real-world accelerometry-derived diurnal patterns of physical activity to evaluate disability in multiple sclerosis. J Rehabil Assist Technol Eng. 2022;9. [CrossRef] [Medline]
  10. Abbadessa G, Lavorgna L, Miele G, et al. Assessment of multiple sclerosis disability progression using a wearable biosensor: a pilot study. J Clin Med. Mar 10, 2021;10(6):1-8. [CrossRef] [Medline]
  11. Ezeugwu V, Klaren RE, Hubbard EA, Manns PT, Motl RW. Mobility disability and the pattern of accelerometer-derived sedentary and physical activity behaviors in people with multiple sclerosis. Prev Med Rep. 2015;2:241-246. [CrossRef] [Medline]
  12. Cederberg KLJ, Mathison BG, Schuetz ML, Motl RW. Discrepancies between self-reported and device-measured sleep parameters in adults with multiple sclerosis. J Clin Sleep Med. Feb 1, 2022;18(2):415-421. [CrossRef] [Medline]
  13. Gulde P, Rieckmann P. The association between actigraphy-derived behavioral clusters and self-reported fatigue in persons with multiple sclerosis: cross-sectional study. JMIR Rehabil Assist Technol. Mar 17, 2022;9(1):e31164. [CrossRef] [Medline]
  14. Neal WN, Cederberg KL, Jeng B, Sasaki JE, Motl RW. Is symptomatic fatigue associated with physical activity and sedentary behaviors among persons with multiple sclerosis? Neurorehabil Neural Repair. Jun 2020;34(6):505-511. [CrossRef] [Medline]
  15. Jones CD, Jeng B, Silic P, Motl RW. Do device-measured physical activity and sedentary behavior differ by depression symptom status in persons with multiple sclerosis? Mult Scler Relat Disord. Jul 2022;63:103889. [CrossRef] [Medline]
  16. Macdonald E, Buchan D, Cerexhe L, Renfrew L, Sculthorpe N. Accelerometer measured physical activity and sedentary time in individuals with multiple sclerosis versus age matched controls: a systematic review and meta-analysis. Mult Scler Relat Disord. Jan 2023;69:104462. [CrossRef] [Medline]
  17. User Guide: ActiGraph GT9X Link + ActiLife, Revision 6. Amazon S3. Feb 18, 2020. URL: https:/​/s3.​amazonaws.com/​actigraphcorp.com/​wp-content/​uploads/​2020/​03/​05155628/​ActiGraph_Link_UserGuide_E.​200.​6001_Revision6_FINAL.​pdf [Accessed 2023-09-03]
  18. Choi L, Liu Z, Matthews CE, Buchowski MS. Validation of accelerometer wear and nonwear time classification algorithm. Med Sci Sports Exerc. Feb 2011;43(2):357-364. [CrossRef] [Medline]
  19. Di J, Leroux A, Urbanek J, et al. Patterns of sedentary and active time accumulation are associated with mortality in US adults: the NHANES study. bioRxiv. Preprint posted online on Aug 31, 2017. [CrossRef]
  20. Schrack JA, Kuo PL, Wanigatunga AA, et al. Active-to-sedentary behavior transitions, fatigability, and physical functioning in older adults. J Gerontol A Biol Sci Med Sci. Mar 14, 2019;74(4):560-567. [CrossRef] [Medline]
  21. Leroux A, Xu S, Kundu P, et al. Quantifying the predictive performance of objectively measured physical activity on mortality in the UK Biobank. J Gerontol A Biol Sci Med Sci. Jul 13, 2021;76(8):1486-1494. [CrossRef] [Medline]
  22. Welk GJ, Lamoureux NR, Zeng C, et al. Equating NHANES monitor-based physical activity to self-reported methods to enhance ongoing surveillance efforts. Med Sci Sports Exerc. Jun 1, 2023;55(6):1034-1043. [CrossRef] [Medline]
  23. Montoye AHK, Clevenger KA, Pfeiffer KA, et al. Development of cut-points for determining activity intensity from a wrist-worn ActiGraph accelerometer in free-living adults. J Sports Sci. Nov 2020;38(22):2569-2578. [CrossRef] [Medline]
  24. Rhudy MB, Dreisbach SB, Moran MD, Ruggiero MJ, Veerabhadrappa P. Cut points of the Actigraph GT9X for moderate and vigorous intensity physical activity at four different wear locations. J Sports Sci. Mar 2020;38(5):503-510. [CrossRef] [Medline]
  25. Xiao Q, Lu J, Zeitzer JM, Matthews CE, Saint-Maurice PF, Bauer C. Rest-activity profiles among U.S. adults in a nationally representative sample: a functional principal component analysis. Int J Behav Nutr Phys Act. Mar 24, 2022;19(1):32. [CrossRef] [Medline]
  26. Yeung CHC, Lu J, Soltero EG, Bauer C, Xiao Q. U.S. adolescent rest-activity patterns: insights from functional principal component analysis (NHANES 2011-2014). Int J Behav Nutr Phys Act. Oct 13, 2023;20(1):125. [CrossRef] [Medline]
  27. Stein MJ, Baurecht H, Sedlmeier AM, et al. Association between circadian physical activity patterns and mortality in the UK Biobank. Int J Behav Nutr Phys Act. Sep 1, 2023;20(1):102. [CrossRef] [Medline]
  28. Di J, Spira A, Bai J, et al. Joint and individual representation of domains of physical activity, sleep, and circadian rhythmicity. Stat Biosci. Jul 2019;11(2):371-402. [CrossRef] [Medline]
  29. Lock EF, Hoadley KA, Marron JS, Nobel AB. Joint and individual variation explained (JIVE) for integrated analysis of multiple data types. Ann Appl Stat. Mar 1, 2013;7(1):523-542. [CrossRef] [Medline]
  30. Salomon A, Galperin I, Buzaglo D, et al. Fragmentation, circadian amplitude, and fractal pattern of daily-living physical activity in people with multiple sclerosis: is there relevant information beyond the total amount of physical activity? Mult Scler Relat Disord. Dec 2022;68:104108. [CrossRef] [Medline]
  31. Ghosal R, Varma VR, Volfson D, et al. Distributional data analysis via quantile functions and its application to modeling digital biomarkers of gait in Alzheimer’s disease. Biostatistics. Jul 14, 2023;24(3):539-561. [CrossRef] [Medline]
  32. Ghosal R, Ghosh SK, Schrack JA, Zipunnikov V. Distributional outcome regression and its application to modelling continuously monitored heart rate and physical activity. arXiv. Preprint posted online on Jan 26, 2023. [CrossRef]
  33. Ziemssen T, Bhan V, Chataway J, et al. Secondary progressive multiple sclerosis: a review of clinical characteristics, definition, prognostic tools, and disease-modifying therapies. Neurol Neuroimmunol Neuroinflamm. Jan 2023;10(1):381-384. [CrossRef] [Medline]


ASTP: active-to-sedentary transition probability
CR: circadian rhythm
DARE: daytime activity ratio estimate
EDSS: Expanded Disability Status Scale
fPC: functional principal component
HEAL-MS: home-based evaluation of actigraphy to predict longitudinal function in multiple sclerosis
IPAQ: International Physical Activity Questionnaire
IS: interdaily stability
IV: intradaily variability
JIVE: Joint and Individual Variation Explained
L5: average log acceleration during the least active 5 hours of the day
LIPA: light intensity physical activity
M10: average log acceleration during the most active 10 hours of the day
MESOR: midline estimating statistic of rhythm
MS: multiple sclerosis
MVPA: moderate to vigorous physical activity
NHANES: National Health and Nutrition Examination Survey
PA: physical activity
PMS: progressive multiple sclerosis
RA: relative amplitude
REDCap: Research Electronic Data Capture
RRMS: relapsing-remitting multiple sclerosis
SATP: sedentary-to-active transition probability
TAC: total activity count
TLAC: total log-transformed activity count


Edited by Lorraine Buis; submitted 01.03.24; peer-reviewed by Chorong Park, Igor Vujovic; final revised version received 14.12.24; accepted 04.02.25; published 31.03.25.

Copyright

© Nicole Bou Rjeily, Muraleetharan Sanjayan, Pratim Guha Niyogi, Blake E Dewey, Alexandra Zambriczki Lee, Christy Hulett, Gabriella Dagher, Chen Hu, Rafal D Mazur, Elena M Kenney, Erin Brennan, Anna DuVal, Peter A Calabresi, Vadim Zipunnikov, Kathryn C Fitzgerald, Ellen M Mowry. Originally published in JMIR mHealth and uHealth (https://mhealth.jmir.org), 31.3.2025.

This is an open-access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work, first published in JMIR mHealth and uHealth, is properly cited. The complete bibliographic information, a link to the original publication on https://mhealth.jmir.org/, as well as this copyright and license information must be included.