Next Article in Journal
Knowledge, Attitude, and Practice towards COVID-19 among Patients Attending Phuentsholing Hospital, Bhutan: A Cross-Sectional Study
Previous Article in Journal
Volumetric Assessment of Apical Periodontitis Using Cone-Beam Computed Tomography—A Systematic Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling the Dynamics of Children’s Musculoskeletal Fitness

1
Instituto Superior Manuel Teixeira Gomes (ISMAT), 8500-724 Portimão, Portugal
2
Centre of Research, Education, Innovation and Intervention in Sport (CIFI2D), Faculty of Sport, University of Porto, 4099-002 Porto, Portugal
3
Academic Department of Physical Education, Federal University of Technology of Paraná, Curitiba 80230-901, Brazil
4
Centro de Investigação em Desporto, Educação Física, Exercício e Saúde (CIDEFES), Lusófona University, 1749-024 Lisboa, Portugal
5
School of Physical Education and Sport, University of São Paulo, São Paulo 13560-970, Brazil
6
Physical Education, University of South Carolina, Columbia, SC 29208, USA
7
Department of Public Health Sciences, University of Chicago, Chicago, IL 60637, USA
8
College of Kinesiology, University of Saskatchewan, Saskatoon, SK S7N 5B5, Canada
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2023, 20(4), 2938; https://doi.org/10.3390/ijerph20042938
Submission received: 13 December 2022 / Revised: 4 February 2023 / Accepted: 6 February 2023 / Published: 8 February 2023

Abstract

:
This study models children’s musculoskeletal fitness (MSF) developmental trajectories and identifies individual differences related to effects of time-invariant, as well as time-varying covariates. A total of 348 Portuguese children (177 girls) from six age cohorts were followed for three years. MSF tests (handgrip strength, standing long jump and shuttle run), age, body mass index (BMI), socioeconomic status (SES), gross motor coordination (GMC) and physical activity (PA) were assessed. Data were analyzed using multilevel models. Between 5 and 11 years of age, boys outperformed girls in all three MSF tests (p < 0.05). Birth weight was positively associated with shuttle run performance (β = −0.18 ± 0.09, p < 0.05). BMI was positively associated with handgrip strength (β = 0.35 ± 0.04, p < 0.001) and shuttle run performance (β = 0.06 ± 0.01, p < 0.001), but negatively associated with standing long jump performance (β = −0.93 ± 0.23, p < 0.001). GMC was positively associated (p < 0.001) with all three MSF tests, while PA was associated with standing long jump (β = 0.08 ± 0.02, p < 0.05) and shuttle run (β = −0.003 ± 0.002, p < 0.05) performance only. No school environment effects were found, and SES was not related to any MSF tests. Children’s MSF development showed a curvilinear shape with increasing age, with boys outperforming girls. Weight status and physical behavior characteristics predicted MSF development, while environmental variables did not. Examining potential longitudinal predictors of MSF across multiple dimensions is important to gain a more holistic understanding of children’s physical development as well as to future interventions.

1. Introduction

Musculoskeletal fitness (MSF), a multidimensional construct comprising muscle strength, muscle endurance, and muscle power, is an important marker of health status [1]. Differences in children’s MSF are explained by a variety of genetic and environmental factors [2]. Cross-sectional studies of children’s MSF have shown positive associations with factors such as physical activity (PA), birth weight, motor competence, gross motor coordination (GMC), and socioeconomic status (SES) [3,4]. School environments also influence MSF [5]. However, few studies have investigated how all these factors combined impact MSF development.
Longitudinal studies of human growth indicate the pattern of curvilinear growth is consistent, but the timing and magnitude of growth is highly individualized [6,7]. This is also true for MSF outcomes [2]. To adequately distinguish the independent effects of sex, birth weight, SES, school environment, body mass index, GMC, and PA on MSF development, time-invariant and time-varying factors are required; this can only be achieved using longitudinal data and appropriate statistical analysis. This approach accounts for the wide variation shown amongst children’s growth parameters at any given age, and in the velocity of these parameters across age, and requires hierarchical (or multilevel) modeling [8].
A longitudinal study by Ruedl et al. [9] reported that physical fitness development in overweight boys and girls was lower than their non-overweight peers. In addition, a study by Rodrigues et al. [10], describing different velocities, found that children with low or average rates of change in MSF and motor competence were more likely to become overweight or obese by the end of childhood, regardless of their sex or initial body mass index (BMI). Also, Haugen and Johansen [11] found that the initially high gross motor competence-grouped children performed significantly better on a multidimensional physical fitness test, compared to the initially low gross motor competence-grouped children, at all time-points. Additionally, Bai et al. [12] reported that higher SES and school size were associated with better MSF profiles. One limitation of these studies is that they did not consider other potential MSF correlates, nor did they consider the hierarchical data structure [9] of the school context [10], or the time dependent effects of PA or GMC [12].
The growth, motor development and cognition study (GMDC) [11] was initiated to address these limitations. The objectives of this paper were: (i) to model intra-individual MSF developmental trajectories, (ii) identify the magnitude of inter-individual differences and their dependence on school contexts, and (iii) investigate the effects of time-invariant (e.g., sex, birth weight and SES) and time-varying (e.g., age, BMI, GMC and PA) covariates on children’s MSF development. We hypothesized that MFS development would show a curvilinear trend (positive for hand grip strength and for standing long jump, but negative for shuttle run), and the trajectories would be greater in boys. Furthermore, we hypothesized that school contexts would influence trajectories, and that time-variant and invariant predictors would have independent significant effects on MFS development.

