Skip to main content
Sports Medicine - Open
Download PDF
  • Original Research Article
  • Open access
  • Published:

Reference Values of Cardiopulmonary Exercise Test Parameters in the Contemporary Paediatric Population

Sports Medicine - Open volume 9, Article number: 68 (2023) Cite this article

Abstract

Background

The evaluation of health status by cardiopulmonary exercise test (CPET) has shown increasing interest in the paediatric population. Our group recently established reference Z-score values for paediatric cycle ergometer VO2max, applicable to normal and extreme weights, from a cohort of 1141 healthy children. There are currently no validated reference values for the other CPET parameters in the paediatric population. This study aimed to establish, from the same cohort, reference Z-score values for the main paediatric cycle ergometer CPET parameters, apart from VO2max.

Results

In this cross-sectional study, 909 healthy children aged 5–18 years old underwent a CPET. Linear, quadratic, and polynomial mathematical regression equations were applied to identify the best CPET parameters Z-scores, according to anthropometric parameters (sex, age, height, weight, and BMI). This study provided Z-scores for maximal CPET parameters (heart rate, respiratory exchange ratio, workload, and oxygen pulse), submaximal CPET parameters (ventilatory anaerobic threshold, VE/VCO2 slope, and oxygen uptake efficiency slope), and maximum ventilatory CPET parameters (tidal volume, respiratory rate, breathing reserve, and ventilatory equivalent for CO2 and O2).

Conclusions

This study defined paediatric reference Z-score values for the main cycle ergometer CPET parameters, in addition to the existing reference values for VO2max, applicable to children of normal and extreme weights. Providing Z-scores for CPET parameters in the paediatric population should be useful in the follow-up of children with various chronic diseases. Thus, new paediatric research fields are opening up, such as prognostic studies and clinical trials using cardiopulmonary fitness outcomes.

Trial registration NCT04876209—Registered 6 May 2021—Retrospectively registered, https://clinicaltrials.gov/ct2/show/NCT04876209.

Key Points

  • The dissemination of CPET in paediatrics requires valid paediatric reference values to define the upper and lower normal limits for the main parameters, in addition to the existing reference values for VO2max.

  • The paediatric cycle ergometer CPET reference values were generated using the best mathematical model for the Z-score construction, according to the main anthropometric predictors (sex, age, height, weight, and BMI).

  • Providing Z-scores for the main cycle ergometer CPET parameters in the paediatric population will be useful in the follow-up of children with chronic diseases.

Background

There has been increasing interest in evaluating the health status of children by cardiopulmonary exercise test (CPET) [1]. In many paediatric chronic diseases, impaired physical capacity assessed by CPET correlates with lower health-related quality of life scores and may indicate the early onset of physical deconditioning, such as in congenital heart disease (CHD) [2], cancer [3], asthma [4], or kidney disease [5]. Therefore, CPET stands as a key examination to evaluate cardiopulmonary fitness in healthy and chronically ill children [6], as well as to promote physical activity and cardiovascular rehabilitation from a young age [7].

The VO2max is the main CPET parameter to evaluate the level of physical capacity, and an independent predictor of cardiovascular risk [8]. The VO2max reference values defined by Cooper et al. in 1984, are based on linear equations from a small cohort of 109 healthy children with normal weight and are currently less adapted to the contemporary paediatric population [9, 10]. Our group recently established reference Z-score values for paediatric cycle ergometer VO2max, applicable to normal and extreme weights, from a cohort of 1141 healthy children, including 909 children in the development cohort and 232 children in the validation cohort [11]. For both sexes, the Z-score equations were modelled with a logarithmic function of VO2max, height, and BMI.

While reference values have been established for VO2max, there are currently no validated reference values for the other main CPET parameters. As recently stated by Takken et al., “there is still a lot of progress to be made” in validating reference equations in paediatric CPET [12]. Yet, the other CPET parameters are complementary to VO2max to assess cardiac, muscular and respiratory functions. In particular, greater attention has been given to submaximal parameters in paediatric CPET: the ventilatory anaerobic threshold (VAT) reflects the level of muscular deconditioning and has been used to determine the intensity of physical activity in paediatric cardiovascular rehabilitation programs [7]; the ventilatory efficiency slope (e.g. VE/VCO2 slope) increases in paediatric heart failure, especially in the most complex heart diseases [13]; and the oxygen uptake efficiency slope (OUES) may be used as a surrogate of VO2max in children unable to perform a maximal exercise test [14].

Therefore, the dissemination of CPET in children requires contemporary paediatric reference values, based on valid mathematical models to define the upper and lower normal limits for each CPET parameter.

In this study, we aimed to establish paediatric reference Z-score values for the main cycle ergometer CPET parameters, apart from VO2max, from a large cohort of healthy children representative of the contemporary paediatric population, including extreme weights.

Methods

Study Design and Population

To elaborate paediatric CPET parameters Z-scores, we used a French CPET database, initially generated by pooling all subjects aged less than 18 years enrolled in previous prospective controlled studies carried out in two paediatric CPET laboratories (centre 1: M3C Regional Paediatric and Congenital Cardiology Centre, Montpellier University Hospital, France; centre 2: Paediatric Cardiology and Rehabilitation Centre, Saint-Pierre Institute, Palavas-Les-Flots, France) [2, 7, 13,14,15,16,17,18]. We identified all subjects who underwent a complete CPET with a high-quality score (≥ 10 points) according to the ATS/ACCP statement [19], over a period of 10 years (from November 2010 to March 2020).

Children referred for CPET with the following clinical criteria were selected: non-severe functional symptoms at rest (murmur, palpitation) or during exercise (chest pain or dyspnoea), and completely normal cardiological check-up, including physical examination, electrocardiogram, and echocardiography. Children with overweight or obesity who were referred to the CPET laboratory for physical fitness check-up were also eligible for the study. However, children with any other chronic disease, medical condition (cardiac, neurological, respiratory, muscular, or renal), or medical treatment and those requiring any further specialized medical consultation were not eligible.

