Skip to main content

Accelerometer-derived movement features as predictive biomarkers for muscle atrophy in neurocritical care: a prospective cohort study

Abstract

Background

Physical inactivity and subsequent muscle atrophy are highly prevalent in neurocritical care and are recognized as key mechanisms underlying intensive care unit acquired weakness (ICUAW). The lack of quantifiable biomarkers for inactivity complicates the assessment of its relative importance compared to other conditions under the syndromic diagnosis of ICUAW. We hypothesize that active movement, as opposed to passive movement without active patient participation, can serve as a valid proxy for activity and may help predict muscle atrophy. To test this hypothesis, we utilized non-invasive, body-fixed accelerometers to compute measures of active movement and subsequently developed a machine learning model to predict muscle atrophy.

Methods

This study was conducted as a single-center, prospective, observational cohort study as part of the MINCE registry (metabolism and nutrition in neurointensive care, DRKS-ID: DRKS00031472). Atrophy of rectus femoris muscle (RFM) relative to baseline (day 0) was evaluated at days 3, 7 and 10 after intensive care unit (ICU) admission and served as the dependent variable in a generalized linear mixed model with Least Absolute Shrinkage and Selection Operator regularization and nested-cross validation.

Results

Out of 407 patients screened, 53 patients (age: 59.2 years (SD 15.9), 31 (58.5%) male) with a total of 91 available accelerometer datasets were enrolled. RFM thickness changed − 19.5% (SD 12.0) by day 10. Out of 12 demographic, clinical, nutritional and accelerometer-derived variables, baseline RFM muscle mass (beta − 5.1, 95% CI − 7.9 to − 3.8) and proportion of active movement (% activity) (beta 1.6, 95% CI 0.1 to 4.9) were selected as significant predictors of muscle atrophy. Including movement features into the prediction model substantially improved performance on an unseen test data set (including movement features: R2 = 79%; excluding movement features: R2 = 55%).

Conclusion

Active movement, as measured with thigh-fixed accelerometers, is a key risk factor for muscle atrophy in neurocritical care patients. Quantifiable biomarkers reflecting the level of activity can support more precise phenotyping of ICUAW and may direct tailored interventions to support activity in the ICU. Studies addressing the external validity of these findings beyond the neurointensive care unit are warranted.

Trial registration

DRKS00031472, retrospectively registered on 13.03.2023.

Background

Intensive care unit acquired weakness (ICUAW) describes a neuromuscular dysfunction secondary to critical illness and its treatment with consecutive generalized weakness. Data on prevalence for ICUAW show considerable variation due to diverse patient demographics and heterogenous methodology. However, with a systematic review pinpointing the median prevalence at 43% [1], its ubiquity in critical care is evident. Moreover, the impact resulting from ICUAW is profound and long-lasting, with patient outcomes significantly compromised for up to five years after discharge [2,3,4,5,6]. Therefore, ICUAW is acknowledged as a key component of post intensive care syndrome (PICS), highlighting its importance in the continuum of long-term recovery following critical care [7, 8].

ICUAW needs to be recognized as a clinical syndrome, rather than a specific disease entity. As such, it exhibits great heterogeneity and partially overlapping pathologies, which has diluted research findings and made the identification of treatable targets challenging in the past [9,10,11,12]. Relevant and common entities include critical illness myopathy (CIM), critical illness polyneuropathy (CIP) as well as critical illness polyneuromyopathy (CIPNM) as an overlap syndrome [9, 11, 12]. Electrophysiological methods including nerve conduction studies (NCS), electromyography and direct muscle stimulation have been successfully used to establish biomarkers for CIM, CIP and CIPMN [9, 13, 14]. Muscle atrophy due to mechanical unloading is also being recognized as a critical component of ICUAW. However, measurable biomarkers to assess the extent of inactivity of muscles are lacking.

In this regard, it is important to note that activity arises from active movement, as opposed to passive movement during mobilization without active patient participation. Hence, we postulate that establishing a proxy for activity can be achieved by applying non-invasive, body-fixed accelerometers to the lower extremities of critically ill patients while prospectively excluding episodes with passive mobilization such as intrahospital transports, physiotherapy and patient positioning. By introducing these biomarkers as continuous measures of active movement and incorporating these variables into a machine learning model, we aimed to predict rectus femoris muscle atrophy, as measured by ultrasound up to day 10 of intensive care unit (ICU) treatment. Based on the hypothesis that neurocritical care patients exhibit a higher prevalence of inactivity due to disorders of consciousness and motor deficits, we specifically included patients with acute brain injury in this trial.

Methods

Study design, setting and clinical management

This study was designed as a single-center, prospective, observational cohort study as part of the MINCE registry (metabolism and nutrition in neurointensive care, DRKS-ID: DRKS00031472, retrospectively registered on 13.03.2023) at a tertiary academic center (LMU University Hospital, Munich, Germany). Reporting follows the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) reporting guidelines. This study was approved by the local ethics committee (LMU Munich, project number 22-0173, 11.04.2022). Written consent was obtained from all participants or their next of kin. The study recruited from April 2022 to March 2024 and included patients within 48 h after ICU admission with age ≥ 18 years, neurologic disease as admitting diagnosis, and expected ICU length of stay ≥ 10 days. Patients with pre-existing neuromuscular disease, renal replacement therapy, pregnancy, pre-existing neoplastic disease, recent hospitalization (hospital stays longer than three days in the last three months, ICU treatment within the last three months), pre-existing confinement to bed, and pre-existing frailty (Clinical Frailty Scale > 3) were excluded.