2. Materials and Methods

The GMDC was carried out in Vouzela, a central region of Portugal, between 2013 and 2016. The mixed-longitudinal design has been described in detail elsewhere [13]. In brief, data were collected in six age cohorts, followed for three consecutive years (4–6, 5–7, 6–8, 7–9, 8–10 and 9–11 years). All children from the 19 schools of the region were invited (participation rate ~90%). Children with special needs (n = 24) were excluded from the analysis as were 4-year-old’s (n = 88) whose GMC was not assessed. Written informed consent was obtained from parents or legal guardians, and the project was approved by local authorities and the ethics committee of the Faculty of Sport, University of Porto (CEFADE 01.2016).
Data from all cohorts (n = 348, 177 girls) with complete data were used in the analysis (Table 1). In addition, differences in a set of variables (sex, age, PF tests, birth weight, BMI, GMC, PA, and SES) were analyzed to identify putative differences between those with complete and those with missing data (n = 150, 75 girls). Significant differences were found in BMI, shuttle-run test and SES, favoring those included in the analysis (p < 0.05).
Information on the school context was obtained via a questionnaire developed by the authors in conjunction with the Vouzela city hall Education Department, namely: (i) school characterization (location; size and setting—rural, semi-urban and urban as determined by the Portuguese National Statistics Institute); (ii) policies and practices to stimulate physical activity and adequate nutrition; (iii) physical infrastructures of the school (playground, multi-sports roofed facilities and equipment for physical education); (iv) physical education (PE) classes; and (v) human resources (number of teachers in each school, academic degree of teachers, and teaching time).
SES was identified using the school social support classification system (Portuguese Ministry of Education) [14]. Three SES groupings were identified: (i) level A: up to EUR 2.934 per year−1, children received books and school lunches; (ii) level B: from EUR 2.934 to EUR 5.896 per year−1, with half of level A support; (iii) level C: ≥EUR 5.897 per year−1 implies no support.
Height (cm) was measured, without shoes, using a portable stadiometer (Holtain, UK) holding the child’s head in the Frankfurt plane, (accuracy = 0.1 cm). Body mass (kg) was measured with a portable bioelectrical impedance scale (TANITA BC-418 MA Segmental Body Composition Analyzer; Tanita, Corporation, Japan) with children wearing light clothing and without shoes (precision = 0.1 kg). Body mass index (BMI) was calculated using the standard formula [BMI = body weight (kg)/height (m)2]. Birth weight (kg) was obtained retrospectively from children’s health booklets.
MSF (isometric strength and explosive strength/agility) [1] was assessed by: (i) Grip strength (HG in kg) via hand held dynamometry (Takei Physical Fitness Test GRIP-D, Japan). In all measurements the dynamometer was well fitted to the hand of the child, making the child as comfortable as possible. The most comfortable and stable positions were adopted to achieve optimal conditions for maximum muscle strength efforts. The child stood with the dynamometer beside their body, without touching the body. The orientation of each desired action was explained to the child while the investigator supported the movement. Verbal encouragement was given during the test. The dominant hand was previously established, as the children were evaluated with respect to their laterality [13]. We report the mean of the two trials performed with the dominant hand. The time set between the two trials was 2 min; (ii) Standing long jump (SLJ in cm) assessed lower body explosive strength via maximum jumping distance. The greatest distance from two trials was recorded; (iii) A shuttle run (SR in s) to measure explosive agility. Children ran from a starting line to a line of wooden blocks nine meters away, picked up a block, returned and placed the block on the starting line, before repeating the procedure. The quickest time of two trials, assessed by chronometry (Géonaute ON START 300), was recorded. The order of the testing was as follows: (i) anthropometric measurements, (ii) handgrip strength; (iii) SLJ and SH were assessed in either a multi-sports roofed or playground area of the schools. Tests were performed after a 10-min break between tests for the children to rest. All children participated in a warm-up period of 5 min, conducted by a member of the research team.
The Köperkoordinationstest für Kinder (KTK) test battery [15] was used to assess GMC: (i) walking backwards—children walked backwards on three different balance beams, decreasing in width, three times. The number of successful steps was counted and recorded; (ii) hopping—children hop, with one leg, over an increasing number of foam squares, and successful completions were recorded (three attempts are allowed at each height); (ii) jumping sideways—with feet together, children jumped sideways over a wooden slat for 15 s. This was repeated twice, and the number of jumps was summed; (iv) moving sideways—children stood with both feet on one platform and had to place both hands on an adjacent platform, thus passing to the next platform. The total number of relocations over 20′s in two trials was recorded and summed. The sum of scores from four tests was used as a measure of total GMC (GMCTS) [15]. The KTK battery tests were arranged in the circuit manner in the order described above, with a 5-min break between each test.
PA was objectively measured using an Actigraph GT3X+ accelerometer (Actigraph, Pensacola, FL, USA) over 7 consecutive days (5 weekdays and 2 week-end days). Children were instructed to only remove the device during water activities and while sleeping at night. The accelerometer was placed on the iliac crest and held in place by an elastic belt with an adjustable clip. Each child received individual instructions and all parents received written instructions. ActiLife® software 6.5.4 was used to process the recorded data immediately upon retrieval of each accelerometer. Children’s data, to be considered valid, had to have a minimum of at least 4 days (with at least one weekend day) and a minimum of 10 h of daily wear time. Any sequence of at least 20 consecutive minutes of zero activity counts was considered non-wear time [16]. Moderate-to-vigorous PA (MVPA) was derived using cut-points developed by Evenson et al. [17], defined as all activities greater than 574 counts/15 s epochs, and expressed in min∙day−1.
All measures were taken by trained team members of the Kinanthropometry and Motor Learning Labs (Faculty of Sport, University of Porto). Each year, a sample of ~60 children was randomly selected from each cohort and retested. Technical error of measurement (TEM) and ANOVA-based intraclass correlation (R) were used for reliability estimates: TEM = 0.17 cm in height, and 0.10 kg in weight; R-values ranged from 0.71 to 0.97 for PF tests, and from 0.78 to 0.91 in GMC tests. Finally, a systematic check was performed regarding errors in data entry and normality of variables’ distributions.
Descriptive statistics by sex across age were calculated. A multilevel modeling approach was used [18]. As advocated [19], the time metric was anchored at 5 years of age (time 0), to ensure intercepts occurred within the data, such that age centers of 0, 1, 2, 3, 4, 5, and 6 years, correspond to 5, 6, 7, 8, 9, 10 and 11 years of age. All analyses were performed in SuperMix and model fit was simultaneously estimated by maximum likelihood [19]. Models were tested sequentially: Model 1 described the best developmental trajectory for each MSF test using age, age2, sex, age-by-sex and age-by-sex2 interactions as predictors; parameters were dropped from the final models if they were non-significant (p > 0.05); Model 2 included the cohort effect, non-varying (sex, birth weight, SES and school environment), and time-varying (BMI, GMC, MVPA) covariates. Individual models were developed for each of the three MSF markers (grip strength, standing long jump, and shuttle run). The significance level was set at p < 0.05.

