PLoS ONE plos plosone PLOS ONE 1932-6203 Public Library of Science San Francisco, CA USA 10.1371/journal.pone.0159678 PONE-D-16-21094 Research Article People and placesPopulation groupingsAge groupsAdolescentsPeople and placesPopulation groupingsAge groupsBiology and life sciencesPhysiologyRespiratory physiologyMedicine and health sciencesPhysiologyRespiratory physiologyPeople and placesPopulation groupingsEthnicitiesGerman peopleMedicine and health sciencesPulmonologyAsthmaMedicine and health sciencesPulmonologyBronchitisMedicine and health sciencesPulmonologyRespiratory infectionsLower respiratory tract infectionsMedicine and health sciencesPulmonologyRespiratory infections Age Dependency of GLI Reference Values Compared with Paediatric Lung Function Data in Two German Studies (GINIplus and LUNOKID) Age Dependency of GLI Reference Values for Paediatric Lung Function Hüls Anke 1 * Krämer Ursula 1 Gappa Monika 2 Müller-Brandes Christine 3 Schikowski Tamara 1 von Berg Andrea 2 Hoffmann Barbara 1 4 Schuster Antje 5 Wisbauer Matthias 5 Flexeder Claudia 6 Heinrich Joachim 6 7 Schulz Holger 6 7 Berdel Dietrich 2 IUF-Leibniz Research Institute for Environmental Medicine, Düsseldorf, Germany Marien Hospital Wesel, Children’s Hospital and Research Institute, Wesel, Germany Medical School of Hanover, Department of Anaesthesiology and Intensive Care Medicine, Hanover, Germany Medical Faculty, Heinrich-Heine University of Düsseldorf, Düsseldorf, Germany Heinrich-Heine University, University Children’s Hospital, Düsseldorf, Germany Institute of Epidemiology I, Helmholtz Zentrum München - erman Research Center for Environmental Health, Munich, Germany Comprehensive Pneumology Center Munich (CPC-M), Member of the German Center for Lung Research, Munich, Germany Morty Rory Edward Editor University of Giessen Lung Center, GERMANY

The study received funding from GlaxoSmithKline GmbH & Co. KG, Aerocrine AB, MSD Sharp & Dohme GmbH, AstraZeneca GmbH, Novartis Pharma GmbH, Astellas Pharma GmbH, Munich, Germany, Ndd Medizintechnik AG, Zürich, Switzerland, Lions Club Wesel, Mead Johnson and Nestlé companies. This does not alter the authors' adherence to all the PLOS ONE policies on sharing data and materials.

Conceived and designed the experiments: UK AvB DB JH HS. Performed the experiments: MG CMB AvB AS MW JH HS DB. Analyzed the data: AH. Wrote the paper: AH UK MG CMB TS AvB BH AS MW CF JH HS DB.

 * E-mail: Anke.Huels@IUF-Duesseldorf.de 20 7 2016 2016 11 7 e0159678 25 5 2016 6 7 2016 2016 Hüls et al This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

A hallmark of the newly published GLI (Global Lungs Initiative) spirometric reference values is their "all-age" (3-95yr) predictive power, accomplished by incorporating non-linear age dependencies into modelling parameters. This modelling strategy is especially promising for the age range of puberty; however, the performance of GLI-values for adolescents is currently unknown. We calculated GLI-based z-scores for children/adolescents without apparent respiratory diseases from two different German studies, LUNOKID (N = 1943, 4–19 years) and GINIplus (N = 1042, 15 years) and determined the goodness of fit for specific age groups. We defined fit sufficient if the absolute mean of z-scores was <0.5. For children (<10yr) the mean GLI-based z-scores for FEV1 and FVC reached a good fit with mean z-scores for FEV1 between -0.11 and 0.01 and mean z-scores for FVC between 0.01 and 0.16, but larger deviations were observed in adolescents, especially boys (mean z-score -0.58 for FEV1 and -0.57 for FVC in GINIplus). The fit for FEV1/FVC was sufficient. GLI reference values provided reasonable estimates for the individuals enrolled in our studies, which span the age range of lung growth and development. However, we found that GLI-predictions overestimated lung volumes, especially those for German adolescent boys, which may, left unrecognised, lead to erroneous diagnosis of lung disease. Caution should be taken when applying these reference values to epidemiologic studies.

 LUNOKID was supported by GlaxoSmithKline GmbH & Co. KG, Munich, Germany; Aerocrine AB, Solna, Sweden; MSD Sharp & Dohme GmbH, Haar, Germany; AstraZeneca GmbH, Wedel, Germany; Novartis Pharma GmbH, Nuernberg, Germany; Astellas Pharma GmbH, Munich, Germany, Deutsche Atemwegsliga; Ndd Medizintechnik AG, Zürich, Switzerland, Lions Club Wesel. The GINIplus study was mainly supported for the first 3 years by the Federal Ministry for Education, Science, Research and Technology (interventional arm) and Helmholtz Zentrum Munich (former GSF) (observational arm). The 4-year, 6-year and 10-year follow-up examinations of the GINIplus study were covered from the respective budgets of the 5 study centres (Helmholtz Zentrum Munich (former GSF), Marien-Hospital Wesel, LMU Munich, TU Munich and from 6 years onward also from IUF - Leibniz Research-Institute for Environmental Medicine) and a grant from the Federal Ministry for Environment (IUF, FKZ 20462296). The 15-year follow-up examination of the GINIplus study was supported by the Commission of the European Communities, the 7th Framework Program (MeDALL project) and the Mead Johnson and Nestlé companies. The publication of this article was funded by the Open Access fund of the Leibniz Association. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Data Availability The authors confirm that, for approved reasons, some access restrictions apply to the data underlying the findings. The informed consent given by LUNOKID and GINIplus study participants does not cover posting individual data in public databases. All data underlying the findings described in the manuscript are archived after publication. Data are available upon request from the Steering Committee (LUNOKID: Prof. Dr. med. Dietrich Berdel, berdel.vonberg@t-online.de; GINIplus: Dr. Marie Standl, marie.standl@helmholtz-muenchen.de). All authors were affiliated with the cohort and had free access to the data underlying the findings described in the manuscript. Co-authorship is not necessary to gain access to the data, but a cooperation contract is required. Introduction