Anthropometric parameters were collected (sex, age, weight, and height) and body mass index (BMI) percentiles were used to define patient groups: underweight (BMI < 5th percentile), normal weight (BMI between 5 and 84th), and overweight or obesity (BMI ≥ 85th) [20].

CPET Procedures

CPET procedures for children enrolled in the study were harmonized and similar in both laboratories [2, 13, 14, 16]. We used a single CPET paediatric cycle ergometer protocol adapted to children to obtain a homogeneous incremental overall duration between 8 to 12 min including (a) 1-min baseline; (b) 3-min warm-up (10 or 20 watts), (c) fixed increments of 10, 15, or 20 watts each minute, (d) pedalling rate of 60–80 revolutions per minute, (e) 3-min active recovery (20 watts); f) 2-min rest [21]. The CPET was considered maximal when the child was unable to maintain a pedalling rate above 60 despite verbal encouragement. When the VO2max did not reach a plateau, the peak VO2 value was collected, as usual in paediatrics [22, 23].

Spirometry (flow volume curve) was performed at rest to measure the forced expiratory volume in 1 s (FEV1), the forced vital capacity (FVC), and the FEV1/FVC ratio (FEV1/FVC%) [24].

From 2010 to 2015, the same technical devices were used in both CPET laboratories: paediatric face masks (Hans Rudolph, Shawnee, KS, USA), calibrated gas analyser (Oxycon Pro, Jaeger, Erich Jaeger GmbH, Hoechberg, Germany), breath-to-breath measurement software (Windows 98, Jaeger), 12-lead ECG equipment (CardioSoft, GE Healthcare, Little Chalfont, UK), pulse oximeter (Nellcor, Medtronic, Fridley, MN, USA), and manual sphygmomanometer with adapted paediatric cuffs. From 2015, the centre 1 used the following technical devices: calibrated gas analyser (Quark CPET, Cosmed Srl, Pavonna di Albano, Italy), breath-by-breath measurement software (Windows 7–10, Omnia, Cosmed), 12-lead ECG equipment (Norav, Germany), and pulse oximeter (Nonin Medical Inc, Plymouth, MN 55441 USA).

CPET Parameters

Apart from VO2max, the CPET parameters were grouped into 3 categories: maximal parameters, submaximal parameters, and ventilatory parameters.

The maximal CPET parameters included: (1) the maximal heart rate; (2) the maximal respiratory exchange ratio (RERmax), which corresponds to the ratio of CO2 elimination (VCO2) and oxygen uptake (VO2); (3) the maximum workload in Watts (maintained for at least 30 s by the patient); (4) and the maximum oxygen pulse (O2 pulsemax), e.g. a surrogate of stroke volume at peak exercise, which corresponds to the ratio between VO2max and maximal heart rate [19].

The submaximal CPET parameters included: (1) the VAT, e.g. the point at which minute ventilation increases disproportionally relative to VO2, which reflects muscular response to exercise [2, 25], and manually calculated by a single investigator using V-slope method [26] and expressed as a percentage of the predicted VO2max (%-predicted VAT) [11]; (2) the VE/VCO2 slope, e.g. an indicator of ventilation-perfusion ratio during exercise [13, 27], calculated from breath-by-breath data using linear regression of minute ventilation (VE) and CO2 production (VCO2) obtained during incremental exercise (VE = [VE/VCO2 slope] × VCO2 + b) [13], and measured from the beginning of incremental exercise (after the warm-up period) to maximum exercise (or respiratory compensation point when present) [2, 13]; (3) and the OUES, e.g. a submaximal surrogate parameter of VO2max, calculated from breath-by-breath data using linear regression of VO2 on logarithmically converted VE (VO2 = OUES × log10 VE + b), measured from the beginning of the exercise test to its maximal point, and expressed as a weight-normalized value (OUESKg) [14].

The ventilatory parameters included: (1) the maximum tidal volume (VTmax), which corresponds to the patient's inspiratory volume for each breath at maximum effort averaged over 30 s [19]; (2) the maximum respiratory rate (RRmax), which corresponds to the respiratory rate averaged over 30 s at maximum effort [19]; (3) the breathing reserve (BR), calculated using the formula BR (%) = [MVV − VEmax] ÷ MVV*100 (with MVV, maximal voluntary ventilation = 35*FEV1) [19]; (4) the maximum ventilatory equivalent for CO2 (VEqCO2max) which corresponds to the ratio between VE and VCO2 averaged over 30 s at maximum effort [28]; and (5) the maximum ventilatory equivalent for O2 (VEqO2max) which corresponds to the ratio between the VE and the VO2 averaged over 30 s at maximum effort [19].

Statistics

The study population was described using means and standard deviations (SD) for quantitative variables and with frequencies for qualitative variables. Quantitative variables were compared using Student's t-test when the distribution was Gaussian and using the Mann–Whitney test otherwise. For qualitative variables, groups were compared using the chi-squared test or Fisher’s exact test.