Patients were mobilized at the treating physicians’ discretion. If indicated, patients received physiotherapy for 20–40 min/day on six days of the week and were repositioned and transferred from bed to chair regularly by the nursing staff. Nutritional therapy was conducted according to the European Society of Parenteral and Enteral Nutrition (ESPEN) guidelines, with caloric and protein targets of 25 kcal/kg/day and 1.3 g/kg/day, respectively [15]. As a reference, body weight as measured with bed scales was used for non-obese patients, and ideal body weight was used for patients with a body mass index (BMI) > 30 kg/m2 [15]. During the acute phase of illness (days 1–3), hypocaloric nutrition (70% of energy expenditure (EE)) was aimed for. From day 4 on, isocaloric (100% of EE) nutrition was implemented.

Data collection

Clinical data prospectively collected on the ICU included age, sex, body mass index (BMI), admission diagnosis, cumulative protein and calorie deficit, duration of mechanical ventilation, ICU length of stay (LOS), daily Sepsis-related Organ Failure Assessment score (SOFA) and SOFA without Glasgow Coma Scale (GCS) score (mSOFA), Acute Physiology and Chronic Health Evaluation (APACHE II) score on ICU admission, Nutrition Risk in Critically ill score (NUTRIC) on admission, premorbid modified Rankin Scale (pmRS) and Glasgow Outcome Scale Extended (GOSE) at ICU discharge.

Ultrasound of the upper thigh (rectus femoris muscle, RFM) and temporalis muscle (TM) was performed bilaterally using a 20 MHz linear probe (MyLabOmega, Esaote, Genoa, Italy) upon admission, and on days 3, 7 and 10. As previously described [16, 17], the site of measurement for RFM was marked in the lower third of the connecting line between the anterior superior iliac spine and the upper edge of the patella with a permanent marker to ensure reproducibility between measurements (Supplementary Fig. 1). Measurements were conducted according to a local protocol that emphasized minimal compression during RFM sonography and called for individual adjustments of depth and gain to optimally visualize the surface of the femur and to delineate fascial borders, respectively. Measurement of TM followed the protocol as described by Maskos et al. [17] Three repeated measurements were performed by one of six raters (TP, LG, LR, AM, JB, SI), using the built-in software of the ultrasound machine to measure muscle thickness. The mean value of the repeated measurements was used for further analysis. Reliability of repeated ultrasound measurements is reported in Supplementary Table 1.

Tri-axial accelerometers (range ± 16 g; sampling rate 12.5 Hz; Axivity Ltd., Newcastle upon Tyne, UK) were attached within 48 h of admission to both upper thighs with transparent adhesive tape. Skin inspections and minor adjustments to the sensor placement were performed every third day to prevent any pressure damage. To exclude any passive movement not contributing to the patient’s activity, episodes with physiotherapy, intrahospital transports, or repositioning by nursing staff were prospectively documented and excluded from the recorded data. The sensor data was extracted and analyzed by an author (MW) not involved in the patients’ clinical management and blinded for the ultrasound measurements.

As a positive control, accelerometers were attached to healthy individuals (n = 3). Static placement of sensors served as a negative control (Supplementary Table 2).

Accelerometer data processing

OmGUI version 1.0.0.11 software was used to download the raw acceleration data from the devices. MATLAB (Mathworks Inc., Natick, USA, release R2022a, version 9.12.0) was used for data processing. Periods of active movement were identified following a previously established procedure that has been already used for activity recognition in ICU settings [18, 19]. Accordingly, recorded time series from each axial component (x, y, z) were first down-sampled to 10 Hz and subsequently high-pass filtered (4th order Butterworth filter, cutoff frequency at 0.2 Hz) to remove baseline offset and low-frequency effects reflecting static postural orientation. Filtered time series were segmented into non-overlapping 5 s windows for subsequent motion feature extraction. Signal magnitude area (SMA) was then computed for every window [19] to identify activity bouts (AB) using a defined threshold of SMA ≥ 0.135 g [18]. Across all identified AB, the mean intensity (AB-intensity) and duration (AB-duration) as well as the variability (standard deviation, SD) of these features were calculated. The distribution of motion features across ABs is log-normal, which required estimating mean and SD via a maximum likelihood technique [20]. The overall movement intensity, the proportion of active movement (%active) and ABs per hour were calculated based on the entire duration of the recording (Fig. 1).

Fig. 1
figure 1

Accelerometer-derived features. Body motion was monitored using tri-axial accelerometers bilaterally attached to the upper thigh (1). Raw triaxial accelerometer recordings were first offset eliminated, and the time series were segmented into non-overlapping 5 s windows. The signal magnitude area was computed for every window (2). Bouts of dynamic activity were identified based on the threshold ≥ 0.135 g (3) and a set of motion features was computed for every bout of activity (4). Finally, the average and distribution of motion features across all bouts of activity were computed (5). acc = acceleration

Predictive modeling and statistical analysis

