Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

MRI-based anatomical characterisation of lower-limb muscles in older women

  • Erica Montefiori ,

    Contributed equally to this work with: Erica Montefiori, Barbara M. Kalkman

    Roles Formal analysis, Methodology, Software, Validation, Visualization, Writing – original draft, Writing – review & editing

    e.montefiori@sheffield.ac.uk

    Affiliations Department of Mechanical Engineering, University of Sheffield, Sheffield, United Kingdom, INSIGNEO Institute for in silico Medicine, University of Sheffield, Sheffield, United Kingdom

  • Barbara M. Kalkman ,

    Contributed equally to this work with: Erica Montefiori, Barbara M. Kalkman

    Roles Data curation, Formal analysis, Methodology, Writing – review & editing

    Affiliations Department of Mechanical Engineering, University of Sheffield, Sheffield, United Kingdom, INSIGNEO Institute for in silico Medicine, University of Sheffield, Sheffield, United Kingdom

  • William H. Henson,

    Roles Methodology, Writing – review & editing

    Affiliations Department of Mechanical Engineering, University of Sheffield, Sheffield, United Kingdom, INSIGNEO Institute for in silico Medicine, University of Sheffield, Sheffield, United Kingdom

  • Margaret A. Paggiosi,

    Roles Conceptualization, Data curation, Methodology, Resources, Supervision, Writing – review & editing

    Affiliations INSIGNEO Institute for in silico Medicine, University of Sheffield, Sheffield, United Kingdom, Department of Oncology and Metabolism, Centre for Integrated research in Musculoskeletal Ageing, Mellanby Centre for Bone Research, University of Sheffield, Sheffield, United Kingdom

  • Eugene V. McCloskey,

    Roles Conceptualization, Funding acquisition, Investigation, Project administration, Resources, Supervision, Writing – review & editing

    Affiliations INSIGNEO Institute for in silico Medicine, University of Sheffield, Sheffield, United Kingdom, Department of Oncology and Metabolism, Centre for Integrated research in Musculoskeletal Ageing, Mellanby Centre for Bone Research, University of Sheffield, Sheffield, United Kingdom

  • Claudia Mazzà

    Roles Conceptualization, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Writing – review & editing

    Affiliations Department of Mechanical Engineering, University of Sheffield, Sheffield, United Kingdom, INSIGNEO Institute for in silico Medicine, University of Sheffield, Sheffield, United Kingdom

Abstract

The ability of muscles to produce force depends, among others, on their anatomical features and it is altered by ageing-associated weakening. However, a clear characterisation of these features, highly relevant for older individuals, is still lacking. This study hence aimed at characterising muscle volume, length, and physiological cross-sectional area (PCSA) and their variability, between body sides and between individuals, in a group of post-menopausal women. Lower-limb magnetic resonance images were acquired from eleven participants (69 (7) y. o., 66.9 (7.7) kg, 159 (3) cm). Twenty-three muscles were manually segmented from the images and muscle volume, length and PCSA were calculated from this dataset. Personalised maximal isometric force was then calculated using the latter information. The percentage difference between the muscles of the two lower limbs was up to 89% and 22% for volume and length, respectively, and up to 84% for PCSA, with no recognisable pattern associated with limb dominance. Between-subject coefficients of variation reached 36% and 13% for muscle volume and length, respectively. Generally, muscle parameters were similar to previous literature, but volumes were smaller than those from in-vivo young adults and slightly higher than ex-vivo ones. Maximal isometric force was found to be on average smaller than those obtained from estimates based on linear scaling of ex-vivo-based literature values. In conclusion, this study quantified for the first time anatomical asymmetry of lower-limb muscles in older women, suggesting that symmetry should not be assumed in this population. Furthermore, we showed that a scaling approach, widely used in musculoskeletal modelling, leads to an overestimation of the maximal isometric force for most muscles. This heavily questions the validity of this approach for older populations. As a solution, the unique dataset of muscle segmentation made available with this paper could support the development of alternative population-based scaling approaches, together with that of automatic tools for muscle segmentation.

Introduction

The characterisation of the ability of individual muscles to produce force is of particular relevance in older individuals, for whom ageing-associated muscle loss (sarcopenia) can significantly affect the ability of a muscle to produce strength [13]. The ability of a muscle to generate force depends on its fibre composition and characteristics, and on its structural and morphological features [4]. The loss of muscle strength at older ages has been explained by a reduction in muscle mass [5], an increase in slower muscle fibres [4], a higher percentage of intramuscular fat, in combination with a smaller physiological cross-sectional area (PCSA) [4,6].

Being able to quantify lower-limb muscle forces during dynamic tasks can help to understand the capability of an individual to control a movement, with relevant application in the prediction of risk of fall and fractures in older people. Musculoskeletal (MSK) models have been increasingly adopted for this purpose [710]. These models, however, often rely on a number of assumptions about anatomical features of muscles, neglecting possible variabilities between subjects and within subject (i.e. body asymmetry) which knowingly affect the accuracy of their outputs [1115]. Additionally, muscle properties such as the maximal isometric force (Fmax) are often derived from cadaver-based dissection studies [16,17], hence neglecting population-specific features like reduction of muscle strength in the elderly.

