Source document: https://docs-lawep.s3.us-east-2.amazonaws.com/1698059115121.pdf (Liberia)
OPEN Soft tissue radiodensity parameters mediate the relationship between self‑reported physical activity and lower extremity function in AGES‑Reykjavík participants Kyle J. Edmunds1,2,3*, Ozioma C. Okonkwo1,2,3, Sigurdur Sigurdsson5, Sarah R. Lose1,2,3, Vilmundur Gudnason5,6, Ugo Carraro7 & Paolo Gargiulo4,8 Although previous studies have highlighted the association between physical activity and lower extremity function (LEF) in elderly individuals, the mechanisms underlying this relationship remain debated.
Our recent work has recognized the utility of nonlinear trimodal regression analysis (NTRA) parameters in characterizing changes in soft tissue radiodensity as a quantitative construct for sarcopenia in the longitudinal, population‑based cohort of the AGES‑Reykjavík study.
For the present work, we assembled a series of prospective multivariate regression models to interrogate whether NTRA parameters mediate the 5‑year longitudinal relationship between physical activity and LEF in AGES‑Reykjavík participants.

Healthy elderly volunteers from the AGES‑Reykjavík cohort underwent mid‑thigh X‑ray CT scans along with a four‑part battery of LEF tasks: normal gait speed, fastest‑comfortable gait speed, isometric leg strength, and timed up‑and‑go.
These data were recorded at two study timepoints which were separated by approximately 5 years: AGES‑I (n = 3157) and AGES‑II (n = 3098).

Participants in AGES‑I were likewise administered a survey to approximate their weekly frequency of engaging in moderate‑to‑vigorous physical activity (PAAGES‑I).
Using a multivariate mediation analysis framework, linear regression models were assembled to test whether NTRA parameters mediated the longitudinal relationship between PAAGES‑I and LEFAGES‑II; all models were covariate‑adjusted for age, sex, BMI, and baseline LEF, and results were corrected for multiple statistical comparisons.
Our first series of models confirmed that all four LEF tasks were significantly related to PAAGES‑I; next, modelling the relationship between PAAGES‑I and NTRA AGES‑II identified muscle amplitude (Nm) and location (μm) as potential mediators of LEF to test.
Finally, adding these two parameters into our PAAGES‑I → LEFAGES‑II models attenuated the prior effect of PAAGES‑I; bootstrapping confirmed Nm and μm as significant partial mediators of the PAAGES‑I → LEFAGES‑II relationship, with the strongest effect found in isometric leg strength.

This work describes a novel approach toward clarifying the mechanisms that underly the relationship between physical activity and LEF in aging individuals.
Identifying Nm and μm as significant partial mediators of this relationship provides strong evidence that physical activity protects aging mobility through the preservation of both lean tissue quantity and quality.

1Department of Medicine, University of Wisconsin School of Medicine and Public Health, Madison, WI, USA.

2Wisconsin Alzheimer’s Disease Research Center, University of Wisconsin School of Medicine and Public Health, Madison, WI, USA.
3Wisconsin Alzheimer’s Institute, University of Wisconsin School of Medicine and Public Health, Madison, WI, USA.

4Institute for Biomedical and Neural Engineering, Reykjavík University, Reykjavík, Iceland.

5Icelandic Heart Association (Hjartavernd), Kópavogur, Iceland.

6Faculty of Medicine, University of Iceland, Reykjavík, Iceland.

7Department of Biomedical Sciences, University of Padova, Padua, Italy.

8Department of Science, Landspítali, Reykjavík, Iceland.
*email: kedmunds@medicine.wisc.edu Scientific Reports | (2021) 11:20173 | https://doi.org/10.1038/s41598-021-99699-7 1 Vol.

:(0123456789)www.nature.com/scientificreports Sarcopenia is consistently identified in literature as an independent risk factor for a host of deleterious health outcomes in the elderly1–5.
Although clinical definitions of the syndrome remain debated, there is growing consensus that sarcopenia is characterized primarily by the age-related loss of skeletal muscle strength and function, followed by the commensurate loss of lean tissue mass6–8.
The loss of mobility incurred by these char- acteristic changes has been directly associated with increased risk for injurious falls, physical frailty, disability, and eventual mortality9, and a growing body of evidence supports the relationship between sarcopenia and other aging-related diseases or comorbidities.
For example, decreased muscle strength and mass have been found to progress alongside advancing stages of aging-related cognitive decline10, Alzheimer’s disease (AD) dementia11,12, and brain atrophy in early preclinical AD13.

Characteristic inflammatory cytokine production and myosteatosis in sarcopenia have likewise been linked with the onset of type-2 diabetes14, hypertension15, dyslipidemia16, coronary heart disease (CHD)17, and all-type cardiovascular disease (CVD)18–21.
Taken together, these findings suggest that sarcopenia’s etiology or pathophysiology is entangled with those of other age-related diseases—a notion further reflected by sarcopenia’s high and growing global prevalence22.

Physical activity may play a critical role in not only protecting against muscle degeneration and mobility loss from sarcopenia, but in attenuating its downstream effects on other health outcomes22.
Indeed, the many advantages of maintaining an active lifestyle throughout one’s lifecourse have been repeatedly demonstrated in literature.

Regarding direct effects on skeletal muscle health, resistance exercise training in elderly individu- als has been correlated with greater muscle strength, increased lean mass, and improved myofibrillar protein accretion23.
Aerobic exercise has also been shown to improve insulin sensitivity24 and protect against a host of metabolic, cardiovascular, and oncological diseases9.

A growing body of literature further cites physical activity as an important protective factor against AD pathogenesis by promoting brain plasticity and neurotrophic factors25, improving cerebral glucose metabolism26, and reducing cerebrospinal fluid AD biomarkers27.
To interrogate the role of skeletal muscle health on the relationship between physical activity and mobility necessitates the development of standardized quantitative methods for assessing and monitoring soft tissue changes in aging individuals.

In this regard, recent investigations have realized the potential of X-ray computed tomography (CT) image analysis to quantify changes in muscle volume and quality28,29.
This is performed by the linear transformation of radiodensitometric absorption values, represented in Hounsfield units (HU), which can then be used to illustrate changes in muscle mass through segmented cross-sectional areas30–34, while average HU values can be used as a construct for muscle quality29,35,36.
Our team has recently extended this literature through the computational modelling of entire radiodensitometric distributions from CT cross-sections of the upper thigh, featuring the novel nonlinear trimodal regression analysis (NTRA) method37,38.
This method results in eleven subject-specific distribution parameters that altogether present a unique index for skeletal muscle health and have shown differential sensitivity to nutritional factors, lower extremity function (LEF), cardiovascular pathophysiology, and type-2 diabetes mellitus in healthy aging adults from the longitudinal AGES-Reykjavík study38–40.
While these works have highlighted the potential of NTRA parameters in quantitatively describing sarcopenia and its association with aging health, the mediating role of NTRA-based soft tissue radiodensity on the relation- ship between self-reported physical activity and LEF was not investigated.

