Does maturity status explain training intensity and physical fitness variations of youth football players?

DOI: https://doi.org/10.21203/rs.3.rs-1635086/v1

Abstract

Background: Different predictive maturity status methods are available for coaches to apply during the training process. However, understanding the potential use of such methods to explain both physical fitness and training intensity (TI) measures, can give coaches new insights to better adjust training according to each individual needs.

Objectives: The purpose of the present study was two-fold: (1) To analyse TI and physical fitness variations of youth football players after a full-season; and (2) to examine whether training intensity and physical fitness variations are explainable by estimated maturation status.

Methods: Twenty-seven youth elite Under-15 football players were daily monitored for TI measures during 38 weeks. At the beginning and at the end of the season, all players were assessed for physical fitness. Also, different methods of maturity status estimation were collected at the beginning of the season.

Results: Significant differences were found for all physical fitness and TI measures after the season. The 2-digit (2D) and 4-digit (4D) of left and right hands had negative moderate correlations with change of direction (COD) variations (r=-.39 to -.45 | p=.05 to .02). The right-fingers 2D:4D (RF2D:4D) had positive moderate-to-large correlations with all TI measures (r=-.40 to -.56 | p=.04 to .001). From the reported significant correlations, the RF2D:4D measure explained the greatest proportion of all TI variations (r=.40 to .62 | p=.04 to .001).

Results: The use of the 2D:4D method of maturity status estimation may be a potential predictor of TI variations. Furthermore, the maturity offset and 2D:4D methods may not be useful to predict physical fitness variations after an Under-15 football season. 

Background