Estimated muscle forces obtained from MSK models are known to be sensitive to variations in architectural musculotendon parameters [18]. Moderate and muscle-specific sensitivity to Fmax has been previously reported in young adults, both for generic-scaled [19] and subject-specific models [13]. In these studies, Fmax was only made to vary within small ranges [19] or proportionally for all muscles [13], and the effects of muscle- or subject-specificity have not been investigated [20]. Additionally, MSK models do not typically account for loss of muscle strength typically associated to ageing [3,21], which can be both subject- and muscle-specific [22].

In an attempt to overcome the above limitations, clinical measurements of muscle strength, such as those from hand-held dynamometer measurements of grip strength, can be integrated in MSK models [23]. These, however, provide an overall indication of strength rather than muscle-specific strength properties. Medical imaging such as Magnetic Resonance Imaging (MRI) has been successfully adopted for deriving individual muscle volume and muscle length through image segmentation [24]. The ratio of muscle volume and length is proportional to the muscle PCSA [25] and hence to the maximal isometric force a muscle can generate [26].

Except for tendon slack length, which cannot yet be quantified with any routine non-invasive techniques, a full characterisation of the muscle parameters based on MRI is certainly feasible [27], but not commonly pursued due to the time and repeatability challenges associated with image processing. As a result, very little is known about specific characteristics of these parameters, especially in older individuals. The main aim of this study was to investigate lower-limb muscle anatomical characteristics, including volume, length, and PCSA, in a group of post-menopausal women. When enough knowledge about individual muscles parameters is available, this could be used to either build population-based statistical models [28] or, as recently proposed by Handsfield et al. [24] for young individuals, establish their relationship with the body mass or length, overcoming the need for segmenting individual muscles in future applications. The second aim of this paper was to verify the suitability of this approach and to provide the community with a fully characterised database including 3D muscle and bone geometries as obtained from lower-limb MRI from a group of post-menopausal women, in the attempt to foster the community efforts towards the development of automatic image processing and modelling tools.

Methods

Participants and data acquisition

Eleven post-menopausal women (mean (standard deviation, SD): 69 (7) y. o., 66.9 (7.7) kg, 159 (3) cm) with no movement limitations were recruited by the Metabolic Bone Centre, Northern General Hospital in Sheffield, UK as part of larger studies (Multisim and Multisim 2, EP/K03877X/1 and EP/S032940/1, https://epsrc.ukri.org). Inclusion criteria were having a bone mineral density T-score at the lumbar spine or total hip (whichever was the lower value) less than or equal to -1. Bone mineral density was measured by dual energy x-ray absorptiometry using a Discovery A densitometer (Hologic Inc., Bedford, MA, USA). Exclusion criteria were: body mass index (BMI) <18 or >35, history of or current conditions known to affect bone metabolism and bone mineral density, history of or current neurological disorders, prescription of oral corticosteroids for more than three months within the last year, history of any long term immobilization (>3 months), conditions that prevent the acquisition of musculoskeletal images, use of medications or treatment known to affect bone metabolism other than calcium/vitamin D supplementation and alcohol intake greater than 21 units per week. The study was approved by the East of England—Cambridgeshire and Hertfordshire Research Ethics Committee and the Health Research Authority and was conducted in accordance with the Declaration of Helsinki (October 2000). Written informed consent was obtained from all participants.

During a hospital visit, full lower-limb MRI was collected using a Magnetom Avanto 1.5 T scanner (Siemens, Erlangen Germany). A T1-weighted scanning sequence was used with an echo time of 2.59 ms, a repetition time of 7.64 ms, flip angle of 10 degrees and voxel sizes of 1.1x1.1x5.0 mm for the long bones and 1.1x1.1x3.0 mm for the joints. In this occasion participants’ lower-limb dominance was determined asking them “If you kicked a football which foot would you use?” [29].

Data processing

Muscle segmentation.

Lower-limb bones were segmented within the MRI scans using Mimics 20.0 (Materialise, Leuven, Belgium). In each limb, 30 muscles were segmented, initially using the automated muscle segmentation toolbox (Mimics Research 20.0, Materialise, Belgium), followed by manual adjustments when necessary. Inter-operator repeatability of the muscle segmentation procedure was assessed by calculating the ratio between SD and mean (referred to as coefficient of variation, CoV) of the muscle volumes (VM) calculated by three different operators on a subset of three participants. According to literature suggestions [30,31], values of CoV can be considered as acceptable when below 10%. Using a conservative approach, for those muscles where inter-operator CoV was higher than 5% we also performed an intra-operator analysis, asking the same operator to repeat the segmentation three times on the same dataset. Following the latter analysis, we discarded all the muscles with non-acceptable repeatability (CoV > 10%). The Psoas major muscle was removed from the repeatability study since it was partially cut off from the MRI field of view in some cases. Similarly, the foot extensors and flexors were not evaluated, since their external boundaries were not identifiable in many of the MRI datasets.

Calculation of the maximal isometric force.

Two different approaches were used to calculate Fmax. Firstly, a linear scaling of Fmax based on lower-limb mass [32], which is typically used in MSK models when individual muscle geometries are not available (Lower-limb mass-based scaling, LLMS). Secondly, Fmax was calculated as a function of muscle PCSA, calculated from individual muscle volumes and length (Volume and length-based scaling, VLS).

In the LLMS approach [32], Fmax was linearly scaled to the lower-limb mass according to (1): (1) where mLL is the mass of the lower limbs of the subject, calculated as a product of the volume of the lower limbs (estimated from the MRI) and the density of the tissue [33]), mLLGen is the mass of the lower limbs of the generic OpenSim model gait2392 [17] and FmaxGen is the default Fmax of each muscles in the gait2392 model. An equivalent estimate of Fmax could be theoretically obtained in the absence of MRI by estimating mLL after a scaling procedure (e.g. using the Scaling Tool in OpenSim [34]).