This was the primary motivation for the present work, where NTRA mediation is interrogated using prospective multivariate regression models of LEF performance constructed from longitudinal AGES-Reykjavík data (Fig.
1).

Results AGES‑Reykjavik dataset covariates, LEF tasks, and NTRA parameters.

Prior to the construc- tion of multivariate regression models for the present mediation analysis, descriptive statistics, mean LEF task performance, and mean NTRA parameters from mid-thigh CT scans were assembled from AGES-I and AGES- II data, and pairwise t tests were performed between datasets to illustrate 5-year longitudinal changes.

Table 1 summarizes these results.


= 3157, n As shown, from the AGES-I sample size of n 3098 of these individuals remained as full partici- pants in AGES-II.

Aside from age, there were no significant differences between AGES-I and AGES-II covariates, but LEF task performance was significantly worse in AGES-II across all four tasks, suggesting a generalized loss of mobility in the cohort.
Significant between-dataset changes in mid-thigh NTRA parameters likewise describe a general radiodensitometric shift toward characteristic sarcopenia, illustrated primarily in lean muscle; in AGES- II, subjects’ muscle distributions were significantly wider and lower in amplitude, with an increase in positive skewness and a shift in location toward 0 HU.

This location shift also occurred in fat, while loose connective tissue distributions were significantly wider with a more negative mean location.
= NTRA parameter mediation analysis on the relationship between PAAGES‑I and LEFAGES‑II.

As the first step in mediation analysis, multivariate multiple regression models were assembled on the longitudinal relationship between physical activity recorded at AGES-I (PAAGES-I) and LEF performance recorded at AGES- II, approximately 5 years later (LEFAGES-II).

Model covariates were age, sex, BMI, and baseline LEFAGES-I defined respectively by each of the four LEF tasks.

Table 2 shows these results.


As shown in Table 2, all four PAAGES-I β coefficients were significantly related to LEFAGES-II ( p < 0.05).

Fur- thermore, all four coefficients were less than their ranked thresholds for HB correction, suggesting that all four LEFAGES-II tasks were significantly related to PAAGES-I, even after correcting for multiple comparisons.

The relationship between PAAGES‑I and hypothesized NTRA AGES‑II mediators.
Next, the relation- ship between PAAGES-I and hypothesized NTRA AGES-II mediators was tested, again using multivariate multiple regression modelling.

Results from these analyses are given in Table 3.

Scientific Reports | (2021) 11:20173 | https://doi.org/10.1038/s41598-021-99699-7 2 Vol:.

(1234567890)www.nature.com/scientificreports/ Figure 1.
Graphical summary of the present NTRA-based multiple mediation study, showing the theorized role of soft tissue radiodensity changes (NTRA parameters) on the 5-year longitudinal relationship between self- reported physical activity and LEF in elderly AGES-Reykjavík subjects.
After controlling for multiple comparisons, PAAGES-I was found to be a significant predictor ( p < 0.05) for muscle amplitude ( Nm ) and location ( µm ), suggesting their potential to serve as mediators of the relationship between PAAGES-I and LEFAGES-II.

The relationship between LEFAGES‑II and NTRA AGES‑II parameters significantly related to PAAGES‑I.
The penultimate step in mediation analysis involved testing the relationship between LEFAGES-II and Nm and µm as the two NTRA AGES-II parameters found to be significantly related to PAAGES-I.

Results from these analyses are presented in Table 4.

Here, across all four LEF tasks, Nm and µm were significant predictors ( p < 0.05) even after HB correction for multiple comparisons.
These results confirm the relationship between Nm and µm and each LEF task, suggesting the particular utility of muscle radiodensity as a quantitative correlate of LEF performance.

Testing the mediating roles of Nm and µm.

As prior regression models implicated Nm and µm as poten- tial NTRA-based mediators of the relationship between PAAGES-I and LEFAGES-II, the final step was to test their indirect mediation effects along with any direct effects from PAAGES-I on LEFAGES-II.
This was performed first by adding Nm and µm as mediators to our prior PAAGES-I → LEFAGES-II regression models (see Table 2) and testing the attenuation of PAAGES-I β coefficients.

Results from these analyses are shown in Table 5.

The addition of Nm and µm NTRA parameters attenuated prior PAAGES-I β coefficients for all four LEF tasks.

This attenuation effect differed according to task, ranging from − 18.3% (GSF) to − 32.1% (STR).
Taken together, Nm and µm are significant partial mediators of the 5-year longitudinal relationship between self-reported PA and LEF, with the strongest mediating effect seen in STR performance.

Finally, to further test Nm and µm mediation, a bootstrapping procedure using the Hayes and Preacher SPSS macro41 was performed to interrogate the significance of both indirect (mediation) effects and direct effects from PAAGES-I.

Results from this procedure are shown in Table 6.


As indicated in Table 6, bootstrapped indirect effects estimates yielded 95% confidence intervals (CI’s) whose range did not encompass zero, or the line of null effect, for any of the four LEF models.

Therefore, the indirect effects of PAAGES-I on all four LEF tasks through both Nm and µm are all significant.
Discussion Prior works have highlighted the quantitative potential of NTRA parameters from mid-thigh CT imaging in associating soft tissue radiodensity changes with aging health outcomes37–40.