3. Results

Boys’ and girls’ descriptive statistics (mean ± SD or percentages) are shown in Table 2; PF, BMI, and GMCTS increased with increasing age. No differences in birth weights were found between sexes. Girls’ MVPA decreased until 9 years and increased thereafter; in boys this decrease was seen until 10 years and then increased. SES level C (≥ EUR 5.870 per year−1) was the most frequent category noted.
The number of students per school varied, from a minimum of 5 to a maximum of 87 (Table 3). Most schools were located in rural areas and had policies and practices for PA. All schools had playgrounds without obstacles (any obstacle that could interfere with the children’s playground space) and their dimensions were generally greater than 70 m2. Only five schools had a sport center, three of them with more than 50 m2. All schools had equipment for PA; however, fifteen schools had only one infrastructure space for PA. The duration of PE classes ranged from 45 to 60 min and the time children spent active varied from 30 to 50 min. The same PE opportunities were offered to all children since schools employed PE specialists.
In model 1, no significant effects were identified for the interactions of age-by-sex or age2-by-sex, or for the school random effects. These factors were not included in the second model (Table 4). Regarding MSF tests, the averages for girls at 5 years (the centered age) were: HG = 7.74 kg, SLJ = 88.95 cm and SR = 15.49 s. Boys significantly outperformed girls at this age (+1.13 kg in HG, +7.88 cm in SLJ, and −0.33 s in SR) and this effect did not vary with age. The HG instantaneous velocity at 5 years was 1.43 ± 0.19 kg and the average acceleration (Age2) was also statistically significant (β = 0.06 ± 0.03, p < 0.05). A similar trend was observed for SLJ (velocity, β = 12.29 ± 1.20, p < 0.001; acceleration, β = 1.08 ± 0.17, p < 0.001), and SR (velocity, β = 0.93 ± 0.09 s, p < 0.001; acceleration, β = 0.07 ± 0.01, p < 0.001). Cohort effects were all significant (p < 0.05) in all MSF markers. Birth weight was only significantly related to SR, i.e., children born with higher birth weight demonstrated, on average, better agility during childhood (β = −0.18 ± 0.09, p < 0.05), i.e., faster SR test time. Further, SES was not significantly linked to any PF developmental trajectory. In contrast, children with higher BMI were stronger in HG (β = −0.35 ± 0.04, p < 0.001); less agile (β = 0.06 ± 0.01, p < 0.001), i.e., slower SR test time, and showed lower explosive strength (β = −0.93 ± 0.23, p < 0.001). Children with higher GMCTS demonstrated better PF developmental trajectories (HG, β = 0.02 ± 0.003, p < 0.001; SLJ, β = 0.23 ± 0.02, p < 0.001; SR, β = −0.02 ± 0.001, p < 0.001).
Additionally, more physically active children displayed a higher rate of increase in SLJ and SR, but not in the rate of HG development. Finally, children’s MSF developmental trajectories showed significant inter-individual differences, i.e., there was significant variation (σ2) in the longitudinal age-related trends across individuals (HG, σ2 = 0.21 ± 0.10, p < 0.05; SLJ, σ2 = 9.39 ± 3.42, p < 0.05; SR, σ2 = 0.06 ± 0.02, p < 0.001), meaning that children’s MSF developmental trajectories were significantly heterogeneous.