In the VLS approach, muscle segmentations were used to calculate the muscle volume (VM) and the anatomical muscle length (lM) was calculated as the length of the centreline from the 3D muscle segmentation. This was generated as the line connecting the points representing the topological skeleton of each muscle cross section in the 2D MRI slices. A smooth curve was fitted to the centreline using a moving average filter, with the span of the filter being selected individually for each muscle. Values for lM were then denoted as the arc length of the fitted smoothed curve constituting the centreline of the 3D segmentations. All above computations were performed in MATLAB R2019b (The Mathworks Inc., Natick, MA, USA). VM and lM were then used to calculate the muscle PCSA according to (2): (2) where k is the ratio between a muscle optimal fibre length () and length, as taken from the literature [25].

Values of VM and PCSA were compared to those available in the literature for healthy young adults [24] and cadavers [25,27].

Fmax was calculated as a product of the PCSA described in Eq (1) and the specific tension (σ = 61 N/cm2, [16,35]), [26]: (3) For the Glutei and Adductor magnus, 1/3 of the total Fmax value was attributed to each of the three bundles constituting the muscle and used for comparison to the values obtained with the LLMS method.

Statistical analysis.

All variables were tested for normality using the one-sample Kolmogorov-Smirnov test in MATLAB and null hypothesis were then consistently tested using either a student’s t test in the case of normally distributed data or a Wilcoxon signed-rank test in the case of non-normally distributed data. To discard the hypothesis of anatomical symmetry, VM, lM and PCSA of the muscles belonging to the right and left limb were compared. The percentage difference between the values in the right and left limb was also quantified for all the muscles and all the subjects. CoV was calculated for each muscle to quantify the inter-subject variability.

Linear regressions were computed between total lower limb muscle volume (VTOT equal to the sum of the muscles whose segmentation resulted repeatable) and lower-limb mass, body mass, height, and BMI.

The effect of accounting for individual muscle geometry on the calculated Fmax was quantified by comparing the Fmax values obtained using the LLMS and VLS approaches. Percentage difference between Fmax estimated with the two methods was calculated. Significance level α was set to 0.05 for all statistical tests.

Results

Muscle segmentation

The inter-operator analysis provided higher CoV than the intra-operator analysis (Table 1) for all the muscles tested. The Gastrocnemii and Vastus medialis were easily identifiable and led to very high inter-operator repeatability. The Peronei had the worst inter-operator CoV (close to 50%). Even though better results were found for the intra-operator analysis for the Peroneus brevis (CoV = 7.6%), this was not the case for the Peroneus longus (CoV = 10.9%), which was removed from further analysis together with the Gluteus minimus (CoV = 21.6%).

In light of the high inter-operator differences, only muscle segmentations generated by the same single expert operator were used for the following analyses.

Muscle anatomical parameters

From the dominance test, all participants resulted right limb dominant.

All investigated parameters were not normally distributed; therefore, non-parametric tests were selected for the statistical analysis. An evident intra- and inter-subject variability was observed for VM and lM, as depicted by the bar plots in Figs 1 and 2 (individual VM and lM values are available as Supplementary material).

thumbnail
Fig 1. Muscle volume variability.

Median (minimum, maximum) of muscle volume for the right and left limb (significant difference between limbs: * p<0.05, **p<0.01). Individual percentage difference between the limbs is reported as a bar plot where each bar represents a participant: blue positive (red negative) values show that the right leg is bigger (smaller). Minimum and maximum percentage difference across the subjects is reported for each muscle.

https://doi.org/10.1371/journal.pone.0242973.g001

thumbnail
Fig 2. Muscle length variability.

Median (minimum, maximum) of muscle length for the right and left limb (significant difference between limbs: **p<0.01). Individual percentage difference between the limbs is reported as a bar plot where each bar represents a participant: blue positive (red negative) values show that the right leg is bigger (smaller). Minimum and maximum percentage difference across the subjects is reported for each muscle.

https://doi.org/10.1371/journal.pone.0242973.g002