However, this is the first large study to interrogate the mediating role of these radiodensity parameters as an underlying mechanism in the relation- ship between self-reported physical activity and LEF.
This was achieved by assembling prospective multivariate multiple regression models of LEF performance from longitudinal AGES-Reykjavík data, in accordance with an Scientific Reports | (2021) 11:20173 | https://doi.org/10.1038/s41598-021-99699-7 3 Vol.
:(0123456789)www.nature.com/scientificreports/ AGES-I: n = 3157 74.9 (4.8) AGES-II: n = 3098 80.0 (4.8)a Variable Age (years), mean (SD) Women, n (%) BMI (kg/m2), mean (SD) Underweight ( < 18.5 kg/m2), n (%) Normal (18.5–25 kg/m2), n (%) Overweight (25–30 kg/m2), n (%) Obese ( > 30 kg/m2), n (%) PAAGES-I, mean score (SD) (1) ‘Never’, n (%) (2) ‘Rarely’, n (%) (3) ‘Occasionally’, n (%) (4) ‘Moderate’, n (%) (5) ‘High’, n (%) LEF task 1832 (58.0) 27.3 (4.2) 18 (0.6) 949 (30.1) 1454 (46.1) 736 (23.3) 2.6 (1.6) 1258 (39.8) 527 (16.7) 243 (7.7) 578 (18.3) 551 (17.5) AGES-I Fast gait speed (GSF; m/s), mean (SD) 1.34 (0.25) Normal gait speed (GSN; m/s), mean (SD) 0.99 (0.19) Isometric leg strength (STR; N), mean (SD) 337.9 (115.9) Timed up-and-go (TUG; s), mean (SD) 11.5 (2.7) Tissue Fat Conn Muscle NTRA parameter Amplitude: N , mean (SD) Location: µ , mean (SD) Width: σ , mean (SD) Skewness: α , mean (SD) Amplitude: N , mean (SD) Location: µ , mean (SD) Width: σ , mean (SD) Amplitude: N , mean (SD) Location: µ , mean (SD) Width: σ , mean (SD) Skewness: α , mean (SD) 1807 (58.3) 26.8 (4.4) 37 (1.2) 1090 (35.2) 1309 (42.3) 662 (21.4) N/A AGES-II 1.26 (0.26)a 0.92 (0.20)a 279.5 (101.4)a 12.7 (3.8)a AGES-I 62.0 (33.1) − 117.8 (3.3) 8.2 (6.0) − 2.5 (2.1) 41.6 (8.4) − 24.1 (28.6) 25.1 (5.7) 78.0 (17.9) 61.5 (2.7) 8.6 (2.2) 2.8 (0.8) AGES-II 62.3 (33.2) − 117.2 (4.7)a 8.1 (5.7) − 2.5 (2.0) 41.8 (9.3) − 25.9 (28.2)a 24.6 (5.6)a 72.6 (17.4)a 60.9 (2.9)a 9.2 (2.6)a 3.0 (0.8)a Table 1.
AGES-I and AGES-II covariate statistics, LEF tasks, and NTRA parameters.

a Significant difference ( P < 0.05) between AGES-I and AGES-II.
Fast gait speed (GSF) Normal gait speed (GSN) Isometric leg strength (STR) Timed up-and-go (TUG) Variable β SE p-value β SE p-value β SE p-value β SE p-value PAAGES-I 0.1686 0.003 < 0.001 0.1755 0.002 < 0.001 0.1060 0.934 < 0.001 − 0.1406 0.038 < 0.001 Age Sex BMI − 0.3454 0.001 < 0.001 − 0.4046 0.001 < 0.001 − 0.2603 0.314 < 0.001 0.3467 0.014 < 0.001 − 0.3032 0.009 < 0.001 − 0.1619 0.007 < 0.001 − 0.5895 2.967 < 0.001 0.0497 0.130 0.003 − 0.1292 0.001 < 0.001 − 0.1466 0.001 < 0.001 0.0889 0.342 < 0.001 0.1454 0.015 < 0.001 LEFAGES-I 0.0180 0.018 0.279 0.0154 0.017 0.335 0.0094 0.013 0.512 0.0180 0.023 0.285 Table 2.
Multivariate multiple regression models of the relationship between PAAGES-I and LEFAGES-II.

All four PAAGES-I hypothesis tests surpassed their respective thresholds for Holm–Bonferroni (HB) correction.

established framework for mediation analysis.


From the first step in this procedure, models of the 5-year longi- tudinal relationship between PAAGES-I and LEFAGES-II indicated that all four LEF tasks were significantly related to PAAGES-I after adjusting for covariates and controlling for multiple statistical comparisons.

Next, modelling the relationship between PAAGES-I and NTRA AGES-II parameters as theorized mediators indicated Nm and µm as potential mediators of LEF.
Finally, adding Nm and µm to prior PAAGES-I → LEFAGES-II regression models resulted in the attenuation of PAAGES-I β coefficients, and bootstrapping solidified Nm and µm as significant partial mediators.

This work altogether presents a novel approach toward clarifying the nature of the relationship between physical activity and LEF in aging populations.

With respect to PAAGES-I, increased physical activity consistently incurred better LEF performance.
However, while all LEF regression models were highly significant, the variance captured by each model differed widely according to LEF task.

In this regard, results were consistently strongest in STR models and weakest with TUG.

This may be due to isometric strength tests having less propensity for Scientific Reports | (2021) 11:20173 | https://doi.org/10.1038/s41598-021-99699-7 4 Vol:.
(1234567890)www.nature.com/scientificreports/ Tiss Fat Muscle Connective Amplitude: N β SE − 0.0164 − 0.0165 0.5544 0.5884 0.0841 − 0.2152 − 0.5082 0.5135 − 0.0001 − 0.0437 − 0.3507 0.6716 0.214 0.070 0.676 0.077 0.130 0.043 0.410 0.047 0.071 0.023 0.225 0.026 Variable PAAGES-I Age Sex BMI PAAGES-I Age Sex BMI PAAGES-I Age Sex BMI Location: µ Width: σ Skewness: α p-value β 0.108 0.107 < 0.001 < 0.001 < 0.001 < 0.001 < 0.001 < 0.001 0.991 < 0.001 < 0.001 < 0.001 − 0.0233 0.0320 0.0476 0.0761 0.1076 − 0.2434 − 0.0717 0.1242 0.0243 − 0.0392 − 0.6422 − 0.3107 SE 0.054 0.018 0.170 0.020 0.031 0.010 0.099 0.011 0.224 0.074 0.709 0.081 p-value β 0.203 0.081 0.008† < 0.001 < 0.001 < 0.001 < 0.001 < 0.001 0.054 0.002 < 0.001 < 0.001 0.0409 0.0354 − 0.4222 − 0.2000 − 0.0408 0.1381 0.2243 0.4722 0.0155 − 0.1437 0.1876 − 0.0581 SE 0.058 0.019 0.184 0.021 0.025 0.008 0.081 0.009 0.063 0.021 0.198 0.022 p-value 0.011† β − 0.0323 0.028 − 0.0703 0.2820 0.2587 − 0.0276 0.1822 0.1920 0.1689 < 0.001 < 0.001 0.008† < 0.001 < 0.001 < 0.001 0.384 < 0.001 < 0.001 0.001 N/A SE 0.021 0.007 0.066 0.008 0.009 0.003 0.028 0.003 p-value 0.054 < 0.001 < 0.001 < 0.001 0.114 < 0.001 < 0.001 < 0.001 Table 3.
Multivariate multiple regression models of the relationship between PAAGES-I and hypothesized NTRA AGES-II mediators.

† p-value did not surpass HB correction for multiple comparisons.
Fast gait speed (GSF) Normal gait speed (GSN) Isometric leg strength (STR) Timed up-and-go (TUG) Variable Nm µm Age Sex BMI LEFAGES-I β 0.2990 0.0983 − 0.2805 − 0.1530 − 0.3099 0.0244 SE 3.9E−4 0.002 0.001 0.011 0.001 0.017 p-value β < 0.001 < 0.001 < 0.001 < 0.001 < 0.001 0.136 0.2722 0.1226 − 0.3385 − 0.0268 − 0.3148 0.0132 SE 2.9E−4 0.001 0.001 0.008 0.001 0.017 p-value β < 0.001 < 0.001 < 0.001 0.184 < 0.001 0.401 0.3265 0.1023 − 0.1783 − 0.4202 − 0.0978 0.0092 SE 0.126 0.522 0.318 3.61 0.408 0.012 p-value β < 0.001 < 0.001 < 0.001 < 0.001 < 0.001 0.502 − 0.0482 − 0.1606 0.2874 − 0.0615 0.2843 0.0129 SE 0.006 0.024 0.014 0.164 0.019 0.023 p-value < 0.001 < 0.001 < 0.001 0.004 < 0.001 0.437 Table 4.
Multivariate multiple regression models of the relationship between LEFAGES-II and NTRA AGES-II parameters significantly related to PAAGES-I.