4. Discussion

The results show that isometric and explosive strength, as well as explosive agility, exhibit a curvilinear developmental trajectory (Age2 was significant) during childhood, which are analogous to other longitudinal studies [20]. When controlling for various covariates (birth weight, SES, BMI, GMCTS, and MVPA), boys outperformed girls in all PF developmental trajectories, corroborating other studies’ results [21]. Such sex-differences have been equated with biological and sociocultural aspects of development [2]. For example, it has been suggested that parents are more prone to let boys be involved in more intense rough play, whereas girls are treated as being more fragile [22] and participated more often and in more skill-based activities [23]. Boys also tend to demonstrate higher performance in GMC tasks [2] which is linked to the development of advanced neuromuscular function (e.g., inter- and intra-muscular coordination, motor unit recruitment/firing rate, proprioceptor sensitivity) that also promotes the development of strength and explosive strength/agility (i.e., musculoskeletal fitness) [24].
No relationship was found between children’s SES and PF trajectories, which is in contrast to other available data [12] that reported that school children attending schools with higher SES showed better profiles in aerobic capacity and BMI. This is probably related to the fact that Vouzela children, regardless of their SES inequality, generally have the same opportunities to develop their fundamental motor skills. Consequently, their MSF may not be related to SES because all children have mandatory physical education classes and equal access to city hall sporting facilities.
BMI increases as children age, and the net expression of inter-individual differences as they increase in height and weight reflect factors linked to motor performance and PF levels [8]. We show that children’s changes in static strength and agility are positively related to increases in BMI; in contrast, SLJ is negatively associated. These results may be interpreted as follows: (i) with increasing age, children static strength increases and is related to improved neuromotor function and increases in absolute muscle mass [12]; (ii) agility development is linked to improved neuromotor function, GMC and more effective motor control, which all contribute to the effective acceleration and deceleration of the center of mass noted in repeated change of direction tasks, which was better expressed in children with lower BMI values; and (iii) although there is an assumption that vertical jumps are theoretically independent of body dimensions [25], and SLJ performance may be proportional to height [26], empirical data contradict these suggestions [25,26]. In contrast, force production is proportional to muscle physiological cross-sectional area to the power of 2, and body size is expected to be proportional to the power of 3, which means that greater body mass relative to height may impair the SLJ performance; thus, potentially explaining the negative relationship found between increases in BMI and decreased SLJ distance.
Physiological mechanisms have previously linked birth weight with fat-free mass and muscle fiber development, and birth weight has been linked with physical performance [27]. In the present study we found no association of birth weight with static or explosive strength, or a positive relationship with agility. Robič Pikel et al. [28] found that although low birth weight children underperformed in explosive leg strength and speed tests, their agility was linked to GMC. Overall, the mechanisms of why specific relationships are found speaks to the importance of understanding both biological and experiential mechanisms for various aspects of physical performance and its development across time.
Children’s MVPA was not associated with HG strength but was positively interrelated to SLJ and SR performance. Although the longitudinal links between children’s PA and MSF trajectories are rarely investigated, Augste et al. [29] showed that PA affected not only children’s composite physical fitness baseline levels, but also their trajectory trends. Further, Breu et al. [30] indicated that children’s physical fitness was associated with their MVPA, and that they met the WHO guidelines at follow-up. Children’s participation in a variety of locomotor physical activities, specifically in early childhood, provides a logical link to tasks linked to lower extremity fitness and GMC performance tests. Finally, we speculate that since HG is specifically linked to static strength, which in turn is linked with size, it may not be as strongly linked to children’s general PA levels.
It has been shown in primary school children that motor competence is positively related to physical fitness, and that change in motor competence development is linked to overweight or obesity [10]. We found that more coordinated children displayed better MSF developmental trajectories. Motor tasks in the GMC and MSF tests include multiplanar and multiaxial movements, and the combinations of different muscle contractions, integrated across different phases of movements, require a high degree of inter- and intramuscular coordination and control [24]. Similarly, the structure of explosive MSF tasks requires high levels of multijoin coordination and control inherent within GMC tasks.
Finally, we found that the school environment did not significantly explain differences in MSF development, probably because there was a high degree of similarity across schools. For example, the school-based curriculum in all Portuguese primary schools is applied in the same manner and rural areas have similar PA/sports facilities and physical education specialists. Further, school effects remain non-significant even when cohort effects were modelled and tested [31], although it is possible that children’s lives and educational histories within each age cohort could have impacted their MSF developmental trends.
In summary, our results illustrate the importance of understanding the relationships between various aspects of physical growth, SES and physical activity on MSF developmental during mid-childhood—an important time window that can be considered as a “sensitive” and transient period marked by a higher degree of plasticity in neuromuscular function. Our results may also increase physical educators’ and coaches’ awareness of these windows of opportunity for children’s MSF development within a life-long perspective. Even parents themselves need to always consider the links between physical growth and maturation that shape children’s MSF development. Ideally, it is important not only to focus on MSF outcomes, but to understand “how,” “why,” and “when” children are able, or not, to perform MSF tasks in terms of standardized movements and skills. In this way, such professionals and parents will be better equipped to provide appropriate learning environments, and conditions, linked to children’s intra- and inter-individual differences, allowing for more appropriate feedback that can impact children’s healthy trajectories.
A few limitations are acknowledged: (i) the GMDC-Vouzela study does not cover all Portuguese regions, and generalizability is questioned; (ii) a pure longitudinal study would have provided a more robust design. However, comprehensive testing of multiple intra-individual and environmental factors linked to physical development has never been conducted and the mixed-longitudinal approach is a highly viable and reliable design.