The percentage difference of VM between the two limbs was above 85% for the Gracilis in one subject and for the Rectus femoris in another subject. A significant difference between the two limbs was found for the VM of the Sartorius, Gluteus maximus, Adductor magnus, and Vastus lateralis, with lower values in the left limb. Between-subject CoV (see Supplementary material) ranged between 14% (Vastus medialis) and 36% (Sartorius).

The percentage difference of lM between the two limbs was up to 22% (Adductor brevis). A significant difference between the two limbs was observed for the lM of the Gluteus medius and Vastus lateralis, with lower values in the left limb. Between-subject CoV (see Supplementary material) ranged between 3% (Sartorius) and 13% (Gastrocnemius lateralis).

Mean and SD of the VM are reported in Table 2 for the sake of comparison with literature data. Overall, our values were higher than dissection-based muscle volumes from elderly cadavers [25] but smaller than muscle volumes from mixed-age cadavers [27] and MRI-based muscle volumes from healthy young adults [24] both of mixed sexes and females only.

thumbnail
Table 2. Comparison of muscle volume to literature values.

https://doi.org/10.1371/journal.pone.0242973.t002

Among the tested anthropometric parameters height and lower-limb mass did not significantly correlate with VTOT (Fig 3), whereas BMI and body mass showed significant correlations, with coefficient of determination, R2 = 0.44 (p = 0.003) and R2 = 0.50 (p = 0.004), respectively.

thumbnail
Fig 3. Linear regression between muscle volume and anthropometric parameters.

Linear regression and coefficients of determination (R2) between total lower-limb muscle volume and body mass (R2 = 0.50, p = 0.003, left), height (R2 = 0.02, p > 0.05, middle-left), lower-limb mass (R2 = 0.14, p > 0.05, middle-right), BMI (R2 = 0.44, p = 0.004¸ right).

https://doi.org/10.1371/journal.pone.0242973.g003

The percentage difference in PCSA (Fig 4) between right and left limb ranged between -84% for the Gracilis and Rectus femoris (with smaller PCSA in the right limb) and 38% for the Gracilis (with bigger PCSA in the right limb). Only for Gastrocnemius medialis, Gluteus maximus and medius and Soleus between-limb variations were below 20% for all the subjects. The values found for the participants in this study were similar to those reported by other authors for mixed-age/sex cadavers [25,27], but smaller than those from healthy young adults of both sexes [24] (Fig 5).

thumbnail
Fig 4. PCSA variability.

Median (minimum, maximum) of the physiological cross-sectional area (PCSA) for 23 lower-limb muscles for eleven subjects in our study (n = number of limbs). PCSA are derived from the segmented VM and lM and using the average optimal fibre length to muscle length ratio proposed by Ward et al. [25]; *PCSA of the Tensor fasciae latae was calculated setting the optimal fibre length to muscle length ratio equal to 1 (as proposed by Handsfield et al. [24]) since the actual values were not available from the literature source. Minimum and maximum percentage difference across the subjects is reported for each muscle.

https://doi.org/10.1371/journal.pone.0242973.g004

thumbnail
Fig 5. PCSA distribution and comparison to literature data.

Distribution of the PCSA for the 22 limbs analysed in this study (grey violin plots) compared to PCSA values from literature. Red circles represent individual data points for three cadavers as calculated by Charles et al. [27]. Blue diamonds represent mean PCSA values for twenty-one cadavers as calculated by Ward et al. [25] and divided by the cosine of the mean pennation angle reported by the same authors. Green squares with deviation error bars represent PCSA values estimated by Handfield et al. [24] from MRI segmentation of thirty-two healthy young adults.

https://doi.org/10.1371/journal.pone.0242973.g005

Maximal isometric force

The Fmax calculated from the VLS approach differed from that of the LLMS by up to 400% (Biceps femoris short head) for individual subjects (Fig 6), with overall smaller estimates of Fmax with the VLS model. On average, the percentage difference between the two approaches was between -176% (for the Iliacus where Fmax was smaller in the VLS) and 36% (for the Adductor magnus II where Fmax was bigger in the VLS). Differences were found significant for all muscles except for Gluteus maximus I and III, Adductor magnus III, Biceps femoris long head, Semimembranosus, Rectus femoris, and Peroneus brevis.

thumbnail
Fig 6. Maximal isometric force calculated with the VLS and LLMS approach.

Median (minimum, maximum) of the maximal isometric force for the VLS and LLMS approaches with p values representing the statistical significance of Wilcoxon test. Individual percentage difference of Fmax between VLS and LLMS reported as a bar plot where each bar represents a participant: green positive (orange negative) bars show that the value is bigger (smaller) with the VLS approach. Minimum and maximum percentage difference across the subjects is reported for each muscle. * These values correspond to one third of the total muscle Fmax.

https://doi.org/10.1371/journal.pone.0242973.g006

Discussion