As the dataset includes multiple observations (both legs) per patient and exhibits linearity as evaluated by exploratory data analysis, a generalized linear mixed model (GLMM) to account for intra-patient correlation by using individual patients as a random effect was chosen. Given the numerous independent variables of interest, including demographic, clinical, and activity-related features, a rigorous approach to model selection and validation to prevent overfitting was required. Therefore, we employed regularization with Least Absolute Shrinkage and Selection Operator (LASSO), which penalizes the GLMM model via L1-norm and in effect shrinks the weight of non-contributing features to zero.

First, multicollinearity among predictors was mitigated by excluding variables with a variance inflation factor (VIF) exceeding 5 (removing AB per hour and SOFA) [21]. Next, standardization (z-score normalization) of the remaining prediction variables (age, sex, baseline RFM muscle mass, mSOFA, calorie deficit, protein deficit, overall intensity, %active, AB-intensitylog-mean, AB-durationlog-mean, AB-intensitylog-SD, AB-durationlog-SD) was performed to ensure equal weights and comparable units. To allow testing on unseen data, a stratified split was executed to divide the data into training (80%) and test (20%) sets. The training set was further used for optimizing the hyperparameter of GLMM-LASSO using a machine learning approach with nested cross-validation (Fig. 2) [22]. Model performance was evaluated on the test set using mean squared error, root mean squared error, mean absolute error, R-squared (squared correlation method, R2) and a plot depicting actual versus predicted values.

Fig. 2
figure 2

Nested-cross validation of a regularized GLMM model. After standardization and a stratified 80/20 split, the training data set was partitioned into 4 folds (outer loop). Within each outer fold, an inner loop of 2 folds was used for hyperparameter tuning. The hyperparameter (lambda) that minimized the mean squared error in the inner loop was selected. The model with this optimal lambda was then evaluated on the validation fold of the outer loop. This process was repeated for all 4 outer folds, resulting in an optimal lambda for each fold. The final model was chosen using the average of the optimal lambdas from all outer folds. Finally, the performance of this final model was assessed using the unseen test set. GLMM = generalized mixed effects model; lasso = least absolute shrinkage and selection operator

To illustrate the level of uncertainty of the model coefficients, bootstrapping with 10,000 resamples was performed to estimate the bias-corrected and accelerated (BCa) confidence intervals for the model coefficients.

Additionally, and to test the relevance of leg movement on TM as a muscle group unaffected by thigh movement, we used a linear regression model for the prediction of TM atrophy at day 10 including all demographic, clinical and nutritional variables as well as %active as independent variables. To compare %active between healthy individuals and the neurointensive care unit (NICU) cohort, a two-tailed t-test was performed. Further, we identified patients with unilateral upper motor neuron damage and corresponding motor deficits to investigate the contribution of upper motor neuron lesion on muscle atrophy. To compare the magnitude of muscle atrophy and to account for within-subject correlation, we used Generalized Estimating Equations (GEE) with post hoc pairwise comparisons and Bonferroni adjustment.

Summary statistics for continuous variables are presented as means with standard deviation (SD) for normally distributed data and as medians with interquartile ranges (IQR) for non-normally distributed data, with normality assessed using Quantile–Quantile plots and Shapiro–Wilk test. Categorical variables are summarized as frequencies and percentages.

All analyses were performed using R (2023.06.1 + 524) using ‘stats’, ‘psych’, ‘ggplot2’, ‘dplyr’, ‘lme4’, ‘nlme’, ‘geepack’, ‘multcomp’, ‘emmeans’, ‘MASS’, ‘svglite’ ‘glmmLasso’,‘caret’ and ‘boot’, packages. ChatGPT (version 4) was used for error handling, repetitive programming, and overall optimization of code in R.

Results

Patient characteristics

Out of 407 patients screened, 53 with a total of 91 available accelerometer datasets were enrolled in this study (Fig. 3). Clinical baseline characteristics are presented in Table 1. Of all patients included in the analysis, mean age was 59.2 years (SD 15.9) and 31 (58.5%) were male. Cerebrovascular diseases were the most frequent ICU admission diagnoses (86.8%, 46/53). Mean ICU length of stay was 17.0 days (IQR 8.0), while the mean duration of mechanical ventilation was 15.6 days (SD 9.2). During the observation period, patients met 62.6% (SD 18.4) of the caloric goals and 57.9% (SD 21.6) of protein goals according to the ESPEN guidelines.

Fig. 3
figure 3

Screening and study inclusion. ICU = Intensive Care Unit;

Table 1 Baseline characteristics

Active movement and muscle atrophy during ICU treatment

Muscular atrophy as measured with ultrasound was more pronounced in RFM compared to TM (-19.5% (SD 12.0) versus -15.3% (SD 11.1) at day 10) (Fig. 4A). Active movement of NICU patients as indicated by proportion of active movement (%active) over time is infrequent, particular at early stages of the ICU stay. While mean %active stays low over the entire time, some patients exhibit higher activity starting around day 3. (Fig. 4B). Compared to healthy individuals, NICU patients demonstrate a significant reduction in active movement (%active: healthy individuals 13.3 (SD 0.8) vs. NICU patients 0.84 (SD 1.08), p < 0.001) (Supplementary Tables 2 and 4). No adverse events were observed in association with the placement of accelerometers in the ICU setting (Supplementary Table 3).

Fig. 4
figure 4