The physiological and physical characteristics of academy football players are well described at different age-categories [13]. For instance, youth football players may present VO2max values as high as approximately 67 ml/kg/min and lactate concentrations that can be between 2 to 12 11.9 mmol·L-1 during a football match [4]. During a youth football match, players spent a great proportion of playing time above 80% of their individual maximal heart rate (measured by a cardiorespiratory test until exhaustion), for all age-categories [5, 6]. Also, the anaerobic power values (measured by the Wingate Anaerobic Test (WAnT), can reach to approximately 11 W·kg-1 in youth football players [7].

To cope with football match physiological and physical demands, youth players must gradually develop great physical fitness levels [8]. While a focus on the development of technical skills, agility and running speed seem to be more important in the U13 and U14 age-categories, the development of cardiorespiratory capacity can be crucial in U15 and U16 academy football players [9]. Also, at those age-categories (U15 and U16), the development and demonstration of strength, lower-limbs power and sprint performance assume imperative roles during match performance [10]. In fact, it was previously demonstrated that higher levels of lower-body strength showed strong associations with better sprint and jump performance in youth football players [10]. Furthermore, match running performance and the physical capacities described above seem to be related to each other, although its magnitude tends to be different when categorizing the players into the different on-field playing positions. Moreover, field tests aerobic performances such as during the Yo-Yo Intermittent Recovery (YYIR) test, are related with the capacity to spend more time executing high-intensity activities during a football match [12].

A football season is usually divided by three different moments (e.g., pre-season, in-season, and off-season), where it is expected that both youth and adult players present significant variations in terms of physical fitness during those different periods of the season [13]. However, it seems to persist some incongruences in literature regarding physical fitness changes throughout a football season [14, 15]. Despite that, it is well described that the pre-season is considered an important period where it is supposed to occur significant physical fitness improvements [16]. Those improvements can be even more pronounced if a training plan during the off-season period was not well prescribed for players of the same team [13]. Furthermore, during the in-season period slight improvements can be observed in some players, although there is a tendency for a maintenance of physical status during the in-season, followed by decreases during the final stages of the season [17]. The training-dose imposed to youth football players, and their accumulated perceived intensity of training, explain some magnitude of the observed physical fitness changes [18].

Furthermore, players’ maturational status might also influence youth players physical fitness changes throughout a full-season [19]. The physical development of males during the adolescence, usually starts peaking at approximately 14/15 of their chronological age [20]. However, this can vary significantly if it is considered the biological age of players. That is, players of the same chronological age can be at different levels of maturation [21]. Players that are more advanced in maturity status compared to their same-age colleagues, usually present greater physical performance mainly in strength manifestation as a consequence of greater increases in muscle mass [22]. On the other hand, more mature players in the U14 and U15 age-categories, can express the so-called “motor-awkwardness”, that may limit their performance in other tasks, such as technical skills [23].

In this sense, a recent study conducted on 88 youth football players from different age-categories (U12, 13, 14 and 15), revealed that the accumulated training, maturity and initial physical fitness status explained only small and inconsistent proportions of the observed physical fitness changes after a full-season [24]. Another study conducted on 68 youth football players from different age-categories revealed that player’s maturation status have a moderate effect on match work rate [25]. Interestingly, a recent study revealed that the observed physical fitness changes after a full-season seemed to be influenced by both accumulated training intensity and maturation status [26].

Despite that, few studies focused on the influence that maturity status has on both training intensity physical fitness variations of U15 teams after a full-season [24, 26]. For those reasons, the present study aims to analyse TI and physical fitness variations of youth football players after a full-season and to examine whether training intensity and physical fitness variations are explainable by estimated maturation status.

Materials And Methods

Design and procedures

A prospective cohort study design was used in this study. An under-15 youth elite football team was analyzed throughout 33 weeks. Anthropometric and body composition measures were conducted. Also, the maturity status of each player was assessed using anthropometric data. Physical assessments were carried out in August 2021, and in April 2022. The participants were assessed during three days. Anthropometry, body compositions, flexibility and change of direction (COD) assessments were carried out in the first day. The assessment of anaerobic performance was made during the second day, using the Wingate test. In the third day of assessments, the 30-15 intermittent fitness test (30-15 IFT) was conducted to assess the participant’s aerobic performance. Only the COD and the 30-15 IFT assessments were conducted outdoors, on a synthetic turf soccer field. 

Participants

Twenty-seven youth elite football players (age: 15.0 ± 0.4 years old; height: 175 ± 0.6 cm; body mass: 62.1 ± 7.0 kg) from the same team participated in this study. The inclusion criteria were: (i) all players had to participate in at least 90% of training sessions throughout the season; (ii) for each week, the players had to participate in all training sessions; and (iii) not be injured during the observations and assessments. The goalkeepers were excluded from the sample. Before the beginning of this study, all participants and their parents or their legal representants signed a written informed consent form. Thus, all the advantages and disadvantages of the study procedures were well explained to all the involved. The present study followed the ethical recommendations for the study in humans as suggested by the Declaration of Helsinki (updated version from 2013). 

Internal intensity quantification

Internal intensity was collected using the rate of perceived exertion (RPE) based on the CR-10 Borg scale [27]. Thus, the RPE values were collected approximately 10–30 min after each training session, as recommended in previous research [28]. Based on the CR-10 scale, 1 means “very light activity” and 10 means “maximal exertion”. All players answered to the question “How intense was your session?”. Their responses were given in an individual way, and without the influence of their colleagues. Additionally, the duration of the training sessions, in minutes, was recorded. After obtaining each player RPE value, the session-rate of perceived exertion (s-RPE) was used [29]. To obtain the s-RPE values, the duration of each training session was recorded and multiplied by the RPE value attributed by each player, and was presented as an arbitrary unit (A.U.).

Moreover, from the s-RPE values, the weekly training intensity (wTI) (sum of the intensity of all sessions and match), mean training intensity (mTI) (mean of the intensity of all sessions and match), 5-day average (5d-AVG) (mean of the intensity of five training sessions without match), the training monotony (TM) (mean of training intensity of 7 days divided by the standard deviation), and the training strain (TS) (sum of the intensity of all sessions and match multiplied by training monotony per week) were calculated.  

Physical Fitness Assessments

Anthropometry and body composition

Each participant height was measured using a stadiometer (Seca model 213, Germany) with an accuracy of ±5 mm. While, the participants weight was measured using a balance (Seca model 813, UK) with a precision of 0.1 per kilogram. All players were assessed without shoes and with their lower back as close to the stadiometer as possible. For measuring body composition, three-point skinfolds (chest, abdominal and thigh) were conducted to measure the the body fat percentage (BF%). All skinfold measures were assessed using a Lafayette caliper (Lafayette, IN, USA) with an accuracy of 0.1 mm. The skinfold measurements were applied twice on the right side of the athlete’s body, and the final score recorded were the mean of two measurements. If the measurement error was high (>5%), the measurements had to be performed again and the median of the three repetitions were used for analysis. All measurements were performed by an ISAK credited person. Thus, the calculations of body density and body fat% were calculated based on the Jackson and Pollock formula [30]. 

Sit and Reach test

For conducting the Sit and Reach test, all participants had to sit on the floor with their bare feet against the sit-and-reach equipment and with their middle fingers stacked on top of one another. Participants were informed to stretch as far as possible without bending their knees. The final outcome to be used was the distance between the tip of the middle fingers and the toe line, as previously recommended [31]. 

Modified 5-0-5 COD test

For measuring the participant's ability to change directions, the modified 505 test was used as in elsewhere [32]. Three cones were placed at 5-meters apart from each other, and a pair of photocells with a digital timer connected to it was placed at cone B. The photocell system used was the Newtest Power timer 300-series, that was adjusted to each player’s hip height. Each participant started the test 70 cm before the cone A (starting line). After a beep sound, each participant had to run as quickly as possible until reaching cone C, turn on the cone C line and return as quickly as possible through the photocells (cone B). Test time was measured to the nearest 0.01 s with the fastest value obtained from 2 maximal trials. After each trial, the players had a 3-minute recovery. 

Wingate test

For measuring the anaerobic performance of each participant, the Wingate Anaerobic Test (WAnT) was performed on a cycle ergometer (Monark model 894-E, Vansbro, Sweden). During 5 seconds, the participants had to pedal at maximum speed to determine the repetition per minute (RPM) in the ergometer monitor. After that, a braking force was determined by the product of body mass in kg by 0.075. The participants had to pedal at their maximum effort during 30 seconds with verbal encouragement from the coach and/or colleagues. The peak power (PP) and fatigue index (FI) measures were used for further analysis [33]. 

30-15 Intermittent Fitness test

For measuring the aerobic performance, the 30-15 Intermittent Fitness test (30-15 IFT) was applied. The test initial velocity was set at 8 km.h−1 during the first run, and was increased by 0.5 km/h-1 after each running sequence. All participants had to run back and forth within a 40-meter straight line. Each shuttle consists of 30-second runs interspersed with 15-seconds of walking. 3-meter zones were delineated in both extremities and at the middle of the test setup. Each participant had to complete as many stages as possible, and the test ended when the players could not maintain the running speed demanded, or could not reach the 3-meter zone before the beep during three times. As a final outcome to be analyzed, the velocity of intermittent fitness test (VIFT) score of each participant was recorded. The VIFT consists of the final velocity recorded during the last stage [34, 35]. Also, the Vo2max was estimated by the following equation for each player [36]: Estimated Vo2max = 28.3 −(2.15 × 1) −(0.741 × age) −(0.0357 × mass) + (0.0586 × age × VIFT) + (1.03 × VIFT). 

Maturity Status

Maturity Offset and Age at PHV

To determine the age at peak heigh velocity (PHV) of each player, the maturity offset was calculated using the chronological age, standing height, sitting height, leg length and body weight measures, according to the following equation [37]: Maturity Offset = −9.236 + 0.0002708 (leg length × sitting height) −0.001663 (chronological age × leg length) + 0.007216 (chronological age × sitting height) + 0.02292 (mass by height ratio). For measuring the sitting height, the athletes were asked to sit on the 50 cm height box, facing forwards. Then the height between the highest point of the head and the bottom of the box the player was sitting in, was measured. For measuring the leg length, the standing height minus the sitting height was calculated for each athlete. Finally, to obtain the age at PHV, the chronological age was subtracted by the maturity offset score of each player. 

2D:4D Ratio

The digit two (2D) and digit four (4D) length fingers were measured [38]. Each player placed the right and left-hand palm on a scanner with the fingers kept 2 cm apart from each other. The image of player’s palms in scanner was transferred to a computer and the Kinovea software was used to analyse fingers’ length. The 2D and 4D fingers length were measured from the proximal phalanx to the distal phalanx. The ratio of both fingers was calculated as follows: 2D/4D. The model of the scanner was the Scanjet (5590 HP Scanjet, USA) with an accuracy of 0.01 cm measurement of second and fourth finger to the tip of the finger. The difference of right-fingers 2D:4D ratio (RF2D:4D) minus left-fingers 2D:4D ratio (LF2D:4D) was calculated [39]. The intra-observer reliability was assessed by the same observer two times a week apart. The intra-class correlation (ICC) for 2D:4D ratio was 0.93 and 0.95, respectively. 

Statistical Analysis

Tests of normal distribution and homogeneity (Kolmogorov–Smirnov and Levene’s, respectively) were conducted on all data before analysis. A repeated-measures ANOVA was used to analyse TI measures from W1 to W38. Effect size was indicated with partial eta squared for Fs. A paired sample t-test was used to determine differences as a repeated measures analysis in two conditions (Pre – Post) for physical fitness variables. Cohen d was used as the effect size indicator. To interpret the magnitude of the effect size, we adopted the following criteria [40]: d = 0.20, small; d = 0.50, medium; and d = 0.80, large. Posteriorly, the percentage of change of internal training intensity measures were calculated [100-(Pre*100)/Post] considering Pre as the mean values from W1 to W4 and Post as the mean values from W35 to W38. In addition, percentage of changes were calculated for physical fitness assessments [100-(Pre*100)/Post]. 

Also, the Pearson’s correlation coefficient r was used to examine the relationship between the maturity status [i) maturity-offset and age at PHV; and ii) 2D:4D finger length], and training intensity and physical fitness variations. To interpret the magnitude of these correlations, the following criteria was adopted [41]: r ≤ 0.1, trivial; 0.1 < r ≤ 0.3, small; 0.3 < r ≤ 0.5, moderate; 0.5 < r ≤ 0.7, large; 0.7 < r ≤ 0.9, very large; and r > 0.9, almost perfect. A regression analysis was used to model the prediction of percentage of change of maturity status assessment from the remaining variables with significant correlations. All data were analysed using the software Statistica (version 13.1; Statsoft, Inc., Tulsa, OK, USA) and the significance level was set at p<0.05.