This study aimed to quantify lower-limb muscle anatomical characteristics from medical images in a group of post-menopausal women. To this purpose, the 3D geometries of 23 lower-limb muscles segmented from MRI from a cohort of eleven post-menopausal women were used to assess inter- and intra-individual differences and compared to existing literature data. The use of image segmentation for the calculation of muscle parameters is complicated by the time and repeatability challenges associated with this technique. However, broadening the knowledge of muscle anatomical characteristics could support the development of tools (e.g. population-based statistical models [28] or regression models [24]) to overcome the need for segmentation.

This is, to our knowledge, the first study providing a quantification of lower-limb muscle volumes and lengths in older women, and a thorough assessment of the differences observable both between body sides and across individuals. An ultrasound-based study quantified up to 24% of muscle thickness asymmetry in abdominal muscles in healthy individuals of different ages [37], suggesting that analogous results could be expected in the lower limbs. When comparing the two limbs of each subject in our cohort, we observed differences of up to 85% for VM and of up to 22% for lM (Figs 1 and 2). Except for very few muscles (Sartorius, Gluteus maximus, Adductor magnus, Vastus lateralis) which were significantly bigger on the right side, no recognisable pattern was observed across the cohort to be associated with limb dominance. In fact, both muscle volumes and lengths were notably variable in the population. This clearly indicates that care should be taken in assuming limb symmetry when assigning musculotendon parameters, even in healthy populations.

Even though different approaches to the image segmentation may have affected the estimate of the muscle parameters, the comparison to MRI-based values from the literature [24] led to valuable insights. Despite the average height and weight of our participants being smaller than those previously reported for an ex-vivo cohort [25], slightly larger VM were found (Table 2). This could be explained by the loss in muscle mass in cadavers [36]. On the contrary, our VM was smaller than that estimated in-vivo from MRI in healthy young adults (25.5 (11.1) y. o.) [24], which might be explained by both younger age and mixed-sex participants. In fact, when isolating the female component from the young population, smaller average VM and SD were still observed in our cohort for all the muscles. This explains the smaller inter-subject variability (as quantified by CoV) found in our study, i.e. between 14% (Vastus medialis) and 36% (Sartorius), compared to literature values for healthy young mixed-sex adults (quantified between 20% and 40% from the reported mean and SD) [25] and even more when isolating the female component (except for the Tensor fasciae latae muscle). The VM calculated from our cohort remained consistently smaller to those from young females, except for the Vastus intermedius, likely due to ageing-related muscles volume loss [3,21].

In order to overcome the need for individual muscle segmentations to estimate muscle-specific parameters in MSK models, Handsfield et al. [24] proposed a series of regression equations linearly correlating individual muscle volume to participants’ total lower-limb muscle volume, body mass and height. Lower correlations were quantified in this study (Fig 3), likely due to having included only 23 instead of 35 lower-limb muscles. This discrepancy could also be preferential weakening or atrophy of certain muscles caused by ageing [38], an hypothesis which seems to be confirmed by the lower volumes found in our cohort when compared to younger females. Surprisingly, VTOT correlated more strongly with total body mass than with lower-limb mass, suggesting that scaling muscle forces based on lower-limb mass (LLMS) [32,39] might not be a suitable approach in an older population, and a simple scaling to body mass should be preferred in the absence of MRI.

The maximal force that a muscle can produce is highly affected by its PCSA [26]. Since optimal fibre length could not be calculated from available MRI data, the PCSA was here calculated by scaling lM according to ex-vivo literature values from an older population [25]. This led to PCSA values in agreement with literature [24,25,27], except for bigger values for the Gluteus maximus and smaller values for the Iliacus (Fig 5). The PCSA of the Sartorius muscle presented a 37% of CoV between the subjects, due to high variability in its volume and small variability in its length. This was also the only muscle showing significantly different Fmax between the body sides at group level, with larger values in the dominant limb. Previous studies highlighted intra-subject variability in the tendon-to-muscle belly length ratio as well as in the location of the widest part of the muscle along its axis [40], therefore confirming our findings.

The specific tension (σ) of a muscle also contributes to the estimate of Fmax. The choice of setting σ to 61 N/cm2 was suggested by previous literature where this value was proposed for elderly populations [16,35]. Sensitivity of models to this parameter was previously tested by Valente et al. [13], finding a moderate effect on the model output. In the effort of maximally personalizing muscle parameters, individual values for the specific tension should be obtained for different subjects and different muscles, however such a measure is not currently available in-vivo. The use of dynamometer could provide further insight in the specific tension of muscle groups and overcome this limitation.

Estimated Fmax were overall significantly smaller when based on VM than when linearly scaled to lower-limb mass (Fig 6), except for the Adductor magnus, Vastus intermedius and medialis, Gastrocnemius medialis, and Soleus, that, on the contrary, presented significantly higher values for the LLMS approach. Declining muscle strength has been observed from the age of fifty [3] and a reduction by 20% of Fmax has been quantified in older people aged seventy [41]. This could explain the smaller Fmax obtained from individual VM (when volume loss associated with ageing was taken into account) compared to a scaling approach. This also confirms previous literature suggesting that a scaling approach might only be appropriate if starting from values from a sex- and age-matched population [28].