All eight PAAGES-I hypothesis tests surpassed their respective thresholds for HB correction.
Fast gait speed (GSF) Normal gait speed (GSN) Isometric leg strength (STR) Timed up-and-go (TUG) β SE p-value β SE p-value β SE p-value β SE p-value Variable PAAGES-I Nm µm Age Sex BMI LEFAGES-I Prior effect PAAGES-I % Change 0.1377 0.2781 0.0853 − 0.2716 − 0.1522 − 0.2853 0.0224 0.1686 − 18.3% 0.003 3.9E−4 0.002 0.001 0.011 0.001 0.017 < 0.001 < 0.001 < 0.001 < 0.001 < 0.001 < 0.001 0.166 0.1431 0.2494 0.1092 − 0.3264 − 0.0261 − 0.2876 0.0135 0.002 2.8E−4 0.001 0.001 0.008 0.001 0.017 < 0.001 < 0.001 < 0.001 < 0.001 0.191 < 0.001 0.384 0.0720 0.3162 0.0953 − 0.1718 − 0.4192 − 0.0848 0.0104 0.894 0.126 0.522 0.317 3.60 0.410 0.012 < 0.001 < 0.001 < 0.001 < 0.001 < 0.001 < 0.001 0.443 0.003 < 0.001 0.1755 0.002 < 0.001 0.1060 0.934 < 0.001 − 18.5% − 32.1% − 0.1118 − 0.0443 − 0.1464 0.2779 − 0.0622 0.2630 0.0139 − 0.1406 − 20.3% 0.041 0.006 0.023 0.014 0.162 0.019 0.023 < 0.001 < 0.001 < 0.001 < 0.001 0.003 < 0.001 0.399 0.038 < 0.001 Table 5.
Prior PAAGES-I → LEFAGES-II regression models with Nm and µm added as theorized mediators.

All twelve hypothesis tests surpassed their respective thresholds for HB correction.

influence from external modifying factors and inter-rater variability, as reflected by conflicting discourse on the comparative reliability of gait speed measurement and TUG 42.
Nevertheless, a comparative assessment of each LEF task and its relationship with physical activity is beyond the scope of the present work, but future investiga- tion in this regard could help to clarify these findings.

The identification of Nm and µm as significant partial mediators of this relationship is strong evidence that physical activity promotes mobility in aging through the preservation of skeletal muscle quantity and quality.
The directionality of this relationship is evidenced by the signs of β coefficients for Nm and µm in each of the present regression models.

In relation to PAAGES-I, increased physical activity incurred larger lean muscle distributions Scientific Reports | (2021) 11:20173 | https://doi.org/10.1038/s41598-021-99699-7 5 Vol.
:(0123456789)www.nature.com/scientificreports/ Fast gait speed (GSF) Normal gait speed (GSN) Isometric leg strength (STR) Timed up-and-go (TUG) Effect 95% CI Effect 95% CI Effect 95% CI Effect 95% CI Indirect TotalInd Nm µm 0.0051 0.0036 0.0012 [0.0035, 0.0067] 0.0041 [0.0030, 0.0053] 2.1748 [1.5200, 2.9146] − 0.0691 [0.0023, 0.0050] 0.0026 [0.0017, 0.0036] 1.5060 [0.9527, 2.1404] − 0.0398 [0.0006, 0.0020] 0.0012 [0.0006, 0.0019] 0.5414 [0.2747, 0.8589] − 0.0233 [− 0.0905, − 0.0497] [− 0.0562, − 0.0255] [− 0.0366, − 0.0123] Bootstrap coef- ficients with 95% CIs Direct PAAGES-I 0.0228 [0.0174, 0.0281] 0.0182 [0.0142, 0.0221] 4.6179 [2.8640, 6.3717] − 0.2676 [− 0.3473, − 0.1878] Table 6.
Bootstrapped indirect and direct effects analyses for each LEFAGES-II task.

95% CI values for all indirect and direct effects did not include zero; therefore, all model effects are significant.

with mean radiodensities further from 0 HU, conferring denser, healthier lean tissue.

This relationship was likewise indicated in prospective models of LEFAGES-II, which consistently found that improved performance across all LEF tasks was related to higher Nm and µm values.
In this regard, it is important to further outline how these NTRA parameters are associated with body composition.

As noted, the NTRA method takes advantage of standard mathematical parameters that define the Gaussian distributions of radiodensity values from three canonical soft tissue types: adipose tissue, loose connective or water-equivalent tissue, and lean muscle.
This procedure generates eleven subject-specific NTRA parameters that altogether serve as a multivariate quantitative construct for capturing changes in soft tissue quantity and lean tissue quality.

These parameters can be grouped by both their relationship with soft tissue and their mathematical origin; for example, the amplitude parameter, N , is linked with soft tissue quantity through its direct association with the count of voxel elements for each tissue type.
In contrast, greater lean tissue qual- ity, reflected by limited partial volume effects from characteristic myosteatosis38, is indicated by narrower (↓σ ), less-skewed (↓α ) muscle peaks located further from 0 HU (↑µ ).
To our knowledge, this is first time that changes in soft tissue radiodensity have been shown to mediate the relationship between self-reported PA and LEF, and these findings are especially robust given both the prospective nature of this work and the large sample size offered by AGES-Reykjavík.

Another key advantage of the present NTRA-based methodology stems from its derivation from CT images.
As an imaging modality widely reported for its standardizability in diagnostic applications and pathophysiologi- cal monitoring43,44, CT-derived distributions of soft-tissue radiodensity can be directly compared across clinical contexts.

As such, NTRA-based classification is a highly reproducible method for assessing changed in body composition that can be readily built into existing CT analysis frameworks.
This could have numerous implica- tions for the field of aging research, which typically relies on estimates of body composition, such as BMI.

For example, aging-related changes in BMI have been associated with cognitive aging; Arvanitakis et al.

found that higher BMI in late-life was protective against cognitive decline, likely due to less sarcopenic muscle wasting12.

In contrast, the recent work of Kivimaki et al.


shows that midlife obesity is tied to cognitive decline and an increased risk for developing dementia10.