Results

Descriptive statistics were calculated for each variable (Table 1). 

Table 1

Maturity status and physical fitness before (Pre) and after (Post) of seasons (mean ± SD).

 

Pre 

Post

 

Maturity status

 

Maturity-offset 

1.6 ± 0.5

-

 

Age at PHV

13.4 ± 0.3

-

 

2D:4D finger length

Left 2D

7.6 ± 0.5

-

 

Left 4D

7.9 ± 0.5

-

 

LF2D:4D

0.9 ± 0.0

-

 

Right 2D

7.5 ± 0.5

-

 

Right 4D

7.9 ± 0.5

-

 

RF2D:4D

0.9 ± 0.0

-

 

Dif L:R-Ratio

0.0 ± 0.0

-

 

Physical Fitness

t-Test (p)

Sit and -Reach Test (cm)

37.8 ± 8.3

40.7 ± 9.0

p=0.001** | d=-0.32

505 test (sec)

2.1 ± 0.2

1.9 ± 0.1

p=0.001** | d=-0.47

Peak Power (AU)

708.0 ± 128.4

821.0 ± 106.4

p=0.001** | d= 1.40

Fatigue Index (%)

37.1 ± 2.2

42.0 ± 2.7

p=0.001** | d=-0.95

30-15 IFT [V02max (mL·kg−1·min−1)] 

45.0 ± 4.3

47.3 ± 5.1

p=0.001** | d=-2.02

PHV: peak height velocity; D: digit; AU: arbitrary units; RF2D:4D: right finger 2-digit:4-digit; Dif L:R-ratio: difference lef:right-ratio; IFT: intermitent fitness test; * Denotes significance at p<0.05, and ** denotes significance at p<0.01.

 A paired measures t-test with participants’ physical fitness assessment revealed significant differences in Sit and Reach, 505 test, PP, FI and 30-15 IFT (= 0.001, = -0.32; = 0.001, = -0.47; = 0.001, =1.40; = 0.001, = -0.95 and = 0.001, = -2.02, respectively). For more information, see Figure 1. 

A repeated-measures ANOVA with mean data of TI measures revealed a significant main effect of TM: (2.89) = 19.22, = 0.001, η2partial= 0.44; TS: (6.99) = 9.67, = 0.001, η2partial = 0.29; mTI: F (6.43) = 15.24, = 0.001, η2partial = 0.39; wTI: (7.24) = 9.39, = 0.001, η2partial = 0.28, and 5d-AVG, (7.08) = 9.78, = 0.001, η2 = 0.29. For more information, see Figure 2.

A correlation analysis was performed between the maturity status [i) maturity-offset and age at PHV, and ii) Left 2, Left 4, L F Ratio, Right 2, Right 4, R:F Ratio and Dif LR Ratio] and percentage of change of physical fitness assessment [Pre – Post (Sit and Reach test, 5-0-5, Peak Power, Fatigue Index and 30-15 IFT)] and the TI measures (TM, TS, mTI, wTI and 5d-AVG). Only a negative moderate correlation was found between percentage of change of 505 test and maturity offset (r =-.40, p= .04). See table 2. 

Table 2

 Correlations between maturity status (maturity-offset and age at PHV), and the variations of TI and physical fitness measures.  

 

Internal training intensity variations


TM 

TS 

mTI 

wTI 

5d-AVG 

Maturity offset

r=-.27 | p=.18

r=-.12 | p=.56

r=-.18 | p=.38

r=-.17 | p=.41

r=-.18 | p=.38

Age at PHV

r=.21 | p=.30

r=.05 | p=.78

r=.20 | p=.33

r=.17 | p=.41

r=.22 | p=.28

 

Physical fitness changes


S&R

30-15 IFT

505

PP

FI

Maturity offset

r=-.07 | p=.71

r=-.09 | p=.64

r=-.40 | p=.04*

r=.22 | p=.28

r=.18 | p=.37

Age at PHV

r=-.02 | p=.89

r=.10 | p=.61

r=.07 | p=.73

r=-.01 | p=.93

r=-.19 | p=.35

PHV: peak height velocity; TM: training monotony; TS: training strain; mTI: mean training intensity; wTI: weekly training intensity; 5d-AVG: 5 day average; S&R: sit and reach test; 30-15 IFT: 30-15 intermittent fitness test; PP: peak power; FI: fatigue index; * Denotes significance at p<0.05, and ** denotes significance at p<0.01.

In reference to correlation analysis performed between 2D:4D fingers. A negative moderate correlation was found between Right 4 and TM (r =-.40, p= .04) and 505 test (r =-.41, p= .04). In the same line, a negative moderate correlation between Right 2 and 505 test (r =-.45, p=.02)also between Right 4 and TM (r =-.40, p=.04).

On the other hand, positive moderate correlations between LF Ratio and TM, mTI, and wTI (r =.53, p=.01; r =.43, p=03. and r =.41, p= .04, respectively) were found. In the same direction, a positive moderate correlation between RF Ratio and TM and wTI (r =.40, p=.04, and r =.50, p=.01, respectively). Positive large correlations between RF ratio and TS, mTI, and 5d-AVG (r =.56, p=.001; r =.54, p=.01, and r =.51, p=01) were found. Lastly, a positive moderate correlation was found between Dif LR Ratio and TS and FI (r =.43, p=.03, and r =.40, p=05, respectively). See Table 3, for more information. 

Table 3

 Correlations between maturity status (2D:4D finger length) and the variations of TI and physical fitness measures.  

 

Internal Training Intensity


TM 

TS 

mTI 

wTI 

5d-AVG 

Left 2D

r=-.15 | p=.46

r=-.01 | p=.96

r=-.03 | p=.87

r=-.05 | p=.80

r=.04 | p=.83

Left 4D

r=-.38 | p=.06

r=-.15 | p=.46

r=-.14 | p=.48

r=-.23 | p=.26

r=-.10 | p=.60

LF2D:4D

r=.53 | p=.01*

r=.34 | p=.08

r=.43 | p=.03*

r=.41 | p=.04*

r=.36 | p=.07

Right 2D

r=-.18 | p=.36

r=-.01 | p=.97

r=.02 | p=.92

r=-.06 | p=.74

r=.02 | p=.89

Right 4D

r=-.40 | p=.04*

r=-.29 | p=.15

r=-.25 | p=.22

r=-.33 | p=.10

r=-.23 | p=.26

RF2D:4D

r=.40| p=.04*

r=.56 | p=.001**

r=.54 | p=.01*

r=.50 | p=.01*

r=.51 | p=.01*

Dif L:R

r=-.02| p=.91

r=.43 | p=.03*

r=.29 | p=.15

r=.26 | p=.20

r=.34 | p=.09

 

Physical Fitness


S&R

IFT

505

PP

FI

Left 2D

r=-.07 | p=.72

r=.01 | p=.96

r=-.39 | p=.05*

r=.04 | p=.82

r=.01 | p=.92

Left 4D

r=-.11 | p=.57

r=.11 | p=.58

r=-.41 | p=.04*

r=.01 | p=.99

r=.12 | p=.54

LF2D:4D

r=.10 | p=.60

r=-.24 | p=.23

r=.03 | p=.88

r=.10 | p=.61

r=-.28 | p=.17

Right 2D

r=-.02 | p=.92

r=.01 | p=.98

r=-.45 | p=.02*

r=.05 | p=.79

r=.02 | p=.89

Right 4D

r=-.04 | p=.83

r=.05 | p=.79

r=-.44 | p=.03*

r=.01 | p=.94

r=.01 | p=.99

RF2D:4D

r=.05| p=.78

r=-.09 | p=.66

r=-.06 | p=.74

r=.07 | p=.72

r=.05 | p=.80

Dif L:R

r=-.04| p=.84

r=.15 | p=.46

r=-.13 | p=.52

r=-.01 | p=.95

r=.40 | p=.05*

D: digit; RF2D:4D: right finger 2-digit:4-digit; LF2D:4D: left finger 2-digit:4-digit Dif L:R-ratio: difference lef:right-ratio; TM: training monotony; TS: training strain; mTI: mean training intensity; wTI: weekly training intensity; 5d-AVG: 5 day average; S&R: sit and reach test; 30-15 IFT: 30-15 intermittent fitness test; PP: peak power; FI: fatigue index; * Denotes significance at p<0.05, and ** denotes significance at p<0.01.

A multilinear regression analysis was performed to verify which variable of maturity status [i) maturity-offset and age at PHV, and ii) Left 2, Left 4, L F Ratio, Right 2, Right 4, R F Ratio and Dif LR Ratio] could be used to better explain the percentage of change of physical fitness [Pre – Post (Sit and Reach test, 5-0-5, Peak Power, Fatigue Index and 30-15 IFT)] and TI measures (TM, TS, mTI, wTI and 5d-AVG). For more information, see the table 4. 

Table 4

Values of regression analysis explaining percentage of change of the maturity status on the remaining variables.

 

 

b*

SE of b*

R

R2

Adjusted R2

F

p

Maturity offset

505

.40

.19

.40

.16

.13

4.62

.04*

Left 4D

505

-.41

.18

.41

.17

.14

4.87

.04*

LF2D:4D

TM

.53

.17

.57

.28

.25

9.33

.01*

 

mTI

.43

.18

.43

.18

.15

5.29

.03*

 

wTI

.41

.18

.41

.17

.13

4.85

.04*

Right 2D

505

-.45

.18

.45

.20

.17

6.07

.02*

Right 4D

TM

-.40

.19

.40

.16

.13

4.59

.04*

 

505

-.44

.18

.44

.19

.16

5.66

.03*

RF2D:4D

TM

.40

.19

.40

.16

.12

4.57

.04*

 

TS

.56

.17

.56

.32

.29

11.01

.001**

 

mTI

.54

.17

.54

.29

.26

9.50

.01*

 

wTI

.50

.18

.50

.25

.22

8.07

.01*

 

5d-AVG

.62

.16

.62

.38

.36

14.51

.01*

Dif L:R

TS

.43

.18

.43

.18

.15

5.35

.03*

 

FI

.53

.18

.53

.28

.26

8.85

.01*

D: digit; RF2D:4D: right finger 2-digit:4-digit; LF2D:4D: left finger 2-digit:4-digit Dif L:R-ratio: difference lef:right-ratio; TM: training monotony; TS: training strain; mTI: mean training intensity; wTI: weekly training intensity; 5d-AVG: 5 day average; FI: fatigue index; * Denotes significance at p<0.05, and ** denotes significance at p<0.01.

Discussion

The aims of the present study were to analyse the training intensity and physical fitness variations of youth football players after a full-season, and to examine whether training intensity and physical fitness variations are explainable by different estimated maturation status methods. The main findings were that from the analysed maturity status estimations (maturity offset, age at PHV and 2D:4D), only the 2D:4D method revealed moderate to large associations with different TI measures and with COD variations. On the other hand, neither maturity offset nor age at PHV revealed significant relationships with both TI and physical fitness variations, except a moderate association between maturity offset and COD variations. The RF2D:4D measure explained the variations of all TI measures.

There are several studies reporting physical fitness seasonal variations of both juvenile and adult football players [15, 42]. The present study revealed significant improvements of all the analysed physical fitness measures, from the beginning to the end of the season. Indeed, Meckel et al. [15] also reported significant improvements in different physical fitness measures from the beginning to the end of a professional adult football season. The above-mentioned study [15] was conducted on adult elite population, and the physical fitness measures were not similar to the measures used in the present study. Given that, a direct comparison with our findings is not straightforward. However, a recent study [26] conducted on 23 under-16 football players, revealed that VO2max, PP and FI measures had significant variations between pre- and post-assessments, which is in concordance with the findings of the present study. Although some studies reported significant improvements in physical fitness after a football season, others revealed that these changes are not so straightforward as they may be dependent on baseline values [43, 44].

Regarding the different TI measures variations, our findings are congruent with some studies conducted on adult football players that showed significant between-week variations in different perceived TI measures throughout the season [45]. However, the overall literature has focused mainly on adult players and that there is a lack of studies conducted on youth [45]. From the available literature, a recent study conducted on under-16 football players analysed the between-week variations of perceived TI measures (wTI, TM and TS), revealed reduced between-week variability (CV < 6%) for all the analysed TI measures [46]. This is in contrast with our findings, which revealed significant main effects for all the analysed TI measures. These discrepancies in between-week TI variations in youth football may be related to different training methodologies and coach's ideologies, as this can influence team’s training and matches demands [47].

The available research regarding the biological maturity status of young athletes and its associations with physical fitness and TI measures lacks consistency. To the best of the authors knowledge, there are only few studies focusing on the relationships between estimated maturity status, physical performance and TI [2426]. While, only one study focused on the influence of maturity on physical fitness and TI variations [42]. Furthermore, while some studies used the Mirwald’s biological maturity predictive equation [24, 25], other studies used the Khamis-Roche method [48] and the 2D:4D finger length method [42]. In the present study, the 2D:4D finger length method and the Mirwald’s maturity offset and age at PHV were used as predictor variables of both physical fitness and TI measures variations throughout a football season. However, the methodological discrepancies between the above-mentioned studies and the present study, make comparisons more difficult to generalize findings.

The Mirwald’s maturity offset and predicted adult height methods are the most non-invasive methods to predict the maturity status of young athletes [49]. However, the 2D:4D method was previously reported to be a potential predictor of physical fitness variations [42]. Indeed, the 2D:4D ratio of the left and right hands of 24 under-17 football players, had negative large correlations with VO2max variations (r = -0.55, p = 0.005; r = -0.50, p = 0.013) [42], which contrast with our findings. In fact, in the present study, the 2D and 4D of the left and right hands had negative moderate correlations with COD variations, while the maturity offset method presented a negative moderate correlation only with COD variations. Another study determined that the predictors used (initial fitness, maturity offset and training time) to explain the observed physical fitness variations were impacted by other factors that were not measured in their model [24]. However, the study of King et al. [24] used only the maturity offset method as a predictor variable. Although it seems evident that exists inconsistencies regarding the use of maturity offset method as a predictor of physical fitness changes, the use of 2D:4D method seems to be more promising, as showed by our model.

The present study had some limitations. The main limitation refers to the small sample size used. The fact that only one male team was included in the sample is another main limitation. However, in professional youth football competitions, the use of more than one team is a major concern for both coaches and practitioners. Future studies should use larger sample sizes and examine the associations between different predictive maturity status and objective internal TI variations, such as heart rate-based measures.

Conclusions

The overall physical fitness and TI measures revealed significant changes from the beginning to the end of the season. The 2D:4D method showed to be a promising predictor of subjective internal TI variations throughout the season. Specifically, the use of RF2D:4D measure explained the variations of all the analysed TI measures. Despite that, maturity offset, age at PHV and 2D:4D methods explain only a small proportion physical fitness measures, namely, COD performance. Coaches and practitioners can benefit from the use of the 2D:4D method in young football players, to better monitor and adjust weekly TI and ensure greater biological individualization and talent selection.

Abbreviations

COD: change-of-direction; 30-15 IFT: 30-15 intermittent fitness test; s-RPE: session-rate of perceived exertion; TM: training monotony; TS: training strain; wTI: weekly training intensity; mTI: mean training intensity; 5d-AVG: 5-day average; PP: peak power; FI: fatigue index; PHV: peak height velocity; D: digit; RF2D:4D: right finger 2-digit:4-digit; LF2D:4D: left finger 2-digit:4-digit Dif L:R-ratio: difference lef:right-ratio; 

Declarations

Ethics approval and consent to participate

Parents/guardians informed consents were obtained and signed regarding all players involved in this study. The study was conducted according to the Declaration of Helsinki (2013), and was approved by the institute’s research ethics committee (the Scientific Council of Polytechnic Institute of Viana do Castelo), with the code CTC-ESDL-CE005-2021.

 Consent for publication

No individual or indemnifiable data is being published as part of this manuscript. 

Competing interests

The authors declare no competing interests. 

Authors’ contributions

Conceptualization, R.S., F.M.C., and J.M.C.C.; methodology, R.S. and F.M.C.; formal analysis, F.T.G.F.; writing—original draft preparation, R.S., F.M.C., and J.M.C.C.; writing—review and editing, R.S., F.M.C., F.T.G.F., H.N., H.H., and J.M.C.C.; supervision, F.M.C. and J.M.C.C. All authors have read and agreed to the published version of the manuscript. All authors contributed equally to the manuscript and read and approved the final version of the manuscript. 

Funding

This research received no external funding. 

Availability of data and materials

The datasets generated and analysed during the current study are not publicly available due to ethical restrictions, however are available from the corresponding author on reasonable request.

References

  1. Deprez D, Fransen J, Boone J, Lenoir M, Philippaerts R, Vaeyens R. Characteristics of high-level youth soccer players: variation by playing position. J Sports Sci. 2015;33:243–54.
  2. Murr D, Raabe J, Höner O. The prognostic value of physiological and physical characteristics in youth soccer: A systematic review. Eur J Sport Sci. 2018;18:62–74.
  3. Beboucha W, Belkadi A, Benchehida A, Bengoua A. The Anthropometric and Physiological Characteristics of Young Algerian Soccer Players. Acta Fac Educ Phys Univ Comenianae. 2021;61:35–51.
  4. Aslan A, Açikada C, Güvenç A, Gören H, Hazir T, Özkara A. Metabolic demands of match performance in young soccer players. J Sport Sci Med. 2012;11:170–9.
  5. Mendez-Villanueva A, Buchheit M, Simpson B, Bourdon PC. Match play intensity distribution in youth soccer. Int J Sports Med. 2013;34:101–10.
  6. Fitzpatrick JF, Hicks KM, Hayes PR. Dose–Response Relationship Between Training Load and Changes in Aerobic Fitness in Professional Youth Soccer Players. Int J Sports Physiol Perform. 2018;13:1365–70.
  7. Slimani M, Nikolaidis PT. Anthropometric and physiological characteristics of male soccer players according to their competitive level, playing position and age group: A systematic review. J Sports Med Phys Fitness. 2019;:141–63.
  8. Bujnovsky D, Maly T, Ford K, Sugimoto D, Kunzmann E, Hank M, et al. Physical Fitness Characteristics of High-level Youth Football Players: Influence of Playing Position. Sports. 2019;7:46.
  9. Vaeyens R, Malina RM, Janssens M, Van Renterghem B, Bourgois J, Vrijens J, et al. A multidisciplinary selection model for youth soccer: The Ghent Youth Soccer Project. Br J Sports Med. 2006;40:928–34.
  10. Comfort P, Stewart A, Bloom L, Clarkson B. Relationships between strength, sprint, and jump performance in well-trained youth soccer players. J Strength Cond Res. 2014;28:173–7.
  11. Buchheit M, Mendez-Villanueva A, Simpson BM, Bourdon PC. Match running performance and fitness in youth soccer. Int J Sports Med. 2010;31:818–25.
  12. Rebelo A, Brito J, Seabra A, Oliveira J, Krustrup P. Physical match performance of youth football players in relation to physical capacity. Eur J Sport Sci. 2014;14:148–56.
  13. Silva JR, Nassis GP, Rebelo A. Strength training in soccer with a specific focus on highly trained players. Sport Med - Open. 2015;1.
  14. Magal M, Smith R, Dyer JJ, Hoffman JR. Seasonal Variation in Physical Performance-Related Variables in Male NCAA Division III Soccer Players. J Strength Cond Res. 2006;20:874–81.
  15. Meckel Y, Doron O, Eliakim E, Eliakim A. Seasonal Variations in Physical Fitness and Performance Indices of Elite Soccer Players. Sports. 2018;6:14.
  16. Clemente F, Ramirez-Campillo R, Sarmento H. Detrimental Effects of the Off-Season in Soccer Players: A Systematic Review and Meta-analysis. Sport Med. 2021;51:795–814.
  17. Mohr M, Krustrup P, Bangsbo J. Match performance of high-standard soccer players with special reference to development of fatigue. J Sports Sci. 2003;21:519–28.
  18. Gil-Rey E, Lezaun A, Los Arcos A. Quantification of the perceived training load and its relationship with changes in physical fitness performance in junior soccer players. J Sports Sci. 2015;33:2125–32.
  19. Nobari H, Silva AF, Clemente FM, Siahkouhian M, García-Gordillo MÁ, Adsuar JC, et al. Analysis of Fitness Status Variations of Under-16 Soccer Players Over a Season and Their Relationships With Maturational Status and Training Load. Front Physiol. 2021;11:1–11.
  20. Söğüt M, Luz LGO, Kaya ÖB, Altunsoy K, Doğan AA, Kirazci S, et al. Age-and maturity-related variations in morphology, body composition, and motor fitness among young female tennis players. Int J Environ Res Public Health. 2019;16.
  21. Malina RM, Eisenmann JC, Cumming SP, Ribeiro B, Aroso J. Maturity-associated variation in the growth and functional capacities of youth football (soccer) players 13-15 years. Eur J Appl Physiol. 2004;91:555–62.
  22. Gastin PB, Bennett G, Cook J. Biological maturity influences running performance in junior Australian football. J Sci Med Sport. 2013;16:140–5.
  23. John C, Rahlf AL, Hamacher D, Zech A. Influence of biological maturity on static and dynamic postural control among male youth soccer players. Gait Posture. 2019;68:18–22.
  24. King M, Ball D, Weston M, McCunn R, Gibson N. Initial fitness, maturity status, and total training explain small and inconsistent proportions of the variance in physical development of adolescent footballers across one season. Res Sport Med. 2021;00:1–12.
  25. Francini L, Rampinini E, Bosio A, Connolly D, Carlomagno D, Castagna C. Association between Match Activity, Endurance Levels and Maturity in Youth Football Players. Int J Sports Med. 2019;40:576–84.
  26. Nobari H, Silva AF, Clemente FM, Siahkouhian M. Analysis of Fitness Status Variations of Under-16 Soccer Players Over a Season and their Relationships with Maturational Status and Training Load. Front Physiol. 2020;:1–21.
  27. Borg G. Borg’s perceived exertion and pain scales. 1998.
  28. Foster C, Florhaug JA, Franklin J, Gottschall L, Hrovatin LA, Parker S, et al. A new approach to monitoring exercise training. J Strength Cond Res. 2001;15:109–15.
  29. Foster C. Monitoring training in athletes with reference to overtraining syndrome. Med Sci Sports Exerc. 1998;30:1164–8.
  30. Jackson AS, Pollock ML. Generalized equations for predicting body density of men. Br J Nutr. 1978;40:497–504.
  31. Tabacchi G, Sanchez GFL, Sahin FN, Kizilyalli M, Genchi R, Basile M, et al. Field-based tests for the assessment of physical fitness in children and adolescents practicing sport: A systematic review within the ESA program. Sustainability. 2019;11.
  32. Gabbett T, Kelly JN, Sheppard JM. Speed, change of direction speed, and reactive agility of rugby league players. J Strength Cond Res. 2008;22:174–81.
  33. Nikolaïdis P. Anaerobic power across adolescence in soccer players. Hum Mov. 2011;12:342–7.
  34. Buchheit M, Haddad H Al, Millet GP, Lepretre PM, Newton M, Ahmaidi S. Cardiorespiratory and cardiac autonomic response to 30-15 IFT in team sport players. J strength Cond Res. 2009;23:93–100.
  35. Grgic J, Lazinica B, Pedisic Z. Test-retest reliability of the 30–15 Intermittent Fitness Test (IFT): A systematic review. J Sport Heal Sci. 2020;:0–20.
  36. Buchheit M. The 30-15 Intermittent Fitness Test : 10 year review. 2010;1:1–9.
  37. Mirwald RL, Baxter-Jones ADG, Bailey DA, Beunen GP. An assessment of maturity from anthropometric measurements. Med Sci Sports Exerc. 2002;34:689–94.
  38. Manning JT, Scutt D, Wilson J, Lewis-Jones DI. The ratio of 2nd to 4th digit length: A predictor of sperm numbers and concentrations of testosterone, luteinizing hormone and oestrogen. Hum Reprod. 1998;13:3000–4.
  39. Bennett M, Manning JT, Cook CJ, Kilduff LP. Digit ratio (2D:4D) and performance in elite rugby players. J Sports Sci. 2010;28:1415–21.
  40. Cohen J. Statistical power analysis for the behavioral sciences. 2nd edition. New York: Routledge; 1988.
  41. Hopkins W, Marshall S, Batterham A, Hanin J. Progressive statistics for studies in sports medicine and exercise science. Med Sci Sports Exerc. 2009;41:3–12.
  42. Nobari H, Alves AR, Clemente FM, Pérez-Gómez J. Influence of 2D:4D ratio on fitness parameters and accumulated training load in elite youth soccer players. BMC Sports Sci Med Rehabil. 2021;13:1–12.
  43. Los Arcos A, Castillo D, Martínez-Santos R. Influence of initial performance level and tactical position on the aerobic fitness in soccer players after preseason period. Sci Med Footb. 2018;2:294–8.
  44. Silva R, Lima R, Camões M, Leão C, Matos S, Pereira J, et al. Physical fitness changes among amateur soccer players: effects of the pre-season period. Biomed Hum Kinet. 2021;13:63–72.
  45. Teixeira JE, Forte P, Ferraz R, Leal M, Ribeiro J, Silva AJ, et al. Monitoring accumulated training and match load in football: A systematic review. Int J Environ Res Public Health. 2021;18:1–41.
  46. Nobari H, Aquino R, Clemente FM, Khalafi M, Adsuar JC, Pérez-Gómez J. Description of acute and chronic load, training monotony and strain over a season and its relationships with well-being status: A study in elite under-16 soccer players. Physiol Behav. 2020;225.
  47. Machado G, Bach Padilha M, González Víllora S, Clemente FM, Teoldo I. The effects of positional role on tactical behaviour in a four-a-side small-sided and conditioned soccer game. Kinesiology. 2019;51:261–70.
  48. Radnor JM, Staines J, Bevan J, Cumming SP, Kelly AL, Lloyd RS, et al. Maturity has a greater association than relative age with physical performance in english male academy soccer players. Sports. 2021;9:1–13.
  49. Salter J, Cumming S, Hughes JD, De Ste Croix M. Estimating somatic maturity in adolescent soccer players: Methodological comparisons. Int J Sport Sci Coach. 2022;17:11–7.