The choice of Fmax highly impacts the output of MSK models [42,43], since a change in an individual muscle ability to produce force alters the solution of the static optimisation problem [42], affecting both individual muscle force estimates and the resulting joint contact force. A previous study found limited sensitivity of muscle forces and joint contact forces to Fmax [32] estimating its values based on scaling of literature values or using Handsfield’s regression equations. Ackland et al. [19] studied the effect of variation between +10% and −10% of Fmax nominal value, reporting no significant changes in the model output. However, in their study, they did not account for actual muscle geometry to estimate Fmax, which proved to cause variation up to 400% in our study when compared to scaling approaches. This suggests that calculating individual Fmax from MRI-segmentations could affect the estimates of muscle forces and joint contact forces on a larger scale than reported in the literature and lead to more accurate estimates. This supports the conclusions from Arnold et al. [16] that tuning individual muscle parameters might provide estimates of internal forces that compare better to experimental measurements [16]. Further studies are needed to confirm this hypothesis.

This study had some limitations. Out of the 35 muscles commonly included in lower-limb MSK models, only 23 were included in the study, as these were not significantly affected by operator-related error in the segmentation. Muscle segmentation is a time-consuming (10 hours per subject on average for this study) and operator-dependent procedure, therefore further effort should be put into developing automated algorithms based on machine learning [28] for the segmentation of individual muscles or statistical shape modelling-based approaches for the extraction of muscle volume and muscle centreline/length. The dataset associated with this paper is publicly available, which will likely foster advances in this field, i.e. acting as a reference atlas.

The cohort enrolled for this study included eleven participants; a larger sample size would be needed to ensure generalisability of the results observed here. Our results suggest that muscle asymmetry could be higher in older adults due to age-related processes. However, this finding is based on comparison to literature [24], where data were obtained following a slightly different methodology. Therefore, a wider study, including a control group of younger women, should be designed to prove our hypothesis.

In the attempt of preserving a degree of subject-specificity in the muscle parameters, PSCA was calculated from muscle volume and length. Nonetheless, due to the impossibility of estimating the optimal fibre length from the implemented MRI sequence, the required ratio between optimal fibre length and muscle length was taken from cadaveric data. Diffusion Tensor imaging recently proved to be a valuable option to enable both muscle segmentation and the estimate of fibre length [27]. Further studies are needed to understand whether this technique might be included within an MSK modelling imaging protocol to overcome this limitation.

In conclusion, this study uniquely proved the existence of significantly large muscle- and subject-specific asymmetry in muscle volume, length, and PCSA. This suggests that individual differences in muscle geometry must not be neglected, and inter-limb symmetry cannot be assumed in older women. Personalised muscle characteristics should be accounted for in MSK models aiming at investigating dynamic tasks such as walking, where strength asymmetry plays an important role in older women. This could be of substantial relevance when internal forces are used in clinical contexts, such as prediction of osteoporotic risk of fracture.

Supporting information

S2 Table. Right-limb muscle volumes segmented by three operators for three randomly selected subjects.

Maximum coefficient of variation (CoV) across the three datasets is reported.

https://doi.org/10.1371/journal.pone.0242973.s002

(DOCX)

S3 Table. Right-limb muscle volumes segmented three times by one operator.

Coefficient of variation (CoV) across the three repetitions is reported.

https://doi.org/10.1371/journal.pone.0242973.s003

(DOCX)

S4 Table. Right and left volume of the muscles segmented in the lower limbs of the eleven subjects enrolled in the study Mean, standard deviation (SD) and coefficient of variation (CoV) are reported.

https://doi.org/10.1371/journal.pone.0242973.s004

(DOCX)

S5 Table. Right and left length of the muscles segmented in the lower limbs of the eleven subjects enrolled in the study Mean, standard deviation (SD) and coefficient of variation (CoV) are reported.

https://doi.org/10.1371/journal.pone.0242973.s005

(DOCX)

S6 Table. Physiological cross-sectional areas (PCSAs) measured for the eleven subjects enrolled in our study (for right and left muscles) and three cadavers included in Charles et al., 2019.

Mean and mean and SD PCSA are reported for Ward et al. 2009 and Handsfield et al. 2014, respectively.

https://doi.org/10.1371/journal.pone.0242973.s006

(DOCX)

Acknowledgments

The authors would like to acknowledge Dr Geoffrey Handsfield for sharing relevant data from [24] and Dr Enrico Dall’Ara for his valuable input around image processing.

We are particularly grateful to the participants who volunteered for the study.