5. Conclusions

In conclusion, children’s MSF development exhibited a curvilinear trend with boys outperforming girls. The complex nature of interactions among MSF attributes, GMC and bodyweight status are difficult to untangle, specifically with the additional influence of biological growth and maturation variables. However, the unique nature of this investigation provides new novel insight on the overall physical development of children, supplying parents, physical education teachers and coaches with a more holistic view of children´s MSF development for future interventions.

Author Contributions

Conceptualization, J.M.; methodology, A.R., R.C., J.M. and A.B.-J.; formal analysis, D.H. and J.M.; investigation, A.R., R.C., O.V., S.P. and G.T.; data curation, J.M.; writing—original draft preparation, A.R., J.M. and A.B.-J.; writing—review and editing, R.C., O.V., S.P., G.T. and D.S.; project administration, J.M.; All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki and approved by the ethics committee of the Faculty of Sport, University of Porto, approved the research project (CEFADE 01.2016).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data used in this study can be requested from one of the authors (J.M.) who is the curator of the broad data set of this study.

Acknowledgments

We thank the Vouzela City Hall, School authorities, Health Center nurses, psychologists, school teachers, children and parents for their support in participating in the study across the three years.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Pate, R.; Oria, M.; Pillsbury, L. Fitness Measures and Health Outcomes in Youth. In Fitness Measures and Health Outcomes in Youth; Pate, R., Oria, M., Pillsbury, L., Eds.; National Academies Press (US): Washington, DC, USA, 2012. [Google Scholar]
  2. Malina, R.M.; Bouchard, C.; Bar-Or, O. Growth, Maturation, and Physical Activity, 2nd ed.; Human Kinetics: Champaign, IL, USA, 2004. [Google Scholar]
  3. Kaioglou, V.; Dania, A.; Kambas, A.; Venetsanou, F. Associations of Motor Competence, Cardiorespiratory Fitness, and Physical Activity. The Mediating Role of Cardiorrespiratory Fitness. Res. Q. Exerc. Sport 2022, 30, 1–7. [Google Scholar] [CrossRef] [PubMed]
  4. Utesch, T.; Bardid, F.; Büsch, D.; Strauss, B. The Relationship between Motor Competence and Physical Fitness from Early Childhood to Early Adulthood: A Meta-Analysis. Sports Med. 2019, 49, 541–551. [Google Scholar] [CrossRef]
  5. Lo, K.-Y.; Wu, M.-C.; Tung, S.-C.; Hsieh, C.C.; Yao, H.-H.; Ho, C.-C. Association of School Environment and After-School Physical Activity with Health-Related Physical Fitness among Junior High School Students in Taiwan. Int. J. Environ. Res. Public Health 2017, 14, 83. [Google Scholar] [CrossRef] [PubMed]
  6. Roche, A.; Sun, S. Human Growth: Assessment and Interpretation; Cambridge University Press: Cambridge, UK, 2003. [Google Scholar]
  7. Fraser, B.J.; Blizzard, L.; Buscot, M.J.; Schmidt, M.D.; Dwyer, T.; Venn, A.J.; Magnussen, C.G. Muscular strength across the life course: The tracking and trajectory patterns of muscular strength between childhood and mid-adulthood in an Australian cohort. J. Sci. Med. Sport 2021, 24, 696–701. [Google Scholar] [CrossRef] [PubMed]
  8. Goldstein, H. Multilevel Statistical Models, 2nd ed.; Arnold: London, UK, 1995. [Google Scholar]
  9. Ruedl, G.; Niedermeier, M.; Posch, M.; Kirschner, W.; Wirnitzer, K.; Cocca, A.; Greier, K. Association of modifiable factors with the development of physical fitness of Austrian primary school children: A 4-year longitudinal study. J. Sport. Sci. 2022, 40, 920–927. [Google Scholar] [CrossRef] [PubMed]
  10. Rodrigues, L.P.; Stodden, D.F.; Lopes, V.P. Developmental pathways of change in fitness and motor competence are related to overweight and obesity status at the end of primary school. J. Sci. Med. Sport 2016, 19, 87–92. [Google Scholar] [CrossRef]
  11. Haugen, T.; Johansen, B.T. Difference in physical fitness in children with initially high and low gross motor competence: A ten-year follow-up study. Hum. Movment Sci. 2018, 62, 143–149. [Google Scholar] [CrossRef]
  12. Bai, Y.; Saint-Maurice, P.F.; Welk, G.J. Fitness Trends and Disparities among School-Aged Children in Georgia, 2011–2014. Public Health Rep. 2017, 132, 39S–47S. [Google Scholar] [CrossRef]
  13. Reyes, A.C.; Chaves, R.; Baxter-Jones, A.D.G.; Vasconcelos, O.; Tani, G.; Maia, J. A mixed-longitudinal study of children’s growth, motor development and cognition. Design, methods and baseline results on sex-differences. Ann. Hum. Biol. 2018, 45, 376–385. [Google Scholar] [CrossRef]
  14. Ministério da Segurança Social e do Trabalho. Decreto-Lei, n.o 176/2003 de 2 de Agosto; Ministério da Segurança Social e do Trabalho: Lisboa, Portugal, 2003. [Google Scholar]
  15. Kiphard, E.; Schilling, F. Körperkoordinationstest für Kinder; Beltz Test GmbH: Weinheim, Germany, 1974. [Google Scholar]
  16. Barreira, T.V.; Schuna, J.M., Jr.; Mire, E.F.; Katzmarzyk, P.T.; Chaput, J.-P.; Leduc, G.; Tudor-Locke, C. Identifying Children’s Nocturnal Sleep Using 24-h Waist Accelerometry. Med. Sci. Sport. Exerc. 2015, 47, 937–943. [Google Scholar] [CrossRef]
  17. Evenson, K.R.; Catellier, D.J.; Gill, K.; Ondrak, K.S.; McMurray, R.G. Calibration of two objective measures of physical activity for children. J. Sport. Sci. 2008, 26, 1557–1565. [Google Scholar] [CrossRef]
  18. Hedeker, D.; Gibbons, R. Longitudinal Data Analysis; Wiley Interscience: New York, NY, USA, 2006. [Google Scholar]
  19. Hedeker, D.; Gibbons, R.; du Toit, M.; Cheng, Y. Supermix: Mixed Effects Models, Version 1.0.; Scientific Software International: Lincolnwood, IL, USA, 2008. [Google Scholar]
  20. Branta, C.; Haubenstricker, J.; Seefeldt, V. Age changes in motor skills during childhood and adolescence. Exerc. Sport Sci. Rev. 1984, 12, 467–520. [Google Scholar] [CrossRef]
  21. Guthold, R.; Stevens, G.A.; Riley, L.M.; Bull, F.C. Global trends in insufficient physical activity among adolescents: A pooled analysis of 298 population-based surveys with 1·6 million participants. Lancet Child Adolesc. Health 2020, 4, 23–35. [Google Scholar] [CrossRef]
  22. Barnett, L.M.; van Beurden, E.; Morgan, P.J.; Brooks, L.O.; Beard, J.R. Gender Differences in Motor Skill Proficiency from Childhood to Adolescence. Res. Q. Exerc. Sport 2010, 81, 162–170. [Google Scholar] [CrossRef]
  23. Nyquist, A.; Moser, T.; Jahnsen, R. Fitness, Fun and Friends through Participation in Preferred Physical Activities: Achievable for Children with Disabilities? Int. J.Disabil. Dev. Educ. 2016, 63, 334–356. [Google Scholar] [CrossRef]
  24. Stodden, D.F.; Gao, Z.; Goodway, J.D.; Langendorfer, S.J. Dynamic relationships between motor skill competence and health-related fitness in youth. Pediatr. Exerc. Sci. 2014, 26, 231–241. [Google Scholar] [CrossRef]
  25. Asmussen, E. Growth in Muscular Strength and Power, in Physical Activity, Human Growth and Development, 1st ed.; Academic Press: New York, NY, USA, 1973. [Google Scholar]
  26. Corlett, J.T. Power function analysis of physical performance by Tswana children. J. Sport. Sci. 1984, 2, 131–137. [Google Scholar] [CrossRef]
  27. Béghin, L.; Vanhelst, J.; Drumez, E.; Kersting, M.; Molnar, D.; Kafatos, A.; De Henauwn, S.; Wildhalm, K.; Karaglani, E.; Moreno, L.A.; et al. Birth weight and breastfeeding are differentially associated with physical fitness components. Eur. J. Clin. Nutr. 2022, 76, 871–878. [Google Scholar] [CrossRef]
  28. Robič Pikel, T.; Starc, G.; Strel, J.; Kovač, M.; Babnik, J.; Golja, P. Impact of prematurity on exercise capacity and agility of children and youth aged 8 to 18. Early Hum. Dev. 2017, 110, 39–45. [Google Scholar] [CrossRef]
  29. Augste, C.; Lämmle, L.; Künzell, S. Does current behaviour predict the course of children’s physical fitness? Eur. J. Sport Sci. 2015, 15, 429–435. [Google Scholar] [CrossRef]
  30. Breau, B.; Brandes, M.; Veidebaum, T.; Tornaritis, M.; Moreno, L.A.; Molnár, D.; Lissner, L.; Eiben, G.; Lauria, F.; Kaprio, J.; et al. Longitudinal association of childhood physical activity and physical fitness with physical activity in adolescence: Insights from the IDEFICS/I.Family study. Int. J. Behaviral Nutr. Phys. Act. 2022, 9, 147. [Google Scholar] [CrossRef] [PubMed]
  31. Prahl-Andersen, B.; Kowalski, C.J. A mixed longitudinal, interdisciplinary study of the growth and development of Dutch children. Growth 1973, 37, 281–295. [Google Scholar] [PubMed]