Active movement and muscle atrophy during ICU treatment. Muscle atrophy at days 3, 5 7 and 10 relative to day 0 for RFM and TM (A). Proportion of active movement (%active) over time (B). ICU = Intensive Care Unit;

Predictive models for muscle atrophy

The machine learning model based on a total of 12 demographic, clinical, nutritional and accelerometer-derived variables selected baseline RFM muscle mass (beta − 5.1, 95% confidence interval (95% CI) − 7.9 to − 3.8) and %active (beta 1.6, 95% CI 0.1 to 4.9) to explain 79% (R2 = 79%) of the occurring variance in muscle wasting in an unseen test data set (Fig. 5, A and B). Thus, for every standard deviation increase in baseline RFM (2.6 mm), RFM thickness at day 10 is estimated to decrease by another 5.1 percentage points (49.5% relative change). In contrast, a standard unit increase in %active (1.1%) is projected to result in 1.6 percentage points less RFM atrophy at day 10 (relative change 15.5%). Ignoring movement features as predictors for RFM muscle atrophy results in substantially worse model performance (R2 = 55%) (Fig. 5C, Supplementary Table 6). RFM atrophy was not significantly different between immobile limbs with upper motor neuron lesions (UMNL) and immobile limbs without UMNL. However, muscle atrophy was markedly decreased in limbs with active movement (Supplementary Fig. 2). Thigh-fixed accelerometer data did not contribute significantly to a model predicting TM atrophy (Supplementary Table 5).

Fig. 5
figure 5

Prediction of MRF muscle atrophy with and without movement features. Standardized coefficients and 95% confidence intervals (asterisks indicate significant predictors) of the regularized regression models with (modelmovement+, A) and without movement features (modelmovement−, C). Out of all demographic (age, sex), clinical (baseline RFM muscle mass, mSOFA), nutritional (calorie deficit, protein deficit) and movement variables (intensity, %active, AB-intensitylog-mean, AB-durationlog-mean, AB-intensitylog-SD, AB-durationlog-SD), the depicted 10/12 independent variables for model modelmovement+ and 4/6 independent variables for modelmovement− were selected for the final models, respectively. Significant predictors in modelmovement+ included baseline RFM muscle mass (beta − 5.1, 95% confidence interval (95% CI) − 7.9 to − 3.8) and %active (beta 1.6, 95% CI 0.1 to 4.9). For modelmovement-, only baseline RFM muscle was found as a statistically significant predictor (beta − 4.6, 95% CI − 7.6 to − 3.9). Scatter plots with regression line of predicted versus actual muscle wasting (grey dots: training data; black dots: unseen test data) for modelmovement+ (B) and modelmovement- (D), respectively (R2: 0.79 vs. 0.55, RMSE: vs. 8.4 vs. 10.7 mm; MAE: 6.2 vs. 8.0 mm). mSOFA = SOFA without GCS;

Discussion

In this prospective cohort study, we used thigh-fixed accelerometers to establish movement features as predictive biomarkers for muscle atrophy in neurocritical care patients. Proportion of active movement (%active) demonstrated a significant protective effect against muscle wasting and improved the precision of muscle atrophy prediction in an unseen test data set. To the best of our knowledge, this is the first quantifiable and validated measure that provides information on the relative importance of inactivity for muscle atrophy in critically ill patients.

It is crucial to distinguish between immobility and inactivity, especially within the ICU context, as inactivity can occur despite mobilization efforts due to a lack of active patient participation (passive movement). To address this, we excluded periods such as physiotherapy, intrahospital transports, and patient positioning, from our analysis. Therefore, we consider the movement features as surrogates for activity rather than measures of mobility. Importantly, and as a limitation of this approach, movement sensors are unable to capture any muscle activity without movement via isometric contractions (active immobility).

The relevance of inactivity, as compared to immobility, as the variable of interest in this context is further exemplified by clinical trials on electrical muscle stimulation (EMS) and interventions focusing on early (passive) mobilization. The current evidence highlights the efficacy of EMS [23, 24], whereas mobilization trials demonstrated limited efficacy and raised safety concerns [25,26,27,28,29,30,31]. While the latter often involve passive mobilization without genuine patient activity, EMS generates muscle activity without requiring mobility. Considering the data, a reasonable strategy to prevent muscle atrophy in critically ill patients may involve first measuring the extent of active movement with accelerometers to identify those at risk, and subsequently promoting activity (with or without mobilization) based on the patient's stability.

The pathophysiology of mechanical unloading leading to atrophy has so far only been systematically studied and quantified outside the ICU. Studies with cast immobilization of the lower extremity for two weeks in healthy adults and examination of astronauts after 8 days of space flight revealed a 5% and 6% decrease in quadriceps muscle mass, respectively [32, 33]. A recent meta-analysis analyzing the general ICU population estimated muscle atrophy to be around 16% at day 10 for RFM [34]. In comparison, our data demonstrate a more pronounced rate of RFM atrophy, showing a 19.5% decrease by day 10. While additional factors such as CIP, CIM, and CIPNM certainly contribute to the higher rate of atrophy in ICU patients, the residual activity in cast immobilization (via isometric contractions) and during space flight (via active movement with reduced muscle activity) may also account for the observed differences.