References

  1. 1. Cruz-Jentoft AJ, Baeyens JP, Bauer JM, Boirie Y, Cederholm T, Landi F, et al. Sarcopenia: European consensus on definition and diagnosis. Report of the European working group on sarcopenia in older people. Age and Ageing. 2010;39(4):412–23. pmid:20392703
  2. 2. Häkkinen K, Häkkinen A. Muscle cross-sectional area, force production and relaxation characteristics in women at different ages. European Journal of Applied Physiology and Occupational Physiology. 1991;62(6):410–4. pmid:1893904
  3. 3. Lindle R, Metter E, Lynch N, Fleg J, Fozard J, Tobin J, et al. Age and gender comparisons of muscle strength in 654 women and men aged 20–93 yr. Journal of applied physiology. 1997;83(5):1581–7. pmid:9375323
  4. 4. Rahemi H, Nigam N, Wakeling JM. The effect of intramuscular fat on skeletal muscle mechanics: implications for the elderly and obese. Journal of The Royal Society Interface. 2015;12(109):20150365. pmid:26156300
  5. 5. Larsson L, Grimby G, Karlsson J. Muscle strength and speed of movement in relation to age and muscle morphology. Journal of Applied Physiology. 1979;46(3):451–6. pmid:438011
  6. 6. Yoshiko A, Hioki M, Kanehira N, Shimaoka K, Koike T, Sakakibara H, et al. Three-dimensional comparison of intramuscular fat content between young and old adults. BMC medical imaging. 2017;17(1):1–8. pmid:28056868
  7. 7. Jonkers I, Sauwen N, Lenaerts G, Mulier M, Van der Perre G, Jaecques S. Relation between subject-specific hip joint loading, stress distribution in the proximal femur and bone mineral density changes after total hip replacement. Journal of Biomechanics. 2008;41(16):3405–13. pmid:19019372
  8. 8. Steele KM, Seth A, Hicks JL, Schwartz MS, Delp SL. Muscle contributions to support and progression during single-limb stance in crouch gait. Journal of Biomechanics. 2010;43(11):2099–105. pmid:20493489
  9. 9. Liu MQ, Anderson FC, Schwartz MH, Delp SL. Muscle contributions to support and progression over a range of walking speeds. Journal of Biomechanics. 2008;41(15):3243–52. pmid:18822415
  10. 10. Viceconti M, Taddei F, Cristofolini L, Martelli S, Falcinelli C, Schileo E. Are spontaneous fractures possible? An example of clinical application for personalised, multiscale neuro-musculo-skeletal modelling. Journal of biomechanics. 2012;45(3):421–6. pmid:22204893
  11. 11. Bahl JS, Zhang J, Killen BA, Taylor M, Solomon LB, Arnold JB, et al. Statistical shape modelling versus linear scaling: effects on predictions of hip joint centre location and muscle moment arms in people with hip osteoarthritis. Journal of Biomechanics. 2019. pmid:30770197
  12. 12. Martelli S, Valente G, Viceconti M, Taddei F. Sensitivity of a subject-specific musculoskeletal model to the uncertainties on the joint axes location. Computer Methods in Biomechanics and Biomedical Engineering. 2015;18(14):1555–63. pmid:24963785
  13. 13. Valente G, Pitto L, Testi D, Seth A, Delp SL, Stagni R, et al. Are subject-specific musculoskeletal models robust to the uncertainties in parameter identification? PLoS One. 2014;9(11):e112625. pmid:25390896
  14. 14. Wesseling M, De Groote F, Bosmans L, Bartels W, Meyer C, Desloovere K, et al. Subject-specific geometrical detail rather than cost function formulation affects hip loading calculation. Computer Methods in Biomechanics and Biomedical Engineering. 2016;19(14):1475–88. pmid:26930478
  15. 15. Gerus P, Sartori M, Besier TF, Fregly BJ, Delp SL, Banks SA, et al. Subject-specific knee joint geometry improves predictions of medial tibiofemoral contact forces. Journal of Biomechanics. 2013;46(16):2778–86. pmid:24074941
  16. 16. Arnold EM, Ward SR, Lieber RL, Delp SL. A model of the lower limb for analysis of human movement. Annals of Biomedical Engineering. 2010;38(2):269–79. pmid:19957039
  17. 17. Delp SL, Loan JP, Hoy MG, Zajac FE, Topp EL, Rosen JM. An interactive graphics-based model of the lower extremity to study orthopaedic surgical procedures. IEEE Transactions On Biomedical Engineering. 1990;37(8):757–67. pmid:2210784
  18. 18. Bujalski P, Martins J, Stirling L. A Monte Carlo analysis of muscle force estimation sensitivity to muscle-tendon properties using a Hill-based muscle model. Journal of Biomechanics. 2018;79:67–77. pmid:30146173
  19. 19. Ackland DC, Lin Y-C, Pandy MG. Sensitivity of model predictions of muscle function to changes in moment arms and muscle–tendon properties: a Monte-Carlo analysis. Journal of Biomechanics. 2012;45(8):1463–71. pmid:22507351
  20. 20. Buchanan TS, Lloyd DG, Manal K, Besier TF. Neuromusculoskeletal modeling: estimation of muscle forces and joint moments and movements from measurements of neural command. Journal of Applied Biomechanics. 2004;20(4):367–95. pmid:16467928
  21. 21. Brooks SV, Faulkner JA. Skeletal muscle weakness in old age: underlying mechanisms. Medicine and science in sports and exercise. 1994;26(4):432–9. pmid:8201898
  22. 22. Maden-Wilkinson T, Degens H, Jones D, McPhee J. Comparison of MRI and DXA to measure muscle size and age-related atrophy in thigh muscles. Journal of Musculoskeletal Neuronal Interactions. 2013;13(3):320–8. pmid:23989253
  23. 23. Kainz H, Goudriaan M, Falisse A, Huenaerts C, Desloovere K, De Groote F, et al. The influence of maximum isometric muscle force scaling on estimated muscle forces from musculoskeletal models of children with cerebral palsy. Gait & Posture. 2018;65:213–20. pmid:30558934
  24. 24. Handsfield GG, Meyer CH, Hart JM, Abel MF, Blemker SS. Relationships of 35 lower limb muscles to height and body mass quantified using MRI. Journal of Biomechanical Engineering. 2014;47(3):631–8.
  25. 25. Ward SR, Eng CM, Smallwood LH, Lieber RL. Are current measurements of lower extremity muscle architecture accurate? Clinical Orthopaedics and Related Research. 2009;467(4):1074–82. pmid:18972175
  26. 26. Fick R. Handbuch der anatomie und mechanik der gelenke: Рипол Классик; 2012.
  27. 27. Charles JP, Moon C-H, Anderst WJ. Determining subject-specific lower-limb muscle architecture data for musculoskeletal models using diffusion tensor imaging. Journal of Biomechanical Engineering. 2019;141(6).
  28. 28. Yeung S, Fernandez J, Handsfield G, Walker C, Besier T, Zhang J. Rapid muscle volume prediction using anthropometric measurements and population-derived statistical models. Biomechanics and modeling in mechanobiology. 2019:1–11. pmid:31667655
  29. 29. Greenberger HB, Paterno MV. Relationship of knee extensor strength and hopping test performance in the assessment of lower extremity function. Journal of Orthopaedic & Sports Physical Therapy. 1995;22(5):202–6. pmid:8580947
  30. 30. Campbell MJ, Machin D, Walters SJ. Medical statistics: a textbook for the health sciences: John Wiley & Sons; 2010.
  31. 31. Wang C-Y, Olson SL, Protas EJ. Test-retest strength reliability: hand-held dynamometry in community-dwelling elderly fallers. Archives of Physical Medicine and Rehabilitation. 2002;83(6):811–5. pmid:12048660
  32. 32. Modenese L, Montefiori E, Wang A, Wesarg S, Viceconti M, Mazzà C. Investigation of the dependence of joint contact forces on musculotendon parameters using a codified workflow for image-based modelling. Journal of Biomechanics. 2018;73:108–18. pmid:29673935
  33. 33. White D, Woodard H, Hammond S. Average soft-tissue and bone models for use in radiation dosimetry. The British Journal Of Radiology. 1987;60(717):907–13. pmid:3664185
  34. 34. Delp SL, Anderson FC, Arnold AS, Loan P, Habib A, John CT, et al. OpenSim: open-source software to create and analyze dynamic simulations of movement. IEEE Transactions On Biomedical Engineering. 2007;54(11):1940–50. pmid:18018689
  35. 35. Wickiewicz TL, Roy RR, Powell PL, Edgerton VR. Muscle architecture of the human lower limb. Clinical Orthopaedics and Related Research. 1983;(179):275–83. pmid:6617027
  36. 36. Ward SR, Lieber RL. Density and hydration of fresh and fixed human skeletal muscle. Journal of Biomechanics. 2005;38(11):2317–20. pmid:16154420
  37. 37. Rankin G, Stokes M, Newham DJ. Abdominal muscle size and symmetry in normal subjects. Muscle & Nerve: Official Journal of the American Association of Electrodiagnostic Medicine. 2006;34(3):320–6. pmid:16775833
  38. 38. Stock MS, Oranchuk DJ, Burton AM, Phan DC. Age, sex, and region-specific differences in skeletal muscle size and quality. Applied Physiology, Nutrition, and Metabolism. 2020;(ja).
  39. 39. Montefiori E, Modenese L, Di Marco R, Magni-Manzoni S, Malattia C, Petrarca M, et al. Linking joint impairment and gait biomechanics in patients with juvenile idiopathic arthritis. Annals of Biomedical Engineering. 2019:1–13. pmid:31111329
  40. 40. Dziedzic D, Bogacka U, Ciszek B. Anatomy of sartorius muscle. Folia Morphologica. 2014;73(3):359–62. pmid:25242250
  41. 41. Thelen DG. Adjustment of muscle mechanics model parameters to simulate dynamic contractions in older adults. Journal of Biomechanical Engineering. 2003;125(1):70–7. pmid:12661198
  42. 42. Carbone V, van der Krogt M, Koopman HF, Verdonschot N. Sensitivity of subject-specific models to Hill muscle–tendon model parameters in simulations of gait. Journal of Biomechanics. 2016;49(9):1953–60. pmid:27131851
  43. 43. De Groote F, Van Campen A, Jonkers I, De Schutter J. Sensitivity of dynamic simulations of gait and dynamometer experiments to hill muscle model parameters of knee flexors and extensors. Journal of Biomechanics. 2010;43(10):1876–83. pmid:20392450