These age-related differences in the impact of body composition on cognitive health are problematic and cannot be resolved using traditional anthropometric estimates of body composition like BMI—a notion which further motivates the pursuit of novel soft tissue assessment methods.
Nevertheless, assessment of the present findings necessitates the consideration of several key study limitations.

Firstly, while both NTRA and LEF parameters were repeated measures available from both AGES timepoints, self-reported physical activity was only reported at AGES-I.

As such, the present study cannot conclude whether the observed decline in LEF measures is associated with potential changes in physical activity levels.
There were also a limited number of covariates considered for this work; in addition to age, sex, BMI, and LEFAGES-I it may have been of interest to test other comorbidities, biometrics, or psychosocial factors in our models45.

Likewise, as the AGES-Reykjavík study was almost entirely comprised of Caucasian subjects46, these findings may not be generalizable to other races or ethnicities.
Finally, in order to further generalize the NTRA method as a complimentary tool for studies involving physi- cal activity assessment in elderly individuals, it is essential to discuss existing measurement constructs for physical activity and some of their potential limitations.

Questionnaires, such as the one utilized in AGES-Reykjavík, have Scientific Reports | (2021) 11:20173 | https://doi.org/10.1038/s41598-021-99699-7 6 Vol:.
(1234567890)www.nature.com/scientificreports/ been validated for estimating physical activity in community health research at the population-level47,48.

However, self-reported measures have shown limited reliability and validity compared to objective methods of physical activity measurement, such as accelerometry49.
Nevertheless, systematic review findings report significant data heterogeneity from accelerometer-based studies due to differences in accelerometer design50 and disagreement on accelerometer cut-points for certain types of physical activity51.

Furthermore, wearing accelerometers and similar devices could artificially motivate study participants to perform more activity than usual49.
However, self-reporting can also be vulnerable to measurement bias from social desirability or faulty recall52, which may be particularly problematic in elderly subjects with significant memory limitations53.

Older individuals are likewise more likely to engage in light- to moderate-intensity physical activity, which are generally the most difficult types of activity to recall accurately54.
As such, exploring alternative constructs for measuring the impact of physical activity on skeletal muscle health could further situate its role in healthy aging while clarifying the underlying pathophysiological mechanisms of sarcopenia.

This study has shown the potential of NTRA-based soft tissue radiodensity analysis in this regard.

Materials and methods The AGES‑I and AGES‑II databases.
The AGES-Reykjavík study began in 2002 with the recruitment of 3316 healthy volunteer subjects between 66 and 98 years of age (mean: 77.46) to participate in a series of two multimetric health assessments separated by approximately 5 years: AGES-I and AGES-II29,55.

Written and informed consent was obtained from all participants, and ethical approval for patient data acquisition was obtained by the Icelandic Science and Ethics Committee (RU Code of Ethics, cf.
Paragraph 3 in Article 2 of the Higher Education Institution Act no.

63/2006).

All present methods were carried out in accordance with the relevant guidelines and regulations outlined in this approval.
In addition to receiving CT scans (see “CT acquisition and segmentation”) and a battery of four lower extremity function (LEF) tests (see “LEF biometrics measurement”), subjects were tasked with self-reporting their level of moderate to vigorous physical activity in the past 12-months (PAAGES-I) using a five-point Likert scale ranging from 1 (‘Never’) to 5 (‘High’).
Of the full AGES-Reykjavík recruitment, n = 3157 subjects in AGES-I had complete datasets for the purpose of this work, and n = 3098 of these subjects remained full participants in AGES-II.

Summary statistics from AGES-I and AGES-II datasets were assembled for the present study (see Table 1).

LEF biometrics measurement.


A battery of four LEF tasks were assessed as part of the AGES-Reykjavik study as previously described29,55,56: normal and fastest-comfortable gait speed (GSN and GSF, respectively), timed up-and-go (TUG), and isometric leg strength (STR).

A minimum of two LEF test leaders were tasked with recording these biometrics, and all study personnel were trained on requisite measurement procedures.
Subjects were provided with detailed instructions on each LEF task to ensure the reliability of recorded data.

Gait measurement involved participants walking a straight-line distance of 6 m in accordance with established protocol57, and speed was recorded in meters per second (m/s) over two trials, which were averaged together for each subject.

Individuals who could not complete these two trials were excluded from analyses.
The TUG test measured the time (in seconds) each subject required to stand from a seated position (height = 45.5 cm58), walk three meters, turn around, walk back to the chair, and sit down again56.

Those who could not rise from the chair by themselves or walk unaided were excluded from the test, and the duration of each subject’s first complete trial was recorded for analyses.
Finally, STR was measured in Newtons (N) via knee extension using an adjustable digital dynamometer on a fixed chair (Good Strength, Metitur, Palokka.

Finland)59.

Here, subjects’ ankles were fastened to a strain-gauge transducer to record maximal sustained force, and all trials utilized a fixed knee angle of 60° from full extension; each subject performed two 4-s trials separated by a 30-s rest period, and STR was selected as the greater of the two sustained forces in these trials.
Participants who had received surgery on their leg or hand were excluded from the tests, as well as individuals diagnosed with any ischemic heart condition two months prior to the study.

Physical activity and covariate assessment.

Physical activity was assessed in AGES-I by a self-reported questionnaire; here, subjects were asked to report their weekly participation in moderate–vigorous intensity physical activity over the previous year.
Provided examples of activities that qualified as moderate–vigorous in intensity included badminton, golf, biking, swimming, heavy gardening, weightlifting, hiking/mountain climb- ing, fast walking/heavy housework, rowing, aerobics, jogging, and running57.

Predefined physical activity cat- egories were presented on a five-level Likert scale as: (1) never, (2) rarely, (3) occasionally (weekly, but < 1 h), (4) moderate (1–3 h per week), and (5) high ( > 4 h per week).
Age, sex, and body mass index (BMI) were assessed in AGES-I and AGES-II and selected as covariates for the purpose of the present work; age and sex were assessed by questionnaire, while BMI was calculated by dividing body weight, measured in kilograms (kg), by height squared, measured in meters (m2), in accordance with established practice60.
Finally, baseline LEF performance obtained at AGES-I (LEFAGES-I) was likewise included in all LEFAGES-II models as a covariate to adjust for the pos- sibility that a prospective LEF association may be partially explained by the cross-sectional association between PA and LEF recorded during AGES-I.

CT acquisition and segmentation.


All AGES-Reykjavík participants were scanned with a 4-row CT detector system (Sensation; Siemens Medical Systems, Erlangen, Germany) using a standardized helical scan- ning protocol (120 kVp, 140 mAs, 1 mm thickness, Pitch = 1, and 4 mm collimation), as previously described61.

Standard calibration and daily quality assurance procedures were performed by multiple CT technicians to ensure the reliability and validity of acquired CT data.
The mean effective dose of radiation to participants was estimated at 1.61 mSv.