Spirometry is the most frequently used method with which to routinely assess lung function in clinical and research settings. However, the interpretation of results is crucially dependent on the use of appropriate reference data [1,2]. The use of historic reference values e.g. those developed by Zapletal et al. [3] in the 80s, although still in clinical use, are no longer the most appropriate [4,5].

Recently, the Global Lung Initiative (GLI) published multi-ethnic reference values with which to interpret spirometric results. These provide mean values and values for the lower limit of normal (LLN) from preschool to old age [5], achieved by including non-linear age dependencies, which are especially important for adolescents where the relationship between lung volume and height depends on maturity [6].

Global comparisons of GLI reference values with observations made from national studies have been published [7–10]. Using absolute mean z-scores of less than 0.5 to indicate sufficient fit, Hall and Thompson [7–10] showed that the newly published reference values were applicable to Australasian Caucasians. Likewise, overall mean GLI-based z-scores for the datasets (N>150) included in its reference population did not exceed 0.5 [11], indicating sufficient fit.

Global comparisons may however be misleading; a small global absolute mean of GLI based z-scores across ages does not necessarily preclude higher deviations for certain age groups. Systematic deviations from the proposed age relationship in GLI reference equations have already been investigated for adults, with two Scandinavian [12,13], and one Japanese report [14]. Thus far, only one study has investigated systematic z-score trends over age, for children and adolescents. Unfortunately that study predated the introduction of the GLI indices, and used historic reference equations instead [6]. Consequently, equivalent investigations for children and adolescents using GLI reference equations are still lacking.

Therefore, we aimed to investigate the fit between GLI derived means and LLN values during specific age periods in childhood and adolescence, by taking advantage of two independent German studies. The first incorporated children and adolescents from the LUNOKID-study (LUng function NOrmal values for KIDs in Germany), part of the original GLI reference population [11]. The second involved 15-year old adolescents from the GINIplus study (German Infant Nutritional Intervention plus environmental and genetic influences on allergy development study), which was not included in the GLI reference population.

 Materials and Methods Study design and study population

The first study analysed is the LUNOKID-study which aimed to generate new German-specific reference values for spirometry in children and adolescents. Details on the design, recruitment, lung function measurements, and quality criteria have been published [4,15–17]. Briefly, Caucasian children and adolescents from randomized kindergartens and schools, aged 4 to 19, were recruited from three cities in Germany (Wesel, Düsseldorf, Hannover) in 2008/2009. The children were tested in their schools and kindergartens, with parents completed a questionnaire on the respiratory health of their child and influencing factors (e.g. smoking, preterm birth).

The second study is the 15-year follow-up of GINIplus, which was designed to prospectively investigate influences on allergy development in a population based birth-cohort study, with nutritional intervention included for predisposed infants in the first months of life. Details on the design, recruitment and follow-up of this study have been published [18,19]. Briefly, a total of 5991 new-borns were recruited in obstetric clinics in Munich and Wesel, Germany, between September 1995 and July 1998. Follow-up occurred at the age of one, two, three, four, six, ten, and 15 years of age. Our analyses are based on lung function data from the 15-year follow-up.

Both studies were approved by the relevant ethics committees (Ethikkommission der Ärztekammer Nordrhein and Ethikkommision der Bayerischen Landesärztekammer) with written informed consent obtained from the parents of all participants.

 Lung function measurement

Detailed descriptions of the lung function measurements performed in LUNOKID and GINIplus have been published elsewhere [15,20]. Briefly, in both studies spirometry was carried out using an EasyOne handheld device (ndd-Medizintechnik AG, Zürich, Switzerland). Only non-smoking children with acceptable lung function according to the ATS/ERS criteria including visual control were considered for the analysis (N = 3205 in LUNOKID and N = 1628 in GINIplus) [1,2,15,20]. Furthermore we restricted our analyses to children/adolescents without apparent respiratory disease. Children with physician diagnosed asthma, spastic bronchitis or an infection on the day of investigation were excluded. Additionally, we excluded children/adolescents with a lower respiratory tract infection in the preceding 6 weeks from the LUNOKID reference population [17] and adolescents with a respiratory tract infection in the preceding week from the GINIplus population. In a sensitivity analysis, adolescents with a respiratory tract infection in the preceding 2–4 weeks were additionally excluded from the GINIplus dataset. Further we performed a sensitivity analysis in the GINIplus cohort, in which we excluded participants with a nutritional intervention in the first months of life. In additional sensitivity analyses we excluded from both cohorts participants who were exposed to passive smoke at home, who had atopic eczema or hayfever and who were born preterm (<37 weeks).

 Statistical analyses

