- Original research article
- Open Access
Ventilatory response to exercise of elite soccer players
Multidisciplinary Respiratory Medicinevolume 9, Article number: 20 (2014)
The purpose of this study was to evaluate the role of ventilatory parameters in maximal exercise performance in elite soccer players.
From September 2009 to December 2012, 90 elite soccer players underwent evaluation of lung function test and ergospirometry by means of an incremental symptom-limited treadmill test. Results were analyzed according to i) maximal exercise velocity performed (Hi-M: high-performers, >18.65 km/h; Lo-M: low-performers, <18.65 km/h) and ii) usual role in the team.
Hi-M showed higher peak minute ventilation (
: 158.3 ± 19.5 vs 148.0 ± 18.54 L/min, p = 0.0203), and forced expiratory volume at first second (5.28 ± 0.50 vs 4.89 ± 0.52 liters, p < 0.001) than Lo-M, independently of playing role. Moreover, a significant correlation between peak oxygen uptake and
(r = 0.57, p < 0.001) was found.
Ventilatory response plays a role in the assessment of exercise capacity in elite soccer players.
Playing soccer is the result of combination of technical, tactical, psychological and athletic skills. To assess soccer player performance model, studies have used match analysis systems [1–5], and laboratory and field tests [6–9]. Other studies addressed different roles [10, 11], different data tests [12, 13], and relationships between endurance field tests and match analysis data .
Soccer is a sport characterized by more than 1,000 unpredictable and acyclic changes in activity, each occurring every 3 to 5 seconds, involving up to 40 sprints, tackles and jumps per match . Decelerations, kicking, dribbling, and tackling are additional actions required as well . Therefore, the physical effort imposed on the players is elicited by all these gestures: as a consequence soccer is a highly physiologically demanding sport, with additional stress resulting by frequent matches and high load training sessions are quite high [17–21].
Computerized systems for time-motion analysis have shown that during a match elite soccer players are able to perform 2 to 3 km of high-intensity running (at a speed >15 km/h) and about 0.6 km sprinting (>20 km/h): these distances are 28% and 58% respectively greater than those covered by intermediate level professional players . Technical and tactical skills in soccer are highly dependent on the player physical capacity [22, 23].
Relevant effort has been devoted to the evaluation of muscular strength, power, maximal speed, aerobic and anaerobic endurance, but poor attention has been addressed to ventilatory capacity. In most studies cardiorespiratory response to maximal exercise was evaluated using peak oxygen uptake () and the corresponding anaerobic threshold value [23–31]. To the best of our knowledge, currently available studies do not focus on ventilatory function as evaluation tool of professional soccer players performance.
In the present study we hypothesized that ventilatory capacity can be a determinant of soccer players exercise capacity. In order to demonstrate our hypothesis the relationship between ventilatory parameters and exercise tolerance was evaluated in a group of elite Italian soccer players using an incremental symptom-limited cardiopulmonary exercise test on a treadmill. Several different measurements reflecting cardiovascular, respiratory, and metabolic response can be recorded during the exercise test. The assessment of ventilatory parameters can contribute to improve training programmes and, consequently, athletic performance.
From 2009 to 2012, ninety professional soccer players from five Italian serie A soccer teams were evaluated in the period September-December of each year after completing pre-seasonal training programme in the frame of routine medical evaluations. Each subject provided informed signed consent to the use of their data for scientific purposes.
As expected, no player reported either smoking habit or any relevant disease, with negative chest physical examination. Players were unaware of the aim of the study and researchers performing analysis of results were blind to players’ identity.
Lung function test
Dynamic lung volumes were assessed by means of a pneumotachograph (V-Max Encore, Yorba Linda, CA, USA). Predicted values were those of American Thoracic Society . Maximal Voluntary Ventilation (MVV) was estimated multiplying Forced Expiratory Volume at first second (FEV1) value by 40 .
Resting and exercise electrocardiography (EKG) was assessed in upright position by means of a 10-lead electrocardiograph (Cardiosoft, GE medical systems, Fairfield, CT, USA) applied on the cardiac screening (six precordial leads) and the posterior wall of the chest (four peripheral leads).
An incremental symptom-limited exercise test was performed on a treadmill (Runrace 900, Technogym, Gambettola, Italy) under EKG and pulse oximetry monitoring. Subjects standing on the treadmill breathed through a mask. A continuous “ramp” protocol at constant grade (1%) (starting from 8 km/h, increasing speed by 1 km/h every 60 seconds) was used. The test was stopped when subjects complained of exhaustion. Exercise tolerance was evaluated as the maximal speed reached (Maximal Exercise Velocity: MEV), adjusted according to a modified Kuiper’s equation (Equation 1) .
where v l represents the speed achieved at the last exercise step and n the number of seconds attained during the last stage.
The following variables were measured at peak exercise through breath by breath analysis of inhaled and exhaled gases, by mass flow meter and fast-responding gas analyzer (V-Max Encore, Yorba Linda, CA, USA): Oxygen uptake () and normalized to body weight, and its relationship with heart rate (HR) (pulse oxygen or ), CO2 production (), the physiological dead space to the tidal volume ratio (V d /V t ), minute ventilation (), maximal () expressed as the highest value recorded either during exercise or at the first time recovery phase, breathing respiratory reserve (BRR% expressed as to MVV ratio). The anaerobic threshold (AT) was estimated by V slope method and ventilatory equivalent method . Anaerobic phase time (APT) was defined as the time spent during > 1. Predicted values of HR were computed according to Tanaka et al. .
Subjects were categorized into four groups according to their role in the team as reported by the technical staff: forwards (F), central midfielders (CM), central defenders (CD), wide players (WP): goalkeepers were excluded from the study. According to MEV performed in the exercise test, players were divided into two groups (Hi-M: high-performers: able to run at a speed greater than median value of all subjects; Lo-M: low-performers: able to run at a speed lower than that median). Furthermore, the subjects were divided into two groups according to median (Hi- and Lo- , respectively).
Linear regression analyses between and , HR and , and HR and MEV, MEV and were computed at peak exercise.
Two way analysis of variance (ANOVA) was carried out as follows: 1) dependent variable: ; source of variation: role, MEV; 2) dependent variable: FEV1; source of variation: role, MEV. All pairwise multiple comparison procedures were carried out using Holm-Sidak method.
Student’s t-test was used for analysing statistical significance of differences in the following parameters: , Body Mass Index (BMI), peak Respiratory Rate (RR peak ), peak tidal volume (Vt peak ), MEV and on Hi-M vs Lo-M and Hi-VE vs Lo-VE, respectively.
Factor analysis was computed on the set of cardiovascular, metabolic, and respiratory variables measured at peak exercise using non-rotated Principal Component Analysis (PCA). PCA is a simple, non-parametric method of extracting relevant information from multivariate datasets. The central idea of PCA is to reduce the dimensionality of a dataset consisting of a large number of interrelated variables, while retaining the variation present in the dataset. This is achieved by transforming into a new variable set, the Principal Components (PCs) which are uncorrelated, and which are ordered so that the first few retain most of the variation present in all the original variables .
Multiple linear regression among variables resulting from PCA analysis was carried out. Correlation analysis between variables of interest was computed using Pearson coefficient.
Statistical analysis was carried out using SigmaStat version 3.5 (Systat Software, Inc., USA), except for factor analysis carried out using SPSS v10.1 (SPSS Inc. USA).
Anthropometric, demographic, and resting physiological characteristics of the whole study population and according to exercise capacity are shown in Table 1. The MEV median value was 18.65 km/h. BMI, height, weight, and age did not show any significant differences between Hi-M and Lo-M. MVV and FEV1 were significantly higher in Hi-M than in Lo-M. There was no significant difference in anthropometric and demographic characteristics among different roles played.
Table 2 shows the physiological parameters at peak exercise according to exercise capacity. Only was significantly higher in Hi-M than Lo-M, resulting from non significantly greater Vtpeak and RRpeak, and non significantly lower Vd/Vt in Hi-M than Lo-M subjects.
During the test, all players reached 97.9 ± 4.7% of their predicted maximal HR without any significant difference between groups.
Table 3 shows resting and at peak exercise physiological characteristics according to . Median was 153.05 L/min. Hi-VE showed significantly higher , , MEV, BRR%, Vtpeak and FEV1 than Lo-VE subjects.
was significantly correlated to (r = 0.619, p = 0.001) (Figure 1) and MEV was significantly correlated to (r = 0.267, p = 0.011). No statistical significant correlation was found between HR peak and MEV.
A weak, although significant correlation between V d /V t and at peak exercise was found in all players (r = −0.282; p = 0.007). No significant differences in , MEV and BRR were found among different roles, whereas mean MVV was significantly lower in F than in CM (193.3 ± 21.2 vs 206.1 ± 19.9 liters respectively, p = 0.043). ANOVA showed that the difference in the mean value of and FEV1 among the different levels of role of players was not significantly different after allowing for the effects of differences in MEV (p = 0.377 and p = 0.543, respectively).
Factor analysis using the PCA on the variables measured at peak exercise and at rest was carried out. Table 4 reports the distribution along the first three components of the following variables: , , (V d /V t )peak, HR peak , BRRpeak, RRpeak, , FEV1. The first three components accounted for 74.33% of cumulative variance. Multiple linear regression on most representative variables on each component deriving from PCA at peak exercise is represented by the following equation:
Equation (2) represents a statistically significant relationship between VO2peak and VE peak , FEV1 and BRR.
The novelty of this study is represented by the analysis of ventilatory parameters as evaluation tool of professional soccer players performance. Our elite soccer players were evaluated during seasonal activity. The main finding of our study is that and FEV1 but not are the main determinants for discriminating high and low performers: nevertheless, neither Hi-M nor Lo-M were ventilatory limited. Therefore, our results are in agreement with available literature about elite athletes: the main limitation to maximal performance is represented by cardiovascular limit. Furthermore, our data highlight the key role of resting dynamic ventilatory parameters like FEV1 as predictive factor of maximal exercise capacity, indicating the need to include lung function test in the routine evaluation of these elite soccer players.
In our study population measured is in agreement with average values reported in previous studies on international elite soccer players . Currently, ventilatory parameters are poorly investigated in these athletes and our study provides a contribution to this field.
Significant relationships between , FEV1 and BRR with exercise capacity were found, and these parameters were significantly different between high and low performers. This finding offers a new insight in this field as it provides a quantitative relationship between dynamic ventilatory parameters and exercise performance in these athletes. Indeed, described by a linear combination of most representative ventilatory parameters resulting by PCA analysis (i.e., , FEV1 and BRR) corresponds to a physiological relationship among these variables. Therefore, the physiological parameters obtained from a simple test like spirometry might be used to distinguish high from low performers.
The significant correlation found between (V d /V t )peak and on the whole study population represents an interesting finding: high values of found in Hi-M elite soccer athletes correspond to low values of (V d /V t )peak. High ventilatory efficiency is more important in achieving better performance than VE. These parameters depend also on anthropometric and demographic characteristics which in our study were not significantly different between groups of players, independently of exercise capacity, ventilatory capacity, or role played.
Our study did not find any significant difference in physiological parameters among different playing roles. This may be related to lack of specific training among the roles.
Limitation of the study
The subjective reason to stop exercise, either dyspnoea, muscular fatigue or both, was not recorded, despite most subjects informally reported muscular fatigue. This is not surprising given the lack of ventilatory limitation to exercise shown by these athletes.
Respiratory parameters can play a determinant role in qualitative and quantitative evaluation of professional soccer players performance as it was confirmed by results discussed so far. Further studies should be addressed to investigate the impact of different training programmes on dynamic ventilatory parameters.
Bradley PS, Sheldon W, Wooster B, Olsen P, Boanas P, Krustrup P: High-intensity running in English FA premier league soccer matches. J Sports Sci. 2009, 27: 159-168.
Bradley PS, Di Mascio M, Peart D, Sheldon B: High-intensity activity profiles of elite soccer players at different performance levels. J Strength Cond Res. 2010, 24: 2343-2351.
Di Salvo V, Pigozzi F, González-Haro C, Laughlin MS, De Witt JK: Match performance comparison in Top English soccer leagues. Int J Sports Med. 2013, 34: 526-532.
Gregson W, Drust B, Atkinson G, Salvo VD: Match-to-match variability of high-speed activities in premier league soccer. Int J Sports Med. 2010, 31: 237-242.
Osgnach C, Poser S, Bernardini R, Rinaldo R, di Prampero PE: Energy cost and metabolic power in elite soccer: a new match analysis approach. Med Sci Sports Exerc. 2010, 42: 170-178.
Hoff J: Training and testing physical capacities for elite soccer players. J Sports Sci. 2005, 23: 573-582.
Kemi OJ, Hoff J, Engen LC, Helgerud J, Wisløff U: Soccer specific testing of maximal oxygen uptake. J Sports Med Phys Fitness. 2003, 43: 139-144.
Oberacker LM, Davis SE, Haff GG, Witmer CA, Moir GL: The Yo-Yo IR2 test: physiological response, reliability, and application to elite soccer. J Strength Cond Res. 2012, 26: 2734-2740.
Svensson M, Drust B: Testing soccer players. J Sports Sci. 2005, 23: 601-618.
Lago-Peñas C, Casais L, Dellal A, Rey E, Dominguez E: Anthropometric and physiological characteristics of young soccer players according to their playing positions: relevance for competition success. J Strength Cond Res. 2011, 25: 3358-3367.
Sporis G, Jukic I, Ostojic SM, Milanovic D: Fitness profiling in soccer: physical and physiologic characteristics of elite players. J Strength Cond Res. 2009, 23: 1947-1953.
Da Silva JF, Guglielmo LG, Bishop D: Relationship between different measures of aerobic fitness and repeated-sprint ability in elite soccer players. J Strength Cond Res. 2010, 24: 2115-2121.
Nassis GP, Geladas ND, Soldatos Y, Sotiropoulos A, Bekris V, Souglis A: Relationship between the 20-m multistage shuttle run test and 2 soccer-specific field tests for the assessment of aerobic fitness in adult semi-professional soccer players. J Strength Cond Res. 2010, 24: 2693-2697.
Castagna C, Manzi V, Impellizzeri F, Weston M, Barbero Alvarez JC: Relationship between endurance field tests and match performance in young soccer players. J Strength Cond Res. 2010, 24: 3227-3233.
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-528.
Bangsbo J: The physiology of soccer-with special reference to intense intermittent exercise. Acta Physiol Scand. 1994, 15 (suppl 619): 1-156.
Bangsbo J, Mohr M, Krustrup P: Physical and metabolic demands of training and match-play in the elite football player. J Sports Sci. 2006, 24: 665-674.
Bangsbo J, Iaia FM, Krustrup P: Metabolic response and fatigue in soccer. Int J Sports Physiol Perform. 2007, 2: 111-127.
Gatterer H, Faulhaber M, Patterson C: Real time VO2 measurements during soccer match-play. J Sports Med Phys Fitness. 2010, 50: 109-110.
Stølen T, Chamari K, Castagna C, Wisløff U: Physiology of soccer: an update. Sports Med. 2005, 35: 501-536.
Ekstrand J, Waldén M, Hägglund M: A congested football calendar and the wellbeing of players: correlation between match exposure of European footballers before the World Cup 2002 and their injuries and performances during that World Cup. Br J Sports Med. 2004, 38: 93-97.
Hoff J, Wisløff U, Engen LC, Kemi OJ, Helgerud J: Soccer specific aerobic endurance training. Br J Sports Med. 2002, 36: 218-221.
Baldari C, Videira M, Madeira F, Sergio J, Guidetti L: Lactate removal during active recovery related to the individual anaerobic and ventilatory thresholds in soccer players. Eur J Appl Physiol. 2004, 93: 224-230.
Boone J, Vaeyens R, Steyaert A, VandenBossche L, Bourgois J: Physical fitness of elite Belgian soccer players by player position. J Strength Cond Res. 2012, 26: 2051-2057.
Casajús JA: Seasonal variation in fitness variables in professional soccer players. J Sports Med Phys Fitness. 2001, 41: 463-469.
Kalapotharakos VI, Ziogas G, Tokmakidis SP: Seasonal aerobic performance variations in elite soccer players. J Strength Cond Res. 2011, 25: 1502-1507.
Marques-Neto SR, Maior AS, MaranhãoNeto GA, Santos EL: Analysis of heart rate deflection points to predict the anaerobic threshold by a computerized method. J Strength Cond Res. 2012, 26: 1967-1974.
Metaxas TI, Koutlianos NA, Kouidi EJ, Deligiannis AP: Comparative study of field and laboratory tests for the evaluation of aerobic capacity in soccer players. J Strength Cond Res. 2005, 19: 79-84.
McMillan K, Helgerud J, Grant SJ, Newell J, Wilson J, Macdonald R, Hoff J: Lactate threshold responses to a season of professional British youth soccer. Br J Sports Med. 2005, 39: 432-436.
Tønnessen E, Hem E, Leirstein S, Haugen T, Seiler S: Maximal aerobic power characteristics of male professional soccer players, 1989–2012. Int J Sports Physiol Perform. 2013, 8: 323-329.
Ziogas GG, Patras KN, Stergiou N, Georgoulis AD: Velocity at lactate threshold and running economy must also be considered along with maximal oxygen uptake when testing elite soccer players during preseason. J Strength Cond Res. 2011, 25: 414-419.
ATS/ERS Task Force Standardisation of Lung Function Testing: General consideration for lung function testing. Eur Respir J. 2005, 26: 153-161.
Campbell SC: A comparison of the maximum volume ventilation with forced expiratory volume in one second: an assessment of subject cooperation. J Occup Med. 1982, 24: 531-533.
Kuipers H, Verstappen FT, Keizer HA, Geurten P, van Kranenburg G: Variability of aerobic performance in the laboratory and its physiologic correlates. Int J Sports Med. 1985, 6: 197-201.
Wasserman K, Hansen JE, Sue DY, Stringer WW, Whipp BJ: Principles of Exercise Testing and Interpretation, Volume chapter 4. 2005, Philadelphia: Lippincott Williams and Wilkins, 87-90. Four
Tanaka H, Monahan KD, Seals DR: Age-predicted maximal heart rate revisited. J Am Coll Cardiol. 2001, 37: 153-156.
Jolliffe IT: Introduction. Principal Component Analysis. Edited by: Bickel P, Diggle P, Fienberg S, Krickeberg K, Olkin I, Wermuth N, Zege S. 2002, New York: Springer, 1-9. 2
Special thanks to ACF Fiorentina, Genoa CFC, Lecce Calcio, Napoli SSC, AC Milan, US Città di Palermo soccer teams.
The authors declare that they have no competing interests.
ADP conceived the study, participated in its design, collected data and helped to draft the manuscript. SM participated in the design of the study, performed the statistical analysis and helped to draft the manuscript. GAC and MLM participated in the design of the study. GV participated in the design of the study and performed the statistical analysis. NA coordinated the study and helped to draft the manuscript. All authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.