To identify the best model for the Z-score construction, the values of each CPET parameter were modelized according to the main anthropometric predictors: sex, age, height, weight, and BMI. Three different mathematical regression models were successively applied: (1) a linear model, (2) a quadratic model, and (3) a polynomial model of degree 2 [29]. Models with several potential predictors were also evaluated. The models maximizing Pearson's coefficient of determination (R2, e.g. the proportion of variance of each CPET parameter explained by the model) were identified. To generate a model respecting the condition of homoscedasticity required for the Z-score construction, we also evaluated all linear, quadratic, and polynomial models with log–log relationships (e.g., modelling the natural logarithm of the CPET parameters with the natural logarithm of the anthropometric predictor(s)). Homoscedasticity was assessed by comparing the standard deviation of the estimated model residuals between quintiles of the predicted value. Finally, all the models selected met this condition. The data were first analysed for all sexes. A sex effect and an interaction between sex and anthropometric predictor(s) were then tested in the selected model. When one of these effects was significant, separate models were performed for boys and girls. For each gender, we finally selected the model with the highest R2 and the lowest heteroscedasticity in the residuals [30]. When the complexity of the model (quadratic, polynomial or multivariate) only allowed a gain of less than 0.01 point of R2, the simplest model was chosen. When no anthropometric predictor explained more than 5% of the variability of the parameter of interest (e.g. R2 < 0.05 for all models), a simple Z-score, obtained by subtracting the mean and dividing by the standard deviation, was constructed after checking graphically that the distribution was close to the Gaussian distribution. The final Z-score distribution was described using median, 5th, and 95th percentiles. Under the condition of standard normal distribution, the expected values of the Z-score are 0, -1.64, and 1.64, respectively.

Statistical analyses were performed at the conventional two-tailed α level of 0.05 using SAS version 9.4 (SAS Institute, Cary, NC).

Results

Population

During the 10-year study period, a total of 909 children aged 5–18 (mean age at 11.4 ± 2.7 years, 477 boys and 432 girls), following the eligibility criteria, were selected from the CPET database, and included in the study. The numbers of patients per age group were N = 142 (i.e. 60 girls and 82 boys) [5–8 years], N = 333 (i.e. 156 girls and 177 boys) [9–11 years], N = 317 (i.e. 153 girls and 164 boys) [12–14 years], and N = 117 (i.e. 63 girls and 54 boys) [15–18 years]. After stratification on BMI, we identified 38 underweight children, 639 children with normal weight, and 232 children with overweight or obesity (of which 81 children with overweight ≥ 85th percentile and 151 children with obesity ≥ 95th percentile). No significant differences between boys and girls were observed in terms of age, height, weight, and BMI. The main anthropometric data of the population are summarized in Table 1.

Table 1 Main anthropometric and CPET data
Full size table

CPET Parameters

Overall, this study fulfilled 11 of the 14 criteria from the ATS/ACCP statement on cardiopulmonary exercise testing (Additional file 1: Table S1).

Significant sex differences were observed in most CPET parameters, except for VE/VCO2 slope, VEqCO2max, maximum heart rate, and %-predicted VAT (Table 1). VAT, OUES, O2 pulsemax, and VTmax were mainly influenced by height and weight, the maximum workload was influenced by height and age, and the VE/VCO2 slope was slightly influenced by age. Maximum heart rate, breathing reserve, and VEqCO2max were not influenced by any anthropometric parameters. The CPET parameters with the highest coefficients of determination (R2) according to anthropometric data were reported in Additional file 1: Table S2.

CPET Parameters Z-Scores

Table 2 summarizes the equations for each CPET parameter Z-score and reported the variability of each parameter logarithm (R2), including median, 5th, and 95th Z-score percentiles. Figures 1, 2, and 3 summarize the correlation between observed and predicted values, distinguishing the different BMI groups (underweight, normal weight, and overweight/obesity). The median CPET values with their ranges [5th percentile; 95th percentile] in the normal weight population were presented in Additional file 1: Table S3.

Table 2 Z-score equations of CPET parameters
Full size table
Fig. 1
figure 1

Correlation between observed and predicted values of CPET maximal parameters using the Z-score model. The “underweight” group was represented by blue points, the “normal weight” group by green points, and the “overweight/obesity” group by red points. The correlation between measured and predicted values using the Z-score model for RERmax (panel A), workloadmax (panel B) and O2 pulsemax (panel C)

Full size image
Fig. 2
figure 2

Correlation between observed and predicted values of CPET submaximal parameters using the Z-score model. The “underweight” group was represented by blue points, the “normal weight” group by green points, and the “overweight/obesity” group by red points. The correlation between measured and predicted values using the Z-score model for VAT (panel A), VE/VCO2 slope (panel B) and OUES (panel C)

Full size image
Fig. 3
figure 3

Correlation between observed and predicted values CPET ventilatory parameters using the Z-score model. The “underweight” group was represented by blue points, the “normal weight” group by green points, and the “overweight/obesity” group by red points. The correlation between measured and predicted values using the Z-score model for VTmax (panel A) and RRmax (panel B)

Full size image

Z-Scores of Maximal CPET Parameters

Maximum Heart Rate (HRmax)

A single linear equation valid for both sexes was generated to define HRmax Z-scores:

$${\text{HR}}_{\max } Z{\text{ - score}} = \frac{{\left( {{\text{HR}}_{\max } \left( {{\text{bpm}}} \right) - 188.1} \right)}}{8.9}$$

In the overall cohort, a HRmax value of 173 bpm corresponded to the 5th percentile.

Maximum Respiratory Exchange Ratio (RERmax)

For boys, the mathematical model using natural logarithms of RERmax, height, and weight was the best fit for the data. This equation determined 16% of the variability of the RERmax logarithm (R2 = 0.16). For girls, the mathematical model using RERmax without transformation and age was the best fit for the data. This equation determined 14% of the variability of the RERmax (R2 = 0.14). The correlation between predicted and observed RERmax using the Z-score model is illustrated in Fig. 1 (panel A). In children with a normal weight, RERmax values of 1 in girls and 1.01 in boys corresponded to the 5th percentile (Additional file 1: Table S3).

Maximum Workload (Workloadmax)

For both sexes, the mathematical model using natural logarithms of workloadmax, height, and age was the best fit for the data. These equations determined 60% of the variability of the workloadmax logarithm (R2 = 0.60) in girls and 76% (R2 = 0.76) in boys. The correlation between predicted and observed workloadmax using the Z-score model is illustrated in Fig. 1 (panel B).

Maximum Oxygen Pulse (O2 Pulsemax)