The GLI [5] and LUNOKID [4] reference values are both based on regression model LMS. For a detailed description of the LMS model see Cole and Green [21], or Rigby et al. [22]. Due to the definition of the LMS model, we can calculate GLI or LUNOKID based reference values for a participant’s expected mean (M), and additionally for a participant’s expected coefficient of variation (S), and skewness (L) [5]. Given M, S and L, the standardized z-scores can be calculated as z = (YM)L-1L× S where Y is the observed FEV1, FVC, and FEV1/FVC. These z-scores have a standard normal distribution, and are independent of age for the reference populations. The coefficients needed for these calculation are provided in Quanjer et al. [5] (GLI), or Hüls et al. (LUNOKID) [4].

Prior to analyses, we checked datasets for outliers as our statistical tests are sensitive to these [11]. For the LUNOKID reference data we separately checked the FEV1, FVC, and FEV1/FVC values (adjusted for median height) for outliers, across 4 different age ranges; 3-<5 years, 5-<10 years, 10-<15 years, and 15-<20 years. Values more than 1.5 IQR (interquartile range) beneath, or above the lower and upper quartiles respectively, were excluded. The same procedure was also applied to the GINIplus data.

To investigate age-specific fits for the GLI reference values to measurements completed for LUNOKID, we divided our data into 4 different age groups, two were for children aged 3-<5 years, and 5-<10 years, and two were for adolescent groups of 10-<15 years, and 15-<20 years. Should the reference values prove to be a perfect fit for the study population, then the mean would be zero, and the standard deviation one, for every age group. Following the example of Hall [7] and Thompson [10], an absolute mean z-score of less than 0.5 was chosen to indicate sufficient fit. Two one-sided tests (TOST) for equivalence were conducted, and we tested our null hypothesis of a mean z-score outside of the interval [-0.5, 0.5]. If rejected, then the fit is sufficient. Furthermore, we calculated the percentage of observed lung function values lying below the predicted lower limit of normal (LLN). With a perfect fit, 5% of values should lie below the LLN. Additionally, we plotted the GLI based z-scores of the LUNOKID observations with age, together with the smoothed mean and 5th percentile using GAMLSS in R [23]. With a perfect fit to age, we would expect the smoothed mean to be constant zero, and the smoothed 5th percentile to be constant -1.64 (LLN). All analyses where conducted using R 3.1.1 [24].

 Results Description of the study population

The LUNOKID study population comprised 3205 participants with a mean age of 11.2 years (48.6% males), with the GINIplus study population comprising 1628 participants with a mean age of 15.3 years (47.2% males) (Table 1). With the exception of wheezy bronchitis, the prevalence of asthma and respiratory tract infections were comparable for both studies. In total, after excluding participants with apparent respiratory disease, our study population comprised 1943 (1922 without outliers) lung healthy participants from LUNOKID, and 1042 (1032 without outliers) from GINIplus, this is slightly more than 60% of the original study population. A more detailed description of these lung healthy participants is shown in Table 2. Table 2 further gives an overview about additional exclusion criteria that were used in the sensitivity analyses.

10.1371/journal.pone.0159678.t001 Description of study populations (LUNOKID and GINIplus).
LUNOKID GINIplus
Total N* 3205 1628
Age range, mean (min-max) 11.2 (4.0–19.0) 15.3 (14.7–16.8)
Sex, n (%) male 1557 (48.6%) 769 (47.2%)
BMI, mean (sd) 19.4 (4.9) 20.9 (3.2)
Asthma, n (%) 279 (8.7%) 190 (11.7%)
Wheezy bronchitis, n (%) 671 (20.9%) 76 (4.7%)
Airway infection (day), n (%) 440 (13.7%) 183 (11.2%)
Airway infection (last weeks), n (%) 241 (7.5%) 208 (12.8%)
Lung Healthy, n (%) 1943 (60.6%)† 1042 (64.0%)#

*: non-smokers with acceptable lung function according to the ATS/ERS criteria including visual control;

: LUNOKID: last 6 weeks; GINIplus: last week;

†: without asthma, wheezy bronchitis, infection (day and/or last 6 weeks);

#: without asthma, wheezy bronchitis, infection (day and/or last week)

 10.1371/journal.pone.0159678.t002 Description of the lung healthy study populations (LUNOKID and GINIplus).

Further exclusion criteria for sensitivity analyses.

LUNOKID GINIplus
Total N 1943 1042
Age range, mean (min-max) 11.0 (4.0–19.0) 15.3 (14.7–16.8)
Sex, n (%) male 901 (46.4%) 466 (44.7%)
Exposure to passive smoke at home, n (%) 656 (33.8%) 125 (12.0%)
Ever atopic eczema, n (%) 202 (10.4%) 189 (18.1%)
Ever hayfever, n (%) 119 (6.1%) 167 (16.0%)
Neither atopic eczema nor hayfever, n (%) 1622 (83.5%) 721 (69.2%)
Born preterm (<37 weeks), n (%) 181 (9.3%) 23 (2.2%)
Infection (last 2–4 weeks), n (%) n.a. 121 (11.6%)
Intervention group*, n (%) no intervention 490 (47.0%)

*: Intervention group: Nutritional intervention for predisposed infants in the first months of life due to study design;

: Lung healthy participants that were used for the main analysis (compare Table 1)

 LUNOKID data compared to the GLI reference values