Table 1. Number of children per cohort and age group.
Table 1. Number of children per cohort and age group.
CohortsAges of Follow-UpTotal
Cohort 1456 41
Cohort 2 567 65
Cohort 3 678 58
Cohort 4 789 59
Cohort 5 8910 58
Cohort 6 9101167
348
Table 2. Descriptive statistics (mean ± SD or percentage) for girls (n = 177) and boys (n = 171) from 5 to 11 years.
Table 2. Descriptive statistics (mean ± SD or percentage) for girls (n = 177) and boys (n = 171) from 5 to 11 years.
TestsGirls
5 years
(n = 28)
6 years
(n = 58)
7 years
(n = 73)
8 years
(n = 73)
9 years
(n = 86)
10 years
(n = 47)
11 years
(n = 23)
Physical fitness
Handgrip strength (kgf)6.2 ± 1.18.0 ± 2.010.0 ± 2.512.7 ± 3.114.1 ± 2.916.4 ± 3.319.4 ± 3.5
Standing long jump (cm)76.3 ± 16.393.7 ± 17.3102.7 ± 16.5113.8 ± 16.9120.9 ± 18.2131.3 ± 18.3135.1 ± 17.6
Shuttle run (s)16.6 ± 2.014.9 ± 1.514.0 ± 1.213.3 ± 1.112.9 ± 1.112.1 ± 1.011.8 ± 1.0
Gestational Information
Birth weight (kg)3.0 ± 0.63.0 ± 0.53.1 ± 0.53.2 ± 0.53.2 ± 0.53.2 ± 0.53.4 ± 0.4
Anthropometry
BMI (kg·m−2)16.7 ± 1.816.9 ± 2.017.7 ± 2.918.2 ± 3.219.0 ± 3.319.3 ± 3.420.2 ± 3.4
Gross Motor Coordination
GMCTS (points)69.2 ± 22.5103.0 ± 27.9127.9 ± 28.8155.8 ± 31.2173.7 ± 32.4207.1 ± 36.7216.5 ± 42.7
Physical activity
MVPA (min∙day−1)64.6 ± 21.161.6 ± 18.960.3 ± 18.158.9 ± 19.354.0 ± 15.858.2 ± 17.864.2 ± 32.2
Socioeconomic status
A (up to EUR 2.934 per year−1)32.1%15.5%21.9%20.5%22.1%21.3%26.1%
B (EUR 2.934 to EUR 5.896 per year−1)32.1%36.2%37.0%26.0%26.7%21.3%30.4%
C (≥EUR 5.870 per year−1)35.7%48.3%41.1%53.4%51.2%57.4%43.5%
TestsBoys
5 years
(n = 44)
6 years
(n = 62)
7 years
(n = 67)
8 years
(n = 60)
9 years
(n = 60)
10 years
(n = 39)
11 years
(n = 15)
Physical fitness
Handgrip strength (kgf)7.0 ± 2.39.1 ± 2.411.5 ± 2.513.8 ± 3.415.5 ± 3.517.8 ± 4.420.1 ± 4.4
Standing long jump (cm)85.6 ± 21.499.2 ± 19.2114.7 ± 16.8123.6 ± 18.6130.4 ± 16.9139.6 ± 15.3143.5 ± 27.0
Shuttle run (s)15.8 ± 2.114.7 ± 1.713.4 ± 1.113.0 ± 1.112.4 ± 1.111.7 ± 0.711.6 ± 1.1
Gestational Information
Birth weight (kg)3.3 ± 0.63.4 ± 0.53.4 ± 0.53.3 ± 0.53.2 ± 0.53.3 ± 0.63.4 ± 0.5
Anthropometry
BMI (kg·m−2)17.0 ± 1.917.0 ± 2.117.2 ± 2.317.9 ± 3.718.5 ± 3.918.8 ± 4.319.0 ± 3.1
Gross Motor Coordination
GMCTS (points)66.7 ± 27.0103.1 ± 33.4129.9 ± 37.1156.0 ± 33.2174.6 ± 38.2201.5 ± 28.3215.5 ± 44.4
Physical activity
MVPA (min∙day−1)83.7 ± 23.582.6 ± 22.077.2 ± 24.174.8 ± 22.971.7 ± 20.868.5 ± 22.671.4 ± 19.8
Socioeconomic status
A (up to EUR 2.934 per year−1)15.9%14.5%14.9%11.7%15.0%7.7%-
B (EUR 2.934 to EUR 5.896 per year−1)18.2%24.2%23.9%25.0%23.3%33.3%46.7%
C (≥EUR 5.870 per year−1)65.9%61.3%61.2%63.3%61.7%59.0%53.3%
BMI = body mass index; GMCTS = gross motor coordination score; MVPA = moderate-to-vigorous physical activity.
Table 3. Descriptive statistics of the school environment.
Table 3. Descriptive statistics of the school environment.
Schools (n = 19)Mean ± SDMin–Max
School Characterization
School size
 Number of children23 ± 225–87
 Number of teachers4 ± 31–8