Localized scanning regions extended from the iliac crest to the knee joint, and prior to Scientific Reports | (2021) 11:20173 | https://doi.org/10.1038/s41598-021-99699-7 7 Vol.

:(0123456789)www.nature.com/scientificreports/ Figure 2.
Conceptual path diagram for the present NTRA-based multiple mediator model, showing the theorized role of soft tissue radiodensity changes (characterized by eleven NTRA parameters) on the longitudinal relationship between physical activity and our four LEF tasks.

transaxial imaging, correct scanning positions were determined by utilizing the maximum femoral length on an anterior–posterior localizer image followed by femoral long axis localization.
Following each CT image acquisi- tion, a single 10 mm cross-section was taken from the mid-thigh of each subject, located halfway between the hip acetabulum and the knee joint.

Within this scanning field-of-view, a density calibration phantom (Image Analysis, Columbia, KY, USA) was included for all subject to allow for the linear transformation of pixel attenua- tion coefficients to their equivalent radiodensity values.
This transformation resulted in the definition of subject- specific radiodensitometric distributions across the soft tissue range of − 200 to 200 HU, as described38.

Subjects with missing or incomplete CT data were excluded from analyses.

Nonlinear trimodal regression analysis (NTRA).
The computational method utilized to model radi- odensity distributions featured a form of modified nonlinear regression analysis where the distribution of radio- absorption values (measured in HU) across voxel elements is taken to be a quasi-probability density function defined by three summed Gaussian distributions: 3 i 1 = ϕ(x, Ni, µi, σi, αi) 3 1 Ni σi√2π = (x 2 µi ) 2 − 2σi e− erfc αi(x µi) σi√2 − , where N is the amplitude, µ is the location, σ is the width, and α is the skewness of each distribution—all of which are iteratively evaluated at each radiodensity bin, x .
This trimodal definition operationalizes the hypothesis that HU distributions across segmented soft tissue represent the sum of three distinct tissue types whose linear attenuation coefficients occupy approximate HU domains: fat [− 200 to − 10 HU], loose connective, water- equivalent tissue and atrophic muscle [− 9 to 40 HU], and lean muscle [41–200 HU].
The respective inwardly- sloping asymmetries of fat and muscle distributions are captured by their positive and negative skewness terms, while the central ‘connective’ tissue distribution is assumed to be standard, or non-skewed.

Next, to estimate radiodensity distribution parameters, an iterative curve-fitting procedure featuring a generalized reduced gradi- ent algorithm was employed to minimize the sum of standard errors at bin.
This general procedure, cited as the ‘NTRA’ method37, generates eleven subject-specific ‘NTRA parameters’ as a multivariate quantitative construct for capturing changes in soft tissue quantity and lean tissue quality.

Multivariate multiple regression models and statistical analyses.
To interrogate the mediating role of soft tissue NTRA parameters on the relationship between physical activity and LEF performance, we fitted a series of multivariate multiple regression models on 5-year longitudinal data captured between AGES-I and AGES-II timepoints.

These models were adjusted for age, sex, BMI, and baseline LEFAGES-I as relevant covari- ates, and to account for the simultaneous performance of multiple hypothesis tests, HB correction for multiple 0.05 .
Between-dataset differences in covariates, statistical comparisons was employed62, assuming a target αT = LEF, and NTRA parameters were likewise assembled as summary statistics, and significant cohort differences were tested using paired sample t tests.

Significance was set at p < 0.05 (adjusted for multiple comparisons) in all statistical analyses.

NTRA parameter mediation analysis on the relationship between PAAGES‑I and LEFAGES‑II.
The R software package ‘RMediation’ was employed for the present mediation analysis63, which followed a standard prospective procedure in accordance with the conceptual mediation model depicted in Fig.

2.

The first step in this procedure involved modelling the relationship between PAAGES-I and LEFAGES-II.
Next, the relationship between PAAGES-I and all eleven NTRA AGES-II parameters was modelled, and significant parameters from these models, adjusted for Type I errors, were then tested as independent predictors of LEFAGES-II.

Indirect mediation effects from these NTRA parameters were then calculated, along with any direct effects from PAAGES-I on LEFAGES-II.
This was performed by first testing the attenuation of PAAGES-I β coefficients in our prior PAAGES-I → LEFAGES-II regres- Scientific Reports | (2021) 11:20173 | https://doi.org/10.1038/s41598-021-99699-7 8 Vol:.

(1234567890)www.nature.com/scientificreports/ sion models following the addition of implicated NTRA mediators.

Finally, a bootstrapping procedure using the Hayes and Preacher SPSS macro64 was performed to interrogate the significance of indirect and direct effects.
Data availability The AGES I-II dataset cannot be made publicly available, since the informed consent signed by the participants prohibits data sharing on an individual level, as outlined by the study approval by the Icelandic National Bioethics Committee.

Requests for these data may be sent to the AGES-Reykjavik Study Executive Committee, contact: Ms. Camilla Kritjansdottir, Camilla@hjarta.is.

Received: 5 October 2020; Accepted: 8 September 2021 References 1.
Barberi, L., Scicchitano, B. M. & Musaro, A. Molecular and cellular mechanisms of muscle aging and sarcopenia and effects of electrical stimulation in seniors.

Eur.

J. Transl.

Myol.

25(4), 231–236.

https:// doi.

org/ 10.

4081/ ejtm.

2015.

5227. eColl ectio n2015 Aug24.

Review (2015).

2.

Newman, A.

B. et al.

Strength, but not muscle mass, is associated with mortality in the health, aging and body composition study cohort.

J. Gerontol.

A Biol.

Sci.

Med.

Sci.

61(1), 72–77 (2006).

3.


Goodpaster, B. H. et al.

Attenuation of skeletal muscle and strength in the elderly: The Health ABC Study.

J. Appl.

Physiol.

(Bethesda MD: 1985) 90(6), 2157–2165 (2001).

4.

Lynch, G. S. Therapies for improving muscle function in neuromuscular disorders.

Exerc.

Sport Sci.

Rev.

29(4), 141–148 (2001).

5.

Kalyani, R. R., Corriere, M. & Ferrucci, L. Age-related and disease-related muscle loss: The effect of diabetes, obesity, and other diseases.
Lancet Diabetes Endocrinol.

2(10), 819–829 (2014).

6.

Cruz-Jentoft, A. J. et al.

Sarcopenia: Revised European consensus on definition and diagnosis.

Age Ageing 48(1), 16–31 (2019).

7.

Lees, M. J., Wilson, O. J., Hind, K. & Ispoglou, T. Muscle quality as a complementary prognostic tool in conjunction with sarcopenia assessment in younger and older individuals.

Eur.

J. Appl.

Physiol.

119(5), 1171–1181 (2019).

8.

Kara, M., Ata, A. M., Kaymak, B.


& Özçakar, L. Muscle strength and function rather than muscle mass in sarcopenia.

Eur.

J. Appl.

Physiol.