Table 2 shows means and standard deviations of the GLI-based z-scores for FEV1, FVC and FEV1/FVC for our four age categories in LUNOKID. A graphical representation of the GLI based z-scores over age is given in Fig 1 for boys, and in S1 Fig for girls. The mean GLI-based z-scores for FEV1 and FVC reached a good fit for children below 10 years of age and decreased with age. From the age of 10 onwards z-scores for FEV1 reached values of -0.2 to -0.3 and z-scores for FVC reached values of -0.3 to -0.4. The effect was larger for boys than for girls. For the 315 boys in the ≥10 to <15 age group, the mean of the GLI based z-scores for FEV1 was -0.48, with a standard deviation of 0.81 (Table 3). Consequently, the null hypotheses of a deviation ≥ 0.5 could not be rejected (indicating insufficient fit). In comparison, the mean of the GLI based z-scores for the 334 girls in this age group was -0.20, with a standard deviation of 0.85 (indicating a sufficient fit). The fit of the mean GLI based z-scores for FEV1/FVC was sufficient for all age and gender groups.

10.1371/journal.pone.0159678.g001 Comparison of the GLI-z-scores (mean and lower limit of normal (LLN)) with these observed in boys from the LUNOKID data. 10.1371/journal.pone.0159678.t003 LUNOKID: GLI-based z-scores of all and of all male and female healthy children/adolescents grouped by age.
Age group N mean (sd) # %≤LLN
All
 FEV1 ≥4 & <6 290 -0.11 (0.76) *** 2.1%
≥6 & <10 561 0.01 (0.82) *** 2.0%
≥10 & <15 649 -0.33 (0.84) *** 4.9%
≥15 & <19 421 -0.40 (0.77) ** 6.2%
All 1921 -0.21 (0.82) *** 3.9%
 FVC ≥4 & <6 286 0.01 (0.76) *** 1.4%
≥6 & <10 560 0.16 (0.76) *** 0.4%
≥10 & <15 651 -0.23 (0.81) *** 3.5%
≥15 & <19 425 -0.30 (0.81) *** 4.7%
All 1922 -0.09 (0.81) *** 2.5%
 FEV1/FVC ≥4 & <6 291 -0.24 (0.98) *** 7.6%
≥6 & <10 553 -0.31 (0.83) *** 5.2%
≥10 & <15 641 -0.16 (0.83) *** 2.3%
≥15 & <19 424 -0.15 (0.85) *** 2.8%
All 1909 -0.21 (0.86) *** 4.1%
Males
 FEV1 ≥4 & <6 147 -0.11 (0.74) *** 2.0%
≥6 & <10 268 0.02 (0.81) *** 1.5%
≥10 & <15 315 -0.48 (0.81) 6.0%
≥15 & <19 160 -0.41 (0.79) 5.0%
All 890 -0.26 (0.82) *** 3.8%
 FVC ≥4 & <6 146 -0.02 (0.78) *** 1.4%
≥6 & <10 267 0.15 (0.77) *** 0.0%
≥10 & <15 316 -0.37 (0.79) ** 4.1%
≥15 & <19 160 -0.40 (0.83) 5.6%
All 889 -0.16 (0.83) *** 2.7%
 FEV1/FVC ≥4 & <6 147 -0.15 (0.98) *** 4.8%
≥6 & <10 261 -0.28 (0.75) *** 3.8%
≥10 & <15 311 -0.15 (0.83) *** 2.6%
≥15 & <19 160 -0.10 (0.77) *** 2.5%
All 879 -0.18 (0.83) *** 3.3%
Females
 FEV1 ≥4 & <6 143 -0.11 (0.79) *** 2.1%
≥6 & <10 293 0.01 (0.82) *** 2.4%
≥10 & <15 334 -0.20 (0.85) *** 3.9%
≥15 & <19 261 -0.38 (0.76) * 6.9%
All 1031 -0.17 (0.82) *** 4.0%
 FVC ≥4 & <6 140 0.05 (0.73) *** 1.4%
≥6 & <10 293 0.17 (0.75) *** 0.7% **
≥10 & <15 335 -0.09 (0.81) *** 3.0%
≥15 & <19 265 -0.24 (0.79) *** 4.2%
All 1033 -0.04 (0.79) *** 2.4% **
 FEV1/FVC ≥4 & <6 144 -0.34 (0.97) * 10.4%
≥6 & <10 292 -0.34 (0.90) ** 6.5%
≥10 & <15 330 -0.17 (0.82) *** 2.1%
≥15 & <19 264 -0.18 (0.89) *** 3.0%
All 1030 -0.25 (0.89) *** 4.8%

#: Two one-sided tests (TOST) for equivalence. Fit is sufficient, if the null-hypothesis of a mean z-score outside of the interval [-0.5, 0.5] is rejected

(*:p<0.05, **: p<0.01, ***: p<0.001).

For all of the age groups and lung function indices, the standard deviations of the GLI based z-scores were less than one, which influenced the fit of the GLI based LLN. Only 2.1% (1.4%) of the values for 4- to <6-year old children, and 2% (0.4%) of the values for 6- to <10-year olds, were below the LLN for FEV1 (FVC). This effect was not seen in the adolescent groups because the influence of their smaller standard deviation was offset by lower mean values.

 GINIplus data compared to the GLI and LUNOKID reference values

Table 4 shows the means and standard deviations of the GLI-based z-scores for FEV1, FVC, and FEV1/FVC, measured in GINIplus. For FEV1 and FVC, these means were below -0.5 indicating insufficient fit, with deviation more pronounced for boys than girls. Although the standard deviation was smaller than the expected 1, more adolescents than the expected 5% fell beneath the LLN. Again, this effect was more pronounced in males (15.0%≤LLN for FEV1 and 14.3%≤LLN for FVC) than females (both ≈8%). The difference in mean GLI-based z-scores between males and females was significant for FVC (t-test, p = 0.008). In contrast the mean GLI based z-scores for FEV1/FVC values were approximately zero, with standard deviations approximating one, indicating a sufficient fit.

10.1371/journal.pone.0159678.t004 GINIplus: GLI- and LUNOKID-based z-scores of all and of all male and female healthy participants.
GLI LUNOKID
N mean (sd) # %≤LLN mean (sd) # %≤LLN
All
 FEV1 1032 -0.58 (0.87) 11.5% *** 0.02 (1.12) *** 7.3% *
 FVC 1029 -0.57 (0.85) 10.9% *** -0.12 (1.06) *** 8.0% **
 FEV1/FVC 1038 -0.04 (0.97) *** 5.4% 0.28 (1.17) *** 6.2%
Males
 FEV1 461 -0.62 (0.91) 15.0% 0.07 (1.19) *** 9.1%
 FVC 461 -0.65 (0.88) 14.3% -0.13 (1.08) *** 9.1%
 FEV1/FVC 464 -0.03 (0.94) *** 4.1% 0.34 (1.20) ** 5.6%
Females
 FEV1 571 -0.55 (0.84) 8.8% * -0.02 (1.05) *** 5.8%
 FVC 568 -0.51 (0.82) 8.1% * -0.10 (1.04) *** 7.0%
 FEV1/FVC 574 -0.06 (0.98) *** 6.4% 0.23 (1.14) *** 6.6%

# Two one-sided tests (TOST) for equivalence. Fit is sufficient if the null-hypothesis of a mean z-score outside of the interval [-0.5, 0.5] is rejected

(*: p<0.05, **: p<0.01, ***: p<0.001).

Table 4 also shows the data following application of the LUNOKID derived reference values to GINIplus. Means and standard deviations of the LUNOKID based z-scores for observations from GINIplus indicated a sufficient fit for males and females.

 Sensitivity analyses: Use of stricter exclusion criteria

In the sensitivity analyses we excluded beside participants with asthma, wheezy bronchitis or infection (day and/or last 6 weeks) also participants who were exposed to passive smoke at home, who ever had atopic eczema or hayfever and who were born preterm (<37 weeks) (Table 5). These sensitivity analyses showed that our previous findings were not influenced by the exclusion criteria we chose in our main analysis. In a further sensitivity analysis in the GINIplus cohort, we additionally excluded participants with a respiratory tract infection in the last 2–4 weeks and participants with a nutritional intervention in the first months of life (S1 Table). These exclusion criteria had again no impact on our findings thus confirming our previous conclusions.

10.1371/journal.pone.0159678.t005 GLI-z-scores in GINIplus and LUNOKID in all lung healthy children (All) and in subgroups after using stricter exclusion criteria.
LUNOKID GINI
N mean (sd) # %≤LLN N mean (sd) # %≤LLN
All
 FEV1 1921 -0.21 (0.82) *** 3.9% 1032 -0.58 (0.87) 11.5% ***
 FVC 1922 -0.09 (0.81) *** 2.5% *** 1029 -0.57 (0.85) 10.9% ***
 FEV1/FVC 1909 -0.21 (0.86) *** 4.1% 1038 -0.04 (0.97) *** 5.4%
No passive smoke
 FEV1 1267 -0.20 (0.83) *** 3.9% 890 -0.56 (0.87) 11.1% ***
 FVC 1273 -0.10 (0.82) *** 2.4% *** 888 -0.56 (0.85) 10.5% ***
 FEV1/FVC 1266 -0.18 (0.86) *** 3.6% 896 -0.01 (0.96) *** 4.6%
No atopic eczema / hayfever
 FEV1 1600 -0.21 (0.83) *** 4.0% 713 -0.59 (0.87) 10.8% ***
 FVC 1602 -0.09 (0.81) *** 2.5% *** 714 -0.57 (0.85) 10.2% ***
 FEV1/FVC 1596 -0.22 (0.85) *** 4.1% 717 -0.05 (0.96) *** 5.0%
Not born preterm (≥37 weeks)
 FEV1 1725 -0.21 (0.82) *** 3.6% 993 -0.58 (0.87) 11.5% ***
 FVC 1727 -0.10 (0.80) *** 2.7% *** 991 -0.57 (0.85) 11.0% ***
 FEV1/FVC 1715 -0.20 (0.86) *** 3.9% 999 -0.05 (0.96) *** 5.4%

# Two one-sided tests (TOST) for equivalence. Fit is sufficient if the null-hypothesis of a mean z-score outside of the interval [-0.5, 0.5] is rejected

(*: p<0.05, **: p<0.01, ***: p<0.001).

 Discussion

We compared the age dependency of the GLI reference values with observations from the German LUNOKID study, and demonstrated a sufficient fit for children up to the age of ten years. However, for older children and adolescents, we found an ever-increasing divergence in volume measurements, with an insufficient fit for boys of 10 to 15 years of age. The same was true for adolescents in the GINIplus study, with a mean age of 15 years. However, a sufficient fit was found for FEV1/FVC. Furthermore, the standard deviations for the z-scores for FEV1 and FVC observed in the German studies were generally less than one.

Until now, several studies have compared GLI reference values with national reference data. The study conducted by Hall et al. [7] concluded that the GLI indices were applicable for Caucasians from Australia and New Zealand (n = 2066, age: 4–80 years). In our previous study on the applicability of the GLI reference values for longitudinal data of elderly women from the German SALIA study, we came to the same conclusion (n = 1726, mean age of 54 years at baseline examination) [25]. In contrast, Ben Saad et al. found that the GLI reference values were of less value for adult North African Caucasians in Tunisia (n = 489, age: 18–60 years) [8]. Age specific comparison studies have also been completed. For example, Backman et al. [12] evaluated the GLI values for a Swedish reference population (N = 501, age: 22–91) and found that mean GLI based z-scores for lung volume systematically increased with age, with the effect magnified for men versus women. The same trend was found for a population of Finish adults (n = 1000, age: 18–83) [13]. Again the z-scores for volume increased with age, reaching values of approximately 0.8 in men older than 70 years of age. Increasing z-scores between the ages of 40 and 80 were also found in a Japanese reference group (n = 20,336, age = 17–95) [14]. However, in contrast to the Scandinavian studies, the best fit was with the oldest age group, with z-scores of less than -0.5 occurring in the youngest group. For the Japanese study, secular changes in relative leg length (which was higher for the younger age group), partially explained the trend of deviation from the height-based GLI reference values; a similar effect could eventually explain the trends found in the Scandinavian studies, where height and age were also negatively correlated.

Similar age-specific comparison studies for GLI reference data during adolescence are however missing. A single study, which predated GLI and used older reference data, revealed a U-shaped deviation with the lowest z-scores obtained for males during puberty [6].

We found the lowest GLI based z-scores for FEV1 and FVC after the age of 10, particularly for males, which roughly coincides with the beginning of their growth spurt during puberty [26]. In contrast, the fit for FEV1/FVC was found to be sufficient. These data suggest a modifying effect for adolescents derived from a basic physiologic parameter. It is well appreciated that adolescents have a lower lung volume than expected for their height [27–29] during the time when growth is most rapid. However, height development differs even between Caucasian populations, especially for adolescents with respect to final height reached, as for the timing of the beginning of the growth spurt. Whereas median height up until the age of ten does not vary greatly between European populations, thereafter, values for Northern and Southern Europeans diverge. Median height is, for instance, 176 cm in an 18 year old Italian male compared to 182 cm for a Dutch male (and 180 cm for a German male) [30]. The age at puberty on the other hand is earlier in Italy than in the Netherlands [26]. GLI categorises Caucasians as a single ethnic group, and reflects mean age and height related lung development, but without accounting for age related deviations for specific Caucasian groups. This may lead to age and population specific mean deviations, and offers a tentative explanation for why the German LUNOKID based reference values fit the mean values of German adolescents from GINIplus better than the GLI reference values.

Whether the mean differences observed are clinically relevant is controversial. We adopted the criterion introduced by Hall [7] and Thompson [10], who considered a deviation of 0.5 z-scores (1/2 standard deviation) to be clinically relevant. We found that at least 50% of German adolescent reference males presented a lung volume below -0.5 z-scores. However, whether the absolute mean of the z-score is 0.4, 0.5, or 0.6, mostly will not influence diagnostic decision-making. Of more importance is the LLN, which is used as a diagnostic cut-off in clinical practice.

The height of the LLN is dependent on the mean and standard deviation. The standard deviations of our GLI based z-scores were smaller than expected. This may be partly explained by our strict visual quality control, which was applied in order to generate an accurate data set of normal values, suitable for clinical and epidemiological studies. Furthermore our study populations were ethnically homogeneous, with the same spirometric devices used for all measurements. These conditions contrast with those found for the combined reference data set for GLI, which is essentially a coalition of several different data sets. Thus the coefficient of variation from the GLI reference data set reflects inherent differences between centres due to their different measuring devices and population characteristics. The observed standard deviations of z-scores were also generally lower than expected in the Scandinavian reference populations [12,13]. Thus, even if the local and expected GLI means match, the local 5th percentile may still lie above the GLI based LLN, which may have consequences for the diagnoses of lung diseases. Therefore, an impaired lung function of German children younger than 10 years of age, might not be recognized using the LLN of the GLI reference values. However, this is not a problem for children and adolescents older than 10 years of age. For these children the LLN of the GLI reference values fits well to the 5th percentile observed in adolescents from the LUNOKID study.

Furthermore GLI based LLNs are not free from age or gender related bias when applied to German populations. German adolescents and males are more likely to be misdiagnosed with a lung disease when using GLI reference values than German children or females. This should be taken into account when assessing lung function for a German male adolescent. The quotient FEV1/FVC is free from this bias and should be the preferred metric in these cases. For clinical practice however it is important to keep in mind that a diagnosis should never be based on a single lung function test, other clinical parameters should be considered. Therefore, the differences in LLN should not lead to a false diagnosis if clinical decision-making is multifactorial.

The observed gender and age specific bias in GLI based z-scores are especially relevant for the application of reference values in epidemiological analyses. Such analyses sometimes use z-scores instead of the original spirometric measurements of lung volume to reduce the huge age and height related variance. LLNs are often used in epidemiology to objectively define obstructive or restrictive lung disease from lung function measurements. However if we applied GLI reference values to our GINIplus study, we would, for example, erroneously consider males to be more affected than females, and adolescents more so than children. For epidemiological applications it might therefore be better to use Nation/Germany-specific reference values (for Germany those derived from LUNOKID), where this bias is smaller, or to use no reference values at all but instead include age, height and gender in the analysis models.

 Conclusion

GLI reference values provide reasonable estimates for two independent German studies covering the age range of lung growth and development. Whilst universal availability and the continuity of reference values from childhood to adulthood, makes GLI reference data an attractive and valuable resource for clinical practice, we should bear in mind those patients for whom GLI- predictions may be slightly less accurate. In particular, the overestimates of lung volume in German adolescent boys could lead to erroneous diagnosis of lung diseases. We would advise an especially cautious approach when applying these reference values in the context of epidemiological studies.

 Supporting Information Comparison of the GLI-z-scores (mean and lower limit of normal (LLN)) with these observed in girls from the LUNOKID data.

(PDF)

 GINIplus (sensitivity analysis): GLI- and LUNOKID-based z-scores of all lung healthy 15-year old participants (All), participants without a respiratory tract infection in the last 2–4 weeks and participants with no nutritional intervention in the first months of life.

(PDF)

LUNOKID. The authors thank the children and their families for taking part in this study. The authors also thank following people for technical assistance: Christina Beckmann, Julia Bienen, Cornelia Bisdorf, Irene Groß, Christina Müller, Sandra Werth: Children's Hospital and Research Unit, Marien Hospital, Wesel, Germany; Heike Beermann, Marion Kliemt: Paediatric Pulmonology, Allergology and Neonatology, Hannover Medical School, Germany; Sabina Illi; University Children’s Hospital, LMU, Munich, Germany; Özgü Altin, Gisela Bartkowiak, Ursula Pfeiffer, Michaela Strempel: University Hospital, Heinrich-Heine University, Düsseldorf, Germany;

GINIplus. We thank the families for their participation in the study, the obstetric units for allowing recruitment and the GINI study team for its excellent work. The GINI Study Team: Institute of Epidemiology I, Helmholtz Zentrum München, German Research Center for Environmental Health, Neuherberg (J. Heinrich, I. Brüske, H. Schulz, C. Flexeder, C. Zeller, M. Standl, M. Schnappinger, M. Sußmann, E. Thiering, C. Tiesler); Department of Pediatrics, Marien-Hospital, Wesel (D. Berdel, A. von Berg); Ludwig-Maximilians-University of Munich, Dr von Hauner Children’s Hospital (S. Koletzko); Child and Adolescent Medicine, University Hospital rechts der Isar of the Technical University Munich (C.P. Bauer, U. Hoffmann); IUF—Leibniz Research Institute for Environmental Medicine, Düsseldorf (T. Schikowski, E. Link, U. Krämer).

 Abbreviations LLN

lower limit of normal (5th percentile)

 CoV

coefficient of variation

 sd

standard deviation

 FEV1

forced expiratory volume in 1 s

 FVC

Forced vital capacity

 IQR

interquartile range

 References Miller MR, Hankinson J, Brusasco V, Burgos F, Casaburi R, Coates A, et al. Standardisation of spirometry. Eur Respir J. 2005;26: 319–38. doi: 10.1183/09031936.05.00034805 16055882 Beydon N, Davis SD, Lombardi E, Allen JL, Arets HGM, Aurora P, et al. An official American Thoracic Society/European Respiratory Society statement: pulmonary function testing in preschool children. Am J Respir Crit Care Med. 2007;175: 1304–45. doi: 10.1164/rccm.200605-642ST 17545458 Zapletal A, Šamánek M, Paul T. Lung Function in Children and Adolescents: Methods, Reference Values. S Karger Ag; 1987. Hüls A, Krämer U, Gappa M, Müller-Brandes C, Seitner-Sorge G, von Berg A, et al. [New spirometric reference values for children and adolescents in Germany considering height and non-linear age effects: the LUNOKID-Study]. Pneumologie. 2014;68: 393. doi: 10.1055/s-0034-1367575 24914869 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. doi: 10.1183/09031936.00080312 22743675 Quanjer PH, Borsboom GJJM, Kivastik J, Merkus PJFM, Hankinson JL, Houthuijs D, et al. Cross-sectional and longitudinal spirometry in children and adolescents: interpretative strategies. Am J Respir Crit Care Med. 2008;178: 1262–70. doi: 10.1164/rccm.200805-700OC 18931334 Hall GL, Thompson BR, Stanojevic S, Abramson MJ, Beasley R, Coates A, et al. The Global Lung Initiative 2012 reference values reflect contemporary Australasian spirometry. Respirology. 2012;17: 1150–1. doi: 10.1111/j.1440-1843.2012.02232.x 22849658 Ben Saad H, El Attar MN, Hadj Mabrouk K, Ben Abdelaziz A, Abdelghani A, Bousarssar M, et al. The recent multi-ethnic global lung initiative 2012 (GLI2012) reference values don’t reflect contemporary adult's North African spirometry. Respir Med. 2013;107: 2000–8. doi: 10.1016/j.rmed.2013.10.015 24231283 Bonner R, Lum S, Stocks J, Kirkby J, Wade A, Sonnappa S. Applicability of the global lung function spirometry equations in contemporary multiethnic children. Am J Respir Crit Care Med. 2013;188: 515–6. doi: 10.1164/rccm.201212-2208LE 23947526 Thompson BR, Stanojevic S, Abramson MJ, Beasley R, Coates A, Dent A, et al. The all-age spirometry reference ranges reflect contemporary Australasian spirometry. Respirology. 2011;16: 912–7. doi: 10.1111/j.1440-1843.2011.01970.x 21401802 Quanjer PH, Stocks J, Cole TJ, Hall GL, Stanojevic S. Influence of secular trends and sample size on reference equations for lung function tests. Eur Respir J. 2011;37: 658–64. doi: 10.1183/09031936.00110010 20817707 Backman H, Lindberg A, Sovijärvi A, Larsson K, Lundbäck B, Rönmark E. Evaluation of the global lung function initiative 2012 reference values for spirometry in a Swedish population sample. BMC Pulm Med. 2015;15: 26. doi: 10.1186/s12890-015-0022-2 25887050 Kainu A, Timonen KL, Toikka J, Qaiser B, Pitkäniemi J, Kotaniemi JT, et al. Reference values of spirometry for Finnish adults. Clin Physiol Funct Imaging. 2015; [Epub ahead of print]. doi: 10.1111/cpf.12237 Quanjer PH, Kubota M, Kobayashi H, Omori H, Tatsumi K, Kanazawa M, et al. Secular changes in relative leg length confound height-based spirometric reference values. Chest. 2015;147: 792–7. doi: 10.1378/chest.14-1365 25254426 Müller-Brandes C, Krämer U, Gappa M, Seitner-Sorge G, Hüls A, Von Berg A, et al. LUNOKID: Can numerical American Thoracic Society/European Respiratory Society quality criteria replace visual inspection of spirometry? Eur Respir J. 2014;43: 1347–1356. doi: 10.1183/09031936.00058813 24232698 Berdel D, Beckmann C, von Berg A, Gappa M, Illi S, Krämer U, et al. Erhebung von Lungenfunktionsnormalwerten (Spirometrie) bei Kindern und Jugendlichen in Deutschland: Die LUNOKID-Studie. Atemw -Lungenkrkh. 2010; 395–404. Hüls A, Krämer U, Schuster A, Gappa M, Wisbauer M, Müller-Brandes C, et al. Inclusion of children with airway disease for the development of spirometry reference data. Eur Respir J. 2016;47: 1290–2. doi: 10.1183/13993003.01465-2015 27037320 Filipiak B, Zutavern A, Koletzko S, von Berg A, Brockow I, Grübl A, et al. Solid Food Introduction in Relation to Eczema: Results from a Four-Year Prospective Birth Cohort Study. J Pediatr. 2007;151: 352–358. doi: 10.1016/j.jpeds.2007.05.018 17889067 von Berg A, Krämer U, Link E, Bollrath C, Heinrich J, Brockow I, et al. Impact of early feeding on childhood eczema: Development after nutritional intervention compared with the natural course—The GINIplus study up to the age of 6 years. Clin Exp Allergy. 2010;40: 627–636. doi: 10.1111/j.1365-2222.2009.03444.x 20082618 Fuertes E, Bracher J, Flexeder C, Markevych I, Klümper C, Hoffmann B, et al. Long-term air pollution exposure and lung function in 15 year-old adolescents living in an urban and rural area in Germany: The GINIplus and LISAplus cohorts. Int J Hyg Environ Health. 2015;218: 656–65. doi: 10.1016/j.ijheh.2015.07.003 26220828 Cole TJ, Green PJ. Smoothing reference centile curves: the LMS method and penalized likelihood. Stat Med. 1992;11: 1305–1319. 1518992 Rigby RA, Stasinopoulos DM. Generalized additive models for location, scale and shape (with discussion). J R Stat Soc Ser C. 2005;54: 507–554. Stasinopoulos D, Rigby R. Generalized additive models for location scale and shape (GAMLSS) in R. J Stat Softw. 2007;23: Available: http://www.jstatsoft.org/v23/i07. R Development Core Team. R: A language and environment for statistical computing [Internet]. Vienna, Austria: R Foundation for Statistical Computing; 2014. Available: http://www.r-project.org/ Hüls A, Krämer U, Stolz S, Hennig F, Hoffmann B, Ickstadt K, et al. Applicability of the Global Lung Initiative 2012 Reference Values for Spirometry for Longitudinal Data of Elderly Women. PLoS One. 2016;11. doi: 10.1371/journal.pone.0157569 de Muinich Keizer-Schrama SMPF, Mul D. Trends in pubertal development in Europe. Hum Reprod Update. 7: 287–91. 11392375 Schrader PC, Quanjer PH, van Zomeren BC, Wise ME. Changes in the FEV1-height relationship during pubertal growth. Bull Eur Physiopathol Respir. 1984;20: 381–388. 6332654 DeGroodt EG, Van Pelt W, Borsboom GJJM, Quanjer Ph. H, Van Zomeren BC. Growth of lung and thorax dimensions during the pubertal growth spurt. Eur Respir J. 1988;1: 102–108. 3360086 Degroodt EG, Quanjer PH, Wise ME, van Zomeren BC. Changing relationships between stature and lung volumes during puberty. Respir Physiol. 1986;65: 139–153. doi: 10.1016/0034-5687(86)90046-0 3764119 Bonthuis M, van Stralen KJ, Verrina E, Edefonti A, Molchanova E, Hokken-Koelega ACS, et al. Use of national and international growth charts for studying height in European children: Development of up-to-date European height-for-age charts. PLoS One. 2012;7: 1–11. doi: 10.1371/journal.pone.0042506