Ratio teachers/students0.13 ± 0.110.05–0.55
n (%)
School setting
 Rural12 (63.2)
 Semi-urban7 (36.8)
Policies and Practices for PA
 Policies and practices10 (52.6)
 Policies5 (26.3)
 Practices4 (21.1)
Physical Structure of the School
Playground
 Yes19 (100)
 No0 (0)
Playground area
 With obstacles19 (100)
 Without obstacles0 (0)
Playground dimension
 Smaller (10 m2 to 39 m2)2 (10.5)
 Medium (40 m2 to 69 m2)4 (21.1)
 Large (>70 m2)13 (68.4)
Multi-sports roofed
 Yes5 (26.3)
 No 14 (73.7)
Number of Infrastructures
 One Infrastructure15 (78.9)
 Two Infrastructures4 (21.1)
Equipment for PA
 Yes15 (78.9)
 No4 (21.1)
PE Classes
Duration of PE classes
 45 min6 (31.6)
 60 min13 (68.4)
Active in PE classes
 30 min6 (31.6)
 40 min4 (21.1)
 50 min9 (47.3)
Human Resources
Academic Degree (all graduated)19 (100)
PA = physical activity; PE = physical education.
Table 4. Parameter estimates (standard errors) for fixed and random effects of the three multilevel models.
Table 4. Parameter estimates (standard errors) for fixed and random effects of the three multilevel models.
Handgrip Strength (kgf)Standing Long Jump (cm)Shuttle-Run Test (s)
Regression coefficients (fixed effects—β)
Intercept (5 years)7.74 (0.43) ***88.95 (2.72) ***15.49 (0.19) ***
Age (velocity)1.43 (0.19) ***12.29 (1.20) ***−0.93 (0.09) ***
Age2 (acceleration)0.06 (0.03) *−1.08 (0.17) ***0.07 (0.01) ***
Sex (boys)1.13 (0.25) ***7.88 (1.48) ***−0.33 (0.09) ***
CE_c2_1−0.80 (0.32) *−8.81 (2.38) ***0.18 (0.21) ns
CE_c3_2−0.87 (0.44) ns−12.48 (2.91) ***0.62 (0.24) **
CE_c4_3−1.18 (0.49) *−8.18 (2.90) **0.68 (0.22) **
CE_c5_4−1.25 (0.46) **−5.27 (2.52) *0.62 (0.18) ***
CE_c6_5−1.36 (0.35) ***−4.35 (1.73) * 0.50 (0.11) ***
Birth weight (kg)0.43 (0.23) ns0.25 (1.41) ns−0.18 (0.09) *
SES (level B)−0.35 (0.35) ns−0.55 (2.13) ns−0.10 (0.13) ns
SES (Level C)−0.50 (0.32) ns−0.94 (1.93) ns0.02 (0.12) ns
BMI (kg·m−2)0.35 (0.04) ***−0.93 (0.23) ***0.06 (0.01) ***
GMCTS (points)0.02 (0.003) ***0.23 (0.02) ***−0.02 (0.001) ***
MVPA (min·day−1)0.001 (0.003) ns0.08 (0.02) **−0.003 (0.002) *
Variance components (random effects—σ2)
Child Level
Intercept