For both sexes, the mathematical model using natural logarithms of O2 pulsemax, weight, and height was the best fit for the data. These equations determined 69% of the variability of the O2 pulsemax logarithm (R2 = 0.69) in girls and 77% (R2 = 0.77) in boys. The correlation between predicted and observed O2 pulsemax using the Z-score model is illustrated in Fig. 1 (panel C).

Z-Scores of Submaximal CPET Parameters

Ventilatory Anaerobic Threshold (VAT)

For both sexes, the mathematical model using natural logarithms of VAT, height, and weight was the best fit for the data. These equations determined 54% of the variability of the VAT logarithm (R2 = 0.54) in girls and 68% (R2 = 0.68) in boys. The correlation between predicted and observed VAT when using the Z-score model is illustrated in Fig. 2 (panel A). In children with a normal weight, a %-predicted VAT value of 55% corresponded to the 5th percentile Additional file 1: Table S3).

VE/VCO2 Slope

For both sexes, the mathematical model using natural logarithms of VE/VCO2 slope and age was the best fit for the data. These equations determined 6% of the variability of the VE/VCO2 slope logarithm (R2 = 0.06) for the girl and 19% (R2 = 0.19) for the boys. The correlation between predicted and observed VE/VCO2 slopes using the Z-score model is illustrated in Fig. 2 (panel B). In children with a normal weight, a VE/VCO2 slope value of 37.6 corresponded to the 95th percentile, overall (Additional file 1: Table S3), with a decrease with age from a VE/VCO2 slope value of 40.7 in young children aged 5 to 8 years, to a VE/VCO2 slope value of 35.9 in adolescents aged 15 to 18 years.

OUES

For both sexes, the mathematical model using natural logarithms of OUES, height, and weight was the best fit for the data. These equations determined 57% of the variability of the OUES logarithm (R2 = 0.57) for girls and 67% (R2 = 0.67) for boys. The correlation between predicted and observed OUES using the Z-score model is illustrated in Fig. 2 (panel C). In children with a normal weight, OUESkg values of 28.9 in girls and 38.9 in boys corresponded to the 5th percentile (Additional file 1: Table S3).

Z-Scores of Ventilatory CPET Parameters

Maximum Tidal Volume (VTmax)

For both sexes, the mathematical model using natural logarithms of VTmax, height, and weight was the best fit for the data. These equations determined 72% of the variability of the VTmax logarithm (R2 = 0.72) for girls and 78% (R2 = 0.78) for boys. The correlation between predicted and observed VTmax using the Z-score model is illustrated in Fig. 3 (panel A).

Maximum Respiratory Rate (RRmax)

For girls, the mathematical model using RRmax and height was the best fit for the data. This equation determined 5% of the variability of the RRmax logarithm (R2 = 0.05). For boys, the mathematical model using natural logarithms of RRmax and weight was the best fit for the data. This equation determined 7% of the variability (R2 = 0.07). The correlation between predicted and observed RRmax using the Z-score model is illustrated in Fig. 3 (panel B). In children with a normal weight, RRmax values of 67/min for girls and 69/min for boys corresponded to the 95th percentile.

Breathing Reserve

This parameter was not influenced by anthropometric variables. The median value of breathing reserve was 30% [0%; 52%] for girls and 25% [0%; 45%] for boys.

Maximum Ventilatory Equivalent for CO2 (VEqCO2max)

No significant sex differences were found for VEqCO2max and this parameter was not influenced by anthropometric variables. The median value of VEqCO2max was 30.6 [25.1; 37.3].

Maximum Ventilatory Equivalent for O2 (VEqO2max)

This parameter was not influenced by anthropometric variables, except for sex. The median value of VEqO2max was 35.8 [27.9; 46.0] for girls and 34.3 [27.2; 45.0] for boys.

Discussion

From a cohort of 909 healthy children aged 5–18 years who underwent a cycle ergometer CPET, this study defined paediatric reference Z-score values of the main paediatric CPET parameters, apart from VO2max, whose reference values were recently reported by our group [11]. These paediatric CPET reference values were generated using the best mathematical model (linear, quadratic, or polynomial) for the Z-score construction, according to the main anthropometric predictors (sex, age, height, weight, and BMI), to apply to normal and extreme weights. This paediatric cohort is representative of the general contemporary paediatric population, with a balanced sex ratio and a 25%-proportion of overweight or obese children [31], and had recently undergone external validation from a CPET cohort of 232 German and American healthy children [11]. This study intended to satisfy a high methodological quality level and fulfilled 11 of the 14 criteria from the ATS/ACCP statement on cardiopulmonary exercise testing [19]. Unsurprisingly, sex differences were observed in most CPET parameters (apart from maximum heart rate and VEqCO2max), resulting in Z-score models generated for boys and girls separately.

In terms of maximal CPET parameters, the existence of Z-scores for the main maximum parameters will improve discussions on maximality criteria in paediatric CPET. The historic equations on predicted maximum heart rate ([220–age]; or 0.65 × [210–age]) were adapted from adult CPET studies and the commonly used value of 80% of predicted maximum heart rate to define maximum exercise was set arbitrarily [19, 32]. Our data confirmed that maximum heart rate was not significantly influenced by age in the paediatric population. Similarly, different cut-off values of maximum RER have been used to define maximal exercise in paediatrics (> 1.05 or > 1.1) [19, 33], but the youngest healthy children often do not reach an RER of 1.05. In our cohort, 10% of subjects had a maximum RER between 1 and 1.05, and 15% between 1.06 and 1.1, most of which (78%) were aged < 12 years. A similar value of RER around 1 at maximal exercise in children has been previously reported [23]. These data confirm that it is probably not appropriate to use RER in maximality criteria for paediatric exercise testing. This study also provides reference values for maximum oxygen pulse, a determinant of aerobic capacity which can be impaired in children with cardiac disease, as in Fontan circulation with limited preload [34]. The maximum oxygen pulse equation will help the clinician to better define the lower normal limit for this important cardiac parameter.

Moreover, the existence of paediatric reference values for submaximal CPET parameters will be useful to interpret cardiopulmonary fitness in children with serious chronic diseases. We found that the reference values of VAT were determined by natural logarithms of height and weight. The VAT is an indicator of aerobic fitness, useful for exercise prescription, especially in cardiac rehabilitation programs to monitor the effect of physical training [7, 35, 36]. In adult studies, wide range of normal values for VAT from 35 to 80% of the predicted VO2max have been reported [19]. Furthermore, percent-predicted VAT values from 50 to 60% observed in adult sedentary subjects are commonly used to define physical deconditioning and patient eligibility for rehabilitation programs. Interestingly, in our study, the 5th percentile of the VAT expressed as a percentage of the predicted VO2max was at 55%, which is exactly the cut-off value used in the QUALIREHAB randomized controlled trial to define children with CHD eligible for cardiac rehabilitation [7]. Furthermore, in this study, the reference values of VE/VCO2 slope were determined by a natural logarithm of age, with normal mean values of 33 in youngest children and 28 in adolescents. These results are consistent with the VE/VCO2 cut-off value < 28 previously reported to define normal ventilatory efficiency in paediatric CPET [33]. The VE/VCO2 slope increases in pulmonary blood flow maldistribution and its prognostic value has been demonstrated in adult heart failure [37]. The existence of paediatric reference Z-scores for VE/VCO2 slope opens up new research perspectives to define the prognostic value of this submaximal parameter in various paediatric chronic diseases involving ventilation/perfusion mismatch. As for the OUES, it was strongly influenced by the natural logarithms of height and weight. This submaximal parameter is classically strongly correlated with VO2max and may be useful in severe conditions compromising the achievement of a maximal CPET and making it difficult to accurately interpret peak VO2 [14]. The availability of paediatric reference Z-score values for OUES will therefore facilitate the dissemination of this CPET parameter in routine clinical practice. For example, if we focus on children with normal weight, a submaximal CPET with OUESKg values < 28.9 in girls or < 38.9 in boys, e.g. corresponding to the 5th percentile, may be suggestive of impaired aerobic fitness.

Finally, the existence of paediatric reference Z-score values for ventilatory exercise parameters will improve the interpretation of paediatric CPET, as exercise ventilation may be affected in many paediatric chronic diseases [3, 4, 38, 39]. Ventilatory exercise parameters interact with each other and should not be interpreted independently [28]. For example, in paediatric asthma, e.g. the most common paediatric chronic disease worldwide, physical aerobic fitness is impaired in a quarter of children, with patterns of physical deconditioning in a third of them [4]. Previous studies have suggested that the hyperventilation syndrome during exercise, well described in adult patients [40], was also present in paediatrics and reflected dysfunctional breathing and poor asthma control [41]. Moreover, our group recently showed that a lower breathing reserve was a predictor of limited physical capacity in children with asthma [4]. Moreover, children with expiratory flow limitation, as in cystic fibrosis, present an abnormal breathing pattern during exercise with rapid shallow breathing, evidenced by a higher RRmax and a lower VTmax [39]. Nevertheless, overall, exercise ventilatory parameters have been scarcely analysed in children [42, 43]. In addition to the existing Z-scores for paediatric spirometry [44], the Z-scores for ventilatory CPET parameters represent a major advance in the analysis of exercise physiology in children.

Study Limitation

Reference values were not analysed in function of other parameters influencing physical capacity, such as lean mass, pubertal status, or the level of physical activity, as in the study on paediatric VO2max reference Z-scores [11]. The impact of ethnic variation on body fat and muscle mass has not been investigated because of restrictive regulations on the assessment of ethnicity in paediatric clinical research. Moreover, the children included in this study were not community-based, as we were concerned that if we had drawn lots of schools to identify volunteers to perform a CPET, a large proportion of child volunteers would have been athletes. As previously stated by Cumming, “clinic patients without heart defects probably serve as a better normal control group than children obtained from the school system” [45].

Conclusions

This study defined paediatric Z-score reference values for the main cycle ergometer CPET parameters, applicable to children of normal and extreme weights. This study intended to satisfy a high methodological quality level, by enrolling a large number of subjects, covering wide age and weight ranges, and fulfilling most criteria on high-quality CPET assessment. In addition to the existing reference values for VO2max, providing Z-scores reference values of maximal, submaximal, and ventilatory CPET parameters in the paediatric population should be useful in the follow-up of children with various chronic diseases. Thus, new paediatric research fields are opening up, such as prognostic studies and clinical trials using cardiopulmonary fitness parameters as primary or secondary outcomes.

Availability of Data and Materials

The data underlying this article will be shared on reasonable request to the corresponding author.

Abbreviations

ATS:

American Thoracic Society

ACCP:

American College of Chest Physicians

BMI:

Body mass index

BR:

Breathing reserve

CPET:

Cardiopulmonary exercise test

VCO2 :

Carbon dioxide production

FEV1:

Forced expiratory volume in 1 s

FVC:

Forced vital capacity

FEV1/FVC:

Tiffeneau index

HRmax :

Maximum heart rate

O2 pulsemax :

Maximum oxygen pulse

OUES:

Oxygen uptake efficiency slope

OUESkg :

Weight-normalized value of oxygen uptake efficiency slope

RERmax :

Maximum respiratory exchange ratio

RRmax :

Maximum respiratory rate

SD:

Standard deviation

VAT:

Ventilatory anaerobic threshold

%-Predicted VAT:

VAT expressed as a percentage of the predicted VO2max

VE:

Minute ventilation

VEqCO2max :

Maximum ventilatory equivalent for CO2

VEqO2max :

Maximum ventilatory equivalent for O2

VO2max :

Maximum oxygen uptake

VTmax :

Maximum tidal volume