Given that no or passive movement was described in more than 70% of patients in the first 48 h, and still more than 40% after two weeks, in the TEAM trial [26], it is plausible to assume that inactivity also significantly contributes to ICUAW in the general ICU population. Yet, as disorders of consciousness and focal-neurological deficits are major barriers to mobilization and activity [26], this might be even more relevant for neurointensive care patients. Although our accelerometer data are not directly comparable to the ICU mobility scale used in the TEAM trial, it indicates extremely infrequent periods of active movement for most patients over a 10-day observation period, reaching only 6% (0.84/13.3) of the activity level of healthy individuals. These numbers are parallelled by data from González-Seguel et al., who found mechanically ventilated patients to be inactive during the ICU stay in over 96% of the time [35]. This, coupled with the prominence of movement features as predictors of muscle atrophy in our prospective cohort, further strengthens the significance of inactivity in (neuro-) critical care. Other studies within the ICU have investigated accelerometry primarily in the context of sleep, circadian rhythm, and sedation levels. However, these studies exhibit limitations, such as narrow observation periods and the absence of well-defined thresholds for activity measurement [36,37,38,39].

Accelerometer-derived data have also been validated as biomarkers for muscle atrophy outside the ICU setting. In a study with almost 500 elderly participants, Sanchez-Sanchez et al. investigated the association of physical activity as measured with hip-worn accelerometers and sarcopenia. Here, higher physical activity correlated with better performance in sarcopenia-related scores [40]. Similarly, Foong et al. showed a positive association of accelerometer-derived physical activity with muscle mass and muscle strength [41].

The exercise stimulus, as the ultimate determinant for activity, can be delineated into two primary variables: volume and intensity. In exercise physiology, the volume of exercise is traditionally quantified by the number of repetitions performed, while intensity is commonly measured by the force exerted during exercise [42]. In our ICU cohort, we utilized proportion of active movement (%active) and AB-durationlog-mean as proxies for exercise volume. For critically ill patients, the force generated cannot be measured pragmatically. We therefore introduced movement intensity (resultant acceleration magnitude) as a surrogate of exercise intensity. The LASSO regularization used to address the high number and potential co-linearity of parameters revealed %active as an approximation of exercise volume as a relevant predictor, while surrogates of intensity were not selected. Thus, intensity may either be irrelevant considering the uniform force generated by patients moving against gravity and not against resistance, or movement intensity is not a valid biomarker of the exercise intensity.

Besides proportion of active movement, baseline muscle mass was predictive of muscle atrophy. This finding is in line with studies in healthy participants. Here, higher age with lower baseline muscle mass showed significantly less pronounced atrophy. However, older participants with lower baseline muscle mass suffered from greater loss of muscle strength after immobilization [43,44,45]. Furthermore, variations in atrophy can be observed across muscle and fiber types. Anti-gravity muscles with high proportions of slow type 1 fibers, such as RFM, seem to exhibit a selective vulnerability [45], which is in line with our data demonstrating pronounced atrophy of RFM over TM.

The strengths of our study include the selection of a neurointensive care population devoid of frailty, specifically targeting individuals at high risk of muscle atrophy without pre-existing muscle wasting. However, the study's focus on neurocritical care may limit its generalizability and further research is needed to confirm the applicability of our results to more diverse patient populations. The extent of neurogenic atrophy mediated via various damage to the lower motor neuron was not explicitly measured in our study but can be assumed to be minimal given the demographics of our cohort. As that general ICU cohorts also experience lower motor neuron lesions due to critical illness or its treatments, this confounder is likely to be similar across groups. For UMNL, we do not expect neurogenic atrophy, as it does not result in direct muscle denervation. Supporting this pathophysiological hypothesis, we could not observe any difference in the extent of muscle atrophy between immobile limbs with upper motor neuron lesions (UMNL) and immobile limbs without UMNL. However, muscle atrophy was markedly decreased in limbs with active movement, suggesting no role for UMNL as a mediator for atrophy. We ensured high data quality by filtering out passive mobilization and prospectively collecting clinical data. Furthermore, our analysis is underpinned by a strong statistical framework leveraging machine learning to identify the most important predictors and using unseen data to validate these findings. Further limitations of our study are primarily rooted in the fact that muscle morphology does not necessarily equate to function. We used muscle ultrasound, a widely adopted and validated surrogate for ICUAW [34, 46, 47], instead of the Medical Research Council Sum Score (MRC-SS), as the latter is often deemed infeasible in the general ICU cohort, and even more so in neurocritical care. Additionally, we decided against including measures of the upper extremities because muscle volume is challenging to determine via ultrasound due to variability in muscle thickness relative to positioning, difficult anatomical landmarks, and less pronounced atrophy compared to the lower extremities. Instead, we focused on the RFM as the established sonographic gold standard, along with the TM as a muscle group unrelated to the movement captured by the accelerometers.

Conclusion

Active movement, as a surrogate of muscle activity, can be quantified using non-invasive, thigh-fixed accelerometers and adds value for the prediction of muscle atrophy in neurocritical care patients. Establishing movement-derived biomarkers enables better phenotyping of ICUAW, thereby providing a foundation for tailored interventions and should be included as covariates in trials on ICUAW in the future. Studies addressing the external validity of these findings beyond the neurointensive care unit are warranted.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request. Matlab code for analyzing motion features is available at https://github.com/DSGZ-MotionLab/ICU_motion_analysis/.

Abbreviations

AB:

Activity bout

APACHE II:

Acute physiology and chronic health evaluation II

BCa:

Bias-corrected and accelerated

BMI:

Body mass index

CI:

Confidence interval

CIM:

Critical illness myopathy

CIP:

Critical illness polyneuropathy

CIPNM:

Critical illness polyneuromyopathy

EE:

Energy expenditure

EMS:

Electrical muscle stimulation

ESPEN:

European Society of Parenteral and Enteral Nutrition

GCS:

Glasgow Coma Scale

GEE:

Generalized Estimating Equations

GLMM:

Generalized linear mixed model

GOSE:

Glasgow Outcome Scale—Extended

HLF:

High/low pass filter

ICU:

Intensive care unit

ICUAW:

Intensive care unit acquired weakness

LASSO:

Least Absolute Shrinkage and Selection Operator

LOS:

Length of stay

MINCE:

Metabolism and nutrition in neurointensive care

MRC-SS:

Medical Research Council Sum Score

mSOFA:

Modified SOFA

NCS:

Nerve conduction studies

NICU:

Neurointensive care unit

NUTRIC:

Nutrition risk in the critically ill

PICS:

Post intensive care syndrome

pmRS:

Premorbid modified Rankin scale

QQ:

Quantile–Quantile

RFM:

Rectus femoris muscle

SMA:

Signal magnitude area

SOFA:

Sequential organ failure assessment

TM:

Temporalis muscle

UMNL:

Upper motor neuron lesion

References

  1. Fan E, Cheek F, Chlan L, Gosselink R, Hart N, Herridge MS, et al. An official American Thoracic Society Clinical Practice guideline: the diagnosis of intensive care unit–acquired weakness in adults. Am J Respir Crit Care Med. 2014;190:1437–46.

    Article  PubMed  Google Scholar 

  2. Van Aerde N, Meersseman P, Debaveye Y, Wilmer A, Gunst J, Casaer MP, et al. Five-year impact of ICU-acquired neuromuscular complications: a prospective, observational study. Intensive Care Med. 2020;46:1184–93.

    Article  PubMed  Google Scholar 

  3. Kelmenson DA, Held N, Allen RR, Quan D, Burnham EL, Clark BJ, et al. Outcomes of ICU patients with a discharge diagnosis of critical illness polyneuromyopathy: a propensity-matched analysis. Crit Care Med. 2017;45:2055–60.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Saccheri C, Morawiec E, Delemazure J, Mayaux J, Dubé B-P, Similowski T, et al. ICU-acquired weakness, diaphragm dysfunction and long-term outcomes of critically ill patients. Ann Intensive Care. 2020;10:1. https://doi.org/10.1186/s13613-019-0618-4.

    Article  PubMed  PubMed Central  Google Scholar 

  5. De Jonghe B, Sharshar T, Lefaucheur J-P, Authier F-J, Durand-Zaleski I, Boussarsar M, et al. Paresis acquired in the intensive care unit: a prospective multicenter study. JAMA. 2002;288:2859–67.

    Article  PubMed  Google Scholar 

  6. Hermans G, Van Mechelen H, Bruyninckx F, Vanhullebusch T, Clerckx B, Meersseman P, et al. Predictive value for weakness and 1-year mortality of screening electrophysiology tests in the ICU. Intensive Care Med. 2015;41:2138–48.

    Article  PubMed  Google Scholar 

  7. Needham DM, Davidson J, Cohen H, Hopkins RO, Weinert C, Wunsch H, et al. Improving long-term outcomes after discharge from intensive care unit. Crit Care Med. 2012;40:502–9.

    Article  PubMed  Google Scholar 

  8. Hermans G, Van den Berghe G. Clinical review: intensive care unit acquired weakness. Crit Care. 2015;19:274.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Vanhorebeek I, Latronico N, Van den Berghe G. ICU-acquired weakness. Intensive Care Med. 2020;46:637–53. https://doi.org/10.1007/s00134-020-05944-4.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Friedrich O, Reid MB, Van den Berghe G, Vanhorebeek I, Hermans G, Rich MM, et al. The sick and the weak: neuropathies/myopathies in the critically Ill. Physiol Rev. 2015;95:1025–109.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Latronico N, Rasulo FA, Eikermann M, Piva S. Illness weakness, polyneuropathy and myopathy: diagnosis, treatment, and long-term outcomes. Crit Care. 2023;27:439.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Piva S, Fagoni N, Latronico N. Intensive care unit–acquired weakness: unanswered questions and targets for future research. F1000Res. 2019;8:508.

    Article  Google Scholar 

  13. Latronico N, Bertolini G, Guarneri B, Botteri M, Peli E, Andreoletti S, et al. Simplified electrophysiological evaluation of peripheral nerves in critically ill patients: the Italian multi-centre CRIMYNE study. Crit Care. 2007;11:R11.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Kelmenson DA, Quan D, Moss M. What is the diagnostic accuracy of single nerve conduction studies and muscle ultrasound to identify critical illness polyneuromyopathy: a prospective cohort study. Crit Care. 2018;22:1–9.

    Article  Google Scholar 

  15. Singer P, Blaser AR, Berger MM, Alhazzani W, Calder PC, Casaer MP, et al. ESPEN guideline on clinical nutrition in the intensive care unit. Clin Nutr. 2019;38:48–79.

    Article  PubMed  Google Scholar 

  16. Pardo E, El Behi H, Boizeau P, Verdonk F, Alberti C, Lescot T. Reliability of ultrasound measurements of quadriceps muscle thickness in critically ill patients. BMC Anesthesiol. 2018;18:205.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Maskos A, Schmidbauer ML, Kunst S, Rehms R, Putz T, Römer S, et al. Diagnostic utility of temporal muscle thickness as a monitoring tool for muscle wasting in neurocritical care. Nutrients. 2022;14:4498.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Lugade V, Fortune E, Morrow M, Kaufman K. Validity of using tri-axial accelerometers to measure human movement—Part I: posture and movement detection. Med Eng Phys. 2014;36:169–76.

    Article  PubMed  Google Scholar 

  19. Bhattacharyay S, Rattray J, Wang M, Dziedzic PH, Calvillo E, Kim HB, et al. Decoding accelerometry for classification and prediction of critically ill patients with severe brain injury. Sci Rep. 2021;11:23654.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Rochester L, Chastin SFM, Lord S, Baker K, Burn DJ. Understanding the impact of deep brain stimulation on ambulatory activity in advanced Parkinson’s disease. J Neurol. 2012;259:1081–6.

    Article  PubMed  Google Scholar 

  21. James G, Witten D, Hastie T, Tibshirani R. An Introduction to Statistical Learning. 1st ed. New York: Springer; 2013. https://doi.org/10.1007/978-1-4614-7138-7.

    Book  Google Scholar 

  22. Raschka S. Model evaluation, model selection, and algorithm selection in machine learning. arXiv:1811.12808. 2018.

  23. Nakanishi N, Yoshihiro S, Kawamura Y, Aikawa G, Shida H, Shimizu M, et al. Effect of neuromuscular electrical stimulation in patients with critical illness: an updated systematic review and meta-analysis of randomized controlled trials. Crit Care Med. 2023. https://doi.org/10.1097/CCM.0000000000005941.

    Article  PubMed  Google Scholar 

  24. Valenzuela PL, Morales JS, Pareja-Galeano H, Izquierdo M, Emanuele E, de la Villa P, et al. Physical strategies to prevent disuse-induced functional decline in the elderly. Ageing Res Rev. 2018;47:80–8.

    Article  PubMed  Google Scholar 

  25. Kho ME, Berney S, Pastva AM, Kelly L, Reid JC, Burns KEA, et al. Early in-bed cycle ergometry in mechanically ventilated patients. NEJM Evid. 2024. https://doi.org/10.1056/EVIDoa2400137.

    Article  PubMed  Google Scholar 

  26. The TEAM Study Investigators and the ANZICS Clinical Trials Group. Early Active Mobilization during Mechanical Ventilation in the ICU. New Engl J Med. 2022;387:1747–58. https://doi.org/10.1056/NEJMoa2209083

  27. Schaller SJ, Anstey M, Blobner M, Edrich T, Grabitz SD, Gradwohl-Matis I, et al. Early, goal-directed mobilisation in the surgical intensive care unit: a randomised controlled trial. The Lancet. 2016;388:1377–88.

    Article  Google Scholar 

  28. Wright SE, Thomas K, Watson G, Baker C, Bryant A, Chadwick TJ, et al. Intensive versus standard physical rehabilitation therapy in the critically ill (EPICC): a multicentre, parallel-group, randomised controlled trial. Thorax. 2018;73:213–21.

    Article  PubMed  Google Scholar 

  29. Moss M, Nordon-Craft A, Malone D, Van Pelt D, Frankel SK, Warner ML, et al. A randomized trial of an intensive physical therapy program for patients with acute respiratory failure. Am J Respir Crit Care Med. 2016;193:1101–10.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Schweickert WD, Pohlman MC, Pohlman AS, Nigos C, Pawlik AJ, Esbrook CL, et al. Early physical and occupational therapy in mechanically ventilated, critically ill patients: a randomised controlled trial. The Lancet. 2009;373:1874–82.

    Article  Google Scholar 

  31. Morris PE, Berry MJ, Files DC, Thompson JC, Hauser J, Flores L, et al. Standardized rehabilitation and hospital length of stay among patients with acute respiratory failure a randomized clinical trial. JAMA J Am Med Assoc. 2016;315:2694–702.

    Article  Google Scholar 

  32. Jones SW, Hill RJ, Krasney PA, O’Conner B, Peirce N, Greenhaff PL. Disuse atrophy and exercise rehabilitation in humans profoundly affects the expression of genes associated with the regulation of skeletal muscle mass. FASEB J. 2004;18:1025–7.

    Article  PubMed  Google Scholar 

  33. LeBlanc A, Rowe R, Schneider V, Evans H, Hedrick T. Regional muscle loss after short duration spaceflight. Aviat Space Environ Med. 1995;66:1151–4.

    PubMed  Google Scholar 

  34. Fazzini B, Märkl T, Costas C, Blobner M, Schaller SJ, Prowle J, et al. The rate and assessment of muscle wasting during critical illness: a systematic review and meta-analysis. Crit Care. 2023;27:2. https://doi.org/10.1186/s13054-022-04253-0.

    Article  PubMed  PubMed Central  Google Scholar 

  35. González-Seguel F, Camus-Molina A, Leiva-Corvalán M, Mayer KP, Leppe J. Uninterrupted actigraphy recording to quantify physical activity and sedentary behaviors in mechanically ventilated adults. J Acute Care Phys Ther. 2022;13:190–7.

    Article  Google Scholar 

  36. Fazio S, Doroy A, Da Marto N, Taylor S, Anderson N, Young HM, et al. Quantifying mobility in the ICU: comparison of electronic health record documentation and accelerometer-based sensors to clinician-annotated video. Crit Care Explor. 2020;2:E0091.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Gupta P, Martin JL, Needham DM, Vangala S, Colantuoni E, Kamdar BB. Use of actigraphy to characterize inactivity and activity in patients in a medical ICU. Heart Lung. 2020;49:398–406.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Kamdar BB, Kadden DJ, Vangala S, Elashoff DA, Ong MK, Martin JL, et al. Feasibility of continuous actigraphy in patients in a medical intensive care unit. Am J Crit Care. 2017;26:329–35.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Verceles AC, Hager ER. Use of accelerometry to monitor physical activity in critically ill subjects: a systematic review. Respir Care. 2015;60:1330–6.

    Article  PubMed  Google Scholar 

  40. Sánchez-Sánchez JL, Mañas A, García-García FJ, Ara I, Carnicero JA, Walter S, et al. Sedentary behaviour, physical activity, and sarcopenia among older adults in the TSHA: isotemporal substitution model. J Cachexia Sarcopenia Muscle. 2019;10:188–98.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Foong YC, Chherawala N, Aitken D, Scott D, Winzenberg T, Jones G. Accelerometer-determined physical activity, muscle mass, and leg strength in community-dwelling older adults. J Cachexia Sarcopenia Muscle. 2016;7:275–83.

    Article  PubMed  Google Scholar 

  42. Marston KJ, Peiffer JJ, Newton MJ, Scott BR. A comparison of traditional and novel metrics to quantify resistance training. Sci Rep. 2017;7:5606.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Suetta C, Hvid LG, Justesen L, Christensen U, Neergaard K, Simonsen L, et al. Effects of aging on human skeletal muscle after immobilization and retraining. J Appl Physiol. 2009;107:1172–80.

    Article  PubMed  Google Scholar 

  44. Hvid L, Aagaard P, Justesen L, Bayer ML, Andersen JL, Ørtenblad N, et al. Effects of aging on muscle mechanical function and muscle fiber morphology during short-term immobilization and subsequent retraining. J Appl Physiol. 2010;109:1628–34.

    Article  PubMed  Google Scholar 

  45. Bodine SC. Disuse-induced muscle wasting. Int J Biochem Cell Biol. 2013;45:2200–8.

    Article  PubMed  Google Scholar 

  46. Zhang W, Wu J, Gu Q, Gu Y, Zhao Y, Ge X, et al. Changes in muscle ultrasound for the diagnosis of intensive care unit acquired weakness in critically ill patients. Sci Rep. 2021;11:18280.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Mourtzakis M, Parry S, Connolly B, Puthucheary Z. Skeletal muscle ultrasound in critical care: a tool in need of translation. Ann Am Thorac Soc. 2017;14:1495–503.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

The authors thank the staff of the neurological and neurosurgical intensive care units at LMU Munich for their great interest and cooperation.

Funding

Open Access funding enabled and organized by Projekt DEAL. MLS was supported by the Deutsche Forschungsgemeinschaft (TRR 274) and Stiftungen zugunsten der Medizinischen Fakultät (LMU Munich).

Author information

Authors and Affiliations

Authors

Contributions

Conceptualization, methodology, validation: MLS, MW, KD. Data curation, formal analysis, visualization, software: MLS, MW. Investigation: TP, LG, LR, AM. Writing—original draft: MLS, TP, MW. Writing—review and editing: LG, LR, AM, JZ, TW, KD. Resources: MW, TW, KD. Project administration, supervision, funding acquisition: MLS, KD. All authors reviewed the manuscript.

Corresponding author

Correspondence to Moritz L. Schmidbauer.

Ethics declarations

Ethics approval and consent to participate

Approval for this study was granted by the local ethics committee (LMU Munich, approval number: 22-0173, data of approval: 11.04.2022). This study was conducted in accordance with the Declaration of Helsinki and ethical guidelines for medical and health research involving human subjects. Written informed consent was obtained from all patients or next of kin.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

13054_2024_5067_MOESM1_ESM.docx

Supplementary Figure 1. Accelerometer placement. Figure 2. Influence of movement and upper motor neuron lesion on muscle atrophy. Table 1. Reliability of RFM ultrasound measurements. Table 2. Negative and positive controls. Table 3. Adverse events of accelerometer placement. Table 4. Accelerometer data. Table 5. Linear regression model for TM atrophy. Table 6. Performance metrics of models with and without movement features.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Schmidbauer, M.L., Putz, T., Gehri, L. et al. Accelerometer-derived movement features as predictive biomarkers for muscle atrophy in neurocritical care: a prospective cohort study. Crit Care 28, 288 (2024). https://doi.org/10.1186/s13054-024-05067-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13054-024-05067-y

Keywords