1.90 (0.58) **

217.67 (36.43) ***

2.01 (0.27) ***
 Age0.21 (0.10) *9.36 (3.42) **0.06 (0.02) ***
 Covariance (intercept/age)
Residual Level
Intercept
0.09 (0.23) ns

1.95 (0.16) ***
−33.36 (10.79) **

74.79 (6.25) ***
−0.34 (0.06) ***

0.44 (0.03) ***
Model Summary
Number of estimated parameters191919
CE = cohort effects; CE_c2_1 is the overlapping effect of cohort 2 on cohort 1; SES = socioeconomic status; BMI = body mass index; GMCTS = gross motor coordination score; MVPA = moderate-to-vigorous physical activity; ns = non-statistically significant; * p < 0.05; ** p < 0.01; *** p < 0.001.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Reyes, A.; Chaves, R.; Vasconcelos, O.; Pereira, S.; Tani, G.; Stodden, D.; Hedeker, D.; Maia, J.; Baxter-Jones, A. Modeling the Dynamics of Children’s Musculoskeletal Fitness. Int. J. Environ. Res. Public Health 2023, 20, 2938. https://doi.org/10.3390/ijerph20042938

AMA Style

Reyes A, Chaves R, Vasconcelos O, Pereira S, Tani G, Stodden D, Hedeker D, Maia J, Baxter-Jones A. Modeling the Dynamics of Children’s Musculoskeletal Fitness. International Journal of Environmental Research and Public Health. 2023; 20(4):2938. https://doi.org/10.3390/ijerph20042938

Chicago/Turabian Style

Reyes, Ana, Raquel Chaves, Olga Vasconcelos, Sara Pereira, Go Tani, David Stodden, Donald Hedeker, José Maia, and Adam Baxter-Jones. 2023. "Modeling the Dynamics of Children’s Musculoskeletal Fitness" International Journal of Environmental Research and Public Health 20, no. 4: 2938. https://doi.org/10.3390/ijerph20042938

APA Style

Reyes, A., Chaves, R., Vasconcelos, O., Pereira, S., Tani, G., Stodden, D., Hedeker, D., Maia, J., & Baxter-Jones, A. (2023). Modeling the Dynamics of Children’s Musculoskeletal Fitness. International Journal of Environmental Research and Public Health, 20(4), 2938. https://doi.org/10.3390/ijerph20042938

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop