Next Article in Journal
Genetic Signature of Pinctada fucata Inferred from Population Genomics: Source Tracking of the Invasion in Mischief Reef of Nansha Islands
Previous Article in Journal
The Strigolactone Pathway Is a Target for Modifying Crop Shoot Architecture and Yield
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Do Repeated Sprints Affect the Biceps Femoris Long Head Architecture in Football Players with and without an Injury History?—A Retrospective Study

1
CIPER, Faculdade de Motricidade Humana, Universidade de Lisboa, 1649-004 Cruz Quebrada-Dafundo, Portugal
2
Department of Biochemistry, King’s College London, London WC2R 2LS, UK
*
Author to whom correspondence should be addressed.
Biology 2023, 12(1), 96; https://doi.org/10.3390/biology12010096
Submission received: 19 October 2022 / Revised: 5 January 2023 / Accepted: 6 January 2023 / Published: 8 January 2023
(This article belongs to the Special Issue Biomechanics in Rehabilitation Engineering)

Abstract

:

Simple Summary

Football is a sport which imposes an intermittent type of physical exertion, and one of the main physical demands of the game is sprinting, covering mostly short distances on repeated occasions. Sprinting imposes maximum demands, especially on the hamstrings muscle group; this activity frequently contributes to the occurrence of injury to the biceps femoris long head muscle, particularly in the later stages of the game, suggesting that fatigue may play an important role. After the injury, sprint performance and strength might be reduced, and muscle morphology can be altered even after the player returns to competition. Therefore, the aim of the study was to verify the effect of a fatiguing task on sprint performance, the biceps femoris long head muscle morphology and the hamstring muscle force capacity. Moreover, a comparison was performed between injured players and healthy controls and injured versus contralateral limbs. Muscle morphology was assessed through ultrasound techniques, while muscle force was estimated using a dynamometry assessment. It was concluded that injured players tend to be faster than the healthy control group and display shorter biceps femoris long head fascicle length when compared with the contralateral limb and healthy controls.

Abstract

The aim of this study was to compare the biceps femoris long head (BFlh) architecture between football players with (twelve) and without (twenty) history of BFlh injury before and after a repeated sprint task. Fascicle length (FL), pennation angle (PA) and muscle thickness (MT) were assessed at rest and in the active condition before and after the repeated sprint protocol. Athletes with previous BFlh injury showed shorter FL at rest (p = 0.014; η2p = 0.196) and active state (p < 0.001; η2p = 0.413), and greater PA at rest (p = 0.002; η2p = 0.307) and active state (p < 0.001; η2p = 0.368) before and after the task. Intra-individual comparisons showed that injured limbs have shorter FL at rest (p = 0.012; η2p = 0.519) and in the active state (p = 0.039; η2p = 0.332), and greater PA in passive (p < 0.001; η2p = 0.732) and active conditions (p = 0.018; η2p = 0.412), when compared with contralateral limbs. Injured players, at rest and in the active condition, display shorter BFlh FL and greater PA than contralateral and healthy controls after repeated sprints. Moreover, the BFlh of injured players presented a different architectural response to the protocol compared with the healthy controls.

1. Introduction

Hamstring strain injury (HSI) is the most common non-contact muscle injury in high-speed running sports [1,2,3,4,5,6,7], with a particularly high injury rate in football. It accounts for 37% of all muscle injuries [8], with a recurrence of 12–33% [7,9,10], which implies a significant financial cost for clubs. Indeed, the average cost of having a professional first team player injured for 1 month is approximately EUR 500.000 in European football leagues [11]. Due to these implications, HSI has drawn the attention of sports scientists and researchers. However, until today, the most consistent risk factor is recent (same season) HSI and previous HSI, although other factors such as old age, anterior cruciate ligament injury history and calf injury history have been suggested more recently [12]. It should also be noted that injuries occurred in the later stages of each half [8], suggesting that fatigue could also be a potential risk factor.
One way to induce fatigue is through repeated sprints [13], with the sprint having been identified as the most common mechanism for sustaining an HSI [4], particularly in the biceps femoris long head (BFlh) [14]. Despite limited evidence, the FL of the BFlh has been associated with a risk of HSI [12], since previous studies have reported differences in BFlh architecture between previously injured and uninjured players [15,16,17]. The potential physiological hypothesis is that the BFlh and the semimembranosus (SM) are two muscles with shorter fibers and greater PA, indicating a lower excursion capacity [18]. In terms of injury, muscles with greater fiber lengths have a greater lengthening capacity than muscles consisting of shorter fibers [19]. Therefore, when all hamstrings contract as they lengthen (eccentrically), muscles with shorter fibers will undergo higher strain, potentially predisposing the BFlh and SM to greater injury risk [20]. Curiously, these two muscles are the most affected in HSI, with most injuries in the BFlh occurring during sprint in the swing late phase [14,21] and in the SM during stretching-type movements (i.e., during movements with large joint excursions) [22,23]. Therefore, the rationale exists to support the hypothesis that shorter fascicles are a risk factor for HSI. Indeed, it has been reported that having shorter BFlh fascicles increased the risk of a future HSI, and low levels of eccentric knee flexor strength increased the risk of an HSI occurring in the subsequent season [15]. Moreover, Timmins et al. (2016) reported that the increased risk associated with increasing age and a history of HSI can be mitigated by higher levels of eccentric knee flexor strength and longer BFlh fascicles. It should be noted that muscle architecture is altered significantly when the muscle is active [19,24], with the serial sarcomeres number (which determines the FL) playing a role in contraction velocity and shortening capacity, while sarcomeres in parallel (linked to the muscle cross-sectional area and positively correlated with PA) dictate the force output potential [25]. Therefore, players with previous muscle damage, specifically HSI, would be expected to display a shorter FL and greater PA than healthy players with no previous HSI and compared with the contralateral (uninjured) limb due to a reduced sarcomere number in series. However, it should be interesting to analyze the muscle architecture after an ecological task with the most common injury mechanism, the sprint, as task specificity could induce different outcomes [26].
The purpose of this study was to compare the acute effects of a repeated sprint protocol on resting and active state BFlh architecture between professional football players with and without BFlh injury. Architecture was determined using ultrasound B-mode. We hypothesized that injured players would have shorter fascicles in comparison with their contralateral limbs and the uninjured players. Secondly, we hypothesized that the sprint task will induce changes on BFlh architecture with shorter FL after the task in the injured players compared with the healthy controls and contralateral limbs.

2. Materials and Methods

Thirty-two male professional football players (age = 25.03 ± 3.72 years, height = 1.79 ± 0.07 m, body mass = 74.69 ± 8.67 kg) participated in the present study divided in two groups. The healthy control group was composed of twenty professional football players without BFlh strain injury history. The injury group was composed of twelve professional football players with a unilateral BFlh strain injury history within the last 11.7 ± 7.4 months (ranging between 2 and 23 months). Lower limb injury history was registered by a physiotherapist, obtained through interview and by using musculoskeletal injury information recorded in the club’s medical department. Each injury diagnosis was based on MRI, ultrasound or palpation examination. All the players returned to pre-injury levels of training and competition and were recruited to assess the differences in architecture between their previously injured and uninjured BFlh. All participants read and signed an informed consent form before participation in the study. The Ethical Committee at the Faculty of Human Kinetics at the University of Lisbon approved the study (#5/2021).

2.1. Dynamometry

The knee flexor linear torque was measured at a sampling rate of 1000 Hz using a custom-made equipment. Participants were placed in a prone position, with the hips in neutral anatomical position, knees flexed at 30° (0° = full extension) and the ankle in 15° of plantar flexion, as hamstrings are in minimal passive tension [27]. The feet were fixed in a foot holder which contained a force transducer (Model STC, Vishay Precision, Malvern, PA, USA) at the heel level to collect the linear force perpendicular to the leg orientation and with the ankle at 90°. Force data was amplified (Model UA73.202, Sensor Techniques, Cowbridge, UK), digitally converted (USB-230 Series, Measurement Computing Corporation Norton, MA, USA), recorded using the DAQami software (v4.1, Measurement Computing Corporation, Norton, MA, USA) and multiplied by the perpendicular distance between the force transducer center and the femoral lateral condyle in order to estimate the knee torque. Visual feedback of force production was provided to individuals during the assessments.

2.2. BFlh Architecture Assessment

Two identical ultrasound scanners (Aixplorer, v11; Supersonic Imagine, Aix-en-Provence, France) were used to assess the FL, muscle thickness (MT), and PA of the BFlh which were determined from ultrasound images taken along the longitudinal axis of the muscle belly using an ultrasound transducer (SL10-2, 38 mm, 2–10 MHz. Super Linear, Vermon, Tours, France).
The scanning site was determined with the considerations previously reported [28], at approximately 50% of the femur length (distance between the lateral femoral condyle and greater trochanter). To ensure a stable measure of the muscle architecture during the contractions and the measurement in the same region, a plastic cast was fixed to the skin using bi-adhesive tape due to the importance of keeping the same region of interest (ROI) over time for ultrasound imaging, as errors in probe placement could strongly influence the muscle architecture outcomes.
To minimize the time between the sprint task and the measurements of both limbs, the images were acquired by two examiners to enable the assessment of both limbs simultaneously. The most experienced examiner, having demonstrated their reliability in a previous study analyzing BFlh architecture [28], determined the region of interest by placing the casts, minimising the influence of the other examiner.

2.3. Protocol

Participants visited the track and field high performance center of the Jamor indoor facility. Both limbs were tested for neuromuscular and muscle architecture parameters with the participants in the previously described position. Testing began with participants resting for ~30 min to allow for body fluid shift stabilization [29,30]. Then, two videos were recorded for each BFlh with the individuals at rest. A standardized warm-up for force production was performed (i.e., 20 submaximal and 3 maximal knee flexions), followed by two 3-s MVIC trials with 30 s recovery between trials. Based on the MVIC highest peak torque (PT), individuals familiarized themselves with the 20% of MVIC trials using visual feedback. Each trial lasted ~40 s. The videos during active state were taken twice at 20% of MVIC for each BFlh. After active measures, a standardized warm-up protocol was conducted. Immediately after the warm-up, a 10 × of 30 m repeated sprint task with 30 s of rest between sets was performed, with participants positioned 1 m behind the photocells. The same active and rest measurements (in this order) were conducted after the repeated sprint task, followed by two MVIC at the end of the protocol.

2.4. Data Processing

Prior selection of the sonograms matching the criteria explained below was performed for each video recorded. All sonograms were digitized using ImageJ software (NIH, 1.47v, Bethesda, MD, USA), and FL, PA and MT were determined from each sonogram. The researcher was blinded to the outcome during the digitization process. As only part of fascicle could be visualized during static-image assessments, the linear extrapolation method was used to estimate the non-visible portion of the fascicle (15) using the following Equation (1):
FL = L + (h/sinβ)
where L is the observable fascicle length from the mid-muscle aponeurosis to the most visible endpoint, h is the perpendicular distance between the superficial aponeurosis and the fascicle’s visible distal end-point and β is the angle between the fascicle (drawn linearly to the most distal point) and the superficial aponeurosis, as in Figure 1. The MT was measured as the distance between the superficial and the mid-muscle aponeurosis measured at the ROI site. For each condition, two videos were recorded, and a frame was selected from each video. This was conducted in pre- and post-task instants (totaling 8 images for each limb). Three fascicles were selected in each image, meeting the criterion that a reasonable portion of the fascicle (~25% or more of the total estimated length) must have been visible within the ultrasound transducer’s field of view to be included in the present study [31].

2.5. Statistical Analysis

Data are presented as mean ± standard deviation. Data analysis was performed using IBM SPSS Statistics 27.0 (IBM Corporation, Armonk, NY, USA). Normality of data distribution was confirmed using the Shapiro–Wilk test. A paired sample t-test was performed for the descriptive statistics: age, height, and weight were calculated and comparisons between players with HSI history and healthy controls were made. An inter-individual comparison was performed between injured players and healthy controls. With respect to the effects of the sprint task, their interaction was examined by conducting a two-way mixed repeated measures ANOVA [injury × instant (pre and post)] for the variable average maximal sprint speed, FL, PA and MT in rest and active state. To determine the differences in femur length between the two groups, an independent t-test was used. An intra-individual comparison was performed between injured and contralateral limbs to compare the effects of the sprint task, using a two-way repeated measures ANOVA [injury × instant (pre and post)].

3. Results

No significant differences were observed for age (injured players: 25.79 ± 3.70; healthy controls: 25.79 ± 3.56 years; p = 1.000), height (injured players: 1.76 ± 0.06; healthy controls: 1.81 ± 0.09 m; p = 0.208) or weight (injured players: 73.53 ± 9.97; healthy controls: 76.21 ± 9.28 kg; p = 0.472). With respect to the femur length, no significant differences were seen between the previously injured players and healthy controls (injured players: 43.3 ± 1.75 cm; healthy controls: 43.19 ± 2.04 cm; p = 0.86; d = 0.06).
Regarding the sprint performance, a significant decrease was seen for the average maximal sprint speed (p < 0.001; η2p = 0.471) with a significantly higher speed for the injured players (injured: 7.00 ± 0.25 m/s; healthy controls: 6.82 ± 0.19 m/s; p = 0.026; η2p = 0.155); however, no interaction was seen for instant × injury (p = 0.442; η2p = 0.029) (Figure 2). For the comparison between injured players and healthy controls, the sprint task effects on BFlh FL, PA, MT and PT are shown in Table 1, reporting the ANOVA statistical analysis. In the resting state, a shorter FL (p = 0.014; η2p = 0.196) and greater PA (p = 0.002; η2p = 0.307) were seen for the previously injured players, with only the PA showing a significant increase between instants (p = 0.011; η2p = 0.208), as well as an instant × injury interaction (p = 0.041; η2p = 0.141). In the active state, a lower FL (p < 0.001; η2p = 0.413) and PA (p < 0.001; η2p = 0.368) were seen for the injured players. The peak torque produced in the MVIC demonstrated a significant decrease between instants (p < 0.001; η2p = 0.399), but not for the injury factor, nor an instant × injury interaction.
For the comparison between injured vs. contralateral limbs, the sprinting effects on BFlh FL, PA and MT and peak torque are shown in Table 2, reporting the ANOVA statistical analysis. In the resting state, a lower FL was demonstrated for injured limbs (p = 0.012; η2p = 0.519) with a decrease for both limb groups between instants (p < 0.001; η2p = 0.799). The PA also presented a significant increase for both groups between instants (p < 0.001; η2p = 0.732) with higher values for injured limbs (p = 0.006; η2p = 0.582). In the active state, a lower FL was seen for injured limbs compared to contralateral limbs (p = 0.039; η2p = 0.332) with no differences between instants nor instant × injury interaction. A higher PA was also shown for injured limbs (p = 0.018; η2p = 0.412). The peak torque only showed a significant difference for the instant factor (p < 0.017; η2p = 0.419).

4. Discussion

In this study we examined BFlh architecture in passive (i.e., at rest) and active (i.e., at 20% of MVIC) conditions before and after a repeated sprint protocol in professional football players with and without previous HSI in the BFlh. To the best of our knowledge, this is the first retrospective study to examine the acute effects of a repeated sprint protocol on BFlh architecture in professional football players. The main findings were: (i) injury players have shorter FL and greater PA in passive and active conditions in relation to controls; (ii) injury limbs have shorter FL and greater PA in passive and active conditions in relation to their contralateral limb; (iii) the sprint task had limited impact on the BFlh architectural parameters.
Firstly, in relation to sprint performance, significant differences were seen between groups, with the injured football players displaying higher average maximal speed than the healthy control group. These results are in concordance with a previous study that reported higher maximal mean speed in using a protocol with the same distance, although with 6 repetitions and 90 s of rest between sets instead of the 10 repetitions and 30 s of rest between sets in this study [32]. On other hand, Røksund et al. (2017) reported that football players with a previous HSI showed a higher drop in speed during a repeated sprint test (8 × 20 m) with a slightly lower, but not significant, sprint time (5.25 s vs. 5.35 s) during one maximal sprint of 40-m. One possible reason for these contradictory results could be the heterogeneity between studies due to different protocols and samples, recovery, return to play and screening periods. Therefore, this warrants further investigation.
Regarding the PT, no significant differences were seen between injured and uninjured players, which agrees with a previous study that demonstrated (based on effect sizes and confidence intervals) similar short-length (hip = 0°, knee = 90°) isometric strength between previously injured limbs and uninjured contralateral limbs at assessment points beyond 7 days post-injury [33]. These results are similar in a meta-regression analysis of long-length (hip = 0°, knee = 0–15°) isometric strength, with no statistical differences between the two groups at assessment points beyond 7 days post-injury [34]. Despite the different protocol positions compared with our study (the hip at 0° and knee flexed at 30°), the results are consistent, indicating that time is fundamental for the recovery of isometric strength. Meanwhile, significant differences have been demonstrated between previously injured and uninjured contralateral limbs regarding their eccentric contraction strength during Nordic exercises [34]. This was further supported in a prospective study that showed no differences in isometric contraction strength but demonstrated differences in eccentric contraction strength between later-injured and uninjured players [15]. It should be noted that isometric contractions are extremely different from eccentric ones. Therefore, it is possible for the neuromuscular system to not be altered between these groups in isometric contractions and yet demonstrate differences in eccentric contractions. Interestingly, most injuries in the aforementioned prospective study occurred during the sprint (82%) in the BFlh (88%). Sprinting involves eccentric contractions in the hamstring group [15,35], and the BFlh is one of the greatest contributors in the late swing phase (where these injuries tend to occur more often [35]), it being relatively more eccentrically stretched in this phase [36]; this could be an explanation for the increased sprint-related BFlh injury among participants demonstrating statistically lower eccentric strength prior to injury. This warrants further investigation.
Furthermore, during eccentric actions, it was verified that the distance between the aponeuroses (averaged at various angular velocities) at a constant ankle joint angle was significantly greater compared with concentric muscle actions [37], this distance being possibly related to the reported smaller pennation angle and longer fascicles. It might be possible that with greater pennation angles and shorter fascicles, the distance between the aponeuroses would be even greater according to trigonometric equations. If so, the force produced by the fascicles would not be transmitted as efficiently given the orientation between the fascicles relative to the aponeuroses. However, these effects could be muscle dependent, and no differences were seen in the PT between all the groups, which can indicate that BFlh architecture is not a preponderant factor for isometric knee flexion torque production. Indeed, it should be noted that BFlh has a greater hip moment arm compared to the knee moment arm, being more sensitive to mobilization of hip flexion [38]; on the other hand, semitendinosus is being consistently identified as the greatest contributor to knee flexion [39,40,41]. Therefore, it would be interesting to analyze the behavior of the sprint protocol in hamstrings PT measured in hip extension. Future researchers should explore this approach.
Regarding the muscle architecture analysis, and in accordance with our initial hypothesis, shorter FL and greater PA were observed in the injured group when compared with both the healthy controls and contralateral limbs, thus indicating that these results are inter- and intra-individual dependent. The implications of shorter fascicles have also been studied, demonstrating that shorter fascicles can affect the joint angle where optimal force is produced during activity [42,43] and reduce the contraction velocity [44,45]. During repeated sprint running, hamstrings reach relatively long fascicle lengths. Thus, at relatively longer muscle lengths, shorter BFlh fascicles would present lower force-generating capacity compared with longer BFlh fascicles; however, it should be noted that in present study, previously injured players ran at higher speed than the uninjured, and higher FL was seen in the contralateral group. Since the study is retrospective, it can be suggested that the higher speed of the previously injured players (greater excursion and velocity capacity of the sarcomeres) and the higher FL in the contralateral limb favors the hypothesis that the reduced FL in previously injured players is a consequence of the injury.
With respect to the comparison to the contralateral limb, these results are in line with the literature, with previous studies reporting shorter FL in the BFlh of players or athletes with previous injury [15,46,47]. Timmins et al. (2015) showed shorter BFlh fascicles in previously injured muscles compared with the contralateral uninjured muscles using a retrospective design, as in our study. However, the greater contralateral FL in the injured group indicates that it is not possible to determine whether this characteristic is a consequence of the injury or is a factor which predisposes to a greater risk of injury. Lima-E-Silva (2020) also reported that in professional senior and under-20 players, injured (in the previous season) players presented shorter BFlh FL compared with players without injury history. Additionally, injured players showed shorter BFlh FL in both limbs when compared with the threshold of 10.56 cm [15], although BFlh FL was symmetrical between limbs. Ribeiro-Alvares et al. (2020) also studied the prevalence of HSI risk factors in professional and under-20 male football players, finding that almost half of the sample of players had short BFlh relative FL (fascicle length/muscle length < 0.25) [15] and 30% of injury prevalence [47]. Overall, the greater contralateral FL in the present study possibly indicates that the shorter fascicles of football players can be a consequence of the injury and not a predictive measure.
In relation to the specificity of the task, in the present study, there was limited impact on the BFlh architectural parameters; specifically, a significant increase was seen only for PA in (i) the previously injured group in comparison with the healthy control group, and (ii) in the comparison between injured and contralateral limbs, via an exercise-induced decrease in FL for the contralateral with an increase in PA for both limbs in resting conditions. Repeated sprint exercise seems to decrease the net moments around the hip and knee, and the knee and hip peak joint angles [48]. Special consideration must be given to the fatigability associated with the repeated loading, as this has been suggested to increase the contractile compliance (decreasing muscle stiffness) as a result of repeated damage to the muscle with no changes in tendon compliance [49,50], suggesting that exercise induced changes in muscle-tendon function may lead to a greater strain in the BFlh. Moreover, muscle damage might be greater for individuals with shorter fascicle lengths, assuming the same (force) stress applied will induce a higher strain in BFlh. Therefore, these mechanical alterations due to lower FL are possible causes that place previously injured players at a higher injury risk of sustaining other HSI (with previous HSI being the most consistent risk factor), even though PT was similar in the uninjured individuals. However, no muscle and tendon stiffness was not measured and no differences were seen in FL after the repeated sprint task. Since the differences between instants are within the minimal detectable changes of the static technique (0.84 cm for FL and 1.5° for PA) [51] architectural changes following repeated sprints could be only statistically but not physiologically relevant. Moreover, it should be noted that while the task induced increases in PA and decreases in FL, the groups compared were not differentially affected by the task; the magnitude of the differences in FL and PA before the task were not statistically altered. This warrants further investigation.
Furthermore, we assessed the BFlh architecture using static technique; this technique allows for measurement during contractions. However, it is important to understand some methodological considerations when assessing muscle architecture, especially of the BFlh, which is heterogeneous, with the fascicles typically following a concave-to-convex path being differentially oriented along the muscle length. The FL estimation equation does not consider the different orientation of the fascicle, which might result in over- (for concave fascicles) or under-estimation (for convex fascicles) of the FL [28], implying that the assumptions with this type of methodology could affect the interpretation of the results. Additionally, it has been reported that transducer width is important in the linear extrapolation of the FL. Indeed, a 1 cm reduction in the transducer width would imply that approximately 11% of the FL would need to be estimated, which indicates an overestimation of the fascicle [28]. In our case, the probe used had 3.8 cm of width, which could overestimate ~24.2% of the FL when compared with Pimenta et al. (2018). It should also be noted that the visible portion of the fascicle in the linear extrapolation technique represents an extremely important factor in the equation that could increase or decrease the overestimation of the FL.
This study has some limitations. Firstly, these results can only be interpreted for senior professional male football players, since different results have been shown for different age groups [52], and the anatomical and biomechanical characteristics of running in females are different from those in males [53,54]. Secondly, the static technique has some limitations with respect to the extrapolation of the non-visible part of the fascicle, which is determined by the width of the transducer [28] and the equation used, which in this case overestimates FL by 0.5 cm, particularly with real FL above ~7 cm [31]; thus, these results should be analyzed carefully. Fascicle length will also depend on the femur length, where subjects with a larger femur would be expected to present larger fascicles. However, no significant differences were seen between groups, indicating that an anatomical criterion did not explain the differences between groups.

5. Conclusions

The present study provides evidence that injured football players in the BFlh are faster than the healthy control group, without difference in a mechanical parameter as PT. Regarding muscle architecture, in resting condition, injured players display shorter FL and greater PA than contralateral and healthy control limbs after a repeated sprint task. The same architectural differences were detected in the active condition after the repeated sprint task. Regarding the specificity of a repeated sprint exercise bout, the current study does not provide sufficient evidence to suggest the impact of the task on BFlh architectural parameters, nor that the task differentially impacts previously injured and uninjured groups. Overall, while one must take care when comparing architectural parameters between the limbs of different subjects; the FL differences observed between previously injured limbs and uninjured contralateral ones in both resting and active conditions could suggest the shorter FL to be a consequence of the injury rather than a predictive measure, as well as potentially indicating an increased risk of re-injury in the same limb. In fact, a history of HSI has been shown to be the most consistent risk factor for sustaining an HSI.

Author Contributions

R.P. contributed to the conception and design of the study, collected all the data, performed the statistical analysis and interpretation of the data, produced the first draft of the manuscript, reviewed the article critically for important intellectual content and edited subsequent versions of the manuscript. H.A. contributed to the processing of the data, manuscript editing and improved intellectual content. T.L. contributed to the data collection and processed all data, and manuscript editing. A.V. participated in the interpretation of the data, editing, and reviewing of the manuscript critically for important intellectual content. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by CIPER–Centro Interdisciplinar para o Estudo da Performance Humana (unit 447) which is also supported by Fundação para a Ciência e Tecnologia (grant: UID/DTP/ UI447/2019). R.P. is funded by a scholarship from the University of Lisbon (number 793).

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki and approved by the Ethics Committee of Faculty of Human Kinetics at University of Lisbon (#5/2021).

Informed Consent Statement

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

Data Availability Statement

Not applicable.

Acknowledgments

The authors thank the clubs and athletes for their participation, to Sport Lisboa e Benfica for lending one Aixplorer unit and to the Centro de Alto Rendimento Jamor for providing the sports facilities.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gabbe, B.J.; Bennell, K.L.; Finch, C.F.; Wajswelner, H.; Orchard, J.W. Predictors of Hamstring Injury at the Elite Level of Australian Football. Scand. J. Med. Sci. Sport. 2006, 16, 7–13. [Google Scholar] [CrossRef] [PubMed]
  2. Gabbe, B.; Finch, C.; Wajswelner, H.; Bennell, K. Australian Football: Injury Profile at the Community Level. J. Sci. Med. Sport 2002, 5, 149–160. [Google Scholar] [CrossRef] [PubMed]
  3. Feeley, B.T.; Kennelly, S.; Barnes, R.P.; Muller, M.S.; Kelly, B.T.; Rodeo, S.A.; Warren, R.F. Epidemiology of National Football League Training Camp Injuries from 1998 to 2007. Am. J. Sport. Med. 2008, 36, 1597–1603. [Google Scholar] [CrossRef] [PubMed]
  4. Brooks, J.H.M.; Fuller, C.W.; Kemp, S.P.T.; Reddin, D.B. Incidence, Risk, and Prevention of Hamstring Muscle Injuries in Professional Rugby Union. Am. J. Sport. Med. 2006, 34, 1297–1306. [Google Scholar] [CrossRef] [PubMed]
  5. Henderson, G.; Barnes, C.A.; Portas, M.D. Factors Associated with Increased Propensity for Hamstring Injury in English Premier League Soccer Players. J. Sci. Med. Sport 2010, 13, 397–402. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Ekstrand, J.; Hägglund, M.; Waldén, M. Injury Incidence and Injury Patterns in Professional Football: The UEFA Injury Study. Br. J. Sport. Med. 2011, 45, 553–558. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Woods, C.; Hawkins, R.D.; Maltby, S.; Hulse, M.; Thomas, A.; Hodson, A. Football Association Medical Research Programme The Football Association Medical Research Programme: An Audit of Injuries in Professional Football—Analysis of Hamstring Injuries. Br. J. Sport. Med. 2004, 38, 36–41. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Ekstrand, J.; Hägglund, M.; Waldén, M. Epidemiology of Muscle Injuries in Professional Football (soccer). Am. J. Sport. Med. 2011, 39, 1226–1232. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Hägglund, M.; Waldén, M.; Ekstrand, J. Risk Factors for Lower Extremity Muscle Injury in Professional Soccer: The UEFA Injury Study. Am. J. Sport. Med. 2013, 41, 327–335. [Google Scholar] [CrossRef] [Green Version]
  10. Heiderscheit, B.C.; Sherry, M.A.; Silder, A.; Chumanov, E.S.; Thelen, D.G. Hamstring Strain Injuries: Recommendations for Diagnosis, Rehabilitation, and Injury Prevention. J. Orthop. Sport. Phys. Ther. 2010, 40, 67–81. [Google Scholar] [CrossRef]
  11. Ekstrand, J. Keeping Your Top Players on the Pitch: The Key to Football Medicine at a Professional Level. Br. J. Sport. Med. 2013, 47, 723–724. [Google Scholar] [CrossRef] [Green Version]
  12. Green, B.; Bourne, M.N.; van Dyk, N.; Pizzari, T. Recalibrating the Risk of Hamstring Strain Injury (HSI): A 2020 Systematic Review and Meta-Analysis of Risk Factors for Index and Recurrent Hamstring Strain Injury in Sport. Br. J. Sport. Med. 2020, 54, 1081–1088. [Google Scholar] [CrossRef] [PubMed]
  13. Girard, O.; Mendez-Villanueva, A.; Bishop, D. Repeated-Sprint Ability—Part I: Factors Contributing to Fatigue. Sport. Med. 2011, 41, 673–694. [Google Scholar] [CrossRef]
  14. Askling, C.M.; Tengvar, M.; Saartok, T.; Thorstensson, A. Acute First-Time Hamstring Strains during High-Speed Running: A Longitudinal Study Including Clinical and Magnetic Resonance Imaging Findings. Am. J. Sport. Med. 2007, 35, 197–206. [Google Scholar] [CrossRef] [PubMed]
  15. Timmins, R.G.; Bourne, M.N.; Shield, A.J.; Williams, M.D.; Lorenzen, C.; Opar, D.A. Short Biceps Femoris Fascicles and Eccentric Knee Flexor Weakness Increase the Risk of Hamstring Injury in Elite Football (soccer): A Prospective Cohort Study. Br. J. Sport. Med. 2016, 50, 1524–1535. [Google Scholar] [CrossRef] [PubMed]
  16. Porter, K.; Opar, D.; Williams, M.; Shield, A. Biceps Femoris Fascicle Length Is Shorter in a Previously Hamstring Injured Athlete. J. Sci. Med. Sport 2012, 15, S152. [Google Scholar] [CrossRef]
  17. Timmins, R.G.; Shield, A.J.; Williams, M.D.; Lorenzen, C.; Opar, D.A. Biceps Femoris Long Head Architecture: A Reliability and Retrospective Injury Study. Med. Sci. Sport. Exerc. 2015, 47, 905–913. [Google Scholar] [CrossRef] [Green Version]
  18. Kellis, E.; Galanis, N.; Kapetanos, G.; Natsis, K. Architectural Differences between the Hamstring Muscles. J. Electromyogr. Kinesiol. 2012, 22, 520–526. [Google Scholar] [CrossRef]
  19. Lieber, R.L.; Fridén, J. Functional and Clinical Significance of Skeletal Muscle Architecture. Muscle Nerve 2000, 23, 1647–1666. [Google Scholar] [CrossRef]
  20. Butterfield, T.A. Eccentric Exercise in Vivo: Strain-Induced Muscle Damage and Adaptation in a Stable System. Exerc. Sport Sci. Rev. 2010, 38, 51–60. [Google Scholar] [CrossRef]
  21. Silder, A.; Thelen, D.G.; Heiderscheit, B.C. Effects of Prior Hamstring Strain Injury on Strength, Flexibility, and Running Mechanics. Clin. Biomech. 2010, 25, 681–686. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Askling, C.; Saartok, T.; Thorstensson, A. Type of Acute Hamstring Strain Affects Flexibility, Strength, and Time to Return to Pre-Injury Level. Br. J. Sport. Med. 2006, 40, 40–44. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Askling, C.M.; Tengvar, M.; Saartok, T.; Thorstensson, A. Acute First-Time Hamstring Strains during Slow-Speed Stretching: Clinical, Magnetic Resonance Imaging, and Recovery Characteristics. Am. J. Sport. Med. 2007, 35, 1716–1724. [Google Scholar] [CrossRef]
  24. Kwah, L.K.; Pinto, R.Z.; Diong, J.; Herbert, R.D. Reliability and Validity of Ultrasound Measurements of Muscle Fascicle Length and Pennation in Humans: A Systematic Review. J. Appl. Physiol. 2013, 114, 761–769. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Wisdom, K.M.; Delp, S.L.; Kuhl, E. Use It or Lose It: Multiscale Skeletal Muscle Adaptation to Mechanical Stimuli. Biomech. Model. Mechanobiol. 2015, 14, 195–215. [Google Scholar] [CrossRef]
  26. Guex, K.; Degache, F.; Morisod, C.; Sailly, M.; Millet, G.P. Hamstring Architectural and Functional Adaptations following Long vs. Short Muscle Length Eccentric Training. Front. Physiol. 2016, 7, 340. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Silder, A.; Whittington, B.; Heiderscheit, B.; Thelen, D.G. Identification of Passive Elastic Joint Moment-Angle Relationships in the Lower Extremity. J. Biomech. 2007, 40, 2628–2635. [Google Scholar] [CrossRef] [Green Version]
  28. Pimenta, R.; Blazevich, A.J.; Freitas, S.R. Biceps Femoris Long-Head Architecture Assessed Using Different Sonographic Techniques. Med. Sci. Sport. Exerc. 2018, 50, 2584–2594. [Google Scholar] [CrossRef]
  29. Berg, H.E.; Tedner, B.; Tesch, P.A. Changes in Lower Limb Muscle Cross-Sectional Area and Tissue Fluid Volume after Transition from Standing to Supine. Acta Physiol. Scand. 1993, 148, 379–385. [Google Scholar] [CrossRef]
  30. Kositsky, A.; Gonçalves, B.A.M.; Stenroth, L.; Barrett, R.S.; Diamond, L.E.; Saxby, D.J. Reliability and Validity of Ultrasonography for Measurement of Hamstring Muscle and Tendon Cross-Sectional Area. Ultrasound Med. Biol. 2020, 46, 55–63. [Google Scholar] [CrossRef]
  31. Franchi, M.V.; Fitze, D.P.; Raiteri, B.J.; Hahn, D.; Spörri, J. Ultrasound-Derived Biceps Femoris Long Head Fascicle Length: Extrapolation Pitfalls. Med. Sci. Sport. Exerc. 2020, 52, 233–243. [Google Scholar] [CrossRef] [Green Version]
  32. Ishøi, L.; Thorborg, K.; Hölmich, P.; Krommes, K. Sprint Performance In Football (Soccer) Players With And Without A Previous Hamstring Strain Injury: An Explorative Cross-Sectional Study. Int. J. Sport. Phys. Ther. 2020, 15, 947–957. [Google Scholar] [CrossRef] [PubMed]
  33. Reurink, G.; Goudswaard, G.J.; Moen, M.H.; Weir, A.; Verhaar, J.A.N.; Bierma-Zeinstra, S.M.A.; Maas, M.; Tol, J.L. Rationale, Secondary Outcome Scores and 1-Year Follow-up of a Randomised Trial of Platelet-Rich Plasma Injections in Acute Hamstring Muscle Injury: The Dutch Hamstring Injection Therapy Study. Br. J. Sport. Med. 2015, 49, 1206–1212. [Google Scholar] [CrossRef] [PubMed]
  34. Maniar, N.; Shield, A.J.; Williams, M.D.; Timmins, R.G.; Opar, D.A. Hamstring Strength and Flexibility after Hamstring Strain Injury: A Systematic Review and Meta-Analysis. Br. J. Sport. Med. 2016, 50, 909–920. [Google Scholar] [CrossRef] [PubMed]
  35. Chumanov, E.S.; Heiderscheit, B.C.; Thelen, D.G. The Effect of Speed and Influence of Individual Muscles on Hamstring Mechanics during the Swing Phase of Sprinting. J. Biomech. 2007, 40, 3555–3562. [Google Scholar] [CrossRef]
  36. Baumert, P.; Temple, S.; Stanley, J.M.; Cocks, M.; Strauss, J.A.; Shepherd, S.O.; Drust, B.; Lake, M.J.; Stewart, C.E.; Erskine, R.M. Neuromuscular Fatigue and Recovery after Strenuous Exercise Depends on Skeletal Muscle Size and Stem Cell Characteristics. Sci. Rep. 2021, 11, 7733. [Google Scholar] [CrossRef]
  37. Reeves, N.D.; Narici, M.V. Behavior of Human Muscle Fascicles during Shortening and Lengthening Contractions in Vivo. J. Appl. Physiol. 2003, 95, 1090–1096. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Sant, G.L.; Le Sant, G.; Ates, F.; Brasseur, J.-L.; Nordez, A. Elastography Study of Hamstring Behaviors during Passive Stretching. PLoS ONE 2015, 10, e0139272. [Google Scholar]
  39. Bourne, M.N.; Williams, M.D.; Opar, D.A.; Al Najjar, A.; Kerr, G.K.; Shield, A.J. Impact of Exercise Selection on Hamstring Muscle Activation. Br. J. Sport. Med. 2017, 51, 1021–1028. [Google Scholar] [CrossRef] [Green Version]
  40. Bourne, M.N.; Timmins, R.G.; Opar, D.A.; Pizzari, T.; Ruddy, J.D.; Sims, C.; Williams, M.D.; Shield, A.J. An Evidence-Based Framework for Strengthening Exercises to Prevent Hamstring Injury. Sport. Med. 2018, 48, 251–267. [Google Scholar] [CrossRef]
  41. Evangelidis, P.E.; Shan, X.; Otsuka, S.; Yang, C.; Yamagishi, T.; Kawakami, Y. Hamstrings Load Bearing in Different Contraction Types and Intensities: A Shear-Wave and B-Mode Ultrasonographic Study. PLoS ONE 2021, 16, e0251939. [Google Scholar] [CrossRef] [PubMed]
  42. Frisk, R.F.; Lorentzen, J.; Barber, L.; Nielsen, J.B. Characterization of Torque Generating Properties of Ankle Plantar Flexor Muscles in Ambulant Adults with Cerebral Palsy. Eur. J. Appl. Physiol. 2019, 119, 1127–1136. [Google Scholar] [CrossRef]
  43. Gao, F.; Zhang, L.-Q. Altered Contractile Properties of the Gastrocnemius Muscle Poststroke. J. Appl. Physiol. 2008, 105, 1802–1808. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Moreau, N.G.; Falvo, M.J.; Damiano, D.L. Rapid Force Generation Is Impaired in Cerebral Palsy and Is Related to Decreased Muscle Size and Functional Mobility. Gait Posture 2012, 35, 154–158. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Moreau, N.G.; Holthaus, K.; Marlow, N. Differential Adaptations of Muscle Architecture to High-Velocity versus Traditional Strength Training in Cerebral Palsy. Neurorehabil. Neural Repair 2013, 27, 325–334. [Google Scholar] [CrossRef] [PubMed]
  46. De Lima-E-Silva, F.X.; de Lima-E-Silva, F.X.; Oliveira, G.S.; Medeiros, T.M.; Dornelles, M.P.; Ribeiro-Alvares, J.B.A.; Baroni, B.M. Flexibility, Strength, and Fascicle Length of Football Players with and without History of Hamstring Strain Injury in the Prior Season. Sci. Med. Footb. 2020, 4, 322–328. [Google Scholar] [CrossRef]
  47. Ribeiro-Alvares, J.B.; Dornelles, M.P.; Fritsch, C.G.; de Lima-E-Silva, F.X.; Medeiros, T.M.; Severo-Silveira, L.; Marques, V.B.; Baroni, B.M. Prevalence of Hamstring Strain Injury Risk Factors in Professional and Under-20 Male Football (Soccer) Players. J. Sport Rehabil. 2020, 29, 339–345. [Google Scholar] [CrossRef]
  48. Gonçalves, B.A.M.; Meinders, E.; Saxby, D.J.; Barrett, R.S.; Bourne, M.N.; Diamond, L.E. Repeated Sprints Alter Mechanical Work Done by Hip and Knee, but Not Ankle, Sagittal Moments. J. Sci. Med. Sport 2021, 24, 939–944. [Google Scholar] [CrossRef]
  49. Butterfield, T.A.; Herzog, W. Quantification of Muscle Fiber Strain during in Vivo Repetitive Stretch-Shortening Cycles. J. Appl. Physiol. 2005, 99, 593–602. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Kubo, K.; Ikebukuro, T. Changes in Joint, Muscle, and Tendon Stiffness Following Repeated Hopping Exercise. Physiol Rep 2019, 7, e14237. [Google Scholar] [CrossRef]
  51. Freitas, S.R.; Marmeleira, J.; Valamatos, M.J.; Blazevich, A.; Mil-Homens, P. Ultrasonographic Measurement of the Biceps Femoris Long-Head Muscle Architecture. J. Ultrasound Med. 2018, 37, 977–986. [Google Scholar] [CrossRef] [PubMed]
  52. Lee Dow, C.; Timmins, R.G.; Ruddy, J.D.; Williams, M.D.; Maniar, N.; Hickey, J.T.; Bourne, M.N.; Opar, D.A. Prediction of Hamstring Injuries in Australian Football Using Biceps Femoris Architectural Risk Factors Derived From Soccer. Am. J. Sport. Med. 2021, 49, 3687–3695. [Google Scholar] [CrossRef] [PubMed]
  53. Ferber, R.; Davis, I.M.; Williams, D.S., 3rd. Gender Differences in Lower Extremity Mechanics during Running. Clin. Biomech. 2003, 18, 350–357. [Google Scholar] [CrossRef]
  54. Willson, J.D.; Petrowitz, I.; Butler, R.J.; Kernozek, T.W. Male and Female Gluteal Muscle Activity and Lower Extremity Kinematics during Running. Clin. Biomech. 2012, 27, 1052–1057. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Biceps femoris long head (BFlh) architecture, measured at pre-task active condition. L, observable fascicle length; h, height; β, pennation angle; MT, muscle thickness; FL, fascicle length.
Figure 1. Biceps femoris long head (BFlh) architecture, measured at pre-task active condition. L, observable fascicle length; h, height; β, pennation angle; MT, muscle thickness; FL, fascicle length.
Biology 12 00096 g001
Figure 2. Acute effects of a sprint protocol on the average maximal speed injured players vs. healthy controls. Data are presented as mean ± standard deviation. * Significant difference between the current sprint and each subsequent sprint indicated.
Figure 2. Acute effects of a sprint protocol on the average maximal speed injured players vs. healthy controls. Data are presented as mean ± standard deviation. * Significant difference between the current sprint and each subsequent sprint indicated.
Biology 12 00096 g002
Table 1. Acute effects of a sprint protocol on BFlh injured players vs. healthy controls on: fascicle length (FL), pennation angle (PA), muscle thickness (MT) and peak torque (PT) before and after the sprint task. Data are presented as mean ± standard deviation.
Table 1. Acute effects of a sprint protocol on BFlh injured players vs. healthy controls on: fascicle length (FL), pennation angle (PA), muscle thickness (MT) and peak torque (PT) before and after the sprint task. Data are presented as mean ± standard deviation.
Injured Players vs. Healthy Controls
PRE SPRINTPOST SPRINTInstant
(Pre vs. Post)
Injury
(Injured vs. Control)
Instant × Injury Interaction
InjuredControlInjuredControlpη2ppη2ppη2p
RESTFL (cm)7.40 ± 0.988.20 ± 0.807.02 ± 1.208.29 ± 1.380.2980.0390.0140.1960.1030.092
PA (°)21.35 ± 2.6718.47 ± 2.3222.91 ± 3.0918.65 ± 3.000.0110.2080.0020.3070.0410.141
MT (cm)2.61 ± 0.282.58 ± 0.302.66 ± 0.262.60 ± 0.330.2080.0560.7070.0050.6190.009
ACTIVEFL (cm)6.36 ± 1.217.89 ± 0.966.22 ± 1.137.86 ± 0.830.5040.015<0.0010.4130.6770.006
PA (°)25.04 ± 3.3420.28 ± 2.9724.68 ± 4.0620.60 ± 2.700.971<0.001<0.0010.3680.4920.016
MT (cm)2.62 ± 0.272.69 ± 0.30 2.54 ± 0.242.73 ± 0.260.4760.0170.1730.0610.0360.138
PT (N⋅m)141.17 ± 43.90149.39 ± 35.45132.54 ± 37.18136.64 ± 32.27<0.0010.3990.6420.0070.3970.024
Abbreviations: FL, Fascicle length, PA, Pennation angle, MT, Muscle Thickness, PT, Peak Torque, p, p value, η2p, eta squared.
Table 2. Acute effects of a sprint protocol on BFlh injured vs. contralateral limbs on: fascicle length (FL), pennation angle (PA), muscle thickness (MT) and peak torque (PT) before and after sprint task. Data are presented as mean ± standard deviation.
Table 2. Acute effects of a sprint protocol on BFlh injured vs. contralateral limbs on: fascicle length (FL), pennation angle (PA), muscle thickness (MT) and peak torque (PT) before and after sprint task. Data are presented as mean ± standard deviation.
Injured vs. Contralateral Limbs
PRE SPRINTPOST SPRINTInstant
(Pre vs. Post)
Injury
(Injured vs. Healthy)
Instant × Injury Interaction
InjuredContralateralInjuredContralateralpη2ppη2ppη2p
RESTFL (cm)7.40 ± 0.989.14 ± 2.187.02 ± 1.208.36 ± 1.88<0.0010.7990.0120.5190.3940.082
PA (°)21.35 ± 2.6718.81 ± 4.2222.91 ± 3.0920.14 ± 3.620.0010.7320.0060.5820.8340.005
MT (cm)2.61 ± 0.282.77 ± 0.242.66 ± 0.262.75 ± 0.240.6460.0240.0880.2890.2500.144
ACTIVEFL (cm)6.36 ± 1.217.73 ± 2.576.22 ± 1.137.60 ± 1.840.5430.0350.0390.3320.989<0.001
PA (°)25.04 ± 3.3422.06 ± 4.2824.68 ± 4.0621.23 ± 3.460.4540.0520.0180.4120.6910.015
MT (cm)2.62 ± 0.272.70 ± 0.40 2.54 ± 0.242.60 ± 0.350.0730.2630.3710.0730.6850.016
PT (N⋅m)141.17 ± 43.90141.44 ± 36.59132.54 ± 37.18133.76 ± 35.530.0170.4190.8570.0030.7350.011
Abbreviations: FL, Fascicle length, PA, Pennation angle, MT, Muscle Thickness, PT, Peak Torque, p, p value, η2p, eta squared.
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

Pimenta, R.; Antunes, H.; Lopes, T.; Veloso, A. Do Repeated Sprints Affect the Biceps Femoris Long Head Architecture in Football Players with and without an Injury History?—A Retrospective Study. Biology 2023, 12, 96. https://doi.org/10.3390/biology12010096

AMA Style

Pimenta R, Antunes H, Lopes T, Veloso A. Do Repeated Sprints Affect the Biceps Femoris Long Head Architecture in Football Players with and without an Injury History?—A Retrospective Study. Biology. 2023; 12(1):96. https://doi.org/10.3390/biology12010096

Chicago/Turabian Style

Pimenta, Ricardo, Hugo Antunes, Tomás Lopes, and António Veloso. 2023. "Do Repeated Sprints Affect the Biceps Femoris Long Head Architecture in Football Players with and without an Injury History?—A Retrospective Study" Biology 12, no. 1: 96. https://doi.org/10.3390/biology12010096

APA Style

Pimenta, R., Antunes, H., Lopes, T., & Veloso, A. (2023). Do Repeated Sprints Affect the Biceps Femoris Long Head Architecture in Football Players with and without an Injury History?—A Retrospective Study. Biology, 12(1), 96. https://doi.org/10.3390/biology12010096

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