119(7), 1671–1672 (2019).

9.

Marzetti, E. et al.

Sarcopenia: An overview.

Aging Clin.

Exp.

Res.

29(1), 11–17 (2017).

10.

Kivimäki, M. et al.

Body mass index and risk of dementia: Analysis of individual-level data from 1.3 million individuals.

Alzheimer’s Dement.

14(5), 601–609 (2018).

11.

Ogawa, Y. et al.


Sarcopenia and muscle functions at various stages of Alzheimer disease.

Front.

Neurol.

9, 710 (2018).

12.

Arvanitakis, Z., Capuano, A. W., Bennett, D. A.

& Barnes, L. L. Body mass index and decline in cognitive function in older black and white persons.

J. Gerontol.

Ser.

A 73(2), 198–203 (2018).

13.


Burns, J. M., Johnson, D. K., Watts, A., Swerdlow, R. H. & Brooks, W. M. Reduced lean mass in early Alzheimer disease and its association with brain atrophy.

Arch.

Neurol.

67(4), 428–433 (2010).

14.

Newman, A.

B. et al.

Sarcopenia: Alternative definitions and associations with lower extremity function.

J.

Am.

Geriatr.

Soc.

51, 1602–1609.

https:// doi.

org/ 10.

1046/j.

1532- 5415.

2003.

51534.x (2003).

15.

Brooks, S. V. & Faulkner, J.


A. Skeletal muscle weakness in old age: Underlying mechanisms.

Med.

Sci.

Sports Exerc.

26(4), 432 (1994).

16.

Maughan, R. J., Watson, J. S. & Weir, J.

Strength and cross-sectional area of human skeletal muscle.

J. Physiol.

338, 37–49 (1983).

17.

Campos, A. M. et al.

Sarcopenia, but not excess weight or increased caloric intake, is associated with coronary subclinical athero- sclerosis in the very elderly.

Atherosclerosis 258, 138–144 (2017).


18.

Reed, R. L., Pearlmutter, L., Yochum, K., Meredith, K. E. & Mooradian, A. D. The relationship between muscle mass and muscle strength in the elderly.

J.

Am.

Geriatr.

Soc.

39(6), 555–561 (1991).

19.

Jubrias, S. A., Odderson, I. R., Esselman, P. C. & Conley, K. E. Decline in isokinetic force with age: Muscle cross-sectional area and specific force.

Pflügers Arch.

Eur.

J. Physiol.

434(3), 246–253 (1997).

20.


Overend, T. J., Cunningham, D. A., Kramer, J. F., Lefcoe, M. S. & Paterson, D. H. Knee extensor and knee flexor strength: Cross- sectional area ratios in young and elderly men.

J. Gerontol.

47(6), M204–M210 (1992).

21.

Han, P. et al.

The increased risk of sarcopenia in patients with cardiovascular risk factors in Suburb-Dwelling older Chinese using the AWGS definition.

Sci.

Rep. 7(1), 9592 (2017).

22.

Cruz-Jentoft, A. J. et al.


Prevalence of and interventions for sarcopenia in ageing adults: a systematic review.

Report of the Inter- national Sarcopenia Initiative (EWGSOP and IWGS).

Age Ageing 43(6), 748–759 (2014).

23.

Landi, F., Marzetti, E., Martone, A. M., Bernabei, R. & Onder, G. Exercise as a remedy for sarcopenia.

Curr.

Opin.

Clin.

Nutr.

Metab.

Care 17(1), 25–31 (2014).

24.

Landi, F., Onder, G. & Bernabei, R. Sarcopenia and diabetes: Two sides of the same coin.

J.

Am.

Med.
J.

Am.

Med.

Dir.

Assoc.

14(8), 540–541 (2013).

25.

Larson, E. B., Wang, L. I., Bowen, J. D., McCormick, W. C., Teri, L., Crane, P., & Kukull, W. (2006).

Exercise is associated with reduced risk for incident dementia among persons 65 years of age and older.

Annals of internal medicine, 144(2), 73–81.

26.

Dougherty, R. J. et al.

Moderate physical activity is associated with cerebral glucose metabolism in adults at risk for Alzheimer’s disease.

J. Alzheimers Dis.

58(4), 1089–1097 (2017).

27.


Law, L. L., Rol, R. N., Schultz, S. A., Dougherty, R. J., Edwards, D. F., Koscik, R. L., et al.

(2018).

Moderate intensity physical activity associates with CSF biomarkers in a cohort at risk for Alzheimer’s 28.

Goodpaster, B. H., Kelley, D. E., Thaete, F. L., He, J.

& Ross, R. Skeletal muscle attenuation determined by computed tomography is associated with skeletal muscle lipid content.

J. Appl.

Physiol.

89(1), 104–110 (2000).

29.


Goodpaster, B. H. et al.

The loss of skeletal muscle strength, mass, and quality in older adults: The health, aging and body composi- tion study.

J. Gerontol.

A Biol.

Sci.

Med.

Sci.

61(10), 1059–1064 (2006).

30.

Carraro, U. et al.

Persistent muscle fiber regeneration in long term denervation.

Past, present, future.

Eur.

J. Transl.

Myol.

25(2), 77–92 (2015).

31.

Gargiulo, P. et al.


Quantitative color three-dimensional computer tomography imaging of human long-term denervated muscle.

Neurol.

Res.

32(1), 13 (2010).

32.

Helgason, T. et al.

Monitoring muscle growth and tissue changes induced by electrical stimulation of denervated degenerated muscles with CT and stereolithographic 3D modeling.

Artif.

Organs.

29(6), 440–443 (2005).

33.

Snijder, M. et al.


Low subcutaneous thigh fat is a risk factor for unfavourable glucose and lipid levels, independently of high abdominal fat.

The Health ABC Study.

Diabetologia 48(2), 301–308 (2005).

34.

Carraro, U., Edmunds, K. J.

& Gargiulo, P. 3D false color computed tomography for diagnosis and follow-up of permanently denervated human femoral muscles submitted to functional electrical stimulation.

Eur.

J. Transl.

Myol.

25(2), 129–140 (2015).


Scientific Reports | (2021) 11:20173 | https://doi.org/10.1038/s41598-021-99699-7 9 Vol.

:(0123456789)www.nature.com/scientificreports/ 35.

Edmunds, K. J., & Gargiulo, P. Imaging approaches in functional assessment of implantable myogenic biomaterials and engineered muscle tissue.

Eur.

J. Transl.

Myol., 25(2).

(2015).

36.

Hicks, G. E. et al.


Cross-sectional associations between trunk muscle composition, back pain, and physical function in the health, aging and body composition study.

J. Gerontol.

Ser.

A Biol.

Sci.

Med.

Sci.

60(7), 882–887 (2005).

37.

Edmunds, K. J., Árnadóttir, Í., Gíslason, M. K., Carraro, U., Gargiulo, P. Nonlinear trimodal regression analysis of radiodensito- metric distributions to quantify sarcopenic and sequelae muscle degeneration.

Comput.

Math.

Methods Med.


8932950 (2016).

38.

Edmunds, K. J. et al.

Advanced quantitative methods in correlating sarcopenic muscle degeneration with lower extremity function biometrics and comorbidities.

PLoS One 13(3), e0193241 (2018).

39.

Ricciardi, C. et al.

Assessing cardiovascular risks from a mid-thigh CT image: A tree-based machine learning approach using radiodensitometric distributions.

Sci.

Rep. 10(1), 1–13 (2020).

40.

Recenti, M. et al.


Healthy aging within an image: Using muscle radiodensitometry and lifestyle factors to predict diabetes and hypertension.

IEEE J. Biomed.

Health Inform.

25(6), 2103–2112 (2021).

41.

Preacher, K. J.

& Hayes, A. F. SPSS and SAS procedures for estimating indirect effects in simple mediation models.

Behav.

Res.

Methods Instrum.

Comput.

36(4), 717–731 (2004).

42.


Cameron, M. H., & Monroe, L. Physical Rehabilitation-E-Book: Evidence-Based Examination, Evaluation, and Intervention.

(Elsevier Health Sciences, 2007).

43.

Dickerman, B.

A. et al.

Body fat distribution on computed tomography imaging and prostate cancer risk and mortality in the AGES-Reykjavik study.

Cancer 125(16), 2877–2885 (2019).

44.

Yoshizumi, T. et al.

Abdominal fat: Standardized technique for measurement at CT. Radiology 211(1), 283–286 (1999).

45.
45.

Okonkwo, O. C. et al.

Physical activity attenuates age-related biomarker alterations in preclinical AD.

Neurology 83(19), 1753–1760 (2014).

46.

Keyak, J. H. et al.

Hip load capacity and yield load in men and women of all ages.

Bone 137, 115321 (2020).

47.

Speck, B. J.

& Looney, S. W. Self-reported physical activity validated by pedometer: A Pilot Study.

Public Health Nurs.

23(1), 88–94 (2006).

48.

Carter-Nolan, P. L. et al.


Validation of physical activity instruments: Black Women’s Health Study.

Ethn.

Dis.

16(4), 943 (2006).

49.

Field, M. J., Gebruers, N., Shanmuga Sundaram, T., Nicholson, S., & Mead, G. (2013).

Physical activity after stroke: a systematic review and meta-analysis.

International Scholarly Research Notices, 2013(1), 1–13.

50.

Skender, S. et al.

Accelerometry and physical activity questionnaires—A systematic review.

BMC Public Health 16(1), 515 (2016).

51.


Prince, S. A. et al.

A comparison of accelerometer cut-points among individuals with coronary artery disease.

PLoS One 10(9), e0137759 (2015).

52.

Prince, S. A. et al.

A comparison of direct versus self-report measures for assessing physical activity in adults: A systematic review.

Int.

J. Behav.

Nutr.

Phys.

Act.

5, 56 (2008).

53.

Pruitt, L. A. et al.

Use of accelerometry to measure physical activity in older adults at risk for mobility disability.
J.

Aging Phys.

Act.

16(4), 416–434 (2008).

54.

Baranowski, T. Validity and reliability of self report measures of physical activity: An information-processing perspective.

Res.

Q. Exerc.

Sport 59(4), 314–327 (1988).

55.

Harris, T. B. et al.

Age, gene/environment susceptibility–Reykjavik study: Multidisciplinary applied phenomics.

Am.

J. Epidemiol.

165(9), 1076–1087 (2007).

56.

Chang, M. et al.


Midlife physical activity preserves lower extremity function in older adults: Age gene/environment susceptibil- ity–Reykjavik study.

J.

Am.

Geriatr.

Soc.

61(2), 237–242 (2013).

57.

Guralnik, J. M., Ferrucci, L., Simonsick, E. M., Salive, M. E. & Wallace, R. B. Lower-extremity function in persons over the age of 70 years as a predictor of subsequent disability.

N. Engl.

J. Med.

332(9), 556–562 (1995).

58.


Siggeirsdottir, K., Jonsson, B. Y., Jonsson, H. Jr. & Iwarsson, S. The timed ‘Up & Go’is dependent on chair type.

Clin.

Rehabil.

16(6), 609–616 (2002).

59.

Cvecka, J. et al.

Physical activity in elderly.

Eur.

J. Transl.

Myol.

25(4), 249 (2015).

60.

Mijnarends, D. M. et al.

Physical activity and incidence of sarcopenia: the population-based AGES—Reykjavik Study.

Age Ageing 45(5), 614–620 (2016).

61.

Johannesdottir, F. et al.


Distribution of cortical bone in the femoral neck and hip fracture: A prospective case-control analysis of 143 incident hip fractures; the AGES-REYKJAVIK Study.

Bone 48(6), 1268–1276 (2011).

62.

Frewen, P. A., Schmittmann, V. D., Bringmann, L. F. & Borsboom, D. Perceived causal relations between anxiety, posttraumatic stress and depression: Extension to moderation, mediation, and network analysis.

Eur.

J. Psychotraumatol.

4(1), 20656 (2013).

63.


Tofighi, D. & MacKinnon, D. P. RMediation: An R package for mediation analysis confidence intervals.

Behav.

Res.

Methods 43(3), 692–700 (2011).

64.

Preacher, K. J.

& Hayes, A. F. Asymptotic and resampling strategies for assessing and comparing indirect effects in multiple media- tor models.

Behav.

Res.

Methods 40(3), 879–891 (2008).


Acknowledgements The authors wish to thank the University Hospital Landspitali in Reykjavik for its infrastructural support.

Author contributions K.J.E.

performed regression modelling and statistical analyses while leading the literature review and manuscript preparation.

O.C.O.

and S.R.L.

advised the mediation analysis procedure and contributed their knowledge of physical activity in promoting aging health.

S.S. and V.G.
S.S. and V.G.

incorporated their complete knowledge of the AGES- Reykjavík study and coordinated dataset management.

U.C.

contributed his knowledge on muscle physiology and translational myology.

P.G.

supervised and coordinated all aspects of this work, and all authors contributed to editing and revising the manuscript.

Competing interests The authors declare no competing interests.

Additional information Correspondence and requests for materials should be addressed to K.J.E.
Reprints and permissions information is available at www.nature.com/reprints.

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

Scientific Reports | (2021) 11:20173 | https://doi.org/10.1038/s41598-021-99699-7 10 Vol:.
(1234567890)www.nature.com/scientificreports/ 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:// creat iveco mmons.

org/ licen ses/ by/4.

0/.


© The Author(s) 2021 Scientific Reports | (2021) 11:20173 | https://doi.org/10.1038/s41598-021-99699-7 11 Vol.

:(0123456789)www.nature.com/scientificreports/