References

  1. Takken T, Bongers BC, van Brussel M, Haapala EA, Hulzebos EHJ. Cardiopulmonary exercise testing in pediatrics. Ann Am Thorac Soc. 2017;14:S123–8.

    Article  PubMed  Google Scholar 

  2. Amedro P, Gavotto A, Guillaumont S, Bertet H, Vincenti M, De La Villeon G, et al. Cardiopulmonary fitness in children with congenital heart diseases versus healthy children. Heart Br Card Soc. 2018;104:1026–36.

    Google Scholar 

  3. Gavotto A, Dubard V, Avesani M, Huguet H, Picot M-C, Abassi H, et al. Impaired aerobic capacity in adolescents and young adults after treatment for cancer or non-malignant haematological disease. Pediatr Res. 2023;

  4. Moreau J, Socchi F, Renoux MC, Requirand A, Abassi H, Guillaumont S, et al. Cardiopulmonary fitness in children with asthma versus healthy children. Arch Dis Child. 2022;archdischild-2021-323733.

  5. Pattaragarn A, Warady BA, Sabath RJ. Exercise capacity in pediatric patients with end-stage renal disease. Perit Dial Int J Int Soc Perit Dial. 2004;24:274–80.

    Article  Google Scholar 

  6. Verduci E, Bronsky J, Embleton N, Gerasimidis K, Indrio F, Köglmeier J, et al. Role of dietary factors, food habits and lifestyle in childhood obesity development. A position paper from the European Society for paediatric gastroenterology, hepatology and nutrition committee on nutrition. J Pediatr Gastroenterol Nutr. 2021

  7. Amedro P, Gavotto A, Legendre A, Lavastre K, Bredy C, De La Villeon G, et al. Impact of a centre and home-based cardiac rehabilitation program on the quality of life of teenagers and young adults with congenital heart disease: the QUALI-REHAB study rationale, design and methods. Int J Cardiol. 2019;283:112–8.

    Article  PubMed  Google Scholar 

  8. Kodama S, Saito K, Tanaka S, Maki M, Yachi Y, Asumi M, et al. Cardiorespiratory fitness as a quantitative predictor of all-cause mortality and cardiovascular events in healthy men and women: a meta-analysis. JAMA. 2009;301:2024–35.

    Article  CAS  PubMed  Google Scholar 

  9. Cooper DM, Weiler-Ravell D, Whipp BJ, Wasserman K. Aerobic parameters of exercise as a function of body size during growth in children. J Appl Physiol. 1984;56:628–34.

    Article  CAS  PubMed  Google Scholar 

  10. Cooper DM, Weiler-Ravell D. Gas exchange response to exercise in children. Am Rev Respir Dis. 1984;129:S47-48.

    Article  CAS  PubMed  Google Scholar 

  11. Gavotto A, Mura T, Rhodes J, Yin SM, Hager A, Hock J, et al. Reference values of aerobic fitness in the contemporary paediatric population. Eur J Prev Cardiol. 2023;zwad054.

  12. Takken T, Mylius CF, Paap D, Broeders W, Hulzebos HJ, Van Brussel M, et al. Reference values for cardiopulmonary exercise testing in healthy subjects—an updated systematic review. Expert Rev Cardiovasc Ther. 2019;17:413–26.

    Article  CAS  PubMed  Google Scholar 

  13. Gavotto A, Huguet H, Picot M-C, Guillaumont S, Matecki S, Amedro P. The V̇e/V̇co2 slope: a useful tool to evaluate the physiological status of children with congenital heart disease. J Appl Physiol Bethesda Md. 1985;2020(129):1102–10.

    Google Scholar 

  14. Gavotto A, Vandenberghe D, Abassi H, Huguet H, Macioce V, Picot M-C, et al. Oxygen uptake efficiency slope: a reliable surrogate parameter for exercise capacity in healthy and cardiac children? Arch Dis Child. 2020;105:1167–74.

    Article  PubMed  Google Scholar 

  15. Amedro P, Picot MC, Moniotte S, Dorka R, Bertet H, Guillaumont S, et al. Correlation between cardio-pulmonary exercise test variables and health-related quality of life among children with congenital heart diseases. Int J Cardiol. 2016;203:1052–60.

    Article  CAS  PubMed  Google Scholar 

  16. Abassi H, Gavotto A, Picot MC, Bertet H, Matecki S, Guillaumont S, et al. Impaired pulmonary function and its association with clinical outcomes, exercise capacity and quality of life in children with congenital heart disease. Int J Cardiol [Internet]. 2019 [cited 2019 Mar 11]; Available from: http://www.sciencedirect.com/science/article/pii/S0167527318369316

  17. Moreau J, Lavastre K, Romieu H, Charbonnier F, Guillaumont S, Bredy C, et al. Impact of Sophrology on cardiopulmonary fitness in teenagers and young adults with a congenital heart disease: the SOPHROCARE study rationale, design and methods. Int J Cardiol Heart Vasc. 2020;27:100489.

    PubMed  PubMed Central  Google Scholar 

  18. Amedro P, Werner O, Abassi H, Boisson A, Souilla L, Guillaumont S, et al. Health-related quality of life and physical activity in children with inherited cardiac arrhythmia or inherited cardiomyopathy: the prospective multicentre controlled QUALIMYORYTHM study rationale, design and methods. Health Qual Life Outcomes. 2021;19:187.

    Article  PubMed  PubMed Central  Google Scholar 

  19. American Thoracic Society, American College of Chest Physicians. ATS/ACCP Statement on cardiopulmonary exercise testing. Am J Respir Crit Care Med. 2003;167:211–77.

  20. Barlow SE, Expert Committee. Expert committee recommendations regarding the prevention, assessment, and treatment of child and adolescent overweight and obesity: summary report. Pediatrics. 2007;120(Suppl 4):S164-192.

    Article  PubMed  Google Scholar 

  21. Takken T, Blank AC, Hulzebos EH, van Brussel M, Groen WG, Helders PJ. Cardiopulmonary exercise testing in congenital heart disease: (contra)indications and interpretation. Neth Heart J Mon J Neth Soc Cardiol Neth Heart Found. 2009;17:385–92.

    CAS  Google Scholar 

  22. Rowland TW. Does peak VO2 reflect VO2max in children?: evidence from supramaximal testing. Med Sci Sports Exerc. 1993;25:689–93.

    Article  CAS  PubMed  Google Scholar 

  23. Armstrong N, Welsman J, Winsley R. Is peak VO2 a maximal index of children’s aerobic fitness? Int J Sports Med. 1996;17:356–9.

    Article  CAS  PubMed  Google Scholar 

  24. Brusasco V, Crapo R, Viegi G. American Thoracic Society, European Respiratory Society. Coming together: the ATS/ERS consensus on clinical pulmonary function testing. Eur Respir J. 2005;26:1–2.

    Article  CAS  PubMed  Google Scholar 

  25. Stephens P, Paridon SM. Exercise testing in pediatrics. Pediatr Clin North Am. 2004;51:1569–87.

    Article  PubMed  Google Scholar 

  26. Beaver WL, Wasserman K, Whipp BJ. A new method for detecting anaerobic threshold by gas exchange. J Appl Physiol Bethesda Md. 1985;1986(60):2020–7.

    Google Scholar 

  27. Weatherald J, Sattler C, Garcia G, Laveneziana P. Ventilatory response to exercise in cardiopulmonary disease: the role of chemosensitivity and dead space. Eur Respir J. 2018;51.

  28. Watson M, Ionescu MF, Sylvester K, Fuld J. Minute ventilation/carbon dioxide production in patients with dysfunctional breathing. Eur Respir Rev Off J Eur Respir Soc. 2021;30:200182.

    Article  Google Scholar 

  29. DeVore GR. Computing the Z score and centiles for cross-sectional analysis: a practical approach. J Ultrasound Med Off J Am Inst Ultrasound Med. 2017;36:459–73.

    Google Scholar 

  30. Mawad W, Drolet C, Dahdah N, Dallaire F. A review and critique of the statistical methods used to generate reference values in pediatric echocardiography. J Am Soc Echocardiogr Off Publ Am Soc Echocardiogr. 2013;26:29–37.

    Article  Google Scholar 

  31. Nittari G, Scuri S, Petrelli F, Pirillo I, di Luca NM, Grappasonni I. Fighting obesity in children from European World Health Organization member states. Epidemiological data, medical-social aspects, and prevention programs. Clin Ter. 2019;170:e223–30.

    CAS  PubMed  Google Scholar 

  32. Cooper KH, Purdy JG, White SR, Pollock ML, Linnerud AC. Age-fitness adjusted maximal heart rates. Role Exerc Int Med. 1977;10:78–88.

    Google Scholar 

  33. Rhodes J, Ubeda Tikkanen A, Jenkins KJ. Exercise testing and training in children with congenital heart disease. Circulation. 2010;122:1957–67.

    Article  PubMed  Google Scholar 

  34. Mahendran AK, Katz D, Opotowsky AR, Lubert AM. Exercise pathophysiology and testing in individuals with a Fontan circulation. CJC Pediatr Congenit Heart Dis [Internet]. 2023 [cited 2023 Mar 16]; Available from: https://www.sciencedirect.com/science/article/pii/S2772812923000027

  35. Wasserman K. The anaerobic threshold measurement to evaluate exercise performance. Am Rev Respir Dis. 1984;129:S35-40.

    Article  CAS  PubMed  Google Scholar 

  36. Casaburi R. Physiologic responses to training. Clin Chest Med. 1994;15:215–27.

    Article  CAS  PubMed  Google Scholar 

  37. Paolillo S, Veglia F, Salvioni E, Corrà U, Piepoli M, Lagioia R, et al. Heart failure prognosis over time: how the prognostic role of oxygen consumption and ventilatory efficiency during exercise has changed in the last 20 years. Eur J Heart Fail. 2019

  38. Novais ARB, Matecki S, Jaussent A, Picot M-C, Amedro P, Guillaumont S, et al. Hyperventilation during exercise in very low birth weight school-age children may implicate inspiratory muscle weakness. J Pediatr. 2012;160:415-420.e1.

    Article  Google Scholar 

  39. Keochkerian D, Chlif M, Delanaud S, Gauthier R, Maingourd Y, Ahmaidi S. Breathing pattern adopted by children with cystic fibrosis with mild to moderate pulmonary impairment during exercise. Respir Int Rev Thorac Dis. 2008;75:170–7.

    Google Scholar 

  40. Brat K, Stastna N, Merta Z, Olson LJ, Johnson BD, Cundrle I. Cardiopulmonary exercise testing for identification of patients with hyperventilation syndrome. PLoS ONE. 2019;14:e0215997.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. de Groot EP, Duiverman EJ, Brand PLP. Dysfunctional breathing in children with asthma: a rare but relevant comorbidity. Eur Respir J. 2013;41:1068–73.

    Article  PubMed  Google Scholar 

  42. Azab NY, El Mahalawy II, El Aal GAA, Taha MH. Breathing pattern in asthmatic patients during exercise. Egypt J Chest Dis Tuberc. 2015;64:521–7.

    Article  Google Scholar 

  43. Clark CJ, Cochrane LM. Assessment of work performance in asthma for determination of cardiorespiratory fitness and training capacity. Thorax. 1988;43:745–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Quanjer PH, Stanojevic S, Cole TJ, Baur X, Hall GL, Culver BH, et al. Multi-ethnic reference values for spirometry for the 3–95-yr age range: the global lung function 2012 equations. Eur Respir J. 2012;40:1324–43.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Cumming GR. Maximal exercise capacity of children with heart defects. Am J Cardiol. 1978;42:613–9.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We sincerely thank Annie Auer and Amandine Marquina (CPET laboratory technicians).

In Memoriam

This study is dedicated to Dr. Jonathan Rhodes, in memory of his extensive work on exercise physiology in children.

Funding

Not applicable.

Author information

Author notes

  1. Pascal Amedro, Stefan Matecki, Thibault Mura and Arthur Gavotto have contributed equally to this work.

Authors and Affiliations

  1. Department of Paediatric and Congenital Cardiology, M3C National Reference Centre, Bordeaux University Hospital, 1 Avenue Magellan, 33604, Pessac, France

    Pascal Amedro

  2. IHU Liryc, Electrophysiology and Heart Modelling Institute, INSERM 1045, Bordeaux University Foundation, Avenue du Haut Lévêque, 33600, Pessac, France

    Pascal Amedro

  3. PhyMedExp, CNRS, INSERM, University of Montpellier, 371 Avenue du Doyen Giraud, 34295, Montpellier, France

    Stefan Matecki, Johan Moreau, Luc Souilla, Marie Vincenti & Arthur Gavotto

  4. Paediatric Functional Exploration Laboratory, Physiology Department, Montpellier University Hospital, 371 Avenue du Doyen Giraud, 34295, Montpellier, France

    Stefan Matecki & Anne Requirand

  5. Department of Biostatistics, Clinical Epidemiology, Public Health, and Innovation in Methodology, Nimes University Hospital, University of Montpellier, Place du Professeur Debré, 30029, Nimes, France

    Taissa Pereira dos Santos & Thibault Mura

  6. Department of Paediatric and Congenital Cardiology, M3C Regional Reference CHD Centre, Montpellier University Hospital, 371 Avenue du Doyen Giraud, 34295, Montpellier, France

    Sophie Guillaumont, Gregoire De La Villeon, Johan Moreau & Marie Vincenti

  7. Paediatric Cardiology and Rehabilitation Unit, St-Pierre Institute, 371 Avenue de L’Évêché de Maguelone, 34250, Palavas-Les-Flots, France

    Sophie Guillaumont & Gregoire De La Villeon

  8. Department of Cardiology, Boston Children’s Hospital, 300 Longwood Avenue, Boston, MA, 02115, USA

    Jonathan Rhodes & Suellen Moli Yin

  9. Clinic for Paediatric Cardiology and Congenital Heart Diseases, German Heart Centre, Lazarettstrasse 36, 80636, Munich, Germany

    Alfred Hager & Julia Hock

  10. Clinical Research and Epidemiology Unit, Montpellier University Hospital, INSERM-CIC 1411, Clinical Investigation Centre, University of Montpellier, Montpellier, France

    Marie-Christine Picot & Helena Huguet

  11. INSERM, U1061, Neuropsychiatry: Epidemiological and Clinical Research, University of Montpellier, 39 Av. Charles Flahault, 34090, Montpellier, France

    Thibault Mura

  12. Paediatric Intensive Care Unit, Arnaud de Villeneuve Hospital, Montpellier University Hospital, 371 Avenue du Doyen Giraud, 34295, Montpellier, France

    Arthur Gavotto

  13. Department of Paediatric and Adult Congenital Cardiology, M3C National CHD Reference Centre, Bordeaux University Hospital, Haut-Leveque Hospital, Avenue de Magellan, 33604, Pessac Cedex, France

    Pascal Amedro

Authors
  1. Pascal Amedro
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Stefan Matecki
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Taissa Pereira dos Santos
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Sophie Guillaumont
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Jonathan Rhodes
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Suellen Moli Yin
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Alfred Hager
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Julia Hock
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Gregoire De La Villeon
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. Johan Moreau
    View author publications

    You can also search for this author in PubMed Google Scholar

  11. Anne Requirand
    View author publications

    You can also search for this author in PubMed Google Scholar

  12. Luc Souilla
    View author publications

    You can also search for this author in PubMed Google Scholar

  13. Marie Vincenti
    View author publications

    You can also search for this author in PubMed Google Scholar

  14. Marie-Christine Picot
    View author publications

    You can also search for this author in PubMed Google Scholar

  15. Helena Huguet
    View author publications

    You can also search for this author in PubMed Google Scholar

  16. Thibault Mura
    View author publications

    You can also search for this author in PubMed Google Scholar

  17. Arthur Gavotto
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

PA and AG directed, conceptualized, and designed the study, enrolled patients, supervised the data management, drafted the initial manuscript, reviewed, and revised the manuscript. SM, SG, GDLV, JM, AR, LS, MV enrolled patients, reviewed and revised the manuscript. JR, SMY, AH, JH, MCP, HH supervised the data management, reviewed, and revised the manuscript. TM, TDS supervised the data management, performed the statistical analyses, and reviewed and revised the manuscript. All authors approved the final manuscript as submitted and agree to be accountable for all aspects of the work.

Corresponding author

Correspondence to Pascal Amedro.

Ethics declarations

Ethics Approval and Consent to Participate

The study was conducted in compliance with the Good Clinical Practices and Declaration of Helsinki principles. It was approved by the Institutional Review Board of Montpellier University Hospital (2019_IRB-MTP_10-20) and registered on ClinicalTrials.gov (NCT04876209). Informed consent was obtained from all parents or legal guardians.

Consent for Publication

Not applicable.

Competing interests

The authors have no competing interests to declare.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1

. Supplementary Table 1. Enumerates the criteria from the ATS/ACCP statement on high-quality CPET assessment. Supplementary Table 2. Reports the best R2 values according to CPET and anthropometric parameters. Supplementary Table 3. Describes the distribution of CPET parameter values in children of normal weight.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Amedro, P., Matecki, S., Pereira dos Santos, T. et al. Reference Values of Cardiopulmonary Exercise Test Parameters in the Contemporary Paediatric Population. Sports Med - Open 9, 68 (2023). https://doi.org/10.1186/s40798-023-00622-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40798-023-00622-3

Keywords